CN116338812B - 一种气藏储层含水饱和度上限的确定方法 - Google Patents

Abstract

本发明公开了一种气藏储层含水饱和度上限的确定方法，包括：对气藏储层物性下限值区域的岩心进行首次选样，并钻切成设定直径和长度的正圆柱体岩心实验标样；将岩心实验标样烘干冷却至常温；测定岩心实验标样的孔隙度和渗透率后再进行第二次选样；对第二次选样的岩心实验标样做稳定法气水相对渗透率实验并绘制气水相对渗透率曲线；对实验后岩心实验标样做水饱和度为50%条件下的气相体积流量测定，并计算出气相渗透率和气相相对渗透率；根据水饱和度为50%条件下的气相相对渗透率值在气水相对渗透率曲线上内插，得到该岩心物性构成储层的水饱和度上限值。本发明应用时，能提升确定气藏储层含水饱和度上限值的精准度。

Description

一种气藏储层含水饱和度上限的确定方法

技术领域

本发明涉及天然气地质勘探和开发技术，具体是一种气藏储层含水饱和度上限的确定方法。

背景技术

气藏储层含水饱和度上限（含气饱和度下限）是储层划分、分类、储层评价和储量计算及开发方案编制的最基础参数。目前为了获得气藏储层含水饱和度上限，普遍采用预先确定储层的孔隙度下限，再根据测井解释的目的层纵向水饱和度剖面和经过校正的孔隙度剖面建立的数学模型计算出储层水饱和度上限的方法。当前，试油方法仍是确定孔隙度下限和水饱和度上限应用最广泛、最直接的方法。然而，受钻、完井过程中的储层损害和裂缝的影响而给水饱和度上限的确定带来了不确定性，精准度较低，尤其是下限段储层，是储层中的物性最差段，受此影响尤为严重。

发明内容

本发明的目的在于解决现有技术采用优先确定孔隙度下限，再确定气藏储层含水饱和度上限精准度低的问题，提供了一种气藏储层含水饱和度上限的确定方法，其通过使用目标层低孔渗段预测的物性下限值区域的岩心，将水饱和度为50%时测定的气相渗透率所计算的气相相对渗透率在该岩心实测相对渗透率曲线上内插，从而直接确定出气藏水饱和度上限，其为确定储层水饱和度上限的一种新思路和新方法，并能提升确定气藏储层含水饱和度上限的精准度。

本发明的目的主要通过以下技术方案实现：

一种气藏储层含水饱和度上限的确定方法，包括以下步骤：

步骤S1、对被测气藏储层井下目标层物性下限值区域的岩心进行首次选样，并钻切制成设定直径和长度的正圆柱体岩心实验标样；

步骤S2、将岩心实验标样烘干冷却至常温；

步骤S3、测定岩心实验标样的孔隙度和渗透率后再进行第二次选样；

步骤S4、对第二次选样的岩心实验标样做稳定法气水相对渗透率实验并绘制气水相对渗透率曲线；

步骤S5、对实验后岩心实验标样做水饱和度为50%条件下的气相体积流量测定，并计算出气相渗透率和气相相对渗透率；

步骤S6、根据水饱和度为50%条件下所得到的气相相对渗透率值在气水相对渗透率曲线上内插，得到该岩心物性构成储层的水饱和度上限值。

进一步的，所述步骤S1中钻切制成的岩心实验标样为直径25mm、长度40mm的正圆柱体；所述步骤S2中岩心实验标样在105℃温度下烘干8小时，并在干燥器中冷却至常温。

进一步的，所述步骤S1中进行首次选样选取的数量为10-20个，所述步骤S3中进行第二次选样选取的数量为6-10个。

进一步的，所述步骤S5计算气相渗透率和气相相对渗透率包括以下步骤：

根据测定的气体体积流量计算出气相渗透率K_g50，计算公式为：

式中，K_g50水饱和度为50%时的气相渗透率，P₀为测试时大气压，Q_g50水饱和度为50%时的气体体积流量，µ_g为测试气体粘度， L为测试岩心长度， A为测试岩心横截面积， P₁为岩心进口压力， P₂为岩心出口压力；

按照下式对水饱和度为50%时的气相相对渗透率计算：

式中：K_rg50水饱和度为50%时气相相对渗透率，K_g50水饱和度为50%时气相渗透率，K为岩心绝对渗透率。

进一步的，所述步骤S6包括以下步骤：

