CN107239648A - 页岩气井产量构成确定方法及装置 - Google Patents
页岩气井产量构成确定方法及装置 Download PDFInfo
- Publication number
- CN107239648A CN107239648A CN201610180381.5A CN201610180381A CN107239648A CN 107239648 A CN107239648 A CN 107239648A CN 201610180381 A CN201610180381 A CN 201610180381A CN 107239648 A CN107239648 A CN 107239648A
- Authority
- CN
- China
- Prior art keywords
- msub
- mrow
- production
- under
- pseudotime
- 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.)
- Granted
Links
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
- Production Of Liquid Hydrocarbon Mixture For Refining Petroleum (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供一种页岩气井产量构成确定方法及装置，其中方法包括：步骤100，设置当前无因次拟时间下的平均地层压力使得与原始地层压力相同；步骤101，根据无因次产能公式，获取拉式空间下的无因次产量步骤102，根据利用数值反演法计算获得当前拟时间下的无因次产量qD；步骤103，根据qD，计算获得累积产量Gp；步骤104，根据Gp，计算获得平均地层压力步骤105，判断是否成立，σ为预设阈值，若成立，则执行步骤106，若不成立，则使转步骤103执行；步骤106，根据计算获得页岩气井的产量构成。通过本发明中的方法计算出的页岩气井构成情况并非固定值，而是随时间发生变化，计算出的页岩气井产量构成与实际情况更相符。
Description
技术领域
本发明涉及气藏开发技术领域，尤其涉及一种页岩气井产量构成确定方法及装置。
背景技术
页岩气井中生产的天然气有两种来源：吸附气和游离气。解吸压缩系数可用于表征吸附气的供给作用，解吸压缩系数占总压缩系数的比重可看做吸附气天然气总气量的比重。现在普遍认为解吸压缩系数是一个常数，由此计算出的页岩气井的产量构成，吸附气和游离气占总产气量的比重都是固定的。而实际上由于页岩气藏普遍衰竭式开发，在开发过程中地层压力不断下降，此时解吸压缩系数也随地层压力的变化而变化，因此忽略解吸压缩系数随地层压力的变化所计算出来的页岩气井产量构成是不合理的。因此，亟需一种考虑解吸压缩系数随地层压力变化而变化的页岩气井产量构成确定方法。
发明内容
本发明提供一种页岩气井产量构成确定方法及装置，用以解决现有技术中计算出来的页岩气井产量构成不合理的技术问题。
本发明一方面提供一种页岩气井产量构成确定方法，包括：
步骤100，设置当前无因次拟时间下的平均地层压力使得当前无因次拟时间下的平均地层压力与原始地层压力相同；
步骤101，根据无因次产能公式，获取拉式空间下的无因次产量
步骤102，根据拉式空间下的无因次产量利用数值反演法计算获得当前拟时间下的无因次产量qD；
步骤103，根据当前拟时间下的无因次产量qD，计算获得累积产量Gp；
步骤104，根据累积产量Gp，计算获得平均地层压力
步骤105，判断是否成立，σ为预设阈值，若成立，则执行步骤106，若不成立，则使转步骤103执行；
步骤106，根据平均地层压力计算获得页岩气井的产量构成。
进一步的，步骤103具体包括：
计算当前无因次拟时间下的平均地层压力的气体物性参数，所述气体物性参数包括气体压缩因子zg、气体密度ρg、气体粘度μg、气体压缩系数cg；
根据气体物性参数，获取拟时间ta；
根据拟时间，获得真实时间t；
根据当前拟时间下的无因次产量qD，获得真实产量qsc；
根据真实产量qsc，计算获得真实时间t内的累积产量Gp。
进一步的，步骤106具体包括：
其中，为解吸压缩系数，该参数是平均地层压力的函数，psc为标况压力，T为储层温度，为基质孔隙度，Tsc为标况温度，VL为Langmuir体积，pL为Langmuir压力，cpm为基质孔隙压缩系数，表示压力为平均地层压力时的气体压缩因子。
进一步的，步骤104具体包括：
平均地层压力根据如下公式确定：
其中，pi为原始地层压力，G为气井控制储量，zg(pi)表示压力为原始地层压力pi时的气体压缩因子。
进一步的，根据当前拟时间下的无因次产量qD，获得真实产量qsc具体包括：
其中，kF为裂缝渗透率，Acw为泄流面积，ψi为原始地层条件下的拟压力，ψwf为井底压力条件下的拟压力，T为储层温度。
本发明另一方面提供一种页岩气井产量构成确定装置，包括：
初始化模块，用于设置当前无因次拟时间下的平均地层压力使得当前无因次拟时间下的平均地层压力与原始地层压力相同；
拉式空间下的无因次产量获取模块，用于根据无因次产能公式，获取拉式空间下的无因次产量
当前拟时间下的无因次产量获取模块，用于根据拉式空间下的无因次产量利用数值反演法计算获得当前拟时间下的无因次产量qD；
累积产量获取模块，用于根据当前拟时间下的无因次产量qD，计算获得累积产量Gp；
平均地层压力获取模块，用于根据累积产量Gp，计算获得平均地层压力
判断模块，用于判断是否成立，σ为预设阈值，若成立，则触发页岩气井产量构成计算模块，若不成立，则使触发累积产量获取模块；
页岩气井产量构成计算模块，用于根据平均地层压力计算获得页岩气井的产量构成。
进一步的，累积产量获取模块具体用于：
计算当前无因次拟时间下的平均地层压力的气体物性参数，所述气体物性参数包括气体压缩因子zg、气体密度ρg、气体粘度μg、气体压缩系数cg；
根据气体物性参数，获取拟时间ta；
根据拟时间，获得真实时间t；
根据当前拟时间下的无因次产量qD，获得真实产量qsc；
根据真实产量qsc，计算获得真实时间t内的累积产量Gp。
进一步的，页岩气井产量构成计算模块具体用于：
其中，为解吸压缩系数，该参数是平均地层压力的函数，psc为标况压力，T为储层温度，为基质孔隙度，Tsc为标况温度，VL为Langmuir体积，pL为Langmuir压力，cpm为基质孔隙压缩系数，表示压力为平均地层压力时的气体压缩因子。
进一步的，平均地层压力获取模块具体用于：
平均地层压力根据如下公式确定：
其中，pi为原始地层压力，G为气井控制储量，zg(pi)表示压力为原始地层压力pi时的气体压缩因子。
进一步的，根据当前拟时间下的无因次产量qD，获得真实产量qsc具体用于：
其中，kF为裂缝渗透率，Acw为泄流面积，ψi为原始地层条件下的拟压力，ψwf为井底压力条件下的拟压力，T为储层温度。
本发明提供的页岩气井产量构成确定方法及装置，通过根据当前拟时间下的无因次产量，计算获得累积产量；然后根据累积产量，计算获得平均地层压力；最后通过计算获得平均地层压力来确定页岩气井的产量构成，即页岩气井的吸附气构成和游离气构成，从而使计算出的页岩气井构成情况并非固定值，而是随时间发生变化，通过本发明中的方法计算出的页岩气井产量与实际情况更相符。
附图说明
在下文中将基于实施例并参考附图来对本发明进行更详细的描述。其中：
图1为根据本发明实施例一的页岩气井产量构成确定方法的流程示意图；
图2为根据本发明实施例二的页岩气井产量构成确定方法的流程示意图；
图3为根据本发明实施例三的页岩气井产量构成确定装置的流程示意图。
在附图中，相同的部件使用相同的附图标记。附图并未按照实际的比例绘制。
具体实施方式
下面将结合附图对本发明作进一步说明。
实施例一
图1为根据本发明实施例一的页岩气井产量构成确定方法的流程示意图；如图1所示，本实施例提供一种页岩气井产量构成确定方法，包括：
步骤100，设置当前无因次拟时间下的平均地层压力使得当前无因次拟时间下的平均地层压力与原始地层压力相同。
具体的，当前无因次拟时间可在0.001～1000000中任意取值，可根据实际需要选取；平均地层压力指生产进行到某一时间时地层的平均压力，原始地层压力是指开采前的地层压力，原始地层压力的值可通过查阅地质资料或者试井解释报告获得。
步骤101，根据无因次产能公式，获取拉式空间下的无因次产量
具体的，拉式空间是一种虚拟的空间，先把无因次产能公式转换到拉式空间，然后后续步骤中利用数值反演法反演到实际空间，就可以求出无因次产量。
无因次产量可根据如下公式获取：
其中，yeD为无因次缝长，ω为弹性储容比，λcw为窜流系数，s为拉普拉斯算子。
步骤102，根据拉式空间下的无因次产量利用数值反演法计算获得当前拟时间下的无因次产量qD。
具体的，数值反演法包括stefest算法，可利用stefest算法计算获得当前拟时间下的无因次产量qD，“当前拟时间”是针对“当前无因次拟时间”而言的，把当前无因次拟时间变换为有因次，则为当前拟时间。
步骤103，根据当前拟时间下的无因次产量qD，计算获得累积产量Gp。
步骤104，根据累积产量Gp，计算获得平均地层压力
步骤105，判断是否成立，σ为预设阈值，若成立，则执行步骤106，若不成立，则使转步骤103执行。
具体的，σ的取值范围为0到0.005之间，优选取0.001。当前拟时间下的平均地层压力与之间差值小于σ时，转步骤106继续执行，否则，需要再次转步骤103，重新计算的值。
步骤106，根据平均地层压力计算获得页岩气井的产量构成。
本实施例提供的页岩气井产量构成确定方法，通过根据当前拟时间下的无因次产量，计算获得累积产量；然后根据累积产量，计算获得平均地层压力；最后通过计算获得平均地层压力来确定页岩气井的产量构成，即页岩气井的吸附气构成和游离气构成，从而使计算出的页岩气井构成情况并非固定值，而是随时间发生变化。
实施例二
本实施例是在上述实施例的基础上进行的补充说明。
图2为根据本发明实施例二的页岩气井产量构成确定方法的流程示意图；如图2所示，本实施例提供一种页岩气井产量构成确定方法，包括：
步骤200，设置当前无因次拟时间下的平均地层压力使得当前无因次拟时间下的平均地层压力与原始地层压力相同。
步骤201，根据无因次产能公式，获取拉式空间下的无因次产量
步骤202，根据拉式空间下的无因次产量利用数值反演法计算获得当前拟时间下的无因次产量qD。
步骤200-步骤202与实施例一中的步骤100-步骤102相同，具体可参见实施例一中相应的描述，在此不再赘述。
步骤2031，计算当前无因次拟时间下的平均地层压力的气体物性参数，所述气体物性参数包括气体压缩因子zg、气体密度ρg、气体粘度μg、气体压缩系数cg。
具体的，气体压缩因子zg可通过如下公式进行计算获得：
其中，Tpr为对比温度。
气体密度ρg的确定方法可根据如下公式获得：
其中，γg表示气体相对密度，T为储层温度，当的值还未确定时，
气体粘度μg的确定方法可根据如下公式获得：
其中，Mg为气体分子量，T为储层温度。
气体压缩系数cg的确定方法可根据如下公式获得：
步骤2032，根据气体物性参数，获取拟时间ta。
具体的，拟时间ta可根据如下公式确定：
其中，pi为原始地层压力，cpm为基质孔隙压缩系数；Acw为泄流面积；tDa为当前无因次拟时间；cpF为裂缝压缩系数，kF为裂缝渗透率，φm为基质孔隙度，φF为裂缝孔隙度。
上式中表示对zg(pi)中的pi求导数，T代表储层温度，表明该气体压缩系数是在等温条件下获得的。
步骤2033，根据拟时间ta，获得真实时间t。
真实时间的确定方法：
其中，
步骤2034，根据当前拟时间下的无因次产量qD，获得真实产量qsc。
其中，kF为裂缝渗透率，Acw为泄流面积，ψi为原始地层条件下的拟压力，ψwf为井底压力条件下的拟压力，T为储层温度。
步骤2035，根据真实产量qsc，计算获得真实时间t内的累积产量Gp。
具体的，累积产量Gp的确定方法为：
步骤204，根据累积产量Gp，计算获得平均地层压力
平均地层压力可根据如下公式确定：
其中，pi为原始地层压力，G为气井控制储量，zg是压力的函数，zg(pi)表示压力为原始地层压力pi时的气体压缩因子zg。
步骤205，判断是否成立，σ为预设阈值，若成立，则执行步骤206，若不成立，则使转步骤2031执行。
步骤206，根据平均地层压力计算获得页岩气藏产井的产量构成。
其中，为解吸压缩系数，该参数是平均地层压力的函数，psc为标况压力，T为储层温度，为基质孔隙度，Tsc为标况温度，VL为Langmuir体积，pL为Langmuir压力，cpm为基质孔隙压缩系数，表示压力为平均地层压力时的气体压缩因子。
本实施例提供的页岩气井产量构成确定方法，通过根据当前拟时间下的无因次产量，计算获得累积产量；然后根据累积产量，计算获得平均地层压力；最后通过计算获得平均地层压力来确定页岩气井的产量构成，即页岩气井的吸附气构成和游离气构成，从而使计算出的页岩气井构成情况并非固定值，而是随时间发生变化。本实施例中的求取吸附气构成和游离气构成的方法，考虑了解吸压缩系数随地层压力变化而产生的变化，从而使计算出的页岩气井构成情况随时间发生变化。
实施例三
本实施例为装置实施例，用于执行上述实施例一中的方法。
图3为根据本发明实施例三的页岩气井产量构成确定装置的结构示意图；如图3所示，本实施例提供一种页岩气井产量构成确定装置，包括初始化模块301、拉式空间下的无因次产量获取模块302、当前拟时间下的无因次产量获取模块303、累积产量获取模块304、平均地层压力获取模块305、判断模块306和页岩气井产量构成计算模块307。
其中，初始化模块301，用于设置当前无因次拟时间下的平均地层压力使得当前无因次拟时间下的平均地层压力与原始地层压力相同；
拉式空间下的无因次产量获取模块302，用于根据无因次产能公式，获取拉式空间下的无因次产量
当前拟时间下的无因次产量获取模块303，用于根据拉式空间下的无因次产量利用数值反演法计算获得当前拟时间下的无因次产量qD；
累积产量获取模块304，用于根据当前拟时间下的无因次产量qD，计算获得累积产量Gp；
平均地层压力获取模块305，用于根据累积产量Gp，计算获得平均地层压力
判断模块306，用于判断是否成立，σ为预设阈值，若成立，则触发页岩气井产量构成计算模块307，若不成立，则使触发累积产量获取模块304；
页岩气井产量构成计算模块307，用于根据平均地层压力计算获得页岩气井的产量构成。
本实施例是与方法实施例一对应的装置实施例，具体可参见实施例一中的描述，在此不再赘述。
实施例四
本实施例是在实施例三的基础上进行的补充说明，用于执行上述实施例二中的方法。
本实施例提供一种页岩气井产量构成确定装置，包括初始化模块301、拉式空间下的无因次产量获取模块302、当前拟时间下的无因次产量获取模块303、累积产量获取模块304、平均地层压力获取模块305、判断模块306和页岩气井产量构成计算模块307。
其中，累积产量获取模块304具体用于：
计算当前无因次拟时间下的平均地层压力的气体物性参数，所述气体物性参数包括气体压缩因子zg、气体密度ρg、气体粘度μg、气体压缩系数cg；
根据气体物性参数，获取拟时间ta；
根据拟时间，获得真实时间t；
根据当前拟时间下的无因次产量qD，获得真实产量qsc；
根据真实产量qsc，计算获得真实时间t内的累积产量Gp。
页岩气井产量构成计算模块307具体用于：
其中，为解吸压缩系数，该参数是平均地层压力的函数，psc为标况压力，T为储层温度，为基质孔隙度，Tsc为标况温度，VL为Langmuir体积，pL为Langmuir压力，cpm为基质孔隙压缩系数，表示压力为平均地层压力时的气体压缩因子。
进一步的，平均地层压力获取模块305具体用于：
平均地层压力根据如下公式确定：
其中，pi为原始地层压力，G为气井控制储量，zg(pi)表示压力为原始地层压力pi时的气体压缩因子。
进一步的，根据当前拟时间下的无因次产量qD，获得真实产量qsc具体用于：
其中，kF为裂缝渗透率，Acw为泄流面积，ψi为原始地层条件下的拟压力，ψwf为井底压力条件下的拟压力，T为储层温度。
本实施例是与方法实施例二对应的装置实施例，具体可参见实施例二中的描述，在此不再赘述。
虽然已经参考优选实施例对本发明进行了描述，但在不脱离本发明的范围的情况下，可以对其进行各种改进并且可以用等效物替换其中的部件。尤其是，只要不存在结构冲突，各个实施例中所提到的各项技术特征均可以任意方式组合起来。本发明并不局限于文中公开的特定实施例，而是包括落入权利要求的范围内的所有技术方案。
Claims (10)
1.一种页岩气井产量构成确定方法，其特征在于，包括：
步骤100，设置当前无因次拟时间下的平均地层压力使得当前无因次拟时间下的平均地层压力与原始地层压力相同；
步骤101，根据无因次产能公式，获取拉式空间下的无因次产量
步骤102，根据拉式空间下的无因次产量利用数值反演法计算获得当前拟时间下的无因次产量qD；
步骤103，根据当前拟时间下的无因次产量qD，计算获得累积产量Gp；
步骤104，根据累积产量Gp，计算获得平均地层压力
步骤105，判断是否成立，σ为预设阈值，若成立，则执行步骤106，若不成立，则使转步骤103执行；
步骤106，根据平均地层压力计算获得页岩气井的产量构成。
2.根据权利要求1所述的页岩气井产量构成确定方法，其特征在于，步骤103具体包括：
计算当前无因次拟时间下的平均地层压力的气体物性参数，所述气体物性参数包括气体压缩因子zg、气体密度ρg、气体粘度μg、气体压缩系数cg；
根据气体物性参数，获取拟时间ta；
根据拟时间，获得真实时间t；
根据当前拟时间下的无因次产量qD，获得真实产量qsc；
根据真实产量qsc，计算获得真实时间t内的累积产量Gp。
3.根据权利要求1所述的页岩气井产量构成确定方法，其特征在于，步骤106具体包括：
其中，为解吸压缩系数，该参数是平均地层压力的函数， psc为标况压力，T为储层温度，为基质孔隙度，Tsc为标况温度，VL为Langmuir体积，pL为Langmuir压力，cpm为基质孔隙压缩系数，表示压力为平均地层压力时的气体压缩因子。
4.根据权利要求1所述的页岩气井产量构成确定方法，其特征在于，步骤104具体包括：
平均地层压力根据如下公式确定：
<mrow>
<mfrac>
<msub>
<mover>
<mi>p</mi>
<mo>&OverBar;</mo>
</mover>
<mi>m</mi>
</msub>
<mrow>
<msub>
<mi>z</mi>
<mi>g</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>p</mi>
<mo>&OverBar;</mo>
</mover>
<mi>m</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<msub>
<mi>p</mi>
<mi>i</mi>
</msub>
<mrow>
<msub>
<mi>z</mi>
<mi>g</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>p</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<msub>
<mi>G</mi>
<mi>P</mi>
</msub>
<mi>G</mi>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
其中，pi为原始地层压力，G为气井控制储量，zg(pi)表示压力为原始地层压力pi时的气体压缩因子。
5.根据权利要求2所述的页岩气井产量构成确定方法，其特征在于，根据当前拟时间下的无因次产量qD，获得真实产量qsc具体包括：
<mrow>
<msub>
<mi>q</mi>
<mi>sc</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>k</mi>
<mi>F</mi>
</msub>
<msqrt>
<msub>
<mi>A</mi>
<mi>cw</mi>
</msub>
</msqrt>
<mrow>
<mo>(</mo>
<msub>
<mi>&psi;</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>&psi;</mi>
<mi>wf</mi>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>q</mi>
<mi>D</mi>
</msub>
</mrow>
<mrow>
<mn>1.291</mn>
<mo>&times;</mo>
<msup>
<mn>10</mn>
<mrow>
<mo>-</mo>
<mn>3</mn>
</mrow>
</msup>
<mi>T</mi>
</mrow>
</mfrac>
</mrow>
其中，kF为裂缝渗透率，Acw为泄流面积，ψi为原始地层条件下的拟压力，ψwf为井底压力条件下的拟压力，T为储层温度。
6.一种页岩气井产量构成确定装置，其特征在于，包括：
初始化模块，用于设置当前无因次拟时间下的平均地层压力使得当前无因次拟时间下的平均地层压力与原始地层压力相同；
拉式空间下的无因次产量获取模块，用于根据无因次产能公式，获取拉式空间下的无因次产量
当前拟时间下的无因次产量获取模块，用于根据拉式空间下的无因次产量利用数值反演法计算获得当前拟时间下的无因次产量qD；
累积产量获取模块，用于根据当前拟时间下的无因次产量qD，计算获得累积产量Gp；
平均地层压力获取模块，用于根据累积产量Gp，计算获得平均地层压力
判断模块，用于判断是否成立，σ为预设阈值，若成立，则触发页岩气井产量构成计算模块，若不成立，则使触发累积产量获取模块；
页岩气井产量构成计算模块，用于根据平均地层压力计算获得页岩气井的产量构成。
7.根据权利要求6所述的页岩气井产量构成确定装置，其特征在于，累积产量获取模块具体用于：
计算当前无因次拟时间下的平均地层压力的气体物性参数，所述气体物性参数包括气体压缩因子zg、气体密度ρg、气体粘度μg、气体压缩系数cg；
根据气体物性参数，获取拟时间ta；
根据拟时间，获得真实时间t；
根据当前拟时间下的无因次产量qD，获得真实产量qsc；
根据真实产量qsc，计算获得真实时间t内的累积产量Gp。
8.根据权利要求6所述的页岩气井产量构成确定装置，其特征在于，页岩气井产量构成计算模块具体用于：
其中，为解吸压缩系数，该参数是平均地层压力的函数， psc为标况压力，T为储层温度，为基质孔隙度，Tsc为标况温度，VL为Langmuir体积，pL为Langmuir压力，cpm为基质孔隙压缩系数，表示压力为平均地层压力时的气体压缩因子。
9.根据权利要求6所述的页岩气井产量构成确定装置，其特征在于，平均地层压力获取模块具体用于：
平均地层压力根据如下公式确定：
<mrow>
<mfrac>
<msub>
<mover>
<mi>p</mi>
<mo>&OverBar;</mo>
</mover>
<mi>m</mi>
</msub>
<mrow>
<msub>
<mi>z</mi>
<mi>g</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>p</mi>
<mo>&OverBar;</mo>
</mover>
<mi>m</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<msub>
<mi>p</mi>
<mi>i</mi>
</msub>
<mrow>
<msub>
<mi>z</mi>
<mi>g</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>p</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<msub>
<mi>G</mi>
<mi>P</mi>
</msub>
<mi>G</mi>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
其中，pi为原始地层压力，G为气井控制储量，zg(pi)表示压力为原始地层压力pi时的气体压缩因子。
10.根据权利要求7所述的页岩气井产量构成确定装置，其特征在于，根据当前拟时间下的无因次产量qD，获得真实产量qsc具体用于：
<mrow>
<msub>
<mi>q</mi>
<mi>sc</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>k</mi>
<mi>F</mi>
</msub>
<msqrt>
<msub>
<mi>A</mi>
<mi>cw</mi>
</msub>
</msqrt>
<mrow>
<mo>(</mo>
<msub>
<mi>&psi;</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>&psi;</mi>
<mi>wf</mi>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>q</mi>
<mi>D</mi>
</msub>
</mrow>
<mrow>
<mn>1.291</mn>
<mo>&times;</mo>
<msup>
<mn>10</mn>
<mrow>
<mo>-</mo>
<mn>3</mn>
</mrow>
</msup>
<mi>T</mi>
</mrow>
</mfrac>
</mrow>
其中，kF为裂缝渗透率，Acw为泄流面积，ψi为原始地层条件下的拟压力，ψwf为井底压力条件下的拟压力，T为储层温度。
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 true CN107239648A (zh) | 2017-10-10 |
| CN107239648B 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) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108798654A (zh) * | 2018-04-27 | 2018-11-13 | 中国石油天然气股份有限公司 | 确定页岩气井的井底压力与时间的对应关系的方法及装置 |
| CN109958429A (zh) * | 2017-12-21 | 2019-07-02 | 中国石油化工股份有限公司 | 一种判断页岩气吸附气产出的方法 |
| CN111506865A (zh) * | 2020-07-01 | 2020-08-07 | 西南石油大学 | 一种页岩气井生产拟合与预测方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20120101787A1 (en) * | 2010-10-26 | 2012-04-26 | Wentao Zhou | 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 |
| CN104389594A (zh) * | 2014-10-13 | 2015-03-04 | 成都创源油气技术开发有限公司 | 页岩气井产能评价预测方法 |
| CN104632187A (zh) * | 2013-11-14 | 2015-05-20 | 中国石油化工股份有限公司 | 一种确定产水煤层气井动态储量的方法 |
| CN104948163A (zh) * | 2014-03-24 | 2015-09-30 | 中国石油化工股份有限公司 | 一种页岩气井产能测定方法 |
-
2016
- 2016-03-25 CN CN201610180381.5A patent/CN107239648B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20120101787A1 (en) * | 2010-10-26 | 2012-04-26 | Wentao Zhou | 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 |
| 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)
| Title |
|---|
| ADAM MICHAEL LEWIS ET AL: ""Production Data Analysis of Shale Gas Reservoirs"", 《SPE ANNUAL TECHNICAL CONFERENCE AND EXHIBITION》 * |
| S. GERAMI ET AL: ""Type Curves for Dry CBM Reservoirs with Equilibrium Desorption"", 《JOURNAL OF CANADIAN PETROLEUM TECHNOLOGY》 * |
| SALMAN A. MENGAL ET AL: ""Accounting For Adsorbed Gas in Shale Gas Reservoirs"", 《SPE MIDDLE EAST OIL AND GAS SHOW AND CONFERENCE》 * |
| 丁志文等: ""考虑吸附作用的页岩气井生产数据分析方法"", 《石油化工应用》 * |
| 程子洋等: ""基于试井技术的页岩气藏非稳态产量分析"", 《第三届全国特殊气藏开发技术研讨会优秀论文集》 * |
| 艾爽等: ""页岩气体积压裂水平井非稳态产能评价模型"", 《天然气地球科学》 * |
| 郑学锐等: ""一种裂缝性油藏产能预测新方法"", 《断块油气田》 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109958429A (zh) * | 2017-12-21 | 2019-07-02 | 中国石油化工股份有限公司 | 一种判断页岩气吸附气产出的方法 |
| CN109958429B (zh) * | 2017-12-21 | 2021-10-22 | 中国石油化工股份有限公司 | 一种判断页岩气吸附气产出的方法 |
| CN108798654A (zh) * | 2018-04-27 | 2018-11-13 | 中国石油天然气股份有限公司 | 确定页岩气井的井底压力与时间的对应关系的方法及装置 |
| CN108798654B (zh) * | 2018-04-27 | 2021-07-02 | 中国石油天然气股份有限公司 | 确定页岩气井的井底压力与时间的对应关系的方法及装置 |
| CN111506865A (zh) * | 2020-07-01 | 2020-08-07 | 西南石油大学 | 一种页岩气井生产拟合与预测方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107239648B (zh) | 2020-07-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104389594B (zh) | 页岩气井产能评价预测方法 | |
| CN101560879B (zh) | 用于低渗透气藏的试井分析控制系统及方法 | |
| CN107239648A (zh) | 页岩气井产量构成确定方法及装置 | |
| CN102980828B (zh) | 一种测量单管岩心泡沫驱气相饱和度的装置的工作方法 | |
| US20220049604A1 (en) | Shale gas well dynamic production allocating method | |
| WO2012115997A3 (en) | Apparatus and method for increasing the ultimate recovery of natural gas contained in shale and other tight gas reservoirs | |
| RU2013132014A (ru) | Система и способ для моделирования потока текучей среды в имеющем трещины пласте | |
| CN105893679B (zh) | 低产水平井续流修正试井解释方法 | |
| Manik et al. | Development and validation of a compositional coalbed simulator | |
| CN101806224A (zh) | 井下煤层气抽采产能预测方法 | |
| CN104612635A (zh) | 煤层群瓦斯联合抽采达标预评判方法 | |
| CN107387068B (zh) | 一种用于确定页岩气储层游离气含量的方法及系统 | |
| CN104949899B (zh) | 一种驱油用聚合物在多孔介质中有效粘度的测定方法 | |
| CN107387051A (zh) | 低渗透非均质油藏多段压裂水平井重复压裂选井的方法 | |
| CN104343445B (zh) | 识别气层与水层、高气油比储层及含气饱和度的方法 | |
| CN107451311A (zh) | 一种油水相对渗透率曲线计算方法及装置 | |
| CN106529184A (zh) | 倾斜有水气藏产水气井产能的计算方法 | |
| CN106156439A (zh) | 一种具有泡沫油现象的稠油油藏溶解气驱数值模拟方法 | |
| CN106202748A (zh) | 一种基于渗透率和井距的煤层气采收率预测方法 | |
| Zhang et al. | Numerical simulation of a coupled gas flow and geomechanics process in fractured coalbed methane reservoirs | |
| CN106127604A (zh) | 一种动态储量计算方法及装置 | |
| CN105930667B (zh) | 一种计算气井煤层气吸附量的方法 | |
| CN107725044A (zh) | 基于阵列感应、侧向测井的砂岩含气储层产水率预测的方法 | |
| Morad | Selected topics in coalbed methane reservoirs | |
| Yu | Developments in modeling and optimization of production in unconventional oil and gas reservoirs |
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 |