一种研究多孔介质中水气分散体系运移机制的仿真模拟新方法.pdf

第30卷第5期油气地质与采收率Vol.30，No.52023年9月Petroleum Geology and Recovery EfficiencySep.2023引用格式：王富琼，王秀宇，单学军，等.一种研究多孔介质中水气分散体系运移机制的仿真模拟新方法J.油气地质与采收率，2023，30（5）：130-138.WANG Fuqiong，WANG Xiuyu，SHAN Xuejun，et al.A new simulation method to study transport mechanism of water-gas dispersion system in porous mediaJ.Petroleum Geology and Recovery Efficiency，2023，30（5）：130-138.一种研究多孔介质中水气分散体系运移机制的仿真模拟新方法王富琼1，王秀宇1，单学军2，温志远1，陈秀林1，杨胜来1（1.中国石油大学（北京）油气资源与探测国家重点实验室，北京 102249；2.中国石化集团国际石油勘探开发公司，北京 100029）摘要：水气分散体系可有效改变地层油渗流通道，抑制驱替流体的窜流，扩大水气分散体系在油层中的波及体积，提高剩余油的动用效果。由于难以直接观察多孔介质中水气分散体系的流动，因此对于水气分散体系在多孔介质中的运移机理研究较少。笔者提出一种研究水气分散体系在多孔介质中运移的仿真模拟新方法，基于水平集方法和气液两相流，使用多物理场耦合数值计算软件COMSOL Multiphysics进行研究。对水气分散体系在孔喉中的贾敏效应、聚并和卡断机理进行了仿真研究，探讨了水气分散体系运移过程中孔喉压力的影响因素，分析分散相气泡聚并过程的形态变化，研究水气分散体系在孔喉中发生卡断现象的动态特征，结果表明数值模拟可以直观的体现水气分散体系在多孔介质中的运移机制。关键词：水气分散体系；水平集方法；气液两相流；多孔介质；运移机制文章编号：1009-9603（2023）05-0130-09DOI：10.13673/j.pgre.202207012中图分类号：TE341文献标识码：AA new simulation method to study transport mechanism ofwater-gas dispersion system in porous mediaWANG Fuqiong1，WANG Xiuyu1，SHAN Xuejun2，WEN Zhiyuan1，CHEN Xiulin1，YANG Shenglai1（1.State Key Laboratory of Petroleum Resources and Exploration，China University of Petroleum（Beijing），Beijing City，102249，China；2.SINOPEC International Petroleum Exploration and Production Corporation，Beijing City，100029，China）Abstract：The water-gas dispersion system can effectively change flow channels of formation oil，inhibit the channeling of displacement fluid，expand the swept volume of the water-gas dispersion system in reservoirs，and enhance the producing effect of theremaining oil.Since it is difficult to directly observe the flow of the water-gas dispersion system in porous media，there are few studies on the transport mechanism of the water-gas dispersion system in porous media.The author proposes a new simulation method tostudy the transport of the water-gas dispersion system in porous media.The level set method and the gas-liquid two-phase flow areadopted，and the coupled numerical calculation software COMSOL Multiphysics is employed for research.The mechanisms of theJamin effect，coalescence，and division of the water-gas dispersion system in the pore throat are simulated and studied.The influencing factors of dynamic pressures during the transport of the water-gas dispersion system are discussed，and the morphologicalchanges of the dispersed-phase bubble during the coalescence process are analyzed.The dynamic characteristics of the water-gas dispersion system during division in the pore throat are studied.The result shows that numerical simulation can intuitively reflect thetransport mechanism of the water-gas dispersion system in the porous media.Key words：water-gas dispersion system；level set method；gas-liquid two-phase flow；porous media；transport mechanism油藏进入高含水或特高含水期后，储层动用性差，层间矛盾显著，高渗透层注入水窜流严重1-6。针收稿日期：2022-07-20。作者简介：王富琼（1998），女，湖北恩施人，在读硕士研究生，从事油气田开发工程研究。E-mail：。通信作者：王秀宇（1976），女，辽宁庄河人，副教授，博士。E-mail：。基金项目：国家自然科学基金项目“超深层碎屑岩油气藏渗流物理基础研究”（51774300）。第30卷第5期王富琼等.一种研究多孔介质中水气分散体系运移机制的仿真模拟新方法 131 对这一开发问题，许多学者提出了一种利用水气分散可调阻力驱替体系开发水驱油田的新理念7-10，水气分散可调阻力驱替体系可分为3类：针对低渗透油藏的水气分散体系。适用于中高渗透油藏的泡沫驱体系。适用于高含水油藏的耐温、耐盐乳状液驱替体系。水气分散驱油体系通过气泡变形以及渗流阻力增加等方式有效改变渗流通道11，提高剩余油的动用效果。众多学者通过实验对水气分散体系提高原油采收率机理开展了研究工作。俞宏伟等通过岩心驱替实验观察水气分散体系的驱替特征，评价其驱油效率，并研究水气分散体系提高原油采收率的能力8。端祥刚研发渗流阻力可调的系列分散体系，证明可调阻力驱油体系的剖面调整性能，能够同时提高中低渗透层和高渗透层采收率，实现剖面均匀推进12。可调阻力水气分散体系及其设计方法的建立，为进一步提高水（气）驱波及效率，从而提高原油采收率提供了新的技术思路和理论基础。LIU等通过孔板喷射法制备水气分散体系，并分析不同因素对气泡数量和大小的影响13。陈兴隆等建立以气泡直径、气泡上升速度等为关键参数的水气分散体系性能评价方法，从理论上评价了超声波振荡生成方法生成的微米级气泡的特点11。根据长度低渗透岩心驱油实验，证明水气分散体系可通过气泡形变及调节渗流阻力等方式有效扩大波及体积，提高剩余油动用效果。然而，对于水气分散体系的实验操作难度很大，且难以直接观察孔隙尺度下两相流运移过程，而基于数值模拟方法可以更好地分析该体系在多孔介质中的运移机理。在孔隙尺度下，水气分散体系是涉及气液界面运动的两相流，其气液界面张力对流动的影响无法忽略，对水气分散体系流动过程模拟难点在于对界面运动的描述。在孔隙尺度上研究气液界面运动的方法有流体体积（Volume of Fluid，简称VOF）方法13、水平集（Level Set）方法14-17以及相场方法（Phase-Field method）18-19。VOF方法最早是由HIRT等在 1981年提出的，基本思想是通过在全流场定义每个欧拉网格单元内某种流体的体积比函数F来确定界面的位置20。部分学者通过VOF方法研究气泡在孔隙和喉道中的运移过程21-23。PAN等利用VOF方法研究了微通道中气泡传输的动力学24。KANG等利用VOF方法对微蛇形通道中的气泡生成和流动进行建模，通过与大量实验结果比较，讨论和验证了由气液流动、表面张力和黏度影响引起的气泡产生过程25。FENG等通过VOF方法研究了低雷诺数下2个串联气泡之间的聚并过程，观察气液两相的速度和液膜变化，聚并过程聚结可分为排液（薄膜变薄）和合并（破膜）26。相场方法以金兹堡-朗道理论为基础，利用一个有一定厚度的两相间的界面代替实际模型中的尖锐界面，由于流体间扩散效应的存在，相界面的变化为对流与扩散的共同作用27。MUKHERJEE 等使用 SIMPLER 方法求解完整的Navier-Stokes方程以及连续性和能量方程，利用水平集方法捕获气液界面，并进行实验以验证数值模型，且数值与实验结果之间显示出很好的一致性28。梁猛等通过相场方法，研究微观尺度下液滴变形和运动规律，分析表面张力和液滴大小对液滴变形的影响，结果表明当表面张力变小，较大半径的液滴更容易发生较大的形变29。水平集方法考虑了两相的接触，已成功应用于模拟液体中的气体流动。朱前林等基于Navier-Stokes方程与守恒的水平集方法，对孔隙介质中气泡尺寸对流动阻力的影响进行了耦合，证明气泡尺寸是影响流体阻力的重要因素30。杨银采用水平集方法对带表面张力的二维两相流问题进行数值求解，采用包含界面的小支集的光滑函数近似函数以及二阶精度的投影算法求解Navier-Stokes方程组，通过数值模拟，证明该方法可以很好地反映界面拓扑变化31。根据超声波振荡原理，结合孔板微孔的剪切作用，陈兴隆等在实验室已经生成了微米级水气分散体系，将气态CO2在超声波作用下分散于水中，由高速摄像机及体视显微镜获取图像并测量微米气泡半径的平均值约为2.5 m，远小于常规孔板喷射法生成的气泡半径（50 m），气泡的均匀度、分散性及稳定性均大幅提高11。根据长度为100 cm、直径为3.8cm的低渗透岩心的驱油实验，微米级水气分散体系在水驱结束后可继续提高采出程度超过10%，证明水气分散体系可通过气泡形变及调节渗流阻力等方式有效扩大波及体积，提高剩余油动用效果。由于水气分散体系中分散相气泡对孔隙的封堵具有累加效应，简化孔喉下的单个分散相气泡运移特征的研究将有助于加深对水气分散体系流动的认识。笔者以分散相为CO2的水气分散体系作为研究对象，其中微米气泡的半径为 2.5 m，基于 COMSOL Multiphysics开展简化孔喉中的水气分散体系运移过程数值模拟，研究多孔介质中水气分散体系运移机制，以期可以为解决狭窄空间、多孔介质等复杂结构中水气分散体系的微观渗流规律研究提供一定的参考。1计算模型1.1 水平集方法OSHER 等在 1988 年提出水平集方法，其广泛 132 油气地质与采收率2023年9月应用于多孔介质流体力学14。水平集方法将平面闭合曲线表示为点集，利用水平集函数曲面的转化来求解运动程32。这种转化虽然在一定程度上使问题复杂化，但在问题的求解上具有非常大的优势，可以使曲线的拓扑变化处理十分自然33。为了描述两相流气液界面的位置和形态变化，将气液界面转化为更高维的Level Set函数34，在模拟过程中，可以实时观察界面的变化过程。该方法是一种零等值面函数法，能有效地处理拓扑结构改变，从而在处理大多数相对复杂的多介质问题中得到很好的效果35。1.2 流动控制方程水平集函数为光滑连续函数，在界面上其值为0.5。在界面附近的过渡层，函数由0平滑地变为1。在充满水的区域 0.5。水平集函数一般与流场相结合，描述水平集函数的对流：t+u=0（1）但在模拟移动边界时，流体参数在界面处发生变化后会出现不连续性，通过定义固定的界面厚度可以解决。在这一固定界面上，参数可以平滑地变化。在COMSOL Multiphysics的CFD模块中，求解以下水平集方程：t+u=-(1-)|（2）（2）式包含稳定项和重新初始化项，其左侧定义了界面运动，右侧描述了数值稳定和重新初始化。其中为界面厚度控制参数，较理想的取值是网格大小控制参数的二分之一，即：=hmax/2（3）由于密度和黏度在界面上的突变会对仿真产生不利影响，为了解决该问题，对水平集函数进行平滑处理36，消除水与连续相、分散相的密度差、黏度差，其密度与黏度的定义为：=l+(g-l)=l+(g-l)（4）水气分散体系为不可压缩流体，流体动量方程用气液两相流不可压缩Navier-Stokes方程表示，在计算过程中需考虑表面张力。Navier-Stokes 方程为：ut+(u)u-pI+()+T=Fst（5）u=0（6）在水平集界面上，表面张力表示为：Fst=n（7）其中：=-n（8）=6|()1-|（9）2水气分散体系运移过程仿真建模根据超声波振荡原理，结合孔板微孔的剪切作用，可在实验室生成分散相为CO2的微米级水气分散体系，其中微米气泡的平均半径为 2.5 m（图1）11。当微米气泡的半径小于6 m时，气泡半径与分散性几乎呈线性关系。笔者针对实验室已生成的微米气泡半径为2.5 m的水气分散体系进行建模研究。图1微米气泡生成效果Fig.1Effect of micron bubble generation2.1 贾敏效应几何模型建立将 Navier-Stokes 方程与水平集函数耦合，在COMSOL Multiphysics 中求解，在狭窄孔喉内模拟水气分散体系两相流的流动。孔喉左侧为半径为R1=2.8 m的孔隙，长度为 15m；右侧为半径为R2=1.4 m的喉道，长度为15m，半径比为2 1，水气分散体系的分散相气泡半径为R=2.5 m。如图2所示，采用二维平面模型。第30卷第5期王富琼等.一种研究多孔介质中水气分散体系运移机制的仿真模拟新方法 133 图2水气分散体系贾敏效应模型Fig.2Jamin effect model of water-gas dispersion system2.2 边界条件设置由于简化的渗流孔喉中流体雷诺数（Re）很小，即Re 1，因此流体流态为层流。针对水气分散体系在简化孔喉中运移的仿真模型，设置边界条件如下：入口边界条件。水气分散体系以一定流速从左侧入口流入，通过狭窄孔喉，使用水平集函数为1的水平集变量的入口边界条件。出口边界条件。水气分散体系流体从固定出口流出时的压力恒定，出口压力为0MPa（此处为相对压力，参考压力为10MPa）。润湿壁边界条件。将简化的孔喉壁面加入润湿性属性37，壁面润湿程度设置接触角为，润湿壁的滑移长度设置为1m。如图3所示为贾敏效应模型边界条件。图3水气分散体系贾敏效应模型边界条件Fig.3Boundary conditions for Jamin effect model ofwater-gas dispersion system2.3 初始条件设置模拟过程中水气分散体系的分散相与连续相的物理参数，在地层条件下（50 C，10 MPa），通过实验测量得到数据（图4）38。其中，连续相（水）的密度为1 000 kg/m3，黏度为 0.562 410-3Pas，矿化度为20 000 mg/L，接触角为90，滑移长度为1m；分散相（CO2）的密度为 377 kg/m3，黏度为 2.779 0810-5Pa s，表面张力为34.9 mN/m，接触角为90，滑移长度为1m。图4界面接触角的定义Fig.4Definition of contact angle at interface2.4 网格单元划分由于分散相气泡在通过孔喉时，气液界面形态发生变化。为精确描述分散相与连续相之间的界面，在仿真过程中使用自适应网格划分。自适应网格细化有许多优势，在不对网格进行全局加密的前提下可以提高计算精度，同时降低求解自由度。在计算过程中，COMSOL Multiphysics 计算占用内存大，对网格尺寸有一定要求。应合理安排网格尺寸，使划分后的网格具有1 000100 000单元、15 00048 000个自由度，且二维模型网格质量应大于0.3。水气分散体系贾敏效应模型网格划分如图 5 所示，采用自由三角形网格进行剖分，包含9 846个三角形网格单元和5 125个网格顶点，其中最小单元质量为0.533 7，平均单元质量为0.934 6。图5水气分散体系贾敏效应模型网格划分Fig.5Meshing of Jamin effect model ofwater-gas dispersion system3模拟结果及讨论水气分散体系在多孔介质中的运移机理包括贾敏效应、聚并机理、卡断机理等，其中驱替体系在贾敏效应的作用下可提高渗流阻力，抑制驱替流体的窜流，扩大流体在油层中的波及体积。3.1 贾敏效应如图6所示，分散相CO2气泡在孔喉中运移时，其直径大于喉道直径，使得分散相气泡在运移过程中存在一定阻力并发生形变。为了克服阻力通过喉道，当分散相气泡运移至0.02 s时开始拉伸变形，随着气泡逐渐运移，在0.07 s时气泡变形程度达到最大，在0.30 s时气液界面方向发生变化，直至完全通过简化孔喉。因此，如果依靠单个水气分散体系CO2分散相气泡的贾敏效应，在地层中不具有较强的封堵能力，但当CO2分散相气泡由较大直径通过喉道且流速较大时，分散相气泡将卡堵在孔喉处，此时流动阻力增大，可以达到封堵效果。为研究水气分散体系分散相气泡通过狭窄孔喉时的压力变化，以图7中点1作为压力分析研究点。当分散相气泡通过狭窄孔喉时，由于之间存在一定半径差，分散相气泡在运移过程中必须改变自身形状以通过狭窄喉道。从图8可以看出，孔喉所受压力为动态变化过程，当分散相气泡变形程度最大时，孔喉受到压力最大。当分散相气泡大部分通过孔喉，此时压力快速降低，气液界面方向发生跳跃变化，压力出现负值。由于孔喉处压力与水气分散体系本身特征具有相关性，因此着重分析表面张力与孔喉所受压力的关系。图9表示当气液界面张力不同时，研究点1处的 134 油气地质与采收率2023年9月图6水气分散体系分散相气泡运移过程Fig.6Transport process of dispersed-phase bubblesin water-gas dispersion system图7孔喉处压力分析研究点Fig.7Pressure analysis point at pore throat图8孔喉处压力随时间变化曲线Fig.8Variation curve of pressure at pore throat with time孔喉压力随时间的变化曲线。在整体情况下，当气液界面张力增大时，压力也会随之增大；在分散相气泡运移过程中，随着表面张力的增大，分散相气泡越不易发生拉伸变形；随着分散相气泡所受压力增大，孔喉处所受压力也随之增大。将气泡在运移过程中变形的最大点作为研究点，通过曲线反映出表面张力与最大孔喉压力之间的关系（图10），表面张力与最大孔喉压力呈正比关系：p=k（10）当水气分散体系界面张力减小时，孔喉压力减小，渗流阻力减小，水气分散体系更易通过多孔介质，但不利于水气分散体系封堵孔隙，改变渗流通道，提高剩余油的动用效果。图9不同表面张力下孔喉压力随时间的变化曲线Fig.9Variation curve of pore throat pressure with timeunder different surface tension图10表面张力与最大孔喉压力关系曲线Fig.10Relationship between pressure at maximum porethroat and surface tension3.2 聚并机理水气分散体系在多孔介质运移过程中，当分散相气泡之间的间距小于一定距离时，两者之间会发生聚并现象。如图11所示，针对该现象建立二维轴对称聚并模型并进行模拟。简 化 渗 流 通 道 半 径 为3.8 m，孔 喉 长 度 为图11二维轴对称聚并模型Fig.11Two-dimensional axisymmetric coalescence model第30卷第5期王富琼等.一种研究多孔介质中水气分散体系运移机制的仿真模拟新方法 135 30m，水为连续相，CO2为分散相，2个分散相气泡初始半径均为2.5m，气液界面张力为34.9 mN/m，入口流速为50 mm/s，出口端压力为0 Pa。通过自由三角形网格对聚并模型进行网格剖分（图12）。图12水气分散体系分散相气泡聚并模型网格划分Fig.12Meshing of coalescence model of dispersed-phasebubbles in water-gas dispersion system由分散相气泡聚并过程模拟结果（图 13）可以看出，2个分散相气泡在0.001 s时逐渐接近并发生接触，其间会形成一薄液膜，相撞的分散相气泡逐渐靠近，使两者间的连续相逐渐向外排出，当连续相液体被排出到临界点时，分散相气泡聚集在一起，界面开始融合，分散相气泡界面逐渐向最小表面能趋势进化。图14为分散相气泡接触后的流场分布，流线沿计算区域的两极朝赤道方向分布，同时2个分散相气泡中间的液体由中线流向外缘，实现排液过程，达到聚并。而在实际地层温度与压力条件下，水气分散体系可压缩，其运移速度远小于数值模拟过程，且水气分散体系所处的温度条件是不断变化的，其黏度与表面张力也在不断变化。3.3 卡断机理水气分散体系在简化渗流通道中流动，通过狭窄孔喉时易发生卡断现象。建立简化孔喉下的二维卡断模型（图15），并对水气分散体系卡断机理进行分析。模拟渗流通道宽度为5.6m，长度为30m，分散相气泡初始半径均为 2.5m，其中水为连续相，CO2为分散相，气液界面张力为34.9 mN/m，入口流速为20 mm/s，出口端压力为0 Pa。如图16所示，采用自由三角网格对水气分散体系卡断模型进行网格划分。如图17所示，由于模拟的渗流通道尺寸存在较大差异，水气分散体系的分散相在通过狭窄喉道前，所受压力逐渐增大，分散相气泡被压缩。当分散相气泡运移至接近0.03 s时，喉道阻力达到一定程度，此时分散相气泡前端通过狭窄喉道，并且产生局部卡断，形成小气泡。由于此时小气泡距离较近，达到了聚并临界距离，从而产生聚并。随后通过的小气泡在喉道右侧与分散相气泡发生聚并形成大气泡，随着液体流出。4结论基于水平集方法与气液两相流模型，通过多物图13水气分散体系分散相气泡聚并过程模拟结果Fig.13Simulation results of coalescence process ofdispersed-phase bubbles in water-gas dispersion system图14水气分散体系分散相聚并过程及对应的流场分布Fig.14Dispersed-phase coalescence process and corresponding flow field distribution of water-gas dispersion system图15二维卡断模型Fig.15Two-dimensional division model理场耦合方法，对水气分散体系在多孔介质简化模型中的贾敏效应、聚并机理、卡断机理进行一系列数 136 油气地质与采收率2023年9月图16水气分散体系卡断模型网格划分Fig.16Meshing of division model of water-gasdispersion system图17水气分散体系卡断过程模拟结果Fig.17Simulation results of division process ofwater-gas dispersion system值模拟，得出以下结论：针对超声波振荡法产生的分散相为 CO2的微米级水气分散体系建立二维模型，对不同孔喉形状下的水气分散体系运移过程进行模拟，分析水气分散体系在孔喉运移过程中的贾敏效应，探讨水气分散体系运移过程中孔喉压力的影响因素。研究发现在孔喉中，分散相气泡的变形程度愈大，孔喉中的压力愈大，渗流阻力越大。建立水气分散体系在孔喉中的聚并、卡断模型，研究分散相气泡聚并、卡断现象的动态特征，结果表明，距离较近的水气分散体系的分散相气泡在运移过程中会逐渐聚集在一起，界面发生融合，分散相气泡界面逐渐向最小表面能趋势进化，其流场表现为流线沿计算区域的两极朝赤道方向分布，分散相气泡中间的液体由中线流向外缘，实现排液过程，达到聚并；分散相气泡在通过狭窄孔喉时会发生卡断现象产生小气泡，当到达聚并临界距离时，会再次产生聚并随着液体流出。通过数值模拟可以直接观察分散相气泡在孔喉运移过程中的变化，明确水平集方法在计算过程中的优势。建立的数值模拟方法在研究微细气泡运移方面具有一定的应用前景，可以为水气分散体系通过气泡形变及调节渗流阻力、有效扩大波及体积、提高剩余油动用效果的数值模拟研究提供一种新的思路与方法。符号解释F在全流场定义每个欧拉网格单元内某种流体的体积比函数；Fst气液界面张力，N/m；hmax网格大小控制参数，m；I单位矩阵；k曲率，m-1；n指向分散相的界面单位法线界向量；p最大孔喉压力，Pa；Re渗流孔喉中的流体雷诺数；t时间，s；u流体速度向量，m/s；u流体速度，m/s；界面初始化参数（等于入口速度），m/s；流体界面处非零的狄拉克函数；界面厚度控制参数，m；界面接触角，（）；比例系数；动力黏度，Pas；gCO2的黏度，Pas；l水的黏度，Pas；密度，kg/m3；gCO2的密度，kg/m3；l水的密度，kg/m3；表面张力系数；水平集函数；t水平集函数对t的时间偏导数；拉普拉斯算子。参考文献 1 熊春明，唐孝芬.国内外堵水调剖技术最新进展及发展趋势J.石油勘探与开发，2007，34（1）：83-88.XIONG Chunming，TANG Xiaofen.Technologies of water shut-off and profile control：An overview J.Petroleum Explorationand Development，2007，34（1）：83-88.2 冯其红，王波，王相，等.高含水油藏细分注水层段组合优选方法研究J.西南石油大学学报：自然科学版，2016，38（2）：103-108.FENG Qihong，WANG Bo，WANG Xiang，et al.Study on lay第30卷第5期王富琼等.一种研究多孔介质中水气分散体系运移机制的仿真模拟新方法 137 er combination optimization method of subdivision water injection in high water-cut reservoirs J.Journal of Southwest Petroleum University：Science and Technology Edition，2016，38（2）：103-108.3 胡文瑞，魏漪，鲍敬伟.中国低渗透油气藏开发理论与技术进展J.石油勘探与开发，2018，45（4）：646-656.HU Wenrui，WEI Yi，BAO Jingwei.Development of the theoryand technology for low permeability reservoirs in China J.Petroleum Exploration and Development，2018，45（4）：646-656.4 张世明，杨勇.特高含水水驱油藏流场边界动态追踪模型研究J.油气地质与采收率，2021，28（2）：91-99.ZHANG Shiming，YANG Yong.Study on the dynamic trackingmodel of flow field boundaries in water drive reservoirs at ultra-high water cut stage J.Petroleum Geology and Recovery Efficiency，2021，28（2）：91-99.5 夏惠芬，王立辉，韩培慧，等.薄膜状剩余油动用条件研究J.特种油气藏，2021，28（3）：106-111.XIA Huifen，WANG Lihui，HAN Peihui，et al.Research onproduction conditions of film remaining oil J.Special Oil&Gas Reservoirs，2021，28（3）：106-111.6 张庆龙.多层非均质层间窜流定量表征及开发对策J.特种油气藏，2021，28（1）：136-141.ZHANG Qinglong.Quantitative characterization and development countermeasures of multi-layer heterogeneous interlayerchanneling flow J.Special Oil&Gas Reservoirs，2021，28（1）：136-141.7 俞宏伟，李实，陈兴隆，等.微分散驱油体系流动性能实验评价J.科学技术与工程，2013，13（22）：6 426-6 429.YU Hongwei，LI Shi，CHEN Xinglong，et al.The evaluation offlow properties for differential oil flooding system J.ScienceTechnology and Engineering，2013，13（22）：6 426-6 429.8 俞宏伟，马德胜，李实，等.水气分散体系提高采收率实验研究J.科学技术与工程，2013，13（33）：9 946-9 949.YU Hongwei，MA Desheng，LI Shi，et al.Experimental studyon enhanced oil recovery by water-gas dispersion flooding system J.Science Technology and Engineering，2013，13（33）：9 946-9 949.9 王琼，王金光.自生泡沫复合驱油体系性能研究J.内蒙古石油化工，2017，43（9）：115-117.WANG Qiong，WANG Jinguang.Research on the performanceof authigenic foam composite flooding system J.Inner Mongolia Petrochemical，2017，43（9）：115-117.10孙哲，吴行才，康晓东，等.连续相与分散相驱油体系驱油机理及其性能对比J.石油勘探与开发，2019，46（1）：116-124.SUN Zhe，WU Xingcai，KANG Xiaodong，et al.Comparisonof oil displacement mechanisms and performances between continuous and dispersed phase flooding agents J.Petroleum Exploration and Development，2019，46（1）：116-124.11陈兴隆，伍家忠，韩海水，等.水气分散体系的超声波振荡生成方法及驱油实验研究J.天然气与石油，2021，39（4）：48-52.CHEN Xinglong，WU Jiazhong，HAN Haishui，et al.Formation of gas bubble-water dispersion system using ultrasonic vibration method and research on its impact on ultimate oil recovery J.Natural Gas and Oil，2021，39（4）：48-52.12端祥刚.可调阻力水气分散体系提高采收率方法与理论研究D.北京：中国石油大学（北京），2015.DUAN Xianggang.The method and theoretical research on EORof water-gas dispersion system with adjustable resistance D.Beijing：China University of Petroleum（Beijing），2015.13LIU N，CHEN X，JU B，et al.Microbubbles generation by anorifice spraying method in a water-gas dispersion flooding system for enhanced oil recovery J.Journal of Petroleum Scienceand Engineering，2021，198：108196.14OSHER S，SETHIAN J A.Fronts propagating with curvature-dependent speed：Algorithms based on Hamilton-Jacobi formulations J.Journal of Computational Physics，1988，79（1）：12-49.15OLSSON E，KREISS G，ZAHEDI S.A conservative level setmethod for two phase flow II J.Journal of Computational Physics，2007，225（1）：785-807.16RAEINI A Q，BLUNT M J，BIJELJIC B.Modelling two-phaseflow in porous media at the pore scale using the volume-of-fluidmethod J.Journal of Computational Physics，2012，231（17）：5 653-5 668.17OSHER S,FEDKIW R.Level set methods and dynamic implicitsurfacesM.New York:Springer,2003:23.18SUN Y，BECKERMANN C.Sharp interface tracking using th