CN103046924B - 基于指数式的凝析气藏绝对无阻流量的获取方法及*** - Google Patents

Abstract

本发明涉及不同地层压力下凝析气藏绝对无阻流量的获取方法及***，该方法基于不同地层压力下凝析气藏内气体组份的变化和气相渗透率的变化对绝对无阻流量的获取；包括：确定不同地层压力下指数式产能方程；根据所述指数式产能方程，由第一地层压力和第二地层压力及第一地层压力对应的第一天然气粘度、第一偏差系数和第一气相渗透率，以及第二地层压力对应的第二天然气粘度、第二偏差系数和第二气相渗透率获取指数式产能方程的第一系数和第二系数之间的关系；根据指数式产能方程的第一系数和第二系数之间的关系、已知的第一凝析气藏绝对无阻流量以及第一地层压力和第二地层压力来获取不同地层压力下凝析气藏绝对无阻流量。

Description

基于指数式的凝析气藏绝对无阻流量的获取方法及***

技术领域

本发明涉及绝对无阻流量获取领域，特别涉及基于指数式产能方程的不同地层压力下凝析气藏绝对无阻流量的获取方法及***。

背景技术

文献1：陈元千，油气藏工程实践[M].北京:石油工业出版社，2005，142-146提出的方法中考虑了凝析油对气井产能的影响，通过公式表示为：

$P_e^2-P_{\mathrm{wf}}^2={\mathrm{aq}}_{\mathrm{sc}}+{\mathrm{bq}}_{\mathrm{sc}}^2$

其中： $a=\frac{\mathrm{TZ}{\mathrm{\μp}}_{\mathrm{sc}}}{\mathrm{\πkh}{T}_{\mathrm{sc}}{Z}_{\mathrm{sc}}}(\ln \frac{r_e}{r_w}+S_c+S_b);$ $b=\frac{\mathrm{\β}{\mathrm{\ρ}}_{\mathrm{sc}}}{2{\mathrm{\π}}^2{h}^2}\frac{p_{\mathrm{sc}}\mathrm{TZ}}{Z_{\mathrm{sc}}{T}_{\mathrm{sc}}}\frac{1}{r_w{e}^{-s}}{S}_b=(\frac{1}{k_{\mathrm{rgc}}}-1)\ln \frac{r_b}{r_w};$ P_e为地层压力，MPa，ρ为天然气相对密度；P_wf为井底流压，MPa，μ为天然气粘度，mpa·s；Z为天然气偏差系数，T_sc为地面标准温度，K；T为地层温度，K；h为有效厚度，m；r_e为气井控制半径，m，r_w为气井井筒半径，m；K为渗透率，10^-3μm²，S_c为完井表皮系数；a和b为二项式系数，S_b为凝析油阻塞表皮系数；r_b为凝析油阻塞半径，m；β为非Darcy渗流系数；K_rgc为临界凝析油饱和度下气相相对渗透率。

文献2：廖华伟，唐海，李宁等.确定不同地层压力下气井产能的方法[J].天然气技术,2008,16(4):30-32提出的方法中考虑了地层压降对气井产能的影响，通过公式表示为：

指数式： $\begin{array}{ccc}{q_{\mathrm{AOF}}}_2=\frac{c_2}{c_1}{\left(\frac{p_{e2}}{p_{e1}}\right)}^{2n}& {q_{\mathrm{AOF}}}_1\frac{Z_1{\mathrm{\μ}}_1}{Z_2{\mathrm{\μ}}_2}{\left(\frac{p_{e2}}{p_{e1}}\right)}^{2n}& {q_{\mathrm{AOF}}}_1\end{array}$

二项式： $q_{\mathrm{AOF}2}=\frac{p_{e2}}{p_{e1}}\sqrt{\frac{Z_1}{Z_2}}\·q_{\mathrm{AOF}1}$

式中：a₁、b₁和a₂、b₂分别为二项式产能方程的系数，且对应的天然气粘度、偏差系数和气相渗透率分别为μ₁、Z₁、K₁和μ₂、Z₂、K₂；P_e1和P_e2分别为不同开采阶段的地层压力。

文献3：郑丽坤，考虑渗透率应力敏感的气井无阻流量预测方法研究[J].石油地质与工程.2010,24(2).

该文章考虑了地层压力下降引起的储层渗透率应力敏感效应，通过二项式产能方程推出不同地层压力条件下气井无阻流量计算方法。

$\frac{k}{k_0}={\mathrm{be}}^{\mathrm{ap}}$

$P_R^2-P_{\mathrm{wf}}^2={\mathrm{Aq}}_g+{\mathrm{Bq}}_g^2$

