【正文】 i)), 2) 18 yfw(i) = Round(100 * (1 10 ^ ((aa + ab * ynp(i) + Log( * ab) / Log(10)))), 2) If yfw(i) 0 Then yfw(i) = 0 Next i mn = 50 Rem 表格輸出 (0) = 1000 39。最高年產(chǎn)量 = tm 39。參數(shù) c = afw 39。參數(shù) a = bb1 39。截距 = rmax 39?？刹蓛α?NR 公式， NR=a*c/(b+1) = a2 39。最高年產(chǎn)量對應(yīng)時(shí)間 tm 公式， tm=(b*c/(b+1))^(1/(b+1)) qmax = aa1 * (bb1 * cc1 / (bb1 + 1)) ^ (bb1 / (bb1 + 1)) * Exp(bb1 / (bb1 + 1)) 39。r 取最大值 rmax = r bmax = b 17 a2 = ((bn + 1) * sxy sx * sy) / ((bn + 1) * sx2 sx ^ 2) a1 = (sx2 * sy sx * sxy) / ((bn + 1) * sx2 sx ^ 2) End If Next b aa1 = 10 ^ a1 bb1 = bmax cc1 = 1 / ( * a2) afw = 100 * (1 10 ^ ((aa + ab * aa1 * cc1 / (bb1 + 1) + Log( * ab) / Log(10)))) 39。最大相關(guān)系數(shù) Rem weibull 參數(shù)計(jì)算 bn = 12 l = ub bn bmax = 0 rmax = 0 16 For b = 1 To 2 Step sx = 0: sy = 0: sxy = 0: sx2 = 0: sy2 = 0 For i = l To l + bn t(i) = i + 1 Next i For 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 + qt(i) ^ 2 Next i r = Abs(((bn + 1) * sxy sx * sy) / (((bn + 1) * sx2 sx ^ 2) ^ * ((bn + 1) * sy2 sy ^ 2) ^ )) 39。截距，參數(shù) A = ab 39。求相關(guān)系數(shù) r,趨勢線 sx = 0: sy = 0: sxy = 0: sx2 = 0: sy2 = 0 l = ub an For i = l To an + l sx = sx + np(i) Next i 15 For i = l To an + l sy = sy + lp(i) Next i For i = l To an + l sxy = sxy + np(i) * lp(i) Next i For i = l To an + l sx2 = sx2 + np(i) ^ 2 Next i For i = l To an + l sy2 = sy2 + lp(i) ^ 2 Next i a2 = ((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) ^ * ((an + 1) * sy2 sy ^ 2) ^ ) 39。課程設(shè)計(jì)同時(shí)又是一門講道課，一門思辨課，給了我許多道，給了我許多思，得到很大的提升空間。 （ 5） 這次課程設(shè)計(jì)，使我認(rèn)識到學(xué)好計(jì)算機(jī)的重要性與實(shí)踐操作的必要性， 在老師的指導(dǎo)下和與 同學(xué)的討論中，不斷地發(fā)現(xiàn)問題、不斷地解決問題、不斷地領(lǐng)悟、不斷地獲取，使自己的程序不斷完善。同時(shí)，提高了使用計(jì)算機(jī)技應(yīng)用的能力以及解決實(shí)際問題的工程應(yīng)用能力，使理論與實(shí)踐相結(jié)合。雖然在課程設(shè)計(jì)的過程中遇到了一些問題，但經(jīng)過反復(fù)的思考與檢查找出原因所在并進(jìn)行改正。 圖 45 含水率與時(shí)間關(guān)系曲線 12 認(rèn)識與結(jié)論 通過本次課程設(shè)計(jì)，得出如下認(rèn)識與結(jié)論： （ 1）通過此次 油田開發(fā)的 課程設(shè)計(jì)， 我 基本掌握了威布爾（ Weibull）模 型與乙型水驅(qū)曲線聯(lián)解法，并能初步對油田進(jìn)行預(yù)測分析。 圖 43 產(chǎn)量與時(shí)間關(guān)系曲線 11 按照（ 319）式的關(guān)系在同一坐標(biāo)系中繪制實(shí)際累計(jì)產(chǎn)量與預(yù)測累 計(jì)產(chǎn)量及其對應(yīng)時(shí)間 t 之間的關(guān)系曲線，繪于圖 44。 預(yù)測的理論產(chǎn)油量、累積產(chǎn)油量和含水率，與實(shí)際產(chǎn)油量、累積產(chǎn)油量和含水率進(jìn)行對比。經(jīng)線性回歸求得直線的截距 ? =；直線的斜率 ? =；參數(shù) a=， b=， c=， 直線的相關(guān)系數(shù) r=。經(jīng)線性回歸求得直線的截距 A=；直線的斜率 B=；直線的相關(guān)系數(shù) r=。 在確定預(yù)測模型常數(shù) a、 b、 c、 NR時(shí)，其值是否正確可靠，要利用（ 318）式、（ 319）式、（ 328）式預(yù)測的理論產(chǎn)油量、累積產(chǎn)油量和含水率，與實(shí)際產(chǎn)油量、累積產(chǎn)油量和含水率進(jìn)行對比加以確定，而達(dá)到最佳擬合效果的參數(shù)才是最準(zhǔn)確、可靠的。它的理論推導(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、 c 以及可采儲量 NR的數(shù)值，對（ 318）式可進(jìn)行如下處理： 7 1lglg ??? bb tcatQ o （ 331） 若設(shè)： alg?? （ 332） 1?? （ 333） 則得： 1lg ??? bb ttQ ??o （ 334） 根據(jù)實(shí)際的開發(fā)數(shù)據(jù)，首先利用（ 334）式進(jìn)行線性試差求解 ， 根據(jù)最大線性相關(guān)系數(shù)求出 b，然后利用最小二乘法求得 α 和 β。因此，可以將（ 33）式再改寫為： 1 ( / )R e tNQt ????? ??? （ 37） 為了確定最高年產(chǎn)量發(fā)生的時(shí)間，有（ 37）式對時(shí)間 t 求導(dǎo)數(shù)得： ? ?2 ( / )Rd 1ed tN