CMG中IMEX中文操作手冊(cè)(一)

1、 目 錄 imex 的介紹1指導(dǎo)段5關(guān)鍵字輸入系統(tǒng)中的數(shù)據(jù)段5如何建立數(shù)據(jù)文件的文檔方式6如何執(zhí)行重新啟動(dòng)運(yùn)行7控制打印輸出文件的內(nèi)容8圖形文件(sr2)的內(nèi)容控制9網(wǎng)格系統(tǒng)描述10局部加密網(wǎng)格描述11雙孔/雙滲模型的使用12死結(jié)點(diǎn)的輸入17水區(qū)選項(xiàng)的使用18擬混相選項(xiàng)的使用19注聚合物模型的使用20程序執(zhí)行時(shí)間長或時(shí)間步長太小問題21單相油藏的模擬23水平井25垂直平衡計(jì)算25定義多個(gè) pvt 區(qū)29井的定義30井的類型定義31如何關(guān)井以及重新開井32虧空填充(voidage replacement)34井筒流動(dòng)模型的使用36操作及監(jiān)測(cè)限制37井指數(shù)的輸入39中止模擬運(yùn)行41在井的列表中使用

2、通配符41關(guān)于井?dāng)?shù)據(jù)段設(shè)置的指導(dǎo)42在循環(huán)數(shù)據(jù)段內(nèi)可用的其他段關(guān)鍵字44關(guān)鍵字?jǐn)?shù)據(jù)輸入系統(tǒng)46關(guān)鍵字系統(tǒng)介紹46注釋行 (可選擇)52空行52包含文件53控制數(shù)據(jù)文件列表54基巖網(wǎng)格性質(zhì)的輸入57裂縫網(wǎng)格性質(zhì)的輸入57加密網(wǎng)格性質(zhì)的輸入58由 i 方向確定 j 和 k 方向的數(shù)據(jù)59常數(shù)值數(shù)組59以 ijk 方式輸入數(shù)組60數(shù)組輸入值沿 i 方向變化61數(shù)組輸入值沿 j 方向變化62數(shù)組輸入值沿 k 方向變化62大多數(shù)或所有網(wǎng)格的值都是不同的63輸入/輸出控制段67基本網(wǎng)格定義99k 坐標(biāo)方向105i 方向的網(wǎng)格尺寸107j 方向的網(wǎng)格尺寸108k 方向的網(wǎng)格尺寸110網(wǎng)格的頂部深度112加密

3、網(wǎng)格的位置123雙孔介質(zhì)minc方法126死結(jié)點(diǎn)標(biāo)識(shí)符128孔隙度130巖石壓縮系數(shù)131滲透率134壓力影響函數(shù)147流體組分性質(zhì)數(shù)據(jù)段153流體模型153油藏溫度154油和氣的 pvt 表155作為壓力函數(shù)的油壓縮系數(shù)168密度173水相粘度179原始?xì)庥捅?80聚合物粘度混合（條件）183巖石流體數(shù)據(jù)段186相對(duì)滲透率表186油水相對(duì)滲透率表187氣液相對(duì)滲透率表194評(píng)價(jià)三相共存時(shí)油相相對(duì)滲透率的方法（任選）203初始條件數(shù)據(jù)段209初始條件標(biāo)識(shí)209油藏初始油相壓力212初始泡點(diǎn)壓力212初始泡點(diǎn)壓力與深度關(guān)系213初始含油飽和度215初始含水飽和度216初始聚合物濃度216參考深度

4、和參考?jí)毫?17油水界面深度218數(shù)值計(jì)算方法數(shù)據(jù)控制段225井?dāng)?shù)據(jù)段241井的改變?nèi)掌?41井的改變時(shí)間242設(shè)定井底流壓初始化的頻率248附錄32065 imex 的介紹簡介 imex 是一個(gè)考慮重力及毛細(xì)管力的三相黑油模擬軟件,網(wǎng)絡(luò)系統(tǒng)可采用直角坐標(biāo),徑向坐標(biāo),變深度/變厚度坐標(biāo),在任何網(wǎng)絡(luò)系統(tǒng)中.都可建立兩維或三維模型.在處理氣相的出現(xiàn)及消失情況時(shí),程序采用了變量替換方法.imex 的一些特征和功能為：自適應(yīng)隱式方法 imex可以在顯示,全隱式以及自適應(yīng)隱式三種方式下運(yùn)行.在大多數(shù)情況下,只有很少一部分網(wǎng)格需要采用全隱式求解,而大部分網(wǎng)格都可采用顯式方法求解.自適應(yīng)隱式方法正是適合于這