$q_{\mathrm{AOF}2}\frac{P_{R2}}{P_{R1}}\sqrt{\frac{Z_1{\mathrm{\μ}}_{g1}{k}_2}{Z_2{\mathrm{\μ}}_{g2}{k}_1}}{q}_{\mathrm{AOF}1}=\frac{P_{R2}}{P_{R1}}\sqrt{\frac{Z_1{\mathrm{\μ}}_{g1}}{Z_2{\mathrm{\μ}}_{g2}}{\mathrm{be}}^{\mathrm{ap}}}{q}_{\mathrm{AOF}1}$

该方法的不足是它只适用于普通的干气气藏，不能应用于凝析气藏，因为改篇文章中渗透率的计算是基于岩石的应力敏感效应，而并非凝析油析出对储层气相渗透率产生影响。虽然都考虑了渗透率对计算结果的影响，但是这两种情况下渗透率的计算方法是有本质的不同。

由以上文献资料可知，根据气顶地层压力的变化情况，目前有两种方程确定气井的绝对无阻流量。第一，如果压力变化不大，则认为气井产能方程的系数在开发过程中是不变的，可以直接利用早期的气井产能方程；第二，考虑地层压力变化对气体粘度、偏差系数等方面的影响，通过气井产能方程推导出不同地层压力下的气井绝对无阻流量，但是目前的推导过程是基于储层渗透率不变的前提条件，而且在计算气体的粘度和偏差系数过程中，认为气体的组成是不变的，即分子量没有发生变化。这种假设对于一般的干气气藏是成立的，但是对于凝析气藏则会出现较大的计算误差。

发明内容

本发明的目的是针对上述问题，提出一种基于指数式的凝析气藏绝对无阻流量的获取方法及***，实现凝析气藏采气井绝对无阻流量的获取精度提高的效果。

为实现上述目的，本发明提供了不同地层压力下凝析气藏绝对无阻流量的获取方法，该方法基于不同地层压力下凝析气藏内气体组份的变化和气相渗透率的变化对绝对无阻流量的获取；包括：

确定不同地层压力下指数式产能方程；

根据所述指数式产能方程，由第一地层压力和第二地层压力及第一地层压力对应的第一天然气粘度、第一偏差系数和第一气相渗透率，以及第二地层压力对应的第二天然气粘度、第二偏差系数和第二气相渗透率获取指数式产能方程的第一系数和第二系数之间的关系；

根据指数式产能方程的第一系数和第二系数之间的关系、已知的第一凝析气藏绝对无阻流量以及第一地层压力和第二地层压力来获取不同地层压力下凝析气藏绝对无阻流量。

可选的，在本发明一实施例中，所述气相渗透率的获取方法包括：

根据等容衰竭实验得出不同地层压力下的凝析油含量；

根据不同地层压力下的凝析油含量和束缚水饱和度得到不同地层压力下的凝析油饱和度；

根据凝析油的析出不改变油气相渗曲线的形态，利用油气相渗曲线和不同地层压力下的凝析油饱和度得到不同含气饱和度下气相相对渗透率；

根据气相相对渗透率和空气绝对渗透率来获取气相渗透率。

可选的，在本发明一实施例中，所述根据不同地层压力下的凝析油含量和束缚水饱和度得到不同地层压力下的凝析油饱和度的步骤包括：

根据下式来获取不同地层压力下的凝析油饱和度；

S_o(p)＝[V_roCVD(p)](1-S_wi)

其中，S_o(p)为不同地层压力下的凝析油饱和度；S_wi为束缚水饱和度；V_roCVD(p)为不同地层压力下的凝析油含量。

可选的，在本发明一实施例中，所述根据气相相对渗透率和空气绝对渗透率来获取气相渗透率的步骤包括：

根据下式来获取气相渗透率；

K_rg=K/K_a

其中，K_rg为气相相对渗透率；K_a为空气绝对渗透率，由岩芯实验测得；K为气相渗透率。

可选的，在本发明一实施例中，所述指数式产能方程的第一系数和第二系数之间的关系表示为：

$\frac{c_1}{c_2}=\frac{Z_1{\mathrm{\μ}}_1{K}_2}{Z_2{\mathrm{\μ}}_2{K}_1}$

其中，c₁为第一指数式产能方程系数，c₂为第二指数式产能方程系数，Z₁为第一地层压力下对应地第一偏差系数，Z₂为第二地层压力下对应地第二偏差系数，μ₁为第一地层压力下对应地第一天然气粘度，μ₂为第二地层压力下对应地第二天然气粘度，K₁为第一地层压力下对应地第一气相渗透率，K₂为第二地层压力下对应地第二气相渗透率。

可选的，在本发明一实施例中，所述不同地层压力下凝析气藏绝对无阻流量计算公式为：

$\begin{array}{ccc}{q_{\mathrm{AOF}}}_2=\frac{c_2}{c_1}{\left(\frac{p_{e2}}{p_{e1}}\right)}^{2n}& {q_{\mathrm{AOF}}}_1\frac{Z_1{\mathrm{\μ}}_1{K}_2}{Z_2{\mathrm{\μ}}_2{K}_1}{\left(\frac{p_{e2}}{p_{e1}}\right)}^{2n}& {q_{\mathrm{AOF}}}_1\end{array}$

