一种曲线拟合的处理方法技术

本发明专利技术公开了一种曲线拟合的处理方法，该方法可基于均质油藏图版特征，根据坐标映射关系利用鼠标位置获得到实测导数曲线的峰值位置，并以此计算获得对应的CDe2S参数；根据获得的CDe2S参数计算并绘制一理论曲线及对应的理论导数曲线；达到了控制光标一次即可实现曲线的快速拟合的目的，大大的提高了工作的效率。

【技术实现步骤摘要】
一种曲线拟合的处理方法
本专利技术涉及计算机技术和石油工程领域，尤其涉及一种曲线拟合的处理方法。
技术介绍
现有的试井分析方法主要包括：直线段分析方法、曲线拟合法及非线性回归方法。其中，现有曲线拟合法的本质是通过试凑法实现曲线拟合；即，不断的更改参数，每更改一次后，查看拟合结果，直至曲线的外形与弧度近似；再移动曲线使其重合。该方法的缺陷主要为：1）它是试凑方法，对解释工程师水平有很高要求；2）需要多次反复的更改参数，并进行相应计算，工作效率较低；3）当参数范围较大时（相关参数的范围10-6～10200）更改的次数非常多，无法快速的进行曲线拟合。
技术实现思路
本专利技术的目的是提供一种曲线拟合的处理方法，提高了工作效率。本专利技术的目的是通过以下技术方案实现的：一种曲线拟合的处理方法，该方法包括：将实测曲线及对应的实测导数曲线绘制在双对数图上，使得所述实测曲线及实测导数曲线的初始部分，与经过原点且斜率为预定值s的斜线重合；且该实测导数曲线的结尾部分与数值为预定值n的水平线重合；接收光标控制指令，根据该指令控制并移动光标，当光标移动至该实测导数曲线的峰值处时，计算该实测导数曲线峰值对应的无因次组合参数CDe2S；根据计算到的CDe2S参数计算并绘制一理论曲线及对应的理论导数曲线，完成曲线拟合。由上述本专利技术提供的技术方案可以看出，基于坐标映射关系，从而利用鼠标位置获得到实测导数曲线的峰值位置，并根据峰值位置计算对应的CDe2S参数；根据获得的CDe2S参数计算并绘制一理论曲线及对应的理论导数曲线；达到了控制光标一次即可实现曲线的快速拟合的目的，大大的提高了工作的效率。与试凑的曲线拟合方法相比，本方法根据均质油藏图版特征，人工获取导数曲线峰值，计算相应参数，达到一次拟合成功的效果。附图说明为了更清楚地说明本专利技术实施例的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本专利技术的一些实施例，对于本领域的普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他附图。图1为本专利技术实施例提供的一种曲线拟合的处理方法的流程图；图2为本专利技术实施例提供的双对数图中绘制辅助线的示意图；图3为本专利技术实施例提供的双对数图中绘制实测曲线及实测导数曲线的示意图；图4为本专利技术实施例提供的曲线拟合过程的示意图；图5为本专利技术实施例提供的完成曲线拟合后的示意图。具体实施方式下面结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本专利技术一部分实施例，而不是全部的实施例。基于本专利技术的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本专利技术的保护范围。图1为本专利技术实施例提供的一种曲线拟合的处理方法的流程图。如图1所示，该方法主要包括如下步骤：步骤11、将实测曲线及对应的实测导数曲线绘制在双对数图上，使得所述实测曲线及实测导数曲线的初始部分，与经过原点且斜率为预定值s的斜线重合；且该实测导数曲线的结尾部分与数值为预定值n的水平线重合。如图2所示，在双对数图（例如，横坐标为时间，其范围为10-2～104，纵坐标为压力，其范围为10-2～102）上绘制两条辅助线：一条斜率为预定值s（例如，1）的斜线，另一条数值为预定值n（例如，0.5）的水平线。所述实测曲线可以是井底压力数据曲线，对应的实测导数曲线为该压力数据求导后的曲线。如图3所示，将两条曲线绘制在双对数图上，并使得所述实测曲线及实测导数曲线的初始部分，与经过原点且斜率为s的斜线重合；且该实测导数曲线的结尾部分与数值为n的水平线重合；即通过辅助线将上述两条曲线的位置进行标定，便于后期的曲线拟合工作。步骤12、接收光标控制指令，根据该指令控制并移动光标，当光标移动至该实测导数曲线的峰值处时，计算该实测导数曲线峰值对应的无因次组合参数CDe2S。由于本实施例的曲线拟合基于无因次组合参数CDe2S，且导数曲线相对于CDe2S参数较为敏感，因此需要计算出实测导数曲线的CDe2S参数；该参数与曲线的峰值存在对应关系，当获知曲线的峰值后，可计算获得对应的CDe2S参数。同理，当得知实测导数曲线峰值对应的CDe2S参数后，则可计算出与实测曲线及对应的实测导数曲线完全重合的理论曲线及对应的理论导数曲线。下面针对实测导数曲线的峰值CDe2S参数的求解过程做详细介绍：示例性的，若实测曲线为均质无限大地层的井底压力数据曲线，则无量纲井底压力表达式为：其中，J0(u)与J1(u)为0阶与1阶实测类Bessel函数（贝塞尔函数）；Y0(u)与Y1(u)为0阶与1阶理论类Bessel函数；为无量纲表达式；为无量纲时间；为无量纲井筒存储常数；S为表皮系数；pwf为井底流动压力（Mpa）；pi为原始地层压力（Mpa）；h为地层有效厚度（m）；μ为流体粘度（Mpa.s）；φ为地层孔隙度；Ct为综合压缩系数（1/Mpa）。求解所述无量纲井底压力表达式获得实测导数峰值对应的时间tD：其中，Δ＝[CDuJ0(u)-(1-CDu2S)J1(u)]2+[CDuY0(u)-(1-CDu2S)Y1(u)]2。对（1）-（2）进行归纳总结后可以得出：Lg(CDe2S)＝2.351187hD+4.73227×10-2；（3）tD/CD＝0.8968505Lg(CDe2S)+3.604613。（4）其中，hD为无量纲化后实测导数曲线的峰值；（4）式表示无量纲化后时间的表达式。由（3）式可知，当得知实测导数曲线的hD后，则可计算出CDe2S参数。本专利技术实施例基于坐标映射关系，利用鼠标位置计算实测导数曲线的hD，进而计算出对应的CDe2S参数。具体方法如下：首先，接收光标控制指令，根据该指令控制并移动光标（即用户点击光标并保持点击状态移动光标时，计算机根据接受到的指令进行相应的操作），当光标移动至该实测导数曲线的峰值处时，利用坐标映射关系计算峰值处的坐标，其公式为：其中，y为光标在屏幕中的纵坐标，y0为双对数图的原点纵坐标对应的屏幕纵坐标，ly为屏幕的高度。再根据峰值坐标计算峰值hD，其公式为：最后结合（3）式计算CDe2S参数。另一方面，为了便于实时观察曲线拟合的结果，在光标移动至实测导数曲线峰值之前，也可以实时计算理论曲线及理论导数曲线，如图4所示。这样，用户就可以根据理论曲线及理论导数曲线是否与实测曲线及实测导数曲线的拟合程度判断，光标的位置是否为实测导数曲线的峰值处。步骤13、根据计算到的CDe2S参数计算并绘制一理论曲线及对应的理论导数曲线，完成曲线拟合。如图5所示，当光标移动至该实测导数曲线的峰值处时，根据计算出的CDe2S参数计算并绘制一理论曲线及对应的理论导数曲线；此时，四条曲线完全重合，完成曲线拟合，可以终止对光标的点击操作，通过上述方式实现了移动光标一次完成曲线拟合的目的，降低了工作的难度及时间，提高了工作的效率。本专利技术实施例基于坐标映射关系，从而利用鼠标位置获得到实测导数曲线的峰值位置，并以此计算获得对应的CDe2S参数；根据获得的CDe2S参数计算并绘制一理论曲线及对应的理论导数曲线；达到了控制光标一次即可实现曲线的快速拟合的目的，大大的提高了工作的效率。通过以上的实施方式的描述，本领域的技术人员本文档来自技高网...