5、種情況的解法,并且在井附近以及層狀油藏的薄層中,開采時(shí)會(huì)產(chǎn)生高速流動(dòng)的錐進(jìn)問題,采用自適應(yīng)隱式處理這類問題是很有效的. 采用自適應(yīng)隱式選項(xiàng)可節(jié)省三分之一到一半的運(yùn)行時(shí)間。計(jì)算時(shí)可采用和全隱式方法同樣大的時(shí)間步長.用戶可以指定采用全隱式方法計(jì)算的網(wǎng)格，可根據(jù)用戶確定的界限或矩陣轉(zhuǎn)換臨界值，動(dòng)態(tài)地選擇采用全隱式計(jì)算的網(wǎng)絡(luò)網(wǎng)格。雙孔/雙滲 雙孔隙度選項(xiàng)允許采用兩種方法對(duì)基巖模型進(jìn)行離散化處理，其中一種為嵌套格式，成為“多重內(nèi)部作用連續(xù)域”(minc)方法，另一種為層狀格式，稱作“子區(qū)域”方法。雙孔隙模型對(duì)裂縫油藏進(jìn)行了理想化的近似處理，認(rèn)為裂隙油藏由兩部分組成：主要孔隙度和次要孔隙度， 主要孔隙度(

6、基巖)代表巖塊中的微小粒間孔隙，次要孔隙度(裂縫)由裂縫，通道和溶洞組成。雙孔隙模型將油藏分為兩個(gè)連續(xù)域，裂隙是流體流動(dòng)的主要通道，只具有很小的儲(chǔ)集性能；而基巖具有較低的流體傳導(dǎo)能力，但具有較大的儲(chǔ)存能力。使用單個(gè)基巖巖塊/一個(gè)裂縫系統(tǒng)就可確定一個(gè)簡單的雙孔隙模型。形狀因子的計(jì)算基于 warren 和 root 或 gilman 和 kszwmi 的工作。 在這種情況下，假設(shè)基巖和裂縫是在半穩(wěn)定流狀態(tài)下進(jìn)行傳輸?shù)?。為了恰?dāng)?shù)乇硎净鶐r與裂縫間的傳導(dǎo)機(jī)理，有必要將基巖進(jìn)一步劃分為較小的體積單元。minc 法通過將基巖劃分成一系列嵌套的體積單元來達(dá)到這種目的，離散化時(shí)基于這樣一種假設(shè)，即體積單元表面

7、為等勢(shì)面，通過一個(gè)裂縫網(wǎng)格時(shí)，主要變量的變化很小。 minc 模型對(duì)于通過基巖一裂縫的非穩(wěn)定流動(dòng)是一種很好的表示方法，但是當(dāng)基巖的高度很大時(shí)，不能適當(dāng)?shù)奶幚碇亓τ绊?，是模型的一個(gè)嚴(yán)重缺陷。一種可以表示重力影響和流體相態(tài)分異現(xiàn)象的改進(jìn)模型是子區(qū)域方法，在這種模型中，對(duì)基巖中的流體和壓力分布進(jìn)行計(jì)算，但假設(shè)重力分異現(xiàn)象只存在于裂縫之中。 除了 minc 模型和子區(qū)域模型外，雙滲模型允許基巖之間的流體傳導(dǎo)，當(dāng)裂縫的連通具有方向性時(shí)，或者當(dāng)基巖之間的連通性顯的重要時(shí)，例如氣一油的重力驅(qū)替過程，采用雙滲模型進(jìn)行計(jì)算是十分有用的。擬混相驅(qū)選項(xiàng) 擬混相驅(qū)模型可模擬一次接觸及多次接觸混相過程。并能模擬水相中的

8、溶劑及驅(qū)動(dòng)氣。聚合物模型選項(xiàng) 聚合物模型選項(xiàng)可模擬巖石吸附，聚合物不可進(jìn)入孔隙體積，阻力系數(shù)，混合粘度的擴(kuò)散情況。斷層油藏選項(xiàng) 可模擬具有一條或多條斷層的層狀油藏，可精確地計(jì)算錯(cuò)斷層間的流體流動(dòng)過程，模擬執(zhí)行過程中不帶有任何不真實(shí)的性質(zhì)平均現(xiàn)象。不平行于坐標(biāo)軸的下滑斷層也能進(jìn)行模擬計(jì)算。井的全隱式處理 采用了一種非常有效的方式處理井的問題，井底壓力以及井網(wǎng)格的變量使用全隱式方法求解。假如一口井完井時(shí)射開一個(gè)以上的層，在求解時(shí)將采用充分偶合方式計(jì)算井底壓力，這樣便消除了層狀油藏射開多層完井情況下的不收斂問題。同時(shí) imex 軟件還提供了完整有效的井控制功能，可對(duì)井底和井口壓力，產(chǎn)量、油氣比等進(jìn)行

9、限制。當(dāng)達(dá)到某一限制時(shí)，可根據(jù)用戶要求選擇新的限制條件。矩陣求解方法 imex 軟件使用一種基于不完全高斯消去法的先進(jìn)求解過程，作為正交加速的預(yù)處理部驟，并且專門開發(fā)了用于自適應(yīng)隱式賈可比矩陣的線性方程解法。如需要得到進(jìn)一步的信息請(qǐng)查閱 aimsol 技術(shù)手冊(cè)。 對(duì)于大多數(shù)模擬應(yīng)用，imex 選擇的缺省控制值選項(xiàng)都能使計(jì)算能有效地執(zhí)行，所以用戶并不需要具有矩陣求解方法方面的詳細(xì)知識(shí)。局部網(wǎng)格加密 imex 中包括局部網(wǎng)格加密功能，用戶可以確定油藏中的一個(gè)區(qū)域進(jìn)行細(xì)分，程序?qū)⒆詣?dòng)計(jì)算出內(nèi)部網(wǎng)格的連通情況及傳導(dǎo)率，這種功能可在油田規(guī)模的模擬中，用于研究井周圍地區(qū)的各種影響，而不需對(duì)井進(jìn)行擬函數(shù)處理

10、。使用這種技術(shù)也可有效地模擬靜態(tài)裂縫情況。加密所產(chǎn)生的額外項(xiàng)都可在矩陣求解程序中得到正確處理。尖滅層 使用尖滅選項(xiàng)可以有效地模擬尖滅情況，程序?qū)⒓鉁缇W(wǎng)格從有效網(wǎng)格表中移走，自動(dòng)地將尖滅層的上下層偶合在一起。油藏初始化 可利用給定油氣和油水界面的方式建立油藏的初始情況。靈活的網(wǎng)格系統(tǒng) 在imex中可使用幾種網(wǎng)格選項(xiàng)，笛卡爾坐標(biāo)，圓柱坐標(biāo)和變厚度/變深度網(wǎng)格，在以上任何一種坐標(biāo)中，都可使用二維或三維系統(tǒng)?？勺兣蔹c(diǎn)使用了一種嚴(yán)密的采用變量替換方法的在可變泡點(diǎn)公式，在低飽和區(qū)內(nèi)，對(duì)使用不同 pvt 數(shù)據(jù)油的混合情況能進(jìn)行適當(dāng)?shù)哪M。水區(qū)模型可采用兩種方式對(duì)水區(qū)進(jìn)行模擬，可采用在邊界增加一定的含水網(wǎng)格的

