专利名称：一种曲线拟合的处理方法
技术领域：
本发明涉及计算机技术和石油工程领域，尤其涉及一种曲线拟合的处理方法。
背景技术：
现有的试井分析方法主要包括:直线段分析方法、曲线拟合法及非线性回归方法。其中，现有曲线拟合法的本质是通过试凑法实现曲线拟合；S卩，不断的更改参数，每更改一次后，查看拟合结果，直至曲线的外形与弧度近似；再移动曲线使其重合。该方法的缺陷主要为:1)它是试凑方法，对解释工程师水平有很高要求；2)需要多次反复的更改参数，并进行相应计算，工作效率较低；3)当参数范围较大时(相关参数的范围10_6 102°°)更改的次数非常多，无法快速的进行曲线拟合。
发明内容
本发明的目的是提供一种曲线拟合的处理方法，提高了工作效率。本发明的目的是通过以下技术方案实现的:一种曲线拟合的处理方法，该方法包括:将实测曲线及对 应的实测导数曲线绘制在双对数图上，使得所述实测曲线及实测导数曲线的初始部分，与经过原点且斜率为预定值s的斜线重合；且该实测导数曲线的结尾部分与数值为预定值n的水平线重合；接收光标控制指令，根据该指令控制并移动光标，当光标移动至该实测导数曲线的峰值处时，计算该实测导数曲线峰值对应的无因次组合参数CDe2S;根据计算到的CDe2S参数计算并绘制一理论曲线及对应的理论导数曲线，完成曲线拟合。由上述本发明提供的技术方案可以看出，基于坐标映射关系，从而利用鼠标位置获得到实测导数曲线的峰值位置，并根据峰值位置计算对应的CDe2S参数；根据获得的CDe2S参数计算并绘制一理论曲线及对应的理论导数曲线；达到了控制光标一次即可实现曲线的快速拟合的目的，大大的提高了工作的效率。与试凑的曲线拟合方法相比，本方法根据均质油藏图版特征，人工获取导数曲线峰值，计算相应参数，达到一次拟合成功的效果。
为了更清楚地说明本发明实施例的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域的普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他附图。图1为本发明实施例提供的一种曲线拟合的处理方法的流程图；图2为本发明实施例提供的双对数图中绘制辅助线的示意图；图3为本发明实施例提供的双对数图中绘制实测曲线及实测导数曲线的示意图4为本发明实施例提供的曲线拟合过程的示意图；图5为本发明实施例提供的完成曲线拟合后的示意图。
具体实施例方式下面结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明的保护范围。图1为本发明实施例提供的一种曲线拟合的处理方法的流程图。如图1所示，该方法主要包括如下步骤:步骤11、将实测曲线及对应的实测导数曲线绘制在双对数图上，使得所述实测曲线及实测导数曲线的初始部分，与经过原点且斜率为预定值s的斜线重合；且该实测导数曲线的结尾部分与数值为预定值n的水平线重合。如图2所示，在双对数图(例如，横坐标为时间，其范围为10_2 104，纵坐标为压力，其范围为10_2 IO2)上绘制两条辅助线:一条斜率为预定值S (例如，I)的斜线，另一条数值为预定值n (例如，0.5)的水平线。所述实测曲线可以是井底压力数据曲线，对应的实测导数曲线为该压力数据求导后的曲线。如图3所示，将两条曲线绘制在双对数图上，并使得所述实测曲线及实测导数曲线的初始部分，与经过原点且斜率为s的斜线重合；且该实测导数曲线的结尾部分与数值为n的水平线重合；即通过辅助线将上述两条曲线的位置进行标定，便于后期的曲线拟合工作。步骤12、接收光标控制指令，根据该指令控制并移动光标，当光标移动至该实测导数曲线的峰值处时，计算该实测导数曲线峰值对应的无因次组合参数CDe2S。由于本实施例的曲线拟合基于无因次组合参数Qe2s，且导数曲线相对于CDe2S参数较为敏感，因此需要计算出实测导数曲线的CDe2S参数；该参数与曲线的峰值存在对应关系，当获知曲线的峰值后，可计算获得对应的CDe2S参数。同理，当得知实测导数曲线峰值对应的CDe2S参数后，则可计算出与实测曲线及对应的实测导数曲线完全重合的理论曲线及对应的理论导数曲线。下面针对实测导数曲线的峰值CDe2S参数的求解过程做详细介绍:示例性的，若实测曲线为均质无限大地层的井底压力数据曲线，则无量纲井底压力表达式为:
权利要求
1.一种曲线拟合的处理方法，其特征在于，该方法包括: 将实测曲线及对应的实测导数曲线绘制在双对数图上，使得所述实测曲线及实测导数曲线的初始部分，与经过原点且斜率为预定值S的斜线重合；且该实测导数曲线的结尾部分与数值为预定值n的水平线重合； 接收光标控制指令，根据该指令控制并移动光标，当光标移动至该实测导数曲线的峰值处时，计算该实测导数曲线峰值对应的无因次组合参数Qe2s ； 根据计算到的CDe2S参数计算并绘制一理论曲线及对应的理论导数曲线，完成曲线拟口 o
2.根据权利要求1所述的方法，其特征在于，所述实测曲线包括: 均质无限大地层的井底压力数据曲线，其无量纲井底压力表达式为:
3.根据权利要求1或2所述的方法，其特征在于，所述当光标移动至该实测导数曲线的峰值处时，计算该实测导数曲线峰值对应的无因次组合参数CDe2S包括: 当光标移动至该实测导数曲线的峰值处时，根据光标在屏幕中的位置计算所述实测导数曲线的峰值坐标，并根据该实测导数曲线的峰值坐标计算所述实测导数曲线峰值及对应的CDe2S参数。
4.根据权利要求3所述的方法，其特征在于，所述根据光标在屏幕中的位置计算所述实测导数曲线的峰值坐标，并根据该实测导数曲线的峰值坐标计算所述实测导数曲线峰值及对应的CDe2S参数包括: 根据光标在屏幕中的位置计算所述实测导数曲线的峰值纵坐标，其公式为:Jad=(H)Ixio40 其中，y为光标在屏幕中的纵坐标，％为双对数图的原点纵坐标对应的屏幕纵坐标，Iy为屏幕的高度； 根据峰值纵坐标计算峰值hD，其公式为:hD = yact / M ; 由峰值hD计算出对应的CDe2S参数，其公式为:lg(CDe2S) = 2.351187hD+4.73227X 10_2。
全文摘要
本发明公开了一种曲线拟合的处理方法，该方法可基于均质油藏图版特征，根据坐标映射关系利用鼠标位置获得到实测导数曲线的峰值位置，并以此计算获得对应的CDe2S参数；根据获得的CDe2S参数计算并绘制一理论曲线及对应的理论导数曲线；达到了控制光标一次即可实现曲线的快速拟合的目的，大大的提高了工作的效率。
文档编号G06F17/50GK103246778SQ20131017769
公开日2013年8月14日 申请日期2013年5月14日 优先权日2013年5月14日
发明者王磊, 卢德唐, 李道伦 申请人:中国科学技术大学