[发明专利]一种含一条直线封闭边界油藏曲线拟合的处理方法

说明书

技术领域

本发明涉及计算机技术和石油工程领域，尤其涉及一种含一条直线封闭边界油藏曲线拟合的处理方法。

背景技术

现有的试井分析方法主要包括：直线段分析方法、曲线拟合法及非线性回归方法。

其中，现有曲线拟合法的本质是通过试凑法实现曲线拟合；即，不断的更改参数，每更改一次后，查看拟合结果，直至曲线的外形与弧度近似；再移动曲线使其重合。该方法的缺陷主要为：1)它是试凑方法，对解释工程师水平有很高要求；2)需要多次反复的更改参数，并进行相应计算，工作效率较低；3)当参数范围较大时更改的次数非常多，无法快速的进行曲线拟合。

发明内容

本发明的技术解决问题：克服现有技术的不足，提供一种含一条直线封闭边界油藏曲线拟合的处理方法，提高工作效率。

本发明的目的是通过以下技术方案实现的：

一种曲线拟合的处理方法，该方法包括：

将实测曲线及对应的实测导数曲线绘制在双对数图版上，使得所述实测曲线及实测导数曲线的初始部分，与经过原点且斜率为预定值s的斜线重合；该实测导数曲线的第一个水平段，与预定值为m的水平线重合；该实测导数曲线的结尾部分与数值为预定值n的水平线重合；

绘制一理论曲线及对应的理论导数曲线，输入无量纲组合参数C_De^2S的值，使理论导数曲线与实测导数曲线峰值重合；

接收光标控制指令，根据该指令控制并移动光标，当光标移动至该实测导数曲线第一个水平段的终点时，计算该位置对应的无量纲组合参数L_De^S；

根据计算出的L_De^S参数计算并绘制一理论曲线及对应的理论导数曲线，完成曲线拟合。

所述实测曲线包括：

含一条直线封闭边界地层的井底压力数据曲线，其无量纲井底压力表达式为：

其中，J₀(u)与J₁(u)为0阶与1阶第一类Bessel函数(贝塞尔函数)；Y₀(u)与Y₁(u)为0阶与1阶第二类Bessel函数；为无量纲压力表达式；为无量纲时间；为无量纲井筒存储常数；S为表皮系数；为无量纲直线封闭边界距离；p_wf为井底流动压力(MPa)；p_i为原始地层压力(MPa)；k为地层渗透率(D)；h为地层有效厚度(m)；q为井产量(m³/d)；μ为流体粘度(mPa.s)；t为流动时间(h)；r_w为井筒半径(m)；φ为地层孔隙度；C_t为综合压缩系数(1/Mpa)；C为井筒存储参数(m³/MPa)；L为直线封闭边界距离(m)；Ei为指数积分函数。

无量纲参数L_De^S与实测导数曲线第一个水平段终点对应的t_D/C_D的关系为：

t_D/C_D＝0.3248(L_De^S)²/(C_De^2S)。

所述当光标移动至该实测导数曲线第一个水平段的终点时，计算该位置对应的无量纲组合参数L_De^S包括：

当光标移动至该实测导数曲线第一个水平段的终点处时，根据光标在屏幕中的位置计算所述实测导数曲线第一个水平段终点对应的无量纲参数t_D/C_D及所述实测导数曲线对应的L_De^S参数；

所述根据光标在屏幕中的位置计算所述实测导数曲线第一个水平段终点对应的无量纲参数t_D/C_D及所述实测导数曲线对应的L_De^S参数包括：

根据光标在屏幕中的位置计算所述实测导数曲线该实测导数曲线第一个水平段的终点对应的t_D/C_D，其公式为：其中，x为光标在屏幕中的横坐标，x₀为双对数图的坐标原点在屏幕中的横坐标，l_x为图在屏幕上的宽度；