【正文】 橫縱坐標(biāo)表示 With .(VtChAxisIdX, 0).AxisTitle = Np .(VtChAxisIdY, 0).AxisTitle = Log(Lp) . = 乙型水驅(qū)規(guī)律曲線 End With = 2 = False = s End Sub 28 結(jié)果 運(yùn)行 界面 如下： 29 數(shù)據(jù)輸入 讀取數(shù)據(jù) 、 計(jì)算數(shù)據(jù) ， 出現(xiàn)界面如下： 30 曲線繪制 Qo 與 t 曲線繪制 、 Np 與 t 曲線繪制 、 fw 與 t 曲線繪制 ， 數(shù)據(jù)輸入 輸出數(shù)據(jù)到表格， 出現(xiàn)界面如下： 31 曲線模型 Weibull曲線 、 乙型 水驅(qū)規(guī)律曲線 ，出現(xiàn)界面如下： 32 XXXX 石油大學(xué)課程設(shè)計(jì)成績評價(jià)表 課程名稱 石油工程課程設(shè)計(jì) 題目名稱 威布爾模型與乙型水驅(qū)曲線的聯(lián)解法 學(xué)生姓名 學(xué)號 指導(dǎo)教師姓名 職稱 講師 序號 評價(jià)項(xiàng)目 指 標(biāo) 滿分 評分 1 工作量、工作態(tài)度和出勤率 按期圓滿的完成了規(guī)定的任務(wù)，難易程度和工作量符合教學(xué)要求，工作努力，遵守紀(jì)律，出勤率高，工作作風(fēng)嚴(yán)謹(jǐn)，善于與他人合作。設(shè)置圖標(biāo)類型為非自動(dòng)選擇，這樣可以設(shè)置為自己需要的類型 (Index).ShowLine = False 39。散點(diǎn)的點(diǎn)是否可見 .Size = 80 39。顯示圖形 = 2 = False 39。設(shè)置圖標(biāo)類型為非自動(dòng)選擇，這樣可以設(shè)置為自己需要的類型 (Index).ShowLine = False 39。顯示圖形 = 2 = False 39。設(shè)置圖標(biāo)類型為非自動(dòng)選擇，這樣可以設(shè)置為自己需要的類型 (Index).ShowLine = False 39。橫縱坐標(biāo)表示 With .(VtChAxisIdX, 0).AxisTitle = t .(VtChAxisIdY, 0).AxisTitle = Qo . = Qo 與 t 關(guān)系曲線 End With = s 39。設(shè)置圖標(biāo)類型為非自動(dòng)選擇，這樣可以設(shè)置為自己需要的類型 (Index).ShowLine = False 39。輸出數(shù)據(jù)在表 格內(nèi) Rem 曲線繪制 Rem 數(shù)據(jù)輸出 End Sub 21 Private Sub c3_Click() 39。極限含水率 = qmax 39。最高年產(chǎn)量 Qmax 公式， Qmax=a*(b*c/(b+1))^(b/(b+1))*exp(b/(b+1)) NR = aa1 * cc1 / (bb1 + 1) 39。取最大值 r If r ar Then ar = r Next an aa = a1 ab = a2 = aa 39。 （ 3）通過計(jì)算機(jī)的編程學(xué)習(xí)與操作，使我對 VB 語言有了更進(jìn)一步的認(rèn)識和了解。 圖 42 該油田的 weibull 曲線 9 在 確定預(yù)測模型常數(shù) a、 b、 c 后，利用（ 318）式、（ 319）式、（ 328）式計(jì)算油田年產(chǎn)油量、計(jì)算油田累積產(chǎn)油量、計(jì)算含水率如下表： 表 41 數(shù)據(jù)表 時(shí)間（ a） 實(shí)際產(chǎn)量 預(yù)測產(chǎn)量 實(shí)際累計(jì)產(chǎn)量 預(yù)測累計(jì)產(chǎn)量 實(shí)際含水率 預(yù)測含水率 1 0 2 0 3 0 4 0 5 0 6 0 7 0 8 0 9 0 10 11 12 13 14 15 16 17 18 19 20 87 21 2747 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 10 41 42 43 44 45 46 47 48 49 50 51 計(jì)算可采儲量： NR =（ 104t） ，計(jì)算最高年產(chǎn)量發(fā)生的時(shí)間： ， 計(jì)算最高年產(chǎn)量： Qmax=(104t)。該模型具有預(yù)測油田產(chǎn)量、累積產(chǎn)量和可采儲量的功能，其基本關(guān)系式分別為： 1o e xpbb tQ atc????????? （ 318） 1p 1 exp ( )1ba c tNbc???? ? ???? ?? （ 319） R 1acN b? ? （ 320） 乙型水驅(qū)曲線法，首先是由我國著名專家童憲章先生 [2]以經(jīng)驗(yàn)公式的形式，于 1978年提出。該模型的分布密度表示為： )/(1)( ?? ??? xexxf ??? (31) 式中 f(x)—— 威布爾分布的分布密度函數(shù)； x—— 分布變量，根據(jù)實(shí)際問題，分布區(qū)間為 0~∞； α—— 控制分布形態(tài)的形狀參數(shù)； β—— 控制分布峰位和峰值的尺度參數(shù)。 本課題要求學(xué)生對 威布爾（ Weibull） 模型及其與乙型水驅(qū)曲線聯(lián)解關(guān)系式進(jìn)行推導(dǎo)，其結(jié)果包括，對油氣田的年產(chǎn)量、累積產(chǎn)量、含水率、可采儲量、最高年產(chǎn)量、最高年產(chǎn)量發(fā)生的時(shí)間的預(yù)測方法。從而將理論知識和實(shí)際問題相結(jié)合，通過該專題設(shè)計(jì)的訓(xùn)練，加強(qiáng)學(xué)生理論知識運(yùn)用能力，計(jì)算機(jī)技術(shù)應(yīng)用能力及解決實(shí)際問題的工程應(yīng)用能力。 若對 (3l)式進(jìn)行積分，在 x 為 0~∞區(qū)間內(nèi)，可以得到 Weibu11 的分布函數(shù)值等于 1，推證如下： 1)/()()(0)/(0)/(00??????????? ?????????????????xxxxxxxxfxFededed )/(1 （ 32） 為將 Weibull 分布模型用于油氣田開發(fā)指標(biāo)的預(yù)測，將 (3l)式改寫為 )/(1 ?? ??? ttCQ ??? e (33) 式中 Q—— 油氣田的年產(chǎn)量， 104t/a（油）或 108m3/a（氣）； t —— 油氣田的開發(fā)時(shí)間， a； C—— 由 Weibull 分布模型轉(zhuǎn)換為油氣田開發(fā)實(shí)用模型的模型轉(zhuǎn)換常數(shù)。它的理論推導(dǎo)由文獻(xiàn) [3]完成，其基本關(guān)系式為： pp BNAL ??lg （ 321） 由 （ 321）式對時(shí)間 t 求導(dǎo)數(shù)得： tNBtLL dddd ppp ??30 1 （ 322） 已知：wopdd tL ??；opdd QtN ?；woow R ? 故由（ 322）式得： 6 )1(3 0 1 wop RBL ?? （ 323） 將（ 323）式代入（ 321）式得： BBNAR 3 0 )1lg ( pwo ???? （ 324） 取經(jīng)濟(jì)極限水 油比 (Rwo) L，由（ 24）式得到預(yù)測油田可采儲量的關(guān)系式： ? ? B BARN )()(1lg ???? LwoR （ 325） 已知水油比與含水率的關(guān)系為： wwwo ffR ??1 （ 326） 將（ 326）式代入（ 324）式得： )(101 BBNAf ????? pw （ 327） 將（ 319）式代入（ 327）式得： ?????????? ????????? ???????? ???????? BctbacBAbf )e x p (111101w （ 328） 當(dāng)由（ 318）式和（ 328）式得到預(yù)測的產(chǎn)油量和含水率之后，可由下面的公式分別預(yù)測油田的產(chǎn)水量和產(chǎn)液量： )1( wwow ff ?? （ 329） )1 1( woL f ?? （ 330） 最高年產(chǎn)量發(fā)生的時(shí)間 tm的計(jì)算公式如下： 11m 1?????????bbbct 最高年產(chǎn)量 Qmax為： 11m a x e1 ????????? ?? bbb bbbcaQ 模型的求解方法 為了確定預(yù)測模型的模型常數(shù) a、 b、