專(zhuān)利名稱(chēng):用于油氣形成建模的系統(tǒng)和方法
用于油氣形成建模的系統(tǒng)和方法
背景技術(shù):
為了使來(lái)自油氣藏的油氣生產(chǎn)最大化,石油和天然氣公司使用油氣藏模型來(lái)模擬油氣藏提取技術(shù)并且隨后基于標(biāo)識(shí)出的結(jié)果來(lái)實(shí)現(xiàn)實(shí)際的提取。由于計(jì)算機(jī)技術(shù)的進(jìn)步以及油氣藏建模技術(shù)的改進(jìn),油氣藏建模的復(fù)雜度和準(zhǔn)確性都已得到提高。在油氣藏建模領(lǐng)域中,在油氣藏模型準(zhǔn)確性與使用油氣藏模型來(lái)運(yùn)行模擬的速度之間存在折衷。較準(zhǔn)確的油氣藏模型較為復(fù)雜,且產(chǎn)生結(jié)果要花費(fèi)較長(zhǎng)時(shí)間。較不復(fù)雜的油氣藏模型可較快速地或以較少的計(jì)算成本產(chǎn)生結(jié)果,但是可能并不充分地計(jì)及地球物理活動(dòng)和反應(yīng)。因此,渴望尋求較準(zhǔn)確地且較快速地執(zhí)行油氣藏建模的任何技術(shù)。附圖簡(jiǎn)述對(duì)于示例性實(shí)施例的詳細(xì)描述,現(xiàn)將對(duì)附圖進(jìn)行參考,附圖中
圖1示出根據(jù)至少一些實(shí)施例的方法;以及圖2示出根據(jù)至少一些實(shí)施例的計(jì)算機(jī)系統(tǒng)。記法和術(shù)語(yǔ)某些術(shù)語(yǔ)在以下的說(shuō)明書(shū)和權(quán)利要求書(shū)中被通篇用以指示特定系統(tǒng)組件。如本領(lǐng)域技術(shù)人員將領(lǐng)會(huì)的,石油和天然氣公司可以按不同名稱(chēng)來(lái)稱(chēng)呼組件。本文不旨在區(qū)分名稱(chēng)不同但功能相同的組件。在以下討論和權(quán)利要求書(shū)中,術(shù)語(yǔ)“包括”和“包含”是以可擴(kuò)充的方式來(lái)使用的, 并且因而應(yīng)被解釋為意味著“包括但不限于……”。同樣,術(shù)語(yǔ)“耦合”旨在意味著非直接或直接的連接。因而,如果第一設(shè)備耦合到第二設(shè)備,則該連接可以通過(guò)直接連接或通過(guò)經(jīng)由其它設(shè)備和連接的非直接連接?!帮柡投取睉?yīng)表示被建模成分的相對(duì)比例。因此,飽和度可以是經(jīng)建模體積內(nèi)水和油氣的飽和度,或者經(jīng)建模體積內(nèi)不同油氣的飽和度。詳細(xì)描述以下討論針對(duì)本發(fā)明的各實(shí)施例。雖然這些實(shí)施例中的一個(gè)或多個(gè)可以是優(yōu)選的,但所公開(kāi)的各實(shí)施例不應(yīng)被解釋為或以其它方式用于限制本公開(kāi)(包括權(quán)利要求書(shū)) 的范圍。另外,本領(lǐng)域技術(shù)人員將理解,以下描述具有廣泛應(yīng)用,并且對(duì)任一實(shí)施例的討論僅意味著是該實(shí)施例的示例,而并非旨在表示本公開(kāi)(包括權(quán)利要求書(shū))的范圍被限于該實(shí)施例。用于對(duì)油氣形成進(jìn)行建模的有關(guān)現(xiàn)有技術(shù)具有固有的限制。例如,一種用于油氣藏建模的技術(shù)被稱(chēng)為有限差分技術(shù)。有限差分技術(shù)將油氣藏建模為多個(gè)特定大小的網(wǎng)格塊。求解預(yù)測(cè)每個(gè)網(wǎng)格塊內(nèi)的油氣和/或水的壓力的差分方程?;谟?jì)算出的壓力,計(jì)算每個(gè)網(wǎng)格塊的每個(gè)面處的流體流動(dòng)速度。然而,在模型不能容易地計(jì)及被建模時(shí)段(即,時(shí)間步長(zhǎng))內(nèi)跨越一個(gè)以上網(wǎng)格塊的油氣和/或水的流動(dòng)的意義上,有限差分技術(shù)是受限制的。取決于網(wǎng)格塊大小和流體運(yùn)動(dòng)速度,用于有限差分技術(shù)的時(shí)間步長(zhǎng)可被限制成極小的大小(例如,一天或更短)。當(dāng)在油氣藏的壽命上對(duì)油氣藏提取進(jìn)行建模時(shí)(這可能在20 年或更長(zhǎng)的數(shù)量級(jí)上),一天或更短的數(shù)量級(jí)上的時(shí)間步長(zhǎng)可以是非常小的。
另一種不具有小時(shí)間步長(zhǎng)限制的油氣藏建模技術(shù)為流線(xiàn)技術(shù)(也被稱(chēng)為歐拉-拉格朗日技術(shù))。流線(xiàn)技術(shù)最初使用有限差分類(lèi)型技術(shù)來(lái)確定網(wǎng)格塊邊界處的壓力和流體流動(dòng)(即,速度)(歐拉部分),但是隨后使用這些速度跨許多網(wǎng)格塊三維地內(nèi)插流體流動(dòng)(拉格朗日部分)。例如,經(jīng)內(nèi)插的流體流動(dòng)可以在被建模時(shí)段期間流經(jīng)許多網(wǎng)格塊,因此得名 “流線(xiàn)”。用于流線(xiàn)技術(shù)的時(shí)間步長(zhǎng)可以跨長(zhǎng)得多的時(shí)間段(例如,60、90、180天),因此可較快速地對(duì)特定提取技術(shù)的油氣藏反應(yīng)進(jìn)行建模。然而,為了使用流線(xiàn)技術(shù),要忽略流體流動(dòng)的某些物理性質(zhì)。例如,流線(xiàn)技術(shù)不容易地計(jì)及重力;相對(duì)滲透性隨水飽和度變化的改變;毛細(xì)壓力如何影響多孔介質(zhì)中的流體流動(dòng);或者與流線(xiàn)流動(dòng)橫向的流體流動(dòng)(橫向通量)O一種結(jié)合流線(xiàn)技術(shù)計(jì)及流體流動(dòng)的物理性質(zhì)中的至少一些物理性質(zhì)的技術(shù)被稱(chēng)為“算子分裂”。算子分裂可概念化為兩步法過(guò)程;首先應(yīng)用傳統(tǒng)的流線(xiàn)技術(shù)并且流體沿模型中的流線(xiàn)“運(yùn)動(dòng)”;隨后在新的位置處對(duì)靜止的流體應(yīng)用流體流動(dòng)的物理性質(zhì)(諸如重力)。在此方面,算子分裂被稱(chēng)為“顯式”技術(shù),其表示與關(guān)于流線(xiàn)技術(shù)的方程順序地求解關(guān)于流體流動(dòng)的物理性質(zhì)的方程,而不是與關(guān)于流線(xiàn)技術(shù)的方程同時(shí)求解關(guān)于流體流動(dòng)的物理性質(zhì)的方程。雖然考慮被流線(xiàn)技術(shù)忽略的流體流動(dòng)的物理性質(zhì)的算子分裂改善了模型準(zhǔn)確性,但是準(zhǔn)確性提高是有限的。在題 為"Timestep Selection During Streamline Simulation Through Transverse Flux Correction(通過(guò)橫向通量校正來(lái)選擇流線(xiàn)模擬期間的時(shí)間步長(zhǎng))”的石油工程師協(xié)會(huì)論文中,作者Osak^Datta-Gupta和King描述了隱式地考慮橫向通量的流線(xiàn)技術(shù)。然而,由于添加了對(duì)橫向通量的考慮,Osaka等人的技術(shù)在時(shí)間步長(zhǎng)大小方面受到限制。為了使Osaka等人的系統(tǒng)保持?jǐn)?shù)值穩(wěn)定,必須如此選擇時(shí)間步長(zhǎng),以使得最快的波不準(zhǔn)在時(shí)間步長(zhǎng)期間橫穿整個(gè)單元。因此,雖然Osaka等人討論了“改進(jìn)的”流線(xiàn)技術(shù),但卻失去了流線(xiàn)技術(shù)的一個(gè)主要益處(使用大時(shí)間步長(zhǎng)的能力)。各種實(shí)施例涉及執(zhí)行油氣藏建模的系統(tǒng)和方法以及存儲(chǔ)指令的計(jì)算機(jī)可讀存儲(chǔ)介質(zhì),它們具有隱式地計(jì)及物理現(xiàn)象(諸如相對(duì)滲透性和毛細(xì)壓力)以及使用大時(shí)間步長(zhǎng)的能力這兩者的益處。本描述將首先給出字面上的概覽,繼以更加數(shù)學(xué)的處理。各種實(shí)施例涉及邏輯地將(油氣)形成劃分成多個(gè)體積或網(wǎng)格塊。在特定實(shí)施例中,網(wǎng)格塊的數(shù)目可以在數(shù)百萬(wàn)個(gè)網(wǎng)格塊的數(shù)量級(jí)上,但是可以等效地使用更多或更少的此類(lèi)網(wǎng)格塊。在一些實(shí)施例中,網(wǎng)格塊體積相等,但是在其他實(shí)施例中,網(wǎng)格塊可以基于油氣和/或水在該網(wǎng)格塊內(nèi)的運(yùn)動(dòng)的活躍性而具有變化的體積。例如,較小的網(wǎng)格塊可在“活躍”區(qū)域中使用,而較大的網(wǎng)格塊可在流體的運(yùn)動(dòng)較少或沒(méi)有運(yùn)動(dòng)的區(qū)域中使用。對(duì)于每個(gè)網(wǎng)格塊并且在計(jì)及固有的形成壓力以及與源(例如,注入井)和匯(例如,生產(chǎn)井)相關(guān)聯(lián)的壓力的情況下,計(jì)算每個(gè)網(wǎng)格塊邊界處的流體壓力。在特定實(shí)施例中,使用有限差分技術(shù)(即,歐拉技術(shù))來(lái)計(jì)算該壓力?;诿總€(gè)網(wǎng)格塊邊界處的壓力或者更確切地說(shuō)跨網(wǎng)格塊邊界考慮的壓力差,確定流動(dòng)速度?;诹鲃?dòng)速度,在時(shí)間步長(zhǎng)上確定流體飽和度(或質(zhì)量)的級(jí)數(shù)。除非另外指出, 確定時(shí)間步長(zhǎng)結(jié)束時(shí)每個(gè)網(wǎng)格塊中的飽和度(或質(zhì)量)。在特定實(shí)施例中,確定飽和度的級(jí)數(shù)使用拉格朗日技術(shù),并且因此經(jīng)建模流體可跨多個(gè)網(wǎng)格塊“流動(dòng)”。然而,如以上所討論的,該過(guò)程的此步驟的拉格朗日技術(shù)不計(jì)及流動(dòng)的許多物理性質(zhì),這些物理性質(zhì)影響每個(gè)網(wǎng)格塊中計(jì)算出的水飽和度的準(zhǔn)確性。例如,計(jì)算出的水飽和度不計(jì)及重力;相對(duì)滲透性隨水飽和度變化的改變;毛細(xì)壓力如何影響多孔介質(zhì)中的流體流動(dòng);或者與流線(xiàn)流動(dòng)橫向的流體流動(dòng)(橫向通量)。根據(jù)各種實(shí)施例,執(zhí)行類(lèi)似于流線(xiàn)技術(shù)的初始步驟表示對(duì)經(jīng)建模形成中的飽和度(例如,水飽和度)遷移的粗略估計(jì)或初步近似,并且該初步近似隨后被修改或校正以計(jì)及以上提及的物理效應(yīng)中的一些或全部。然而,校正此類(lèi)物理效應(yīng)應(yīng)當(dāng)不會(huì)不利地影響時(shí)間步長(zhǎng)的長(zhǎng)度,如同在以上提及的Osaka等人的論文中的技術(shù)所呈現(xiàn)的情形那樣。具體而言,用初步近似的結(jié)果,各種實(shí)施例計(jì)算每個(gè)網(wǎng)格塊內(nèi)的飽和度的變化乘以單元孔隙量除以時(shí)間步長(zhǎng)大小的值。該值指示在一時(shí)間段期間已發(fā)生的流體流動(dòng)。接下來(lái),再次在每個(gè)網(wǎng)格塊內(nèi),確定流體的總速度。并且最后,該方法轉(zhuǎn)向同時(shí)求解巴克利-萊弗里特 (Buckley-Leverett)方程,這些方程被修改成包括諸如網(wǎng)格塊中的油氣與水之間的相對(duì)滲透性、毛細(xì)壓力、重力、或橫向通量之類(lèi)的考慮中的至少一個(gè),但在特定實(shí)施例中包括多個(gè)這樣的考慮。流體流動(dòng)、流體速度和巴克利-萊弗里特的求解的計(jì)算可執(zhí)行多次,直至該值減小(并且在一些情形中最小化),但是在一些情形中,單次迭代就足夠了。當(dāng)求出解時(shí),這些方程提供對(duì)水飽和度確定的校正。不同于Osaka等人,各種實(shí)施例不會(huì)導(dǎo)致數(shù)值不穩(wěn)定。除非另外指出,這些校正不會(huì)施加時(shí)間步長(zhǎng)限制,因?yàn)檫@些校正能夠跨網(wǎng)格塊邊界“移動(dòng)”飽和度。現(xiàn)在轉(zhuǎn)向根據(jù)各種實(shí)施例的對(duì)校正的更加數(shù)學(xué)的處理。在完成了水飽和度的初步近似的情況下,該方法轉(zhuǎn)向基本上使用下式來(lái)為每個(gè)網(wǎng)格塊計(jì)算流體流動(dòng)的殘差值Rwii, j, k) = - ASw"+1 (i, j, h)
AtBw⑴其中Rw(i,j, k)是特定網(wǎng)格塊的殘差值,Sw(i, j, k)n+1是為該特定網(wǎng)格塊計(jì)算的飽和度,Pv(i,j,k)是該網(wǎng)格塊內(nèi)的孔隙量,At是時(shí)間步長(zhǎng)大小,而B(niǎo)w是流體形成體積因子,并且η是時(shí)間步長(zhǎng)。接下來(lái),基本上使用下式來(lái)確定每個(gè)網(wǎng)格塊的接口處的流體總速度Σ^ΑΦ,
(2)其中ut是總速度,Ti是相位i的上游移動(dòng)性的可傳遞性時(shí)間,Δ Φ是每個(gè)網(wǎng)格塊的接口處的潛在梯度,并且其中相位i是油(O)、水(w)和/或氣體(g)。最后,使用從式( 計(jì)算出的總速度來(lái)求解每個(gè)單元的巴克利-萊弗里特方程,其中迭代地確定解直至誤差或殘差值滿(mǎn)足預(yù)定值,諸如最小值。巴克利-萊弗里特方程基本上采取以下形式
_o] j 叫 ▽人=0。(3)其中Sw是飽和度(例如,水飽和度),Ut是總流體速度,fw是分?jǐn)?shù)流體流動(dòng),t是時(shí)間,并且η是時(shí)間步長(zhǎng)。同樣,式C3)充當(dāng)示例,而不是將該技術(shù)限定成僅水飽和度的解。 也可求解其他飽和度和/或成分。根據(jù)各種實(shí)施例,要被考慮的附加物理效應(yīng)包括在關(guān)于分?jǐn)?shù)流體流動(dòng)fw和/或流體速度項(xiàng)Ut的方程中。例如,在其中應(yīng)用校正以處理網(wǎng)格塊內(nèi)的水和油之間的相對(duì)滲透性的二相系統(tǒng)(即油和水)中并且在僅考慮單個(gè)維度的情況下,fw可采取以下形式
權(quán)利要求
1.一種方法,包括基于實(shí)際的地下油氣形成的數(shù)據(jù)來(lái)制定地下油氣形成的邏輯模型,所述模型包括多個(gè)網(wǎng)格塊;在多個(gè)時(shí)間步長(zhǎng)上模擬所述形成對(duì)油氣提取的反應(yīng),每個(gè)時(shí)間步長(zhǎng)內(nèi)的所述模擬通過(guò)以下步驟進(jìn)行作出所述多個(gè)網(wǎng)格塊中的至少一個(gè)網(wǎng)格塊的飽和度遷移的初步近似,其中至少一個(gè)時(shí)間步長(zhǎng)內(nèi)的飽和度遷移橫跨一個(gè)以上網(wǎng)格塊;并且隨后為沒(méi)有在所述飽和度遷移的初步近似中考慮的效應(yīng)來(lái)校正所述初步近似;以及顯示飽和度位置的視覺(jué)描繪。
2.如權(quán)利要求1所述的方法,其特征在于,所述校正還包括校正從由以下各項(xiàng)構(gòu)成的組中選擇的至少一項(xiàng)重力;當(dāng)水飽和度在網(wǎng)格塊內(nèi)變化時(shí)的相對(duì)滲透性;毛細(xì)壓力;以及橫向通量。
3.如權(quán)利要求1所述的方法,其特征在于,所述校正還包括同時(shí)求解多個(gè)巴克利-萊弗里特方程,所述巴克利-萊弗里特方程針對(duì)所述網(wǎng)格塊中的至少一些網(wǎng)格塊,其中每個(gè)巴克利-萊弗里特方程計(jì)及沒(méi)有在所述初步近似中考慮的效應(yīng)。
4.如權(quán)利要求3所述的方法,其特征在于,同時(shí)求解所述多個(gè)巴克利-萊弗里特方程還包括求解對(duì)從由以下各項(xiàng)構(gòu)成的組中選擇的至少一項(xiàng)進(jìn)行校正的巴克利-萊弗里特方程 重力;當(dāng)水飽和度在網(wǎng)格塊內(nèi)變化時(shí)的相對(duì)滲透性;毛細(xì)壓力;以及橫向通量。
5.如權(quán)利要求1所述的方法,其特征在于,建模還包括對(duì)至少一個(gè)油氣提取點(diǎn)的反應(yīng)進(jìn)行建模。
6.如權(quán)利要求1所述的方法,其特征在于,校正所述初步近似還包括作出所述初步近似,以使得所述時(shí)間段期間的所述飽和度遷移橫跨一個(gè)以上網(wǎng)格塊。
7.如權(quán)利要求1所述的方法,其特征在于,校正所述初步近似還包括校正所述初步近似,以使得所述時(shí)間段期間的所述飽和度遷移橫跨一個(gè)以上網(wǎng)格塊。
8.一種存儲(chǔ)程序的計(jì)算機(jī)可讀介質(zhì),所述程序在由一個(gè)或多個(gè)處理器執(zhí)行時(shí)使得所述處理器在多個(gè)時(shí)間步長(zhǎng)上模擬形成模型對(duì)油氣提取的反應(yīng),每個(gè)時(shí)間步長(zhǎng)內(nèi)的所述模擬使所述處理器作出多個(gè)網(wǎng)格塊中的至少一個(gè)網(wǎng)格塊的飽和度遷移的初步近似,并且在所述形成模型的至少一個(gè)區(qū)域中,所述飽和度遷移橫跨一個(gè)以上網(wǎng)格塊;并且隨后為沒(méi)有在所述飽和度遷移的初步近似中考慮的效應(yīng)來(lái)校正所述初步近似。
9.如權(quán)利要求8所述的計(jì)算機(jī)可讀介質(zhì),其特征在于,在所述處理器進(jìn)行校正時(shí),所述程序還使所述處理器校正從由以下各項(xiàng)構(gòu)成的組中選擇的至少一項(xiàng)重力;當(dāng)飽和度在網(wǎng)格塊內(nèi)變化時(shí)的相對(duì)滲透性;毛細(xì)壓力;以及橫向通量。
10.如權(quán)利要求8所述的計(jì)算機(jī)可讀介質(zhì),其特征在于,在所述程序進(jìn)行校正時(shí),所述程序還使所述處理器同時(shí)求解多個(gè)巴克利-萊弗里特方程,所述巴克利-萊弗里特方程各自針對(duì)所述網(wǎng)格塊中的至少一些網(wǎng)格塊,其中所求解的每個(gè)巴克利-萊弗里特方程計(jì)及沒(méi)有在所述初步近似中考慮的效應(yīng)。
11.如權(quán)利要求10所述的計(jì)算機(jī)可讀介質(zhì),其特征在于,在所述處理器同時(shí)求解時(shí),所述程序還使所述處理器同時(shí)求解對(duì)從由以下各項(xiàng)構(gòu)成的組中選擇的至少一項(xiàng)進(jìn)行校正的所述多個(gè)巴克利-萊弗里特方程重力;當(dāng)飽和度在網(wǎng)格塊內(nèi)變化時(shí)的相對(duì)滲透性;毛細(xì)壓力;以及橫向通量。
12.如權(quán)利要求8所述的計(jì)算機(jī)可讀介質(zhì),其特征在于,在所述處理器作出所述初步近似時(shí),所述程序還使所述處理器對(duì)至少一個(gè)油氣提取點(diǎn)的反應(yīng)進(jìn)行建模。
13.如權(quán)利要求8所述的計(jì)算機(jī)可讀介質(zhì),其特征在于,在所述處理器進(jìn)行校正時(shí),所述程序還使所述處理器校正初步近似,以使得由所述校正為至少一個(gè)區(qū)域施加的飽和度遷移變化橫跨一個(gè)以上網(wǎng)格塊。
14.一種計(jì)算機(jī)系統(tǒng),包括處理器;耦合至所述處理器的存儲(chǔ)器,所述存儲(chǔ)器存儲(chǔ)程序,所述程序在由所述處理器執(zhí)行時(shí)使所述處理器在多個(gè)時(shí)間步長(zhǎng)上模擬形成模型對(duì)油氣提取的反應(yīng),每個(gè)時(shí)間步長(zhǎng)內(nèi)的所述模擬使所述處理器作出多個(gè)網(wǎng)格塊中的至少一個(gè)網(wǎng)格塊的飽和度遷移的初步近似,并且在所述形成模型的至少一個(gè)區(qū)域中,所述飽和度遷移橫跨一個(gè)以上網(wǎng)格塊;并且隨后為沒(méi)有在所述飽和度遷移的初步近似中考慮的效應(yīng)來(lái)校正所述初步近似。
15.如權(quán)利要求14所述的計(jì)算機(jī)系統(tǒng),其特征在于,在所述處理器進(jìn)行校正時(shí),所述程序還使所述處理器校正從由以下各項(xiàng)構(gòu)成的組中選擇的至少一項(xiàng)重力;當(dāng)飽和度在網(wǎng)格塊內(nèi)變化時(shí)的相對(duì)滲透性;毛細(xì)壓力;以及橫向通量。
16.如權(quán)利要求14所述的計(jì)算機(jī)系統(tǒng),其特征在于,在所述程序進(jìn)行校正時(shí),所述程序還使所述處理器同時(shí)求解多個(gè)巴克利-萊弗里特方程,所述巴克利-萊弗里特方程各自針對(duì)所述網(wǎng)格塊中的至少一些網(wǎng)格塊,其中所求解的每個(gè)巴克利-萊弗里特方程計(jì)及沒(méi)有在所述初步近似中考慮的效應(yīng)。
17.如權(quán)利要求16所述的計(jì)算機(jī)系統(tǒng),其特征在于,在所述處理器同時(shí)求解時(shí),所述程序還使所述處理器同時(shí)求解對(duì)從由以下各項(xiàng)構(gòu)成的組中選擇的至少一項(xiàng)進(jìn)行校正的所述多個(gè)巴克利-萊弗里特方程重力;當(dāng)飽和度在網(wǎng)格塊內(nèi)變化時(shí)的相對(duì)滲透性;毛細(xì)壓力; 以及橫向通量。
18.如權(quán)利要求14所述的計(jì)算機(jī)系統(tǒng),其特征在于,在所述處理器進(jìn)行模擬時(shí),所述程序還使所述處理器對(duì)至少一個(gè)非油氣注入點(diǎn)和至少一個(gè)油氣提取點(diǎn)的反應(yīng)進(jìn)行建模。
19.如權(quán)利要求14所述的計(jì)算機(jī)系統(tǒng),其特征在于,在所述處理器進(jìn)行校正時(shí),所述程序還使所述處理器校正初步近似,以使得由所述校正為至少一個(gè)區(qū)域施加的飽和度遷移變化橫跨一個(gè)以上網(wǎng)格塊。
全文摘要
油氣形成建模。解說(shuō)性實(shí)施例中的至少一些是包括在多個(gè)時(shí)間步長(zhǎng)上模擬該形成對(duì)油氣提取的反應(yīng)的方法。在每個(gè)時(shí)間步長(zhǎng)上,該解說(shuō)性方法作出多個(gè)網(wǎng)格塊中的至少一個(gè)網(wǎng)格塊的飽和度遷移的初步近似(其中至少一個(gè)時(shí)間步長(zhǎng)內(nèi)的飽和度遷移橫跨一個(gè)以上網(wǎng)格塊),并且隨后為沒(méi)有在飽和度遷移的初步近似中考慮的效應(yīng)來(lái)校正該初步近似。
文檔編號(hào)A62D3/00GK102576374SQ200980161288
公開(kāi)日2012年7月11日 申請(qǐng)日期2009年9月2日 優(yōu)先權(quán)日2009年9月2日
發(fā)明者J·基洛夫 申請(qǐng)人:界標(biāo)制圖有限公司