威布爾模型與乙型水驅(qū)曲線的聯(lián)解法課程設(shè)計(jì)報(bào)告

1、52/57XXXX 石 油 大 學(xué)課 程 設(shè) 計(jì)課 程 課 程 石油工程課程設(shè)計(jì) 題 目 威布爾模型與乙型水驅(qū)曲線的聯(lián)解法 院 系 石油工程學(xué)院 專業(yè)班級 油工 班 學(xué)生姓名 學(xué)生學(xué)號 指導(dǎo)教師 年 月 日XXXX石油大學(xué)課程設(shè)計(jì)任務(wù)書課程 石油工程課程設(shè)計(jì)題目 威布爾模型與乙型水驅(qū)曲線的聯(lián)解法專業(yè) 石油工程 姓名 學(xué)號 要緊內(nèi)容、差不多要求、要緊參考資料等要緊內(nèi)容：推導(dǎo)威布爾模型及其與乙型水驅(qū)曲線聯(lián)解關(guān)系式的預(yù)測模型；依照油氣田實(shí)際生產(chǎn)數(shù)據(jù)，進(jìn)行線性回歸，求得乙型水驅(qū)曲線的截距、斜率；依照相關(guān)公式，確定威布爾（Weibull）模型常數(shù)a、b、c；計(jì)算油田年產(chǎn)油量；計(jì)算油田累積產(chǎn)油量；計(jì)算含水

2、率；計(jì)算可采儲量；計(jì)算最高年產(chǎn)量發(fā)生的時刻；計(jì)算最高年產(chǎn)量；(10)繪制油田實(shí)際年產(chǎn)量與預(yù)測產(chǎn)量對比曲線；(11)繪制油田實(shí)際累積產(chǎn)量與預(yù)測累積產(chǎn)量對比曲線。差不多要求：該專題設(shè)計(jì)最終要求是，學(xué)生通過自選基礎(chǔ)數(shù)據(jù)，利用威布爾模型與乙型水驅(qū)曲線聯(lián)解預(yù)測模型，結(jié)合油田實(shí)際生產(chǎn)資料，預(yù)測選定油氣田的年產(chǎn)量、累積產(chǎn)量、含水率、可采儲量、最高年產(chǎn)量、最高年產(chǎn)量發(fā)生的時刻等開發(fā)指標(biāo)，并編制相應(yīng)軟件，并提交規(guī)范設(shè)計(jì)報(bào)告。要緊參考資料1陳元千.預(yù)測油氣田產(chǎn)量的Weibull模型J.新疆石油地質(zhì)，1995，16（3）：250-255.2童憲章.天然水驅(qū)和人工注水油藏的統(tǒng)計(jì)規(guī)律探討J.石油勘探與開發(fā)，1978，4

3、（6）：38-64.3 陳元千.一種新型水驅(qū)曲線關(guān)系式的推導(dǎo)及應(yīng)用J.石油學(xué)報(bào)，1993，14（2）：65-73.4 陳元千.水驅(qū)曲線關(guān)系式的推導(dǎo)J. 石油學(xué)報(bào)，1985，6（2）：69-78.完成期限 指導(dǎo)教師 專業(yè)負(fù)責(zé)人 年 月 日目 錄 TOC o 1-2 h z HYPERLINK l _Toc268286572 1 前 言 PAGEREF _Toc268286572 h 1 HYPERLINK l _Toc268286573 1.1 設(shè)計(jì)的目的意義 PAGEREF _Toc268286573 h 1 HYPERLINK l _Toc268286574 1.2 設(shè)計(jì)的要緊內(nèi)容 PAGE

4、REF _Toc268286574 h 1 HYPERLINK l _Toc268286575 2 基礎(chǔ)數(shù)據(jù) PAGEREF _Toc268286575 h 2 HYPERLINK l _Toc268286576 3 基礎(chǔ)理論 PAGEREF _Toc268286576 h 3 HYPERLINK l _Toc268286577 3.1威布爾預(yù)測模型的建立 PAGEREF _Toc268286577 h 3 HYPERLINK l _Toc268286578 3.2威布爾模型與乙型水驅(qū)曲線的聯(lián)解法 PAGEREF _Toc268286578 h 4 HYPERLINK l _Toc268286

5、579 3.3模型的求解方法 PAGEREF _Toc268286579 h 5 HYPERLINK l _Toc268286580 4 設(shè)計(jì)結(jié)果 PAGEREF _Toc268286580 h 7 HYPERLINK l _Toc268286581 認(rèn)識與結(jié)論 PAGEREF _Toc268286581 h 8 HYPERLINK l _Toc268286582 參考文獻(xiàn) PAGEREF _Toc268286582 h 9 HYPERLINK l _Toc268286583 附錄：計(jì)算機(jī)程序與結(jié)果 PAGEREF _Toc268286583 h 101 前 言1.1 設(shè)計(jì)的目的意義本課題的目