根据水饱和度为50%时计算的气相相对渗透率，在对应的气水相对渗透率曲线纵坐标上内插气相相对渗透率计算值，做横坐标平行线与气相相对渗透率曲线的交点，然后过交点作横坐标垂线，垂线与横坐标交点的水饱和度值即为该岩心物性构成储层的水饱和度上限值。

综上所述，本发明与现有技术相比具有以下有益效果：本发明采用井下岩心在实验室测定水饱和度为50%时的气相相对渗透率与该岩心的气水相对渗透率曲线相结合以确定储层水饱和度上限值，气藏水饱和度上限值根据实验室实测资料确定，其为确定储层水饱和度上限的一种新思路和新方法，并能提升确定气藏储层含水饱和度上限值的精准度。

附图说明

此处所说明的附图用来提供对本发明实施例的进一步理解，构成本申请的一部分，并不构成对本发明实施例的限定。在附图中：

图1为本发明的一个具体实施例的流程图；

图2为为本发明的一个具体实施例，应用于川广须四段第1次选出的12个岩心的孔隙度-渗透率相关关系曲线；

图3为本发明的一个具体实施例，应用于川广须四段第2次选出的7个岩心的孔隙度-渗透率相关关系曲线；

图4为本发明的一个具体实施例，应用于川广须四段使用实测水饱和度为50%时的气相相对渗透率在该岩心的气水相对渗透率曲线上内插得到的岩心物性构成储层的水饱和度上限值的内插示意图；

图5为本发明的一个具体实施例，应用于川广须四段不同岩心取样的实测气水相对渗透率曲线。

具体实施方式

为使本发明的目的、技术方案和优点更加清楚明白，下面结合实施例和附图，对本发明作进一步的详细说明，本发明的示意性实施方式及其说明仅用于解释本发明，并不作为对本发明的限定。

实施例：

如图1所示，一种气藏储层含水饱和度上限的确定方法，包括以下步骤：步骤S1、对被测气藏储层井下目标层物性下限值区域的岩心进行首次选样，并钻切制成设定直径和长度的正圆柱体岩心实验标样；步骤S2、将岩心实验标样烘干冷却至常温；步骤S3、测定岩心实验标样的孔隙度和渗透率后再进行第二次选样；步骤S4、对第二次选样的岩心实验标样做稳定法气水相对渗透率实验并绘制气水相对渗透率曲线；步骤S5、对实验后岩心实验标样做水饱和度为50%条件下的气相体积流量测定，并计算出气相渗透率和气相相对渗透率；步骤S6、根据水饱和度为50%条件下所得到的气相相对渗透率值在气水相对渗透率曲线上内插，得到该岩心物性构成储层的水饱和度上限值。

本实施例应用时，对被测气藏储层井下目标层首次取心是根据测井资料预测（估计）的物性下限值确定的，岩心取心数量为10-20个。二次选样是在测定孔隙度和渗透率后重新选取6-10个岩心实验标样执行后续步骤。对第二次选出的6-10个岩心实验标样都测定其气水相对渗透率曲线、测定水饱和度为50%时的气相渗透率并计算其气相相对渗透率，再根据实测计算的水饱和度为50%时的气相相对渗透率，在相应岩心的气水相对渗透率曲线上内插。其中，预测的物性下限值根据测井资料并结合试油成果，在试油井段相对低孔渗段按孔隙度1%左右的间距选样。本实施例中第一次选样（取心）是在现场或岩心库根据测井资料并参考试油成果所预测的下限值左右岩心进行选样（全直径岩心），只能目测选样，所以取心数量较多10-20个，以避免通过实验后再次选样数量不够，所以可以多取，但取多了增加了工作量，少了第二次选样数量不够。将第一次选出的岩心（岩心直径大于Φ65mm的全直径岩心）用小的取心钻头（Φ25mm）重新钻取小岩心并制样。然后，按照实验室孔隙度和渗透率的测定方法和步骤以及计算方法对该岩心做孔隙度和渗透率测试，得到其被测岩心的孔隙度和渗透率并拟合孔隙度-渗透率相关曲线，剔除由于裂缝或实验过程中所产生的异常点，然后根据预测的下限值和孔-渗相关曲线，对其下限值左右的被测岩心进行第二次选样：在拟合好的曲线附近进行选样，如：预测的下限值是孔隙度为7%，则具体选样的孔隙度，低端控制在7%、6%、5%、4%、3%，而高端控制在8%、9%、10%，当然，数值并非一定是整数，也不可能完全是整数。然后，将已经选中的样品进行下一步的实验。

