特高含水期剩余油分布及形成机理

目前我国东部老油区普遍进入特高含水阶段,水驱开发效果变差,采收率明显下降,有效挖掘剩余油潜力是当前开发的重点所在.特高含水期油相主要以非连续状态存在于孔隙空间内,基于宏观认识的达西定律很难准确描述地下剩余油的真实流动状态.为摸清特高含水期剩余油分布状况,研究剩余油形成机理并挖掘剩余油潜力,基于N-S方程模拟油水两相在微尺度多孔介质内流动,运用相场方法实时追踪驱替过程中的两相界面,建立微观流动模型.研究特高含水期剩余油分布,并从力学角度分析残余油形成机理.结果表明特高含水期多孔介质内剩余油主要可以分为5种类型:孤立油滴、孔喉残余油、簇状非均质残余油、油膜、盲端残余油.毛管力,孔隙结构以及壁面润湿性是影响特高含水期剩余油的分布与类型的主要因素.孔隙尺度流动模拟得到的水驱曲线与矿场规律基本一致.特高含水期增注与注入表面活性剂均是通过驱替出簇状非均质剩余油,孔喉残余油以及孤立油滴达到提高剩余油采收程度的目的.簇状非均质剩余油的减少是采收程度提高的主要原因.矿场增注时机对经济开发有重要影响.该研究初步揭示了特高含水期剩余油分布及形成机理,为水驱油藏的后期开发提供有效指导.     

磁场-渗流场耦合作用下的铁磁流体多孔介质流动数值模拟

铁磁流体作为一种新型功能材料,由于同时具备磁性及流体性质被广泛应用于各个技术领域.本文从数值模拟理论上研究了利用磁场控制的铁磁流体驱油问题:针对外磁场作用下的铁磁流体多孔介质流动物理过程,将磁力项引入铁磁流体相运动方程从而耦合磁场-渗流场给出其流动方程,并采用全隐式有限体积法形成相应数值离散格式;通过二维均质填砂平板流动实验验证了模型和算法的正确性,在此基础上,对比了平面及纵向非均质多孔介质中注水及注铁磁流体的驱油效果,计算结果表明通过控制铁磁流体的驱替路径可以提高驱替波及范围继而提高采收率,为提高原油采收率提供了新的思路和方法.     

纳米颗粒表面活性剂对多孔介质中黏性指进现象的抑制

多孔介质中不稳定的两相驱替会导致黏性指进现象,对实际工业中的两相流动过程产生不利影响.已有文献证明水相中的表面功能化纳米颗粒和油相中的功能化基团聚合物可以在油水界面发生反应生成纳米颗粒表面活性剂,降低界面张力并诱发界面黏弹性.本文提出使用纳米颗粒表面活性剂抑制多孔介质中黏性指进的新方法,通过微流体可视化实验,研究了纳米颗粒表面活性剂对多孔介质中驱替界面特性的影响规律.结果表明,纳米颗粒表面活性剂可以有效抑制多孔介质中的黏性指进现象,驱替效率对比纯水纯油提高了约2.5倍.通过分析基于不同毛细数与黏度比下的驱替模式分布相图,得到了多孔介质中纳米颗粒表面活性剂稳定驱替和黏性指进两种驱替模式的转变界限,证明纳米颗粒表面活性剂驱替对比纯水纯油大幅提高了黏性指进出现的临界毛细数和黏度比,在不同的参数范围内对黏性指进都有显著的抑制作用,为多孔介质中黏性指进的抑制提供了新的思路.     

用应力回复的实验方法来研究水解聚丙烯酰胺的黏弹特性

聚合物溶液在强化采油工业中有着十分重要的应用, 其独特的黏弹特性在提高微观驱油效率方面发挥了巨大的作用. 采用应力回复的实验方法对水解聚丙烯酰胺溶液的黏弹性进行了较为系统的研究. 研究表明水解聚丙烯酰胺溶液的黏弹性随浓度的增大而不断地增强; 随阴离子表面活性剂十二烷基硫酸钠(SDS)含量的增加, 黏弹性降低; 氯化钠的加入对黏弹性有破坏作用.     

页岩压裂中压裂液返排率低的孔隙尺度模拟与解释

