2024年2月6日发(作者：)

第29卷2022年3月第2期油气地质与采收率PetroleumGeologyandRecoveryEfficiencyVol.29,No.2Mar.2022文章编号：1009-9603（2022）02-0069-08DOI：10.13673/37-1359/te.2022.02.008两相低速非达西渗流模型及相对渗透率曲线求取方法222赵国忠1，，董大鹏1，，肖鲁川1，（1.中国石油大庆油田有限责任公司勘探开发研究院，黑龙江大庆163712；2.黑龙江省油层物理与渗流力学重点实验室，黑龙江大庆163712）摘要：流体在特低渗透-致密油藏中的渗流规律不符合达西定律，传统计算两相相对渗透率曲线的JBN方法并不适用。为了解决此类低速非达西渗流实验数据的处理问题，得到更为准确的相对渗透率曲线，在考虑两相渗流拟启动压力梯度的基础上，统计得到经验关系式，可计算不同空气渗透率岩样任意含水率下的拟启动压力梯度。依据低渗透岩样两相稳定流实验结果，提出两相流过程中不同相的拟启动压力梯度相等的假设，建立两相低速非达西渗流模型，基于该模型得出求取两相非达西渗流相对渗透率曲线的方法。应用大庆油田岩样实验数据，对比JBN方法和所提出的新方法，结果表明考虑两相渗流拟启动压力梯度的影响后，油相和水相相对渗透率都有明显变化。关键词：稳定流；非稳定流；致密油藏；非达西渗流；相对渗透率中图分类号：TE312文献标识码：ACalculationoftwo-phaserelativepermeabilitycurvesbasedonalow-velocitynon-Darcyflowmodel222ZHAOGuozhong1，，DONGDapeng1，，XIAOLuchuan1，（ationandDevelopmentResearchInstituteofDaqingOilfieldCo.，Ltd.，PetroChina，DaqingCity，HeilongjiangProvince，163712，China；gjiangKeyLaboratoryofReservoirPhysicsandFluidMechanicsinPorousMedium，DaqingCity，HeilongjiangProvince，163712，China）Abstract：Thefluidflowinultra-lowpermeability-tightoilreservoirsdoesnotconformtoDarcy’slaw，whichcannotbede⁃terminedbytheconventionalJohnson，Bossler，andNaumann（JBN）eringthepseudothresholdpressuregradientoftwo-phaseflow，piricalrelationwasemployedtocalculatethepseudothresholdpressuregradientatanywatercutofstatisticsforprocessingtheexperimentaldataoflow-velocitynon-Darcyflowandacquiringmoreaccuraterelativepermea⁃ingtotheexperimentalresultsoftwo-phasesteadyflowinlow-permea⁃bilityrocksamples，ahypothesiswasproposed，inwhichthepseudothresholdpressuregradientsofdifferentphasesinthetwo-phaseflowwereequal，nthismodel，amethodwasdevelopposedmethodwerecomparedwiththeJBNmetultsdemon⁃stratethattherelativepermeabilityoftheoilphaseandthewaterphasechangessignificantlyafterconsideringtheinflu⁃ds：steadyflow；unsteadyflow；tightoilreservoirs；non-Darcyflow；relativepermeability—————————————收稿日期：2021-01-21。作者简介：赵国忠（1964—），男，黑龙江海伦人，教授级高级工程师，硕士，从事油气渗流及油藏模拟方面的研究。E-mail：zhaoguozh@。基金项目：中国石油天然气集团公司“十三五”重大科技专项“大庆油气持续有效发展关键技术研究与应用”（2016-E02）。