根据实验仪器岩心室的大小，本实施例中钻切成的岩心实验标样直径为25mm，长度为40mm，即将首次选出的岩心钻切成Φ25mm×L40mm的岩心实验标样。本实施例步骤S2中岩心实验标样在105℃温度下烘干8小时，并在干燥器中冷却至常温。其中，本实施例中岩心实验标样烘干在电热恒温鼓风干燥箱中进行，冷却在加了干燥剂的干燥器中冷却至常温，能避免放在大气中冷却，岩心吸入大气中的水蒸汽而变潮湿，达不到烘干的目的，从而影响实验结果。本实施例将烘干温度设定为105℃，是因为温度高于105℃时，会将矿物结晶水烘干（矿物结晶水是矿物组成的一部分，应该保留的），低于105℃，岩石中的自由水将烘不干，都将影响实验结果的准确性。本实施例将烘干时间设定为8小时，是经反复实验测定，在烘干温度为105℃时，烘干8小时可将岩石中的自由水全部烘干。

本实施例中孔隙度的测定采用饱和乙醇差减法测定，其包括以下步骤：

步骤11、将烘干并冷却至常温的岩心实验标样放在千分之一天平上称取岩心的干重G1；

步骤12、将已称重的岩心垂直放入玻璃杯中，掺入1/3岩心高度的无水乙醇后放入抽真空系统中，让岩心自动吸入乙醇；24小时后启动真空泵抽空1小时后停泵，让岩心继续自动吸入无水乙醇；2小时后，重新启动真空泵对其抽空1小时后停泵，然后向装岩心的玻璃杯中注入无水乙醇，让乙醇高度高于岩心2cm；4小时后再次启动真空泵继续抽空饱和乙醇，待岩心充分饱和后停止抽空并关闭电源，取出岩心称取饱和乙醇后的岩心在乙醇中的重量G2；

步骤13、从乙醇中取出岩心，用半干乙醇的纱布揩干岩心表面多余的乙醇后立刻称重G3；

步骤14、按下式计算出岩心的孔隙度：

岩心孔隙总体积：

岩心体积：

式中：Φ-被测岩心孔隙度，%；G₁-被测岩心干重，克；G₂-饱和乙醇后的岩心在乙醇中的重量，克；G₃-饱和乙醇后的岩心在空气中的重量，克；ρ_ah-无水乙醇的密度，克/cm³。

本实施例中渗透率的测定采用气体稳定流动法测定，其包括以下步骤：

步骤21、测量烘干并冷却至常温的岩心实验标样长度和直径后，再将岩心放入渗透率仪的岩心室中，开围压开关在3MPa的围压下密封岩心侧面；

步骤22、记录测试时大气压和实验温度；

步骤23、开进口氮气管路开关，使氮气在合适的压力下在岩心中流动，待流动稳定后，测定岩心在该流动压力下的气体流量Q₁，连续测定3次，取平均值并记录测试压差和气体流量；

步骤24、提高进口氮气压力，提高氮气流动压差，待流动稳定平衡后，测定岩心在新的流动压力下的气体流量Q₂，连续测定3次，取平均值同时记录新的测试压差和气体流量；

步骤25、第三次提高压差，重复步骤24的测试步骤，结束测试；

步骤26、根据测试时记录的温度查表得到的氮气粘度、记录的大气压、每一次测试压差和相应的平均气体流量计算出岩心的渗透率，其计算渗透率的公式与本实施例步骤S5计算出气相渗透率K_g50的公式相同。

本实施例中采用的稳定法气水相对渗透率实验、绘制气水相对渗透率曲线均可基于现有技术实现，其包括以下步骤：

步骤31、将上面已经测试完孔隙度和渗透率的岩心重新烘干8小时后，放入干燥器中冷却至常温后重新称干重G₄备用；

步骤32、配制模拟地层水，过滤后测定其密度后备用；

步骤33、岩心自动吸入模拟地层水，24小时后，将岩心放入真空系统中进一步抽空饱和模拟地层水，至无气泡后静置备用；

步骤34、用半渗透隔板仪将岩心压力由低至高作气水毛管压力测试，或最高压力法作一点法气水毛管压力实验，待稳定平衡后，取出岩心称重G₅，即可根据下式计算出岩心的束缚水饱和度；

