【正文】 東北石油大學(xué)課程設(shè)計(jì)（報(bào)告） 12 25 1977 26 1978 27 1979 6191 28 1980 29 1981 30 1982 31 1983 32 1984 33 1985 34 1986 35 1987 36 1988 37 1989 38 1990 39 1991 174907 40 1992 41 1993 42 1994 43 1995 44 1996 45 1997 46 1998 47 1999 48 20xx 190 49 20xx 東北石油大學(xué)課程設(shè)計(jì)（報(bào)告） 13 繪制產(chǎn)量、累積產(chǎn)量和儲(chǔ)采比隨時(shí)間的變化曲線： 圖 51產(chǎn)量與時(shí)間關(guān)系曲線 東北石油大學(xué)課程設(shè)計(jì)（報(bào)告） 14 圖 52 累積產(chǎn)量和時(shí)間的關(guān)系曲線 圖 53 儲(chǔ)采比和時(shí)間的關(guān)系曲線 ，如圖 54所示。 y = R178。 0 . 0 0 02 0 0 . 0 0 04 0 0 . 0 0 06 0 0 . 0 0 08 0 0 . 0 0 01 0 0 0 . 0 0 01 2 0 0 . 0 0 01 4 0 0 . 0 0 01 6 0 0 . 0 0 01 8 0 0 . 0 0 00 1 0 0 0 2 0 0 0 3 0 0 0 4 0 0 0 5 0 0 0 6 0 0 0 7 0 0 0 8 0 0 0儲(chǔ)采比時(shí)間產(chǎn)量上升階段儲(chǔ)采比與產(chǎn)量的變化曲線系列 1 圖 59 產(chǎn)量上升階段儲(chǔ)采比與產(chǎn)量的變化曲線 穩(wěn)產(chǎn)階段儲(chǔ)采比與開發(fā)時(shí)間的關(guān)系曲線，如圖 510 所示。符合其理論直線關(guān)系為 lg lgAQ? ?? 。其理論直線關(guān)系為 lg lg Q? ? ??? ，經(jīng)對(duì)比可知理論描述的遞減階段儲(chǔ)采比的變化關(guān)系較為合理。 求威布爾模型中的 ,ab和 c 的值，及利用 ,ab和 c 的值求得可采儲(chǔ)量 RN 值。威布爾預(yù)測(cè)模型預(yù)測(cè)注水開發(fā)油田可采儲(chǔ)量與實(shí)際值相當(dāng)誤差較小，說明此法預(yù)測(cè)的可采儲(chǔ)量較為準(zhǔn)東北石油大學(xué)課程設(shè)計(jì)（報(bào)告） 20 確。因此我們應(yīng)該用威布爾模型和乙型水驅(qū)曲線聯(lián)解來預(yù)測(cè)油田開發(fā)指標(biāo)。 全程： 全程儲(chǔ)采比和時(shí)間呈雙對(duì)數(shù)的直線關(guān)系。 東北石油大學(xué)課程設(shè)計(jì)（報(bào)告） 22 附錄（計(jì)算機(jī)程序） 一、 程序運(yùn)行界面 數(shù)據(jù)導(dǎo)入及計(jì)算結(jié)果，如附圖 1 所示。 附圖 4 穩(wěn)產(chǎn)階段的計(jì)算結(jié)果和圖形輸出 遞減階段的計(jì)算結(jié)果和 圖形輸出，如附圖 5 所示。畫圖使用 XMIN = x(m): YMIN = y(m) 東北石油大學(xué)課程設(shè)計(jì)（報(bào)告） 27 For i = m To n If XMAX x(i) Then XMAX = x(i) If YMAX y(i) Then YMAX = y(i) If XMIN x(i) Then XMIN = x(i) If YMIN y(i) Then YMIN = y(i) Next i For i = m To n ER(1, i) = x(i) ER(2, i) = y(i) ER(3, i) = ER(1, i) * ER(2, i) ER(4, i) = ER(1, i) ^ 2 ER(5, i) = ER(2, i) ^ 2 Print xx。 Format(ER(1, i), ), yy。 Format(ER(2, i), ) If i = 49 Then CurrentX = / 2 CurrentY = 0 ElseIf i 49 Then CurrentX = / 2 End If ER(1, m 1) = ER(1, m 1) + ER(1, i) ER(2, m 1) = ER(2, m 1) + ER(2, i) ER(3, m 1) = ER(3, m 1) + ER(3, i) ER(4, m 1) = ER(4, m 1) + ER(4, i) ER(5, m 1) = ER(5, m 1) + ER(5, i) Next i ER(1, m 1) = ER(1, m 1) / (n m + 1) ER(2, m 1) = ER(2, m 1) / (n m + 1) R = (ER(3, m 1) (n m + 1) * ER(1, m 1) * ER(2, m 1)) / 東北石油大學(xué)課程設(shè)計(jì)（報(bào)告） 28 _ Sqr((ER(4, m 1) (n m + 1) * ER(1, m 1) ^ 2) * (ER(5, m 1) (n m + 1) * ER(2, m 1) ^ 2)) A = ((n m + 1) * ER(1, m 1) * ER(2, m 1) ER(3, m 1)) / ((n m + 1) * ER(1, m 1) ^ 2 ER(4, m 1)) B = ER(2, m 1) A * ER(1, m 1) If i 49 Then 擬合公式為： 39。 =。 Format(B, ) ElseIf B 0 Then Format(B, ) End If 線性相關(guān)系數(shù) R=。 =。 Format(B, ) ElseIf B 0 Then Format(B, ) End If 線性相關(guān)系數(shù) R=。控制按鈕可用性 = True = True = True End Sub Private Sub Command2_Click() 39。繪圖 (2, Int(YMAX) + 2)(Int(XMAX) + 2, 2) 39。輸出橫坐標(biāo)值 (i, 0)(i, ) = i: = : i Next i For i = 0 To Int(YMAX) + 1 39。穩(wěn)產(chǎn)階段 Dim m As Integer, n As Integer = True Print xx = t: yy = w m = InputBox(輸入數(shù)據(jù)的起點(diǎn) , 穩(wěn)產(chǎn)階段 , 18) n = InputBox(輸入數(shù)據(jù)的終點(diǎn) , 穩(wěn)產(chǎn)階段 , 22) ReDim x(m To n), y(m To n) For i = m To n x(i) = t(i) y(i) = w(i) QS = QS + Q(i) Next i QS = QS / (n m + 1) 東北石油大學(xué)課程設(shè)計(jì)（報(bào)告） 33 npi = NP(m 1) t1 = t(m 1) 穩(wěn)產(chǎn)階段的平均產(chǎn)量為 。遞減階段 = True Dim m As Integer, n As Integer Print xx = lgQ: yy = lgw m = InputBox(輸入數(shù)據(jù)的起點(diǎn) , 遞減階段 , 23) n = InputBox(輸入數(shù)據(jù)的終點(diǎn) , 遞減階段 , 49) ReDim x(m To n), y(m To n) For i = m To n x(i) = Log(Q(i)) / Log(10) y(i) = Log(w(i)) / Log(10) Next i Call ercheng(x(), y(), m, n) ((Int(XMAX) + 2) * 2, (Int(YMAX) + 2) * 2)(2, 2) 39。 Format(qbb, ) 遞減階段所求的 n值為 。畫橫坐標(biāo) 東北石油大學(xué)課程設(shè)計(jì)（報(bào)告） 37 (0, 0)(Int(XMAX + 1), 0) 39。畫散點(diǎn) (x(i), y(i))Step(, ), vbRed, BF Next i For lgQ = XMIN To XMAX + Step 39。畫圖用的刻度 (Int(XMIN) 2, Int(YMAX) + 2)(Int(XMAX) + 2, Int(YMIN 2)) = 2 (Int(XMIN 1), Int(YMIN 1))(Int(XMIN 1), Int(YMAX + 1)) (Int(XMIN 1), Int(YMIN 1))(Int(XMAX + 1), Int(YMIN 1)) 東北石油大學(xué)課程設(shè)計(jì)（報(bào)告） 39 = : = Int(YMIN) : xx = Int(XMIN) : = Int(YMAX) + : yy For i = Int(XMIN) 1 To Int(XMAX) + 1 Step buc (i, Int(YMIN) 1)(i, Int(YMIN) 1 + ) = i: = Int(YMIN) : i Next i For i = Int(YMIN) 1 To Int(YMAX) + 1 (Int(XMIN) 1, i)(Int(XMIN) 1 + , i) = Int(XMIN) 2: = i: i Next i For i = m To n (x(i), y(i))Step(, ), vbRed, BF Next i For lgQ = XMIN To XMAX + Step lgw = A * lgQ + B (lgQ, lgw), vbGreen Next lgQ QAa = 10 ^ B c = 1 / ( * Abs(A)) NR = QAa * c / (qb + 1) If i 49 Then 全程求得 a值為 。 Format(QAa, ) 全程求得 c值為 。返回窗體 1 End Sub Private Sub Command2_Click() 39。繪制儲(chǔ)采比與時(shí)間關(guān)系曲線 w(0) = w(1) TMAX = t(1) wMAX = w(1) For i = 1 To m If t(i) TMAX Then TMAX = t(i) If w(i) wMAX Then wMAX = w(i) Next i (2, Int(wMAX) / 300 + 2)(Int(TMAX) / 8 + 2, 2) = 2 (0, 0)(0, Int(wMAX) / 300 + ) (0, 0)(Int(TM