【導(dǎo)讀】假設(shè)油田產(chǎn)量遞減符合雙曲遞減，推導(dǎo)遞減階段儲(chǔ)采比與產(chǎn)量之間的關(guān)系；假設(shè)油田產(chǎn)量可由威布爾預(yù)測(cè)模型給出，推導(dǎo)全程儲(chǔ)采比與開(kāi)發(fā)時(shí)間的關(guān)系；產(chǎn)量、累積產(chǎn)量和儲(chǔ)采比隨時(shí)間的變化曲線；某油田或區(qū)塊的可采儲(chǔ)量和產(chǎn)量歷史數(shù)據(jù)。容，編制相應(yīng)軟件，并提交規(guī)范設(shè)計(jì)報(bào)告。[1]陳元千.油藏工程實(shí)踐[M].石油工業(yè)出版社.油田的穩(wěn)產(chǎn)形式密切相關(guān)。并通過(guò)實(shí)例運(yùn)用此方法和驗(yàn)證此方法的有效性。

　　

【正文】 出縱坐標(biāo)值 (0, i)(, i) = : = i: i Next i For i = m To n 39。畫(huà)散點(diǎn) (x(i), y(i))Step(, ), vbRed, BF Next i For lgQ = XMIN To XMAX + Step 39。畫(huà)擬合線 lgw = A * lgQ + B (lgQ, lgw), vbGreen Next lgQ qb = Abs(A) If n 49 Then 全程求得 b值為 。 Format(qb, ) Else 全程求得 b值為 。 Format(qb, ) End If End Sub 東北石油大學(xué)課程設(shè)計(jì)（報(bào)告） 38 Private Sub Command6_Click() Dim i As Double, j As Integer Dim m As Integer, n As Integer Dim buc As Single Print qb = InputBox(輸入 b值 , 數(shù)據(jù) , ) xx = t^(b + 1): yy = lgQ/t^b m = InputBox(輸入數(shù)據(jù)的起點(diǎn) , 全程 NR, 2) n = InputBox(輸入數(shù)據(jù)的終點(diǎn) , 全程 NR, 49) ReDim x(m To n), y(m To n) For i = m To n x(i) = t(i) ^ (qb + 1) y(i) = Log(Q(i) / (t(i) ^ qb)) / Log(10) Next i Call ercheng(x(), y(), m, n) buc = Val(InputBox(輸入繪圖步長(zhǎng) , 數(shù)據(jù) , 5000)) 39。畫(huà)圖用的刻度 (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值 為 。 Format(c, ) 全程求得 NR值為 。 Format(NR, ) 東北石油大學(xué)課程設(shè)計(jì)（報(bào)告） 40 Else 全程求得 a值為 。 Format(QAa, ) 全程求得 c值為 。 Format(c, ) 全程求得 NR值為 。 Format(NR, ) End If End Sub Private Sub Form_DblClick() End Sub Private Sub Form_Load() = False = False = False = False = False Print 說(shuō)明 Print 將 已知數(shù)據(jù) .txt。復(fù)制到 E盤(pán)根目錄下，才能運(yùn)行程序 Print 雙擊顯示窗體二，繪制產(chǎn)量，累積產(chǎn)量，儲(chǔ)采比與時(shí)間關(guān)系曲線 Print 根據(jù)已知條件可知： Print 產(chǎn)量上升階段為 1到 17 Print 產(chǎn)穩(wěn)產(chǎn)階段為 18 到 22 Print 產(chǎn)量遞減階段為 23 到 49 Print 全程為 2 到 49 End Sub Form2 代碼 ： Dim t() As Single, nian() As Single, Q() As Single, NP() As Single, w() As Single 東北石油大學(xué)課程設(shè)計(jì)（報(bào)告） 41 Private Sub Command1_Click() 39。返回窗體 1 End Sub Private Sub Command2_Click() 39。讀取數(shù)據(jù) Dim m As Integer Dim i As Integer YZNR = InputBox(輸入已知 NR值 , 數(shù)據(jù) , ) m = InputBox(輸入已知數(shù)據(jù)的個(gè)數(shù) , 數(shù)據(jù) , 49) If m 1 Then ReDim t(m), nian(m), Q(m), NP(m), w(m) Open E:\已知數(shù)據(jù) .txt For Input As 1 j = 0 Do While Not EOF(1) j = j + 1 ReDim Preserve t(j), Q(j) Input 1, t(j), nian(j), Q(j) Loop Close 1 For i = 1 To m NP(i) = NP(i 1) + Q(i) w(i) = (YZNR NP(i)) / Q(i) Next i Else Print 輸入已知數(shù)據(jù)的個(gè)數(shù)必須大于一 End If MsgBox 數(shù)據(jù)已讀取 = True 東北石油大學(xué)課程設(shè)計(jì)（報(bào)告） 42 End Sub Private Sub Command3_Click() m = InputBox(數(shù)據(jù)組數(shù) , 數(shù)據(jù) , 49) TMAX = t(1) 39。繪制產(chǎn)量和時(shí)間曲線 QMAX = Q(1) For i = 1 To m If t(i) TMAX Then TMAX = t(i) If Q(i) QMAX Then QMAX = Q(i) Next i (2, Int(QMAX) / 3000 + 2)(Int(TMAX) / 8 + 2, 2) = 2 (0, 0)(0, Int(QMAX) / 3000 + ) (0, 0)(Int(TMAX) / 8, 0) = Int(TMAX) / 8 + 1: = : t/年 = : = Int(QMAX) / 3000 + 1: Q/萬(wàn)噸 For i = 1 To Int(TMAX) / 8 (i, 0)(i, ) = i : = : i * 8 Next i For i = 0 To Int(QMAX) / 3000 + (0, i)(, i) = : = i: i * 3000 東北石油大學(xué)課程設(shè)計(jì)（報(bào)告） 43 Next i For i = 1 To m (t(i 1) / 8, Q(i 1) / 3000)(t(i) / 8, Q(i) / 3000), vbGreen (t(i) / 8, Q(i) / 3000)Step(, ), vbRed, BF Next i TMAX = t(1) 39。繪制累積產(chǎn)量與時(shí)間曲線 NPMAX = NP(1) For i = 1 To m If t(i) TMAX Then 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/萬(wàn)噸 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 * 東北石油大學(xué)課程設(shè)計(jì)（報(bào)告） 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。繪制儲(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(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 東北石油大學(xué)課程設(shè)計(jì)（報(bào)告） 45 Next i For i = 0 To Int(wMAX) / 300 + (0, i)(, i) = : = i: i * 300 Next i