式中：S_iw-岩心束缚水饱和度，%；G₄-岩心干重，克；G₅-气水毛管压力实验后的岩心重量，克；ρ_w-模拟地层水密度，克/cm³；V_p-岩心孔隙总体积，cm³。

步骤35、将岩心放入稳定法气水相对渗透率测试仪的岩心室中并密封岩心；

步骤36、测束缚水状态下的气相渗透率；

步骤37、从束缚水饱和度点开始，从水饱和度增加方向，按预先计算的每一点加水重量，采用逐点加水，稳定平衡、测试气水流量、取出称重，记录该点压力、气水流量、岩心重量、当天大气压、实验温度；加水、稳定平衡、测试气水流量、取出称重，记录该点压力、气水流量、岩心重量，重复上述过程，直至岩心水饱和度达到85%左右为止；

步骤35、将每一点记录的大气压、温度查表得到的气体和水的粘度、气水流量、实验压差分别代入下列气、水渗流的达西公式，计算出相应点的气、水相渗透率和相对渗透率。

气相渗透率：

气相相对渗透率:

水相渗透率：

相相对渗透率：

水饱和度

式中：S_Wi-实验中i点的水饱和度，%；G_i-i点的岩心重量，克；G_r-岩心的干重，克；

ρ_w-模拟地层水密度，克/cm³；V_P-实验岩心孔隙总体积。

本实施例的步骤S5计算气相渗透率和气相相对渗透率包括以下步骤：

根据测定的气体体积流量计算出气相渗透率K_g50，计算公式为：

式中，K_g50水饱和度为50%时的气相渗透率，dc，1dc=1000mdc；P₀为测试时大气压，at,1at=0.101MPa；Q_g50水饱和度为50%时的气体流量，cm³/秒；µ_g为气体粘度，mPa·s；L为测试岩心长度，cm；A为测试岩心横截面积，cm²；P₁为岩心进口压力，at,1at=0.101MPa；P₂为岩心出口压力，at,1at=0.101MPa；

按照下式对水饱和度为50%时的气相相对渗透率计算：

式中：K_rg50水饱和度为50%时的气相相对渗透率，%；K_g50水饱和度为50%时气相渗透率，dc，1dc=1000md；K为岩心绝对渗透率，dc，1dc=1000md。

本实施例步骤S5中对实验后岩心实验标样做水饱和度为50%条件下的气相体积流量测定的具体实验方法为：

在半渗透隔板仪上做气水毛管压力实验，使岩心中的水饱和度达到50%左右后结束实验。

取出岩心，迅速称其重量后将岩心装入相对渗透率仪的岩心室中，开密封压力开关密封岩心，然后根据测定相对渗透率中的方法测定岩心水饱和度为50%左右时的气体流速（体积流量）。

根据当天大气压、测试压差、温度查表获得的气体粘度、气体的体积流量代入气相渗透率K_g50计算公式，即可计算出岩心水饱和度为50%左右时的气相渗透率（K_g50）。

水饱和度为50%，也即是说，地层岩石中天然气和水的含量各占孔隙体积一半的时候，地层中气、水的渗流状态，它是一个特殊的点。由于气层岩石亲水，地层水被吸附在岩石颗粒表面和微孔微喉中，而气则分布于大孔大喉即有利于气体流动的孔隙空间中，显然，此时的水是不流动或者流动速度是很低的，测定水饱和度为50%，即地层含气饱和度为50%的情况下，气体的流动速度，可进一步评价地层的产气情况；根据国内外的经验，在没有含气饱和度下限值的情况下，常常将50%（水饱和度50%）作为含气饱和度的下限值。所以，本实验选择水饱和度为50%左右的条件来测定岩石的气相渗透率。同时也可以检验该区该类储层的水饱和度上限是否就是50%了。

如图4所示，本实施例步骤S6包括以下步骤：根据水饱和度为50%时计算的气相相对渗透率，在对应的气水相对渗透率曲线纵坐标上内插气相相对渗透率计算值，做横坐标平行线与气相相对渗透率曲线的交点，然后过交点作横坐标垂线，垂线与横坐标交点的水饱和度值即为该岩心物性构成储层的水饱和度上限值。

