低渗含水气藏非达西渗流规律及其应用

徐轩，王继平，田姗姗，焦春艳，张宇轩

1.中国石油勘探开发研究院廊坊分院，河北廊坊 065007

2.中国矿业大学(北京)地球科学与工程测绘学院，北京海淀 100083

3.中国石油长庆油田苏里格研究中心，陕西西安7100184.中国石油天然气管道工程有限公司，河北廊坊 065000

5.中国石油新疆油田公司采油一厂地质所，新疆克拉玛依 834000

低渗含水气藏非达西渗流规律及其应用

徐轩1，2*，王继平3，田姗姗4，焦春艳1，张宇轩5

1.中国石油勘探开发研究院廊坊分院，河北廊坊 065007

2.中国矿业大学(北京)地球科学与工程测绘学院，北京海淀 100083

3.中国石油长庆油田苏里格研究中心，陕西西安710018
4.中国石油天然气管道工程有限公司，河北廊坊 065000

5.中国石油新疆油田公司采油一厂地质所，新疆克拉玛依 834000

气藏开发实践和室内实验研究均表明低渗含水气藏气体渗流存在启动压力现象，研究启动压力梯度的影响因素和变化规律对气藏开发具有重要意义。通过长岩芯多测点气藏物理模拟实验方法，研究了低渗含水砂岩气藏气相启动压力梯度变化规律，对该参数与储层渗透率、含水饱和度进行了相关性分析，建立了低渗含水气藏启动压力梯度计算通式。结合气藏渗流理论，采用建立的启动压力梯度表征方法，针对典型气藏进行应用分析，确定了气井有效动用半径，形成了有效动用半径与储层渗透率和含水饱和度的关系图版，为类似气藏开发早期井网部署以及开发中后期井网加密提供了技术指导。

含水气藏；非达西渗流；启动压力梯度；实验研究；定量表征

徐轩，王继平，田姗姗，等.低渗含水气藏非达西渗流规律及其应用[J].西南石油大学学报(自然科学版)，2016，38(5)：9096.

XU Xuan，WANG Jiping，TIAN Shanshan，et al.Non-Darcy Seapeage Law and its Application in Low Permeability Water-bearing Gas Reservoir[J]. Journal of Southwest Petroleum University(Science&Technology Edition)，2016，38(5)：9096.

引言

随着天然气需求量的快速增长，越来越多的低渗气藏投入开发。在开发实践中，低渗气藏的生产特征表现出气体渗流存在启动压力，这些现象的出现使得低渗气藏渗流机理越来越得到科研工作者的重视[1—4]。任晓娟、贺伟等[5—6]针对含水低渗气体低速非达西渗流进行的实验研究表明：与低渗油藏储层中的渗流不同，低渗气藏气体启动压力梯度同时受储层含水饱和度和渗透率影响。大量学者从单因素的角度研究了其变化规律：含水饱和度方面，刘晓旭、依呷等[7—8]通过室内实验发现启动压力梯度与含水饱和度成正比关系。渗透率方面，对于低渗油藏前人研究已明确渗透率对储层启动压力起决定性作用[9—18]，对低渗气藏前人得出了气体启动压力梯度与储集层渗透率的统计规律，即其大小与渗透率成反比。

那么，气体渗流的启动压力梯度究竟如何受含水饱和度和渗透率的共同影响？是否具有统计规律？该如何定量表征？这些问题，都需要进一步的深入研究。因此，笔者通过大量实验，系统研究了启动压力梯度变化规律。在此基础上，分析建立了综合考虑含水饱和度和渗透率影响的启动压力梯度定量表征方法，并针对实际区块进行了应用分析。

1 实验原理与方法

