[发明专利]一种曲线拟合的处理方法

摘要

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

主权项

一种曲线拟合的处理方法，其特征在于，该方法包括：将实测曲线及对应的实测导数曲线绘制在双对数图上，使得所述实测曲线及实测导数曲线的初始部分，与经过原点且斜率为预定值s的斜线重合；且该实测导数曲线的结尾部分与数值为预定值n的水平线重合；接收光标控制指令，根据该指令控制并移动光标，当光标移动至该实测导数曲线的峰值处时，计算该实测导数曲线峰值对应的无因次组合参数CDe2S；根据计算到的CDe2S参数计算并绘制一理论曲线及对应的理论导数曲线，完成曲线拟合；所述实测曲线包括：均质无限大地层的井底压力数据曲线，其无量纲井底压力表达式为：PWD(S,CD,tD)=4π2∫0∞[1-exp(-u2tD)]duu3{[CDuJ0(u)-(1-CDu2S)J1(u)]2+[CDuY0(u)-(1-CDu2S)Y1(u)]2};]]>其中，J0(u)与J1(u)为0阶与1阶第一类Bessel函数；Y0(u)与Y1(u)为0阶与1阶第二类Bessel函数；tD为无量纲时间；CD为无量纲井筒存储常数；S为表皮系数；pwf为井底流动压力；pi为原始地层压力；h为地层有效厚度；μ为流体粘度；φ为地层孔隙度；Ct为综合压缩系数；实测导数峰值对应的时间tD：tD=∫0∞exp(-u2tD)/(uΔ)du∫0∞uexp(-u2tD)/Δdu;]]>其中，Δ＝[CDuJ0(u)‑(1‑CDu2S)J1(u)]2+[CDuY0(u)‑(1‑CDu2S)Y1(u)]2；所述当光标移动至该实测导数曲线的峰值处时，计算该实测导数曲线峰值对应的无因次组合参数CDe2S包括：当光标移动至该实测导数曲线的峰值处时，根据光标在屏幕中的位置计算所述实测导数曲线的峰值坐标，并根据该实测导数曲线的峰值坐标计算所述实测导数曲线峰值及对应的CDe2S参数；所述根据光标在屏幕中的位置计算所述实测导数曲线的峰值坐标，并根据该实测导数曲线的峰值坐标计算所述实测导数曲线峰值及对应的CDe2S参数包括：根据光标在屏幕中的位置计算所述实测导数曲线的峰值纵坐标，其公式为：其中，y为光标在屏幕中的纵坐标，y0为双对数图的原点纵坐标对应的屏幕纵坐标，ly为屏幕的高度；根据峰值纵坐标计算峰值hD，其公式为：由峰值hD计算出对应的CDe2S参数，其公式为：lg(CDe2S)＝2.351187hD+4.73227×10‑2。