其中，为不同地层压力下凝析气藏绝对无阻流量，c₁为第一指数式产能方程系数，c₂为第二指数式产能方程系数，P_e1为第一地层压力，P_e2为第二地层压力，为已知的凝析气藏绝对无阻流量，μ₁为第一地层压力下对应地第一天然气粘度，μ₂为第二地层压力下对应地第二天然气粘度；K₁为第一地层压力下对应地第一气相渗透率，K₂为第二地层压力下对应地第二气相渗透率；n为产能曲线指数，n的取值范围为0.5-1，其中，n=1时，代表了完全Darcy渗流；当n=0.5时，代表了完全非Darcy渗流。

为实现上述目的，本发明还提出不同地层压力下凝析气藏绝对无阻流量的获取***，该***基于不同地层压力下凝析气藏内气体组份的变化和气相渗透率的变化对绝对无阻流量的获取；包括：

指数式产能方程确定单元，用于确定不同地层压力下指数式产能方程；

指数式系数关系确定单元，用于根据所述指数式产能方程，由第一地层压力和第二地层压力及第一地层压力对应的第一天然气粘度、第一偏差系数和第一气相渗透率，以及第二地层压力对应的第二天然气粘度、第二偏差系数和第二气相渗透率获取指数式产能方程的第一系数和第二系数之间的关系；

绝对无阻流量获取单元，用于根据指数式产能方程的第一系数和第二系数之间的关系、已知的第一凝析气藏绝对无阻流量以及第一地层压力和第二地层压力来获取不同地层压力下凝析气藏绝对无阻流量。

可选的，在本发明一实施例中，所述指数式系数关系确定单元包括：

等容衰竭模块，用于根据等容衰竭实验得出不同地层压力下的凝析油含量；

凝析油饱和度获取模块，用于根据不同地层压力下的凝析油含量和束缚水饱和度得到不同地层压力下的凝析油饱和度；

相对渗透率值获取模块，用于根据凝析油的析出不改变油气相渗曲线的形态，利用油气相渗曲线和不同地层压力下的凝析油饱和度得到不同含气饱和度下气相相对渗透率；

气相渗透率获取模块，用于根据气相相对渗透率和空气绝对渗透率来获取气相渗透率。

可选的，在本发明一实施例中，所述凝析油饱和度获取模块根据下式来获取不同地层压力下的凝析油饱和度；

S_o(p)＝[V_roCVD(p)](1-S_wi)

其中，S_o(p)为不同地层压力下的凝析油饱和度；S_wi为束缚水饱和度；V_roCVD(p)为不同地层压力下的凝析油含量。

可选的，在本发明一实施例中，所述气相渗透率获取模块根据下式来获取气相渗透率；

K_rg=K/K_a

其中，K_rg为气相相对渗透率；K_a为空气绝对渗透率，由岩芯实验测得；K为气相渗透率。

可选的，在本发明一实施例中，所述指数式系数关系确定单元获取的指数式产能方程的第一系数和第二系数之间的关系表示为：

$\frac{c_1}{c_2}=\frac{Z_1{\mathrm{\μ}}_1{K}_2}{Z_2{\mathrm{\μ}}_2{K}_1}$

其中，c₁为第一指数式产能方程系数，c₂为第二指数式产能方程系数，Z₁为第一地层压力下对应地第一偏差系数，Z₂为第二地层压力下对应地第二偏差系数，μ₁为第一地层压力下对应地第一天然气粘度，μ₂为第二地层压力下对应地第二天然气粘度，K₁为第一地层压力下对应地第一气相渗透率，K₂为第二地层压力下对应地第二气相渗透率。

可选的，在本发明一实施例中，所述绝对无阻流量获取单元获得的不同地层压力下凝析气藏绝对无阻流量计算公式为：

$\begin{array}{ccc}{q_{\mathrm{AOF}}}_2=\frac{c_2}{c_1}{\left(\frac{p_{e2}}{p_{e1}}\right)}^{2n}& {q_{\mathrm{AOF}}}_1=\frac{Z_1{\mathrm{\μ}}_1{K}_2}{Z_2{\mathrm{\μ}}_2{K}_1}{\left(\frac{p_{e2}}{p_{e1}}\right)}^{2n}& {q_{\mathrm{AOF}}}_1\end{array}$

其中，为不同地层压力下凝析气藏绝对无阻流量，c₁为第一指数式产能方程系数，c₂为第二指数式产能方程系数，P_e1为第一地层压力，P_e2为第二地层压力，为已知的凝析气藏绝对无阻流量，μ₁为第一地层压力下对应地第一天然气粘度，μ₂为第二地层压力下对应地第二天然气粘度；K₁为第一地层压力下对应地第一气相渗透率，K₂为第二地层压力下对应地第二气相渗透率；n为产能曲线指数，n的取值范围为0.5-1，其中，n=1时，代表了完全Darcy渗流；当n=0.5时，代表了完全非Darcy渗流。