11、方法，或者使用 carter 和 tracy 提出的水區(qū)解析模型。當(dāng)了解水體大小和位置，并且對(duì)于油藏可用包含相對(duì)少量的附加網(wǎng)格的方式就可進(jìn)行模擬的情況下，前一種方法是很有用的；對(duì)于非常大甚至無限大的水體，后一種方法更加有用，在這里需要對(duì)進(jìn)入油藏的水侵量做近似計(jì)算，而通過附加油藏邊界網(wǎng)格的描述方法是不可行的。輸入/輸出單位 國際標(biāo)準(zhǔn)單位(si)，礦場單位(英制)和試驗(yàn)室(修改的si)三種單位制可供選擇?？梢浦残?imex 的程序編寫采用了標(biāo)準(zhǔn)的 fortran77 語言，曾在各種各樣的硬件平臺(tái)上運(yùn)行過，這些包括：ibm大型機(jī)，cdc，cray，honeywell，dec，prime， gould

12、，unisys，apollo，sun，hp，fps 和 ibmpc 386 和 486 的兼容機(jī)等。繪圖系統(tǒng)cmg 的繪圖系統(tǒng) results，將 sr2 文件系統(tǒng)用于模擬輸出的后處理。results 也可用于包括網(wǎng)格設(shè)計(jì)在內(nèi)的輸入數(shù)據(jù)準(zhǔn)備工作。 指導(dǎo)段 _ 介紹 _ 這一指導(dǎo)段是用于對(duì)缺乏經(jīng)驗(yàn)的用戶在關(guān)鍵字輸入系統(tǒng)方面的指導(dǎo)，它并不能替代這個(gè)文件中的用戶參考手冊(cè)。在這一指導(dǎo)段中只討論了一些特殊的關(guān)鍵字和問題，雖然當(dāng)你在建立數(shù)據(jù)文件時(shí)可能會(huì)遇到指導(dǎo)段中“如何做”的問題，但用戶手冊(cè)中包含有每個(gè)關(guān)鍵字的詳細(xì)描述。imex 使用你初始建立的數(shù)據(jù)文件，然后生成三個(gè)或四個(gè)其他文件。imex 的每次運(yùn)行產(chǎn)

13、生一個(gè)輸出文件 (out)，一個(gè)結(jié)果索引文件 (irf) 和一個(gè)主結(jié)果文件 (mrf)。此外根據(jù)用戶的選擇可以產(chǎn)生或不產(chǎn)生一個(gè)可回繞結(jié)果文件 (rrf)。 +-+ | | - .out (輸出) | | | | - .irf (索引) 數(shù)據(jù)文件 - | imex | | | - .mrf (主) | | | | - .rrf (可回繞) +-+ (可選擇) 假如要求進(jìn)行重新啟動(dòng)運(yùn)行，那么將需要幾個(gè)現(xiàn)存文件并產(chǎn)生另外三個(gè)文件，說明見下圖： +-+ 數(shù)據(jù)文件 - | | - .out (輸出) | | input.irf - | | - .irf (索引) | imex | input.mrf -

14、 | | - .mrf (主) | | input.rrf - | | - .rrf (可回繞) (可選擇) +-+ (可選擇)關(guān)鍵字輸入系統(tǒng)中的數(shù)據(jù)段 _當(dāng)你使用關(guān)鍵字輸入系統(tǒng)建立數(shù)據(jù)文件時(shí)，有幾點(diǎn)應(yīng)當(dāng)記住：a)在關(guān)鍵字輸入系統(tǒng)中，有7個(gè)不同的數(shù)據(jù)組。b)這7個(gè)數(shù)據(jù)組必須遵循下列順序輸入: 輸入/輸出控制 油藏描述 組份性質(zhì) 巖石，流體數(shù)據(jù) 初始情況 數(shù)值計(jì)算方法控制 井的數(shù)據(jù)和循環(huán)數(shù)據(jù)。c)屬于各組的關(guān)鍵字不能出現(xiàn)在其他組內(nèi)，除非是在某些特殊情況下。通常是在循環(huán)數(shù)據(jù)組中某些別的數(shù)據(jù)需要改變。d)同樣，也應(yīng)注意關(guān)鍵字在組內(nèi)的輸入順序。 如何建立數(shù)據(jù)文件的文檔方式_使用這些關(guān)鍵字建立你的數(shù)據(jù)文

15、件文檔方式： a) *title1, b) *title2, c) *title3 d) *caseid. 它們必須位于輸入/輸出控制段。這些關(guān)鍵字是可選擇的，可以不包括在數(shù)據(jù)文件中，然而它們對(duì)于文件的文檔化以及區(qū)別相似的數(shù)據(jù)文件是十分有用的，至少應(yīng)使用一個(gè)標(biāo)題，所有的標(biāo)題和情況標(biāo)識(shí)必須包括在單引號(hào)之間。 *title1 和 *caseid 均用于 sr2 文件系統(tǒng)，這一系統(tǒng)用于產(chǎn)生模擬結(jié)果的圖形，*title1 最長可以是40個(gè)字符，而 *title2 和 *title3 每個(gè)最多允許為80個(gè)字符。你也可以使用兩個(gè)關(guān)鍵字標(biāo)識(shí)符即“*”在你的數(shù)據(jù)文件中插入注釋，注釋可以出現(xiàn)在數(shù)據(jù)文件中的任何地