6、的是讓學(xué)生通過自選一組數(shù)據(jù)，利用所學(xué)專業(yè)知識在指導(dǎo)教師的指導(dǎo)下獨(dú)立完成對某一油田或區(qū)塊開發(fā)指標(biāo)的預(yù)測。本課題要求學(xué)生對威布爾（Weibull）模型及其與乙型水驅(qū)曲線聯(lián)解關(guān)系式進(jìn)行推導(dǎo)，其結(jié)果包括，對油氣田的年產(chǎn)量、累積產(chǎn)量、含水率、可采儲量、最高年產(chǎn)量、最高年產(chǎn)量發(fā)生的時刻的預(yù)測方法。從而將理論知識和實(shí)際問題相結(jié)合，通過該專題設(shè)計(jì)的訓(xùn)練，加強(qiáng)學(xué)生理論知識運(yùn)用能力，計(jì)算機(jī)技術(shù)應(yīng)用能力及解決實(shí)際問題的工程應(yīng)用能力。1.2 設(shè)計(jì)的要緊內(nèi)容依照已有的基礎(chǔ)數(shù)據(jù)，利用所學(xué)的專業(yè)知識，在指導(dǎo)教師指導(dǎo)下獨(dú)立完成并提交一個油田或一個區(qū)塊油田開發(fā)指標(biāo)預(yù)測結(jié)果，設(shè)計(jì)要緊內(nèi)容如下：（1）推導(dǎo)威布爾模型及其與乙型水驅(qū)曲

7、線聯(lián)解關(guān)系式；（2）依照油氣田實(shí)際生產(chǎn)數(shù)據(jù)，進(jìn)行線性回歸，求得乙型水驅(qū)曲線的截距、斜率；（3）依照相關(guān)公式，確定威布爾（Weibull）模型常數(shù)a、b、c；（4）計(jì)算油田年產(chǎn)油量；（5）計(jì)算油田累積產(chǎn)油量；（6）計(jì)算含水率；（7）計(jì)算可采儲量；（8）計(jì)算最高年產(chǎn)量發(fā)生的時刻；（9）計(jì)算最高年產(chǎn)量；（10）繪制油田實(shí)際年產(chǎn)量與預(yù)測產(chǎn)量對比曲線；（11）繪制油田實(shí)際累積產(chǎn)量與預(yù)測累積產(chǎn)量對比曲線。2 基礎(chǔ)數(shù)據(jù)某油田的開發(fā)數(shù)據(jù)如下：表2-1 某油田開發(fā)數(shù)據(jù)年份時刻（a）Qo（104t/a）Qw（104t/a）Np（104t）Wp（104t）Lp（104t）1968125.73 0.75 25.73

8、0.75 26.47 1969237.80 1.33 63.53 2.08 65.60 1970350.84 1.63 114.37 3.71 118.07 1971472.05 1.85 186.42 5.56 191.97 1972585.44 4.25 271.86 9.80 281.66 1973696.68 7.24 368.54 17.04 385.58 19747115.16 14.67 483.70 31.71 515.42 19758129.69 24.08 613.39 55.79 669.18 19769127.82 26.31 741.21 82.10 823.31 1

9、97710153.83 39.71 895.03 121.81 1016.84 197811181.82 63.01 1076.85 184.82 1261.67 197912189.15 102.85 1266.00 287.67 1553.67 198013179.77 146.89 1445.77 434.57 1880.33 198114171.35 193.05 1617.12 627.62 2244.73 198215169.18 230.47 1786.29 858.09 2644.38 198316159.08 301.06 1945.37 1159.15 3104.51 19

10、8417152.74 381.82 2098.11 1540.97 3639.08 198518136.70 431.01 2234.81 1971.98 4206.79 198619122.07 483.31 2356.88 2455.29 4812.17 198720104.16 480.49 2461.04 2935.78 5396.82 19882194.44 518.30 2555.48 3454.08 6009.56 19892283.86 517.19 2639.34 3971.26 6610.61 19902375.60 540.77 2714.94 4512.04 7226.

11、98 3 基礎(chǔ)理論威布爾（Weibull）預(yù)測模型能夠?qū)ｉT好地預(yù)測油田產(chǎn)油量隨時刻的變化關(guān)系，但卻不能預(yù)測油田的含水率、產(chǎn)水量、產(chǎn)液量及累積產(chǎn)水量和累積產(chǎn)液量，而這些開發(fā)指標(biāo)正是水驅(qū)開發(fā)油田所需要預(yù)測的。乙型水驅(qū)曲線法是油藏工程中重要的預(yù)測方法，但它只能預(yù)測累積產(chǎn)水量與累積產(chǎn)油量之間的關(guān)系，卻不能預(yù)測開發(fā)指標(biāo)與開發(fā)時刻的關(guān)系，而油田開發(fā)指標(biāo)的預(yù)測，都離不開與開發(fā)時刻的聯(lián)系。將Weibull預(yù)測模型和乙型水驅(qū)曲線法相結(jié)合，則能夠預(yù)測水驅(qū)油田的含水率、產(chǎn)油量、產(chǎn)水量、產(chǎn)液量、可采儲量及其相應(yīng)的累積產(chǎn)量隨開發(fā)時刻變化的聯(lián)解法。聯(lián)解法既能保持兩種方法原有的預(yù)測功能，又克服了兩者的局限性。3.1威布爾預(yù)