针对页岩压裂过程中,压裂液返排效率普遍很低的现象,基于高分辨率页岩岩样SEM(scanning electron microscope)电子显微镜扫描图像,利用马尔科夫链蒙特卡罗(Markov Chain Monte Carlo)方法重构得到页岩的三维数字岩心,采用具有高密度比的格子Boltzmann模型,从孔隙尺度来模拟页岩中气水两相驱替过程.首先通过计算表面张力和相对渗透率来验证模型的准确性,然后模拟页岩数字岩心中的油水两相流动,页岩中首先饱和水(气)相,然后从一端注入气(水)相,模拟终止条件为驱替相在出口端发生突破,气水运动黏度比和密度比分别设置为10:1和1:1000.水驱气过程中发生突破时,水相的饱和度为70%,而气驱水过程中发生突破时,气相的饱和度只有4.5%,给出了三维数字岩心中驱替相分布,在气驱水的过程中发生突破时,大部分水被滞留在页岩孔隙中,从而解释了页岩水力压裂中,压裂液返排效率低于10%的现象.     

高含水期水驱特征曲线上翘机理及其影响因素

水驱特征曲线对预测水驱油藏可采储量及开发动态具有重要作用,其理论推导的基础是油水相对渗透率比值(K_(ro)/K_(rw))与含水饱和度(S_w)在半对数坐标下呈线性关系.然而大量矿场实践及岩心实验研究显示,油田开发进入高含水期后, K_(ro)/K_(rw)与S_w在半对数坐标下偏离直线关系,导致水驱特征曲线出现上翘现象而无法准确预测开发动态.为摸清水驱特征曲线上翘的本质原因,本研究基于Navier-Stokes方程模拟油水两相在微观多孔介质中的驱替过程,采用相场方法实时描述油水两相界面变化.通过微观油水两相数值模拟研究了壁面润湿性、油水黏度比等因素对油水在高含水阶段流动状态的影响,并分析了这些因素对K_(ro)/K_(rw)与S_w在半对数坐标下偏离直线关系时临界含水率的影响.研究结果表明,高含水期剩余油由连续相转变为非连续相,油相流动能力快速下降是水驱特征曲线出现上翘的内在原因;临界含水率在不同润湿性多孔介质内存在差异,渗吸作用导致水湿多孔介质中的临界含水率高于油湿多孔介质;高油水黏度比条件下的黏性指进现象会强化注入水的无效驱替,使水驱特征曲线在相对较低含水率时偏离直线.该研究初步揭示了油藏润湿性及油水黏度比对水驱特征曲线上翘拐点的影响,为高含水油藏开发提供有效指导.     

热弹性耦合系统解的渐近性质

一、引言 热弹性耦合理论是一个引入注目的研究课题。在热弹性材料中,物体上的温度不仅取决于热源、有关的热力学物性系数和换热边界条件,而且还受到弹性变形应变率的影响。另一方面,弹性运动也受到温度的影响,因此热传导方程和弹性运动方程必须耦合求解,其数学模     

SDS在水溶液中聚集的自扩散的NMR研究

目前在强化采油工艺中使用复合驱替液是一个有效的措施.因为它能改变多孔岩石中地下水、油和岩石中微细孔道表面间的界面性质,从而达到最大的驱替效果.这种复合驱替液包含着一定量的表面活性剂和水溶性的高聚物.表面活性物在水溶液中有自发聚集的性质,在驱油过程中可产生油-水间的超低界面张力(IFT),从而提高驱替效果.核磁共振(NMR)技术可以同时测量多组分体系中的不同种类的分子在溶液中的Brown运动的信息(即自扩散系数),从而得到聚集体的大小和聚集过程的信息,这对研究超低IFT的产生过程,从而有目的的改进工艺、提高驱油效率、降低成本有很大意义.传统的用NMR测量自扩散系数的方法是用梯度场-自旋回波(PGSE)的方法,由于自旋回波存在以下几个缺点:(ⅰ)由于脉冲期间磁化强度矢量在横向演化,所以自旋弛豫是依赖T_2弛豫过程.(ⅱ)当自旋耦合存在时,回波将受到J-耦合的调制,使谱图发生畸变.(ⅲ)引进第一个梯度脉冲之时会产生涡流,从而使     

页岩微裂缝内气-水两相流动规律

