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