12、測模型的建立Weibull(威布爾)于1939年提出的統(tǒng)計(jì)分布模型已成為生命試驗(yàn)和可靠性理論研究的基礎(chǔ)。該模型的分布密度表示為： (3-1)式中 f(x)威布爾分布的分布密度函數(shù)；x分布變量，依照實(shí)際問題，分布區(qū)間為0；操縱分布形態(tài)的形狀參數(shù)；操縱分布峰位和峰值的尺度參數(shù)。若對(3-l)式進(jìn)行積分，在x為0區(qū)間內(nèi)，能夠得到Weibu11的分布函數(shù)值等于1，推證如下： （3-2）為將Weibull分布模型用于油氣田開發(fā)指標(biāo)的預(yù)測，將(3-l)式改寫為 (3-3)式中 Q油氣田的年產(chǎn)量，104t/a（油）或108m3/a（氣）； t 油氣田的開發(fā)時刻，a； C由Weibull分布模型轉(zhuǎn)換為油氣田開

13、發(fā)有用模型的模型轉(zhuǎn)換常數(shù)。 油氣田的累積產(chǎn)量表達(dá)式為： （3-4）式中 NP油氣田的累積產(chǎn)量，104t或108t(油)；108m(氣)。式中 NP油氣田的累積產(chǎn)量，104t或108t(油)；108 將（3-3）代入（3-4）式并考慮（3-2）式中的變量變換法，t從0到t積分得： （3-5） 當(dāng)時，則,因此（3-5）式又可改寫為： （3-6）在得到上面的結(jié)果之后，便可對模型轉(zhuǎn)換常數(shù)的性質(zhì)和作用做如此的講明：由于Weibull分布模型，在x從0到區(qū)間的分布函數(shù)F（x）=1.0,這相當(dāng)于實(shí)際開發(fā)的油氣田，在t從0到區(qū)間的累積產(chǎn)量，即油氣田的可采儲量。因此，為了能夠得到（3-5）式的結(jié)果，就必須在（3

14、-3）中引入模型轉(zhuǎn)換常數(shù)C。而該模型轉(zhuǎn)換常數(shù)確實(shí)是油氣田的可采儲量。因此，能夠?qū)ⅲ?-3）式再改寫為： （3-7）為了確定最高年產(chǎn)量發(fā)生的時刻，有（3-7）式對時刻t求導(dǎo)數(shù)得： （3-8）當(dāng)=0時，必定有=0，故能夠得到最高年產(chǎn)量發(fā)生的時刻tm為： （3-9）將（3-9）式代入（3-7）式，得到油氣田的最高年產(chǎn)量（Qmax）的表達(dá)式： （3-10）再將（3-9）式代入（3-6）式，得到油氣田最高年產(chǎn)量發(fā)生時的累積產(chǎn)量（Npm）為： （3-11）油氣田的剩余可采儲量NRR表示為： （3-12）將（3-6）式代入（3-12）式得： （3-13）剩余可采儲量的儲采比表示為： （3-14）將（3-7）

15、式和（3-13）式代入（3-14）式得： （3-15）剩余可采儲量的采油速度為儲采比的倒數(shù)，故由（3-15）式得到剩余可采儲量采油速度的表達(dá)式： （3-16）式中以小數(shù)f表示，若改以百分?jǐn)?shù)%表示時，（3-16）式改為下式： （3-17）3.2威布爾模型與乙型水驅(qū)曲線的聯(lián)解法利用數(shù)理統(tǒng)計(jì)學(xué)中的威布爾（Weibull）分布，研究與推導(dǎo)得到了威布爾（Weibull）預(yù)測模型1。該模型具有預(yù)測油田產(chǎn)量、累積產(chǎn)量和可采儲量的功能，其差不多關(guān)系式分不為： （3-18） （3-19） （3-20）乙型水驅(qū)曲線法，首先是由我國聞名專家童憲章先生2以經(jīng)驗(yàn)公式的形式，于1978年提出。它的理論推導(dǎo)由文獻(xiàn)3完成，其

16、差不多關(guān)系式為： （3-21）由（3-21）式對時刻t求導(dǎo)數(shù)得： （3-22）已知：；故由（3-22）式得： （3-23）將（3-23）式代入（3-21）式得： （3-24）取經(jīng)濟(jì)極限水油比(Rwo) L，由（24）式得到預(yù)測油田可采儲量的關(guān)系式： （3-25）已知水油比與含水率的關(guān)系為： （3-26）將（3-26）式代入（3-24）式得： （3-27）將（3-19）式代入（3-27）式得： （3-28）當(dāng)由（3-18）式和（3-28）式得到預(yù)測的產(chǎn)油量和含水率之后，可由下面的公式分不預(yù)測油田的產(chǎn)水量和產(chǎn)液量： （3-29） （3-30）最高年產(chǎn)量發(fā)生的時刻tm的計(jì)算公式如下：最高年產(chǎn)量Qma

17、x為：3.3模型的求解方法為了確定預(yù)測模型的模型常數(shù)a、b、c以及可采儲量NR的數(shù)值，對（3-18）式可進(jìn)行如下處理： （3-31）若設(shè)： （3-32） （3-33）則得： （3-34）依照實(shí)際的開發(fā)數(shù)據(jù)，首先利用（3-34）式進(jìn)行線性試差求解，依照最大線性相關(guān)系數(shù)求出b，然后利用最小二乘法求得和。再由（3-32）式和（3-33）式改寫的下式，分不確定模型的常數(shù)a和c的數(shù)值： （3-35） （3-36）確定出預(yù)測模型參數(shù)a、b、c后，即可依照（3-20）式求解出可采儲量NR。在確定預(yù)測模型常數(shù)a、b、c、NR時，其值是否正確可靠，要利用（3-18）式、（3-19）式、（3-28）式預(yù)測的理論產(chǎn)