气-水两相流动广泛存在于页岩气藏压裂液返排及页岩气生产阶段,直接决定页岩气藏压裂及开发效果.页岩微裂缝是储层流体运移的主要通道.然而由于物理实验尺度及精度的限制,准确模拟页岩微裂缝中的气-水两相流动存在巨大挑战.为探究页岩微裂缝中气-水两相流动规律,本文系统考虑了气-水两相流体多重微尺度效应,包括:(1)气相滑移;(2)水相滑移及边界层黏度变化;(3)气驱水条件下,存在沿壁面流动的水膜.建立气-水两相流体微尺度流动模型,采用侵入逾渗判断两相分布,求解真实页岩微裂缝内气-水相对渗透率,并验证模型的正确性.研究结果表明,微裂缝开度小于3μm时,微尺度效应影响不可忽略;微裂缝开度及水膜对气-水相对渗透率的影响均取决于各相流体所占流动空间的相对大小;气藏开采过程中压力降低,导致气相表观渗透率升高及微裂缝开度减小,均会造成气相相对渗透率减小;水膜的存在对气相流动能力造成影响的微裂缝开度界限为0.65μm.该研究揭示了页岩微裂缝内气-水两相流动规律,为评价页岩气井压裂及生产效果提供了理论基础.     

内陆浅层咸水区渗透性变异过程

选取鲁北平原区浅层粉土、浅层咸水和大气降水为供试材料,进行室内雨水-咸水驱替试验,研究驱替过程中的含水介质渗透性变化特征,探讨内陆浅层咸水区水动力、水文地球化学过程对含水层渗透性变异影响.试验结果表明:雨水驱替咸水饱和土柱过程中,随雨水注入,盐分突变,渗透系数由开始的6.49×10-4 cm/s最终降至2.28×10-4 cm/s,并维持相对稳定;第188天后咸水再次驱替土柱,且完全穿透土柱后,含水介质孔隙度由开始的43.02%降至39.73%;在雨水驱替过程中,驱替液黏滞性变化是影响渗透系数变化的因素之一.研究认为,内陆咸水区含水介质黏土矿物含量高,水化学与水动力学条件发生变化,Ca/Mg-Na的离子交换作用、黏粒大量膨胀和释放、扩散引起孔隙度降低是造成含水介质的渗透性下降的主因,咸雨水驱替过程中渗透性的变化具有不可逆性.含水介质水文地质参数的变异对于内陆浅层咸水形成和演化有重要影响.     

线性剪切流中弹性界面液滴变形的惯性作用

弹性界面液滴(由弹性界面包裹而成的液滴)在流场中的变形等动力学行为是深入理解液滴群悬浮液流变特性以及发展基于液滴的物质输送技术的重要基础.流体惯性、液滴界面弹性、流体黏性等力学作用相互耦合,共同控制弹性界面液滴的动力学行为.已有文献以Stokes流假设进行简化,集中研究了忽略流体惯性条件下弹性界面液滴的变形动力学行为.本文结合Front Tracking Method和Finite Element Method,在传统两相流模拟方法中引入薄膜弹性力学理论,发展了能同时考虑液滴界面弹性和流体惯性的两相流数值模拟方法,研究了颗粒雷诺数(Re)对线性剪切流中弹性界面液滴变形动力学行为的影响.结果表明:雷诺数较高时,弹性界面液滴的瞬态变形出现振荡,且随雷诺数增大,变形振荡的幅值和周期都增大,液滴瞬态变形时间变长;随雷诺数增大,液滴的最大变形和稳态变形均增大,且液滴形状也有显著变化.由此可见,在Re1的条件下,流体惯性对弹性界面液滴变形动力学行为的影响不可忽略.     

潜热型功能热流体层流强化传热分析新模型

建立了分析潜热型功能热流体的层流强化传热的新模型. 该模型的主要特征是: (1) 应用新推导的潜热型功能热流体定压比热公式; (2) 考虑了真实的相变传热过程, 即, 相变微胶囊与单相工作流体间换热, 壁面与两相混合流体间换热; (3) 用新建立的轴对称双倒易边界元和有限差分联合的数值方法求解该数理方程. 新模型的可靠性用已有实验结果进行了验证. 得到了一些新的物理结果, 如定压比热在空间与时间随相变传热过程变化的特征, 捕获了相变界面随时间推进图像, 等等, 还数值分析了影响潜热型功能热流体强化传热的几个主要物理因数, 对该功能流体强化传热机理获得了一些新认识.     

