【文章內容簡介】 9）式、（ 328）式計算油田年產油量、計算油田累積產油量、計算含水率如下表： 表 41 數(shù)據(jù)表 時間（ a） 實際產量 預測產量 實際累計產量 預測累計產量 實際含水率 預測含水率 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） ，計算最高年產量發(fā)生的時間： ， 計算最高年產量： Qmax=(104t)。 預測的理論產油量、累積產油量和含水率，與實際產油量、累積產油量和含水率進行對比。 按照（ 318）式的關系在同一坐標系中繪制實際產量與預測產量及其對應時間 t 之間的關系曲線， 繪于圖 43。 圖 43 產量與時間關系曲線 11 按照（ 319）式的關系在同一坐標系中繪制實際累計產量與預測累 計產量及其對應時間 t 之間的關系曲線，繪于圖 44。 圖 44 累計產量與時間關系曲線 按照（ 328）式的關系在同一坐標系中繪制實際含水率與預測含水率及其對應時間 t之間的關系曲線，繪于圖 45。 圖 45 含水率與時間關系曲線 12 認識與結論 通過本次課程設計，得出如下認識與結論： （ 1）通過此次 油田開發(fā)的 課程設計， 我 基本掌握了威布爾（ Weibull）模 型與乙型水驅曲線聯(lián)解法，并能初步對油田進行預測分析。 （ 2） 通過這次的課程設計，培養(yǎng)了我獨立思考和動手操作的能力，也使我在日常的學習與討論當中增長了許多課外知識，認識到數(shù)學模型對油氣田的產量預測的重要性。雖然在課程設計的過程中遇到了一些問題，但經過反復的思考與檢查找出原因所在并進行改正。 （ 3）通過計算機的編程學習與操作，使我對 VB 語言有了更進一步的認識和了解。同時，提高了使用計算機技應用的能力以及解決實際問題的工程應用能力，使理論與實踐相結合。 （ 4） 現(xiàn)在能熟練掌握、運用威布爾（ Weibull）模型與乙型 水驅曲線聯(lián)解法的原理、推導及公式，可以對油氣田的年產量、累計產量、含水率、可采儲量、最高年產量、最高年產量發(fā)生的時間以及與最高年產量對應的累積產量進行預測。 （ 5） 這次課程設計，使我認識到學好計算機的重要性與實踐操作的必要性， 在老師的指導下和與 同學的討論中，不斷地發(fā)現(xiàn)問題、不斷地解決問題、不斷地領悟、不斷地獲取，使自己的程序不斷完善。 （ 6）課程設計誠然是一門專業(yè)課，它給了我很多的專業(yè)知識以及專業(yè)上的提升。課程設計同時又是一門講道課，一門思辨課，給了我許多道，給了我許多思，得到很大的提升空間。 13 參考文獻 [1]陳元千 .預測油氣田產量的 Weibull 模型 [J].新疆石油地質， 1995， 16（ 3）： 250255. [2]童憲章 .天然水驅和人工注水油藏的統(tǒng)計規(guī)律探討 [J].石油勘探與開發(fā)， 1978， 4（ 6）：3864. [3] 陳元千 .一種新型水驅曲線關系式的推導及應用 [J].石油學報， 1993， 14（ 2）： 6573. [4] 陳元千 .水驅曲線關系式的推導 [J]. 石油學報， 1985， 6（ 2）： 6978. 14 附錄：計算機程序與結果 計算 機 程序 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。求相關系數(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。最大相關系數(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。最大相關系數(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 *