多孔介質(zhì)中流量的變量高保真模擬的制作方法
【專利摘要】覆蓋一區(qū)域的精細(xì)計(jì)算機(jī)模型包括精細(xì)網(wǎng)格。故障遵循精細(xì)網(wǎng)格路徑。故障將該區(qū)域分成第一精細(xì)側(cè)和第二精細(xì)側(cè)。流體流的源的模型位于該區(qū)域的第一精細(xì)側(cè)。流體流的匯的模型位于該區(qū)域的第二精細(xì)側(cè)。計(jì)算機(jī)使模型粗糙化。在該粗糙化的模型中，故障遵循經(jīng)粗糙網(wǎng)格路徑。粗糙網(wǎng)格路徑把該區(qū)域分成第一粗糙側(cè)和第二粗糙側(cè)。源的模型和匯的模型位于該區(qū)域的第一粗糙側(cè)。源或匯之一被移動到該區(qū)域的第二粗糙側(cè)。
【專利說明】多孔介質(zhì)中流量的變量高保真模擬
【背景技術(shù)】
[0001]多孔介質(zhì)(諸如地里的烴源地層)中的流量模擬一般涉及使用某些形式的網(wǎng)格化將該多孔介質(zhì)細(xì)分成較小部分或塊。求解該域的這個(gè)細(xì)分(網(wǎng)格化)的多孔介質(zhì)中的流的方程的最普遍的形式是有限差分、有限體積、有限元。無論解的形式，從數(shù)值誤差估算的角度來看，一般地觀察到較精細(xì)的網(wǎng)格(或較小的塊)產(chǎn)生更準(zhǔn)確的答案。然而，一般而言，較精細(xì)的網(wǎng)格要求更多的計(jì)算時(shí)間來產(chǎn)生答案。并行計(jì)算有助于在一定程度上減少計(jì)算占用時(shí)間；然而，為了捕捉盡可能多的場景或?yàn)榱烁玫亓炕嗫捉橘|(zhì)物理性質(zhì)中的不確定性要求很多模擬。經(jīng)常，減小模型的尺寸以減少運(yùn)行每個(gè)模擬的時(shí)間。減少模型的尺寸經(jīng)常包括“粗糙化化”或“提升”模型。粗糙化模型同時(shí)近似地維持精細(xì)網(wǎng)格的性質(zhì)以致較粗糙的或“經(jīng)提升”的模型能夠近似地再現(xiàn)精細(xì)網(wǎng)格化模型中的物理性質(zhì)(在不對精細(xì)模型的結(jié)果進(jìn)行簡單內(nèi)插的情況下)這是一個(gè)挑戰(zhàn)。
【專利附圖】

