【文章內(nèi)容簡介】 8 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ǎn)量上升階段儲采比與產(chǎn)量的變化曲線系列 1 圖 59 產(chǎn)量上升階段儲采比與產(chǎn)量的變化曲線 穩(wěn)產(chǎn)階段儲采比與開發(fā)時間的關(guān)系曲線，如圖 510 所示。 圖 510 穩(wěn)產(chǎn)階段儲采比與開發(fā)時間的關(guān)系曲線 東北石油大學(xué)課程設(shè)計（報告） 18 遞減階段儲采比與產(chǎn)量的關(guān)系曲線，如圖 511 所示。 圖 511 遞減階段儲采比與產(chǎn)量的關(guān)系曲線 1全程儲采比與開發(fā)時間的關(guān)系曲線，如圖 512 所示。 圖 512 全程儲采比與開發(fā)時間的關(guān)系曲線 設(shè)計結(jié)果及分析 產(chǎn)量上升階段： 產(chǎn)量上升階段儲采比和產(chǎn)量呈現(xiàn)雙對數(shù) 的直線關(guān)系，直線關(guān)系式為lg 1. 06 85 lg 5. 44 02Q? ? ? ?。符合其理論直線關(guān)系為 lg lgAQ? ?? 。 穩(wěn)產(chǎn)階段： 東北石油大學(xué)課程設(shè)計（報告） 19 穩(wěn) 產(chǎn) 階 段 儲 采 比 與 時 間 呈 直 線 遞 減 關(guān) 系 ， 直 線 關(guān) 系 式 為 59 8 34 .30 33t? ? ? ?。其理論直線關(guān)系為 Et???，經(jīng)對比可知理論描述的穩(wěn)產(chǎn)階段儲采比的變化特征較為合理。 遞減階段： 遞減階段 儲采比和 產(chǎn)量呈 現(xiàn)雙對 數(shù)的直線 關(guān)系， 直線關(guān) 系式為lg 0. 48 73 1 lg 1. 02 33Q? ??，利用直線斜率可求得遞減指數(shù) ? 。其理論直線關(guān)系為 lg lg Q? ? ??? ，經(jīng)對比可知理論描述的遞減階段儲采比的變化關(guān)系較為合理。并由此推斷得的遞減指數(shù)應(yīng)為正確值。 全程： 全程儲采 比和時間 呈雙對 數(shù)的直 線關(guān)系， 直線關(guān) 系式可 表示為lg lg ? ? ? ?。其理論直線關(guān)系為 lg lgA B t? ?? 經(jīng)對比可知理論描述的全程儲采比的變化規(guī)律合理有效。 求威布爾模型中的 ,ab和 c 的值，及利用 ,ab和 c 的值求得可采儲量 RN 值。 直線關(guān)系為： 5 2 . 7 2 2 91 . 7 2 2 9l g 8 . 0 4 2 1 0 1 . 9 7 0 3Q tt ? ? ? ?。 根據(jù)擬合的直線的截距和斜率可求出： α 1. 970310 10 ? ? ? 511 5 3 9 9 . 3 12 . 3 0 3 2 . 3 0 3 8 . 0 4 2 1 0c ?? ? ??? 由于 ? ，將 ,ab和 c 的值代入求可采儲量 RN 公式得可采儲量為： 4R 9 3 . 3 9 5 3 9 9 . 3 1 1 8 5 1 7 7 . 5 ( 1 0 t )1 2 . 7 2 2 9acN b ?? ? ?? 由威布爾模型預(yù)測可采儲量 RN 值： 4R (10 t )N ? ，實際給的已知4R (10 t )N ? ，相對誤差為 %，這說明威布爾模型預(yù)測可采儲量有效。 計算結(jié)果的 認(rèn)識 根據(jù)各階段用實際數(shù)據(jù)的擬合的直線關(guān)系與相應(yīng)的理論關(guān)系式進(jìn)行對比分析，可知各階段描述的儲采比與時間的變化關(guān)系合理有效。威布爾預(yù)測模型預(yù)測注水開發(fā)油田可采儲量與實際值相當(dāng)誤差較小，說明此法預(yù)測的可采儲量較為準(zhǔn)東北石油大學(xué)課程設(shè)計（報告） 20 確。 威布爾預(yù)測模型能很好的預(yù)測注水開發(fā)油田產(chǎn)量、累積產(chǎn)量、可采儲量、儲采比隨時間的變化關(guān)系，但是不能預(yù)測油田的含水率、產(chǎn)水量、產(chǎn)液量及累積產(chǎn)水量和累積產(chǎn)液量，而這些開發(fā)指標(biāo)卻 是水驅(qū)開發(fā)油田所需要預(yù)測的。乙型水驅(qū)曲線法是油藏工程中重要的預(yù)測方法，但它只能預(yù)測累積產(chǎn)水量和累積產(chǎn)油量之間的關(guān)系，卻不能預(yù)測開發(fā)指標(biāo)與時間的關(guān)系，而油田開發(fā)指標(biāo)的預(yù)測，都離不開與時間的關(guān)系。因此如果將兩者結(jié)合起來進(jìn)行預(yù)測，不僅可以完善各種的理論，而且更有利于預(yù)測出更多的油田開發(fā)的相關(guān)參數(shù)。因此我們應(yīng)該用威布爾模型和乙型水驅(qū)曲線聯(lián)解來預(yù)測油田開發(fā)指標(biāo)。 結(jié)論 產(chǎn)量上升階段： 產(chǎn)量上升階段儲采比和產(chǎn)量呈現(xiàn)雙對數(shù)的直線關(guān)系。 穩(wěn)產(chǎn)階段： 穩(wěn)產(chǎn)階段儲采比與時間呈直線遞減關(guān)系。 遞減階段： 遞 減階段儲采比和產(chǎn)量呈現(xiàn)雙對數(shù)的直線關(guān)系。 全程： 全程儲采比和時間呈雙對數(shù)的直線關(guān)系。 威布爾模型預(yù)測可采儲量 RN 值： 由威布爾模型預(yù)測可采儲量 RN 值： 4R (10 t )N ? ，實際給的已知4R (10 t )N ? ，相對誤差為 %，這說明威布爾模型預(yù)測可采儲量有效。 東北石油大學(xué)課程設(shè)計（報告） 21 參考文獻(xiàn) 1. 陳元千，雙曲線遞減的簡化及確定 可采儲量的截距法，天然氣工業(yè) [J， （ 4）， 3237。 2. 陳元千，趙慶飛，油氣田儲采比變化關(guān)系的研究，斷塊油氣田 [J]， （ 6），2326； 3. 陳元千，胡建國，預(yù)測油氣田產(chǎn)量和可采儲量 weibull 模型，新疆石油地質(zhì)，（ 3）， 250255； 4. 陳元千，油田可采儲量計算方法，新疆石油地質(zhì) [J]， （ 2）， 130137； 5. 陳元千，油氣藏工程實踐，石油工業(yè)出版社 [M],。 東北石油大學(xué)課程設(shè)計（報告） 22 附錄（計算機(jī)程序） 一、 程序運行界面 數(shù)據(jù)導(dǎo)入及計算結(jié)果，如附圖 1 所示。 附圖 1 數(shù)據(jù)導(dǎo)入及計算結(jié)果 繪制產(chǎn)量、累積產(chǎn)量和儲采比隨時間的變化曲線程序界面，如附圖 2 所示 .。 東北石油大學(xué)課程設(shè)計（報告） 23 附圖 2 繪制產(chǎn)量、累積產(chǎn)量和儲采比隨時間的變化曲線程序界面 產(chǎn)量上升階段的計算結(jié)果和圖形輸出，如附圖 3 所示。 附圖 3 產(chǎn)量上升階段的計算結(jié)果和圖形輸出 東北石油大學(xué)課程設(shè)計（報告） 24 穩(wěn)產(chǎn)階段的計算結(jié)果和圖形輸出，如附圖 4 所示。 附圖 4 穩(wěn)產(chǎn)階段的計算結(jié)果和圖形輸出 遞減階段的計算結(jié)果和 圖形輸出，如附圖 5 所示。 東北石油大學(xué)課程設(shè)計（報告） 25 附圖 5 遞減階段的計算結(jié)果和圖形輸出 全程的計算結(jié)果和圖形輸出，如附圖 6 所示 附圖 6 全程的計算結(jié)果和圖形輸出 利用威布爾模型求可采儲量 RN 的計算結(jié)果和圖形輸出，如附圖 7 所示 東北石油大學(xué)課程設(shè)計（報告） 26 附圖 7 利用威布爾模型求可采儲量 RN 的計算結(jié)果和圖形輸出 二、計算機(jī)程序代碼 Form1 代碼： Dim t() As Single, nian() As Single, Q() As Single, NP() As Single, w() As Single Dim A As Single, B As Single, QS As Single, QAa As Single Dim xx As String, yy As String Dim XMAX As Single, YMAX As Single Dim XMIN As Single, YMIN As Single Dim x() As Single, y() As Single Sub ercheng(x() As Single, y() As Single, m As Integer, n As Integer) 39。定義最小二乘法過程進(jìn)行線性擬合 Dim i As Integer, R As Single Dim ER() As Single ReDim ER(5, m 1 To n) XMAX = x(m): YMAX = y(m) 39。求出橫坐標(biāo)和縱坐標(biāo)的最大，最小值。畫圖使用 XMIN = x(m): YMIN = y(m) 東北石油大學(xué)課程設(shè)計（報告） 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。 (。 i。 ) =。 Format(ER(1, i), ), yy。 (。 i。 ) = 。 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è)計（報告） 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。輸出擬合公式和相關(guān)性系數(shù) If A = 1 Then yy。 =。 xx Else yy。 =。 Format(A, )。 xx。 End If If B 0 Then +。 Format(B, ) ElseIf B 0 Then Format(B, ) End If 線性相關(guān)系數(shù) R=。 Format(Abs(R), ) Else 擬合公式為： If A = 1 Then yy。 =。 xx Else yy。 =。 Format(A, )。 xx。 End If If B 0 Then 東北石油大學(xué)課程設(shè)計（報告） 29 +。 Format(B, ) ElseIf B 0 Then Format(B, ) End If 線性相關(guān)系數(shù) R=。 F