实验采用自主研发的长岩芯多测点气藏开发物理模拟系统进行，该系统特点是岩芯夹持器长度可达100 cm，岩芯胶皮套上每隔一定距离部署测压孔，通过接头连接到压力传感器，可以精确测量开采过程中沿程压力分布。通过对岩芯注气后衰竭开采模拟可以观察到：对于低渗含水砂岩，随着气体采出，各测点压力开始下降，衰竭开采至停止流动时，各测点压力不再变化，此时与常规砂岩不同的是，各测点压力并不直接降为零，而仍然会长时间维持一定压力值。从岩芯远端边界到出口端可测到明显压力剖面(图1)，认为此时岩芯两端临界压差即为气相渗流启动压力。图1为平均孔隙度10.5%，渗透率0.06 mD，长度为52 cm的模型在不同含水饱和度条件下，衰竭开采至停止产气时测得的压力剖面。通过拟合测点压力与出口端距离可得到气体渗流的启动压力梯度，在含水饱和度分别为35%和58%时测得的启动压力梯度分别为0.13 MPa/m和1.64 MPa/m。

图1 衰竭开采至停止产气时岩芯压力剖面Fig.1Pressure profile when there is no gas flow in natural depletion

本次实验采用现场取得的某低渗含水气藏砂岩，岩样渗透率分布在0.03～1.60 mD，孔隙度分布为5.87%～12.85%。具体实验步骤为：(1)根据研究内容，选择物性一致的砂岩串联得到长岩芯(对于含水模型，建立含水饱和度)；(2)将模型饱和气至设计的孔隙压力后关闭气源，静置平衡一段时间，使模型处于统一压力系统；(3)从模型出口端释放气体，模拟气藏衰竭开采直到停止产气且各测压点压力保持不变，此时记录各测点压力；(4)根据压力与测点到出口端距离的关系拟合求得气相渗流启动压力梯度，结束实验，开始下次实验。

2 气体启动压力梯度规律研究

2.1 启动压力梯度与含水饱和度的关系

采用上文所述方法对不同含水饱和度下的气体渗流启动压力梯度进行了一系列的实验测试，统计分析得实验结果见图2。从实验可以看出，储层启动压力梯度具有规律性，总体上随原始含水饱和度增加而快速增大，相同含水饱和度时，渗透率越低启动压力梯度越大。这是由于在较高含水饱和度下，岩芯中的气体并不能形成连续相，而是被分割成许多小气泡进行流动，由于毛管力的作用，这些小气泡在每个喉道处都产生贾敏效应，形成流动阻力，含水饱和度越高，阻力越大。在宏观上就表现为岩芯含水饱和度越高，启动压力梯度值越大。依呷等[8]拟合了启动压力梯度与含水饱和度的关系式，认为启动压力梯度和含水饱和度Sw正相关。为了进一步明确启动压力梯度随含水饱和度的变化规律，并探索其定量表征方法，根据实验测试结果进行了数据拟合。采用试算法对启动压力梯度和含水饱和度Sw进行幂指数拟合，经过多次拟合，比较拟合效果，发现采用含水饱和度Sw的3.5次幂进行拟合效果较好，相关系数最大，见表1，图3。

分析拟合结果可知，启动压力梯度与含水饱和度Sw的变化规律符合幂函数关系式。不同渗透率岩样表达式具有如下共同形式。

图2 不同含水饱和度时的启动压力梯度Fig.2Starting pressure with different water saturation

表1 启动压力梯度与含水饱和度拟合公式Tab.1Fitting formula of starting pressure and water saturation

图3 启动压力梯度与含水饱和度关系曲线Fig.3Relationship between starting pressure and water saturation

式中：λ—启动压力梯度，MPa/m；

Sw—含水饱和度，无因次；

C1，C2—正常数，与储层的渗透率及孔隙结构有关。

对于岩芯S1，S2和S3，拟合得到的C1分别为16.672 0，8.932 4和3.578 3。可见，C1与渗透率密切相关，渗透率越低，启动压力梯度越大，相应的系数C1越大。拟合得到不同岩样C2差别不大，分别为0.019 0，0.018 4和0.025 6，该系数与启动压力梯度对渗流的阻碍作用相反，体现出滑脱效应的影响。

2.2 启动压力梯度与渗透率的关系

实验系统测定了多块岩芯在含水饱和度相同(含水饱和度分别为35%，47%，58%)时的气体渗流启动压力梯度，见图4。研究表明含水饱和度相同的情况下，储层启动压力梯度和平均渗透率倒数呈线性规律，具有很好的统计意义。