18、油量、累積產(chǎn)油量和含水率，與實(shí)際產(chǎn)油量、累積產(chǎn)油量和含水率進(jìn)行對比加以確定，而達(dá)到最佳擬合效果的參數(shù)才是最準(zhǔn)確、可靠的。4 設(shè)計(jì)結(jié)果將表1中的累積產(chǎn)液量（Lp）和相應(yīng)的累積產(chǎn)油量（NP）數(shù)據(jù)，按照（3-21）式的直線關(guān)系繪于圖4-1，得到了一條專門好的直線。經(jīng)線性回歸求得直線的截距A=2.50；直線的斜率B=0.000502；直線的相關(guān)系數(shù)r=0.9999。圖4-1 該油田的乙型水驅(qū)曲線將表1中的Log(Qo/t)和相應(yīng)的t(b+1)數(shù)據(jù)，按照（3-34）式的直線關(guān)系繪于圖4-2，得到了一條專門好的直線。經(jīng)線性回歸求得直線的截距=1.32；直線的斜率=-0.00132；參數(shù)a=21.002，b

19、=1.1，c=328.039，直線的相關(guān)系數(shù)r=0.9990。圖4-2 該油田的weibull曲線在確定預(yù)測模型常數(shù)a、b、c后，利用（3-18）式、（3-19）式、（3-28）式計(jì)算油田年產(chǎn)油量、計(jì)算油田累積產(chǎn)油量、計(jì)算含水率如下表：表4-1 數(shù)據(jù)表時刻（a）實(shí)際產(chǎn)量預(yù)測產(chǎn)量實(shí)際累計(jì)產(chǎn)量預(yù)測累計(jì)產(chǎn)量實(shí)際含水率預(yù)測含水率125.7320.9425.739.952.830237.8044.4463.5342.443.400350.8468.21114.3798.583.100472.0591.26186.42178.112.500585.44112.82271.86279.934.730696.

20、68132.26368.54402.216.9607115.16149.05483.70542.5911.3008129.69162.82613.39698.2115.6609127.82173.31741.21865.9417.07010153.83180.38895.031042.4220.5217.1811181.82184.051076.851224.2525.7432.8912189.15184.431266.001408.0735.2245.7413179.77181.761445.771590.7344.9756.0714171.35176.351617.121769.3452.

21、9864.2715169.18168.61786.291941.3657.6770.7216159.08158.941945.372104.6765.4375.7617152.74147.812098.112257.5771.4379.6918136.70135.662234.812398.8575.9282.7519122.07122.942356.882527.779.8485.1420104.16110.022461.042643.7682.18872194.4497.272555.48274784.5988.472283.8684.972639.342837.7386.0589.612

22、375.6073.352714.942916.5487.7390.522462.592984.1891.232552.793041.5791.82644.023089.7192.242736.293129.6392.592829.583162.3692.872923.853188.993.083019.013210.1893.253114.993227.0693.383211.693240.2993.48339.013250.5693.56346.883258.4393.62355.193264.4193.66363.873268.993.69372.863272.2393.72382.093

23、274.6793.73391.513276.4593.75401.083277.7393.76410.763278.6493.76420.533279.2893.77430.373279.7293.77440.253280.0393.77450.173280.2393.77460.113280.3793.78470.083280.4793.78480.053280.5393.78490.033280.5793.78500.023280.693.78510.013280.6193.78計(jì)算可采儲量：NR =3280.639（104t），計(jì)算最高年產(chǎn)量發(fā)生的時刻：11.897a，計(jì)算最高年產(chǎn)量：Q

24、max=184.302(104t)。預(yù)測的理論產(chǎn)油量、累積產(chǎn)油量和含水率，與實(shí)際產(chǎn)油量、累積產(chǎn)油量和含水率進(jìn)行對比。按照（3-18）式的關(guān)系在同一坐標(biāo)系中繪制實(shí)際產(chǎn)量與預(yù)測產(chǎn)量及其對應(yīng)時刻t之間的關(guān)系曲線，繪于圖4-3。圖4-3 產(chǎn)量與時刻關(guān)系曲線按照（3-19）式的關(guān)系在同一坐標(biāo)系中繪制實(shí)際累計(jì)產(chǎn)量與預(yù)測累計(jì)產(chǎn)量及其對應(yīng)時刻t之間的關(guān)系曲線，繪于圖4-4。 圖4-4 累計(jì)產(chǎn)量與時刻關(guān)系曲線按照（3-28）式的關(guān)系在同一坐標(biāo)系中繪制實(shí)際含水率與預(yù)測含水率及其對應(yīng)時刻t之間的關(guān)系曲線，繪于圖4-5。 圖4-5 含水率與時刻關(guān)系曲線認(rèn)識與結(jié)論通過本次課程設(shè)計(jì)，得出如下認(rèn)識與結(jié)論：（1）通過此次油田