基于SPH方法的液滴撞击可湿润固壁模拟

液滴撞击固体壁面是一种常见的自然现象和工业过程,在这些现象和过程中液滴运动往往受到壁面可湿润性的影响.本文提出了一种可湿润性固壁边界条件的计算方法,即假设壁面粒子的亲水性以及毛细吸附作用统一表现为对支持域内流体粒子的吸附力.吸附力的大小与流体压强、材料饱和度及亲水性有关.基于光滑粒子流体动力学(SPH)方法,模拟了静态液滴在不同可湿润性壁面上由变形至稳定的过程.考虑壁面可湿润性的影响,模拟了液滴撞击疏水壁面的过程.结合试验数据,分析了壁面可湿润性对液滴运动过程的影响.研究表明:根据模拟液滴静态接触角的变化特点,本文可湿润性固壁边界条件可以较好地反映出壁面可湿润性;液滴撞击输水表面的模拟数据与试验结果趋势上吻合良好,液滴回弹时在壁面可湿润性的影响下形成液柱;壁面可湿润性对撞击后液滴的铺展过程影响很小,在该阶段壁面摩擦力起主要作用,回缩和回弹阶段壁面可湿润性对液滴运动特征的影响较明显.     

山东胜利油田沾3油藏微生物组特征

以山东胜利油田沾3区注入水和采出液为研究对象,利用16S r RNA高通量测序技术,分析了油相样本、水相样本以及注入水样本的微生物群落结构.结果发现沾3油藏微生物多样性较高,注入水样本中细菌群落多样性最高,古菌群落多样性最低;而采出液油相样本的微生物群落多样性较采出液水相样本低.采出液微生物组的特征与注入水中微生物组的特征存在较大差异,具有明显的生境特异性.采出液油相样本与水相样本中主要的细菌类群包括γ-变形菌纲、ε-变形菌纲、α-变形菌纲、放线菌纲和梭状芽孢杆菌纲.在属水平上,采出液以假单胞菌占优势,而注入水中的优势菌属为硫单胞菌.采出液中古菌主要类群为甲烷微菌纲、甲烷杆菌纲、嗜盐菌纲、热原体纲、热变形菌纲和热球菌纲,尤其以鬃毛甲烷菌科和甲烷杆菌为主.对注入水和采出液水相中的油藏微生物与环境理化因子进行冗余分析(RDA)发现,注入水细菌和古菌群落均主要受SO_4~(2-)浓度影响,采出液水相样本中的细菌群落主要受Ca~(2+)浓度和pH等影响,而古菌群落则主要受Cl~-, Mg~(2+)和HCO_3~-浓度影响.该研究结果揭示了微生物在高温油藏环境中的空间分布特征及受控因素,为今后定向调控油藏微生物群落、进而提高微生物驱油的采油率奠定了基础.     

二氧化碳气体的流动及反应特性研究: 从微管到多孔介质

用5 μm微管、人造岩心及天然岩心研究了CO₂气体的流动特性及与岩心作用后对渗透率的影响规律. 结果表明, CO₂气体在微管中的流动速度明显比N₂快; CO₂在饱和水的岩心(含碳酸盐)中流动时, 其渗透率随着注入量的增加而增大, 表明CO₂的溶蚀作用导致了岩心孔径的增加, 用扫描电子显微镜(SEM)也验证了溶蚀扩孔作用. CO₂气体在微管中具有较高的流速是由于尺度效应和压缩效应的综合反映, 而其在水中溶解引起的界面层水分子扩散速度增大, 导致孔隙壁面的水膜厚度减小, 其水溶液流动的有效孔径增大. 上述两种结果表明, CO₂在驱油过程中具有良好的注入能力, 是低渗储层开发过程中能量补充的一种很好的驱替流体, 但也造成其在地层中窜流或散逸的可能性增大.     

基于时驱硬球算法与格子玻尔兹曼方法的颗粒流体系统直接数值模拟

