(19)中华 人民共和国 国家知识产权局 (12)发明 专利申请 (10)申请公布号 (43)申请公布日 (21)申请 号 202111433441.7 (22)申请日 2021.11.29 (71)申请人 长江大学 地址 430100 湖北省武汉市蔡甸区大 学路 111号 (72)发明人 许诗婧　王长权　 (74)专利代理 机构 北京圣州专利代理事务所 (普通合伙) 11818 代理人 黄青青 (51)Int.Cl. G06F 30/28(2020.01) G06F 17/11(2006.01) G06F 113/08(2020.01) G06F 119/14(2020.01) (54)发明名称 一种致密砂岩气藏气体相对渗透率计算方 法 (57)摘要 本发明公开了一种致密砂岩气藏气体相对 渗透率计算方法， 包括如下过程： 根据气体相对 渗透率值与残 余水状态下的气相渗透率值， 计算 不同含气饱和度下的气相渗透率； 根据含气饱和 度， 计算不同含气饱和度下的孔隙度； 根据渗透 率与孔隙半径的关系， 计算不同含气饱和度下的 储层品质指数； 根据岩心初始液测渗透率与储层 品质指数的关系， 计算不同含气饱和度下的液体 渗透率， 即将原气体渗透率转化为同样含气饱和 度下的液体渗透率； 修正气体相对渗透率， 并绘 制气体相对渗透率曲线。 本发明采用上述计算方 法， 考虑气体滑脱的影响， 将气相渗透率转化为 不同含气饱和度下的液相渗透率进行相渗计算， 可以更准确的掌握致密砂岩气藏气水相渗变化 规律， 为气藏 开发提供依据。 权利要求书1页 说明书7页 附图2页 CN 114139470 A 2022.03.04 CN 114139470 A 1.一种致密 砂岩气藏气体相对渗透率计算方法， 其特 征在于， 包括如下步骤： S1、 根据气体相对渗透率值与残余水状态下的气相渗透率值， 计算不同含气饱和度下 的气相渗透率； 其中， 气体相 对渗透率值依据JBN公式计算而得， 残余水状态下的气相渗透 率值依据气测渗透率的计算公式计算而得； S2、 根据含气饱和度， 计算 不同含气饱和度下的孔隙度； S3、 根据渗透率与孔隙半径的关系， 计算 不同含气饱和度下的储层品质指数； S4、 根据岩心初始液测渗透率与储层品质指数的关系， 计算当前地区岩心不同含气饱 和度下的液体渗透率， 即将原气体渗透率 转化为同样含气饱和度下的液体渗透率； S5、 修正气体相对渗透率， 绘制气体相对渗透率曲线。 2.根据权利要求1所述的致密气藏气体相对渗透率计算方法， 其特征在于， S1步骤中， 不同含气饱和度下的气相渗透率的计算公式如下： 其中， Krgi为不同含气饱和度下的气体相对渗透率， Krg(swi)为束缚水饱和度下的气体相 对渗透率， Kg(swi)为束缚水饱和度下的气相渗透率。 3.根据权利要求1所述的致密气藏气体相对渗透率计算方法， 其特征在于， S2步骤中， 不同含气饱和度下的孔隙度的计算公式如下： φgi＝φg×Sgi×100 其中， φg为岩心孔隙度， Sgi为不同的含气饱和度。 4.根据权利要求3所述的致密气藏气体相对渗透率计算方法， 其特征在于， S3步骤中， 不同含气饱和度下的储层品质指数的计算公式如下： 其中， Kgi为不同含气饱和度下的气相渗透率， φgi为不同含气饱和度下的孔隙度。 5.根据权利要求4所述的致密气藏气体相对渗透率计算方法， 其特征在于， S4步骤中， 岩心初始液测渗透率为当前地区岩心在完全饱和水后测得的初始水相渗透率， 岩心初始液 测渗透率与储层品质指数的关系如下： Kw＝aRQIb 其中， RQI为储层品质指数， a、 b为拟合 参数。 6.根据权利要求5所述的致密气藏气体相对渗透率计算方法， 其特征在于， S5步骤中， 气体相对渗透率的修 正公式如下： 其中， Kwi为不同含气饱和度下的液测渗透率， Kw为岩心初始液测渗透率。权　利　要　求　书 1/1 页 2 CN 114139470 A 2一种致密砂岩气藏气体相对渗透率计算方 法 技术领域 [0001]本发明涉及油气藏开发技术领域， 尤其是涉及一种致密砂岩气藏气体相对渗透率 计算方法。 背景技术 [0002]相对渗透率通常被用来表征储层的孔 隙结构、 饱和度分布和润湿性， 气水相渗可 以表征储层中气水的相对流动能力。 因此， 准确获得不同条件 下气水相对渗透率， 对于更好 地了解储层中气体流动的内在规 律十分重要。 [0003]现有气体相对渗透率测量方法中， 由于未充分考虑滑脱效应的影响， 气体相对渗 透率会被高估， 而且最终数据也未进行 校正， 容易产生较大的误差 。 发明内容 [0004]为解决上述 技术问题， 本发明提供了如下技 术方案： [0005]一种致密 砂岩气藏气体相对渗透率计算方法， 包括如下步骤： [0006]S1、 根据气体相对渗透率值与残余水状态下的气相渗透率值， 计算不同含气饱和 度下的气相渗透率； 其中， 气体相对渗透率值依据JBN公式计算而得， 残余水状态下的气相 渗透率值依据气测渗透率的计算公式计算而得； [0007]S2、 根据含气饱和度， 计算 不同含气饱和度下的孔隙度； [0008]S3、 根据渗透率与孔隙半径的关系计算 不同含气饱和度下的储层品质指数； [0009]S4、 根据岩心初始液测渗透率与储层品质指数的关系， 计算当前地区岩心不 同含 气饱和度下的液体渗透率， 即将原气体渗透率 转化为同样含气饱和度下的液体渗透率； [0010]S5、 修正气体相对渗透率， 绘制气体相对渗透率曲线。 [0011]优选的， S1步骤中， 不同含气饱和度下的气相渗透率的计算公式如下： [0012] [0013]其中， Krgi为不同含气饱和度下的气体相对渗透率， Krg(swi)为束缚水饱和度下的气 体相对渗透率， Kg(swi)为束缚水饱和度下的气相渗透率。 [0014]优选的， S2步骤中， 不同含气饱和度下的孔隙度的计算公式如下： [0015]φgi＝φg×Sgi×100 [0016]其中， φg为岩心孔隙度， Sgi为不同的含气饱和度。 [0017]优选的， S3步骤中， 不同含气饱和度下的储层品质指数的计算公式如下： [0018] [0019]其中， Kgi为不同含气饱和度下的气相渗透率， φgi为不同含气饱和度下的孔隙度。 [0020]优选的， S4步骤中， 岩心初始液测渗透率为当前地区岩心在完全饱和水后测得的说　明　书 1/7 页 3 CN 114139470 A 3