【正文】 ............................................................................. 1 2 基礎(chǔ)數(shù)據(jù) .................................................................................................................................. 2 3 基礎(chǔ)理論 .................................................................................................................................. 3 威布爾預(yù)測(cè)模型的建立 ................................................................................................... 3 威布爾模型與乙型水驅(qū)曲線(xiàn)的聯(lián)解法 ........................................................................... 5 模型的求解方法 ............................................................................................................... 6 4 設(shè)計(jì)結(jié)果 .................................................................................................................................... 8 認(rèn)識(shí)與結(jié)論 .................................................................................................................................. 12 參考文獻(xiàn) ...................................................................................................................................... 12 附錄：計(jì)算機(jī)程序與結(jié)果 .......................................................................................................... 14 1 1 前 言 設(shè)計(jì)的目的意義 本課題的目的是讓學(xué)生通過(guò)自選一組數(shù)據(jù)，利用所學(xué)專(zhuān)業(yè)知識(shí)在指導(dǎo)教師的指導(dǎo)下獨(dú)立完成對(duì)某一油田或區(qū)塊開(kāi)發(fā)指標(biāo)的預(yù)測(cè)。 2 2 基礎(chǔ)數(shù)據(jù) 某油田的開(kāi)發(fā)數(shù)據(jù)如下： 表 21 某油田開(kāi)發(fā)數(shù)據(jù) 年份 時(shí)間（ a） Qo（ 104t/a） Qw（ 104t/a） Np（ 104t） Wp（ 104t） Lp（ 104t） 1968 1 1969 2 1970 3 1971 4 1972 5 1973 6 1974 7 1975 8 1976 9 1977 10 1978 11 1979 12 1980 13 1981 14 1982 15 1983 16 1984 17 1985 18 1986 19 1987 20 1988 21 1989 22 1990 23 3 3 基礎(chǔ)理論 威布爾（ Weibull）預(yù)測(cè)模型能夠很好地預(yù)測(cè)油田產(chǎn)油量隨時(shí)間的變化關(guān)系，但卻不能預(yù)測(cè)油田的含水率、產(chǎn)水量、產(chǎn)液量及累積產(chǎn)水量和累積產(chǎn)液量，而這些開(kāi)發(fā)指標(biāo)正是水驅(qū)開(kāi)發(fā)油田所需要預(yù)測(cè)的。 威布爾預(yù)測(cè)模型的建立 Weibull(威布爾 )于 1939 年提出的統(tǒng)計(jì)分布模型已成為生命試驗(yàn)和可靠性理論研究的基礎(chǔ)。 式中 NP—— 油氣田的累積產(chǎn)量， 104t 或 108t(油 )； 108m(氣 )。因此，可以將（ 33）式再改寫(xiě)為： 1 ( / )R e tNQt ????? ??? （ 37） 為了確定最高年產(chǎn)量發(fā)生的時(shí)間，有（ 37）式對(duì)時(shí)間 t 求導(dǎo)數(shù)得： ? ?2 ( / )Rd 1ed tNQ ttt ?? ? ?? ????????? ? ????? （ 38） 當(dāng) ddQt =0 時(shí)，必然有 ? ?1 t??????=0，故可以得到最高年產(chǎn)量發(fā)生的時(shí)間 tm為： ? ? 1/1mt ???????? ???? （ 39） 將（ 39）式代入（ 37）式，得到油氣田的最高年產(chǎn)量（ Qmax）的表達(dá)式： ? ? ? ?1/ 1 1 / 1/m a x R 1eQN ? ? ??? ?? ? ?????????????? （ 310） 再將（ 39）式代入（ 36）式，得到油氣田最高年產(chǎn)量發(fā)生時(shí)的累積產(chǎn)量（ Npm）為： ? ?? ?1/p m R 1eNN ?????????? （ 311） 油氣田的剩余可采儲(chǔ)量 NRR 表 示為： RR R PN N N?? （ 312） 將（ 36）式代入（ 312）式得： ? ?/RR R e tNN ? ??? （ 313） 5 剩余可采儲(chǔ)量的儲(chǔ)采比 ω 表示為： RR /NQ?? （ 314） 將（ 37）式和（ 313）式代入（ 314）式得： 1/ ( )t?? ? ? ?? （ 315） 剩余可采儲(chǔ)量的采油速度為儲(chǔ)采比的倒數(shù)，故由（ 315）式得到剩余可采儲(chǔ)量采油速度 o? 的表達(dá)式： 1o ( ) /t?? ? ??? （ 316） 式中 o? 以小數(shù) f 表示，若改以百分?jǐn)?shù) %表示時(shí)，（ 316）式改為下式： 1o (100 / )%t?? ? ??? （ 317） 威布爾模型與乙型水驅(qū)曲線(xiàn)的聯(lián)解法 利用數(shù)理統(tǒng)計(jì)學(xué)中的威布爾（ Weibull）分布， 研究與推導(dǎo)得到了 威布爾（ Weibull）預(yù)測(cè)模型 [1]。 在確定預(yù)測(cè)模型常數(shù) a、 b、 c、 NR時(shí)，其值是否正確可靠，要利用（ 318）式、（ 319）式、（ 328）式預(yù)測(cè)的理論產(chǎn)油量、累積產(chǎn)油量和含水率，與實(shí)際產(chǎn)油量、累積產(chǎn)油量和含水率進(jìn)行對(duì)比加以確定，而達(dá)到最佳擬合效果的參數(shù)才是最準(zhǔn)確、可靠的。經(jīng)線(xiàn)性回歸求得直線(xiàn)的截距 ? =；直線(xiàn)的斜率 ? =；參數(shù) a=， b=， c=， 直線(xiàn)的相關(guān)系數(shù) r=。 圖 43 產(chǎn)量與時(shí)間關(guān)系曲線(xiàn) 11 按照（ 319）式的關(guān)系在同一坐標(biāo)系中繪制實(shí)際累計(jì)產(chǎn)量與預(yù)測(cè)累 計(jì)產(chǎn)量及其對(duì)應(yīng)時(shí)間 t 之間的關(guān)系曲線(xiàn)，繪于圖 44。雖然在課程設(shè)計(jì)的過(guò)程中遇到了一些問(wèn)題，但經(jīng)過(guò)反復(fù)的思考與檢查找出原因所在并進(jìn)行改正。 （ 5） 這次課程設(shè)計(jì)，使我認(rèn)識(shí)到學(xué)好計(jì)算機(jī)的重要性與實(shí)踐操作的必要性， 在老師的指導(dǎo)下和與 同學(xué)的討論中，不斷地發(fā)現(xiàn)問(wèn)題、不斷地解決問(wèn)題、不斷地領(lǐng)悟、不斷地獲取，使自己的程序不斷完善。求相關(guān)系數(shù) r,趨勢(shì)線(xiàn) sx = 0: sy = 0: sxy = 0: sx2 = 0: sy2 = 0 l = ub an For i = l To an + l sx = sx + np(i) Next i 15 For i = l To an + l sy = sy + lp(i) Next i For i = l To an + l sxy = sxy + np(i) * lp(i) Next i For i = l To an + l sx2 = sx2 + np(i) ^ 2 Next i For i = l To an + l sy2 = sy2 + lp(i) ^ 2 Next i a2 = ((an + 1) * sxy sx * sy) / ((an + 1) * sx2 sx ^ 2) a1 = (sx2 * sy sx * sxy) / ((an + 1) * sx2 sx ^ 2) r = ((an + 1) * sxy sx * sy) / (((an + 1) * sx2 sx ^ 2) ^ * ((an + 1) * sy2 sy ^ 2) ^ ) 39。最大相關(guān)系數(shù) Rem weibull 參數(shù)計(jì)算 bn = 12 l = ub bn bmax = 0 rmax = 0 16 For b = 1 To 2 Step sx = 0: sy = 0: sxy = 0: sx2 = 0: sy2 = 0 For i = l To l + bn t(i) = i + 1 Next i For i = l To l + bn qt(i) = Log(q0(i) / (t(i) ^ b)) / Log(10)