图4 启动压力梯度与渗透率倒数的关系曲线Fig.4Relationship between starting pressure and the inverse of permeability

分析拟合结果可知，不同饱和度岩样启动压力梯度与渗透率的变化规律具有如下共同形式

式中：

K—渗透率，mD；

C3，C4—正常数，与储层的含水饱和度及孔隙结构有关。

饱和度分别为35%，47%和58%时，拟合得到的系数C3分别为0.015 1，0.040 4和0.084 5。可见，C3与含水饱和度正相关，含水饱和度越高，启动压力梯度越大，相应的系数C3越大。拟合得到的不同饱和度下C4差别不大，分别为0.019 4，0.019 7和0.025 2，与式(1)中C2类似，该系数与启动压力梯度对渗流的阻碍作用相反，体现出滑脱效应的影响。

2.3 启动压力梯度的计算通式

上文通过对启动压力梯度随含水饱和度和渗透率变化规律的研究，统计拟合得到了单因素时相应的表达式。综合单因素研究成果，不难得出，综合考虑含水饱和度和渗透率影响时，启动压力梯度应可统一到如下通式

式中：

a，b—正常数。

当渗透率为定值时，启动压力梯度只受含水饱和度的影响，式(3)退化为式(1)，此时aK-1即为常数C1，b即为常数C2。当含水饱和度为定值时，启动压力梯度只受渗透率的影响，式(3)退化为式(2)，此时aS3w.5即为常数C3，b即为常数C4。

通过以上分析表明，启动压力梯度计算通式(3)能够同时反映含水饱和度和渗透率对启动压力梯度影响。这一通式中系数a，b可以通过启动压力梯度与含水饱和度的关系式(1)计算，也可以由启动压力梯度和渗透率的关系式(2)计算。对于同一区块或者岩性和孔隙结构相近的储层，通过式(1)或式(2)获得的a，b应该有同样或基本一致的取值。为了验证通式(3)的适用性，分别采用上述两种方法，计算a，b值，计算结果如表2所示。

表2 启动压力梯度计算通式中参数的拟合Tab.2Parameters fitting of starting pressure calculation

通过上述两种方法计算的a，b值基本一致，表明通式是适用的。取a，b平均值，可得研究区储层启动压力梯度计算通式为

值得指出的是，通过实验测试和拟合公式计算均表明，对于渗透率较高且含水饱和度较低的储层(如本例中当渗透率大于0.5 mD且含水饱和度小于30%时)，其启动压力梯度为零或负值，这反映出气体滑脱效应的影响。由于其值较小且滑脱效应不是本文研究重点，所以认为这种情况下，储层不存在启动压力梯度。

综上所述，研究认为，综合考虑含水饱和度和渗透率影响的启动压力梯度定量表征方法具有较强的统计意义，对于同一区块或者岩性和孔隙结构接近的储层通式是适用的。通式可以通过启动压力梯度与含水饱和度的关系式(1)或启动压力梯度和渗透率的关系式(2)计算得到。一旦拟合得到通式，便可方便估算某一研究区域不同储层，不同原始含水饱和度下的启动压力梯度。该方法便于实际操作，实用性较强。

3 应用分析

某含水低渗气藏，现场取芯分析得到气藏83%岩样渗透率小于1.00 mD，主要分布于0.01～2.00 mD，65%以上岩样孔隙度小于9%。气藏气水关系复杂，原始含水饱和度分布在30%～60%，属于典型的低渗含水气藏。

3.1 建立研究区启动压力梯度关系图版

气藏渗透率较低，原始含水分布较广。不同气井，不同部位储层的渗透率和原始含水饱和度均不相同。为获取不同条件下储层启动压力梯度，选取部分岩芯进行启动压力梯度测试。采用上文所述方法，仅开展少量实验测试即拟合得到该区块启动压力梯度计算通式

利用启动压力梯度通式可快速获得不同储层不同含水饱和度下启动压力梯度。不失一般性，选择3个代表性的饱和度值进行计算，并建立研究区不同渗透率和含水条件下的启动压力梯度关系图版，图5。

图5 气藏不同渗透率和含水条件下的启动压力梯度Fig.5Starting pressure of gas reservoirs with different permeability and saturation

