【正文】 TMAX = t(i) If NP(i) NPMAX Then NPMAX = NP(i) Next i (2, Int(NPMAX) / 30000 + 2)(Int(TMAX) / 8 + 2, 2) = 2 (0, 0)(0, Int(NPMAX) / 30000 + 1) (0, 0)(Int(TMAX) / 8, 0) = Int(TMAX) / 8 + : = : t/年 = : = Int(NPMAX) / 30000: NP/萬噸 For i = 1 To Int(TMAX) / 8 (i, 0)(i, ) = i : = : i * 8 Next i For i = 0 To Int(NPMAX) / 30000 + (0, i)(, i) = 2: = i: i * 東北石油大學課程設(shè)計（報告） 44 30000 Next i For i = 1 To m (t(i 1) / 8, NP(i 1) / 30000)(t(i) / 8, NP(i) / 30000), vbGreen (t(i) / 8, NP(i) / 30000)Step(, ), vbRed, BF Next i t(0) = t(1) 39。輸出縱坐標值 (0, i)(, i) = : = i: i Next i For i = m To n 39。畫擬合線 lgw = A * lgQ + B (lgQ, lgw), vbGreen Next lgQ End Sub Private Sub Command3_Click() 39。 End If If B 0 Then 東北石油大學課程設(shè)計（報告） 29 +。 ) = 。 2. 陳元千，趙慶飛，油氣田儲采比變化關(guān)系的研究，斷塊油氣田 [J]， （ 6），2326； 3. 陳元千，胡建國，預測油氣田產(chǎn)量和可采儲量 weibull 模型，新疆石油地質(zhì)，（ 3）， 250255； 4. 陳元千，油田可采儲量計算方法，新疆石油地質(zhì) [J]， （ 2）， 130137； 5. 陳元千，油氣藏工程實踐，石油工業(yè)出版社 [M],。其理論直線關(guān)系為 lg lgA B t? ?? 經(jīng)對比可知理論描述的全程儲采比的變化規(guī)律合理有效。 00 00 00 00 0010 .00 012 .00014 .00016 .00018 19 20 21 22儲采比時間系列 1圖 55 穩(wěn)產(chǎn)階段儲采比與開發(fā)時間的關(guān)系曲線 遞減階段儲采比與產(chǎn)量的關(guān)系曲線，如圖 56所示。 若產(chǎn)量上升階段的產(chǎn)量隨時間的關(guān)系符合如下指數(shù)遞減關(guān)系： DtiQ Qe? （ 35） 則產(chǎn)量上升階段的累積產(chǎn)量可表示為： 東北石油大學課程設(shè)計（報告） 3 p 0tN Qdt?? （ 36） 將（ 36）式代入（ 37）式得： i D t D tip 0 ( ) ( 1 )t N e d D t eDD? ? ?? (37) 將（ 38）式代入儲采比的定義式中得到下式： DtR P R i ( 1 )N N N Q eQ Q D Q? ??? ? ? (38) 對于產(chǎn)量上升階段 i 0Q? ，因此有： RNQ?? (39) 對（ 39）式兩邊取常用對數(shù)得： Rlg lg lgNQ? ?? (310) 則（ 310）式可改寫為： lg lgAQ? ?? (311) 其中： RlgAN? 。 確定 a， b 和 c 常數(shù)方法 ................................................................ 錯誤 !未定義書簽。 4 程序設(shè)計框圖 ............................................................................................ 錯誤 !未定義書簽。 穩(wěn)產(chǎn)階段儲采比變化特征 在穩(wěn)產(chǎn)階段某一開發(fā)時間的儲采比表示為： SR P I PSON N NQ? ??? (312) 在穩(wěn)產(chǎn)階段某一開發(fā)時間的累積產(chǎn)量可表示為： SSP O 1()N Q t t?? (313) 將（ 313）式代入（ 312）式可得，在穩(wěn)產(chǎn)階段儲采比與開發(fā)時間 的關(guān)系式為： SR P I O 1 R P I1SSOO()N N Q t t NN tt? ? ? ? ?? ? ? ? (314) 則（ 314）式可寫成： Et??? (315) 其中： R PI 1SONNEtQ???。 y = R178。 求威布爾模型中的 ,ab和 c 的值，及利用 ,ab和 c 的值求得可采儲量 RN 值。 東北石油大學課程設(shè)計（報告） 22 附錄（計算機程序） 一、 程序運行界面 數(shù)據(jù)導入及計算結(jié)果，如附圖 1 所示。 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)) / 東北石油大學課程設(shè)計（報告） 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=。穩(wěn)產(chǎn)階段 Dim m As Integer, n As Integer = True Print xx = t: yy = w m = InputBox(輸入數(shù)據(jù)的起點 , 穩(wěn)產(chǎn)階段 , 18) n = InputBox(輸入數(shù)據(jù)的終點 , 穩(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) 東北石油大學課程設(shè)計（報告） 33 npi = NP(m 1) t1 = t(m 1) 穩(wěn)產(chǎn)階段的平均產(chǎn)量為 。畫散點 (x(i), y(i))Step(, ), vbRed, BF Next i For lgQ = XMIN To XMAX + Step 39。繪制儲采比與時間關(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(TMAX) / 8, 0) = Int(TMAX) / 8 + 1: = : t/年 = : = Int(wMAX) / 300 + 1: w For i = 1 To Int(TMAX) / 8 (i, 0)(i, ) = i : = : i * 8 東北石油大學課程設(shè)計（報告） 45 Next i For i = 0 To Int(wMAX) / 300 + (0, i)(, i) = : = i: i * 300 Next i F。 Format(qb, ) Else 全程求得 b值為 。 t1 產(chǎn)量上升階段結(jié)束時累積產(chǎn)量為 。導入數(shù)據(jù)，求解累積產(chǎn)量和儲采比 Dim m As Integer Dim i As Integer = False YZNR = InputBox(輸入已知 NR值 , 數(shù)據(jù) , ) Print 時間 , 年份 , 產(chǎn)量 , 累積產(chǎn)量 , 儲采比 m = InputBox(輸入已知數(shù)據(jù)的個數(shù) , 數(shù)據(jù) , 49) If m 1 Then On Error Resume Next ReDim t(m), nian(m), Q(m), NP(m), w(m) Open E:\已知數(shù)據(jù) .txt For Input As 1 Do Until EOF(1) For j = 1 To 49 If j 50 Then ReDim Preserve t(j), Q(j) Input 1, t(j), nian(j), Q(j) End If Next j Loop 東北石油大學課程設(shè)計（報告） 30 Close 1 For i = 1 To m NP(i) = NP(i 1) + Q(i) w(i) = (YZNR NP(i)) / Q(i) Print t(i), nian(i), Q(i), Format(NP(i), ), Format(w(i), ) If i = 49 Th