【附圖說明】
[0002]圖1是精細(xì)網(wǎng)格的視圖。
[0003]圖2是圖1的精細(xì)網(wǎng)格的粗糙化版本的視圖。
[0004]圖3和4示出在粗糙化的網(wǎng)格中將匯(sink)從故障的一側(cè)移動至故障的另一側(cè)。
[0005]圖5是流程圖。
[0006]圖6是系統(tǒng)的框圖。
【具體實(shí)施方式】
[0007]如圖1所示，用于產(chǎn)生將模型特性近似維持在最精細(xì)網(wǎng)格分辨率的變量高保真模型技術(shù)的一個(gè)實(shí)施例從多孔介質(zhì)的精細(xì)模型100開始。在一個(gè)實(shí)施例中，精細(xì)模型100包括N個(gè)精細(xì)網(wǎng)格單元(例如，網(wǎng)格單元105)的網(wǎng)格。在圖1所示示例中，網(wǎng)格被圖示為二維網(wǎng)格?？梢岳斫獾氖蔷W(wǎng)格可以是三維的(即，“3D”)，或網(wǎng)格可包含附加維，諸如時(shí)間。在圖1所示示例中，網(wǎng)格是16X16正方的單元(或塊)，獲得均勻尺寸的256個(gè)塊?？梢岳斫獾氖蔷?xì)模型100的網(wǎng)格可具有其它形狀，諸如非正方矩形、多邊形、非正方菱形、圓形、非圓橢圓形、或其它類似形狀。此外，在圖1中，每個(gè)單元都被圖示為正方形，并且所有單元都具有相同尺寸?？梢岳斫獾氖菃卧槐仨毷钦叫蔚?也就是，它們可以是六邊形、八邊形、或另一個(gè)形狀)，并且它們的形狀或尺寸不必須一致。即，某些單元與其它單元相比，可較大而且具有不同的形狀。
[0008]在一個(gè)實(shí)施例中，N個(gè)精細(xì)網(wǎng)格單元的每一個(gè)代表多孔介質(zhì)的一個(gè)區(qū)域。例如，假定使精細(xì)模型100投射在地球表面的平坦正方形投影上。在該情況下，每一個(gè)單元，例如網(wǎng)格單元105，代表該單元投影在地球表面上的、地球表面的平坦正方形投影的該區(qū)域。
[0009]在一個(gè)實(shí)施例中，如圖1的分解部分所示，N個(gè)精細(xì)網(wǎng)格單元的每一個(gè)，如，精細(xì)網(wǎng)格單元110，通過由精細(xì)網(wǎng)格邊緣135、140、145、150連接的精細(xì)網(wǎng)格節(jié)點(diǎn)115、120、125、130所定義。在一個(gè)實(shí)施例中，可由兩個(gè)精細(xì)網(wǎng)格單元共享精細(xì)網(wǎng)格邊緣135、140、145、150。在一個(gè)實(shí)施例中，精細(xì)網(wǎng)格單元110的所有邊緣被共享。例如,精細(xì)網(wǎng)格單元110和精細(xì)網(wǎng)格單元160共享邊緣150。在一個(gè)實(shí)施例中，只共享網(wǎng)格單元105的兩個(gè)內(nèi)邊緣。
[0010]在一個(gè)實(shí)施例中，N個(gè)精細(xì)網(wǎng)格單元的每一個(gè)相關(guān)聯(lián)于物理性質(zhì)的值。在一個(gè)實(shí)施例中，所述性質(zhì)是多孔性。在一個(gè)實(shí)施例中，所述特性是電阻率。在一個(gè)實(shí)施例中，所述性質(zhì)是另一個(gè)地質(zhì)學(xué)性質(zhì)。
[0011]在一個(gè)實(shí)施例中，由精細(xì)模型100建模的區(qū)域代表包括故障155(在圖1上用虛線圖示)的區(qū)域。在一個(gè)實(shí)施例中，由精細(xì)網(wǎng)格路徑165來表示模型100中的故障，該精細(xì)網(wǎng)格路徑165沿著N個(gè)精細(xì)網(wǎng)格單元(沿故障路徑)的邊緣的故障精細(xì)網(wǎng)格組。在一個(gè)實(shí)施例中，故障代表該區(qū)域的第一精細(xì)側(cè)170 (—般位于故障155的左上)和區(qū)域的第二精細(xì)側(cè)175 (一般位于故障155的右下)之間的結(jié)構(gòu)性不連續(xù)。
[0012]在一個(gè)實(shí)施例中，該模型流體流的源180 (諸如井，在圖1上用實(shí)心圓代表)的模型，與位于該區(qū)域的第一精細(xì)側(cè)上的精細(xì)網(wǎng)格單元相關(guān)聯(lián)，且流體流的匯185 (諸如注入井，在圖1上用空心圓代表)的模型，與位于該區(qū)域的第二精細(xì)側(cè)上的精細(xì)網(wǎng)格單元相關(guān)聯(lián)。因此，在圖1所示示例中，源180和匯185位于代表故障155的精細(xì)網(wǎng)格路徑165的相對側(cè)上。
[0013]在一個(gè)實(shí)施例中，所述技術(shù)接受精細(xì)模型100，并且粗糙化或提升該模型，從而產(chǎn)生M個(gè)粗糙網(wǎng)格單元的粗糙模型，諸如圖2所示的粗糙模型200。在一個(gè)實(shí)施例中，M小于N。即，在一個(gè)實(shí)施例中，粗糙模型200所具有的單元要比精細(xì)模型100少。在一個(gè)實(shí)施例中，M比N小得多。在一個(gè)實(shí)施例中，M比N小數(shù)個(gè)數(shù)量級。正如精細(xì)模型那樣，M個(gè)粗糙網(wǎng)格單元代表多孔介質(zhì)區(qū)域的相應(yīng)部分。在一個(gè)實(shí)施例中，M個(gè)粗糙網(wǎng)格單元的每一個(gè)代表對應(yīng)于由A個(gè)精細(xì)網(wǎng)格單元覆蓋的區(qū)域部分的區(qū)域的一部分，A大于I。例如，圖2中每個(gè)粗糙網(wǎng)格單元代表圖1中四個(gè)精細(xì)網(wǎng)格單元的相同區(qū)域，從而在圖示示例中，A = 4。在一個(gè)實(shí)施例中，粗糙網(wǎng)格單元的尺寸并不一致，從而由每個(gè)粗糙網(wǎng)格單元覆蓋的精細(xì)網(wǎng)格單元的數(shù)量也不同。在一個(gè)實(shí)施例中，精細(xì)網(wǎng)格單元的尺寸、形狀、和其它屬性的上述討論也適用于粗糙網(wǎng)格單元。
[0014]如精細(xì)模型那樣，通過由粗糙網(wǎng)格邊緣連接的粗糙網(wǎng)格節(jié)點(diǎn)來定義每一個(gè)精細(xì)網(wǎng)格單元。
[0015]在一個(gè)實(shí)施例中，通過粗糙網(wǎng)格路徑205來代表粗糙模型200中的故障155，該粗糙網(wǎng)格路徑205沿著M個(gè)粗糙網(wǎng)格單元(它們沿著故障路徑)的邊緣的故障粗糙網(wǎng)格組。在一個(gè)實(shí)施例中，粗糙網(wǎng)格路徑205把該區(qū)域分成一般位于粗糙網(wǎng)格路徑205左上的該區(qū)域的第一粗糙側(cè)210，以及一般位于粗糙網(wǎng)格路徑205左下的粗糙網(wǎng)格路徑的第二粗糙網(wǎng)格側(cè) 215。
[0016]在一個(gè)實(shí)施例中，相當(dāng)詳細(xì)地，精細(xì)模型100對諸如故障155之類的結(jié)構(gòu)性不連續(xù)進(jìn)行計(jì)數(shù)。在一個(gè)實(shí)施例中，粗糙模型200中故障155的重要性取決于故障的透射率。例如，在一個(gè)實(shí)施例中，透射故障被建模為橫跨相鄰單元邊緣的流量系數(shù)的減小。類似地，在一個(gè)實(shí)施例中，密封或非透射故障被建模為橫跨相鄰單元邊緣具有零流量系數(shù)。
[0017]在諸如圖2圖示的情況中，這類處理可導(dǎo)致誤差，在該情況中，在精細(xì)模型100中處于故障相對側(cè)的流體流的源180和流體流的匯185，在粗糙模型200中出現(xiàn)在故障155的相同側(cè)。避免這個(gè)問題的一個(gè)方法是在粗糙模型200中，在故障附近保留精細(xì)模型100精細(xì)網(wǎng)格化。這是 Kefei Wang 和 John E.Killough 在“A New Upscaling Method of RelativePermeability Curves for ReservoirSimulation (油層模擬的相對滲透性曲線的新的提升方法)” (SPE124819)中描述的技術(shù)。
[0018]在一個(gè)實(shí)施例中，如圖3和4所示，通過把流體流的源180的模型或流體流的匯185的模型之一移動到故障的相對側(cè)來避免上述誤差。在一個(gè)實(shí)施例中，這個(gè)動作保留了源180和匯185之間的故障155的透射特性。
[0019]在一個(gè)實(shí)施例中，可對于不止一個(gè)候選粗糙單元做出橫跨故障的源180或匯185的移動。在圖3所示的例子中，源180被移動到單元305，而在圖4中，源被移動到單元405。單元305和405是候選單元。
[0020]在一個(gè)實(shí)施例中，向具有與源180原始駐留在精細(xì)模型100中的該精細(xì)網(wǎng)格單元的物理性質(zhì)值最接近的物理性質(zhì)值的候選單元，做出移動。例如，在一個(gè)實(shí)施例中，物理性質(zhì)是跨越故障的透射率。在一個(gè)實(shí)施例中，在(a)—方面，位于包含源180的精細(xì)網(wǎng)格單元和包含匯185的精細(xì)網(wǎng)格單元之間、由精細(xì)網(wǎng)格路徑165代表的故障155的透射率，(b)位于單元310和單元305之間、由粗糙網(wǎng)格路徑205代表的故障155的透射率、以及(c)位于單元310和單元405之間、由粗糙網(wǎng)格路徑205代表的故障155的透射率，之間進(jìn)行比較。在一個(gè)實(shí)施例中，如果與(c)相比，(b)是(a)的較佳近似，則向單元305做出移動。在一個(gè)實(shí)施例中，如果與(b )相比，(c )是(a)的較佳近似，則向單元405做出移動。
[0021]在一個(gè)實(shí)施例中，如果上述比較沒有提供解，則應(yīng)用另一個(gè)規(guī)則(rule)。在一個(gè)實(shí)施例中，規(guī)則是始終沿相同軸移動。例如，在圖2、3、和4圖示示例中，規(guī)則是可始終沿水平軸移動，在該情況下的移動如圖3所示?？蛇x地，規(guī)則是可始終沿垂直軸移動，在該情況下的移動如圖4所不。在一個(gè)實(shí)施例中，隨機(jī)選擇移動的方向。在一個(gè)實(shí)施例中，在可能移動的方向之間旋轉(zhuǎn)移動方向，即，水平、然后垂直、然后水平，等。
[0022]在一個(gè)實(shí)施例中，規(guī)則是選擇橫跨故障的移動方向，這個(gè)方向盡可能接近垂直于故障的方向。在一個(gè)實(shí)施例中，根據(jù)故障的窗口區(qū)域來確定“故障的方向”。在一個(gè)實(shí)施例中，窗口是粗糙模型200的整個(gè)范圍。
[0023]在一個(gè)實(shí)施例中，規(guī)則是選擇橫跨故障的移動方向，這個(gè)方向最接近精細(xì)模型100中源180和匯185之間的方向。例如，使用圖1-4中圖示示例，包含源180的單元和包含匯185的單元之間的方向是水平的，這會導(dǎo)致選擇圖3所示的水平移動而舍棄圖4所示的垂直移動。
[0024]在一個(gè)實(shí)施例中，確定與粗糙模型200的每個(gè)單元相關(guān)聯(lián)的物理性質(zhì)。在一個(gè)實(shí)施例中，與代表該區(qū)域的第一部分的粗糙網(wǎng)格單元相關(guān)聯(lián)的物理性質(zhì)的值被從代表相同區(qū)域的精細(xì)網(wǎng)格單元的物理性質(zhì)的值中所確定。例如，從精細(xì)網(wǎng)格單元105、190、195、197的物理性質(zhì)的值中確定粗糙網(wǎng)格單元220的物理性質(zhì)的值(見圖2)。在一個(gè)實(shí)施例中，對于流體、壓力和流量系數(shù)的初始分布，使用每個(gè)粗糙網(wǎng)格單元的性質(zhì)的平均或局部單相流建模，直接從精細(xì)模型100中確定粗糙模型200的物理特性的值。
[0025]在一個(gè)實(shí)施例中，確定與每個(gè)粗糙網(wǎng)格單元相關(guān)聯(lián)的物理性質(zhì)包括多相流近似。在一個(gè)實(shí)施例中，使用在Kefei Wang和 John Ε.Killough 的“A New Upscaling Method ofRelative Permeability Curves for Reservoir Simulation (油層模擬的相對滲透性曲線的新的提升方法)”，(SPE124819)中描述的技術(shù)來修改通常所說的相對滲透率函數(shù)從而計(jì)數(shù)粗糙化的網(wǎng)格模型的流的差分。在一個(gè)實(shí)施例中，這個(gè)技術(shù)涉及通過回歸使精細(xì)模型100的精細(xì)網(wǎng)格單元的滲透率與粗糙模型200的粗糙網(wǎng)格單元的滲透率相匹配。在一個(gè)實(shí)施例中，不僅可把這個(gè)技術(shù)應(yīng)用于單元間的流，而且還可以應(yīng)用于個(gè)別源項(xiàng)，以較佳地匹配多孔介質(zhì)的整體的流體生產(chǎn)行為。在一個(gè)實(shí)施例中，已經(jīng)示出這個(gè)技術(shù)不但能夠在模擬周期上匹配精細(xì)模型100，而且還允許粗糙模型200在模擬周期外的可預(yù)測性。
[0026]實(shí)踐中，如圖5所示，在一個(gè)實(shí)施例中，通過執(zhí)行基本精細(xì)模擬以創(chuàng)建精細(xì)模型100 (框510)，就開始增強(qiáng)網(wǎng)格質(zhì)量(框505)。在一個(gè)實(shí)施例中，當(dāng)執(zhí)行提升過程的每個(gè)迭代時(shí)，使用精細(xì)模型100作為基準(zhǔn)。在一個(gè)實(shí)施例中，然后使網(wǎng)格粗糙化(框505)，例如，以形成粗糙模型200。在一個(gè)實(shí)施例中，然后執(zhí)行井修改(框520)從而，例如，相對于故障移動源或匯，以試圖在粗糙模型200中維持精細(xì)模型100的特性。在一個(gè)實(shí)施例中，然后通過平均、局部單相流建模、或如上所述類似處理，使屬性粗糙化(框525)。在一個(gè)實(shí)施例中，如果期望增強(qiáng)經(jīng)提升的粗糙模型200 (框530的“是”分支)，則在精細(xì)模型上執(zhí)行回歸分析，如上所述，以作出多相流近似，并且保存經(jīng)粗糙化的模型(框540)。在一個(gè)實(shí)施例中，如果不期望粗糙模型的增強(qiáng)計(jì)數(shù)多相流(框530的“否”分支)，則保存粗糙模型540。
[0027]在一個(gè)實(shí)施例中，通過重復(fù)框515到540可進(jìn)一步粗糙化模型。
[0028]在一個(gè)實(shí)施例中，以計(jì)算機(jī)程序的形式把執(zhí)行圖5示出的功能的軟件存儲在諸如⑶或DVD之類的計(jì)算機(jī)可讀介質(zhì)605上,如圖6所示。在一個(gè)實(shí)施例中，計(jì)算機(jī)610通過輸入/輸出設(shè)備615從計(jì)算機(jī)可讀介質(zhì)605讀出計(jì)算機(jī)程序，并且把它存儲在存儲器620中準(zhǔn)備執(zhí)行，如果需要的話，通過匯編和鏈接，然后執(zhí)行。在一個(gè)實(shí)施例中，系統(tǒng)通過諸如鍵盤之類的輸入/輸出設(shè)備615接受輸入，并且通過諸如監(jiān)視器或打印機(jī)之類的輸入/輸出設(shè)備615提供輸出。在一個(gè)實(shí)施例中，系統(tǒng)把計(jì)算結(jié)果存儲在存儲器620中，或修改已經(jīng)存儲在存儲器1220中的計(jì)算。
[0029]在一個(gè)實(shí)施例中，通過網(wǎng)絡(luò)625使駐留在存儲器620中的計(jì)算結(jié)果對于遠(yuǎn)程實(shí)時(shí)操作中心630可用。在一個(gè)實(shí)施例中，通過網(wǎng)絡(luò)635使計(jì)算結(jié)果對于遠(yuǎn)程實(shí)時(shí)操作中心630可用，以有助于油井640的計(jì)劃或油井640的鉆挖。
[0030]例如，在一個(gè)實(shí)施例中，使用粗糙模型200確定鉆機(jī)應(yīng)該把鉆柱轉(zhuǎn)移到模型預(yù)測到會具有高滲透率且因此更可能包含有價(jià)值的烴的區(qū)域內(nèi)。相對故障來移動源或匯以便保持粗糙模型的準(zhǔn)確度的能力提高了表示鉆機(jī)將鉆到包含這種有價(jià)值的烴的地下區(qū)域的可能性。
[0031]上文描述了廣泛發(fā)明的一個(gè)或多個(gè)具體實(shí)施例。還可在多種另外的實(shí)施例中進(jìn)行本發(fā)明，因此，不局限于這里上面描述的那些實(shí)施例。為了示意或說明的目的，已經(jīng)提供了本發(fā)明的較佳實(shí)施例的上述描述。并不打算把本發(fā)明無遺漏地或限制于所揭示的精確形式。鑒于上述教導(dǎo)，許多修改或變型都是可能的。旨在使本發(fā)明的范圍不受本詳細(xì)說明的限制，而是由這里所附的權(quán)利要求書所限制。
【權(quán)利要求】
1.一種方法，包括: Ca)計(jì)算機(jī)接受覆蓋一區(qū)域的多孔介質(zhì)的精細(xì)模型，所述模型包括: N個(gè)精細(xì)網(wǎng)格單元的網(wǎng)格，所述N個(gè)精細(xì)網(wǎng)格單元的每一個(gè)代表所述區(qū)域的一部分，通過由精細(xì)網(wǎng)格邊緣連接的精細(xì)網(wǎng)格節(jié)點(diǎn)來定義所述精細(xì)網(wǎng)格單元的每一個(gè)； 具有所述N個(gè)精細(xì)網(wǎng)格單元的每一個(gè)的值的物理性質(zhì)； 故障，遵循沿所述N個(gè)精細(xì)網(wǎng)格單元的邊緣的故障精細(xì)網(wǎng)格組的精細(xì)-網(wǎng)格路徑，所述故障代表所述區(qū)域的第一精細(xì)側(cè)和區(qū)域的第二精細(xì)側(cè)之間的結(jié)構(gòu)性不連續(xù)； 與位于所述區(qū)域的第一精細(xì)側(cè)上的精細(xì)網(wǎng)格單元相關(guān)聯(lián)的流體流的源的模型；以及 與位于所述區(qū)域的第二精細(xì)側(cè)上的精細(xì)網(wǎng)格單元相關(guān)聯(lián)的流體流的匯的模型； (b)通過如下步驟，所述計(jì)算機(jī)粗糙化所述模型: 創(chuàng)建M個(gè)粗糙網(wǎng)格單元的網(wǎng)格，M〈N，以使所述M個(gè)粗糙網(wǎng)格單元的每一個(gè)代表對應(yīng)于A個(gè)精細(xì)網(wǎng)格單元的區(qū)域的一部分，A>1，通過由粗糙網(wǎng)格邊緣連接的粗糙網(wǎng)格節(jié)點(diǎn)來定義所述粗糙網(wǎng)格單元的每一個(gè)； 故障，遵循沿M個(gè)粗糙單元的粗糙網(wǎng)格邊緣的故障粗糙網(wǎng)格組的粗糙網(wǎng)格路徑，所述粗糙網(wǎng)格路徑把所述區(qū)域分成所述區(qū)域的第一粗糙側(cè)和所述區(qū)域的第二粗糙側(cè)； 與流體流的源相關(guān)聯(lián)的精細(xì)網(wǎng)格單元以及與流體流的匯相關(guān)聯(lián)的精細(xì)網(wǎng)格單元，對應(yīng)于所述區(qū)域的第一粗糙側(cè)上的粗糙網(wǎng)格單元； (c)所述計(jì)算機(jī)把所述流體流的源的模型或所述流體流的匯的模型之一從區(qū)域的所述第一粗糙側(cè)上的源粗糙網(wǎng)格單元移動到所述區(qū)域的所述第二粗糙側(cè)上的目的粗糙網(wǎng)格單元; Cd)所述計(jì)算機(jī)使用經(jīng)粗糙化的模型來創(chuàng)建鉆井的計(jì)劃；以及 Ce)所述計(jì)算機(jī)使用所述計(jì)劃鉆井。
2.如權(quán)利要求1所述的方法，其特征在于，還包括: (f)從代表所述區(qū)域的第一部份的所述精細(xì)網(wǎng)格單元的物理性質(zhì)的值中確定與代表所述區(qū)域的第一部份的粗糙網(wǎng)格單元相關(guān)聯(lián)的物理性質(zhì)的值。
3.如權(quán)利要求1所述的方法，其特征在于，還包括: (f)通過平均代表所述區(qū)域的第一部分的精細(xì)網(wǎng)格單元的物理性質(zhì)的值來確定與代表所述區(qū)域的所述第一部份的粗糙網(wǎng)格單元相關(guān)聯(lián)的物理性質(zhì)的值。
4.如權(quán)利要求1所述的方法，其特征在于，還包括: (f)通過代表所述區(qū)域的第一部分的粗糙網(wǎng)格單元的局部單相流建模來確定與代表所述區(qū)域的第一部份的粗糙網(wǎng)格單元相關(guān)聯(lián)的物理性質(zhì)的值。
5.如權(quán)利要求1所述的方法，其特征在于，還包括: (f)通過表示所述區(qū)域的第一部分的粗糙網(wǎng)格單元的多相流建模來確定與代表所述區(qū)域的第一部份的粗糙網(wǎng)格單元相關(guān)聯(lián)的物理性質(zhì)的值。
6.如權(quán)利要求1所述的方法，其特征在于，所述源粗糙網(wǎng)格單元與所述目的粗糙網(wǎng)格單元共享邊緣。
7.如權(quán)利要求1所述的方法，其特征在于，(c)把所述流體流的源的模型或所述流體流的匯的模型之一從源粗糙網(wǎng)格單元移動到位于所述區(qū)域的第二粗糙側(cè)上的目的粗糙網(wǎng)格單元包括:確定在所述區(qū)域的第二粗糙側(cè)上存在與所述源粗糙網(wǎng)格單元共享邊緣的兩個(gè)候選粗糙網(wǎng)格單元；且 從物理性質(zhì)值最接近包含所述流體流的源的模型或所述流體流的匯的模型之一的精細(xì)網(wǎng)格單元的物理性質(zhì)值的兩個(gè)候選粗糙網(wǎng)格單元之一選擇目的粗糙網(wǎng)格單元。
8.如權(quán)利要求1所述的方法，其特征在于，(c)把所述流體流的源的模型或所述流體流的匯的模型之一從源粗糙網(wǎng)格單元移動到位于所述區(qū)域的第二粗糙側(cè)上的目的粗糙網(wǎng)格單元包括: 確定在所述區(qū)域的第二粗糙側(cè)上存在與所述目的粗糙網(wǎng)格單元共享邊緣的兩個(gè)候選粗糙網(wǎng)格單元； 確定所述兩個(gè)候選粗糙網(wǎng)格單元的物理值基本相同；且 應(yīng)用規(guī)則從所述兩個(gè)候選粗糙網(wǎng)格單元中選擇所述目的粗糙網(wǎng)格單元。
9.如權(quán)利要求8所述的方法，其特征在于，所述M個(gè)粗糙網(wǎng)格單元的網(wǎng)格包括軸，且所述規(guī)則包括從所述源粗糙網(wǎng)格單元沿所述軸選擇所述候選粗糙網(wǎng)格單元作為所述目的粗糙網(wǎng)格單元。
10.如權(quán)利要求8所述的方法，其特征在于，所述規(guī)則包括選擇位于最接近垂直于所述粗糙網(wǎng)格路徑的所述源粗糙網(wǎng)格單元的方向上的候選粗糙網(wǎng)格單元作為目的粗糙網(wǎng)格單
11.一種計(jì)算機(jī)程序，存儲在非瞬態(tài)有形計(jì)算機(jī)可讀存儲介質(zhì)中，所述程序包括可執(zhí)行指令，使得計(jì)算機(jī): Ca)接受覆蓋一區(qū)域的多孔介質(zhì)的精細(xì)模型，所述模型包括: N個(gè)精細(xì)網(wǎng)格單元的網(wǎng)格，所述N個(gè)精細(xì)網(wǎng)格單元的每一個(gè)代表所述區(qū)域的一部分，通過由精細(xì)網(wǎng)格邊緣連接的精細(xì)網(wǎng)格節(jié)點(diǎn)來定義所述精細(xì)網(wǎng)格單元的每一個(gè)； 具有所述N個(gè)精細(xì)網(wǎng)格單元的每一個(gè)的值的物理性質(zhì)； 故障，遵循沿所述N個(gè)精細(xì)網(wǎng)格單元的邊緣的故障精細(xì)網(wǎng)格組的精細(xì)-網(wǎng)格路徑，所述故障代表所述區(qū)域的第一精細(xì)側(cè)和區(qū)域的第二精細(xì)側(cè)之間的結(jié)構(gòu)性不連續(xù)； 與位于所述區(qū)域的第一精細(xì)側(cè)上的精細(xì)網(wǎng)格單元相關(guān)聯(lián)的流體流的源的模型；以及 與位于所述區(qū)域的第二精細(xì)側(cè)上的精細(xì)網(wǎng)格單元相關(guān)聯(lián)的流體流的匯的模型； (b)通過如下步驟，粗糙化所述模型: 創(chuàng)建M個(gè)粗糙網(wǎng)格單元的網(wǎng)格，M〈N，以使所述M個(gè)粗糙網(wǎng)格單元的每一個(gè)代表對應(yīng)于A個(gè)精細(xì)網(wǎng)格單元的區(qū)域的一部分，A>1，通過由粗糙網(wǎng)格邊緣連接的粗糙網(wǎng)格節(jié)點(diǎn)來定義所述粗糙網(wǎng)格單元的每一個(gè)； 故障，遵循沿M個(gè)粗糙單元的粗糙網(wǎng)格邊緣的故障粗糙網(wǎng)格組的粗糙網(wǎng)格路徑，所述粗糙網(wǎng)格路徑把所述區(qū)域分成所述區(qū)域的第一粗糙側(cè)和所述區(qū)域的第二粗糙側(cè)； 與流體流的源相關(guān)聯(lián)的精細(xì)網(wǎng)格單元以及與流體流的匯相關(guān)聯(lián)的精細(xì)網(wǎng)格單元，對應(yīng)于所述區(qū)域的第一粗糙側(cè)上的粗糙網(wǎng)格單元； (c)把所述流體流的源的模型或所述流體流的匯的模型之一從所述區(qū)域的所述第一粗糙側(cè)上的源粗糙網(wǎng)格單元移動到所述區(qū)域的所述第二粗糙側(cè)上的目的粗糙網(wǎng)格單元； Cd)使用經(jīng)粗糙化的模型來創(chuàng)建鉆井的計(jì)劃；以及 (e)使用所述計(jì)劃鉆井。
12.如權(quán)利要求11所述的計(jì)算機(jī)程序，其特征在于，還包括使得所述計(jì)算機(jī)執(zhí)行如下步驟的可執(zhí)行指令: (f)從代表所述區(qū)域的第一部份的所述精細(xì)網(wǎng)格單元的物理性質(zhì)的值中確定與代表所述區(qū)域的第一部份的粗糙網(wǎng)格單元相關(guān)聯(lián)的物理性質(zhì)的值。
13.如權(quán)利要求11所述的計(jì)算機(jī)程序，其特征在于，還包括使得所述計(jì)算機(jī)執(zhí)行如下步驟的可執(zhí)行指令: (f)通過平均代表所述區(qū)域的第一部分的精細(xì)網(wǎng)格單元的物理性質(zhì)的值來確定與代表所述區(qū)域的所述第一部份的粗糙網(wǎng)格單元相關(guān)聯(lián)的物理性質(zhì)的值。
14.如權(quán)利要求11所述的計(jì)算機(jī)程序，其特征在于，還包括使得所述計(jì)算機(jī)執(zhí)行如下步驟的可執(zhí)行指令: (f)通過代表所述區(qū)域的第一部分的粗糙網(wǎng)格單元的局部單相流建模來確定與代表所述區(qū)域的第一部份的粗糙網(wǎng)格單元相關(guān)聯(lián)的物理性質(zhì)的值。
15.如權(quán)利要求11所述的計(jì)算機(jī)程序，其特征在于，還包括使得所述計(jì)算機(jī)執(zhí)行如下步驟的可執(zhí)行指令: (f)通過表示所述區(qū)域的第一部分的粗糙網(wǎng)格單元的多相流建模來確定與代表所述區(qū)域的第一部份的粗糙網(wǎng)格單元相關(guān)聯(lián)的物理性質(zhì)的值。
16.如權(quán)利要求11所述的計(jì)算機(jī)程序，其特征在于: 所述源粗糙網(wǎng)格單元與 所述目的粗糙網(wǎng)格單元共享邊緣。
17.如權(quán)利要求11所述的計(jì)算機(jī)程序，其特征在于，當(dāng)(c)把所述流體流的源的模型或所述流體流的匯的模型之一從源粗糙網(wǎng)格單元移動到位于所述區(qū)域的第二粗糙側(cè)上的目的粗糙網(wǎng)格單元時(shí)，所述計(jì)算機(jī): 確定在所述區(qū)域的第二粗糙側(cè)上存在與所述源粗糙網(wǎng)格單元共享邊緣的兩個(gè)候選粗糙網(wǎng)格單元；且 從物理性質(zhì)值最接近包含所述流體流的源的模型或所述流體流的匯的模型之一的精細(xì)網(wǎng)格單元的物理性質(zhì)值的兩個(gè)候選粗糙網(wǎng)格單元之一中選擇目的粗糙網(wǎng)格單元。
18.如權(quán)利要求11所述的計(jì)算機(jī)程序，其特征在于，當(dāng)(c)把所述流體流的源的模型或所述流體流的匯的模型之一從源粗糙網(wǎng)格單元移動到位于所述區(qū)域的第二粗糙側(cè)上的目的粗糙網(wǎng)格單元時(shí)，所述計(jì)算機(jī): 確定在所述區(qū)域的第二粗糙側(cè)上存在與所述目的粗糙網(wǎng)格單元共享邊緣的兩個(gè)候選粗糙網(wǎng)格單元； 確定所述兩個(gè)候選粗糙網(wǎng)格單元的物理值基本相同；且 應(yīng)用規(guī)則從所述兩個(gè)候選粗糙網(wǎng)格單元中選擇所述目的粗糙網(wǎng)格單元。
19.如權(quán)利要求18所述的計(jì)算機(jī)程序，其特征在于，所述M個(gè)粗糙網(wǎng)格單元的網(wǎng)格包括軸，且所述規(guī)則包括從所述源粗糙網(wǎng)格單元沿所述軸選擇所述候選粗糙網(wǎng)格單元作為所述目的粗糙網(wǎng)格單元。
20.如權(quán)利要求18所述的計(jì)算機(jī)程序，其特征在于，所述規(guī)則包括選擇位于最接近垂直于所述粗糙網(wǎng)格路徑的所述源粗糙網(wǎng)格單元的方向上的候選粗糙網(wǎng)格單元作為目的粗糙網(wǎng)格單元。
【文檔編號】G06G7/50GK103477345SQ201180070050
【公開日】2013年12月25日 申請日期:2011年4月12日 優(yōu)先權(quán)日:2011年4月12日
【發(fā)明者】J·E·基洛夫 申請人:界標(biāo)制圖有限公司