收斂與加快計算

1、上周收到的問題大多數(shù)是關(guān)于如何解決模型收斂性問題以及如何加快模型的計算。收斂性問題可以說是所有從事數(shù)模工作的人員都會面臨的問題，本文將以eclipse軟件為例從兩方面介紹收斂性問題。第一方面介紹數(shù)值模擬計算與收斂有關(guān)的一些概念。第二部分介紹如何通過修改模型數(shù)據(jù)來加速計算，解決收斂性問題。一：數(shù)模計算的收斂性： 在了解收斂性之前，應該首先了解幾個基本概念：1。報告步： 一個數(shù)模作業(yè)包括多個報告步，報告步是用戶設(shè)置要求多長時間輸出運行報告，比如可以每個月，每季度或每年輸出運行報告，運行報告包括產(chǎn)量報告和動態(tài)場（重啟）報告。在eclipse軟件中，報告步是通過dates和tstep關(guān)鍵字來設(shè)置的。2

2、。時間步： 一個報告步包括多個時間步，時間步是軟件自動設(shè)置（vip需要用戶設(shè)置）即通過多個時間步的計算來達到下一個報告步，以eclipse為例，假如報告步為一個月，在缺省條件下，eclispe第一個時間步取一天，然后以三倍增加，即第二個時間步取三天，然后取九天，下一個時間步是17天來達到30天的報告步，然后會以每30天的時間步來計算。時間步可以通過tuning關(guān)鍵字來修改。3。非線形迭代：一個時間步包括多次非線形迭代。在缺省情況下，eclipse如果通過12次的非線形迭代沒有收斂，eclipse將對時間步減小10倍。比如下一個時間步應該是30天，如果通過12次的迭代計算不能達到收斂，eclip

3、se將把時間步縮短為3天。下一個時間步將以1.25倍增長，即3.75天，4.68天，。如果在計算過程中經(jīng)常發(fā)生時間步的截斷，計算將很慢。4。線形迭代： 一個非線形迭代包括多次線形迭代。線形迭代是解矩陣。 在eclipse輸出報告prt文件中可以找到時間步，迭代次數(shù)的信息， step 10 time= 100.00 days ( +10.0 days rept 5 its) (1-feb-2008) “step 10” : 說明這是第10個時間步。 “time= 100.00 days”： 說明現(xiàn)在模擬到第100天。 “+10.0 days”： 說明這個時間步是10天。 “rept ： 說明為什

4、么選10天做為時間步，rept是指由于到了下一個報告步。 “5 its: 說明此10天時間步需要5次非線形迭代。 ”(1-feb-2008)“： 現(xiàn)在的模擬時間。 模擬計算的時間取決于時間步的大小，如果模型沒有發(fā)生時間步的截斷而且能保持長的時間步，那表明該模型沒有收斂性問題，反之如果經(jīng)常發(fā)生時間步截斷，那模型計算將很慢，收斂性差。時間步的大小主要取決于非線形迭代次數(shù)。如果模型只用一次非線形迭代計算就可以收斂，那表明模型很容易收斂，如果需要2到3次，模型較易收斂，如果需要4到9次，那模型不易收斂，大于10次的化模型可能有問題，如果大于12次，時間步將截斷。在prt文件中如果看到以下信息： pro

5、blem: at time 200 days (1-feb-2009): non-linear equation convergence failure iteration limit reached - time step chopped from 10 step 20 time= 200.00 days ( +1.0 days chop 5 its) (1-feb-2009) 那表明時間步發(fā)生了截斷。 （注： 另外如果你見到如下信息： warning at time 0.0 days (1-mar-2004): linear equations not fully converged -

