專利名稱：一種曲線擬合的處理方法
技術(shù)領(lǐng)域：
本發(fā)明涉及計(jì)算機(jī)技術(shù)和石油工程領(lǐng)域，尤其涉及一種曲線擬合的處理方法。
背景技術(shù)：
現(xiàn)有的試井分析方法主要包括:直線段分析方法、曲線擬合法及非線性回歸方法。其中，現(xiàn)有曲線擬合法的本質(zhì)是通過(guò)試湊法實(shí)現(xiàn)曲線擬合；S卩，不斷的更改參數(shù)，每更改一次后，查看擬合結(jié)果，直至曲線的外形與弧度近似；再移動(dòng)曲線使其重合。該方法的缺陷主要為:1)它是試湊方法，對(duì)解釋工程師水平有很高要求；2)需要多次反復(fù)的更改參數(shù)，并進(jìn)行相應(yīng)計(jì)算，工作效率較低；3)當(dāng)參數(shù)范圍較大時(shí)(相關(guān)參數(shù)的范圍10_6 102°°)更改的次數(shù)非常多，無(wú)法快速的進(jìn)行曲線擬合。

發(fā)明內(nèi)容
本發(fā)明的目的是提供一種曲線擬合的處理方法，提高了工作效率。本發(fā)明的目的是通過(guò)以下技術(shù)方案實(shí)現(xiàn)的:一種曲線擬合的處理方法，該方法包括:將實(shí)測(cè)曲線及對(duì) 應(yīng)的實(shí)測(cè)導(dǎo)數(shù)曲線繪制在雙對(duì)數(shù)圖上，使得所述實(shí)測(cè)曲線及實(shí)測(cè)導(dǎo)數(shù)曲線的初始部分，與經(jīng)過(guò)原點(diǎn)且斜率為預(yù)定值s的斜線重合；且該實(shí)測(cè)導(dǎo)數(shù)曲線的結(jié)尾部分與數(shù)值為預(yù)定值n的水平線重合；接收光標(biāo)控制指令，根據(jù)該指令控制并移動(dòng)光標(biāo)，當(dāng)光標(biāo)移動(dòng)至該實(shí)測(cè)導(dǎo)數(shù)曲線的峰值處時(shí)，計(jì)算該實(shí)測(cè)導(dǎo)數(shù)曲線峰值對(duì)應(yīng)的無(wú)因次組合參數(shù)CDe2S;根據(jù)計(jì)算到的CDe2S參數(shù)計(jì)算并繪制一理論曲線及對(duì)應(yīng)的理論導(dǎo)數(shù)曲線，完成曲線擬合。由上述本發(fā)明提供的技術(shù)方案可以看出，基于坐標(biāo)映射關(guān)系，從而利用鼠標(biāo)位置獲得到實(shí)測(cè)導(dǎo)數(shù)曲線的峰值位置，并根據(jù)峰值位置計(jì)算對(duì)應(yīng)的CDe2S參數(shù)；根據(jù)獲得的CDe2S參數(shù)計(jì)算并繪制一理論曲線及對(duì)應(yīng)的理論導(dǎo)數(shù)曲線；達(dá)到了控制光標(biāo)一次即可實(shí)現(xiàn)曲線的快速擬合的目的，大大的提高了工作的效率。與試湊的曲線擬合方法相比，本方法根據(jù)均質(zhì)油藏圖版特征，人工獲取導(dǎo)數(shù)曲線峰值，計(jì)算相應(yīng)參數(shù)，達(dá)到一次擬合成功的效果。