上述技术方案具有如下有益效果：由于本申请提出的技术方案考虑了凝析油对气相渗透率的影响以及凝析气组份的变化，对绝对无阻流量的获取更加精准，以便指导凝析气藏的合理有效开发。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。

图1为本发明提出的基于指数式产能方程不同地层压力下凝析气藏绝对无阻流量的获取方法流程图；

图2为本发明提出的基于指数式产能方程不同地层压力下凝析气藏绝对无阻流量的获取***结构图；

图3为本发明提出的基于指数式产能方程不同地层压力下凝析气藏绝对无阻流量的获取***中指数式系数关系确定单元结构图；

图4为混合物分子量随压力的变化曲线图；

图5为气相渗透率随压力的变化曲线图；

图6为指数式凝析气藏绝对无阻流量随地层压力的变化曲线图；

图7为本发明提出的气相渗透率的获取方法流程图。

具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述。显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

凝析气藏是一种特殊、复杂的气藏。在开发过程中，凝析气井随地层压力的下降呈现不同的流动区域，不同的流动区域就伴随不同的渗流方程。要得到不同区域准确的油气分部和相应的油气相渗曲线是非常困难的，本申请在详细分析和不断尝试的基础上，采用近似的气相渗透率计算方法，取得了较好的效果，降低了计算误差，满足了工程计算精度的要求。

如图1所示，为本发明提出的基于指数式产能方程不同地层压力下凝析气藏绝对无阻流量的获取方法流程图。该方法基于不同地层压力下凝析气藏内气体组份的变化和气相渗透率的变化对绝对无阻流量的获取；包括：

步骤301：确定不同地层压力下指数式产能方程；

步骤302：根据所述指数式产能方程，由地层压力P_e1和地层压力P_e2及地层压力P_e1对应的天然气粘度μ₁、偏差系数Z₁和气相渗透率K₁，以及地层压力P_e2对应的天然气粘度μ₂、偏差系数Z₂及和气相渗透率K₂获取指数式产能方程系数c₁和c₂之间的关系；

步骤303：根据指数式产能方程系数c₁和c₂之间的关系、已知的凝析气藏绝对无阻流量以及地层压力P_e1和地层压力P_e2来获取不同地层压力下凝析气藏绝对无阻流量

可选的，在本发明一实施例中，所述指数式产能方程系数c₁和c₂之间的关系表示为：

$\frac{c_1}{c_2}=\frac{Z_1{\mathrm{\μ}}_1{K}_2}{Z_2{\mathrm{\μ}}_2{K}_1}.$

可选的，在本发明一实施例中，所述凝析气藏绝对无阻流量计算公式为：

$\begin{array}{ccc}{q_{\mathrm{AOF}}}_2=\frac{c_2}{c_1}{\left(\frac{p_{e2}}{p_{e1}}\right)}^{2n}& {q_{\mathrm{AOF}}}_1=\frac{Z_1{\mathrm{\μ}}_1{K}_2}{Z_1{\mathrm{\μ}}_2{K}_1}{\left(\frac{p_{e2}}{p_{e1}}\right)}^{2n}& {q_{\mathrm{AOF}}}_1\end{array}.$

如图7所示，为本发明提出的气相渗透率的获取方法流程图。包括：

步骤a：根据等容衰竭实验得出不同地层压力下的凝析油含量；

步骤b：根据不同地层压力下的凝析油含量和束缚水饱和度得到不同地层压力下的凝析油饱和度；

步骤c：根据凝析油的析出不改变油气相渗曲线的形态，利用油气相渗曲线和不同地层压力下的凝析油饱和度得到不同含气饱和度下气相相对渗透率；

步骤d：根据气相相对渗透率和空气绝对渗透率来获取气相渗透率。

可选的，在本发明一实施例中，所述根据不同地层压力下的凝析油含量和束缚水饱和度得到不同地层压力下的凝析油饱和度的步骤包括：

根据下式来获取不同地层压力下的凝析油饱和度；

S_o(p)＝[V_roCVD(p)](1-S_wi)

其中，S_o(p)为不同地层压力下的凝析油饱和度；S_wi为束缚水饱和度；V_roCVD(p)为不同地层压力下的凝析油含量。

可选的，在本发明一实施例中，所述根据气相相对渗透率和空气绝对渗透率来获取气相渗透率的步骤包括：

根据下式来获取气相渗透率；

K_rg=K/K_a

其中，K_rg为气相相对渗透率；K_a为空气绝对渗透率，由岩芯实验测得；K为气相渗透率。

如图2所示，为本发明提出的基于指数式产能方程不同地层压力下凝析气藏绝对无阻流量的获取***结构图。该***基于不同地层压力下凝析气藏内气体组份的变化和气相渗透率的变化对绝对无阻流量的获取；包括：