本实施例提供了岩心从取心处理至气水相对渗透率曲线、水饱和度为50%时的气相渗透率的实验室测定全流程和方法，使用水饱和度为50%时气相相对渗透率在相应岩心的气水相对渗透率曲线图上内插确定储层水饱和度上限值，具有坚实的理论基础，原理清楚，操作简单，结果可靠。

川广气田位于四川省中部广、南市境内，区域上位于川中古隆中斜平缓构造带南充构造群东部。目的层为三叠系上统须家河组须四段，岩性为浅灰色中粒岩屑砂岩、长石岩屑砂岩为主间夹薄层泥岩及煤线。其主要储集空间为残余粒间孔、粒间溶孔、粒内溶孔和胶结物及杂基溶孔。孔隙度主要分布在1%-15%之间。砂岩孔隙度主要分布在2%-12%之间，占84%，其中2%-8%的岩心占63%。渗透率在0.01md-5md之间，储层渗透率主要集中在0.01md-0.2md之间占67%。目的层井深在2300米-2500米之间。地层压力30MPa、地层温度80℃。从2005年-2006年，先后钻获20余口气井，测试获气(0.8-3.0)×10⁸m³。川广气田须四段为川广气田其中一段，本实施例应用于川广气田须四段时，采用本实施例的方法：

在岩心库根据测井资料在相对低孔、渗段首次选出12个Φ65mm的全直径岩心，然后，在岩心制样室用Φ25mm的取心钻头沿全直径岩心平行层面方向钻取Φ25mm小岩心12个，然后在切片机上切制成长度为42mm的圆柱体，最后在磨床上将岩心磨制成Φ25mm，长度为40mm的正圆柱体。其中，本实施例应用时对须四段，低端孔隙度小于3%和高端孔隙度大于10%的岩心除外，在孔隙度3%-10%之间的岩心按孔隙度大约1%的间距选样。

12个岩心经过105℃恒温烘干8小时后放入干燥器中冷却至常温，然后对12个岩心测定其孔隙度和渗透率（见表1）。

根据测井资料预测的须四段孔隙度下限为5%-6%，首次选择的12个岩心的孔隙度在3.52%-9.71%之间，渗透率在0.0031 md-0.5771md之间，拟合的孔隙度-渗透率相关曲线见图2。

按照孔隙度-渗透率相关曲线附近选样的原则，选择其中7个岩心进行后续实验，7个岩心孔隙度在4.22%-7.23%之间，渗透率在0.0217md-0.1610md之间，拟合的孔隙度-渗透率相关曲线见图3。

对比图2和图3表明：二次选样的相关系数高于首次选样的相关系数，选样是成功的，说明二次选样可以用作后续实验。

表1 川广须四段1、2次选样实测孔隙度、渗透率及水饱和度为50%时气相渗透率实测值和内插的水饱和度上限值表

对二次选中的7个岩心作气水相对渗透率实验所作出的相对渗透率曲线见图5。

7个对应岩心水饱和度为50%时的气相渗透率(K_g50)在0.0038md-0.0254md，所计算的气相相对渗透率(K_rg50=K_g50/K)在11.1%-19.4%，平均值是16.5%，在相应岩心气水相对渗透率曲线上内插的水饱和度上限值（S_wmax）在53.6%-58.7%之间，其平均值为56.4%。这就是说，根据本实施例所确定的川广须四段水饱和度上限值是56.4%，如表1所示。

为了使用方便、储量计算的稳妥和水饱和度就低不就高（气饱和度就高不就低）的原则，将川广须四段水饱和度上限值确定为55%，代入储量计算中，所得结论符合生产实际，符合储量规范的要求，所确定的有效储层物性下限参数已应用于川广气田须四段的探明储量计算。

按照水饱和度上限55%（气饱和度下限45%）的研究成果，探明储量670×10⁸m³天然气。

以上所述的具体实施方式，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施方式而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (5)

1.一种气藏储层含水饱和度上限的确定方法，其特征在于，包括以下步骤：

步骤S1、对被测气藏储层井下目标层物性下限值区域的岩心进行首次选样，并钻切制成设定直径和长度的正圆柱体岩心实验标样；

步骤S2、将岩心实验标样烘干冷却至常温；

步骤S3、测定岩心实验标样的孔隙度和渗透率后再进行第二次选样；

步骤S4、对第二次选样的岩心实验标样做稳定法气水相对渗透率实验并绘制气水相对渗透率曲线；

步骤S5、对实验后岩心实验标样做水饱和度为50%条件下的气相体积流量测定，并计算出气相渗透率和气相相对渗透率；