為了更清楚地說(shuō)明本發(fā)明實(shí)施例的技術(shù)方案，下面將對(duì)實(shí)施例描述中所需要使用的附圖作簡(jiǎn)單地介紹，顯而易見(jiàn)地，下面描述中的附圖僅僅是本發(fā)明的一些實(shí)施例，對(duì)于本領(lǐng)域的普通技術(shù)人員來(lái)講，在不付出創(chuàng)造性勞動(dòng)的前提下，還可以根據(jù)這些附圖獲得其他附圖。圖1為本發(fā)明實(shí)施例提供的一種曲線擬合的處理方法的流程圖；圖2為本發(fā)明實(shí)施例提供的雙對(duì)數(shù)圖中繪制輔助線的示意圖；圖3為本發(fā)明實(shí)施例提供的雙對(duì)數(shù)圖中繪制實(shí)測(cè)曲線及實(shí)測(cè)導(dǎo)數(shù)曲線的示意圖4為本發(fā)明實(shí)施例提供的曲線擬合過(guò)程的示意圖；圖5為本發(fā)明實(shí)施例提供的完成曲線擬合后的示意圖。
具體實(shí)施例方式下面結(jié)合本發(fā)明實(shí)施例中的附圖，對(duì)本發(fā)明實(shí)施例中的技術(shù)方案進(jìn)行清楚、完整地描述，顯然，所描述的實(shí)施例僅僅是本發(fā)明一部分實(shí)施例，而不是全部的實(shí)施例。基于本發(fā)明的實(shí)施例，本領(lǐng)域普通技術(shù)人員在沒(méi)有做出創(chuàng)造性勞動(dòng)前提下所獲得的所有其他實(shí)施例，都屬于本發(fā)明的保護(hù)范圍。圖1為本發(fā)明實(shí)施例提供的一種曲線擬合的處理方法的流程圖。如圖1所示，該方法主要包括如下步驟:步驟11、將實(shí)測(cè)曲線及對(duì)應(yīng)的實(shí)測(cè)導(dǎo)數(shù)曲線繪制在雙對(duì)數(shù)圖上，使得所述實(shí)測(cè)曲線及實(shí)測(cè)導(dǎo)數(shù)曲線的初始部分，與經(jīng)過(guò)原點(diǎn)且斜率為預(yù)定值s的斜線重合；且該實(shí)測(cè)導(dǎo)數(shù)曲線的結(jié)尾部分與數(shù)值為預(yù)定值n的水平線重合。如圖2所示，在雙對(duì)數(shù)圖(例如，橫坐標(biāo)為時(shí)間，其范圍為10_2 104，縱坐標(biāo)為壓力，其范圍為10_2 IO2)上繪制兩條輔助線:一條斜率為預(yù)定值S (例如，I)的斜線，另一條數(shù)值為預(yù)定值n (例如，0.5)的水平線。所述實(shí)測(cè)曲線可以是井底壓力數(shù)據(jù)曲線，對(duì)應(yīng)的實(shí)測(cè)導(dǎo)數(shù)曲線為該壓力數(shù)據(jù)求導(dǎo)后的曲線。如圖3所示，將兩條曲線繪制在雙對(duì)數(shù)圖上，并使得所述實(shí)測(cè)曲線及實(shí)測(cè)導(dǎo)數(shù)曲線的初始部分，與經(jīng)過(guò)原點(diǎn)且斜率為s的斜線重合；且該實(shí)測(cè)導(dǎo)數(shù)曲線的結(jié)尾部分與數(shù)值為n的水平線重合；即通過(guò)輔助線將上述兩條曲線的位置進(jìn)行標(biāo)定，便于后期的曲線擬合工作。步驟12、接收光標(biāo)控制指令，根據(jù)該指令控制并移動(dòng)光標(biāo)，當(dāng)光標(biāo)移動(dòng)至該實(shí)測(cè)導(dǎo)數(shù)曲線的峰值處時(shí)，計(jì)算該實(shí)測(cè)導(dǎo)數(shù)曲線峰值對(duì)應(yīng)的無(wú)因次組合參數(shù)CDe2S。由于本實(shí)施例的曲線擬合基于無(wú)因次組合參數(shù)Qe2s，且導(dǎo)數(shù)曲線相對(duì)于CDe2S參數(shù)較為敏感，因此需要計(jì)算出實(shí)測(cè)導(dǎo)數(shù)曲線的CDe2S參數(shù)；該參數(shù)與曲線的峰值存在對(duì)應(yīng)關(guān)系，當(dāng)獲知曲線的峰值后，可計(jì)算獲得對(duì)應(yīng)的CDe2S參數(shù)。同理，當(dāng)?shù)弥獙?shí)測(cè)導(dǎo)數(shù)曲線峰值對(duì)應(yīng)的CDe2S參數(shù)后，則可計(jì)算出與實(shí)測(cè)曲線及對(duì)應(yīng)的實(shí)測(cè)導(dǎo)數(shù)曲線完全重合的理論曲線及對(duì)應(yīng)的理論導(dǎo)數(shù)曲線。下面針對(duì)實(shí)測(cè)導(dǎo)數(shù)曲線的峰值CDe2S參數(shù)的求解過(guò)程做詳細(xì)介紹:示例性的，若實(shí)測(cè)曲線為均質(zhì)無(wú)限大地層的井底壓力數(shù)據(jù)曲線，則無(wú)量綱井底壓力表達(dá)式為:
權(quán)利要求
1.一種曲線擬合的處理方法，其特征在于，該方法包括: 將實(shí)測(cè)曲線及對(duì)應(yīng)的實(shí)測(cè)導(dǎo)數(shù)曲線繪制在雙對(duì)數(shù)圖上，使得所述實(shí)測(cè)曲線及實(shí)測(cè)導(dǎo)數(shù)曲線的初始部分，與經(jīng)過(guò)原點(diǎn)且斜率為預(yù)定值S的斜線重合；且該實(shí)測(cè)導(dǎo)數(shù)曲線的結(jié)尾部分與數(shù)值為預(yù)定值n的水平線重合； 接收光標(biāo)控制指令，根據(jù)該指令控制并移動(dòng)光標(biāo)，當(dāng)光標(biāo)移動(dòng)至該實(shí)測(cè)導(dǎo)數(shù)曲線的峰值處時(shí)，計(jì)算該實(shí)測(cè)導(dǎo)數(shù)曲線峰值對(duì)應(yīng)的無(wú)因次組合參數(shù)Qe2s ； 根據(jù)計(jì)算到的CDe2S參數(shù)計(jì)算并繪制一理論曲線及對(duì)應(yīng)的理論導(dǎo)數(shù)曲線，完成曲線擬口 o
2.根據(jù)權(quán)利要求1所述的方法，其特征在于，所述實(shí)測(cè)曲線包括: 均質(zhì)無(wú)限大地層的井底壓力數(shù)據(jù)曲線，其無(wú)量綱井底壓力表達(dá)式為:
3.根據(jù)權(quán)利要求1或2所述的方法，其特征在于，所述當(dāng)光標(biāo)移動(dòng)至該實(shí)測(cè)導(dǎo)數(shù)曲線的峰值處時(shí)，計(jì)算該實(shí)測(cè)導(dǎo)數(shù)曲線峰值對(duì)應(yīng)的無(wú)因次組合參數(shù)CDe2S包括: 當(dāng)光標(biāo)移動(dòng)至該實(shí)測(cè)導(dǎo)數(shù)曲線的峰值處時(shí)，根據(jù)光標(biāo)在屏幕中的位置計(jì)算所述實(shí)測(cè)導(dǎo)數(shù)曲線的峰值坐標(biāo)，并根據(jù)該實(shí)測(cè)導(dǎo)數(shù)曲線的峰值坐標(biāo)計(jì)算所述實(shí)測(cè)導(dǎo)數(shù)曲線峰值及對(duì)應(yīng)的CDe2S參數(shù)。
4.根據(jù)權(quán)利要求3所述的方法，其特征在于，所述根據(jù)光標(biāo)在屏幕中的位置計(jì)算所述實(shí)測(cè)導(dǎo)數(shù)曲線的峰值坐標(biāo)，并根據(jù)該實(shí)測(cè)導(dǎo)數(shù)曲線的峰值坐標(biāo)計(jì)算所述實(shí)測(cè)導(dǎo)數(shù)曲線峰值及對(duì)應(yīng)的CDe2S參數(shù)包括: 根據(jù)光標(biāo)在屏幕中的位置計(jì)算所述實(shí)測(cè)導(dǎo)數(shù)曲線的峰值縱坐標(biāo)，其公式為:Jad=(H)Ixio40 其中，y為光標(biāo)在屏幕中的縱坐標(biāo)，％為雙對(duì)數(shù)圖的原點(diǎn)縱坐標(biāo)對(duì)應(yīng)的屏幕縱坐標(biāo)，Iy為屏幕的高度； 根據(jù)峰值縱坐標(biāo)計(jì)算峰值hD，其公式為:hD = yact / M ; 由峰值hD計(jì)算出對(duì)應(yīng)的CDe2S參數(shù)，其公式為:lg(CDe2S) = 2.351187hD+4.73227X 10_2。
全文摘要
本發(fā)明公開(kāi)了一種曲線擬合的處理方法，該方法可基于均質(zhì)油藏圖版特征，根據(jù)坐標(biāo)映射關(guān)系利用鼠標(biāo)位置獲得到實(shí)測(cè)導(dǎo)數(shù)曲線的峰值位置，并以此計(jì)算獲得對(duì)應(yīng)的CDe2S參數(shù)；根據(jù)獲得的CDe2S參數(shù)計(jì)算并繪制一理論曲線及對(duì)應(yīng)的理論導(dǎo)數(shù)曲線；達(dá)到了控制光標(biāo)一次即可實(shí)現(xiàn)曲線的快速擬合的目的，大大的提高了工作的效率。
文檔編號(hào)G06F17/50GK103246778SQ20131017769
公開(kāi)日2013年8月14日 申請(qǐng)日期2013年5月14日 優(yōu)先權(quán)日2013年5月14日
發(fā)明者王磊, 盧德唐, 李道倫 申請(qǐng)人:中國(guó)科學(xué)技術(shù)大學(xué)