CN104392068B

CN104392068B - 一种含一条直线封闭边界油藏曲线拟合的处理方法

Info

Publication number: CN104392068B
Application number: CN201410758423.XA
Authority: CN
Inventors: 王磊; 何桃; 程翔; 赵锦波; 胡丽琴
Original assignee: Hefei Institutes of Physical Science of CAS
Current assignee: Hefei Institutes of Physical Science of CAS
Priority date: 2014-12-10
Filing date: 2014-12-10
Publication date: 2017-09-29
Anticipated expiration: 2034-12-10
Also published as: CN104392068A

Abstract

本发明公开了一种含一条直线封闭边界油藏曲线拟合的处理方法，该方法可基于含一条直线封闭边界均质油藏图版特征，根据坐标映射关系利用鼠标位置获得到实测导数曲线第一个水平段终点位置，并以此计算获得对应的L_De^S参数；根据获得的L_De^S参数计算并绘制一理论曲线及对应的理论导数曲线；达到了控制光标一次即可实现曲线的快速拟合的目的，大大的提高了工作的效率。

Description

一种含一条直线封闭边界油藏曲线拟合的处理方法

技术领域

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

背景技术

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

其中，现有曲线拟合法的本质是通过试凑法实现曲线拟合；即，不断的更改参数，每更改一次后，查看拟合结果，直至曲线的外形与弧度近似；再移动曲线使其重合。该方法的缺陷主要为：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为图在屏幕上的宽度；

由峰值t_D/C_D计算出对应的L_De^S参数，其公式为：

本发明与现有技术相比的优点在于：由上述本发明提供的技术方案可以看出，基于坐标映射关系，从而利用鼠标位置获得实测导数曲线第一个水平段的终点坐标，并根据该坐标计算对应的L_De^S参数；根据获得的L_De^S参数计算并绘制一理论曲线及对应的理论导数曲线；达到了控制光标一次即可得到正确的L_De^S，实现了曲线快速拟合的目的，大大提高了工作效率。与试凑的曲线拟合方法相比，本方法根据含一条直线封闭边界的均质油藏图版特征，人工获取实测导数曲线第一个水平段的终点坐标，计算相应参数，达到一次拟合成功的效果。

附图说明

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

图1为本发明实施例提供的一种曲线拟合的处理方法的流程图；

图2为本发明实施例提供的双对数图中绘制辅助线的示意图；

图3为本发明实施例提供的双对数图中绘制实测曲线及实测导数曲线的示意图；

图4为本发明实施例提供的实测导数曲线峰值和理论导数曲线峰值重合的示意图；

图5为本发明实施例提供的曲线拟合过程的示意图；

图6为本发明实施例提供的完成曲线拟合后的示意图。

具体实施方式

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

图1为本发明实施例提供的一种曲线拟合的处理方法的流程图。如图1所示，该方法主要包括如下步骤：

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

如图2所示，在双对数图(例如，横坐标为时间，其范围为10^-2～10⁴，纵坐标为压力，其范围为10^-2～10²)上绘三条辅助线：一条斜率为预定值s(例如，1)的斜线，另一条数值为预定值m(例如，0.5)的水平线，第三条数值为预定值n(例如，1)的水平线。

所述实测曲线可以是井底压力数据曲线，对应的实测导数曲线为该压力数据求导后的曲线。如图3所示，将两条曲线绘制在双对数图上，并使得所述实测曲线及实测导数曲线的初始部分，与经过原点且斜率为s的斜线重合；该实测导数曲线的第一个水平段(图3虚线方框部分)与数值为m的水平线重合；该实测导数曲线的结尾部分与数值为n的水平线重合；即通过辅助线将上述两条曲线的位置进行标定，便于后期的曲线拟合工作。

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

如图4所示，输入一个无量纲组合参数C_De^2S的值(例如500)，在已绘制实测曲线和实测导数曲线的双对数图上绘制理论曲线与理论导数曲线，根据实测导数曲线峰值与理论导数曲线峰值的相对位置调整C_De^2S的值，使实测导数曲线峰值与理论导数曲线峰值重合。

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