6、run may go faster if you increase litmx (=25 - tuning keyword) 你可以不必管。這只是線形方程不收斂） 除了rept,chop外，在rpt文件中還常見以下信息來表明為什么選擇現(xiàn)在的時間步： init: 表明是初始時間步 trnc: 為滿足時間截斷誤差 mins: 最小時間步 maxs: 最大時間步 half: 接近報告步時的時間步取半 diff: 時間步截斷chop之后的增長 （在eclipse技術(shù)手冊的第125頁還會找到更多） 如果模型中有很多chop,diff,mins,那模型有嚴重的收斂性問題。 二：如何解決數(shù)模計算的收斂性問題

7、： 模型不收斂的原因很多，網(wǎng)格參數(shù)，屬性參數(shù)，流體pvt參數(shù)，巖石相滲曲線，毛管壓力曲線，相滲曲線端點標定，初始化，井軌跡，垂直管流表都會造成模型不收斂，下面分別介紹各部分如何造成不收斂及如何解決。1。網(wǎng)格部分：網(wǎng)格正交性差和網(wǎng)格尺寸相差太大是導致不收斂的主要原因之一。正交性差會給矩陣求解帶來困難，而網(wǎng)格尺寸相差大會導致孔隙體積相差很大，大孔隙體積流到小孔隙體積常會造成不收斂。 解決辦法： 網(wǎng)格正交性差通常是在建角點網(wǎng)格時為描述斷層或裂縫的走向而造成的。在此情況下，最好能使邊界與主斷層或裂縫走向平行，這樣一方面網(wǎng)格可以很好地描述斷層或裂縫，另一方面正交性也很好。 在平面上最好讓網(wǎng)格大小能夠較均

8、勻，在沒有井的地方網(wǎng)格可以很大，但最好能夠從大到小均勻過渡。縱向上有的層厚，有的薄，最好把厚層能再細分。在檢查模型時應該每層每層都在三維顯示中檢查。徑向局部網(wǎng)格加密時里面最小的網(wǎng)格不要太小。在eclipse里用minpv關(guān)鍵字可以把小于設(shè)定孔隙體積的網(wǎng)格設(shè)為死網(wǎng)格，這樣通常會有用。2。屬性參數(shù)： 不合理的插值計算會導致屬性分布很差，如果是從地質(zhì)模型粗化為數(shù)模模型，通常問題不大，只是有時候數(shù)模人員自己插值時會有問題。 解決辦法： 有可能盡量用地質(zhì)模型的數(shù)據(jù)，自己插值時可以加一些控制點使屬性合理分布。x,y方向的滲透率最好相等或級差不大。在井連通網(wǎng)格的z方向滲透率不要設(shè)為0，如果想控制垂向流動，可

9、給一個很小的值。3。流體pvt參數(shù) 流體pvt參數(shù)會有兩種可能的問題，一是數(shù)據(jù)不合理導致了負總壓縮系數(shù)，二是壓力或氣油比范圍給的不夠?qū)е履Ｐ蛯vt參數(shù)進行了外插。 解決辦法： 檢查prt文件中的warning信息，如果在油藏壓力范圍內(nèi)有負總壓縮系數(shù)的警告，應該修改pvt參數(shù)，否則的化會有收斂性問題。如果負總壓縮系數(shù)是在油藏壓力范圍之外，可以忽略該警告。此部分的修改主要可以小規(guī)模修改油和氣的fvf和rs。關(guān)于eclipse如何計算負總壓縮系數(shù)，參照eclipse技術(shù)手冊第936頁。 在eclipse中加extrapms關(guān)鍵字可以要求輸出如果發(fā)生pvt插值后的警告信息。在提供pvt表時，壓力應該

10、覆蓋所有范圍，包括注水后的壓力上升。rs值也應該考慮到氣在油中的重新溶解。4。巖石相滲曲線和毛管壓力曲線。 eclipse不會對輸入模型的相滲曲線和毛管壓力曲線進行光滑，將會應用每一個輸入飽和度和相滲值，所以要保證輸入的參數(shù)是合理的。通常的問題有： （1）： 飽和度和相對滲透率的數(shù)據(jù)位數(shù)過多。 （2）： 飽和度值太接近，導致相滲曲線的傾角變化很大。 （3）： 飽和度有很小變化但相對滲透率發(fā)生了很大變化。 解決辦法： （1）：飽和度和相對滲透率最多給兩位小數(shù)就夠了。 （2）：檢查相滲曲線的導數(shù)?梢雜肊clipse中的scal模塊），導數(shù)要光滑。 ?3）：將臨界飽和度和束縛飽和度設(shè)為不同的值。5。