16、方。 e.g. *title1 simulation run #1 - 1989-01-23 *title2 dual porosity problem using the minc option *title3 this is a 12 x 12 x 10 cartesian grid system *caseid run1 * 如果這些標(biāo)題行沒有足夠的空間用于這個(gè)數(shù)據(jù)文件 * 的文檔化，你也可以使用注釋來描述你輸入的數(shù)據(jù)。 * 這樣你就可以在這里或其它任何地方加入附加信息。如何執(zhí)行重新啟動(dòng)運(yùn)行_ 什么是重新啟動(dòng)文件？重新啟動(dòng)文件是一套sr2文件，它至少包括一個(gè)結(jié)果索引文件(irf)和一個(gè)主

17、結(jié)果文件，這套文件可以包括也可以不包括一個(gè)可回繞重啟文件，這取決于用戶使用的選項(xiàng)，關(guān)于模擬變量的重新啟動(dòng)信息（通常稱作重新啟動(dòng)記錄）根據(jù)用戶選擇的頻率寫入這些文件中，只能在寫有重新啟動(dòng)記錄的時(shí)間步長進(jìn)行重新啟動(dòng)運(yùn)行，是否寫重新啟動(dòng)記錄是可選擇的。為什么你需要做重新啟動(dòng)運(yùn)行？由于下述原因你可能想重新啟動(dòng)：a) 做敏感性研究或歷史擬合，b) 改變井的工作制度，c) 在做一個(gè)較大的，時(shí)間較長的作業(yè)之前，執(zhí)行一個(gè)短的模擬運(yùn)行以觀察結(jié)果是否滿意。d) 為在一系列運(yùn)行中節(jié)省執(zhí)行時(shí)間，例如：你已經(jīng)完成了一個(gè)模擬運(yùn)行并且初步結(jié)果看上去不錯(cuò)，現(xiàn)在你想進(jìn)行預(yù)測(cè)運(yùn)行。因?yàn)樵诔跏歼\(yùn)行中你已經(jīng)生成了重新啟動(dòng)記錄，你可以

18、從你的運(yùn)行中選擇一個(gè)時(shí)間步并且重新啟動(dòng)模擬，模擬軟件并不需要在開始時(shí)啟動(dòng)，而是從你選擇的時(shí)間步開始繼續(xù)執(zhí)行。如何進(jìn)行重新啟動(dòng)運(yùn)行？重新啟動(dòng)記錄是可選擇的并不是必須要寫的，但是如果你計(jì)劃做重新啟動(dòng)，你就需要在初始運(yùn)行中生成重新啟動(dòng)記錄。使用 *wrst 和 restart 關(guān)鍵字就可生成重新啟動(dòng)文件， 它們必須位于輸入/輸出控制段中，然而當(dāng)井的情況改變時(shí)，*wrst 也可出現(xiàn)在井的數(shù)據(jù)段中。 *wrst 用于說明寫重新啟動(dòng)記錄的頻率。*restart 說明當(dāng)前的模擬運(yùn)行是一個(gè)重新啟動(dòng)運(yùn)行。假如你希望在最后一個(gè)時(shí)間步開始重啟，那么就可在 *restart 這一行的關(guān)鍵字后其余部分保持空缺，這與缺省

19、情況相同，否則輸入一個(gè)時(shí)間步數(shù)。 例： *restart 30 *wrst 10 為了執(zhí)行一個(gè)重啟運(yùn)行：a) 除了將 *restart 加在數(shù)據(jù)文件中的輸入/輸出控制段之外，不需要對(duì)原來的油藏?cái)?shù)據(jù)做任何修改。b) 如果必要的話，增加時(shí)間步的最大值，或者不與 *maxsteps 一起考慮。c) 你需要輸入文件(在初始運(yùn)行期間生成的)，包括一個(gè) irf 輸入文件，一個(gè) mrf 輸入文件，并且可能還有一個(gè) rrf 輸入文件。d) 你也需要為新的一組輸出文件命名，確保所有的輸入文件都是由相同的模擬運(yùn)行所產(chǎn)生的。輸出文件的控制內(nèi)容控制打印輸出文件的內(nèi)容_控制輸出文件內(nèi)容的打印時(shí)，需使用關(guān)鍵字： a) *

20、wprn b) *outprn. 這些關(guān)鍵字必須出現(xiàn)在輸入/輸出控制段中，也可能在以后參數(shù)改變時(shí)，出現(xiàn)在井的數(shù)據(jù)段中。 *wprn 控制打印頻率，打印內(nèi)容包括網(wǎng)格數(shù)據(jù)，段和租賃面數(shù)據(jù)，井的數(shù)據(jù)和計(jì)算方法控制數(shù)據(jù)，例如牛頓迭代和時(shí)間步的收斂情況。如果不要求在輸出時(shí)打印網(wǎng)格及井的數(shù)據(jù)，則將打印頻率設(shè)置為零。 例： *wprn *well 0 *wprn *grid 0 *wprn *sector 0 假如它們之中的某些不出現(xiàn)在數(shù)據(jù)文件中，缺省情況是在循環(huán)數(shù)據(jù)段中每個(gè) *time 或 *date 關(guān)鍵字所表示的時(shí)間階段打印 *grid 和 *sector 信息。對(duì)于 *well 缺省是在每個(gè)時(shí)間步打