【技术保护点】
一种曲线拟合的处理方法，其特征在于，该方法包括：?将实测曲线及对应的实测导数曲线绘制在双对数图上，使得所述实测曲线及实测导数曲线的初始部分，与经过原点且斜率为预定值s的斜线重合；且该实测导数曲线的结尾部分与数值为预定值n的水平线重合；?接收光标控制指令，根据该指令控制并移动光标，当光标移动至该实测导数曲线的峰值处时，计算该实测导数曲线峰值对应的无因次组合参数CDe2S；?根据计算到的CDe2S参数计算并绘制一理论曲线及对应的理论导数曲线，完成曲线拟合。

【技术特征摘要】
1.一种曲线拟合的处理方法，其特征在于，该方法包括：将实测曲线及对应的实测导数曲线绘制在双对数图上，使得所述实测曲线及实测导数曲线的初始部分，与经过原点且斜率为预定值s的斜线重合；且该实测导数曲线的结尾部分与数值为预定值n的水平线重合；接收光标控制指令，根据该指令控制并移动光标，当光标移动至该实测导数曲线的峰值处时，计算该实测导数曲线峰值对应的无因次组合参数CDe2S；根据计算到的CDe2S参数计算并绘制一理论曲线及对应的理论导数曲线，完成曲线拟合；所述实测曲线包括：均质无限大地层的井底压力数据曲线，其无量纲井底压力表达式为：其中，J0(u)与J1(u)为0阶与1阶第一类Bessel函数；Y0(u)与Y1(u)为0阶与1阶第二类Bessel函数；tD为无量纲时间；CD为无量纲井筒存储常数；S为表皮系数；pwf为井底流动压力；pi为原始地层压力；h为地层有效厚度；μ为流体粘度；φ...

【专利技术属性】
技术研发人员：王磊，卢德唐，李道伦，
申请(专利权)人：中国科学技术大学，
类型：发明
国别省市：