·70·油气地质与采收率2022年3月在将达西定律应用于石油领域的过程中，两相相对渗透率曲线的获取是个重要环节。1959年，JOHNSON据求取两相相对渗透率曲线的方法等建立了利用岩心一维非稳定流驱替数［1-2］，已被石油行业标准所采用，适用于渗透率大于5mD的岩样。流体在特低渗透-致密油藏中的流动并不遵从达西定律［3-4］。闫庆来等依据低渗透岩心单相渗流实验结果推测，当压力梯度较小时渗流速度和压力梯度呈曲线关系［5］，通过将近似直线段延长交于压力梯度轴上而得到拟启动压力梯度，可以构建拟启动压力梯度模型［6］。然而，如何把单相流模型推广到多相流尚未形成广泛认同的做法，通常把单相流拟启动压力梯度直接用于多相流［7-10］，但这种做法缺乏实验和理论依据。开展两相非达西稳定流实验发现，两相流体整体的拟启动压力梯度并不等于单相时的值，且与含水率有关［11-13］。在此认识基础上建立了两相非达西渗流模型，认为两相情况下不同相的拟启动压力梯度相等且受控于空气渗透率和相分流率（或饱和度）［14］。目前为止，基于拟启动压力梯度依赖含水率的两相非达西渗流模型的相对渗透率曲线求取方法尚未严格建立起来［15-17］。以两相低速非达西稳定流实验为依据，分析各相拟启动压力梯度相等假设的合理性，给出油水两相非达西渗流模型，并采用JBN方法的思路，从源头出发严格给出非达西情况下理论推导的条件，依据恒速法非稳定流驱替实验数据，建立求取两相低速非达西渗流相对渗透率曲线的方法。1油水两相低速非达西稳定流实验采用稳态法，对大庆油田特低渗透岩心进行实验，研究油水两相的渗流特征。实验条件包括：实验温度为45℃；实验用油为模拟油，45℃时黏度为9.426液，45mPa·s℃时黏度为；实验用水为质量分数为0.639mPa·s。3%的KCl溶实验步骤如下：①对岩样依次进行抽真空饱和水、称重、测定孔隙体积、油驱造束缚水。②在束缚水饱和度下改变5次不同的油相渗流速度，记录稳定时的压差，测定束缚水时油相拟启动压力梯度。③油水流速比下改变依次改变油水流速比为5次总的渗流速度，4∶1，1∶1，1∶4记录稳定时，在每一个的压差，测定不同流速比下的油水拟启动压力梯度。④水驱至残余油，在残余油饱和度下改变5次不同的水相渗流速度，记录稳定时的压差，测定残余油时水相拟启动压力梯度。应用上述方法测试得到7块岩样不同含水率下的拟启动压力梯度，由岩样参数和实验结果（表1）可见，同一岩样在不同含水率下的拟启动压力梯度明显不同。所以两相情况下拟启动压力梯度不能直接使用单相的测试结果，必须考虑含水率的影响。随着含水率的增加，3号和4号岩样的拟启动压力梯度呈现先增后降的规律；而对于其他样品，虽然中间有个别波动，但总体也是呈现先增后降的变化趋势。Table1表1不同含水率下的拟启动压力梯度Pseudodifferentthresholdwaterpressurecutgradientat岩样渗透孔隙不同含水率下的拟启动编号率/mD度压力梯度（/MPa·cm-1）空气油相%/120.6202.44220%50%2.19280%100%31.140.0131.7410.16341.470.0510.8712.830.51910.440.2762.1911.9600.7000.043752.350.3263.870.9313.200.0190.6571.41974.700.6513.990.02513.110.0950.1930.8160.3730.2660.2160.1480.1220.3120.0046.850.361.2818.500.0160.1320.08850.2860.0010.014990.2830.2150.00430.0120.0169应用表1实验结果，建立最大拟启动压力梯度与空气渗透率的关系式，对每块岩样的拟启动压力梯度进行归一化处理，最终得到油水两相渗流时拟启动压力梯度的经验关系式为：λ=1.5038K(-3.2879fa-1.2189×w2+3.1423fw+0.2195)（1）对于不同空气渗透率岩样，（1）式可用于计算油水两相流体在任意含水率下的拟启动压力梯度。此关系式并非普适的和唯一的，实验数据增多或岩石、流体类型改变，可重复以上步骤得到不同的近似解析关系式，且所得关系式只能用于油水两相低速非达西渗流模型之中。基于不同的渗流模型也会得到不同的关系式，但显然不可互用。曾有学者采用类似的实验方法得到两相启动压力梯度与含水饱和度的线性关系［18-19］，但所用岩样与本文实验不同，且如何从不同流速下两相稳定流实验结果得到两相启动压力梯度的方法不详。因其考虑的是二次非线性渗流，其应用也必须基于这种渗流模型而建立起来的渗流理论，但目前还未形成相应的完整理论体系。