11、端點標定： 在應用端點標定時，有時標定完后的相滲曲線傾角很大，標定后的毛管了很大。 解決辦法： 在三維顯示中檢查標定完的pcw,可以給pcw一個最大值來控制毛管壓力。 輸出每個網(wǎng)格標定后的相滲曲線進行檢查。6。初始化： 初始化最容易發(fā)生的問題是在初始時模型不穩(wěn)定，流體在初始條件下就會發(fā)生流動，這也會導致模型不收斂。造成模型初始不穩(wěn)定的主要有： （1） 手工賦網(wǎng)格飽和度和壓力值。 （2）擬合初始含水飽和度。 解決辦法： （1）盡量不要直接為網(wǎng)格賦壓力和飽和度值，盡量由模型通過油水界面及參考壓力來進行初始化計算。 （2）要想擬合地質(zhì)提供的初始含水飽和度分布，應該進行毛管壓力的端點標定，這樣毛管壓力

12、會穩(wěn)住每個網(wǎng)格的水，在初始條件下不會流動。 （3）可以通過讓模型在沒有任何井的情況下計算十年來檢查初始條件下模型是否穩(wěn)定，如果10年的計算模型壓力和飽和度度沒有變化，說明模型初始是穩(wěn)定的。7。井軌跡： 在進行井處理時井可能以之字型在網(wǎng)格中竄過，有可能發(fā)生井的實際竄過方向與模型關(guān)鍵字定義的方向不符，這也會導致不收斂。 解決辦法： （1） 在三維顯示中檢查井軌跡。 （2） 如果井已經(jīng)關(guān)掉，在模擬時不要給零產(chǎn)量，要用關(guān)鍵字把井關(guān)掉。 （3） 檢查井射孔，井不要射在孤立的網(wǎng)格上。8：垂直管流曲線： 有了垂直管流曲線很容易導致模型不收斂，這有兩種可能： （1）曲線有交叉。 （2）曲線發(fā)生了外插。 解決辦