由于本实施例的曲线拟合基于无量纲组合参数L_De^S，且导数曲线形态对L_De^S参数较为敏感，因此需要计算出实测导数曲线的L_De^S参数；该参数与实测导数曲线第一个水平段的终点位置存在对应关系，当获知曲线第一个水平段的终点位置后，可计算获得对应的L_De^S参数。同理，当得知实测导数曲线第一个水平段终点位置对应的L_De^S参数后，则可计算出与实测曲线及对应的实测导数曲线完全重合的理论曲线及对应的理论导数曲线。

下面针对实测导数曲线的无量纲组合参数L_De^S参数的求解过程做详细介绍：

示例性的，若实测曲线为含一条直线封闭边界地层的井底压力数据曲线，则无量纲井底压力表达式为：

求解所述无量纲井底压力表达式获得第一个水平段终点位置对应的时间为：

t_D/C_D＝a(L_De^S)²/(C_De^2S)； (2)

使用大量数据对(2)进行回归后可以得到：

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

由(4)式可知，当得知实测导数曲线的C_De^2S和第一个水平段的终点时间t_D/C_D后，则可计算出无量纲组合参数L_De^S参数。本发明实施例基于坐标映射关系，利用鼠标位置计算实测导数曲线的t_D/C_D，同时使用步骤12输入的C_De^2S，代入(4)计算出对应的L_De^S参数。具体方法如下：首先，接收光标控制指令，根据该指令控制并移动光标(即用户点击光标并保持点击状态移动光标时，计算机根据接受到的指令进行相应的操作)，当光标移动至该实测导数曲线第一个水平段的终点位置时，利用坐标映射关系计算该处的t_D/C_D，其公式为：其中，x为光标在屏幕中的横坐标，x₀为双对数图的坐标原点在屏幕中的横坐标，l_x为图在屏幕上的宽度。将t_D/C_D和步骤12获得的C_De^2S代入(4)计算出无量纲组合参数L_De^S。

另一方面，为了便于实时观察曲线拟合的结果，在光标移动至第一个水平段终点之前，也可以实时计算理论曲线及理论导数曲线，光标变为左右箭头提示用户左右移动，如图5所示。这样，用户就可以根据理论曲线及理论导数曲线与实测曲线及实测导数曲线的拟合程度判断，光标的位置是否为实测导数曲线第一个水平段终点处。

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

如图6所示，当光标移动至该实测导数曲线第一个水平段终点处时，根据步骤12输入的C_De^2S和步骤13计算到的L_De^S参数计算并绘制一理论曲线及对应的理论导数曲线；此时，四条曲线完全重合，完成曲线拟合，可以终止对光标的点击操作，通过上述方式实现了移动光标一次即得到L_De^S并完成曲线拟合的目的，降低了工作的难度及时间，提高了工作的效率。

本发明实施例基于坐标映射关系，从而利用鼠标位置获得到实测导数曲线第一个水平段的终点位置，并以此计算获得对应的L_De^S参数；根据获得的L_De^S参数计算并绘制一理论曲线及对应的理论导数曲线；达到了控制光标一次即可实现曲线的快速拟合的目的，大大的提高了工作的效率。

通过以上的实施方式的描述，本领域的技术人员可以清楚地了解到上述实施例可以通过软件实现，也可以借助软件加必要的通用硬件平台的方式来实现。基于这样的理解，上述实施例的技术方案可以以软件产品的形式体现出来，该软件产品可以存储在一个非易失性存储介质(可以是CD-ROM，U盘，移动硬盘等)中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述的方法。

以上所述，仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明披露的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应该以权利要求书的保护范围为准。

Claims

1.一种含一条直线封闭边界油藏曲线拟合的处理方法，其特征在于，该方法包括：

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

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

根据计算到的参数计算并绘制一理论曲线及对应的理论导数曲线，完成曲线拟合；

所述实测曲线包括：

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

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

；

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

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

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

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

由实测导数曲线第一个水平段终点对应的无量纲参数计算出对应的参数，其公式为：。