25、開發(fā)的課程設(shè)計(jì)，我差不多掌握了威布爾（Weibull）模型與乙型水驅(qū)曲線聯(lián)解法，并能初步對油田進(jìn)行預(yù)測分析。（2）通過這次的課程設(shè)計(jì)，培養(yǎng)了我獨(dú)立考慮和動手操作的能力，也使我在日常的學(xué)習(xí)與討論當(dāng)中增長了許多課外知識，認(rèn)識到數(shù)學(xué)模型對油氣田的產(chǎn)量預(yù)測的重要性。盡管在課程設(shè)計(jì)的過程中遇到了一些問題，但通過反復(fù)的考慮與檢查找出緣故所在并進(jìn)行改正。（3）通過計(jì)算機(jī)的編程學(xué)習(xí)與操作，使我對VB語言有了更進(jìn)一步的認(rèn)識和了解。同時，提高了使用計(jì)算機(jī)技應(yīng)用的能力以及解決實(shí)際問題的工程應(yīng)用能力，使理論與實(shí)踐相結(jié)合。（4）現(xiàn)在能熟練掌握、運(yùn)用威布爾（Weibull）模型與乙型水驅(qū)曲線聯(lián)解法的原理、推導(dǎo)及公式，能夠

26、對油氣田的年產(chǎn)量、累計(jì)產(chǎn)量、含水率、可采儲量、最高年產(chǎn)量、最高年產(chǎn)量發(fā)生的時刻以及與最高年產(chǎn)量對應(yīng)的累積產(chǎn)量進(jìn)行預(yù)測。（5）這次課程設(shè)計(jì)，使我認(rèn)識到學(xué)好計(jì)算機(jī)的重要性與實(shí)踐操作的必要性，在老師的指導(dǎo)下和與同學(xué)的討論中，不斷地發(fā)覺問題、不斷地解決問題、不斷地領(lǐng)悟、不斷地獵取，使自己的程序不斷完善。（6）課程設(shè)計(jì)誠然是一門專業(yè)課，它給了我專門多的專業(yè)知識以及專業(yè)上的提升。課程設(shè)計(jì)同時又是一門講道課，一門思辨課，給了我許多道，給了我許多思，得到專門大的提升空間。參考文獻(xiàn)1陳元千.預(yù)測油氣田產(chǎn)量的Weibull模型J.新疆石油地質(zhì)，1995，16（3）：250-255.2童憲章.天然水驅(qū)和人工注水油藏

27、的統(tǒng)計(jì)規(guī)律探討J.石油勘探與開發(fā)，1978，4（6）：38-64.3 陳元千.一種新型水驅(qū)曲線關(guān)系式的推導(dǎo)及應(yīng)用J.石油學(xué)報(bào)，1993，14（2）：65-73.4 陳元千.水驅(qū)曲線關(guān)系式的推導(dǎo)J. 石油學(xué)報(bào)，1985，6（2）：69-78.附錄：計(jì)算機(jī)程序與結(jié)果計(jì)算機(jī)程序Option Base 0Const e = 2.708281828Dim a, t(200), qt(200), t0(200), q0(200), qw(200), ye(200), wp(200)Dim b, c, dDim afw, qmax, tm, NR, x, y, a1, a2, r, rmax, bmax,

28、aa1, bb1, cc1, ub, aa, ab, ar, br, yt, an, bn, mnDim s(), px(), py()Dim sx, sy, sxy, sx2, sy2Dim np(200), lp(200)Dim yq0(200), ynp(200), yfw(200), ynr, fw(200)Private Sub c1_Click()yt = 50Rem 源數(shù)據(jù)輸入Open App.Path + cd.txt For Input As #1ub = -1Do Until EOF(1)ub = ub + 1Input #1, ye(ub), t(ub), q0(ub),

29、 qw(ub), np(ub), wp(ub), lp(ub)lp(ub) = Log(lp(ub) / Log(10)LoopClose #1Private Sub c5_Click()Rem 乙型參數(shù)計(jì)算For an = ub To 2 Step -1 求相關(guān)系數(shù)r,趨勢線sx = 0: sy = 0: sxy = 0: sx2 = 0: sy2 = 0l = ub - anFor i = l To an + lsx = sx + np(i)Next iFor i = l To an + lsy = sy + lp(i)Next iFor i = l To an + lsxy = sxy

30、+ np(i) * lp(i)Next iFor i = l To an + lsx2 = sx2 + np(i) 2Next iFor i = l To an + lsy2 = sy2 + lp(i) 2Next ia2 = (an + 1) * sxy - sx * sy) / (an + 1) * sx2 - sx 2)a1 = (sx2 * sy - sx * sxy) / (an + 1) * sx2 - sx 2)r = (an + 1) * sxy - sx * sy) / (an + 1) * sx2 - sx 2) 0.5 * (an + 1) * sy2 - sy 2) 0

31、.5)取最大值rIf r ar Then ar = rNext anaa = a1ab = a2Form3.Text1.Text = aa 截距，參數(shù)AForm3.Text2.Text = ab 斜率，參數(shù)BForm3.Text3.Text = ar 最大相關(guān)系數(shù)Form3.ShowRem weibull參數(shù)計(jì)算bn = 12l = ub - bnbmax = 0rmax = 0For b = 1 To 2 Step 0.05 sx = 0: sy = 0: sxy = 0: sx2 = 0: sy2 = 0 For i = l To l + bn t(i) = i + 1 Next i Fo