指数式产能方程确定单元401，用于确定不同地层压力下指数式产能方程；

指数式系数关系确定单元402，用于根据所述指数式产能方程，由地层压力P_e1和地层压力P_e2及地层压力P_e1对应的天然气粘度μ₁、偏差系数Z₁和气相渗透率K₁，以及地层压力P_e2对应的天然气粘度μ₂、偏差系数Z₂及和气相渗透率K₂获取指数式产能方程系数c₁和c₂之间的关系；

绝对无阻流量获取单元403，用于根据指数式产能方程系数c₁和c₂之间的关系、已知的凝析气藏绝对无阻流量以及地层压力P_e1和地层压力P_e2来获取不同地层压力下凝析气藏绝对无阻流量

如图3所示，为本发明提出的基于指数式产能方程不同地层压力下凝析气藏绝对无阻流量的获取***中指数式系数关系确定单元结构图。所述指数式系数关系确定单元402包括：

等容衰竭模块4021，用于根据等容衰竭实验得出不同地层压力下的凝析油含量；

凝析油饱和度获取模块4022，用于根据不同地层压力下的凝析油含量和束缚水饱和度得到不同地层压力下的凝析油饱和度；

相对渗透率值获取模块4023，用于根据凝析油的析出不改变油气相渗曲线的形态，利用油气相渗曲线和不同地层压力下的凝析油饱和度得到不同含气饱和度下气相相对渗透率；

气相渗透率获取模块4024，用于根据气相相对渗透率和空气绝对渗透率来获取气相渗透率。

可选的，在本发明一实施例中，所述凝析油饱和度获取模块4022根据下式来获取不同地层压力下的凝析油饱和度；

S_o(p)＝[V_roCVD(p)](1-S_wi)

其中，S_o(p)为不同地层压力下的凝析油饱和度；S_wi为束缚水饱和度；V_roCVD(p)为不同地层压力下的凝析油含量。

可选的，在本发明一实施例中，所述气相渗透率获取模块4024根据下式来获取气相渗透率；

K_rg=K/K_a

其中，K_rg为气相相对渗透率；K_a为空气绝对渗透率，由岩芯实验测得；K为气相渗透率。

可选的，在本发明一实施例中，所述指数式系数关系确定单元402获取的所述指数式产能方程系数c₁和c₂之间的关系表示为：

$\frac{c_1}{c_2}=\frac{Z_1{\mathrm{\μ}}_1{K}_2}{Z_2{\mathrm{\μ}}_2{K}_1}.$

可选的，在本发明一实施例中，所述绝对无阻流量获取单元403获得的凝析气藏绝对无阻流量计算公式为：

$\begin{array}{ccc}{q_{\mathrm{AOF}}}_2=\frac{c_2}{c_1}{\left(\frac{p_{e2}}{p_{e1}}\right)}^{2n}& {q_{\mathrm{AOF}}}_1=\frac{Z_1{\mathrm{\μ}}_1{K}_2}{Z_2{\mathrm{\μ}}_2{K}_1}{\left(\frac{p_{e2}}{p_{e1}}\right)}^{2n}& {q_{\mathrm{AOF}}}_1\end{array}.$

本发明提出的装置对于利用通用处理器，数字信号处理器，专用集成电路（ASIC），现场可编程门阵列（FPGA）或其它可编程逻辑装置，离散门或晶体管逻辑，离散硬件部件，计算机装置，或上述任何组合的设计来实现的配电网自动化装置中线路负荷分配，不需要主站及子站及复杂的通讯装置。

本领域技术人员还可以了解到本发明的装置中列出的各种说明性逻辑块（illustrativelogicalblock）、单元和步骤可以通过电子硬件、电脑软件，或两者的结合进行实现。为清楚展示硬件和软件的可替换性（interchangeability），上述的各种说明性部件（illustrativecomponents）、单元和步骤已经通用地描述了它们的功能。这样的功能是通过硬件还是软件来实现取决于特定的应用和整个装置的设计要求。本领域技术人员可以对于每种特定的应用，可以使用各种方法实现所述的功能，但这种实现不应被理解为超出本发明实施例保护的范围。

本发明的装置中所描述的各种说明性的逻辑块，或单元都可以通过通用处理器，数字信号处理器，专用集成电路（ASIC），现场可编程门阵列（FPGA）或其它可编程逻辑装置，离散门或晶体管逻辑，离散硬件部件，或上述任何组合的设计来实现或操作所描述的功能。通用处理器可以为微处理器，可选地，该通用处理器也可以为任何传统的处理器、控制器、微控制器或状态机。处理器也可以通过计算装置的组合来实现，例如数字信号处理器和微处理器，多个微处理器，一个或多个微处理器联合一个数字信号处理器核，或任何其它类似的配置来实现。