21、印信息，這種缺省將會(huì)產(chǎn)生非常大的輸出文件，能很快地將計(jì)算機(jī)設(shè)備的可用空間添滿。*outprn 用于限制打印什么樣的井?dāng)?shù)據(jù)，網(wǎng)格數(shù)據(jù)，油藏?cái)?shù)據(jù)，以及多少性質(zhì)表。實(shí)際上你可以列出你想打印的網(wǎng)格數(shù)據(jù)類型。井?dāng)?shù)據(jù)的處理差別相當(dāng)大，你可以打印出關(guān)于井的全部數(shù)據(jù)或僅僅打印井的總結(jié)，使用 *outprn *well *layer 可打印出所有井的各層信息，使用 *outprn *well *reservoir 打印油藏條件的生產(chǎn)情況，使用 *outprn *well *brief 可對(duì)每口井印出一行總結(jié)，這也是缺省條件。 圖形文件(sr2)的內(nèi)容控制 _ 為控制 sr2 文件系統(tǒng)的內(nèi)容，使用： a) *ws

22、rf b) *outsrf. 這些關(guān)鍵字必須出現(xiàn)在輸入/輸出控制段中，也可能在以后參數(shù)改變時(shí)，出現(xiàn)在井的數(shù)據(jù)段中。 *wsrf 控制網(wǎng)格數(shù)據(jù)，段數(shù)據(jù)和井?dāng)?shù)據(jù)的輸出頻率，這些信息可以以一定的頻率輸出。如果不要求將網(wǎng)格及井的數(shù)據(jù)輸出到 sr2 文件，則將頻率設(shè)置為零。這可以用于減小文件的輸出量，然而你也可以在隨后的井改變時(shí)對(duì)輸出頻率進(jìn)行修改。如果 *wsrf 不存在于數(shù)據(jù)段內(nèi)，那么缺省條件為在每個(gè)時(shí)間步打印井的信息，當(dāng)在井和循環(huán)數(shù)據(jù)段遇到 *time 關(guān)鍵字時(shí)則輸出 *grid 和 *sector 信息。可通過使用 *outsrf 關(guān)鍵字控制輸出信息量。關(guān)鍵字組合 *wsrf *grid 0 確保

23、了在模擬結(jié)果文件中不寫任何網(wǎng)格信息，在以后井的數(shù)據(jù)改變時(shí)可對(duì)其進(jìn)行修改。*outsrf 用于限制輸出那種網(wǎng)格數(shù)據(jù)，你可以要求在給定的網(wǎng)格區(qū)域輸出確定的變量，對(duì)于網(wǎng)格信息和井的信息可使用分別的變量表。 網(wǎng)格系統(tǒng)描述_對(duì)網(wǎng)格系統(tǒng)進(jìn)行描述時(shí)，你需要使用： a) *grid, b) *kdir, c) *di, d) *dj, e) *dk, 使用 f) *depth 或 g) *paydepth 或是 h) *dtop, 和 i) *dip.在上述關(guān)鍵字中，只有 *kdir 和 *dip 是可選擇項(xiàng)，可以不出現(xiàn)在數(shù)據(jù)段內(nèi)。上述列出的關(guān)鍵字必須出現(xiàn)在油藏描述數(shù)據(jù)段中，并且位于 *null 和 *por

24、 關(guān)鍵字之前。*grid 描述使用的網(wǎng)格系統(tǒng)類型，有三種選擇：規(guī)則直角坐標(biāo)，變深度/變厚度坐標(biāo)，徑向圓柱坐標(biāo)。每種選擇都要求輸入各方向的網(wǎng)格數(shù)。 例： *grid *cart 10 10 6 例： *grid *vari 10 10 6 例： *grid *radial 10 1 151)描述一個(gè)直角坐標(biāo)，網(wǎng)格為10x10x62)描述一個(gè)變深度/變厚度坐標(biāo)，網(wǎng)格與1)相同。3)徑向坐標(biāo)系統(tǒng)用于錐進(jìn)研究，網(wǎng)格為10x1x15。*di, *dj 和 *dk 是要求的關(guān)鍵字，用于輸入網(wǎng)格尺寸，在輸入時(shí)必須使用數(shù)組讀人選項(xiàng)。 例： *grid *cart 10 10 12 *di *con 100.0

25、 *dj *con 100.0 *dk *kvar 25.0 2*50.0 3*40.0 75.0 3*40 2*50例中采用的是規(guī)則的笛卡爾網(wǎng)格系統(tǒng)，i 和 j 方向的每個(gè)網(wǎng)格都是100米寬，而 k 方向的每個(gè)層都具有相同的厚度，而各層的厚度互不相同。應(yīng)當(dāng)注意到，當(dāng)使用 *kdir *up 關(guān)鍵字時(shí)，你的數(shù)據(jù)是從最底層開始輸入的。 局部加密網(wǎng)格描述_ 使用 *refine 關(guān)鍵字描述加密網(wǎng)格的位置，也可和 *hybrid 一起使用描述雜交加密網(wǎng)格，雜交網(wǎng)格是在笛卡爾網(wǎng)格系統(tǒng)中用徑向網(wǎng)格細(xì)分。關(guān)鍵字 *refine 必須出現(xiàn)在油藏描述段中，并位于 *null 和 *por 關(guān)鍵字之前。*ref