32、r i = l To l + bn qt(i) = Log(q0(i) / (t(i) b) / Log(10) t0(i) = t(i) (b + 1) Next i For i = l To l + bn sx = sx + t0(i) Next i For i = l To l + bn sy = sy + qt(i) Next i For i = l To l + bn sxy = sxy + qt(i) * t0(i) Next i For i = l To l + bn sx2 = sx2 + t0(i) 2 Next i For i = l To l + bn sy2 = sy2

33、 + qt(i) 2 Next i r = Abs(bn + 1) * sxy - sx * sy) / (bn + 1) * sx2 - sx 2) 0.5 * (bn + 1) * sy2 - sy 2) 0.5)最大相關(guān)系數(shù) If r rmax Then r取最大值 rmax = r bmax = b a2 = (bn + 1) * sxy - sx * sy) / (bn + 1) * sx2 - sx 2) a1 = (sx2 * sy - sx * sxy) / (bn + 1) * sx2 - sx 2) End IfNext baa1 = 10 a1bb1 = bmaxcc1

34、= -1 / (2.303 * a2)afw = 100 * (1 - 10 (-(aa + ab * aa1 * cc1 / (bb1 + 1) + Log(2.303 * ab) / Log(10) 極限含水率fw公式，fw=1-10(-A+B(a*c/(b+1)(1-exp(-t(b+1)/c)tm = (bb1 * cc1 / (bb1 + 1) (1 / (bb1 + 1)最高年產(chǎn)量對應(yīng)時刻tm公式，tm=(b*c/(b+1)(1/(b+1)qmax = aa1 * (bb1 * cc1 / (bb1 + 1) (bb1 / (bb1 + 1) * Exp(-bb1 / (bb1 +

35、 1)最高年產(chǎn)量Qmax公式，Qmax=a*(b*c/(b+1)(b/(b+1)*exp(-b/(b+1)NR = aa1 * cc1 / (bb1 + 1)可采儲量NR公式，NR=a*c/(b+1)Form3.Text4.Text = a2 斜率Form3.Text5.Text = a1 截距Form3.Text6.Text = rmax 最大相關(guān)系數(shù)Form3.Text7.Text = aa1 參數(shù)aForm3.Text8.Text = bb1 參數(shù)bForm3.Text9.Text = cc1 參數(shù)cForm3.Text10.Text = afw 極限含水率Form3.Text11.Te

36、xt = qmax 最高年產(chǎn)量Form3.Text12.Text = tm 最高年產(chǎn)量對應(yīng)時刻Form3.Text13.Text = NR 可采儲量Rem 預(yù)測數(shù)據(jù)For i = 0 To yt t(i) = i + 1 yq0(i) = Round(aa1 * (t(i) bb1) * (e (-t(i) (bb1 + 1) / cc1), 2) ynp(i) = Round(aa1 * cc1 * (1 - (e (-t(i) (bb1 + 1) / cc1) / (bb1 + 1), 2) If i = ub Then fw(i) = Round(100 * qw(i) / (qw(i)

37、 + q0(i), 2) yfw(i) = Round(100 * (1 - 10 (-(aa + ab * ynp(i) + Log(2.303 * ab) / Log(10), 2) If yfw(i) 0 Then yfw(i) = 0 Next i mn = 50 Rem 表格輸出MSFlexGrid1.ColWidth(0) = 1000 表格行寬MSFlexGrid1.ColWidth(1) = 1200MSFlexGrid1.ColWidth(2) = 1200MSFlexGrid1.ColWidth(3) = 1200MSFlexGrid1.ColWidth(4) = 1200

38、MSFlexGrid1.ColWidth(5) = 1200MSFlexGrid1.ColWidth(6) = 1200MSFlexGrid1.Rows = 201MSFlexGrid1.Cols = 7MSFlexGrid1.Width = MSFlexGrid1.ColWidth(1) * MSFlexGrid1.Cols + MSFlexGrid1.Cols * 12MSFlexGrid1.Height = MSFlexGrid1.RowHeight(1) * MSFlexGrid1.Rows + MSFlexGrid1.Rows * 12MSFlexGrid1.TextMatrix(0

39、, 0) = 時刻/aMSFlexGrid1.TextMatrix(0, 1) = 實(shí)際產(chǎn)量MSFlexGrid1.TextMatrix(0, 2) = 預(yù)測產(chǎn)量MSFlexGrid1.TextMatrix(0, 3) = 實(shí)際累計(jì)產(chǎn)量MSFlexGrid1.TextMatrix(0, 4) = 預(yù)測累計(jì)產(chǎn)量MSFlexGrid1.TextMatrix(0, 5) = 實(shí)際含水率%MSFlexGrid1.TextMatrix(0, 6) = 預(yù)測含水率%For i = 0 To mnMSFlexGrid1.TextMatrix(i + 1, 0) = Format(t(i), 0.#)MSFl

40、exGrid1.TextMatrix(i + 1, 1) = Format(q0(i), 0.#)MSFlexGrid1.TextMatrix(i + 1, 2) = Format(yq0(i), 0.#)MSFlexGrid1.TextMatrix(i + 1, 3) = Format(np(i), 0.#)MSFlexGrid1.TextMatrix(i + 1, 4) = Format(ynp(i), 0.#)MSFlexGrid1.TextMatrix(i + 1, 5) = Format(fw(i), 0.#)MSFlexGrid1.TextMatrix(i + 1, 6) = Fo

41、rmat(yfw(i), 0.#)Next i橫縱坐標(biāo)With MSChart1.ColumnCount = 2 .Plot.SeriesCollection(1).Pen.Width = 15 .Plot.SeriesCollection(1).Pen.Style = VtPenStyleSolid .Plot.Axis(VtChAxisIdX).ValueScale.Auto = True .Plot.Axis(VtChAxisIdY).ValueScale.Auto = True .Plot.Axis(VtChAxisIdY).ValueScale.Minimum = 0 .Plot.A

42、xis(VtChAxisIdX).ValueScale.Minimum = 0 .Plot.Axis(VtChAxisIdX).ValueScale.MajorDivision = 6 .Plot.Axis(VtChAxisIdY).ValueScale.MajorDivision = 6 .Plot.Axis(VtChAxisIdX).ValueScale.MinorDivision = 0 .Plot.Axis(VtChAxisIdY).ValueScale.MinorDivision = 0 .Plot.Axis(VtChAxisIdX).AxisGrid.MajorPen.Style

43、= VtPenStyleDotted .Plot.Axis(VtChAxisIdY).AxisGrid.MajorPen.Style = VtPenStyleDotted MSChart1.Plot.AutoLayout = False MSChart1.Plot.UniformAxis = False MSChart1.chartType = VtChChartType2dXY 二維坐標(biāo) .Column = 1End WithWith MSChart2.ColumnCount = 2 .Plot.SeriesCollection(1).Pen.Width = 15 .Plot.SeriesC

44、ollection(1).Pen.Style = VtPenStyleSolid .Plot.Axis(VtChAxisIdX).ValueScale.Auto = True .Plot.Axis(VtChAxisIdY).ValueScale.Auto = True .Plot.Axis(VtChAxisIdY).ValueScale.Minimum = 0 .Plot.Axis(VtChAxisIdX).ValueScale.Minimum = 0 .Plot.Axis(VtChAxisIdX).ValueScale.MajorDivision = 6 .Plot.Axis(VtChAxi

45、sIdY).ValueScale.MajorDivision = 6 .Plot.Axis(VtChAxisIdX).ValueScale.MinorDivision = 0 .Plot.Axis(VtChAxisIdY).ValueScale.MinorDivision = 0 .Plot.Axis(VtChAxisIdX).AxisGrid.MajorPen.Style = VtPenStyleDotted .Plot.Axis(VtChAxisIdY).AxisGrid.MajorPen.Style = VtPenStyleDotted MSChart1.Plot.AutoLayout

46、= False MSChart1.Plot.UniformAxis = False MSChart1.chartType = VtChChartType2dXY .Column = 1 End With With MSChart3.ColumnCount = 2 .Plot.SeriesCollection(1).Pen.Width = 15 .Plot.SeriesCollection(1).Pen.Style = VtPenStyleSolid .Plot.Axis(VtChAxisIdX).ValueScale.Auto = True .Plot.Axis(VtChAxisIdY).Va

47、lueScale.Auto = True .Plot.Axis(VtChAxisIdY).ValueScale.Minimum = 0 .Plot.Axis(VtChAxisIdX).ValueScale.Minimum = 0 .Plot.Axis(VtChAxisIdX).ValueScale.MajorDivision = 6 .Plot.Axis(VtChAxisIdY).ValueScale.MajorDivision = 6 .Plot.Axis(VtChAxisIdX).ValueScale.MinorDivision = 0 .Plot.Axis(VtChAxisIdY).Va

48、lueScale.MinorDivision = 0 .Plot.Axis(VtChAxisIdX).AxisGrid.MajorPen.Style = VtPenStyleDotted .Plot.Axis(VtChAxisIdY).AxisGrid.MajorPen.Style = VtPenStyleDotted MSChart1.Plot.AutoLayout = False MSChart1.Plot.UniformAxis = False MSChart1.chartType = VtChChartType2dXY .Column = 1End WithEnd SubPrivate

49、 Sub c2_Click()MSFlexGrid1.Visible = True 輸出數(shù)據(jù)在表格內(nèi)Rem 曲線繪制Rem 數(shù)據(jù)輸出End SubPrivate Sub c3_Click()文件輸出Open App.Path + 3.txt For Output As #2Print #2, 年份, 時刻, 預(yù)測產(chǎn)量, 預(yù)測累計(jì)產(chǎn)量, 預(yù)測含水率For i = 0 To ytPrint #2, ye(0) + i, i + 1, yq0(i), ynp(i), yfw(i)Next iClose #2PrintEnd SubPrivate Sub c4_Click()EndEnd SubPri

50、vate Sub d1_Click()理論產(chǎn)油量Qo與時刻t的關(guān)系曲線繪制MSChart1.Visible = TrueReDim px(yt), py(yt), s(yt, 3)st = yt / 500For i = 1 To mnpx(i) = ipy(i) = aa1 * (px(i) bb1) * (e (-px(i) (bb1 + 1) / cc1) 理論產(chǎn)油量Qo公式實(shí)際散點(diǎn)s(i, 0) = px(i)s(i, 1) = q0(i)預(yù)測曲線s(i, 2) = px(i)s(i, 3) = py(i)Next i散點(diǎn) Dim Index As Integer Index = 1

51、設(shè)置圖標(biāo)類型為非自動選擇，如此能夠設(shè)置為自己需要的類型 MSChart1.Plot.SeriesCollection(Index).ShowLine = False 去除連線！ MSChart1.Plot.SeriesCollection.Item(Index).SeriesMarker.Auto = False With MSChart1.Plot.SeriesCollection.Item(Index).DataPoints.Item(-1).Marker .Visible = True 散點(diǎn)的點(diǎn)是否可見 .Size = 120 散點(diǎn)的大小 .Style = VtMarkerStyleFi

52、lledDiamond 選擇菱形的 .FillColor.Automatic = False 顏色不自動填充 .FillColor.Set 400, 15, 400 設(shè)置為紫紅色 .Pen.VtColor.Set 400, 15, 400 同樣那個也必須設(shè)置，否則中間有默認(rèn)的顏色 End With 橫縱坐標(biāo)表示W(wǎng)ith Form1.MSChart1 .Plot.Axis(VtChAxisIdX, 0).AxisTitle = t .Plot.Axis(VtChAxisIdY, 0).AxisTitle = Qo .Title.Text = Qo與t關(guān)系曲線End WithMSChart1.Ch

53、artData = s 顯示圖形End SubPrivate Sub d2_Click()累積產(chǎn)油量Np與時刻t的關(guān)系曲線MSChart2.Visible = TrueReDim px(yt), py(yt), s(yt, 3)st = yt / 500For i = 1 To mnpx(i) = ipy(i) = aa1 * cc1 * (1 - (e (-px(i) (bb1 + 1) / cc1) / (bb1 + 1) 累計(jì)產(chǎn)油量Np公式實(shí)際散點(diǎn)s(i, 0) = px(i)s(i, 1) = np(i)預(yù)測曲線s(i, 2) = px(i)s(i, 3) = py(i)Next i

54、散點(diǎn) Dim Index As Integer Index = 1 設(shè)置圖標(biāo)類型為非自動選擇，如此能夠設(shè)置為自己需要的類型 MSChart2.Plot.SeriesCollection(Index).ShowLine = False 去除連線！ MSChart2.Plot.SeriesCollection.Item(Index).SeriesMarker.Auto = False With MSChart2.Plot.SeriesCollection.Item(Index).DataPoints.Item(-1).Marker .Visible = True 散點(diǎn)的點(diǎn)是否可見 .Size =

55、120 散點(diǎn)的大小 .Style = VtMarkerStyleFilledDiamond 選擇菱形的 .FillColor.Automatic = False 顏色不自動填充 .FillColor.Set 400, 15, 400 設(shè)置為紫紅色 .Pen.VtColor.Set 400, 15, 400 同樣那個也必須設(shè)置，否則中間有默認(rèn)的顏色 End With With Form1.MSChart2 .Plot.Axis(VtChAxisIdX, 0).AxisTitle = t .Plot.Axis(VtChAxisIdY, 0).AxisTitle = Np .Title.Text =

56、 Np與t關(guān)系曲線End With 顯示圖形MSChart2.Column = 2MSChart2.Plot.UniformAxis = False 自動調(diào)節(jié)圖形大小MSChart2.ChartData = sEnd SubPrivate Sub d3_Click()含水率fw與時刻t的關(guān)系曲線MSChart3.Visible = TrueReDim px(yt), py(yt), s(yt, 3)st = yt / 500For i = 1 To mnpx(i) = ipy(i) = 100 * (1 - 10 (-(aa + ab * ynp(i) + Log(2.303 * ab) /

57、Log(10) 含水率fw公式If py(i) 0 Then py(i) = 0實(shí)際散點(diǎn)s(i, 0) = px(i)s(i, 1) = fw(i)預(yù)測曲線s(i, 2) = px(i)s(i, 3) = py(i)Next i散點(diǎn) Dim Index As Integer Index = 1 設(shè)置圖標(biāo)類型為非自動選擇，如此能夠設(shè)置為自己需要的類型 MSChart3.Plot.SeriesCollection(Index).ShowLine = False 去除連線！ MSChart3.Plot.SeriesCollection.Item(Index).SeriesMarker.Auto =

58、False With MSChart3.Plot.SeriesCollection.Item(Index).DataPoints.Item(-1).Marker .Visible = True 散點(diǎn)的點(diǎn)是否可見 .Size = 120 散點(diǎn)的大小 .Style = VtMarkerStyleFilledDiamond 選擇菱形的 .FillColor.Automatic = False 顏色不自動填充 .FillColor.Set 400, 15, 400 設(shè)置為紫紅色 .Pen.VtColor.Set 400, 15, 400 同樣那個也必須設(shè)置，否則中間有默認(rèn)的顏色 End With Wi

59、th Form1.MSChart3 .Plot.Axis(VtChAxisIdX, 0).AxisTitle = t .Plot.Axis(VtChAxisIdY, 0).AxisTitle = fw .Title.Text = fw與t關(guān)系曲線End With 顯示圖形MSChart3.Column = 2MSChart3.Plot.UniformAxis = False 自動調(diào)節(jié)圖形大小MSChart3.ChartData = sEnd SubPrivate Sub z1_Click()Weibull曲線Form2.MSChart1.Visible = TrueReDim s(ub, 3)

60、s(0, 0) = 0s(0, 1) = a1s(1, 0) = t(ub) (bb1 + 1)s(1, 1) = a1 + a2 * t(ub) (bb1 + 1)For i = 0 To ub Step 1s(i, 0) = t(i) (bb1 + 1)s(i, 1) = Log(q0(i) / (t(i) bb1) / Log(10) 實(shí)際散點(diǎn)s(i, 2) = t(i) (bb1 + 1)s(i, 3) = a1 + a2 * t(i) (bb1 + 1)Next i橫縱坐標(biāo)表示W(wǎng)ith Form2.MSChart1 .Plot.Axis(VtChAxisIdX, 0).AxisTit