13、法： （1） 用前處理軟件（eclipse中的vfpi)檢查曲線。 （2）在eclipse中加extrapms關(guān)鍵字可以要求輸出如果發(fā)生vfp插值后的警告信息。 （3）曲線應該覆蓋所有井口壓力，含水，油氣比及產(chǎn)量。 (4）在eclipse用wvfpexp。9：其他解決辦法：如果模型數(shù)據(jù)沒有問題，可以調(diào)整模擬器的收斂計算參數(shù)，對于eclipse，可以做以下調(diào)整：（1）調(diào)整tuning中的最大時間步。如果模型每計算到30天就會截斷時間步，可以將最大時間步調(diào)整為20天，這樣計算會快很多。（2）調(diào)整tuning中的最大線形迭代次數(shù)到70次。（3）降低tuning中的線形收斂誤差標準（4）對于組分模型用

14、freezepc和dpcdt問題00415：謝謝前兩次的回答，非常感謝，這次又出了新問題，請看看，萬分感謝，在prt文件中出現(xiàn)以下提示，不知道什么原因，請指教-warning at time 1318.0 days (11-mar-1994): linear equations not fully converged - run may go faster if you (a) increase nstack (= 15 - runspec section) 18 or (b) make nstack negative ( extends stack on disk ) -problem at

15、 time 1318.0 days (11-mar-1994): non-linear equation convergence failure iteration limit reached - time step chopped from 10.00 nstack表示什么啊，模型前面收斂很好，以后就不行了，速度非常慢，請問可能是什么原因阿謝謝您的回答回答：我在網(wǎng)站內(nèi)寫了如何檢查模型和加快計算的辦法。你好好讀讀。問題00382：老師你好：-problem at time 12484.0 days ( 1-aug-1994): un-physical value for 1/bg in pvt

16、 table 1 well 3-d54 pg = -4.5045e+01 rv = 0.0000e+00 1/bg = -4.7583e+01這是什么呢問題?怎么解決呀?不勝感謝! 回答：這是模型不收斂造成的。檢查你的pvt部分有沒有數(shù)據(jù)外插，檢查你的模型的收斂性。 昨天的問題經(jīng)過今天查看后，線性不收斂問題還是大量存在，又出現(xiàn)了新問題，如果您看的話就不用看千前面兩封郵件了：）這里我再問問一個問題，當昨天我給的pvto為下面的情況時，出現(xiàn)了如下警告：warning：-warning at time 0.0 days ( 1-jan-1981): unphysical value -4.10162

17、 for oil fvf found in oil pressure table 1 at sample pressure 25.75862 due to extrapolation above the maximum entered bubble point pressure 3.03000 今天我將fvf整體調(diào)小，雖然仍然告訴我這是可疑數(shù)據(jù)，但警告已經(jīng)沒有了，為什么警告消失了呢？這里什么數(shù)據(jù)可疑呢？我應該如何修改呢？還有，結(jié)果提示我sgfn和sof3也可疑： 這是為什么呢？：）謝謝：）回答： sof3應該沒關(guān)系，主要是你的swfn應該從0.362開始。 你的pvto壓力數(shù)據(jù)增加了十倍，結(jié)果

18、當然不同了問題00091: 我在算一個ecl模型，沒有任何提示（沒有error,warning problem 也沒有超過限制），算到93年就停止了，我從93年作了個重啟，又可以算下去了，不知道這是什么原因，勞駕月光姐姐幫忙。還有另一個問題，我的模型的有效網(wǎng)格是23萬，算得比較慢，我想讓他算得快一點，可不可以改一下tuning 的設(shè)置，改tuning 的哪個關(guān)鍵詞比較合適，該的幅度與多大，改后和不改的差別主要在什么地方（除了計算速度）?；卮穑?你如果用linux算，有可能會發(fā)生計算停止但沒有錯誤信息。你這種情況不好說，你再算一遍，看看是不是還在93年停止。 我計劃寫如何加速的辦法，不過一直沒有

19、時間。tuning關(guān)鍵字你只能改第二項，不要改別的。你檢查你的模型收斂性如何，如果收斂性好的化，不用修改。檢查所有的warning信息。問題00104: eclipse 計算太慢，您看一下是什么原因? 剛開始很快，前三年用3分鐘，但是后來越來越慢，不知道是什么原因？自己覺得可能是： 1 是不是不收斂了？ 2 動態(tài)數(shù)據(jù)中是一次性把所有的井都設(shè)開，然后給產(chǎn)量0，是不是對動態(tài)數(shù)據(jù)有什么影響？ 3 從移植過來是不是靜態(tài)有問題在vip上用2.5小時就夠了，用小時？一般遇到不收斂的問題怎么解決？回答：1.你可以定義所有井的射孔（compdat)和產(chǎn)量控制(wconhist)，然后你要用welopen將井關(guān)

20、掉，不要給產(chǎn)量0，這樣會快很多，而且結(jié)果也正確。你給0其實結(jié)果是不對的，因為會產(chǎn)出其他相。2。你有很多網(wǎng)格孔隙度很大但滲透率為0，這些網(wǎng)格你可以設(shè)為死網(wǎng)格。3。你需要把schedule部分tuning關(guān)鍵字改為： tuning 1* 15 / / /4。你可以在grid部分加minpv關(guān)鍵字為： minpv 10 /5。你可以把wconinje改為wconinjh關(guān)鍵字。問題00156:對不起，上封郵件忘加附件了我定的timestep是一個月，可在運行時是 基本上以0.1day的速度前進，我的模擬時間是20年，得算到什么時候啊。很可能是收斂性的問題，可看不懂具體問題出在哪。我把運行時的提示放在