本发明的装置中所描述的方法或算法的步骤可以直接嵌入硬件、处理器执行的软件模块、或者这两者的结合。软件模块可以存储于RAM存储器、闪存、ROM存储器、EPROM存储器、EEPROM存储器、寄存器、硬盘、可移动磁盘、CD-ROM或本领域中其它任意形式的存储媒介中。示例性地，存储媒介可以与处理器连接，以使得处理器可以从存储媒介中读取信息，并可以向存储媒介存写信息。可选地，存储媒介还可以集成到处理器中。处理器和存储媒介可以设置于ASIC中，ASIC可以设置于用户终端中。可选地，处理器和存储媒介也可以设置于用户终端中的不同的部件中。

在一个或多个示例性的设计中，本发明的装置中所描述的上述功能可以在硬件、软件、固件或这三者的任意组合来实现。如果在软件中实现，这些功能可以存储与电脑可读的媒介上，或以一个或多个指令或代码形式传输于电脑可读的媒介上。电脑可读媒介包括电脑存储媒介和便于使得让电脑程序从一个地方转移到其它地方的通信媒介。存储媒介可以是任何通用或特殊电脑可以接入访问的可用媒体。例如，这样的电脑可读媒体可以包括但不限于RAM、ROM、EEPROM、CD-ROM或其它光盘存储、磁盘存储或其它磁性存储装置，或其它任何可以用于承载或存储以指令或数据结构和其它可被通用或特殊电脑、或通用或特殊处理器读取形式的程序代码的媒介。此外，任何连接都可以被适当地定义为电脑可读媒介，例如，如果软件是从一个网站站点、服务器或其它远程资源通过一个同轴电缆、光纤电缆、双绞线、数字用户线（DSL）或以例如红外、无线和微波等无线方程传输的也被包含在所定义的电脑可读媒介中。所述的碟片（disk）和磁盘（disc）包括压缩磁盘、镭射盘、光盘、DVD、软盘和蓝光光盘，磁盘通常以磁性复制数据，而碟片通常以激光进行光学复制数据。上述的组合也可以包含在电脑可读媒介中。

对于凝析气藏，当地层压力下降到露点压力以下后，凝析气中的重质组分会逐渐反凝析到地层中，造成凝析气的组成会发生变化，同时反凝析出来的凝析油会吸附岩石表面，影响储层渗流动能力。本发明实施例中探讨了不同地层压力下凝析气藏绝对无阻流量获取的技术方案，通过研究地层压力下降对气相渗透率和凝析气组成的影响，利用指数式产能方程，可以推导出凝析气藏在不同地层压力下的气井产能方程，以便指导凝析气藏的合理有效开发。

对于凝析气藏与干气气藏的主要差别在于以下两个方面：第一，气体组份的变化。如图4所示，为混合物分子量随压力的变化曲线图；第二，气相渗透率的变化。如图5所示，为气相渗透率随压力的变化曲线图。其中，气体组成变化对产能的影响主要表现在气体混合物分子量的变化上。根据已有试气数据，利用数学方法便可以得出不同地层压力下气体的分子量。当凝析油从气体中析出时，气相的渗透率会发生改变，根据实验得出的相渗曲线便可以得出不同含气饱和度下的气相渗透率。

在本申请中，不同地层压力下气相渗透率的获取方法：

气相渗透率K的改变是由于凝析油的析出对气相流动产生影响，这是凝析气藏的特定性质决定的，而非应力敏感引起的渗透率改变。

根据CVD实验可以得出不同地层压力下的凝析油含量，这时岩石孔隙中会出现油、气、水三相，束缚水饱和度S_wi可以由实测数据得到，由公式S_o(p)＝[V_roCVD(p)](1-S_wi)可以得到不同地层压力下的凝析油饱和度。假设凝析油的析出不改变油气相渗曲线的形态，根据油气相渗曲线和不同地层压力下的凝析油饱和度，可以得到不同含气饱和度下油相和气相的相对渗透率值。

在现有的众多获取气相渗透率的方法中，对于凝析气藏来说具有如下特点：

A）：很难得到实际地层状态下油气相渗曲线；

B）：根据实际的岩心资料，在实验室内可以得到油气相渗曲线；

C）：该相渗曲线反映了油气两相渗流特点，再由CVD实验得出不同地层压力下的凝析油含量，就可以近似计算出不同地层压力下的气相渗透率。总之，该方法可以利用有限的条件得到近似于真实地层状态下油气两相的流动效果，大幅降低了凝析气藏气井无阻流量的计算误差。

实施例：

根据指数式产能方程，当井底流压为大气压时，有：

$q_{\mathrm{AOF}}={\mathrm{cP}}_e^{2n}$

同理假设在两个不同的开采阶段，地层压力分别为P_e1和P_e2，对应的产能方程系数分别为c₁和c₂，可以得到：