26、ine 要求給出加密網(wǎng)格在三個(gè)坐標(biāo)方向的細(xì)分網(wǎng)格數(shù)。使用 *hybrid 關(guān)鍵字，允許你將一個(gè)直角坐標(biāo)下的井網(wǎng)格，使用徑向坐標(biāo)進(jìn)行細(xì)分。在網(wǎng)格加密時(shí)，必須使用 *range 關(guān)鍵字指出加密網(wǎng)格的范圍。例如：在一個(gè) 10x13x3 的規(guī)則笛卡爾坐標(biāo)系統(tǒng)中，要將網(wǎng)格(1, 1, 3 )進(jìn)行加密，在三個(gè)坐標(biāo)方向的加密網(wǎng)格數(shù)分別為2,3,2.其表示方式如下: *refine 2 3 2 *range 1:1 1:1 3:3對(duì)一個(gè)基本網(wǎng)格在任一坐標(biāo)方向最多可劃分為4個(gè)加密網(wǎng)格，如果你對(duì)不同范圍內(nèi)的網(wǎng)格采取不同的細(xì)分形式，可順序輸入 *refine 關(guān)鍵字，并指明基本網(wǎng)格的正確位置。對(duì)于使用 *hybri

27、d 關(guān)鍵字進(jìn)行雜交網(wǎng)格加密，如果在不同范圍采用不同的細(xì)分形式，兩個(gè)細(xì)分區(qū)域的相鄰關(guān)系將受到某種限制(詳見 *refine 關(guān)鍵字)。在某些情況下，細(xì)分程度不同的區(qū)域之間至少要有一個(gè)不細(xì)分的基本網(wǎng)格將兩個(gè)區(qū)域分開。同時(shí)也要注意到，在雙孔隙度條件下，不允許進(jìn)行網(wǎng)格局部加密。在變深度，變厚度坐標(biāo)下可以進(jìn)行局部網(wǎng)格加密，然而，對(duì)于每個(gè)網(wǎng)格來說，加密網(wǎng)格的劃分厚度對(duì)于每個(gè)基礎(chǔ)網(wǎng)格是相同的。如果加密網(wǎng)格的性質(zhì)與對(duì)應(yīng)的基本網(wǎng)格性質(zhì)不一致，使用關(guān)鍵字 *rg 描述加密網(wǎng)格的性質(zhì)，否則，認(rèn)為兩者具有相同的性質(zhì)。雙孔/雙滲模型的使用_可以通過使用雙孔隙度模型進(jìn)行天然裂縫性油藏的模擬，在這些模型中，油藏被離散化為兩

28、個(gè)搭配在一起的連續(xù)域（位于同一空間的兩套網(wǎng)格，其中一個(gè)稱作基巖，另一個(gè)叫裂縫。假設(shè)基巖連續(xù)域是由在空間被裂縫隔離的基巖組成，這些基巖的尺寸在整個(gè)油藏中是變化的，并且是裂縫間隔，方位和寬度的函數(shù)。在計(jì)算的網(wǎng)格中，基巖與裂縫之間的傳導(dǎo)由一個(gè)單獨(dú)的流動(dòng)項(xiàng)表示，通過油藏的流體流動(dòng)由裂縫系統(tǒng)產(chǎn)生，這是流體流動(dòng)的主要通道，在一個(gè)有裂縫的網(wǎng)格，井始終是與裂縫連接的。imex有四種不同類型的雙孔隙度模型可以使用。當(dāng)使用任何一種雙孔隙度選項(xiàng)時(shí)，基巖和裂縫的性質(zhì)都必須確定，這通過在主關(guān)鍵字后使用 *matrix 或 *fracture 修正因子進(jìn)行。例如：確定基巖和裂縫的孔隙度和滲透率，基巖孔隙度是 0.10 而

29、裂縫孔隙度為 0.05，基巖滲透率為 50.0 md，kv/kh = 0.10；裂縫滲透率為 500 md，kv/kh = 0.10。 在油藏描述段輸入如下: *por *matrix *con 0.10 *por *fracture *con 0.05 *permi *matrix *con 50.0 *permj *matrix *con 50.0 *permj *matrix *con 5.0 *permi *fracture *con 500.0 *permj *fracture *con 500.0 *permj *fracture *con 50.0 以相似的方式，可以完成對(duì)于基巖

30、和裂縫網(wǎng)格設(shè)定相對(duì)滲透率表或采用impes/全隱式。例：在開始將基巖網(wǎng)格設(shè)定為impes方式而裂縫網(wǎng)格采用全隱式方式，在井和循環(huán)數(shù)據(jù)段第一個(gè) *date 關(guān)鍵字之后，輸入如下： *aimset *matrix *con 0 *aimset *fracture *con 1 此外在確定使用的雙孔隙度選項(xiàng)類型時(shí)，也必須確定計(jì)算形狀因子的方法，imex 有兩種有效方法用于計(jì)算形狀因子，它們分別基于 warren 和 root 模型或 gilman-kazemi 模型，請(qǐng)參考用戶手冊(cè)和技術(shù)手冊(cè)的有關(guān)段落近一步了解形狀因子的兩種計(jì)算方法。在數(shù)據(jù)文件的油藏描述段中確定雙孔隙度選項(xiàng)類型之后設(shè)定 *shape

31、 *wr 以激活 warren 和 root 方法；相似地使用 *shape *gk 激活 gilman-kazemi 方法。為了計(jì)算形狀因子，imex 需要裂縫的幾何形狀信息，這通過確定裂縫在三個(gè)主要方向(i,j,k)的間隔得到，裂縫間隔在整個(gè)油藏的空間內(nèi)可以是變化的，油藏中的某些部分可以是雙孔隙度而其他部分可通過以設(shè)置裂縫孔隙度為零的方式變?yōu)閱慰紫抖?。在?shù)據(jù)文件的油藏描述段中使用關(guān)鍵字 *difrac, *djfrac 和 *dkfrac 確定三個(gè)主要方向上的裂縫間隔，注意如果確定一個(gè)非常大的裂縫間隔則說明形狀因子，sigma 較小，因而基巖裂縫的傳導(dǎo)將會(huì)降低；確定小的裂縫間隔值表明基巖含