实现了一种直接数值模拟颗粒流体系统的耦合算法, 颗粒间相互作用由时驱硬球算法描述, 而流体的控制方程采用格子玻尔兹曼方法求解, 流固耦合用浸入运动边界法实现.该方法使用欧拉网格求解流场, 拉格朗日网格跟踪颗粒, 避免了非结构化贴体网格方法需要重新划分网格的问题. 通过模拟两个圆形颗粒在黏性流体中的沉降过程, 成功地复现了经典的Drafting-Kissing-Tumbling(DKT)过程, 验证了耦合算法的有效性.     

油基水基及部分分散体系水合物堵管机理

多相体系水合物堵管规律作为实现水合物风险控制的关键,一直备受流动安全保障行业的关注.本文结合众多研究者的实验研究及理论分析,总结了部分实验装置及实验方式,概述了油基体系、水基体系及部分分散体系等多相流动体系的水合物堵管规律,并对后续研究的开展提出建议.本文指出不同体系水合物的生长及堵管机理不同,油基体系堵管主要受聚并、壁面黏附及着床等机理的控制;水基体系堵管主要受着床机理控制,堵管趋势最低;部分分散体系堵管主要受膜生长机理控制,堵管趋势最高.建议下一步工作可围绕以下几点展开:油基体系壁面黏附及着床的理论建模,水基体系水合物生长模型的建立,部分分散体系自由水的量化表征以及考虑不同多相流流型对堵管规律的影响.     

考虑非弹性碰撞的低浓度固液两相流动理学模型

用修正的BGK模型模拟粒间非弹性碰撞效应, 用Fokker-Planck扩散算子描述湍流-颗粒作用, 建立了低浓度固液两相流颗粒相的模型动理学方程和宏观水动力学方程组. 运用 Chapman-Enskog迭代法获得动理学方程的二阶近似解, 建立了颗粒相脉动应力和脉动能传导通量的显式本构关系. 建立的动理学模型在颗粒近乎弹性条件下简化为已有的模型, 为固液两相流模拟提供了一个实用模型.     

电场作用下电流变液滴的变形及力学行为

外加电场作用下的液滴变形及力学行为是电流体动力学研究的重要内容.基于volume-of-fluid(VOF)方法,将电场力作为源项添加到流体运动的Navier-Stokes方程中,提出了电流体流场与电场双向耦合的数值方法,分别研究了外加均匀电场与非均匀电场作用下,中性漏电液滴和带电液滴的变形/运动及其力学行为.小变形条件下,数值模拟得到的外加均匀电场中性漏电液滴的变形系数接近于理论值,验证了本文方法的正确性.研究结果表明,电场作用下液滴与周围介质物理属性的差异导致液滴表面不同的自由电荷再分配形式及其变形状态;相对于电场方向,中性漏电液滴可能会发生"扁长型"或"扁平型"变形.均匀电场下,中性漏电液滴内部形成稳定的回转运动,液滴不会发生宏观运动;对于存在净电荷的液滴,由于库仑力的作用,液滴不仅会发生变形,同时也会沿电场线方向运动.非均匀电场下,中性漏电液滴与带电液滴都会沿电场线运动,并发生不同的变形.本文所提出电流体流场与电场双向耦合的数值方法为静电喷雾、电泳等复杂的工程电流体动力学研究提供了理论基础.     

模拟运动流体中柔性体耦合振动的数值方法

介绍了一种模拟流体诱发柔性体运动的数值方法. 该方法采用二维面元法计算流体作用力, 用模态叠加法求解欧拉梁的运动控制方程, 通过迭代来求解流固耦合问题. 利用该方法模拟了单个以及两个并行、串行排列的柔性梁在运动流体中的耦合运动. 计算结果表明, 均匀来流中单个柔性梁的稳定边界和流体诱发振动都与前人结果一致. 对于两并行排列的柔性体, 当无量纲间距H<0.25时, 两柔性体同向摆动; H>0.25时, 两柔性体反向摆动, 并且摆动频率在模式转换时发生突变; 当间距很大时(H>1.0), 两并排柔性体之间的耦合作用明显变弱, 两个物体运动呈解耦状态, 各自按照单个柔性梁在来流中的运动形态运动. 相同均匀来流中两串行柔性梁的数值模拟结果显示, 当排列间距很小时, 上游物体受到的阻力减小, 下游物体阻力增大. 本文计算的所有结果都与前人实验结果定性一致, 证明本文提出的方法可以用来模拟不同方式排列的多个柔性体在流体中的耦合运动.