步骤S6、根据水饱和度为50%条件下所得到的气相相对渗透率值在气水相对渗透率曲线上内插，得到该岩心物性构成储层的水饱和度上限值；

所述孔隙度的测定采用饱和乙醇差减法测定，其包括以下步骤：

步骤11、将烘干并冷却至常温的岩心实验标样放在千分之一天平上称取岩心的干重G1；

步骤12、将已称重的岩心垂直放入玻璃杯中，掺入1/3岩心高度的无水乙醇后放入抽真空系统中，让岩心自动吸入乙醇；24小时后启动真空泵抽空1小时后停泵，让岩心继续自动吸入无水乙醇；2小时后，重新启动真空泵对其抽空1小时后停泵，然后向装岩心的玻璃杯中注入无水乙醇，让乙醇高度高于岩心2cm；4小时后再次启动真空泵继续抽空饱和乙醇，待岩心充分饱和后停止抽空并关闭电源，取出岩心称取饱和乙醇后的岩心在乙醇中的重量G2；

步骤13、从乙醇中取出岩心，用半干乙醇的纱布揩干岩心表面多余的乙醇后立刻称重G3；

步骤14、按下式计算出岩心的孔隙度：

岩心孔隙总体积：

岩心体积：

式中：Φ-被测岩心孔隙度，%；G₁-被测岩心干重，克；G₂-饱和乙醇后的岩心在乙醇中的重量，克；G₃-饱和乙醇后的岩心在空气中的重量，克；ρ_ah-无水乙醇的密度，克/cm³；

所述渗透率的测定采用气体稳定流动法测定，其包括以下步骤：

步骤21、测量烘干并冷却至常温的岩心实验标样长度和直径后，再将岩心放入渗透率仪的岩心室中，开围压开关在3MPa的围压下密封岩心侧面；

步骤22、记录测试时大气压和实验温度；

步骤23、开进口氮气管路开关，使氮气在合适的压力下在岩心中流动，待流动稳定后，测定岩心在该流动压力下的气体流量Q₁，连续测定3次，取平均值并记录测试压差和气体流量；

步骤24、提高进口氮气压力，提高氮气流动压差，待流动稳定平衡后，测定岩心在新的流动压力下的气体流量Q₂，连续测定3次，取平均值同时记录新的测试压差和气体流量；

步骤25、第三次提高压差，重复步骤24的测试步骤，结束测试；

步骤26、根据测试时记录的温度查表得到的氮气粘度、记录的大气压、每一次测试压差和相应的平均气体流量计算出岩心的渗透率。

2.根据权利要求1所述的一种气藏储层含水饱和度上限的确定方法，其特征在于，所述步骤S1中钻切制成的岩心实验标样为直径25mm、长度40mm的正圆柱体；所述步骤S2中岩心实验标样在105℃温度下烘干8小时，并在干燥器中冷却至常温。

3.根据权利要求1所述的一种气藏储层含水饱和度上限的确定方法，其特征在于，所述步骤S1中进行首次选样选取的数量为10-20个，所述步骤S3中进行第二次选样选取的数量为6-10个。

4.根据权利要求1所述的一种气藏储层含水饱和度上限的确定方法，其特征在于，所述步骤S5计算气相渗透率和气相相对渗透率包括以下步骤：

根据测定的气体体积流量计算出气相渗透率K_g50，计算公式为：

；式中，K_g50水饱和度为50%时的气相渗透率，P₀为测试时大气压，Q_g50水饱和度为50%时的气体体积流量，µ_g为测试气体粘度， L为测试岩心长度， A为测试岩心横截面积， P₁为岩心进口压力， P₂为岩心出口压力；

按照下式对水饱和度为50%时的气相相对渗透率计算：

；式中：K_rg50水饱和度为50%时气相相对渗透率，K_g50水饱和度为50%时气相渗透率，K为岩心绝对渗透率。

5.根据权利要求1～4中任意一项所述的一种气藏储层含水饱和度上限的确定方法，其特征在于，所述步骤S6包括以下步骤：

根据水饱和度为50%时计算的气相相对渗透率，在对应的气水相对渗透率曲线纵坐标上内插气相相对渗透率计算值，做横坐标平行线与气相相对渗透率曲线的交点，然后过交点作横坐标垂线，垂线与横坐标交点的水饱和度值即为该岩心物性构成储层的水饱和度上限值。