【文章內(nèi)容簡介】 9）式、（ 328）式計算油田年產(chǎn)油量、計算油田累積產(chǎn)油量、計算含水率如下表： 表 41 數(shù)據(jù)表 時間（ a） 實(shí)際產(chǎn)量 預(yù)測產(chǎn)量 實(shí)際累計產(chǎn)量 預(yù)測累計產(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 計算可采儲量： NR =（ 104t） ，計算最高年產(chǎn)量發(fā)生的時間： ， 計算最高年產(chǎn)量： Qmax=(104t)。 預(yù)測的理論產(chǎn)油量、累積產(chǎn)油量和含水率，與實(shí)際產(chǎn)油量、累積產(chǎn)油量和含水率進(jìn)行對比。 按照（ 318）式的關(guān)系在同一坐標(biāo)系中繪制實(shí)際產(chǎn)量與預(yù)測產(chǎn)量及其對應(yīng)時間 t 之間的關(guān)系曲線， 繪于圖 43。 圖 43 產(chǎn)量與時間關(guān)系曲線 11 按照（ 319）式的關(guān)系在同一坐標(biāo)系中繪制實(shí)際累計產(chǎn)量與預(yù)測累 計產(chǎn)量及其對應(yīng)時間 t 之間的關(guān)系曲線，繪于圖 44。 圖 44 累計產(chǎn)量與時間關(guān)系曲線 按照（ 328）式的關(guān)系在同一坐標(biāo)系中繪制實(shí)際含水率與預(yù)測含水率及其對應(yīng)時間 t之間的關(guān)系曲線，繪于圖 45。 圖 45 含水率與時間關(guān)系曲線 12 認(rèn)識與結(jié)論 通過本次課程設(shè)計，得出如下認(rèn)識與結(jié)論： （ 1）通過此次 油田開發(fā)的 課程設(shè)計， 我 基本掌握了威布爾（ Weibull）模 型與乙型水驅(qū)曲線聯(lián)解法，并能初步對油田進(jìn)行預(yù)測分析。 （ 2） 通過這次的課程設(shè)計，培養(yǎng)了我獨(dú)立思考和動手操作的能力，也使我在日常的學(xué)習(xí)與討論當(dāng)中增長了許多課外知識，認(rèn)識到數(shù)學(xué)模型對油氣田的產(chǎn)量預(yù)測的重要性。雖然在課程設(shè)計的過程中遇到了一些問題，但經(jīng)過反復(fù)的思考與檢查找出原因所在并進(jìn)行改正。 （ 3）通過計算機(jī)的編程學(xué)習(xí)與操作，使我對 VB 語言有了更進(jìn)一步的認(rèn)識和了解。同時，提高了使用計算機(jī)技應(yīng)用的能力以及解決實(shí)際問題的工程應(yīng)用能力，使理論與實(shí)踐相結(jié)合。 （ 4） 現(xiàn)在能熟練掌握、運(yùn)用威布爾（ Weibull）模型與乙型 水驅(qū)曲線聯(lián)解法的原理、推導(dǎo)及公式，可以對油氣田的年產(chǎn)量、累計產(chǎn)量、含水率、可采儲量、最高年產(chǎn)量、最高年產(chǎn)量發(fā)生的時間以及與最高年產(chǎn)量對應(yīng)的累積產(chǎn)量進(jìn)行預(yù)測。 （ 5） 這次課程設(shè)計，使我認(rèn)識到學(xué)好計算機(jī)的重要性與實(shí)踐操作的必要性， 在老師的指導(dǎo)下和與 同學(xué)的討論中，不斷地發(fā)現(xiàn)問題、不斷地解決問題、不斷地領(lǐng)悟、不斷地獲取，使自己的程序不斷完善。 （ 6）課程設(shè)計誠然是一門專業(yè)課，它給了我很多的專業(yè)知識以及專業(yè)上的提升。課程設(shè)計同時又是一門講道課，一門思辨課，給了我許多道，給了我許多思，得到很大的提升空間。 13 參考文獻(xiàn) [1]陳元千 .預(yù)測油氣田產(chǎn)量的 Weibull 模型 [J].新疆石油地質(zhì)， 1995， 16（ 3）： 250255. [2]童憲章 .天然水驅(qū)和人工注水油藏的統(tǒng)計規(guī)律探討 [J].石油勘探與開發(fā)， 1978， 4（ 6）：3864. [3] 陳元千 .一種新型水驅(qū)曲線關(guān)系式的推導(dǎo)及應(yīng)用 [J].石油學(xué)報， 1993， 14（ 2）： 6573. [4] 陳元千 .水驅(qū)曲線關(guān)系式的推導(dǎo) [J]. 石油學(xué)報， 1985， 6（ 2）： 6978. 14 附錄：計算機(jī)程序與結(jié)果 計算 機(jī) 程序 Option Base 0 Const e = Dim a, t(200), qt(200), t0(200), q0(200), qw(200), ye(200), wp(200) Dim b, c, d Dim afw, qmax, tm, NR, x, y, a1, a2, r, rmax, bmax, aa1, bb1, cc1, ub, aa, ab, ar, br, yt, an, bn, mn Dim s(), px(), py() Dim sx, sy, sxy, sx2, sy2 Dim np(200), lp(200) Dim yq0(200), ynp(200), yfw(200), ynr, fw(200) Private Sub c1_Click() yt = 50 Rem 源數(shù)據(jù)輸入 Open + \ For Input As 1 ub = 1 Do Until EOF(1) ub = ub + 1 Input 1, ye(ub), t(ub), q0(ub), qw(ub), np(ub), wp(ub), lp(ub) lp(ub) = Log(lp(ub)) / Log(10) Loop Close 1 Private Sub c5_Click() Rem 乙型參數(shù)計算 For an = ub To 2 Step 1 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。取最大值 r If r ar Then ar = r Next an aa = a1 ab = a2 = aa 39。截距，參數(shù) A = ab 39。斜率，參數(shù) B = ar 39。最大相關(guān)系數(shù) Rem weibull 參數(shù)計算 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。最大相關(guān)系數(shù) If r rmax Then 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 *