21、附件里，肯請您幫我分析一下?；卮穑耗愕哪Ｐ褪怯惺諗啃詥栴}。你把模型發(fā)給我，我?guī)湍憧纯?。我在tuning把時步變大也不管用，在report里都是一些看不太明白的警告信息，沒有出錯信息，很著急啊，愁死我了，領(lǐng)導總吹我。我在把我的原程序在附件里發(fā)過去（怕中午發(fā)的您沒收到），請您幫我看一看?；卮穑褐饕悄愕哪Ｐ蛦栴}，井附近的很多網(wǎng)格有交錯，這及其容易導致不收斂。這是由于在你建模時對斷層處理不合理造成的。你可以在flogrid中編輯你的網(wǎng)格。另外你為什么把你的第七和八平衡區(qū)的equil值都設(shè)為0問題00161:首先對您的幫助表示感謝！我有一個問題要問的是：我在做一區(qū)塊的數(shù)值模擬時，出現(xiàn)不收斂的情況。手冊

22、上說，不收斂通常由于pvt，相滲曲線不光滑引起的，我把這些都修改了，足夠光滑，還是不收斂，請問可能是什么原因引起的？回答：收斂性是一個復雜問題，很多因數(shù)都會導致不收斂，你的模型最好正交性要好一些，網(wǎng)格的尺寸差別不要很大，屬性分布也應該盡量均勻，模型的初始化要穩(wěn)定，不應該在初始是發(fā)生流體流動。井的走向也要檢查。你最好把所有的warning信息先檢查一下。我不久會寫如何解決收斂性問題，你現(xiàn)在如果解決不了，可以把模型發(fā)給我。上周收到的問題大多數(shù)是關(guān)于如何解決模型收斂性問題以及如何加快模型的計算。收斂性問題可以說是所有從事數(shù)模工作的人員都會面臨的問題，本文將以eclipse軟件為例從兩方面介紹收斂性問

23、題。第一方面介紹數(shù)值模擬計算與收斂有關(guān)的一些概念。第二部分介紹如何通過修改模型數(shù)據(jù)來加速計算，解決收斂性問題。一：數(shù)模計算的收斂性： 在了解收斂性之前，應該首先了解幾個基本概念：1。報告步： 一個數(shù)模作業(yè)包括多個報告步，報告步是用戶設(shè)置要求多長時間輸出運行報告，比如可以每個月，每季度或每年輸出運行報告，運行報告包括產(chǎn)量報告和動態(tài)場（重啟）報告。在eclipse軟件中，報告步是通過dates和tstep關(guān)鍵字來設(shè)置的。2。時間步： 一個報告步包括多個時間步，時間步是軟件自動設(shè)置（vip需要用戶設(shè)置）即通過多個時間步的計算來達到下一個報告步，以eclipse為例，假如報告步為一個月，在缺省條件下，

24、eclispe第一個時間步取一天，然后以三倍增加，即第二個時間步取三天，然后取九天，下一個時間步是17天來達到30天的報告步，然后會以每30天的時間步來計算。時間步可以通過tuning關(guān)鍵字來修改。3。非線形迭代：一個時間步包括多次非線形迭代。在缺省情況下，eclipse如果通過12次的非線形迭代沒有收斂，eclipse將對時間步減小10倍。比如下一個時間步應該是30天，如果通過12次的迭代計算不能達到收斂，eclipse將把時間步縮短為3天。下一個時間步將以1.25倍增長，即3.75天，4.68天，。如果在計算過程中經(jīng)常發(fā)生時間步的截斷，計算將很慢。4。線形迭代： 一個非線形迭代包括多次線形

