CN107239648B - 页岩气井产量构成确定方法及装置 - Google Patents

页岩气井产量构成确定方法及装置 Download PDF

Info

Publication number
CN107239648B
CN107239648B CN201610180381.5A CN201610180381A CN107239648B CN 107239648 B CN107239648 B CN 107239648B CN 201610180381 A CN201610180381 A CN 201610180381A CN 107239648 B CN107239648 B CN 107239648B
Authority
CN
China
Prior art keywords
yield
dimensionless
pressure
formation pressure
gas well
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201610180381.5A
Other languages
English (en)
Other versions
CN107239648A (zh
Inventor
艾爽
张同义
毛军
杜娟
吴琼
刘欢乐
邸德家
贺英
庞伟
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Petroleum and Chemical Corp
Sinopec Research Institute of Petroleum Engineering
Original Assignee
China Petroleum and Chemical Corp
Sinopec Research Institute of Petroleum Engineering
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China Petroleum and Chemical Corp, Sinopec Research Institute of Petroleum Engineering filed Critical China Petroleum and Chemical Corp
Priority to CN201610180381.5A priority Critical patent/CN107239648B/zh
Publication of CN107239648A publication Critical patent/CN107239648A/zh
Application granted granted Critical
Publication of CN107239648B publication Critical patent/CN107239648B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
    • G16Z99/00—Subject matter not provided for in other main groups of this subclass

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Production Of Liquid Hydrocarbon Mixture For Refining Petroleum (AREA)

Abstract

本发明提供一种页岩气井产量构成确定方法及装置，其中方法包括：步骤100，设置当前无因次拟时间下的平均地层压力
Figure DDA0000950734180000011
使得
Figure DDA0000950734180000012
与原始地层压力相同；步骤101，根据无因次产能公式，获取拉式空间下的无因次产量
Figure DDA0000950734180000013
步骤102，根据
Figure DDA0000950734180000014
利用数值反演法计算获得当前拟时间下的无因次产量qD；步骤103，根据qD，计算获得累积产量Gp；步骤104，根据Gp，计算获得平均地层压力
Figure DDA0000950734180000015
步骤105，判断
Figure DDA0000950734180000016
是否成立，σ为预设阈值，若成立，则执行步骤106，若不成立，则使
Figure DDA0000950734180000017
转步骤103执行；步骤106，根据
Figure DDA0000950734180000018
计算获得页岩气井的产量构成。通过本发明中的方法计算出的页岩气井构成情况并非固定值，而是随时间发生变化，计算出的页岩气井产量构成与实际情况更相符。

Description

页岩气井产量构成确定方法及装置
技术领域
本发明涉及气藏开发技术领域，尤其涉及一种页岩气井产量构成确定方法及装置。
背景技术
页岩气井中生产的天然气有两种来源：吸附气和游离气。解吸压缩系数可用于表征吸附气的供给作用，解吸压缩系数占总压缩系数的比重可看做吸附气天然气总气量的比重。现在普遍认为解吸压缩系数是一个常数，由此计算出的页岩气井的产量构成，吸附气和游离气占总产气量的比重都是固定的。而实际上由于页岩气藏普遍衰竭式开发，在开发过程中地层压力不断下降，此时解吸压缩系数也随地层压力的变化而变化，因此忽略解吸压缩系数随地层压力的变化所计算出来的页岩气井产量构成是不合理的。因此，亟需一种考虑解吸压缩系数随地层压力变化而变化的页岩气井产量构成确定方法。
发明内容
本发明提供一种页岩气井产量构成确定方法及装置，用以解决现有技术中计算出来的页岩气井产量构成不合理的技术问题。
本发明一方面提供一种页岩气井产量构成确定方法，包括：
步骤100，设置当前无因次拟时间下的平均地层压力
Figure BDA0000950734160000011
使得当前无因次拟时间下的平均地层压力
Figure BDA0000950734160000012
与原始地层压力相同；
步骤101，根据无因次产能公式，获取拉式空间下的无因次产量
Figure BDA0000950734160000013
步骤102，根据拉式空间下的无因次产量
Figure BDA0000950734160000014
利用数值反演法计算获得当前拟时间下的无因次产量qD；
步骤103，根据当前拟时间下的无因次产量qD，计算获得累积产量Gp；
步骤104，根据累积产量Gp，计算获得平均地层压力
Figure BDA0000950734160000015
步骤105，判断
Figure BDA0000950734160000016
是否成立，σ为预设阈值，若成立，则执行步骤106，若不成立，则使
Figure BDA0000950734160000021
转步骤103执行；
步骤106，根据平均地层压力
Figure BDA0000950734160000022
计算获得页岩气井的产量构成。
进一步的，步骤103具体包括：
计算当前无因次拟时间下的平均地层压力
Figure BDA0000950734160000023
的气体物性参数，所述气体物性参数包括气体压缩因子zg、气体密度ρg、气体粘度μg、气体压缩系数cg；
根据气体物性参数，获取拟时间ta；
根据拟时间，获得真实时间t；
根据当前拟时间下的无因次产量qD，获得真实产量qsc；
根据真实产量qsc，计算获得真实时间t内的累积产量Gp。
进一步的，步骤106具体包括：
Figure BDA0000950734160000024
其中，
Figure BDA0000950734160000025
为解吸压缩系数，该参数是平均地层压力
Figure BDA0000950734160000026
的函数，
Figure BDA0000950734160000027
psc为标况压力，T为储层温度，
Figure BDA0000950734160000028
为基质孔隙度，Tsc为标况温度，VL为Langmuir体积，pL为Langmuir压力，cpm为基质孔隙压缩系数，
Figure BDA0000950734160000029
表示压力为平均地层压力
Figure BDA00009507341600000210
时的气体压缩因子。
进一步的，步骤104具体包括：
平均地层压力
Figure BDA00009507341600000211
根据如下公式确定：
Figure BDA00009507341600000212
其中，pi为原始地层压力，G为气井控制储量，zg(pi)表示压力为原始地层压力pi时的气体压缩因子。
进一步的，根据当前拟时间下的无因次产量qD，获得真实产量qsc具体包括：
Figure BDA00009507341600000213
其中，kF为裂缝渗透率，Acw为泄流面积，ψi为原始地层条件下的拟压力，ψwf为井底压力条件下的拟压力，T为储层温度。
本发明另一方面提供一种页岩气井产量构成确定装置，包括：
初始化模块，用于设置当前无因次拟时间下的平均地层压力
Figure BDA0000950734160000031
使得当前无因次拟时间下的平均地层压力
Figure BDA0000950734160000032
与原始地层压力相同；
拉式空间下的无因次产量获取模块，用于根据无因次产能公式，获取拉式空间下的无因次产量
Figure BDA0000950734160000033
当前拟时间下的无因次产量获取模块，用于根据拉式空间下的无因次产量
Figure BDA0000950734160000034
利用数值反演法计算获得当前拟时间下的无因次产量qD；
累积产量获取模块，用于根据当前拟时间下的无因次产量qD，计算获得累积产量Gp；
平均地层压力获取模块，用于根据累积产量Gp，计算获得平均地层压力
Figure BDA0000950734160000035
判断模块，用于判断
Figure BDA0000950734160000036
是否成立，σ为预设阈值，若成立，则触发页岩气井产量构成计算模块，若不成立，则使
Figure BDA0000950734160000037
触发累积产量获取模块；
页岩气井产量构成计算模块，用于根据平均地层压力
Figure BDA0000950734160000038
计算获得页岩气井的产量构成。
进一步的，累积产量获取模块具体用于：
计算当前无因次拟时间下的平均地层压力
Figure BDA0000950734160000039
的气体物性参数，所述气体物性参数包括气体压缩因子zg、气体密度ρg、气体粘度μg、气体压缩系数cg；
根据气体物性参数，获取拟时间ta；
根据拟时间，获得真实时间t；
根据当前拟时间下的无因次产量qD，获得真实产量qsc；
根据真实产量qsc，计算获得真实时间t内的累积产量Gp。
进一步的，页岩气井产量构成计算模块具体用于：
Figure BDA00009507341600000310
其中，
Figure BDA00009507341600000311
为解吸压缩系数，该参数是平均地层压力
Figure BDA00009507341600000312
的函数，
Figure BDA0000950734160000041
psc为标况压力，T为储层温度，
Figure BDA0000950734160000042
为基质孔隙度，Tsc为标况温度，VL为Langmuir体积，pL为Langmuir压力，cpm为基质孔隙压缩系数，
Figure BDA0000950734160000043
表示压力为平均地层压力
Figure BDA0000950734160000044
时的气体压缩因子。
进一步的，平均地层压力获取模块具体用于：
平均地层压力
Figure BDA0000950734160000045
根据如下公式确定：
Figure BDA0000950734160000046
其中，pi为原始地层压力，G为气井控制储量，zg(pi)表示压力为原始地层压力pi时的气体压缩因子。
进一步的，根据当前拟时间下的无因次产量qD，获得真实产量qsc具体用于：
Figure BDA0000950734160000047
其中，kF为裂缝渗透率，Acw为泄流面积，ψi为原始地层条件下的拟压力，ψwf为井底压力条件下的拟压力，T为储层温度。
本发明提供的页岩气井产量构成确定方法及装置，通过根据当前拟时间下的无因次产量，计算获得累积产量；然后根据累积产量，计算获得平均地层压力；最后通过计算获得平均地层压力来确定页岩气井的产量构成，即页岩气井的吸附气构成和游离气构成，从而使计算出的页岩气井构成情况并非固定值，而是随时间发生变化，通过本发明中的方法计算出的页岩气井产量与实际情况更相符。
附图说明
在下文中将基于实施例并参考附图来对本发明进行更详细的描述。其中：
图1为根据本发明实施例一的页岩气井产量构成确定方法的流程示意图；
图2为根据本发明实施例二的页岩气井产量构成确定方法的流程示意图；
图3为根据本发明实施例三的页岩气井产量构成确定装置的流程示意图。
在附图中，相同的部件使用相同的附图标记。附图并未按照实际的比例绘制。
具体实施方式
下面将结合附图对本发明作进一步说明。
实施例一
图1为根据本发明实施例一的页岩气井产量构成确定方法的流程示意图；如图1所示，本实施例提供一种页岩气井产量构成确定方法，包括：
步骤100，设置当前无因次拟时间下的平均地层压力
Figure BDA0000950734160000051
使得当前无因次拟时间下的平均地层压力
Figure BDA0000950734160000052
与原始地层压力相同。
具体的，当前无因次拟时间可在0.001～1000000中任意取值，可根据实际需要选取；平均地层压力指生产进行到某一时间时地层的平均压力，原始地层压力是指开采前的地层压力，原始地层压力的值可通过查阅地质资料或者试井解释报告获得。
步骤101，根据无因次产能公式，获取拉式空间下的无因次产量
Figure BDA0000950734160000053
具体的，拉式空间是一种虚拟的空间，先把无因次产能公式转换到拉式空间，然后后续步骤中利用数值反演法反演到实际空间，就可以求出无因次产量。
无因次产量
Figure BDA0000950734160000054
可根据如下公式获取：
Figure BDA0000950734160000055
其中，yeD为无因次缝长，ω为弹性储容比，λcw为窜流系数，s为拉普拉斯算子。
步骤102，根据拉式空间下的无因次产量
Figure BDA0000950734160000056
利用数值反演法计算获得当前拟时间下的无因次产量qD。
具体的，数值反演法包括stefest算法，可利用stefest算法计算获得当前拟时间下的无因次产量qD，“当前拟时间”是针对“当前无因次拟时间”而言的，把当前无因次拟时间变换为有因次，则为当前拟时间。
步骤103，根据当前拟时间下的无因次产量qD，计算获得累积产量Gp。
步骤104，根据累积产量Gp，计算获得平均地层压力
Figure BDA0000950734160000057
步骤105，判断
Figure BDA0000950734160000058
是否成立，σ为预设阈值，若成立，则执行步骤106，若不成立，则使
Figure BDA0000950734160000061
转步骤103执行。
具体的，σ的取值范围为0到0.005之间，优选取0.001。当前拟时间下的平均地层压力
Figure BDA0000950734160000062
与
Figure BDA0000950734160000063
之间差值小于σ时，转步骤106继续执行，否则，需要再次转步骤103，重新计算
Figure BDA0000950734160000064
的值。
步骤106，根据平均地层压力
Figure BDA0000950734160000065
计算获得页岩气井的产量构成。
本实施例提供的页岩气井产量构成确定方法，通过根据当前拟时间下的无因次产量，计算获得累积产量；然后根据累积产量，计算获得平均地层压力；最后通过计算获得平均地层压力来确定页岩气井的产量构成，即页岩气井的吸附气构成和游离气构成，从而使计算出的页岩气井构成情况并非固定值，而是随时间发生变化。
实施例二
本实施例是在上述实施例的基础上进行的补充说明。
图2为根据本发明实施例二的页岩气井产量构成确定方法的流程示意图；如图2所示，本实施例提供一种页岩气井产量构成确定方法，包括：
步骤200，设置当前无因次拟时间下的平均地层压力
Figure BDA0000950734160000066
使得当前无因次拟时间下的平均地层压力
Figure BDA0000950734160000067
与原始地层压力相同。
步骤201，根据无因次产能公式，获取拉式空间下的无因次产量
Figure BDA0000950734160000068
步骤202，根据拉式空间下的无因次产量
Figure BDA0000950734160000069
利用数值反演法计算获得当前拟时间下的无因次产量qD。
步骤200-步骤202与实施例一中的步骤100-步骤102相同，具体可参见实施例一中相应的描述，在此不再赘述。
步骤2031，计算当前无因次拟时间下的平均地层压力
Figure BDA00009507341600000610
的气体物性参数，所述气体物性参数包括气体压缩因子zg、气体密度ρg、气体粘度μg、气体压缩系数cg。
具体的，气体压缩因子zg可通过如下公式进行计算获得：
Figure BDA0000950734160000071
其中，Tpr为对比温度。
气体密度ρg的确定方法可根据如下公式获得：
Figure BDA0000950734160000072
其中，γg表示气体相对密度，T为储层温度，当
Figure BDA0000950734160000073
的值还未确定时，
Figure BDA0000950734160000074
气体粘度μg的确定方法可根据如下公式获得：
Figure BDA0000950734160000075
其中，Mg为气体分子量，T为储层温度。
气体压缩系数cg的确定方法可根据如下公式获得：
Figure BDA0000950734160000076
步骤2032，根据气体物性参数，获取拟时间ta。
具体的，拟时间ta可根据如下公式确定：
Figure BDA0000950734160000077
其中，pi为原始地层压力，cpm为基质孔隙压缩系数；Acw为泄流面积；tDa为当前无因次拟时间；cpF为裂缝压缩系数，kF为裂缝渗透率，φm为基质孔隙度，φF为裂缝孔隙度。
Figure BDA0000950734160000081
Figure BDA0000950734160000082
上式中
Figure BDA0000950734160000083
表示对zg(pi)中的pi求导数，T代表储层温度，表明该气体压缩系数是在等温条件下获得的。
步骤2033，根据拟时间ta，获得真实时间t。
真实时间的确定方法：
Figure BDA0000950734160000084
其中，
Figure BDA0000950734160000085
Figure BDA0000950734160000086
步骤2034，根据当前拟时间下的无因次产量qD，获得真实产量qsc。
Figure BDA0000950734160000087
其中，kF为裂缝渗透率，Acw为泄流面积，ψi为原始地层条件下的拟压力，ψwf为井底压力条件下的拟压力，T为储层温度。
步骤2035，根据真实产量qsc，计算获得真实时间t内的累积产量Gp。
具体的，累积产量Gp的确定方法为：
Figure BDA0000950734160000088
步骤204，根据累积产量Gp，计算获得平均地层压力
Figure BDA0000950734160000089
平均地层压力
Figure BDA00009507341600000810
可根据如下公式确定：
Figure BDA0000950734160000091
其中，pi为原始地层压力，G为气井控制储量，zg是压力的函数，zg(pi)表示压力为原始地层压力pi时的气体压缩因子zg。
步骤205，判断
Figure BDA0000950734160000092
是否成立，σ为预设阈值，若成立，则执行步骤206，若不成立，则使
Figure BDA0000950734160000093
转步骤2031执行。
步骤206，根据平均地层压力
Figure BDA0000950734160000094
计算获得页岩气藏产井的产量构成。
Figure BDA0000950734160000095
其中，
Figure BDA0000950734160000096
为解吸压缩系数，该参数是平均地层压力
Figure BDA0000950734160000097
的函数，
Figure BDA0000950734160000098
psc为标况压力，T为储层温度，
Figure BDA0000950734160000099
为基质孔隙度，Tsc为标况温度，VL为Langmuir体积，pL为Langmuir压力，cpm为基质孔隙压缩系数，
Figure BDA00009507341600000910
表示压力为平均地层压力
Figure BDA00009507341600000911
时的气体压缩因子。
本实施例提供的页岩气井产量构成确定方法，通过根据当前拟时间下的无因次产量，计算获得累积产量；然后根据累积产量，计算获得平均地层压力；最后通过计算获得平均地层压力来确定页岩气井的产量构成，即页岩气井的吸附气构成和游离气构成，从而使计算出的页岩气井构成情况并非固定值，而是随时间发生变化。本实施例中的求取吸附气构成和游离气构成的方法，考虑了解吸压缩系数随地层压力变化而产生的变化，从而使计算出的页岩气井构成情况随时间发生变化。
实施例三
本实施例为装置实施例，用于执行上述实施例一中的方法。
图3为根据本发明实施例三的页岩气井产量构成确定装置的结构示意图；如图3所示，本实施例提供一种页岩气井产量构成确定装置，包括初始化模块301、拉式空间下的无因次产量获取模块302、当前拟时间下的无因次产量获取模块303、累积产量获取模块304、平均地层压力获取模块305、判断模块306和页岩气井产量构成计算模块307。
其中，初始化模块301，用于设置当前无因次拟时间下的平均地层压力
Figure BDA0000950734160000101
使得当前无因次拟时间下的平均地层压力
Figure BDA0000950734160000102
与原始地层压力相同；
拉式空间下的无因次产量获取模块302，用于根据无因次产能公式，获取拉式空间下的无因次产量
Figure BDA0000950734160000103
当前拟时间下的无因次产量获取模块303，用于根据拉式空间下的无因次产量
Figure BDA0000950734160000104
利用数值反演法计算获得当前拟时间下的无因次产量qD；
累积产量获取模块304，用于根据当前拟时间下的无因次产量qD，计算获得累积产量Gp；
平均地层压力获取模块305，用于根据累积产量Gp，计算获得平均地层压力
Figure BDA0000950734160000105
判断模块306，用于判断
Figure BDA0000950734160000106
是否成立，σ为预设阈值，若成立，则触发页岩气井产量构成计算模块307，若不成立，则使
Figure BDA0000950734160000107
触发累积产量获取模块304；
页岩气井产量构成计算模块307，用于根据平均地层压力
Figure BDA0000950734160000108
计算获得页岩气井的产量构成。
本实施例是与方法实施例一对应的装置实施例，具体可参见实施例一中的描述，在此不再赘述。
实施例四
本实施例是在实施例三的基础上进行的补充说明，用于执行上述实施例二中的方法。
本实施例提供一种页岩气井产量构成确定装置，包括初始化模块301、拉式空间下的无因次产量获取模块302、当前拟时间下的无因次产量获取模块303、累积产量获取模块304、平均地层压力获取模块305、判断模块306和页岩气井产量构成计算模块307。
其中，累积产量获取模块304具体用于：
计算当前无因次拟时间下的平均地层压力
Figure BDA0000950734160000109
的气体物性参数，所述气体物性参数包括气体压缩因子zg、气体密度ρg、气体粘度μg、气体压缩系数cg；
根据气体物性参数，获取拟时间ta；
根据拟时间，获得真实时间t；
根据当前拟时间下的无因次产量qD，获得真实产量qsc；
根据真实产量qsc，计算获得真实时间t内的累积产量Gp。
页岩气井产量构成计算模块307具体用于：
Figure BDA0000950734160000111
其中，
Figure BDA0000950734160000112
为解吸压缩系数，该参数是平均地层压力
Figure BDA0000950734160000113
的函数，
Figure BDA0000950734160000114
psc为标况压力，T为储层温度，
Figure BDA0000950734160000115
为基质孔隙度，Tsc为标况温度，VL为Langmuir体积，pL为Langmuir压力，cpm为基质孔隙压缩系数，
Figure BDA0000950734160000116
表示压力为平均地层压力
Figure BDA0000950734160000117
时的气体压缩因子。
进一步的，平均地层压力获取模块305具体用于：
平均地层压力
Figure BDA0000950734160000118
根据如下公式确定：
Figure BDA0000950734160000119
其中，pi为原始地层压力，G为气井控制储量，zg(pi)表示压力为原始地层压力pi时的气体压缩因子。
进一步的，根据当前拟时间下的无因次产量qD，获得真实产量qsc具体用于：
Figure BDA00009507341600001110
其中，kF为裂缝渗透率，Acw为泄流面积，ψi为原始地层条件下的拟压力，ψwf为井底压力条件下的拟压力，T为储层温度。
本实施例是与方法实施例二对应的装置实施例，具体可参见实施例二中的描述，在此不再赘述。
虽然已经参考优选实施例对本发明进行了描述，但在不脱离本发明的范围的情况下，可以对其进行各种改进并且可以用等效物替换其中的部件。尤其是，只要不存在结构冲突，各个实施例中所提到的各项技术特征均可以任意方式组合起来。本发明并不局限于文中公开的特定实施例，而是包括落入权利要求的范围内的所有技术方案。

Claims (8)

1.一种页岩气井产量构成确定方法，其特征在于，包括：
步骤100，设置当前无因次拟时间下的平均地层压力
Figure FDA0002407612580000011
使得当前无因次拟时间下的平均地层压力
Figure FDA0002407612580000012
与原始地层压力相同；
步骤101，根据无因次产能公式，获取拉式空间下的无因次产量
Figure FDA0002407612580000013
步骤102，根据拉式空间下的无因次产量
Figure FDA0002407612580000014
利用数值反演法计算获得当前拟时间下的无因次产量qD；
步骤103，根据当前拟时间下的无因次产量qD，计算获得累积产量Gp；
步骤104，根据累积产量Gp，计算获得平均地层压力
Figure FDA0002407612580000015
步骤105，判断
Figure FDA0002407612580000016
是否成立，σ为预设阈值，若成立，则执行步骤106，若不成立，则使
Figure FDA0002407612580000017
转步骤103执行；
步骤106，根据平均地层压力
Figure FDA0002407612580000018
计算获得页岩气井的产量构成，
步骤103具体包括：
计算当前无因次拟时间下的平均地层压力
Figure FDA0002407612580000019
的气体物性参数，所述气体物性参数包括气体压缩因子zg、气体密度ρg、气体粘度μg、气体压缩系数cg；
根据气体物性参数，获取拟时间ta；
根据拟时间，获得真实时间t；
根据当前拟时间下的无因次产量qD，获得真实产量qsc；
根据真实产量qsc，计算获得真实时间t内的累积产量Gp。
2.根据权利要求1所述的页岩气井产量构成确定方法，其特征在于，步骤106具体包括：
Figure FDA00024076125800000110
其中，
Figure FDA00024076125800000111
为解吸压缩系数，该参数是平均地层压力
Figure FDA00024076125800000112
的函数，
Figure FDA00024076125800000113
psc为标况压力，T为储层温度，
Figure FDA00024076125800000114
为基质孔隙度，Tsc为标况温度，VL为Langmuir体积，pL为Langmuir压力，cpm为基质孔隙压缩系数，
Figure FDA0002407612580000021
表示压力为平均地层压力
Figure FDA0002407612580000022
时的气体压缩因子。
3.根据权利要求1所述的页岩气井产量构成确定方法，其特征在于，步骤104具体包括：
平均地层压力
Figure FDA0002407612580000023
根据如下公式确定：
Figure FDA0002407612580000024
其中，pi为原始地层压力，G为气井控制储量，zg(pi)表示压力为原始地层压力pi时的气体压缩因子。
4.根据权利要求1所述的页岩气井产量构成确定方法，其特征在于，根据当前拟时间下的无因次产量qD，获得真实产量qsc具体包括：
Figure FDA0002407612580000025
其中，kF为裂缝渗透率，Acw为泄流面积，ψi为原始地层条件下的拟压力，ψwf为井底压力条件下的拟压力，T为储层温度。
5.一种页岩气井产量构成确定装置，其特征在于，包括：
初始化模块，用于设置当前无因次拟时间下的平均地层压力
Figure FDA0002407612580000026
使得当前无因次拟时间下的平均地层压力
Figure FDA0002407612580000027
与原始地层压力相同；
拉式空间下的无因次产量获取模块，用于根据无因次产能公式，获取拉式空间下的无因次产量
Figure FDA0002407612580000028
当前拟时间下的无因次产量获取模块，用于根据拉式空间下的无因次产量
Figure FDA0002407612580000029
利用数值反演法计算获得当前拟时间下的无因次产量qD；
累积产量获取模块，用于根据当前拟时间下的无因次产量qD，计算获得累积产量Gp；
平均地层压力获取模块，用于根据累积产量Gp，计算获得平均地层压力
Figure FDA00024076125800000210
判断模块，用于判断
Figure FDA00024076125800000211
是否成立，σ为预设阈值，若成立，则触发页岩气井产量构成计算模块，若不成立，则使
Figure FDA00024076125800000212
触发累积产量获取模块；
页岩气井产量构成计算模块，用于根据平均地层压力
Figure FDA00024076125800000213
计算获得页岩气井的产量构成，
累积产量获取模块具体用于：
计算当前无因次拟时间下的平均地层压力
Figure FDA0002407612580000031
的气体物性参数，所述气体物性参数包括气体压缩因子zg、气体密度ρg、气体粘度μg、气体压缩系数cg；
根据气体物性参数，获取拟时间ta；
根据拟时间，获得真实时间t；
根据当前拟时间下的无因次产量qD，获得真实产量qsc；
根据真实产量qsc，计算获得真实时间t内的累积产量Gp。
6.根据权利要求5所述的页岩气井产量构成确定装置，其特征在于，页岩气井产量构成计算模块具体用于：
Figure FDA0002407612580000032
其中，
Figure FDA0002407612580000033
为解吸压缩系数，该参数是平均地层压力
Figure FDA0002407612580000034
的函数，
Figure FDA0002407612580000035
psc为标况压力，T为储层温度，
Figure FDA0002407612580000036
为基质孔隙度，Tsc为标况温度，VL为Langmuir体积，pL为Langmuir压力，cpm为基质孔隙压缩系数，
Figure FDA0002407612580000037
表示压力为平均地层压力
Figure FDA0002407612580000038
时的气体压缩因子。
7.根据权利要求5所述的页岩气井产量构成确定装置，其特征在于，平均地层压力获取模块具体用于：
平均地层压力
Figure FDA0002407612580000039
根据如下公式确定：
Figure FDA00024076125800000310
其中，pi为原始地层压力，G为气井控制储量，zg(pi)表示压力为原始地层压力pi时的气体压缩因子。
8.根据权利要求5所述的页岩气井产量构成确定装置，其特征在于，根据当前拟时间下的无因次产量qD，获得真实产量qsc具体用于：
Figure FDA0002407612580000041
其中，kF为裂缝渗透率，Acw为泄流面积，ψi为原始地层条件下的拟压力，ψwf为井底压力条件下的拟压力，T为储层温度。
CN201610180381.5A 2016-03-25 2016-03-25 页岩气井产量构成确定方法及装置 Active CN107239648B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610180381.5A CN107239648B (zh) 2016-03-25 2016-03-25 页岩气井产量构成确定方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610180381.5A CN107239648B (zh) 2016-03-25 2016-03-25 页岩气井产量构成确定方法及装置

Publications (2)

Publication Number Publication Date
CN107239648A CN107239648A (zh) 2017-10-10
CN107239648B true CN107239648B (zh) 2020-07-10

Family

ID=59982593

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610180381.5A Active CN107239648B (zh) 2016-03-25 2016-03-25 页岩气井产量构成确定方法及装置

Country Status (1)

Country Link
CN (1) CN107239648B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109958429B (zh) * 2017-12-21 2021-10-22 中国石油化工股份有限公司 一种判断页岩气吸附气产出的方法
CN108798654B (zh) * 2018-04-27 2021-07-02 中国石油天然气股份有限公司 确定页岩气井的井底压力与时间的对应关系的方法及装置
CN111506865B (zh) * 2020-07-01 2020-09-15 西南石油大学 一种页岩气井生产拟合与预测方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104389594A (zh) * 2014-10-13 2015-03-04 成都创源油气技术开发有限公司 页岩气井产能评价预测方法
CN104632187A (zh) * 2013-11-14 2015-05-20 中国石油化工股份有限公司 一种确定产水煤层气井动态储量的方法
CN104948163A (zh) * 2014-03-24 2015-09-30 中国石油化工股份有限公司 一种页岩气井产能测定方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8788252B2 (en) * 2010-10-26 2014-07-22 Schlumberger Technology Corporation Multi-well time-lapse nodal analysis of transient production systems
US20140083687A1 (en) * 2012-09-27 2014-03-27 Schlumberger Technology Corporation Production in fractured systems

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104632187A (zh) * 2013-11-14 2015-05-20 中国石油化工股份有限公司 一种确定产水煤层气井动态储量的方法
CN104948163A (zh) * 2014-03-24 2015-09-30 中国石油化工股份有限公司 一种页岩气井产能测定方法
CN104389594A (zh) * 2014-10-13 2015-03-04 成都创源油气技术开发有限公司 页岩气井产能评价预测方法

Non-Patent Citations (7)

* Cited by examiner, † Cited by third party
Title
"Accounting For Adsorbed Gas in Shale Gas Reservoirs";Salman A. Mengal et al;《SPE Middle East Oil and Gas Show and Conference》;20111231;全文 *
"Production Data Analysis of Shale Gas Reservoirs";Adam Michael Lewis et al;《SPE Annual Technical Conference and Exhibition》;20081231;全文 *
"Type Curves for Dry CBM Reservoirs with Equilibrium Desorption";S. GERAMI et al;《Journal of Canadian Petroleum Technology》;20071231;第47卷(第7期);全文 *
"一种裂缝性油藏产能预测新方法";郑学锐等;《断块油气田》;20151130;第22卷(第6期);第744-746、751页 *
"基于试井技术的页岩气藏非稳态产量分析";程子洋等;《第三届全国特殊气藏开发技术研讨会优秀论文集》;20141031;全文 *
"考虑吸附作用的页岩气井生产数据分析方法";丁志文等;《石油化工应用》;20131231;第32卷(第12期);全文 *
"页岩气体积压裂水平井非稳态产能评价模型";艾爽等;《天然气地球科学》;20141031;第25卷(第10期);全文 *

Also Published As

Publication number Publication date
CN107239648A (zh) 2017-10-10

Similar Documents

Publication Publication Date Title
CN106351651B (zh) 页岩气井产能的预测方法及装置
CN104196524B (zh) 一种欠饱和煤储层开发的气水产出动态相渗曲线测量方法
Wang et al. Improved apparent permeability models of gas flow in coal with Klinkenberg effect
CN111911115B (zh) 一种页岩气井动态配产方法
CN104612635A (zh) 煤层群瓦斯联合抽采达标预评判方法
CN107239648B (zh) 页岩气井产量构成确定方法及装置
CN108166963B (zh) 一种海上油气井压裂效果评价方法
CN110702585B (zh) 一种岩石压缩系数计算方法
CN107435528B (zh) 火山岩气藏气井配产的方法
CN108133086A (zh) 一种应力敏感储层中产水压裂气井裂缝半长反演方法
CN108952657B (zh) 一种水平井平台压裂裂缝长度确定方法及装置
Peng et al. Benchmark assessment of coal permeability models on the accuracy of permeability prediction
CN110348590A (zh) 一种页岩吸附气含量计算方法
RU2016102698A (ru) Способ и система для калибровки статической модели геологической среды с использованием исследования проницаемости
Zhang et al. A fluid‐solid coupling method for the simulation of gas transport in porous coal and rock media
CN110032744B (zh) 一种非常规油气储量计算方法及系统
Ebeltoft et al. Parameter based scal-analysing relative permeability for full field application
CN111626530B (zh) 确定压裂页岩气藏井控储量的方法及装置
Mengal Accounting for Adsorbed gas and its effect on production bahavior of Shale Gas Reservoirs
CN108153944A (zh) 一种应力敏感储层产水气井天然气单井控制储量确定方法
CN108952697B (zh) 压回法试验模拟系统及模拟方法
CN109426657B (zh) 一种油藏井间连通系数的计算方法及系统
CN103590813B (zh) 两层合采气井动态储量的测试方法、测试管柱及测试装置
Fattah et al. New inflow performance relationship for solution-gas drive oil reservoirs
US8146659B2 (en) Method of modelling enhanced recovery by polymer injection

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant