[发明专利]高压致密气藏射孔总表皮和压实带表皮系数模拟计算方法

权利要求书

1.一种高压致密气藏射孔总表皮系数模拟计算方法，其特征在于，包括以下步骤：

S1：采用局部加密网格，构建有射孔压实带的地质模型一；

S2：为所述地质模型一的每一网格单元赋储层和流体物性的值；

S3：建立气体稳定渗流控制方程，射孔孔眼作为内边界条件，设为井底流压，地质模型径向流半径外的网格单元作为外边界条件，设为地层压力；采用有限体积法对所述气体稳定渗流控制方程进行离散化，相邻网格间采用Two-Point-Flux-Approximation方法对气体流动能力进行调和平均，得到线性方程组；

S4：采用代数多重网格方法对所述线性方程组进行求解，得到所述地质模型一的稳态压力场分布；

S5：根据所述稳态压力场分布计算射孔井流量；

S6：根据所述射孔井流量，计算所述地质模型一的总表皮系数。

2.根据权利要求1所述的高压致密气藏射孔总表皮系数模拟计算方法，其特征在于，所述储层和流体物性包括渗透率、孔隙度和动力粘度，压实带的渗透率和孔隙度均低于储层的渗透率和孔隙度。

3.根据权利要求1所述的高压致密气藏射孔总表皮系数模拟计算方法，其特征在于，在地层压力大于21MPa的情况下，所述气体稳定渗流控制方程具体为：

式中：

为劈形算符；

K为气体有效渗透率，mD；

p为压力，MPa。

4.根据权利要求1所述的高压致密气藏射孔总表皮系数模拟计算方法，其特征在于，所述射孔井流量的计算方法为：

首先，通过达西公式计算得到每个孔眼网格单元边界的速度，所述达西公式具体为：

式中：

u为速度，m/s；

μ为孔眼处天然气动力粘度，mPa·s；

然后，将每个孔眼网格的速度乘以边界的面积得到每个孔眼网格的流量；

最后，计算所有流量之和，即可得到所述射孔井流量。

5.根据权利要求1-4中任意一项所述的高压致密气藏射孔总表皮系数模拟计算方法，其特征在于，所述总表皮系数的计算公式为：

式中：

S为总表皮系数，无量纲；

Δp为地层与孔眼间的压差，MPa；

为储层平均渗透率，mD；

h为储层厚度，m；

q为射孔井流量，m³/s；

为平均地层压力下天然气动力粘度，mPa·s；

r_e为径向流半径，m；

r_w为井筒半径，m。

6.一种高压致密气藏射孔压实带表皮系数模拟计算方法，其特征在于，包括以下步骤：

利用权利要求1-5中任意一所述的高压致密气藏射孔总表皮系数模拟计算方法，计算得到地质模型一的总表皮系数；

将所述地质模型一的压实系数设为1，其余赋值情况不变，即可得到无射孔压实带的地质模型二；

然后采用与所述地质模型一的总表皮系数相同的计算方法，计算得到地质模型二的总表皮系数；

计算所述地质模型一的总表皮系数与所述地质模型二的总表皮系数之差，即可得到所述射孔压实带表皮系数。