第29卷第2期赵国忠等.两相低速非达西渗流模型及相对渗透率曲线求取方法·71·2油水两相低速非达西渗流模型两相低速非达西稳定流动规律可描述为：两相不可压缩流体按一定比例通过岩石介质时总流速与压力梯度的关系，其表达式为：vìt=ïím(fwwï)∇p|-λ(f)î0[|]|∇p由于油、水按一定比例通过，等价于水按一定|∇p||>≤λλ((ffw)（2）w)分流率通过，所以水相流速与压力梯度的关系为：vìfw=ïíwm(fwï)|î0[|∇p|-λ(fw)]由（2）和（3）式知，2条直线在|∇|∇∇pp|p|轴上的截距，|>≤λλ((ffw)（3）w)即拟启动压力梯度相等。类似的关系也适用于油相。这就是说，在一定的分流率下，各相的拟启动压力梯度都与两相总拟启动压力梯度相等，只是流动直线关系的斜率不同。既然通过实验可以在一定的分流率条件下建立两相低速非达西稳定流，若把实验岩心看成油藏内的单元体，考虑实验持续时间相对油藏尺度下流体流动时间很短，可以认为油藏内任一点在一小段时间内的两相流动可以达到稳定。油藏内流体的非稳定流是由这样的无数稳定流构成的。其实达西条件下的两相流模型就是在此思想下建立的。因此，有理由假设两相低速非达西渗流各相的拟启动压力梯度相等，但其值与发生渗流位置的渗透率和当时的含水率有关。有学者在研究两相非线性渗流启动压力梯度时并未区分各相的启动压力梯度，已隐含了两相启动压力梯度相等的假设［18］。本文的研究使用的是相拟启动压力梯度这个参量，并假设其相等，这与文献［18］隐含假设是类似的，具有同样的道理。在该假设条件下，一维均质不考虑毛管压力情况下油水两相不互溶流体低速非达西渗流模型可表示为：ìïKKri(Sw)é||v=ïí-μê|ë|∂p||i|∂x|+λ(fw)ùúû|∂iïï|∂px|||>λ(fw)î0||∂∂px|||≤λ(fw)（4）选取基准渗透率为束缚水时非达西流条件下的油相有效渗透率，等价于定义束缚水时油相相对渗透率为1，即：K=μom(0)（5）多孔介质的渗透率是由达西定律定义的。须指出的是，在非达西模型下流动不再满足达西定律，渗透率的定义也该相应作出关联（3）或（4）式且形如（5）式的调整，才能构成完整理论体系，单相的情况也类似。这样定义的渗透率的量纲不变，且在忽略拟启动压力梯度时即为原达西流时的定义。3油水两相非达西渗流情况下JBN方法的适用条件方程和威尔奇积分的基础之上，JBN方法是建立在巴克利-列维里特前缘推进利用岩心一维非稳定流驱替数据获取两相相对渗透率曲线的方法［1］。由低速非达西渗流模型（4）式，可得出发生两相流动条件下含水率的表达式为：fvw=vwKroμwo+vw=(1+Krwμo)-1=fw(Sw)（6）这表明非达西渗流情况下水相分流率仍然可表示为含水饱和度的函数，而黏度比只是参数，对于确定的流体其为常数。考虑不可压条件，由相连续性方程和饱和度约束方程仍然可以得到巴克利-列维里特前缘推进方程［20］：(ddxt)t)7）sw=v(ϕfw′(Sw)（（7）式表征了一维驱替过程中t时刻含水饱和度为S度为Lw的推进规律。可将其用于横截面积为A、长的一维均质岩心内，忽略毛管压力的油水两相不互溶流体低速非达西驱替实验。忽略毛管压力是在特定条件下研究渗流单一影响因素的常用做法，也是建立基于非稳定流驱替实验数据求取相对渗透率曲线方法的前提条件。在使用驱替实验数据求取相对渗透率曲线的过程中，所涉及的出口端饱和度范围处于前缘饱和度和残余油时水相饱和度之间（这是非稳定流法求取相对渗透率曲线的局限性，达西流时也是如此），而入口端饱和度一直等于残余油时水相饱和度。日常由特低渗透岩样退汞曲线可转换得到油水两相毛管压力曲线，在此饱和度区间的毛管压力一般小于0.1以用于研究特低渗透岩心驱替的模型可忽略毛管MPa，远低于特低渗透驱替实验的驱替压差，所压力。在应用本文方法所得相对渗透率曲线开展

·72·油气地质与采收率2022年3月油藏工程研究或数值模拟计算时，可根据具体情况在模型中随时考虑毛管压力、重力等其他项的作用。对（7）式积分可得t时刻含水饱和度为S进距离为：w的推tx(vst)=∫(t)d(t)fw′(Sw)=V(t)w0ϕϕAfw′(Sw)（8）由（8）式可得t时刻含水饱和度分布剖面，此剖面存在一个前缘饱和度Sw∗，满足的条件是［20］fw′(Sw∗)=fSww∗(-Sw∗S)：（9）wc而剖面上的饱和度范围为：Sw∗≤Swe≤Sw≤1-Sor（10）由于JBN方法求取的相对渗透率与含水饱和度的关系曲线是针对驱替实验出口端，由（10）式知，理论上讲JBN方法只能得到前缘饱和度Sw∗后的数据点。这一点是无论达西流与否都应该注意的适用范围，在此之外都没有理论依据，是不可信的。由（8）式还可看出，在水驱前缘到达出口后的任意时刻，含水饱和度、含水率、含水率导数和坐标x具有一一对应（映射）关系，即它们之间都是两两单值函数关系。（8）式用于出口端可得［1］：L=x(t)s(t)=VweϕAfw′(Swe)⇔Wi-1=ϕAL(=fw′(Swe)=fwe′Vt)（11）为后续推演方便，（8）式可进一步写成［1］：x=Lffwe′w′（12）基于低速非达西流模型即（4）式，由于以上诸条件的成立，威尔奇积分仍然成立，亦即存在岩心出口端含水饱和度与平均含水饱和度的关系为［2］Swe=Sw-1f-f：-wewe′=-Sw-Wi(1-fwe)=-Sw+fwe′平均含水饱和度可依据驱替数据确定：dd-Sfwwe′（13）-S=SVwc+ϕALo(t)w（14）而岩心出口端水相分流率或含水率可表示为：fwe=1-foe=1-dVo(t)=1d[ϕALd(--SV(wt-)Swc)]dV(t)=1-d-1+(dWSwi=fwe′)2dd-Sfwwe′（15）4油水两相低速非达西流相对渗透率公式及应用实例替，通过积分得到岩心两端压差与出口端流动特征JBN方法的关键是把流动方程运用于岩心驱的关系，并由此定义一个可称为相对注入能力的时间变量，然后再对注入孔隙体积倍数的倒数，即fwe′求导，从而得到出口端相流度与相对注入能力的关系。在此基础上进行了公式推导及应用。4.1相对渗透率曲线计算公式由（4）式可把非达西渗流模型在流动发生时刻t的水相方程（亦可选油相方程）表示为：-∂∂px=μwv(t)fΚwΚrw(S（16）式沿岩心长度积分可得岩心两端压差(Sw)w)+λ(fw)（16）对为：Δp(t)=∫Léêμwv(t)fw0êëΚΚrw(Sw)+λ(fùw)úúdx（17）利用（12）式对（17）式换元可得：(Sw)ûΔp(=Léêμwv(t)fwe′fwfwe′êëΚ∫fw0′Κrw(S(Sw)t)w)dfw′+∫fwe′λ(f)df′úùfw0′wwúû（18）把（4）式用于驱替开始时的油相，再考虑（5）式，得：Κμv(0)=Δp(0)-oLLλ(0)（19）若驱替实验采用恒速法，即v(t)=v(0)且为常数，则：μwv(t)μΚ=wv(0)μΚ=w[Δp(0)-Lλ(0)]μoL=λr（20）

第29卷第2期赵国忠等.两相低速非达西渗流模型及相对渗透率曲线求取方法·73·定义相对注入能力为：Ω=Κμwv(t)ΔpL(t)=Δλpr(Lt)（21）则（18）式可表示为：∫fwe′ffww0′Κrw(f(Swe′Swww′fw)（22）式两)dfw′=Ωfwe′-λ1r∫w0′λ(f)df（22）边对ffwe′微分并注意到fw0′=w′|1-Sor=0，整理可得：Krwe=Krw(Swe)=fd(wΩf(Swe)=dfwe′we′)-λ[fwλ(Swer)](fdΩf′we（23）dfwewe)we′)-λ(fλr将（23）式代入（6）式，可求得：Kroe=Kro(Swe)=μo(μ1-fwewfwe)Krwe（24）4.2相对渗透率曲线求取步骤综合两相非达西渗流模型及相对渗透率计算公式，计算步骤可归纳为：①由（14）式利用驱替过程中的产油量数据计算并回归得到平均含水饱和度与注入孔隙体积倍数倒数的解析关系，()即-S-wSwf=we′。②由（15）式求得出口端含水率与注入孔隙体积倍数倒数的解析关系，即fwe=fwe(fwe′)。③由20）式利用前期实验得到的拟启动压力梯度经验公式λ=λ(fw)和驱替实验起始压差计算参考压力梯度λr。④由（21）式利用驱替过程中的压差数据计算并回归得到相对注入能力与注入孔隙体积倍数倒数的解析关系，即Ω=Ω(fwe′)。⑤由（23）和24）式求得出口端两相相对渗透率与注入孔隙体积倍数倒数的解析关系，(即K′)rwe=Krwe(fwe′)和Kroe=Kroefwe。⑥利用步骤①的结果由（13）式求得出口端含水饱和度与注入孔隙体积倍数倒数的解析关系，即Swe=Swefwe′。⑦依据步骤⑤和⑥得到的解析关系算出一系列注入孔隙体积倍数倒数对应的出()口端含水饱和度和相对渗透率，作图得到两相相对渗透率曲线。以上步骤与达西流情况的经典JBN方法相比，步骤③是增加的，步骤④和⑤是不同的。需要指出的是，不同于达西流的情况，如果未采用恒速法驱替实验，则不能由（22）式对fwe′微分而推导出（23）式。所以非达西流的相对渗透率求取须使用恒速法驱替实验数据。4.3相对渗透率曲线求取实例以大庆油田某特低渗透岩心为例验证算法的可行性。实验基本参数包括：岩心长度为8.08cm，直径为2.52cm，孔隙体积为4.92mL，孔隙度为12.20%42.07%，气测渗透率为1.131mD，束缚水饱和度为mD文所选基准渗透率为，对应的束缚水时非达西流油相有效渗透率即本，束缚水时达西流油相有效渗透率为0.04760.0516mD，油的黏度为5.20mPa·s见表2，。水的黏度为0.57mPa·s。所记录的实验数据Table2表2Experimental大庆油田某特低渗透岩心实验数据corefromdataDaqingofaultra-lowOilfieldpermeability累积时产油体产液体注入孔隙平均含水压差间/min积/mL积/mL体积倍数饱和度/%00000.420MPa/6.000.600.542729.57.160.600.720.1220.552727.114.509.530.650.700.950.14500.563826.422.500.751.450.19350.573025.732.000.802.250.29470.45770.583124.947.500.850.593324.067.500.903.200.65030.603523.30.954.750.613622.7140.0097.501.006.750.96540.624822.4200.001.0514.009.751.37241.98100.634022.31.1020.002.84570.644122.2290.001.1529.004.06550.654322.1410.001.205.89400.664422.0530.001.2041.008.333330.664621.9650.001.2053.00770.001.2065.0010.7720.664621.9990.001.2077.0013.211499.0015.650420.122400.664621.90.6646621.921.9分别采用本文非达西流算法和达西流JBN方法，计算得到相对渗透率曲线，结果如图1所示。在计算过程中，平均含水饱和度与注入孔隙体积倍数倒数的拟合函数关系为：-Sw=a1+a2fwe′+a3e-a4fwe′（25）（（

·74·油气地质与采收率2022年3月其中，a1非达西流情况下，=0.592，a2=-0.005相对注入能力与注入孔隙体55，a3=0.0818，a4=1.5。积倍数倒数的拟合函数关系为：Ω=b′b1+b2fwe3+b4e-fwe′（26）其中，b1=7.34，b2=0.36，b3(=0.766)对于达西流，考虑与非达西流情况的对比需，b4=-0.0683。要，基准渗透率仍然选择非达西流条件下的油相有效渗透率，则相对注入能力应为：Ω=ΚΚοΔp(t)×Δvp((t0)=ΚΚο×ΔΔpp((0t))（27）v(0))达西流时相对注入能力与注入孔隙体积倍数倒数的拟合函数关系式与非达西流情况相同，但对应的b1b2平均含水饱和度及相对注入能力与注入孔隙=0.6845，=0.0336，b3=0.766，b4=-0.00637。体积倍数倒数的拟合函数关系有个共同的特点，即当注入孔隙体积倍数倒数趋于0或注入孔隙体积倍数无穷大时，拟合函数的极限存在。这对相对渗透率端点值残余油饱和度及其对应的水相相对渗透率的确定提供了合理可行的办法，因为残余油饱和度是极限驱替条件下的含油饱和度［21-26］，实验达不到这个理想状态［27-28］。计算结果表明，对于特低渗透岩心，考虑两相渗流拟启动压力梯度的影响后，油相和水相相对渗透率在等渗点之前都变高了；在等渗点之后油相近乎相同，而水相略微降低。这里应指出，非达西流相对渗透率曲线只可以用在基于非达西流模型的相关研究或应用中；而达西流相对渗透率曲线只可以用在基于达西流模型的相关研究或应用中。二者不应混淆。Fig.1图1Comparison非达西流与达西流相对渗透率曲线对比non-DarcybetweenflowrelativeandDarcypermeabilityflowcurvesof5结论特低渗透岩心油水两相渗流时，总流速与压力梯度之间表现出明显的非达西渗流特征，也存在启动压力梯度。拟启动压力梯度随着含水率的增加，呈现出先增加后降低的变化规律，依据两相低速非达西稳定流实验结果，对拟启动压力梯度与含水率的关系曲线进行归一化，为二次多项式关系。通过对两相低速非达西稳定流实验结果的进一步分析，论证了一定的分流率下各相拟启动压力梯度相等假设的合理性，建立了油水两相低速非达西渗流模型。导出了类似于JBN方法的油水两相低速非达西流非稳定驱替相对渗透率曲线求取方法。新方法考虑了两相流拟启动压力梯度的影响，给出了相对注入能力的新定义及油水相对渗透率新的计算公式。该方法仅适用于恒速法非稳定流驱替数据。与传统基于达西流的JBN方法相比，所得油相和水相相对渗透率均有明显不同。符号解释A——横截面积，cm2；ffoe——出口端含油率，%；w——含水率，%；fwe——出口端含水率，%；fwe′——注入孔隙体积倍数倒数；f下标w0——入口端含水率，%；i——o,w，分别表示油和水；K——基准渗透率，mD；Ka——空气渗透率，mD；KKri(Sw)——相对渗透率；ro——油相相对渗透率；Krw——水相相对渗透率；K—出口端的水相相对渗透率；L—rwe——长度，cm；m(fw)——受控于fw的系数，表示直线斜率，即有效流度；p——压力，MPa；∇Sp——压力梯度，MPa/cm；or——残余油饱和度；Sw∗——前缘饱和度；Sw——含水饱和度；Swc——束缚水饱和度；Swe——水驱前缘到达出口后出口端的含水饱和度；t——时间，s；

第29卷第2期赵国忠等.两相低速非达西渗流模型及相对渗透率曲线求取方法·75·vi——相流速，cm/s；vo——油相流速，cm/s；vt——总流速，cm/s；v(t)——t时刻总流速，cm/s；vVw——水相流速，cm/s；(t)——t时刻累积注水量，cm3；Vo(t)——累积出油量，cm3；Wi——注入孔隙体积倍数；x——横坐标，cm；xsw(t)——推进距离，cm；xswe——对应于出口端的推进距离，cm；λ——拟启动压力梯度，MPa/cm；λ(fw)——受控于fw的系数，表示直线在|∇p|轴上的截距，即拟启动压力梯度，MPa/cm；λμr——参考压力梯度，常数，MPa/cm；i——相黏度，mPa·s；μo——油相黏度，mPa·s；μϕ—w——水相黏度，mPa·s；—孔隙度，%；Ω——相对注入能力。参考文献1］JOHNSONEF，BOSSLERDP，ationofrelativeTransactionspermeability，AIME，1959from，displacement216（1）：ment［J］.Petroleum2］rive［J］.JournalofPetroleumTechnology，1952，4（4）：3］赵国忠，尹芝林，匡铁，等.大庆油田水驱油藏模拟特色技术［J］.大庆石油地质与开发，2019，38（5）：mulatingGuozhong，YINZhilin，KUANGOilfield［J］eumofGeologythewater&OilfieldfloodedTie，teristicininDaqinging，2019，38（5）：⁃4］李承龙，苗志国，李照永，等.大庆长垣外围特低渗透油藏开发效果综合评价方法［J］.特种油气藏，2019，26（4）：iveChenglongultra-lowevaluation，MIAOofdevelopmentZhiguo，LIZhaoyongperformance，enewperipheralcompre⁃field［J］.SpecialpermeabilityOil尉立岗，&Gasoilreservoirs等Reservoirsin.低渗透油层中单相液体渗流特征，2019Changyuan，26（4）：Oil⁃5］闫庆来，何秋轩，的实验研究［J］.西安石油学院学报：自然科学版，1990，5（2）：inglai，HEQiuxuan，WEILigang，atorystudybilitypercolationSciencereservoirscharacteristicsofsinglephaseflowinlow-permea⁃Edition［，1990J］.Journal，5（2）of：’anPetroleumInstitute：Natural6］肖鲁川，甄力，郑岩.特低渗透储层非达西渗流特征研究［J］.大庆石油地质与开发，2000，19（5）：27-28，chuan，ZHENLi，-Darcyseepageflowcharacteristicology&OilfieldinDevelopmentultra-lowpermeableinDaqingreservoirs，2000，［19J（］.Petroleum5）：27-28，Ge⁃［7］宋付权，刘慈群.低渗油藏的两相相对渗透率计算方法［J］.30.西安石油学院学报：自然科学版，2000，15（1）：10-12，rmeabilitiesFuquan，JournalofXi’eumandAoilmethodInstituteinlow：permeabilityforcalculatingNaturalSciencereservoirstherelativeEdition［J］，.［8］邓英尔，2000，15（王允成，1）：10-12刘慈群，，20.等.低渗非达西渗流相对渗透率计算方法及特征［J］.西南石油学院学报，2000，22（3）：lativeYingerpermeability，WANGofnon-DarcyYuncheng，flowLIUinCiqunlowpermeability，ating［J］.JournalofSouthwestPetroleumInstitute，2000，22（3）media：34-［9］董大鹏，36.冯文光，赵俊峰，等.考虑启动压力梯度的相对渗透率计算［J］.天然气工业，2007，27（10）：sureforgradientrelativeDapeng［Jpermeability，FENGWenguang］.NaturalGascalculation，ZHAOIndustry，2007consideringJunfeng，，27（10）⁃：⁃［10］周英芳，王晓冬，李斌会，等.低渗油藏油水相对渗透率非稳态计算方法［J］.大庆石油地质与开发，2010，29（3）：permeabilitiestheunsteady-stateYingfang，WANGinlow-permeabilitycalculationXiaodongmethod，LIBinhuiofoil，esearchesrelativeGeology&OilfieldDevelopmentinDaqingoilreservoirs，2010，29［（J］3）.Petroleum［11］肖鲁川，曹维政.特低渗透储层油水两相非达西渗流实验研究：93-97.［R］.大庆：大庆油田有限责任公司勘探开发研究院，2002：epageLuchuan，noil-waterNon-DarcyDaqingingOilfield：ExplorationexperimentsCo.，Ltd.，andin2002Developmentultra-lowpermeability：chInstitutereservoirsof［Daq⁃R］.［12］曹维政，肖鲁川，曹维福，等.特低渗透储层油水两相非达西渗流特征［J］.大庆石油地质与开发，2007，26（5）：spermeabilityofWeizhengnon-Darcy，XIAOreservoirsflowLuchuan［Jvelocity，］.PetroleumofCAOoil-waterWeifu，teris⁃Geologyphases&Oilfieldinextra-lowmentinDaqing，2007，26（5）：p⁃［13］石京平.低渗透油田油水两相低速非达西渗流规律研究［D］.成都：成都理工大学，2006：noil-waternon-Darcyseepagerulesinlow赵国忠nology，2006.变启动压力梯度三维三相渗流数值模拟方法：oilfield75-79.［D］.Chengdu：ChengduUniversityofTech⁃［14］［J］.石油学报，2006，27（z1）：119-123，lestart-upNumericalsimulationof3Dandthree-phase15］杨小平nica，2006.精确计算相对渗透率的方法，27（z1）：119-123pressure，nt［J］.ActaPetroleiSi⁃［［J］.石油勘探与开发，1998YANG，25（Xiaoping.6）：dtocalculatetherelativepermeability（accurately6）：63-66.［J］.PetroleumExplorationandDevelopment，1998，25［［［［［［

·76·油气地质与采收率2022年3月［16］陈忠，殷宜平，陈浩.非稳态法计算油水相对渗透率的方法探讨［J］.断块油气田，2005，12（1）：ong，YINYiping，culationmethodofOilandGasField，2005，12（1）：-waterrelativepermeabilitybyunsteady-state［J］.Fault-Block［17］杨正明，李治硕，王学武，等.特低渗透油田相对渗透率曲线测［23］李金宜，陈丹磬，朱文森，等.一种对残余油饱和度合理值范围的预判方法［J］.断块油气田，2016，23（3）：i，CHENDanqing，ZHUWensen，dgewayforBlockOilandGasField，2016，23（3）：ablevaluerangeofirreducibleoilsaturation［J］.Fault-［24］张伟，曹仁义，罗东红，等.南海珠江口盆地海相砂岩油藏高倍试新方法［J］.石油学报，2010，31（4）：voirsfor［JtestingZhengming］.Actarelative，LIPetroleipermeabilityZhishuo，WANGSinica，2010of，31ultra-lowXuewu，（4）：th⁃reser⁃［18］李爱芬，刘敏，张化强，等.低渗透油藏油水两相启动压力梯度变化规律研究［J］.西安石油大学学报：自然科学版，2010，25（6）：47-50，54．LIonAifenlowthreshold，LIUpermeabilitypressureMin，ZHANGreservoirsgradientHuaqiang［J］.Journalofoil-water，-phaseExperimentalseepagestudyin赵国忠，ty：Natural杨清彦，Science唐文锋，Edition等，2010.大庆长垣外围低渗透油田开发机，25（6）of：47-50Xi’an，si⁃［19］理［J］.大庆石油地质与开发，2009，28（5）：mentGuozhong，YANGvoir［J］.mechanismPetroleumGeologyofDaqingQingyan&peripheral，TANGOilfieldDevelopmentlow-permeabilityWenfeng，⁃inDaqingreser⁃2009，28（5）：126-133.，［20］科林斯RE.流体通过多孔材料的流动［M］.陈钟祥，吴望一，译.北京：石油工业出版社，1984：SZhongxiangRPress，1984，：yioffluids，gporous：materialsPetroleum［M］［21］雷霄，查玉强，姜平，等.一种改进的油水相渗实验末端效应校正方法［J］.中国海上油气，2016，28（5）：nXiaomethod，ZHAforYuqiangend-effect，JIANGinoil-waterPing，ovedpermeabilitycorrec⁃ex⁃［22］许家峰，periment［张金庆，J］.China安桂荣，Offshore等Oil.利用水驱曲线动态求取残余油水andGas，2016，28（5）：49-53．相渗透率的新方法［Ｊ］.中国海上油气，2014，26（1）：ulateJiafeng，ZHANGJinqing，ANGuirong，thodtoingandthetherelativepermeabilityofwaterinresidualoilbyus⁃Gasperformance，2014，26（1）of：rivecurve［J］.ChinaOffshoreOil数水驱驱替特征［J］.油气地质与采收率，2018，25（2）：haracteristicsWei，CAOreservoirsintheofPearlhigh-multipleRenyi，LUORiverMouthwaterDonghongBasindrive，，cementsandstonePetroleumGeologyandRecoveryEfficiency，2018South，25China（2）：Sea64-71.［J］.25］吕晶，刘显太，孙业恒，等.水相分流量全过程有理式模型的建立及应用［J］.油气地质与采收率，2020，27（6）：100-105.LÜplicationJing，eumwholeXiantaiprocess，SUNGeologyrationalYehengandformula，ishmentRecoverymodelEfficiencyofwater，2020fraction⁃andap⁃（al6）flow：100-105.［J］，2726］秦国伟，蒲春生，罗明良，等.不同粘弹性驱替液下毛管数对驱油效果的影响［J］.石油天然气学报，2007，29（2）：illaryGuowei，PUChunsheng，LUOMingliang，idsnologywith，2007differentnumbers，29（nceof2）：97-100.［J］.JournaleffectusingofOildisplacingandGasTech⁃flu⁃27］董大鹏.非稳态相渗实验数据的处理方法［J］.西南石油大学学报：自然科学版，2014，36（6）：singUniversityunsteadystatemethodtothedataofrelativepermea⁃：Sciencedisplacement&Technology［J］.JournalEditionof，Southwest2014，36（6Pe⁃m）：28］高慧梅，姜汉桥，陈民锋，等.储集层微观参数对油水相对渗透率影响的微观模拟研究［J］.石油勘探与开发，2006，33（6）：ulationHuimeiwaterrelativeof，theJIANGHanqiao，CHENMinfeng，copicpermeabilityeffectofmicroscopic［J］.PetroleumreservoirExplorationparametersandonDevel⁃oil-opment，2006，33（6）：734-737.编辑刘北羿［［［［