【正文】 (311) 其中： RlgAN? 。 穩(wěn)產(chǎn)階段儲(chǔ)采比變化特征 在穩(wěn)產(chǎn)階段某一開(kāi)發(fā)時(shí)間的儲(chǔ)采比表示為： SR P I PSON N NQ? ??? (312) 在穩(wěn)產(chǎn)階段某一開(kāi)發(fā)時(shí)間的累積產(chǎn)量可表示為： SSP O 1()N Q t t?? (313) 將（ 313）式代入（ 312）式可得，在穩(wěn)產(chǎn)階段儲(chǔ)采比與開(kāi)發(fā)時(shí)間 的關(guān)系式為： SR P I O 1 R P I1SSOO()N N Q t t NN tt? ? ? ? ?? ? ? ? (314) 則（ 314）式可寫成： Et??? (315) 其中： R PI 1SONNEtQ???。 東北石油大學(xué)課程設(shè)計(jì)（報(bào)告） 4 遞減階段儲(chǔ)采比變化關(guān)系 進(jìn)入遞減階段之后的儲(chǔ)采比可表示為： SdR P I P PON N N NQ? ? ? ?? (316) 進(jìn)入遞減階段累積產(chǎn)量與產(chǎn)量之 間的關(guān)系如下： S S nd 1 nOOPOii()(1 ) (1 )NQD n D n? ? ??? (317) 則（ 317）式可以簡(jiǎn)化為： d 1nPON a bQ?? (318) 其中：S S nOOii(),(1 ) (1 )abD n D n????。 將（ 322）式代入（ 320）式可得： S 1 nR P I P OO()N N N a b ? ? ? ? ?? (319) 由于（ 322）式可知： PmaxNa? (320) 注水開(kāi)發(fā)油田的可采儲(chǔ)量可表示為： R P I P S P m a xN N N N? ? ? (321) 由（ 320）式和（ 321）式可得： P m a x R P I P SN N N N a? ? ? ? (322) 將（ 322）式代入（ 323）式可得： nObQ?? (323) 將（ 323）式兩邊取常用對(duì)數(shù)得： Olg lg lgb n Q? ?? (324) 則（ 324）式可改寫為： Olg lg Q? ? ??? (325) 其中： lg ,bn????。 東北石油大學(xué)課程設(shè)計(jì)（報(bào)告） 5 全程儲(chǔ)采比變化關(guān)系 假設(shè)全程產(chǎn)量、累積產(chǎn)量、和可采儲(chǔ)量隨時(shí)間的關(guān)系可由威布爾預(yù)測(cè)模型表示： b+1b (t /c)Q at e? (326) b+ 1t( )cPR1N N e???????? (327) R 1acN b? ? (328) 將（ 326）、（ 327）、（ 328）式代入儲(chǔ)采比的定義式中得到下式： R P RbN N NQ at? ??? (329) 將（ 329）式兩邊取常用對(duì)數(shù)得： Rlg lg lgN bta? ???????? (330) 則（ 330）式可簡(jiǎn)化為： lg lgA B t? ?? (331) 其中： Rlg ,NA B ba????????。 確定 ,ab和 c 常數(shù)方法 由產(chǎn)量隨時(shí)間的關(guān)系可由威布爾預(yù)測(cè)模型表示： b+1b (t /c)Q at e? (332) 將（ 332）式改寫成： b+1(t /c)bQ aet ? (333) 將（ 333）是等號(hào)兩端取常用對(duì)數(shù)得： b + 1b 1lg lg 2 .3 0 3Q attC?? (334) 則（ 334）式可簡(jiǎn)化為： b+1blg Q tt ???? (335) 其中： 1lg , 2 .3 0 3a c???? 東北石油大學(xué)課程設(shè)計(jì)（報(bào)告） 6 那么： α10a? (336) ?? (337) 根據(jù)油田開(kāi)發(fā)的實(shí)際產(chǎn)量 Q 和相應(yīng)的時(shí)間 t 的數(shù)據(jù)，給出不同的 b 值，利用（ 335）式進(jìn)行線性迭代試差，能得到相關(guān)系數(shù)最高的 b 值，即為所求的 b 值。由線性回歸得到直線的斜率 ? 和截距 ? 值后，再由（ 336）式和（ 337）式求出a 和 c 的數(shù)值。 東北石油大學(xué)課程設(shè)計(jì)（報(bào)告） 7 目錄 4 程序設(shè)計(jì)框圖 計(jì)算機(jī)程序 設(shè)計(jì)框圖 最小二乘法程序 設(shè)計(jì)框圖 ，見(jiàn)圖 41。 輸入 x和 y 值 分別計(jì)算 2 2 2 2, , , , , , ,x y x y x y x y x y? ? ? ? ? 計(jì)算 ,xy 計(jì)算直線的斜率 22( ) / ( )A x y n x y x n x? ? ??? 計(jì)算直線的截距 B y Ax?? 計(jì)算線性相關(guān)系數(shù)2222( ) ( )( ) ( )xyxynRxynn??????????? 輸出擬合公式 y Ax B?? 開(kāi)始 結(jié)束 東北石油大學(xué)課程設(shè)計(jì)（報(bào)告） 8 圖 41 最小二乘法程序 設(shè)計(jì)框圖 產(chǎn)量上升階段程序設(shè)計(jì)框圖，如圖 42 所示 圖 42 產(chǎn)量上升階段程序設(shè)計(jì)框圖 穩(wěn)產(chǎn)階段程序設(shè)計(jì)框圖，如圖 43 所示 開(kāi)始 輸入產(chǎn)量上升階段起始點(diǎn) 計(jì)算 lg ,lgQ? 調(diào)用最先二乘法進(jìn)行直線 lg lgAQ? ?? 擬合 繪出散點(diǎn)圖和擬合直線 圖 利用直線截距求出可采儲(chǔ)量 RN 結(jié)束 開(kāi)始 輸入穩(wěn)產(chǎn)階段起始點(diǎn) 調(diào)用最先二乘法進(jìn)行直線 Et???擬合 繪出散點(diǎn)圖和擬合直線圖 結(jié)束 東北石油大學(xué)課程設(shè)計(jì)（報(bào)告） 9 圖 43 穩(wěn)產(chǎn)階段程序設(shè)計(jì)框圖 遞減階段程序設(shè)計(jì)框圖，如圖 44 所示。 圖 44 遞減階段程序設(shè)計(jì)框圖 開(kāi)始 輸入產(chǎn)量遞減階段起始點(diǎn) 計(jì)算 lg ,lgQ? 調(diào)用最先二乘法進(jìn)行直線 lg lg Q? ? ??? 擬合 繪出散點(diǎn)圖和擬合直線圖 利用直線斜率求出遞減指數(shù) n 結(jié)束 東北石油大學(xué)課程設(shè)計(jì)（報(bào)告） 10 全程程序設(shè)計(jì)框圖，如圖 45 所示。 圖 45 全程程序設(shè)計(jì)框圖 計(jì)算機(jī)程序 開(kāi)始 輸入全程起始點(diǎn) 計(jì)算 lg ,lgt? 調(diào)用最先二乘法進(jìn)行直線 lg lgA B Q? ?? 擬合 繪出散點(diǎn)圖和擬合直線圖 利用直線斜率求出 b 值 輸入 b 值 調(diào)用最先二乘法進(jìn)行直線 b+1blg Q tt ????擬合 繪出散點(diǎn)圖和擬合直線圖 利用直線截距與斜率求出 ,ab值 結(jié)束 相關(guān)度好不好 否 是 東北石油大學(xué)課程設(shè)計(jì)（報(bào)告） 11 計(jì)算機(jī)程序代碼，見(jiàn)附錄。 5 設(shè)計(jì)結(jié)果及分析 利用已知數(shù)據(jù)求得累積產(chǎn)量和儲(chǔ)采比，見(jiàn)表 51。 表 51利用已知數(shù)據(jù)求得累積產(chǎn)量和儲(chǔ)采比 可采儲(chǔ)量（萬(wàn)噸） 序號(hào) 年 年產(chǎn)量（萬(wàn)噸） Np（萬(wàn)噸） w 1 1953 2 1954 3 1955 1113 4 1956 5 1957 6 1958 7 1959 8 1960 9 1961 3545 10 1962 11 1963 12 1964 13 1965 14 1966 15 1967 16 1968 17 1969 7195 18 1970 19 1971 20 1972 21 1973 22 1974 23 1975 24 1976 東北石油大學(xué)課程設(shè)計(jì)（報(bào)告） 12 25 1977 26 1978 27 1979 6191 28 1980 29 1981 30 1982 31 1983 32 1984 33 1985 34 1986 35 1987 36 1988 37 1989 38 1990 39 1991 174907 40 1992 41 1993 42 1994 43 1995 44 1996 45 1997 46 1998 47 1999 48 20xx 190 49 20xx 東北石油大學(xué)課程設(shè)計(jì)（報(bào)告） 13 繪制產(chǎn)量、累積產(chǎn)量和儲(chǔ)采比隨時(shí)間的變化曲線： 圖 51產(chǎn)量與時(shí)間關(guān)系曲線 東北石油大學(xué)課程設(shè)計(jì)（報(bào)告） 14 圖 52 累積產(chǎn)量和時(shí)間的關(guān)系曲線 圖 53 儲(chǔ)采比和時(shí)間的關(guān)系曲線 ，如圖 54所示。 y = 91 x + 51 8R178。 = 00 00 0020 40 60 80 10 012 014 016 00 .00 018 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ǔ)采比產(chǎn)量產(chǎn)量上升階段儲(chǔ)采比與產(chǎn)量的變化曲線系列 1 線性 ( 系列 1)圖 54 產(chǎn)量上升階 段儲(chǔ)采比與產(chǎn)量的關(guān)系曲線 東北石油大學(xué)課程設(shè)計(jì)（報(bào)告） 15 穩(wěn)產(chǎn)階段儲(chǔ)采比與開(kāi)發(fā)時(shí)間的關(guān)系曲線，如圖 55所示。 00 00 00 00 0010 .00 012 .00014 .00016 .00018 19 20 21 22儲(chǔ)采比時(shí)間系列 1圖 55 穩(wěn)產(chǎn)階段儲(chǔ)采比與開(kāi)發(fā)時(shí)間的關(guān)系曲線 遞減階段儲(chǔ)采比與產(chǎn)量的關(guān)系曲線，如圖 56所示。 y = R178。 = 3 2 1 01012345LogWLogQ遞減階段儲(chǔ)采比與產(chǎn)量的關(guān)系曲線系列 1線性 ( 系列 1) 圖 56 遞減階段儲(chǔ)采比與產(chǎn)量的關(guān)系曲線 ?東北石油大學(xué)課程設(shè)計(jì)（報(bào)告） 16 全程儲(chǔ)采比與開(kāi)發(fā)時(shí)間的關(guān)系曲線，如圖 57所示。 圖 57 全程儲(chǔ)采比與開(kāi)發(fā)時(shí)間的關(guān)系曲線 求 威布爾模型中的 ,ab和 c 的值關(guān)系曲線 ，如圖 58 所示。 圖 58 求威布爾模型中的 ,ab和 c 的值關(guān)系曲線 b+1t東北石油大學(xué)課程設(shè)計(jì)（報(bào)告） 17 產(chǎn)量上升階段儲(chǔ)采比與產(chǎn)量的變化曲線，如圖 59所示。 0 . 0 0 02 0 0 . 0 0 04 0 0 . 0 0 06 0 0 . 0 0 0