3.2 计算有效动用半径

该气藏原始地层压力约为32 MPa，最小井底流压为2 MPa，井筒半径取0.1 m，气井物理边界半径为800 m。由于启动压力梯度的存在，气井有效动用半径往往达不到物理边界，不能完全控制储量，或者虽然能完全控制储量，但开发效益不一定很好。应用储层启动压力梯度规律，可计算出不同储层，不同含水饱和度下单井的有效动用半径，为开发提供指导。

气藏开发过程中，随着井底流动压力的降低，供给边缘不断扩大，当井底流动压力下降到最小值时，供给边缘达到最大，此时供给边缘处压力梯度应等于其启动压力梯度。根据气体渗流理论，可得均质气藏平面径向稳定渗流时供给边缘处的压力梯度[1—3]，代入式(5)，即可得到研究区单井有效动用半径计算公式

式中：

pe—供给边缘的地层压力，MPa；

pwmin—最小井底流动压力，MPa；

re—供给边缘半径，m；

rw—井筒半径，m。

通过式(6)计算得到研究区气层不同渗透率和含水条件下的有效动用半径，计算结果如图6所示。图6直观反映出气藏不同条件下有效动用半径特征，由图可知，对于该气藏而言：

(1)有效动用半径与气层渗透率和含水饱和度关系密切。渗透率越低，含水饱和度越高，动用半径越小。

(2)含水饱和度较低的气层，储层渗透率仅在一定范围内对井控储量有较大影响，当有效动用半径达到物理边界，渗透率对增大控制储量不再发挥影响，如含水45%时，在渗透率达到0.85 mD以后，储层有效动用半径即可达到800 m，渗透率继续增大对增大控制储量不再发挥影响。

(3)相较于渗透率，原始含水饱和度对有效动用半径影响更为显著，动用半径随原始含水饱和度增加急剧降低。含水饱和度35%时，渗透率达到0.45 mD的储层有效动用半径即可达到800 m，而当含水饱和度达到55%以上时，只有渗透率达到1.80 mD的储层其有效动用半径才能达到800 m。因此，对于低渗含水气藏，在进行井网部署，开发方案制定时，除考虑储层渗透率以外，应同时重视储层原始含水饱和度。

图6 气藏不同渗透率和含水条件下的有效动用半径Fig.6Effective development radius of gas reservoirs with different permeability and saturation

4 结论

(1)通过长岩芯多测点气藏物理模拟实验，系统测量了低渗含水砂岩气藏气相启动压力梯度，研究表明与低渗油藏储层中的渗流不同，低渗气藏气体启动压力梯度同时受储层渗透率和含水饱和度控制。

(2)单因素研究表明，启动压力梯度与含水饱和度Sw的3.5次幂正相关，与储层渗透率倒数呈线性规律。据此，分析得到了综合考虑含水饱和度和渗透率影响的启动压力梯度计算通式。对于同一区块或者岩性和孔隙结构接近的储层，通式是统一，适用的。

(3)应用气藏启动压力梯度的定量表征方法，可方便估算研究区内不同储层，不同原始含水饱和度下的启动压力梯度，确定气井有效动用半径。该方法便于实际操作，为类似气藏开发早期井网部署以及开发中后期井网加密提供了新的研究思路和计算方法，具有重要的工程实用价值。

[1]冯文光.天然气非达西低速不稳定渗流[J].天然气工业，1986，6(3)：4148.

FENG Wenguang.Non-Darcy low-velocity unsteadystateseepage flowofnatural gas[J].Natural GasIndustry，1986，6(3)：4148.

[2]闫庆来，何秋轩，尉立岗，等.低渗透油层中单相液体渗流特征的实验研究[J].西安石油学院学报，1990，5(2)：16.

YAN Qinglai，HE Qiuxuan，WEI Ligang，et al.A laboratory study on percolation characteristics of single phase flow in low-permeability reservoirs[J].Journal of Xi'an Petroleum Institute，1990，5(2)：16.

[3]杨正明，霍凌静，张亚蒲，等.含水火山岩气藏气体非线性渗流机理研究[J].天然气地球科学，2010，21(3)：371374.