25、迭代。線形迭代是解矩陣。 在eclipse輸出報告prt文件中可以找到時間步，迭代次數(shù)的信息， step 10 time= 100.00 days ( +10.0 days rept 5 its) (1-feb-2008) “step 10” : 說明這是第10個時間步。 “time= 100.00 days”： 說明現(xiàn)在模擬到第100天。 “+10.0 days”： 說明這個時間步是10天。 “rept ： 說明為什么選10天做為時間步，rept是指由于到了下一個報告步。 “5 its: 說明此10天時間步需要5次非線形迭代。 ”(1-feb-2008)“： 現(xiàn)在的模擬時間。 模擬計算的時間

26、取決于時間步的大小，如果模型沒有發(fā)生時間步的截斷而且能保持長的時間步，那表明該模型沒有收斂性問題，反之如果經(jīng)常發(fā)生時間步截斷，那模型計算將很慢，收斂性差。時間步的大小主要取決于非線形迭代次數(shù)。如果模型只用一次非線形迭代計算就可以收斂，那表明模型很容易收斂，如果需要2到3次，模型較易收斂，如果需要4到9次，那模型不易收斂，大于10次的化模型可能有問題，如果大于12次，時間步將截斷。在prt文件中如果看到以下信息： problem: at time 200 days (1-feb-2009): non-linear equation convergence failure iteration li

27、mit reached - time step chopped from 10 step 20 time= 200.00 days ( +1.0 days chop 5 its) (1-feb-2009) 那表明時間步發(fā)生了截斷。 （注： 另外如果你見到如下信息： warning at time 0.0 days (1-mar-2004): linear equations not fully converged - run may go faster if you increase litmx (=25 - tuning keyword) 你可以不必管。這只是線形方程不收斂） 除了rept,

28、chop外，在rpt文件中還常見以下信息來表明為什么選擇現(xiàn)在的時間步： init: 表明是初始時間步 trnc: 為滿足時間截斷誤差 mins: 最小時間步 maxs: 最大時間步 half: 接近報告步時的時間步取半 diff: 時間步截斷chop之后的增長 （在eclipse技術(shù)手冊的第125頁還會找到更多） 如果模型中有很多chop,diff,mins,那模型有嚴重的收斂性問題。 二：如何解決數(shù)模計算的收斂性問題： 模型不收斂的原因很多，網(wǎng)格參數(shù)，屬性參數(shù)，流體pvt參數(shù)，巖石相滲曲線，毛管壓力曲線，相滲曲線端點標定，初始化，井軌跡，垂直管流表都會造成模型不收斂，下面分別介紹各部分如何造

29、成不收斂及如何解決。1。網(wǎng)格部分：網(wǎng)格正交性差和網(wǎng)格尺寸相差太大是導致不收斂的主要原因之一。正交性差會給矩陣求解帶來困難，而網(wǎng)格尺寸相差大會導致孔隙體積相差很大，大孔隙體積流到小孔隙體積常會造成不收斂。 解決辦法： 網(wǎng)格正交性差通常是在建角點網(wǎng)格時為描述斷層或裂縫的走向而造成的。在此情況下，最好能使邊界與主斷層或裂縫走向平行，這樣一方面網(wǎng)格可以很好地描述斷層或裂縫，另一方面正交性也很好。 在平面上最好讓網(wǎng)格大小能夠較均勻，在沒有井的地方網(wǎng)格可以很大，但最好能夠從大到小均勻過渡。縱向上有的層厚，有的薄，最好把厚層能再細分。在檢查模型時應該每層每層都在三維顯示中檢查。徑向局部網(wǎng)格加密時里面最小的網(wǎng)

