[发明专利]一种高阶煤岩气相渗透率动态变化的预测方法

权利要求书

1.一种高阶煤岩气相渗透率动态变化的预测方法，其特征在于，包括以下步骤：

在实验室条件下，设置实验温度T为储层温度；

确定原位储层应力大小；

在所述实验温度T和所述原位储层应力条件下，测试n个气体压力下的渗透率k_i数据，记做

在所述气体压力和所述实验温度T下，确定气体密度ρ_i和粘度μ_i值；

计算气体压力为时的渗透率f_i；

$f_i=\frac{2{\mathrm{\μ}}_{i}M(a{\stackrel{\‾}{P}}_i+b)}{3\×{10}^3{\mathrm{RT\ρ}}_i^2}{\left(\frac{8RT}{\mathrm{\π}M}\right)}^{0.5}+\[1+{\left(\frac{8\mathrm{\π}RT}{M}\right)}^{0.5}\frac{{\mathrm{\μ}}_i}{(a{\stackrel{\‾}{P}}_i+b){\stackrel{\‾}{P}}_i}(\frac{2}{c{\stackrel{\‾}{P}}_i+d}-1)\]\frac{{(a{\stackrel{\‾}{P}}_i+b)}^2}{8{\mathrm{\ρ}}_i},i=1...n$

式中，M为分子质量；为平均气体压力；R为气体常数；T为温度；μ_i为压力为时的粘度；ρ_i为压力为时的密度；a、b、c和d为拟合系数；

预测气体压力降低过程中煤岩气相渗透率k_g动态变化。

2.如权利要求1所述的高阶煤岩气相渗透率动态变化的预测方法，其特征在于：运用最小二乘法，求解出a、b、c和d。

3.如权利要求1-2任一项所述的高阶煤岩气相渗透率动态变化的预测方法，其特征在于，煤岩气相渗透率k_g表述为：

$k_g=\frac{2\mathrm{\μ}M(a\stackrel{\‾}{P}+b)}{3\×{10}^3{\mathrm{RT\ρ}}^2}{\left(\frac{8RT}{\mathrm{\π}M}\right)}^{0.5}+\[1+{\left(\frac{8\mathrm{\π}RT}{M}\right)}^{0.5}\frac{\mathrm{\μ}}{(a\stackrel{\‾}{P}+b)\stackrel{\‾}{P}}(\frac{2}{c\stackrel{\‾}{P}+d}-1)\]\frac{{(a\stackrel{\‾}{P}+b)}^2}{8\mathrm{\ρ}}$

式中，M为分子质量；为平均气体压力；R为气体常数；T为温度；μ为压力为时的粘度；ρ为压力为时的密度；a、b、c和d为拟合系数。