YANG Zhengming，HUO Lingjing，ZHANG Yapu，et al. Mechanismofnon-linergasseepageflowinwater-bearing volcanic gas reservoir[J].Natural Gas Geoscience，2010，21(3)：371374.

[4]刘高波，刘荣和，冯文光，等.低渗透气藏非达西渗流研究[J].大庆石油地质与开发，2007，26(1)：6567.

LIU Gaobo，LIU Ronghe，FENG Wenguang，et al.Research on non-Darcy flow in low permeability gas reservoirs[J].Petroleum Geology&Oilfield Development in Daqing，2007，26(1)：6567.

[5]任晓娟，阎庆来，何秋轩，等.低渗气层气体的渗流特征实验研究[J].西安石油学院学报，1997，12(3)：2225.

REN Xiaojuan，YAN Qinglai，HE Qiuxuan，et al.The experimental study of characteristics of gas flow in tight formation[J].Journal of Xi'an Petroleum Institute，1997，12(3)：2225.

[6]贺伟，冯曦，钟孚勋.低渗储层特殊渗流机理和低渗透气井动态特征探讨[J].天然气工业，2002，22(增刊)：9194.

HE Wei，FENG Xi，ZHONG Fuxun.Discussion on the special percolation mechanism of low permeability reservoir and the performance characteristics of low permeability gas well[J].Natural Gas Industry，2002，22(S)：9194.

[7]刘晓旭，胡勇，李宁，等.低渗砂岩气藏气体特殊渗流机理实验研究与分析[J].特种油气藏，2007，14(1)：8083.

LIU Xiaoxu，HU Yong，LI Ning，et al.Experimental study and analysis of special percolation mechanism in low permeability sandstone gas reservoir[J].Special Oil and Gas Reservoirs，2007，14(1)：8083.

[8]依呷，唐海，吕栋梁.低渗气藏启动压力梯度研究与分析[J].海洋石油，2006，26(3)：5154.

YI Ga，TANG Hai，LÜ Dongliang.The study and analysis of starting pressure gradient in low permeability gas reservoirs[J].Offshore Oil，2006，26(3)：5154.

[9]邱勇松.低渗透油藏渗流机理及开发技术研究[D].北京：中国科学院研究生院，2004.

QIU Yongsong.A study on percolation mechanism of low permeability reservoir and recovery technology[D].Beijing：Chinese Academy of Science，2004.

[10]杨正明，张英芝，郝明强，等.低渗透油田储层综合评价方法[J].石油学报，2006，27(2)：6467.

YANG Zhengming，ZHANG Yingzhi，HAO Mingqiang，et al.Comprehensive evaluation of reservoir in lowpermeability oilfields[J].Acta Petrolei Sinica，2006，27(2)：6467.

[11]徐轩，杨正明，刘学伟，等.特低渗透大型露头模型流场测量技术及分布规律研究[J].岩土力学，2012，33(11)：33313337.

XU Xuan，YANG Zhengming，LIU Xuewei，et al.Study of flow field measurement technique and distribution regularity of ultra-low permeability large-scale outcrop model[J].Rock and Soil Mechanics，2012，33(11)：33313337.

[12]熊伟，雷群，刘先贵，等.低渗透油藏拟启动压力梯度[J].石油勘探与开发，2009，36(2)：232236.

XIONGWei，LEIQun，LIUXiangui，etal.Pseudothreshold pressure gradient to flow for low permeability reser-voirs[J].Petroleum Exploration and Development，2009，36(2)：232236.

[13]葛家理.油气层渗流力学[M].北京：石油工业出版社，1982.

[14]赵政璋，杜金虎，邹才能，等.致密油气[M].北京：石油工业出版社，2012.

ZHAOZhengzhang，DUJinhu，ZOUCaineng，etal.Tight oil&gas[M].Petroleum Industry Press，2012.

[15]徐轩，刘学伟，杨正明，等.特低渗透砂岩大型露头模型单相渗流特征实验[J].石油学报，2012，33(3)：453458.

XU Xuan，LIU Xuewei，YANG Zhengming，et al.An experimental study on single-phase seepage characteristics with a large-scale model made of ultra-low permeability sandstoneoutcrops[J].ActaPetroleiSinica，2012，33(3)：453458.