30、格不要太小。在eclipse里用minpv關(guān)鍵字可以把小于設(shè)定孔隙體積的網(wǎng)格設(shè)為死網(wǎng)格，這樣通常會有用。2。屬性參數(shù)： 不合理的插值計算會導致屬性分布很差，如果是從地質(zhì)模型粗化為數(shù)模模型，通常問題不大，只是有時候數(shù)模人員自己插值時會有問題。 解決辦法： 有可能盡量用地質(zhì)模型的數(shù)據(jù)，自己插值時可以加一些控制點使屬性合理分布。x,y方向的滲透率最好相等或級差不大。在井連通網(wǎng)格的z方向滲透率不要設(shè)為0，如果想控制垂向流動，可給一個很小的值。3。流體pvt參數(shù) 流體pvt參數(shù)會有兩種可能的問題，一是數(shù)據(jù)不合理導致了負總壓縮系數(shù)，二是壓力或氣油比范圍給的不夠?qū)е履Ｐ蛯vt參數(shù)進行了外插。 解決辦法：

31、檢查prt文件中的warning信息，如果在油藏壓力范圍內(nèi)有負總壓縮系數(shù)的警告，應該修改pvt參數(shù)，否則的化會有收斂性問題。如果負總壓縮系數(shù)是在油藏壓力范圍之外，可以忽略該警告。此部分的修改主要可以小規(guī)模修改油和氣的fvf和rs。關(guān)于eclipse如何計算負總壓縮系數(shù)，參照eclipse技術(shù)手冊第936頁。 在eclipse中加extrapms關(guān)鍵字可以要求輸出如果發(fā)生pvt插值后的警告信息。在提供pvt表時，壓力應該覆蓋所有范圍，包括注水后的壓力上升。rs值也應該考慮到氣在油中的重新溶解。4。巖石相滲曲線和毛管壓力曲線。 eclipse不會對輸入模型的相滲曲線和毛管壓力曲線進行光滑，將會應用

32、每一個輸入飽和度和相滲值，所以要保證輸入的參數(shù)是合理的。通常的問題有： （1）： 飽和度和相對滲透率的數(shù)據(jù)位數(shù)過多。 （2）： 飽和度值太接近，導致相滲曲線的傾角變化很大。 （3）： 飽和度有很小變化但相對滲透率發(fā)生了很大變化。 解決辦法： （1）：飽和度和相對滲透率最多給兩位小數(shù)就夠了。 （2）：檢查相滲曲線的導數(shù)?梢雜肊clipse中的scal模塊），導數(shù)要光滑。 ?3）：將臨界飽和度和束縛飽和度設(shè)為不同的值。5。端點標定： 在應用端點標定時，有時標定完后的相滲曲線傾角很大，標定后的毛管了很大。 解決辦法： 在三維顯示中檢查標定完的pcw,可以給pcw一個最大值來控制毛管壓力。 輸出每個網(wǎng)

33、格標定后的相滲曲線進行檢查。6。初始化： 初始化最容易發(fā)生的問題是在初始時模型不穩(wěn)定，流體在初始條件下就會發(fā)生流動，這也會導致模型不收斂。造成模型初始不穩(wěn)定的主要有： （1） 手工賦網(wǎng)格飽和度和壓力值。 （2）擬合初始含水飽和度。 解決辦法： （1）盡量不要直接為網(wǎng)格賦壓力和飽和度值，盡量由模型通過油水界面及參考壓力來進行初始化計算。 （2）要想擬合地質(zhì)提供的初始含水飽和度分布，應該進行毛管壓力的端點標定，這樣毛管壓力會穩(wěn)住每個網(wǎng)格的水，在初始條件下不會流動。 （3）可以通過讓模型在沒有任何井的情況下計算十年來檢查初始條件下模型是否穩(wěn)定，如果10年的計算模型壓力和飽和度度沒有變化，說明模型初始是穩(wěn)定的。7。井軌跡： 在進行井處理時井可能以之字型在網(wǎng)格中竄過，有可能發(fā)生井的實際竄過方向與模型關(guān)鍵字定義的方向不符，這也會導致不收斂。 解決辦法： （1） 在三維顯示中檢查井軌跡。 （2） 如果井已經(jīng)關(guān)掉，在模擬時不要給零產(chǎn)量，要用關(guān)鍵字把井關(guān)掉。