基于随机变径毛细管束模型的渗吸-驱替机理分析.pdf

第30卷第5期油气地质与采收率Vol.30，No.52023年9月Petroleum Geology and Recovery EfficiencySep.2023引用格式：高涛，王香增，王锰，等.基于随机变径毛细管束模型的渗吸-驱替机理分析J.油气地质与采收率，2023，30（5）：100-109.GAO Tao，WANG Xiangzeng，WANG Meng，et al.Analysis of imbibition-displacement mechanism based on random variable-diameter capillary bundle modelJ.Petroleum Geology and Recovery Efficiency，2023，30（5）：100-109.基于随机变径毛细管束模型的渗吸-驱替机理分析高 涛1，2，王香增2，3，王 锰1，党海龙1，2，刘芳娜1（1.陕西延长石油（集团）有限责任公司 研究院，陕西 西安 710065；2.陕西省特低渗透油气田勘探与开发工程技术研究中心，陕西 西安 710065；3.陕西延长石油（集团）有限责任公司，陕西 西安 710065）摘要：针对等直径毛细管束模型的不足，结合分段法推导得到考虑压差及固-液壁面作用的液-液系统下变径毛细管多段式渗吸理论公式，研究不同毛细管几何结构和流体性质的变化对变径毛细管渗吸-驱替的影响，并通过Python编制随机变径毛细管束模型，建立基于该模型的采出程度及含水率计算方法。结果表明：压差的增大导致变径毛细管段间渗吸速度差异减小，变径毛细管渗吸-驱替具有方向性；基于随机变径毛细管束模型计算的采出程度及含水率与实验真实岩心驱替数据符合较好，且渗吸作用对总采出程度的贡献量达19.97%。关键词：变径毛细管束；渗吸-驱替模型；温和注水；特低渗透油藏；致密储层文章编号：1009-9603（2023）05-0100-10DOI：10.13673/j.pgre.202210027中图分类号：TE357.6文献标识码：AAnalysis of imbibition-displacement mechanism basedon random variable-diameter capillary bundle modelGAO Tao1，2，WANG Xiangzeng2，3，WANG Meng1，DANG Hailong1，2，LIU Fangna1（1.Research Institute of Shaanxi Yanchang Petroleum（Group）Co.，Ltd.，Xi an City，Shaanxi Province，710065，China；2.Shaanxi Extra Low Permeability Oil and Gas Field Exploration and Development Engineering TechnologyResearch Center，Xi an City，Shaanxi Province，710065，China；3.Shaanxi Yanchang Petroleum（Group）Co.，Ltd.，Xi an City，Shaanxi Province，710065，China）Abstract：In view of the shortcomings of the equal-diameter capillary bundle model，this paper theoretically derives the multi-segment imbibition formula of the variable-diameter capillary under the liquid-liquid system considering pressure difference and solid-liquid wall effect by using the sectioning method.In addition，the paper studies the influence of changes in different capillary geometry and fluid properties on the imbibition-displacement of the variable-diameter capillaries and develops a random variable-diameter capillary bundle model through Python，so as to deduce calculation methods of the recovery and water cut based on this model.The results show that the increase in pressure differences leads to the decrease in the differences of the imbibition velocities betweenvariable-diameter capillary segments，and the imbibition-displacement of the variable-diameter capillaries is directional.The recovery and water cut calculated based on the random variable-diameter capillary bundle model are in good agreement with the real coredisplacement data from experiments，and the contribution of the imbibition effects to the total recovery reaches 19.97%.Key words：variable-diameter capillary bundle；imbibition-displacement model；moderate water injection；ultra-low permeabilityreservoir；tight reservoir通过理论研究和矿场实践认为，利用基质与裂缝之间的油水渗吸置换作用可较好地采出基质内常规注水无法动用的原油，基于该开发理念建立的“控速增渗”温和注水技术可有效改善低渗透致密油藏收稿日期：2022-10-28。作者简介：高涛（1988），男，陕西延安人，高级工程师，硕士，从事特低渗透致密油气开采理论与工程技术攻关工作。E-mail：。通信作者：王香增（1968），男，河南滑县人，教授级高级工程师，博士。E-mail：。基金项目：陕西省科技统筹创新项目“延长难采储量有效动用开发技术研究”（2016KTCL01-12），陕西省创新能力支撑计划项目“非常规油藏低碳环保提高采收率技术创新团队”（2023-CX-TD-31）。第30卷第5期高 涛等.基于随机变径毛细管束模型的渗吸-驱替机理分析 101 水驱开发效果。基质与裂缝之间的油水渗吸置换作用主要利用毛细管模型表征，进而对致密油藏注水开发方式进行优化。微观渗吸理论上，等直径毛细管1-5、非等径毛细管6、双毛细管7-8等都具有一定的研究基础。对于等直径毛细管，LUCAS等最早分析了简单毛细管中影响渗吸的动力学因素，提出了经典LW方程用于计算毛细管内液面上升高度9-11，但由于该方程的假设前提为简单毛细管牛顿流体层流，从而限制了其在复杂微观多孔介质中的应用。部分学者基于LW方程又开展了大量的研究工作。BENAVENTE等通过在岩石微观结构中引入各种校正因子，将理论研究扩展到非均质多孔介质12；TIAN等通过引入分子动力学理论研究了动态接触角对双液相系统下毛细管自发渗吸的影响13-14；WANG等分析了分子相互作用及入口效应对纳米多孔介质渗吸的影响15。对于非等径毛细管，ERICKSON等通过有限元方法模拟了气-液系统在具有变径结构毛细管中的自发渗吸过程16；YOUNG推导了非等径毛细管气-液渗吸公式17；REYSSAT等通过实验和理论的方法研究了2种扩径方式（锲型和幂律型）下的气-液渗吸公式18；BERTHIER等研究了在半开放式毛细管中气-液渗吸规律19。对于双毛细管，NABIZADEH等研究了双毛细管系统中液-液非混相渗吸规律7。综上可知，大部分研究主要分析简单非等径毛细管在气-液系统下的渗吸规律，而对液-液系统复杂毛细管的渗吸规律研究不够深入。为此，笔者针对液-液系统下非等径毛细管的渗流特点建立了渗吸-驱替双重作用模型，并考虑渗吸-驱替过程中压差及固-液壁面作用对渗吸-驱替的影响，结合真实岩心孔径分布曲线建立了表征真实岩心渗流特征的随机变径毛细管束模型，计算该模型中不同毛细管束的含水率和采出程度，并与室内水驱油实验测试结果进行对比，从而验证了该模型的可靠性。研究成果丰富了延长油田“控速增渗”温和注水技术，通过发挥渗吸-驱替协同作用对低渗透致密储层提高采收率具有一定指导意义。1基础模型1.1 模型结构设计2个具有不规则横截面的圆柱变径毛细管（图1），其几何形状复杂程度不同。毛细管完全饱和非润湿相并处于静态，入口和出口边界分别与润湿相和非润湿相相连。变径毛细管模型假设条件如下：流体为不可压缩流体，两相流体不存在混溶，流体流动属于层图1变径毛细管模型示意Fig.1Variable-diameter capillary model流。不考虑相界面厚度对位移距离和渗吸速度的影响20。由于研究的毛细管直径为微米级别，模型推导忽略重力影响。1.2 模型推导1.2.1 单段变径毛细管模型在渗吸过程中驱替相流体主要受液柱两端压差、毛细管力、管壁施加的黏性力、重力以及惯性力的影响，关系式可表示为：Fs+Fc=Fv+Fg+Fi（1）对于流动方向为水平方向的扩径毛细管模型（图1a），忽略重力与惯性力在渗吸方向上的影响，则（1）式简化为：Fs+Fc=Fv（2）假设驱替相流体两端压差等于毛细管两端压差，则有：Fs=r2xps=r2x(psi-pso)（3）驱替相流体前缘到达任意毛细管任意截面时的毛细管力可表示为：Fc=r2xpc=2rxcos（4）驱替相流体受到的壁面摩擦力等于润湿相流体摩擦力与非润湿相流体摩擦力之和，运用牛顿黏性定理可得：Fv=0lx-2rxwdvdrdl-lxl2rxnwdvdrdl（5）将（3）（5）式代入（2）式可得驱替相流体任意截面处流速分布，其表达式为：V(r)=ps+pc41w0lx1r2x-r2dlx+nwlxl1r2x-r2dlx（6）102 油气地质与采收率2023年9月根据（6）式可得驱替相渗吸时的体积流量，其表达式为：Q=(ps+pc)81w0lx1r4xdl+nwlxl1r4xdl（7）将（3）和（4）式代入（7）式，渗吸速度可表示为：v=(psi-pso)rx+2cos8r3x1w0lx1r4xdl+nwlxl1r4xdl（8）将v=dlx/dt代入（8）式，渗吸时间的表达式为：t=0lx8r3x(psi-pso)rx+2cos()w0lx1r4xdlx+nwlxl1r4xdlxdlx（9）1.2.2 多段式变径毛细管模型对于多段式变径毛细管模型（图1b），渗吸时间可用分段法形式表示19-23，首先引入中间积分变量：Ci=0li1r4idlx,i（10）将（8）和（9）式分别进行变换，其表达式依次为：vk,x=(psi-pso)rk,x+2cos8r3k,x1ek,x（11）t=0lx8r3k,x(psi-pso)rk,x+2cosek,xdlx（12）其中：ek,x=w()1k-1Ci+0lk,x1r4k,xdlk,x+nw()k+1nCi+lk,xlk1r4k,xdlk,x（13）1.2.3 考虑固-液壁面作用的多段式变径毛细管模型多段式变径毛细管模型中驱替相流体黏度与被驱替相流体黏度均假设为定值，这一假设在宏观流动中较为合理，但在微纳米尺度中却存在一定问题24-25。由于黏度是流体流动过程中流体分子所受阻力的度量，其本质是对分子间吸引力的表征，在微观流动中由于通道尺寸较小，部分流体分子受管壁固体分子范德华力作用的影响从而产生附加黏性26-35，其计算公式为：=0+bAsAw-Awr（14）从（14）式可以看出附加黏度与流动通道尺寸相关，流动通道尺寸越小，附加的黏性效应越明显，在微纳米尺度的流动中，这一附加黏度对流体流动的影响不可忽略。在多段式毛细管渗吸速度求解过程中，由于各段毛细管存在尺寸差异，则附加黏度也各不相同，当计算段数较多或尺寸差异较明显时，累加计算的误差将导致计算结果出现明显偏差，所以在计算渗吸速度时，需对每段的黏度进行修正，修正后的结构参数表示为：ek,x=()w0+bwAsAw-Awr()1k-1Ci+0lk,x1r4k,xdlk,x+()nw0+bnwAsAnw-Anwr()k+1nCi+lk,xlk1r4k,xdlk,x（15）将（15）式代入（11）式即可计算考虑附加黏度后的渗吸速度。1.3 模型验证选取模型的几何参数和流体参数与 ERICKSON等对简单扩径-缩径毛细管和缩径-扩径毛细管的渗吸研究中一致16。所有毛细管均有2个级别的半径，分别为 25和 50 m，每根毛细管内各段长度l1，l12，l2，l23和l3分别为20，3，40，3和20 mm。流体参数设置为：表面张力为 0.03 N/m，润湿相黏度为 1mPas，非润湿相黏度为0.018 mPas，接触角为30。采用文献中相关的复杂毛细管渗吸数据16，19验证新模型的可靠性。将新模型计算的扩径-缩径毛细管和缩径-扩径毛细管渗吸速度与ERICKSON数值解16和BERTHIER解析解19的预测结果进行对比（图2），结果显示新模型在2种复杂毛细管中的预测结果与ERICKSON数值解拟合效果更好。但从图2b可以看出经过2次变径后，扩径-缩径毛细管模型变径后半段的BERTHIER解析解的拟合效果偏差较大，对于更为复杂的变径毛细管，将会产生较大的计算误差，从而影响预测结果的准确性。相较于BERTHIER等提出的解析模型，新模型在变径段的拟合效果更好，更符合ERICKSON数值模拟的结果。2渗流特征的影响因素2.1 压差毛细管压力与压差均是渗吸-驱替的动力，但由于圆管中的毛细管压力受毛细管半径的影响，所以结构相同但尺寸不同的变径毛细管对压差与毛细管压力的敏感性不同。为了研究不同压差对变径毛细管渗吸-驱替的影响，选取图2所示2种毛细管结构，分别将其在直径方向上按等比例缩放得到6种不同尺寸变径毛细管（图3a3b为缩径-扩径毛细管，图3c3f为扩径-缩径毛细管），各段长度与图2所示相第30卷第5期高 涛等.基于随机变径毛细管束模型的渗吸-驱替机理分析 103 同，毛细管压差设置为010 kPa。由图3可以看出，随着压差的增大，渗吸时间曲线各段斜率趋于一致，表明各段渗吸速度趋于相同，故毛细管两端压差的增大有助于缩小变径毛细管段间渗吸速度的差距。且对于相同结构不同尺寸的变径毛细管，大尺寸结构的渗吸速度趋同现象更为明显。变径毛细管段间渗吸速度差异主要受毛细管压力与流动阻力的控制，若将（6）式拆分成压力项与毛细管压力项，则可以看出压差对渗吸速度的影响受毛细管半径的二次方控制，与毛细管半径呈线性关系，可简化为：v pr2+pcr（16）对于某一确定的变径毛细管结构，毛细管内的小半径段受压力的影响比大半径段更强，所以随着图22种复杂毛细管不同计算模型下渗吸速度与渗吸距离的关系Fig.2Relationships between imbibition velocities and imbibition distances under differentcalculation models of two complex capillaries图3不同尺寸变径毛细管渗吸时间随压差的变化Fig.3Imbibition time with pressure differences of different variable-diameter capillaries 104 油气地质与采收率2023年9月压差的增大，段间渗吸速度差距逐渐消失。而在毛细管整体尺寸变化时，压力的影响将以二次方的形式扩大或缩小，这将导致当毛细管半径足够大时，毛细管压力对渗吸速度的影响减弱，压力的作用占据主导。2.2 毛细管结构为探究半径分布对渗吸时间的影响，利用固定段长的方法研究三段式半径结构对渗吸-驱替作用的影响。设计6种三段式变径毛细管结构以分析气-液系统和液-液系统条件下渗吸时间的变化。6种三段式变径毛细管内各段长度l1，l12，l2，l23和l3均分别为30，10，30，10和 30 mm。每根变径毛细管均由 10，20，30 m三种级别的半径组成，其排列顺序如图4所示。流体参数设置为：表面张力为30 mN/m，接触角为 30，润湿相黏度为 1 mPas，固相和液相 Hamaker常数均为4.210-20J。图4不同半径结构三段式变径毛细管结构示意Fig.4Three-segment variable-diameter capillarystructures with different radii从图5可以看出，对于包含相同尺寸及数量直管段的变径毛细管结构，半径排列顺序对最终渗吸时间产生影响。6种不同半径结构的三段式变径毛细管在气-液系统下的渗吸时间曲线趋于同一点，曲线各段斜率不随排列顺序变化而变化，曲线终点对应的最终渗吸时间基本一致。而液-液系统下的最终渗吸时间具有明显差异，且曲线各段斜率因排列顺序不同而变化显著。图5a中气-液系统在最终渗吸时间上的差异主要来源于变径段结构的差异。而在图5b的液-液系统中，流动阻力受润湿相黏度、非润湿相黏度及毛细管半径共同影响，且非润湿相黏度大于润湿相黏度，所以非润湿相的流动阻力贡献量更大，易随出口段毛细管尺寸的减小而增大。结合图4可以看出，图5b中变径毛细管所需的渗吸时间大于变径毛细管，这一现象说明变径毛细管在液-液系统下渗吸-驱替会表现出一定的方向性。因此，在变径毛细管中的渗吸过程不仅要考虑孔隙尺寸、段长的分布，还需考虑各段之间的位置关系以及流体流动方向。?图5不同半径结构三段式变径毛细管在气-液系统及液-液系统下渗吸时间的变化Fig.5Imbibition time of three-segment variable-diametercapillary structures with different radii undergas-liquid and liquid-liquid systems3随机变径毛细管束模型3.1 模型建立储层是具有复杂孔隙结构的多孔介质，以往的第30卷第5期高 涛等.基于随机变径毛细管束模型的渗吸-驱替机理分析 105 研究中通常将天然岩石中不规则孔隙简化为等直径毛细管组成的理想流动模型，并基于此模型建立采出程度及含水率的计算公式，但该类模型无法描述重要的变径结构信息。而在变径毛细管模型中，对于毛细管半径分布广、段数多、各段长度不定的变径毛细管，无法对所有可能出现的毛细管结构进行单独分析。为了反映孔隙结构的弯曲程度，需为单根毛细管引入迂曲度的定义，其表达式为：=ltL（17）在多种变量作用下，单独分析某一变径结构意义不大。为了研究复杂结构变径毛细管的渗流特征，利用Python编制变径毛细管生成程序，随机批量生成具有一定半径分布特征及段数分布特征的毛细管，并通过对直管段添加随机长度的方式，生成固定迂曲度条件下任意结构的弯曲随机变径毛细管（图6），将获得的毛细管结构信息代入模型计算，即可得到任意结构弯曲随机变径毛细管的渗吸时间。图6随机变径毛细管示意Fig.6Random variable-diameter capillaries基于上述毛细管生成程序，进一步构建以真实岩心总孔隙度、孔隙分布实验数据及迂曲度等参数为约束条件的随机变径毛细管束模型（图7），与经图7不同孔隙度和迂曲度下随机变径毛细管束模型Fig.7Random variable-diameter capillary bundle modelunder different porosity and tortuosity典毛细管束模型相比，其更真实地反映孔隙微观结构。通过建立基于随机变径毛细管束模型的采出程度及含水率计算方法分析真实岩心尺度下渗吸-驱替过程。3.2 采出程度及含水率表征对于上述模型中的单根变径毛细管，某一时间t下，毛细管内油水相界面位于该管第k段中的lkt处，且该处毛细管半径为rkt，则在该时间下第i根变径毛 106 油气地质与采收率2023年9月细管的水相总体积表达式为：Vwi=j=1kt-1r2jlj+0lktr2ktdlkt（18）由于水相总体积等于油相采出体积，则求解所有变径毛细管在t时间下的水相体积并求和，即可得到任意时间下变径毛细管束的油相采出体积。油相采出体积与毛细管束总体积的比值即为采出程度，表达式为：R=i=1nj=1kt-1r2jlj+0lktr2ktdlkti=1nj=1kr2jlj（19）同理，出口端含水率的表达式为：fw=i=1nBir2endi,tvendi,ti=1nr2endi,tvendi,t（20）其中：Bi=1tendi t（21）根据连续性定理（20）式可变为：fw=i=1nBir2i,ktvi,kti=1nr2i,ktvi,kt（22）Bi为第i根毛细管末端的状态参数，是由第i根毛细管出口端见水时间tendi与当前渗吸时间相对大小确定，当Bi=1时表明第i根毛细管末端产水，当Bi=0时表明第i根毛细管末端产油。4实例应用选取延长油田西部油区主力开发层位长6油层岩心 4 块，其平均孔隙度为 10.9%，平均渗透率为1.53 mD，孔隙半径分布如图8所示。设置随机变径毛细管束模型的段数为 50100，岩心长度为 0.05m，迂曲度为12，毛细管个数为10 000。将图8中孔隙半径分布作为毛细管束1的孔隙半径分布，将得到的毛细管束代入模型中进行计算。同时，通过对原始孔隙半径分布添加微孔隙的方式模拟部分死孔隙，得到毛细管束2孔隙半径并代入模型计算。由随机变径毛细管束模型计算的不同毛细管束采出程度及含水率变化（图9）可以看出，毛细管束1计算得到的采出程度在高含水阶段与实验值偏差较大，由于毛细管束1采用的半径分布中各级孔隙半径较大，所以毛细管束1中的所有孔隙均图8延长油田西部油区长6油层岩心孔隙半径分布Fig.8Pore radius distribution of Chang 6 core inwestern Yanchang Oilfield图9不同毛细管束采出程度及含水率的变化Fig.9Change of water content and recovery degree ofdifferent capillary bundles为连通孔隙，其计算得到的最终采出程度基本达100%。而毛细管束2由于添加了微孔隙，计算结果与实验值更加符合，模型计算最终采出程度为65.6%，实验值为 61.86%。同时含水率的变化不是平滑地上升，而是存在台阶式变化，这与最新的水驱油实验测试结果一致。而传统水驱油实验测试由于计量间隔时间长，数据点是一段时间累积效应并不是瞬时值，故本模型计算的含水率变化特征更能真实地反映储层含水率变化特征。由渗吸-驱替作用在总采出程度中的贡献量（图10）可以看出驱替作用图10渗吸-驱替作用在总采出程度中的贡献量Fig.10Contribution of imbibition-displacementeffect to total recovery第30卷第5期高 涛等.基于随机变径毛细管束模型的渗吸-驱替机理分析 107 采出程度为52.5%，渗吸作用采出程度为13.1%，总采出程度为65.6%。在毛细管束2条件下，渗吸作用对总采出程度的影响可达19.97%，即对于延长油田长 6 油 层 渗 吸 作 用 对 总 采 出 程 度 的 贡 献 量 达19.97%。但渗吸作用对总采出程度的贡献是由生产压差、孔隙半径、润湿性的强弱共同决定，孔隙半径越小渗吸作用贡献量越大，同时通过改变储层的亲水性以及合理控制生产压差可以进一步提高渗吸作用的贡献。5结论通过多段式变径毛细管模型，分析了压差与毛细管结构对渗流特征的影响，并建立了符合真实岩心孔隙结构的随机变径毛细管束模型。对于某一确定的变径毛细管结构，毛细管内大半径段的压力敏感性比小半径段更强，当压差高于某一临界值时，变径毛细管的段间渗吸速度差异不明显，且临界压差与毛细管尺寸密切相关。对于包含相同尺寸及数量直管段的变径毛细管结构，半径分布的顺序会对最终渗吸时间产生影响，其中气-液系统受影响较小，而液-液系统的渗吸-驱替过程对变径毛细管结构显示出一定的方向性。采用真实岩心孔隙半径分布数据，构建了随机变径毛细管束模型，计算不同毛细管束的采出程度及含水率，结果表明渗吸作用对总采出程度的贡献量达19.97%。符号解释Anw非润湿相流体Hamaker常数，J；As管壁Hamaker常数，J；Aw润湿相流体Hamaker常数，J；b管壁与流体分子作用的黏度系数，Pasm；bw管壁与润湿相流体分子作用的黏度系数，Pasm；bnw管壁与非润湿相流体分子作用的黏度系数，Pasm；Bi第i根毛细管末端的状态参数；Ci中间积分变量，m-3；ek，x多段式变径毛细管内第k段距出口x处的结构参数，Pas/m2；fw含水率，%；Fc毛细管力，N；Fi惯性力，N；Fg重力，N；Fs压力，N；Fv黏性力，N；i毛细管编号；j直管段编号；k任意多段式变径毛细管直管段总数；ktt时间下含润湿相的直管段总数；l毛细管长度，m；l1，l2，ln多段式变径毛细管内第1，2，n段直管段长度，m；l12，l23多段式变径毛细管内第 1和第 2，第 2和第 3直管段间变径段长度，m；li第i段毛细管长度，m；lj第j段直管段长度，m；lk多段式变径毛细管内第k段直管段长度，m；lkt多段式变径毛细管内t时间下相界面所在第k段长度，m；lk，x多段式变径毛细管内第k段距入口x处直管段长度，m；lt毛细管路径长度，m；lx毛细管内距入口x处的长度，m；lx，i毛细管第i段内距离该段入口x处的长度，m；L毛细管直线长度，m；n毛细管根数；pc毛细管压力，Pa；psi进口压力，Pa；pso出口压力，Pa；Q体积流量，m3/s；r半径，m；r0单段扩径毛细管入口半径，m；r1，r2，r3，rn多段式变径毛细管内第1，2，3，n段直管段半径，m；rendi，tt时间下第i根毛细管出口半径，m；ri第i段毛细管半径，m；ri，kt第i根毛细管t时间下相界面所在第k段半径，m；rj第j段毛细管半径，m；rk多段式变径毛细管内第k段半径，m；rkt多段式变径毛细管内t时间下相界面所在第k段半径，m；rk，x多段式变径毛细管内第k段距入口x处的半径，m；rl长度为l的单段扩径毛细管入口半径，m；rlx长度为l的单段扩径毛细管距入口x处的半径，m；rn多段式变径毛细管出口半径，m；rx变径毛细管内距入口x处的半径，m；R采出程度，%；t渗吸时间，s；tendi第i根毛细管出口见水时间，s；v渗吸速度，m/s；vendi,tt时间下第i根毛细管出口渗吸速度，m/s；vi,kt第i根毛细管t时间下第k段渗吸速度，m/s；vk,x多段式变径毛细管内第k段距入口x处的渗吸速度，m/s；108 油气地质与采收率2023年9月V（r）流速分布函数，m/s；Vw毛细管束的水相总体积，m3；Vwi第i根毛细管的水相体积，m3；x相界面距毛细管入口距离，m；p压差，Pa；ps毛细管两端压差，Pa；界面张力，N/m；接触角，（）；总黏度，Pas；0初始黏度，Pas；nw非润湿相黏度，Pas；nw0非润湿相初始黏度，Pas；w润湿相黏度，Pas；w0润湿相初始黏度，Pas；迂曲度。参考文献 1 XIAO Junfeng，LIU Xiang，LUO Youming，et al.Oscillationsof free surface at the edge of short capillary tubes J.Colloidsand Surfaces A：Physicochemical and Engineering Aspects，2020，591：124572.2 黎明，廖晶，王肃，等.鄂尔多斯盆地超低渗透油藏渗吸特征及其影响因素以渭北油田三叠系延长组三段储层为例J.石油实验地质，2022，44（6）：971-980.LI Ming，LIAO Jing，WANG Su，et al.Imbibition characteristics and influencing factors of reservoirs with ultra-low permeability of Ordos Basina case study of third member of TriassicYanchangFormationinWeibeiOilFieldJ.PetroleumGeology&Experiment，2022，44（6）：971-980.3 XIAO Junfeng，LUO Youming，NIU Muyuan，et al.Study ofimbibition in various geometries using phase field method J.Capillarity，2019，2（4）：57-65.4 PETER L L W，GRGOIRE D，JAMES C B，et al.Capillarydisplacement of viscous liquids J.Langmuir，2016，32（13）：3 186-3 190.5 GORCE J B，HEWITT I J，VELLA D.Capillary imbibition intoconverging tubes：Beating Washburn s law and the optimal imbibition of liquids J.Langmuir，2016，32（6）：1 560-1 567.6 李颖，李茂茂，李海涛，等.水相渗吸对页岩储层的物化作用机理研究J.油气藏评价与开发，2023，13（1）：64-73.LI Ying，LI Maomao，LI Haitao，et al.Physicochemical mechanism of water phase imbibition in shale reservoirs J.ReservoirEvaluation and Development，2023，13（1）：64-73.7 NABIZADEH A，ADIBIFARD M，HASSANZADEH H，et al.Computational fluid dynamics to analyze the effects of initialwetting film and triple contact line on the efficiency of immiscible two-phase flow in a pore doublet model J.Journal of Molecular Liquids，2019，273：248-258.8 CHATZIS I，DULLIEN F.Dynamic immiscible displacementmechanisms in pore doublets：Theory versus experiment J.Journal of Colloid&Interface Science，1983，91（1）：199-222.9 LUCAS R.Ueber das Zeitgesetz des kapillaren aufstiegs vonFlssigkeiten J.Kolloid-Zeitschrift，1918，23：15-22.10CAI Jianchao，JIN Tingxu，KOU Jisheng，et al.Lucas-Washburn equation-based modeling of capillary-driven flow in poroussystems J.Langmuir，2021，37（5）：1 623-1 636.11WASHBURN E W.The dynamics of capillary flow J.PhysicalReview，1921，17（3）：273-283.12BENAVENTE D，LOCK P，MNGELES G D，et al.Predicting the capillary imbibition of porous rocks from microstructureJ.Transport in Porous Media，2002，49：59-76.13TIAN Weibing，WU Keliu，CHEN Zhangxin，et al.Effect of dynamic contact angle on spontaneous capillary-liquid-liquid imbibition by molecular kinetic theory J.SPE Journal，2021，26（4）：1-16.14TIAN Weibing，WU Keliu，CHEN Zhangxin，et al.Dynamicwetting of solid-liquid-liquid system by molecular kinetic theoryJ.Journal of Colloid and Interface Science，2020，579：470-478.15WANG Han，SU Yuliang，WANG Wendong.Investigations onwater imbibing into oil-saturated nanoporous media：Couplingmolecular interactions，the dynamic contact angle，and the entrance effect J.Industrial&Engineering Chemistry Research，2021，60（4）：1 872-1 883.16ERICKSON D，LI D，PARK C B.Numerical simulations of capillary-driven flows in nonuniform cross-sectional capillaries J.Journal of Colloid and Interface Science，2002，250（2）：422-430.17YOUNG W B.Analysis of capillary flows in non-uniform cross-sectional capillaries J.Colloids and Surfaces A：Physicochemical and Engineering Aspects，2004，234（1/3）：123-128.18REYSSAT M，COURBIN L，REYSSAT E，et al.Imbibition ingeometries with axial variations J.Journal of Fluid Mechanics，2008，615：335-344.19BERTHIER Jean，GOSSELIN David，PHAM Andrew，et al.Spontaneous capillar