$\begin{array}{ccc}{q_{\mathrm{AOF}}}_2=\frac{c_2}{c_1}{\left(\frac{p_{e2}}{p_{e1}}\right)}^{2n}& {q_{\mathrm{AOF}}}_1=\frac{Z_1{\mathrm{\μ}}_1{K}_2}{Z_2{\mathrm{\μ}}_2{K}_1}{\left(\frac{p_{e2}}{p_{e1}}\right)}^{2n}& {q_{\mathrm{AOF}}}_1\end{array}$

通常情况下，气井流动指数n的变化是可以忽略的，但是当气井进行了酸化或者压裂等措施以后，方程的流动指数n会有较大的变化，这时必须通过重新试井建立产能方程。

如图6所示，为指数式凝析气藏绝对无阻流量随地层压力的变化曲线图。图6中的A曲线没有考虑凝析油对储层渗透率和气体分子量的影响，B曲线同时考虑地层压降对气体分子量和储层渗透率的影响。从图中可以看出，当地层压力降到露点压力以下时，考虑凝析油对气相渗透率的影响以及凝析气组份的变化，计算结果更为准确。

由下表1中计算结果值可知，依据某油田的实际试气数据和该井的检测分析报告，采用通常的干气气藏计算方法和基于指数式的不同地层压力下凝析气藏计算方法分别计算气井的绝对无阻流量。其计算结果显示，采用干气气藏的计算方法来计算凝析气藏的气井产量都会出现较大的误差，超过了工程计算的精度要求（0-10%），而基于指数式的不同地层压力下凝析气藏计算方法则大大降低了计算误差，满足了工程计算的精度要求，故该方法可以很好的指导油田的生产。

表12007年Q1井指数式的计算结果和误差

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

Claims (10)

1.不同地层压力下凝析气藏绝对无阻流量的获取方法，其特征在于，该方法基于不同地层压力下凝析气藏内气体组份的变化和气相渗透率的变化对绝对无阻流量的获取；包括：

确定不同地层压力下指数式产能方程；

根据所述指数式产能方程，由第一地层压力和第二地层压力及第一地层压力对应的第一天然气粘度、第一偏差系数和第一气相渗透率，以及第二地层压力对应的第二天然气粘度、第二偏差系数和第二气相渗透率获取指数式产能方程的第一系数和第二系数之间的关系；

根据指数式产能方程的第一系数和第二系数之间的关系、已知的第一凝析气藏绝对无阻流量以及第一地层压力和第二地层压力来获取不同地层压力下凝析气藏绝对无阻流量；

其中，所述气相渗透率的获取方法包括：

根据等容衰竭实验得出不同地层压力下的凝析油含量；

根据不同地层压力下的凝析油含量和束缚水饱和度得到不同地层压力下的凝析油饱和度；

根据凝析油的析出不改变油气相渗曲线的形态，利用油气相渗曲线和不同地层压力下的凝析油饱和度得到不同含气饱和度下气相相对渗透率；

根据气相相对渗透率和空气绝对渗透率来获取气相渗透率。

2.根据权利要求1所述的方法，其特征在于，所述根据不同地层压力下的凝析油含量和束缚水饱和度得到不同地层压力下的凝析油饱和度的步骤包括：

根据下式来获取不同地层压力下的凝析油饱和度；

S_o(p)＝[V_roCVD(p)](1-S_wi)

其中，S_o(p)为不同地层压力下的凝析油饱和度；S_wi为束缚水饱和度；V_roCVD(p)为不同地层压力下的凝析油含量。

3.根据权利要求1所述的方法，其特征在于，所述根据气相相对渗透率和空气绝对渗透率的来获取气相渗透率的步骤包括：

根据下式来获取气相渗透率；

K_rg＝K/K_a

其中，K_rg为气相相对渗透率；K_a为空气绝对渗透率，由岩芯实验测得；K为气相渗透率。

4.根据权利要求1～3任一权利要求所述的方法，其特征在于，所述指数式产能方程的第一系数和第二系数之间的关系表示为：

$\frac{c_1}{c_2}=\frac{Z_1{\mathrm{\μ}}_1{K}_2}{Z_2{\mathrm{\μ}}_2{K}_1}$

其中，c₁为第一指数式产能方程系数，c₂为第二指数式产能方程系数，Z₁为第一地层压力下对应的第一偏差系数，Z₂为第二地层压力下对应的第二偏差系数，μ₁为第一地层压力下对应的第一天然气粘度，μ₂为第二地层压力下对应的第二天然气粘度，K₁为第一地层压力下对应的第一气相渗透率，K₂为第二地层压力下对应的第二气相渗透率。

5.根据权利要求1～3任一权利要求所述的方法，其特征在于，所述不同地层压力下凝析气藏绝对无阻流量计算公式为：

$q_{{\mathrm{AOF}}_2}=\frac{c_2}{c_1}{\left(\frac{p_{e2}}{p_{e1}}\right)}^{2n},q_{{\mathrm{AOF}}_1}=\frac{Z_1{\mathrm{\μ}}_1{K}_2}{Z_2{\mathrm{\μ}}_2{K}_1}{\left(\frac{p_{e2}}{p_{e1}}\right)}^{2n}{q}_{{\mathrm{AOF}}_1}$