[16]吴柏志.低渗透油藏高效开发理论与应用[M].北京：石油工业出版社，2009.

WU Bozhi.Efficient development theory and application of the low permeability reservoir[M].Beijing：Petroleum Industry Press，2009.

[17]郝斐，程林松，李春兰，等.特低渗油藏启动压力梯度研究[J].西南石油学院学报，2006，28(6)：2932.

HAO Fei，CHENG Linsong，LI Chunlan，et al.Study on threshold pressure gradient in ultra-low permeability reservoir[J].Journal of Southwest Petroleum Institute，2006，28(6)：2932.

[18]许建红，程林松，钱俪丹，等.低渗透油藏启动压力梯度新算法及应用[J].西南石油大学学报，2007，29(4)：6466，78.

XU Jianhong，CHENG Linsong，QIAN Lidan，et al.A new method and the application on calculating start-up pressure gradient in low permeable reservoirs[J].Journal of Southwest Petroleum University，2007，29(4)：6466，78.

王继平，1978年生，男，汉族，陕西白水人，高级工程师，博士，主要从事储层地质及油气田开发方面的研究。E-mail：wangjp_cq@petrochina.com.cn

田姗姗，1984年生，女，汉族，河北保定人，工程师，硕士，主要从事地质与勘查研究工作。E-mail：cppetianss@cnpc.com.cn

焦春艳，1984年生，女，汉族，天津武清人，工程师，博士，主要从事油气藏渗流实验与开发研究工作。E-mail：jiaochunyan69@petrochina.com.cn

张宇轩，1990年生，男，汉族，新疆克拉玛依人，助理工程师，主要从事油藏动态生产管理工作。E-mail：zyxxuxuan@petrochina.com.cn

编辑：牛静静

编辑部网址：http：//zk.swpuxb.com

Non-Darcy Seapeage Law and its Application in Low Permeability Water-bearing Gas Reservoir

XU Xuan1，2*，WANG Jiping3，TIAN Shanshan4，JIAO Chunyan1，ZHANG Yuxuan5
1.Langfang Branch，Research Institute of Exploration and Development，PetroChina，Langfang，Hebei 065007，China 2.College of Geoscience and Surveying Engineering，China University of Mining&Technology，Haidian，Beijing 100083，China 3.Research Center of Sulige Gas Field，Changqing Oilfield Company，PetroChina，Xi'an，Shaanxi 710018，China 4.Pipeline Company，PetroChina，Langfang，Hebei 065000，China 5.Institute of Geology，No.1 Oil Production Plant of Xinjiang Oilfield Company，Karamay，Xinjiang 834000，China

Gas development practices and laboratory experiments show that there is starting pressure phenomenon in gas seepage in low permeability water-bearing gas reservoir.Therefore，it is very important to study the influence factors and change rules of starting pressure gradient in gas reservoir development.Through experimental methods，the change law of starting pressure in gas seepage in low permeability water-bearing gas reservoir was studied，the correlation analysis of this parameter with water saturation and reservoir permeability was carried out，and quantitative characterization method of starting pressure gradient was researched.An applied analysis was carried out on a typical reservoir with the starting pressure gradient calculation method and gas reservoir percolation theory.Effective development radius of wells was calculated，and the relationship of effective development radius with water saturation and reservoir permeability was established.This method provides a technological guidance for the similar gas reservoirs development especially on well pattern deployment in the early period and well pattern infilling in the middle and later periods.

water-bearing gas reservoir；non-Darcy flow；starting pressure gradient；experiment research；quantitative characterization

10.11885/j.issn.16745086.2014.09.09.02

16745086(2016)05009007

TE312

A

http：//www.cnki.net/kcms/detail/51.1718.TE.20160928.1448.002.html

徐轩，1984年生，男，汉族，湖北荆门人，工程师，博士，主要从事油气藏渗流实验与开发方面的研究工作。E-mail：xuxuan69@petrochina.com.cn

20140909

时间：20160928

徐轩，E-mail：xuxuan69@petrochina.com.cn

国家科技重大专项(2011ZX05013002)。