油藏數(shù)值模擬原理黑油模型(參考版)
2025-03-25 07:41本頁面
【正文】 79 二 、 動態(tài)模型數(shù)據(jù)主要有： 1. 模型控制數(shù)據(jù)： （ 1） 離散方程求解控制參數(shù) （ 2） 線性代數(shù)方程組解法控制參數(shù) （ 3） 時間步長控制參數(shù) 2. 輸出控制參數(shù) : （ 1） 控制井資料報告 （ 2） 控制全油田 、 分區(qū)及分層報告 （ 3） 控制數(shù)據(jù)場報告：壓力 、 飽和度 80 3. 井?dāng)?shù)據(jù)： （ 1） 井位數(shù)據(jù)及射孔數(shù)據(jù) （ 2） 井的分類及生產(chǎn)方式 81 。 77 3． 網(wǎng)格屬性數(shù)據(jù)是用數(shù)組表述的屬性數(shù)據(jù) ,主要來自地質(zhì)建模 ,它們是： （ 1） 油層深度數(shù)據(jù) （ 2） 孔隙度數(shù)據(jù) （ 3） 滲透率數(shù)據(jù) （ 4） 厚度數(shù)據(jù) （ 5） 有效厚度 78 初始化程序輸入以下幾類卡片： （ 1） NIT卡 （ 2） 一般描述卡片 （ 3） PVT數(shù)據(jù)卡片 （ 4） 巖石類型數(shù)據(jù)卡片 （ 5） 數(shù)組數(shù)據(jù) （ 6） END卡 每一類型卡片有引領(lǐng)關(guān)鍵字 ， 標志者該類數(shù)據(jù)開始輸入。 （ 2） 油水系統(tǒng)以及油氣系統(tǒng)的相對滲透率曲線和毛管壓力曲線數(shù)據(jù) 。 76 2． 巖石數(shù)據(jù)包括油水系統(tǒng)和油氣系統(tǒng)的相對滲透率曲線和毛管壓力曲線數(shù)據(jù) ， 以及其它的巖石性質(zhì)數(shù)據(jù) 。 （ 2） 平衡區(qū)數(shù)據(jù) ， 如油水界面 、 油氣界面 。 8. 線性代數(shù)方程組的求解 72 黑油模型油藏數(shù)值模擬器流程 二、黑油模型簡介 73 初始程序 INIT文件 輸出結(jié)果 動態(tài)數(shù)據(jù) 動態(tài)程序 主模型程序 輸出結(jié)果 輸出 重啟文件 輸入重啟文件 拷貝文件 靜態(tài)數(shù)據(jù) 74 黑油模型輸入數(shù)據(jù)： 主要數(shù)據(jù)流 ? ?? ?????????????????? ??( P V T )(油藏流體性質(zhì) 數(shù)據(jù)初始化數(shù)據(jù) 油藏描述 巖性）數(shù)據(jù)網(wǎng)格數(shù)據(jù)模型控制數(shù)據(jù)主模型數(shù)據(jù) 輸出控制數(shù)據(jù)時間相關(guān)（井）數(shù)據(jù) 75 一 、 初始化數(shù)據(jù)主要有： 1． PVT數(shù)據(jù)包括流體性質(zhì)常數(shù) ， 高壓物性試驗數(shù)據(jù)和平衡區(qū)數(shù)據(jù) 。 BO軟件中采用 RSVP方法 。 因此 ， 研制快速求解大型稀疏矩陣的新方法成為油藏數(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? ? ?? ? ? ? ? ? ?? ? ? ? ? ?? ? ? 69 全隱式方法形成的線性代數(shù)方程組的結(jié)構(gòu) 7. 差分方程的線性化 x b A 矩陣結(jié)構(gòu)示意圖 70 油藏數(shù)值模擬的線性代數(shù)方程組是大型稀疏矩陣， 其求解效率在油藏數(shù)值模擬計算中舉足輕重 。 變量進行運算后也保留線性部分 ， 舍去高階項 。 其中： x代表求解變量 n—— 時間迭代步 —— 牛頓迭代步 6. 差分方程 1nnx x x? ???1llx x x? ???11n l lx x x x???? ? ?l 67 求解變量： ， ， 對于二相狀態(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??? 68 對方程組進行線性化處理 。 用一維非線性方程表示： 給定初始近似 對函數(shù) 在點處 進行一階泰勒展開 , 可近 似地表示為 6. 差分方程 0 0 0f ( ) f ( ) f ( ) ( ) 0x x x x x?? ? ? ?f( ) 0x ? 0xf( ) 0x ? 0x 65 半隱式： 全隱式方法： 6. 差分方程 10()()nnnnfxxxfxx?? ??? ???? 為初值000()()fxxxfx?? ? 66 下面介紹全隱式方法。 6. 差分方程 64 全隱式方法和半隱式方法的根本差別在于： 半隱式方法 是將非線性方程直接線性化形成 線性代數(shù)方程組，來求解 n+1時刻的值。 全隱式方法 主要采用了 NewtonRaphson迭代技術(shù) , 對方程中各項系數(shù)，包括井系數(shù)全部隱式處理。因此時間步長不能太大。對流度中與時間有關(guān)的量均取上一時間步的值。 根據(jù)隱式程度的不同，數(shù)值模擬的算法有IMPES (Implicit Pressure Explicit Saturation) 方法、 半隱式 (SemiImplicit Method)方法 和 全隱式方法 (Fully Implicit Method )。非線性最強的參數(shù)是相對滲透率和毛管力。該方法在研究之中，最早是 Heinemann將動態(tài)局部網(wǎng)格用于水驅(qū)油藏數(shù)值模擬 。該技術(shù)已經(jīng)成熟，在多加軟件公司的軟件產(chǎn)品中都有這一功能。 能夠正確的模擬沿著斷層面層之間的接觸關(guān)系 容易建立徑向網(wǎng)格模型 若沒有前處理軟件，則很難建立徑向網(wǎng)格模型。 ZCORN 這個關(guān)鍵字定義網(wǎng)格拐點的深度，數(shù)據(jù)按照 X (或 R) 方向循環(huán)最快，然后是 Y (或 THETA) 方向，最后是 Z方向的順序排列。 ? 層厚度從上至下分別是 32， 22， 20， 4， 32， 4， 26， 26， 4， 28ft 38 :塊中心網(wǎng)格的數(shù)據(jù)結(jié)構(gòu) 前 20個 TOPS定義網(wǎng)格 (1, 1, 1) 到 (20, 1, 1) 下面 20個 TOPS定義網(wǎng)格 (1, 2, 1) 到 (20, 2, 1) 下面 20個
教學(xué)課件相關(guān)推薦
鄂ICP備17016276號-1