32、有許多裂縫，因而基巖裂縫的傳導(dǎo)將會(huì)增加。確定裂縫間隔為零說明沒有裂縫面垂直于那些坐標(biāo)軸。注意裂縫間隔僅影響基巖裂縫的傳導(dǎo)項(xiàng)，而不影響基巖與裂縫的連接數(shù)。為了將油藏的某一區(qū)域設(shè)定為單孔隙度，必須將網(wǎng)格的裂縫孔隙度置為零。注意為了將基巖或裂縫設(shè)為死結(jié)點(diǎn)，兩者孔隙度都必須設(shè)為零，下面的用法是不可能的： *null *matrix *con 0 或者 *null *fracture *con 0這是因?yàn)?*null 應(yīng)用于網(wǎng)格空間(在這種情況下為基巖和裂縫兩者)而不是分別對(duì)基巖和裂縫網(wǎng)格(有效網(wǎng)格)。 例：一個(gè)由 10x10x5 個(gè)網(wǎng)格組成的油藏模型。所有層都有裂縫，裂縫間隔在 i，j 方向?yàn)?50

33、英尺，而在 k 方向?yàn)?10 英尺。在油藏描述段中輸入： *difrac *con 50.0 *djfrac *con 50.0 *dkfrac *con 10.0層 3 和 4 不包含裂縫(可用單孔隙度系統(tǒng)代表)，相應(yīng)地對(duì)數(shù)據(jù)進(jìn)行修改，假設(shè)裂縫孔隙度為 0.05。 *por fracture *ijk 1:10 1:10 1:5 0.05 1:10 1:10 3:4 0.0四種不同的雙孔隙度選項(xiàng)描述如下。在數(shù)據(jù)文件中只能使用一種類型的雙孔隙度選項(xiàng)。 a) 標(biāo)準(zhǔn)雙孔隙度 -油藏的流體流動(dòng)通過裂縫網(wǎng)絡(luò)產(chǎn)生，基巖網(wǎng)格的基本作用是作為源匯項(xiàng)，在油藏描述段通過使用關(guān)鍵字 *dualpor 激活這一模型

34、。例如：使用標(biāo)準(zhǔn)雙孔隙度模型并用 gilman kazemi 方法計(jì)算形狀因子。在油藏描述數(shù)據(jù)段中輸入如下： *dualpor *shape *gk b) 雙滲透率 - 這個(gè)模型與標(biāo)準(zhǔn)的雙孔隙度模型的不同在于每個(gè)基巖網(wǎng)格與裂縫網(wǎng)格和周圍基巖網(wǎng)格兩者都連接，那么流體流動(dòng)既通過裂縫網(wǎng)絡(luò)也通過基巖網(wǎng)格。雙滲透率選項(xiàng)存在連續(xù)的毛細(xì)管現(xiàn)象的情況可能是重要的。在這些情況下，雙滲透率對(duì)于存在自由氣，油并且沿深度有很大變化的油藏經(jīng)常是重要的，通常只有在垂直 k 方向，基巖對(duì)基巖的傳導(dǎo)才是重要的。通過在油藏描述段使用 *dualperm 激活雙滲透率選項(xiàng)，通過使用基巖傳導(dǎo)率乘子可減小基巖對(duì)基巖的流動(dòng)或?qū)⑵湓O(shè)為零

35、。例如：使用雙滲透率模型以及 gilman kazemi 方法計(jì)算形狀因子，僅在垂直方向設(shè)置基巖對(duì)基巖的傳導(dǎo)。 在油藏描述數(shù)據(jù)段輸入如下數(shù)據(jù)： *dualperm *shape *gk *transi *matrix *con 0.0 *transj *matrix *con 0.0 *transk *matrix *con 1.0 c) 子區(qū)域 - 這個(gè)模型是標(biāo)準(zhǔn)雙孔隙度選項(xiàng)的一變形。它允許用戶在垂直方向?qū)鶐r網(wǎng)格進(jìn)行細(xì)分，其目的是為了更精確地表示由基巖到裂縫的重力排驅(qū)過程，這在垂直方向是很重要的。每個(gè)網(wǎng)格以確定數(shù)目進(jìn)行細(xì)分，以此更精確地表示基巖網(wǎng)格內(nèi)的流體壓力和飽和度。對(duì)于裂縫并不需要更高

36、的分辨率，因此僅對(duì)基巖網(wǎng)格進(jìn)行細(xì)分。在油藏描述段通過 *subdomain ndiv 激活子區(qū)域選項(xiàng)，這里的 ndiv 是每個(gè)基巖網(wǎng)格的細(xì)分層數(shù)。注意 ndiv 在這個(gè)油藏內(nèi)是常數(shù)。例如：將上面雙孔隙度的例子修改為子區(qū)域方式，假設(shè)垂直細(xì)分?jǐn)?shù)為 3。 在油藏描述數(shù)據(jù)段： *subdomain 3 *shape *wr 假設(shè)油藏模型由 10x10x5 個(gè)網(wǎng)格組成，所有層都有裂縫，裂縫間隔在 i，j 方向?yàn)?50 英尺，而在 k 方向?yàn)?10 英尺。將下述數(shù)據(jù)加入油藏描述段： *difrac *ijk 1:10 1:10 1:5 50.0 *djfrac *ijk 1:10 1:10 1:5 50.

