注水開發儲采比研究報告——東北石油大學最終版(編輯修改稿)

　

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