其中，为不同地层压力下凝析气藏绝对无阻流量，c₁为第一指数式产能方程系数，c₂为第二指数式产能方程系数，Ρ_e1为第一地层压力，Ρ_e2为第二地层压力，为已知的凝析气藏绝对无阻流量，μ₁为第一地层压力下对应的第一天然气粘度，μ₂为第二地层压力下对应的第二天然气粘度；K₁为第一地层压力下对应的第一气相渗透率，K₂为第二地层压力下对应的第二气相渗透率；n为产能曲线指数，n的取值范围为0.5-1。

6.不同地层压力下凝析气藏绝对无阻流量的获取***，其特征在于，该***基于不同地层压力下凝析气藏内气体组份的变化和气相渗透率的变化对绝对无阻流量的获取；包括：

指数式产能方程确定单元，用于确定不同地层压力下指数式产能方程；

指数式系数关系确定单元，用于根据所述指数式产能方程，由第一地层压力和第二地层压力及第一地层压力对应的第一天然气粘度、第一偏差系数和第一气相渗透率，以及第二地层压力对应的第二天然气粘度、第二偏差系数和第二气相渗透率获取指数式产能方程的第一系数和第二系数之间的关系；

绝对无阻流量获取单元，用于根据指数式产能方程的第一系数和第二系数之间的关系、已知的第一凝析气藏绝对无阻流量以及第一地层压力和第二地层压力来获取不同地层压力下凝析气藏绝对无阻流量；

其中，所述指数式系数关系确定单元包括：

等容衰竭模块，用于根据等容衰竭实验得出不同地层压力下的凝析油含量；

凝析油饱和度获取模块，用于根据不同地层压力下的凝析油含量和束缚水饱和度得到不同地层压力下的凝析油饱和度；

相对渗透率值获取模块，用于根据凝析油的析出不改变油气相渗曲线的形态，利用油气相渗曲线和不同地层压力下的凝析油饱和度得到不同含气饱和度下气相相对渗透率；

气相渗透率获取模块，用于根据气相相对渗透率和空气绝对渗透率来获取气相渗透率。

7.根据权利要求6所述的***，其特征在于，所述凝析油饱和度获取模块根据下式来获取不同地层压力下的凝析油饱和度；

S_o(p)＝[V_roCVD(p)](1-S_wi)

其中，S_o(p)为不同地层压力下的凝析油饱和度；S_wi为束缚水饱和度；V_roCVD(p)为不同地层压力下的凝析油含量。

8.根据权利要求6所述的***，其特征在于，所述气相渗透率获取模块根据下式来获取气相渗透率；

K_rg＝K/K_a

其中，K_rg为气相相对渗透率；K_a为空气绝对渗透率，由岩芯实验测得；K为气相渗透率。

9.根据权利要求6～8任一权利要求所述的***，其特征在于，所述指数式系数关系确定单元获取的指数式产能方程的第一系数和第二系数之间的关系表示为：

$\frac{c_1}{c_2}=\frac{Z_1{\mathrm{\μ}}_1{K}_2}{Z_2{\mathrm{\μ}}_2{K}_1}$

其中，c₁为第一指数式产能方程系数，c₂为第二指数式产能方程系数，Z₁为第一地层压力下对应的第一偏差系数，Z₂为第二地层压力下对应的第二偏差系数，μ₁为第一地层压力下对应的第一天然气粘度，μ₂为第二地层压力下对应的第二天然气粘度，K₁为第一地层压力下对应的第一气相渗透率，K₂为第二地层压力下对应的第二气相渗透率。

10.根据权利要求6～8任一权利要求所述的***，其特征在于，所述绝对无阻流量获取单元获得的不同地层压力下凝析气藏绝对无阻流量计算公式为：

$q_{{\mathrm{AOF}}_2}=\frac{c_2}{c_1}{\left(\frac{p_{e2}}{p_{e1}}\right)}^{2n},q_{{\mathrm{AOF}}_1}=\frac{Z_1{\mathrm{\μ}}_1{K}_2}{Z_2{\mathrm{\μ}}_2{K}_1}{\left(\frac{p_{e2}}{p_{e1}}\right)}^{2n}{q}_{{\mathrm{AOF}}_1}$

其中，为不同地层压力下凝析气藏绝对无阻流量，c₁为第一指数式产能方程系数，c₂为第二指数式产能方程系数，Ρ_e1为第一地层压力，Ρ_e2为第二地层压力，为已知的凝析气藏绝对无阻流量，μ₁为第一地层压力下对应的第一天然气粘度，μ₂为第二地层压力下对应的第二天然气粘度；K₁为第一地层压力下对应的第一气相渗透率，K₂为第二地层压力下对应的第二气相渗透率；n为产能曲线指数，n的取值范围为0.5-1。