37、0 *dkfrac *ijk 1:10 1:10 1:5 10.0 d) minc 方法 - 對(duì)于標(biāo)準(zhǔn)的雙孔隙度概念的另一種擴(kuò)展是由 pruess 和 narasimhan 提出的多重作用連續(xù)域(minc)方法。minc 方法的主要優(yōu)點(diǎn)是使用基巖網(wǎng)格的嵌套式離散化進(jìn)行基巖與裂縫的傳導(dǎo)計(jì)算。這可以十分有效地描述流體的非穩(wěn)定狀況，而這在標(biāo)準(zhǔn)的雙孔隙度模型中是經(jīng)常被忽略的。這種嵌套式離散方法是一維的，可描述壓力，粘滯力和毛細(xì)管力，然而在這種一維的基巖嵌套式細(xì)分方法中不考慮重力。通過在油藏描述段設(shè)定 *minc ndiv 激活 minc 選項(xiàng)，這里的 ndiv 是每個(gè)基巖網(wǎng)格的嵌套層數(shù)。注意 ndiv

38、 的值在整個(gè)油藏中是常數(shù)。例如：將上面雙孔隙度的例子修改為 minc 方式，假設(shè)嵌套細(xì)分?jǐn)?shù)為 3。在油藏描述段中加入： *minc 3 *shape *wr 假設(shè)油藏模型由 10x10x5 個(gè)網(wǎng)格組成，所有層都有裂縫，裂縫間隔在 i，j 方向?yàn)?50 英尺，而在 k 方向?yàn)?10 英尺。 在油藏描述段中加入下述數(shù)據(jù)： *difrac *ijk 1:10 1:10 1:5 50.0 *djfrac *ijk 1:10 1:10 1:5 50.0 *dkfrac *ijk 1:10 1:10 1:5 10.0 總結(jié)：為了激活雙孔/雙滲選項(xiàng)，你可以使用： a) *dualpor b) *minc,

39、c) *subdomain, d) *dualperm, e) *shape, f) *difrac, g) *djfrac, and h) *dkfrac 這些關(guān)鍵字必須出現(xiàn)在油藏描述段，這些不同的選項(xiàng)都是有效的，但在任何給定的數(shù)據(jù)文件中只能使用其中一種。 如果使用了其中任何一種關(guān)鍵字，應(yīng)將它們放在 *null 和 *por 關(guān)鍵字之前。 死結(jié)點(diǎn)的輸入_在給定的網(wǎng)格系統(tǒng)中，有兩種方式可指出存在的死結(jié)點(diǎn):a) *null b) *por它們必須出現(xiàn)在油藏描述數(shù)據(jù)段中， null 是可選擇項(xiàng)，而 por 無論如何必須出現(xiàn)在數(shù)據(jù)段中。*por 是用于確定孔隙度的關(guān)鍵字。使用 *null 時(shí)，用 0

40、 表示死結(jié)點(diǎn)，用 1 表示有效結(jié)點(diǎn)。在下面的例子中，表示 i 方向的網(wǎng)格 1 到 4 ，j 方向的網(wǎng)格 1 到 3，k 方向的網(wǎng)格1 到 3 范圍內(nèi)是死結(jié)點(diǎn)，其余為有效結(jié)點(diǎn)。對(duì)于這個(gè)例子可以使用 *ijk 數(shù)組讀人選項(xiàng)。 例: *null *ijk 1:10 1:10 1:3 1 1:4 1:3 1:3 0可以看出上述例子中的第二行覆蓋了第一行。 *null 是可選擇項(xiàng)，如果不輸入，則假設(shè)所有結(jié)點(diǎn)為有效結(jié)點(diǎn)。*null 的設(shè)定始終覆蓋 *por 的設(shè)定，假如對(duì)于一個(gè)被 *null 設(shè)定的死結(jié)點(diǎn)，即使用 *por 對(duì)此結(jié)點(diǎn)設(shè)置一個(gè)非零的孔隙度，也將被程序忽略。當(dāng)有效結(jié)點(diǎn)設(shè)為非零值時(shí)，也可通過 *

41、por 用零孔隙度指出死結(jié)點(diǎn)。前面的例子同樣可以使用 *por 代替: 例: *por *ijk 1:10 1:10 1:3 0.3 1:4 1:3 1:3 0.0在雙孔/雙滲模型情況下，需要對(duì)基巖和裂縫分別輸入孔隙度值，必須先輸入基巖的孔隙度值，然后輸入裂縫系統(tǒng)的孔隙度，對(duì)于其他數(shù)據(jù)的輸入過程也采用同樣的順序。 例: *por *matrix *ijk 1:10 1:10 1:3 0.3 1:4 1:3 1:3 0.0 *por *fracture *ijk 1:10 1:10 1:3 0.0 8 7:9 1:2 0.4 *mod 8 7:8 1 = 0.45這個(gè)例子同時(shí)也說明了 *mod

42、關(guān)鍵字的使用，它用來修改某些網(wǎng)格的性質(zhì)，將孔隙度由0.40 改變?yōu)?0.45。在雙孔隙系統(tǒng)中，只有當(dāng)基巖和裂縫的孔隙度都為零時(shí)，才將它作為死結(jié)點(diǎn)處理，在一般情況下，基巖和裂縫的孔隙度可以一個(gè)設(shè)為零值，而另一個(gè)不為零。水區(qū)選項(xiàng)的使用_在 imex 中模擬水區(qū)的公式來自 carter 和 tracy 的“一種計(jì)算水侵量的改進(jìn)方法”。這些公式是對(duì) everdingen 和 hurst 公式的一種近似，請(qǐng)參考 van everdingen 和 w. hurst 的“對(duì)于油藏流動(dòng)問題的拉普拉斯變換應(yīng)用”和 l. p. dake 的“油藏工程基礎(chǔ)”第九章。 建立水區(qū)模型時(shí)，可以使用下述關(guān)鍵字: a) *aquifer, b) *aqprop, c) *aqleak d) *aqfunc. 這些關(guān)鍵字放在油藏描述數(shù)據(jù)段中，它們完全是可選擇的，如果使用這些選項(xiàng)，必