數(shù)值模擬原理黑油模型(參考版)
2025-02-24 10:26本頁面
【正文】 INIT卡必須放在卡片疊的第一張 ， END放在卡片疊的最后一張 ， 數(shù)組數(shù)據(jù)放在其它輸入數(shù)據(jù)的后面 。 （ 3） 毛管壓力端點值與束縛水飽和度關(guān)系 。 主要是： （ 1） 輸入巖石壓縮系數(shù)及其參考壓力 。 （ 3） 油 PVT表 ， 氣 PVT表 。 主要是： （ 1） 油 、 氣 、 水的地面密度和比重等常數(shù) ， 以及計算地層水性質(zhì)的有關(guān)參數(shù) 。 該方法計算快 ， 精度高 。 因此 ， 研制快速求解大型稀疏矩陣的新方法成為油藏數(shù)值模擬發(fā)展的一個重要方向 。 7. 差分方程的線性化 1ll T T TT T P S w S xP S w S x? ? ?? ? ? ?? ? ? ?? ? ?1llw P S w S xP S w S x? ? ?? ? ? ? ? ? ?? ? ? ? ? ?? ? ? 45 全隱式方法形成的線性代數(shù)方程組的結(jié)構(gòu) 7. 差分方程的線性化 x b A 矩陣結(jié)構(gòu)示意圖 46 油藏數(shù)值模擬的線性代數(shù)方程組是大型稀疏矩陣， 其求解效率在油藏數(shù)值模擬計算中舉足輕重 。 變量進(jìn)行運算后也保留線性部分 ， 舍去高階項 。 其中： x代表求解變量 n—— 時間迭代步 —— 牛頓迭代步 6. 差分方程 1nnx x x? ???1llx x x? ???11n l lx x x x???? ? ?l 43 求解變量： ， ， 對于二相狀態(tài)： 對于三相狀態(tài)： 6. 差分方程 1 1 11111 1 1 1111( ) [ ( ) ( ) ( ) ]( ) [ ( ) ( ) ( ) ]( ) ( )[ ( ) ( ) (l l l lnw w w w w w w w wlol ol o l o no o o o o o o o ol l l lgl glo o g g o gg l g n go o g g o o g g oVT q S S StVT q S S StT T q qVS S S St? ? ? ? ? ? ?? ? ? ? ? ? ?? ? ? ? ? ? ? ? ? ? ?? ? ????? ? ? ????? ? ? ? ? ? ? ??? ? ? ? ? ? ? ??? ? ? ? ? ? ? ? ? ? ?? ? ? ? ??)]o g gSS ??????????????P? Sw? Sx?S x S o???S x P s??? 44 對方程組進(jìn)行線性化處理 。 用一維非線性方程表示： 給定初始近似 對函數(shù) 在點處 進(jìn)行一階泰勒展開 , 可近 似地表示為 6. 差分方程 0 0 0f ( ) f ( ) f ( ) ( ) 0x x x x x?? ? ? ?f( ) 0x ? 0xf( ) 0x ? 0x 41 半隱式： 全隱式方法： 6. 差分方程 10()()nnnnfxxxfxx?? ??? ???? 為初值000()()fxxxfx?? ? 42 下面介紹全隱式方法。 6. 差分方程 40 全隱式方法和半隱式方法的根本差別在于： 半隱式方法 是將非線性方程直接線性化形成 線性代數(shù)方程組，來求解 n+1時刻的值。 全隱式方法 主要采用了 NewtonRaphson迭代技術(shù) , 對方程中各項系數(shù)，包括井系數(shù)全部隱式處理。因此時間步長不能太大。對流度中與時間有關(guān)的量均取上一時間步的值。 根據(jù)隱式程度的不同，數(shù)值模擬的算法有IMPES (Implicit Pressure Explicit Saturation) 方法、 半隱式 (SemiImplicit Method)方法 和 全隱式方法 (Fully Implicit Method )。非線性最強(qiáng)的參數(shù)是相對滲透率和毛管力。 能夠正確的模擬沿著斷層面層之間的接觸關(guān)系 容易建立徑向網(wǎng)格模型 若沒有前處理軟件，則很難建立徑向網(wǎng)格模型。 有限差分方法大多空間變量采用 網(wǎng)格塊中心 （ Block Center)差分 , 時間變量按 向前差分方法 。 離散化方法主要有：有限差分 （ Finite Difference)、 有限元 (Finite Element)以及流線法 (Stream Line)。 網(wǎng)格剖分方法有很多 ， 像直角坐標(biāo)網(wǎng)格 、 徑向網(wǎng)格 、 角點網(wǎng)格 、 PEBI網(wǎng)格等 。 △ X，△ Y，△ Z—— 網(wǎng)格節(jié)點 x， y和 z方向步長 PI叫采油指數(shù)， 通常寫成 KH,可以直接輸入 2ln ewKZPIrr? ??erwr0 . 2er x y? ? ??BHP BHzKZ?井?dāng)?shù)據(jù)又叫動態(tài)數(shù)據(jù) 26 4 .初邊值條件 ? 外邊界 封閉邊界 :尖滅、斷層、圈閉 流動邊界：邊、底水 ? 內(nèi)邊界，即 井的生產(chǎn)條件 定產(chǎn)邊界條件（ rate cons
教學(xué)課件相關(guān)推薦
鄂ICP備17016276號-1