用于在儲(chǔ)層模擬系統(tǒng)中進(jìn)行粗化的系統(tǒng)和方法
【專利摘要】描述了用于實(shí)現(xiàn)儲(chǔ)層模擬系統(tǒng)的系統(tǒng)和方法。一個(gè)實(shí)施例是粗化包括多個(gè)精細(xì)網(wǎng)格塊的精細(xì)網(wǎng)格的計(jì)算機(jī)實(shí)現(xiàn)的方法，該精細(xì)網(wǎng)格表示其中具有至少一個(gè)間斷點(diǎn)的地理模型。該方法包括：將若干精細(xì)網(wǎng)格塊編組在一起以形成粗糙網(wǎng)格塊，其中粗糙網(wǎng)格塊中的至少一個(gè)是非標(biāo)準(zhǔn)形狀網(wǎng)格塊；并計(jì)算每對(duì)毗鄰的粗糙網(wǎng)格塊的傳導(dǎo)能力，其中該粗糙網(wǎng)格塊對(duì)中的至少一個(gè)網(wǎng)格塊是非標(biāo)準(zhǔn)形狀網(wǎng)格塊。該計(jì)算包括：計(jì)算每對(duì)毗鄰精細(xì)網(wǎng)格塊的傳導(dǎo)能力；將精細(xì)網(wǎng)格塊對(duì)的每一對(duì)映射到粗糙網(wǎng)格塊對(duì)；并對(duì)于每個(gè)網(wǎng)格塊對(duì)，對(duì)映射至此的精細(xì)網(wǎng)格塊對(duì)的傳導(dǎo)能力求和。
【專利說(shuō)明】用于在儲(chǔ)層模擬系統(tǒng)中進(jìn)行粗化的系統(tǒng)和方法
【背景技術(shù)】
[0001]儲(chǔ)層模擬是儲(chǔ)層工程領(lǐng)域，其采用計(jì)算機(jī)模型來(lái)預(yù)測(cè)儲(chǔ)層中諸如石油、水、和氣體之類的流體的傳輸。儲(chǔ)層模擬器由石油生產(chǎn)商使用以確定如何最佳地開(kāi)發(fā)新油田、以及生成生產(chǎn)預(yù)報(bào)，可基于該生產(chǎn)預(yù)報(bào)得到與開(kāi)發(fā)的油田相關(guān)聯(lián)的投資決策。
[0002]儲(chǔ)層模擬模型一般是使用若干離散的塊實(shí)現(xiàn)的，在本文中可互換地稱其為“塊”、“網(wǎng)格塊”、或“單元”。模型的尺寸可從幾個(gè)塊變化至幾億個(gè)塊。通常，儲(chǔ)層模擬工作流程開(kāi)始于創(chuàng)建包括許多“精細(xì)”網(wǎng)格塊的高分辨率模型，在這一點(diǎn)處模型的尺寸可減少為“粗糙”網(wǎng)格塊以使能在合理的時(shí)間段內(nèi)運(yùn)行模擬。該過(guò)程遠(yuǎn)非自動(dòng)，且此外，取決于所建模的儲(chǔ)層中的斷層和井的數(shù)目、位置、和取向，經(jīng)受模型中的間斷點(diǎn)?？赏ㄟ^(guò)對(duì)模型重新劃分網(wǎng)格來(lái)解決一些問(wèn)題；然而，如果用于創(chuàng)建模型的底層數(shù)據(jù)和/或軟件不可用，則重新劃分網(wǎng)格可能是困難的。
[0003]假設(shè)人們開(kāi)始時(shí)使用高分辨率地質(zhì)單元精細(xì)網(wǎng)格模型，存在用于減小模型尺寸的已有過(guò)程，這一般被稱為提升(upscaling)和/或粗化。提升和粗化這兩者均涉及采樣精細(xì)標(biāo)度模型進(jìn)并創(chuàng)建較粗糙的模型，以嘗試遵守(honor)原始模型的流體性質(zhì)(諸如孔隙體積、傳導(dǎo)能力、和飽和率)。從本質(zhì)上說(shuō)，提升和粗化是求平均過(guò)程并且目的之一在于維持模型的流體特性。如此處所用，粗化是其中通過(guò)移除網(wǎng)格節(jié)點(diǎn)而不改變網(wǎng)格的其他部分來(lái)將網(wǎng)格塊合并到較大的塊中的過(guò)程。提升類似于粗化，然而使用提升的話，網(wǎng)格可被改變并可在較粗糙的網(wǎng)格上被重新采樣。
[0004]本文描述的所有方法在三維(“3D”)是有效的；然而，為簡(jiǎn)單起見(jiàn)，下面將參照二維(“2D”)對(duì)這些方法進(jìn)行描述以使附圖和說(shuō)明不過(guò)分復(fù)雜。圖1A-1E—起示出網(wǎng)格粗化的概念。圖1A示出包括多個(gè)網(wǎng)格塊的8X4網(wǎng)格100，由附圖標(biāo)記102圖示地標(biāo)示。圖1B示出精細(xì)網(wǎng)格塊102的編組以執(zhí)行網(wǎng)格100的2X2粗化，其中沿X方向的每?jī)蓚€(gè)網(wǎng)格塊和沿y方向的每?jī)蓚€(gè)網(wǎng)格塊被編組在一起以形成更大、或粗糙的網(wǎng)格塊，由附圖標(biāo)記104圖示地標(biāo)示。圖1C示出在完成2X2粗化并去除粗糙網(wǎng)格塊104的內(nèi)部節(jié)點(diǎn)之后包括網(wǎng)格100的網(wǎng)格100’。圖1D示出網(wǎng)格100”，其包括在其上已執(zhí)行4X2粗化后的網(wǎng)格100’，得到均勻的甚至更粗糙的網(wǎng)格塊，由附圖標(biāo)記106圖示地標(biāo)示。圖1E示出網(wǎng)格100’’’，它包括在其上已執(zhí)行2X1粗化后的網(wǎng)格100”，得到單個(gè)網(wǎng)格塊108。
[0005]如從圖1A-1E中容易觀察地那樣，每個(gè)粗糙網(wǎng)格塊104-108由若干精細(xì)網(wǎng)格塊102構(gòu)成。特定地，每個(gè)粗糙網(wǎng)格塊104包括4個(gè)精細(xì)網(wǎng)格塊102 ;每個(gè)粗糙網(wǎng)格塊106包括16個(gè)精細(xì)網(wǎng)格塊102 ;而粗糙網(wǎng)格塊108包括32個(gè)精細(xì)網(wǎng)格塊102。
[0006]如圖1A-1E所示，當(dāng)結(jié)合非復(fù)雜儲(chǔ)層模型執(zhí)行時(shí)，粗化是相當(dāng)簡(jiǎn)單和直接的。如果將復(fù)雜性——尤其是間斷點(diǎn)——加入到模型中，該過(guò)程變得更為復(fù)雜。例如，圖2A示出與圖1A的網(wǎng)格100—致的網(wǎng)格200，除了網(wǎng)格200包括間斷點(diǎn)202。圖2B示出與圖2A相同的網(wǎng)格，但不是平坦表面。在任何情況下，為了本文的示例和說(shuō)明，將假設(shè)間斷點(diǎn)202是結(jié)構(gòu)性間斷點(diǎn)，諸如斷層。圖3示出在已經(jīng)執(zhí)行圖1B所示的相同2X2編組后的網(wǎng)格200。如圖3所示，2 X 2粗化導(dǎo)致對(duì)某些塊或單元的自然編組；然而，對(duì)于其它塊或單元，在粗化的單元內(nèi)部可能有斷層，這是不允許的，因?yàn)閱卧沓鲈谄渖蠄?zhí)行流體計(jì)算的均質(zhì)體。如果允許內(nèi)部特征，則要求將這些單元細(xì)分成更小的單元，這與粗化的目的相悖。
[0007]當(dāng)前，為應(yīng)對(duì)圖3所示的情況有四種常見(jiàn)的方法:
[0008]1.無(wú)論如何進(jìn)行粗化，并舍棄出現(xiàn)在粗化的單元內(nèi)任何地方的斷層信息(“無(wú)論如何粗化”);
[0009]2.對(duì)模型重新進(jìn)行網(wǎng)格化以在期望的4X2網(wǎng)格上重新采樣屬性和斷層(“對(duì)模型重新網(wǎng)格化”);
[0010]3.粗化能到的任何地方，但不粗化在單元內(nèi)部有斷層的塊(“粗化能到的任何地方，，);以及
[0011]4.粗化能到的任何地方并然后將這些單元邏輯編組成矩形粗化塊，但仍然留下保持?jǐn)鄬铀璧木?xì)標(biāo)度塊(“粗化能到的任何地方并然后重新編組”)。
[0012]圖4A-7B示出將前面四種方法應(yīng)用到網(wǎng)格200。特定地，圖4A示出在應(yīng)用“無(wú)論如何粗化”方法后包括網(wǎng)格200的網(wǎng)格200’，其中粗糙網(wǎng)格內(nèi)的節(jié)點(diǎn)以及間斷點(diǎn)信息被去除。如圖4A所示，在這種情況下，間斷點(diǎn)202來(lái)來(lái)往往(come and go);它沿某些粗化網(wǎng)格邊的一些部分出現(xiàn)并在其它部分消失。圖4B示出圖4A的所得模型的三維(“3D”)視圖。
[0013]圖5A示出在應(yīng)用“對(duì)模型重新網(wǎng)格化”方法之后包括網(wǎng)格200的網(wǎng)格200”，其中在粗化的4X2網(wǎng)格上對(duì)斷層和屬性重新采樣。該方法保持了間斷點(diǎn)；然而，粗糙塊內(nèi)部的間斷點(diǎn)202的一些部分被重新采樣并被移動(dòng)到最近的邏輯網(wǎng)格塊邊緣以保持間斷點(diǎn)，雖然被重新定位。為清楚起見(jiàn)，從中已去除間斷點(diǎn)的原始位置在圖5A中表示為虛線500。圖5B示出圖5A的結(jié)果模型的3D視圖。
[0014]圖6A示出在應(yīng)用“粗化能到的任何地方并維持精細(xì)標(biāo)度”方法之后包括網(wǎng)格200的網(wǎng)格200’’’，其中粗糙塊內(nèi)不存在間斷點(diǎn)的網(wǎng)格被粗化。然而，對(duì)于其內(nèi)有間斷點(diǎn)的塊來(lái)說(shuō)，粗化被跳過(guò)。這恰好地將間斷點(diǎn)維持在其之前所在處，但限制了模型尺寸可被減少的量。圖6B示出圖6A的所得模型的3D視圖。
[0015]圖7A示出在應(yīng)用“粗化能到的任何地方并然后重新編組”方法之后包括網(wǎng)格200的網(wǎng)格200””，這是圖6A和圖6B所示的方法的邏輯延伸。特定地，在該方法中，嘗試對(duì)一些較精細(xì)(即非粗化)塊進(jìn)行編組。圖7B示出圖7A的所得模型的3D視圖。
[0016]前述方法中的每一種方法均具有缺陷。例如，“無(wú)論如何粗化”方法導(dǎo)致間斷點(diǎn)信息的喪失，這可能顯著地影響不同塊的連接性?！皩?duì)模型重新網(wǎng)格化”方法導(dǎo)致間斷點(diǎn)被重新定位。對(duì)于較小的提升因子，這是可接受的；然而，重新定位可能具有顯著的影響。例如，如果重新網(wǎng)格化發(fā)生在井的附近，間斷點(diǎn)的重新定位可導(dǎo)致井從間斷點(diǎn)的一側(cè)位移至另一偵牝藉此還影響了錯(cuò)誤的斷層塊中的井并且那些層被建模為由井穿孔?！按只艿降娜魏蔚胤健狈椒ū苊饬饲懊鎯煞N方法的缺陷；然而，該過(guò)程的縮放性自然是有限的。例如，如圖6A所示，間斷點(diǎn)得以維持，但網(wǎng)格200的尺寸僅從32塊減小至14塊。對(duì)于較小的粗化因子而言這是可接受的，但對(duì)于較大的模型和較大的粗化因子，保持幾乎所有的精細(xì)標(biāo)度細(xì)節(jié)實(shí)際上與粗化的目的相悖。
[0017]剩下的方法，其中在任何可能的地方粗化網(wǎng)格并將其余的塊邏輯編組為矩形粗化塊，其中精細(xì)標(biāo)度塊留存于任何必要的地方以維持?jǐn)鄬?，盡管該方法可能是四種方法中最好的方法，但就縮放性而言仍然是受限的。此外，它是不唯一的方法，這意味著存在許多不同的方式來(lái)對(duì)單元進(jìn)行編組并且目前是一種冗長(zhǎng)乏味的手動(dòng)過(guò)程。
[0018]附圖簡(jiǎn)述
[0019]可通過(guò)結(jié)合附圖參照下面的說(shuō)明來(lái)獲得對(duì)本說(shuō)明書及其優(yōu)勢(shì)的更完整理解，在附圖中:
[0020]圖1A-1E —起示出對(duì)包含儲(chǔ)層模型的網(wǎng)格進(jìn)行粗化的概念。
[0021]圖2A和圖2B示出表示其中具有間斷點(diǎn)的儲(chǔ)層的一部分的網(wǎng)格。
[0022]圖3示出在已執(zhí)行對(duì)其單元的2X2編組之后圖2的網(wǎng)格。
[0023]圖4A-4B、5A-5B、6A-6B、和圖7A-7B示出用于對(duì)包括其中存在間斷點(diǎn)的儲(chǔ)層模型的網(wǎng)格進(jìn)行粗化的四種不同方法。
[0024]圖8A是適于實(shí)現(xiàn)示例性實(shí)施例的儲(chǔ)層模擬系統(tǒng)的計(jì)算機(jī)系統(tǒng)的框圖。
[0025]圖SB是示出用于執(zhí)行本文描述的實(shí)施例的由圖8A的儲(chǔ)層模擬系統(tǒng)實(shí)現(xiàn)的方法的流程圖。
[0026]圖9A示出表示儲(chǔ)層模型的一部分的網(wǎng)格。
[0027]圖9B示出根據(jù)在非結(jié)構(gòu)化網(wǎng)格中使用的實(shí)施例的特征的鋸齒形網(wǎng)格塊的概念。
[0028]圖9C示出根據(jù)在結(jié)構(gòu)化網(wǎng)格中使用的實(shí)施例的特征的鋸齒形網(wǎng)格塊的概念。
[0029]圖10是根據(jù)一個(gè)實(shí)施例的特征的鋸齒形網(wǎng)格塊的概念的另一示圖。
[0030]圖11示出鋸齒形網(wǎng)格塊和毗鄰的網(wǎng)格塊之間的連接和傳導(dǎo)能力。
[0031]圖12示出根據(jù)一個(gè)實(shí)施例的特征的一對(duì)網(wǎng)格塊之間的傳導(dǎo)能力和滲透性之間的關(guān)系。
[0032]圖13不出本發(fā)明的一個(gè)實(shí)施例,其中基于一對(duì)網(wǎng)格塊之間的半傳導(dǎo)能力來(lái)確定傳導(dǎo)能力。
[0033]圖14示出本發(fā)明的一個(gè)實(shí)施例，其中基于一對(duì)網(wǎng)格塊之間的平均滲透性來(lái)確定傳導(dǎo)能力。
【具體實(shí)施方式】
[0034]為了克服當(dāng)前方法的前面提到的和其它限制，本文描述的一個(gè)或多個(gè)實(shí)施例包括粗化包括儲(chǔ)層模擬模型的網(wǎng)格的方法。根據(jù)一個(gè)實(shí)施例的特征，下面描述一種粗化網(wǎng)格的方法，其中指定對(duì)塊或單元的期望編組并創(chuàng)建各個(gè)單元(在本文中稱其為“非標(biāo)準(zhǔn)”形狀單元或塊)，這些單元在任一平面內(nèi)可潛在地具有四個(gè)以上的邊。計(jì)算經(jīng)編組的單元的孔隙體積，因?yàn)檫@是所有經(jīng)編組的單元之間的傳導(dǎo)能力。為了本說(shuō)明書目的，在任一平面內(nèi)具有多于或小于四個(gè)邊的各個(gè)單元在本文中被稱為“非標(biāo)準(zhǔn)”單元。特定非標(biāo)準(zhǔn)單元也被稱為“鋸齒形”單元，意思是這些非標(biāo)準(zhǔn)單元在任一平面內(nèi)的所有邊取向?yàn)楸舜嘶敬怪?。?duì)于嚴(yán)格矩形的塊，這些單元彼此以90°角相交。
[0035]圖8A是適于實(shí)現(xiàn)本文所述的發(fā)明的儲(chǔ)層模擬系統(tǒng)的示例性計(jì)算機(jī)系統(tǒng)800的框圖。在一個(gè)實(shí)施例中，計(jì)算機(jī)系統(tǒng)800包括至少一個(gè)處理器802、存儲(chǔ)器804、任選的I/O設(shè)備806、以及任選的顯示器808，這些部件全部經(jīng)由系統(tǒng)總線809互連?？捎商幚砥?02執(zhí)行的用于實(shí)現(xiàn)根據(jù)本文描述的實(shí)施例的儲(chǔ)層模擬系統(tǒng)810的軟件指令，可被存儲(chǔ)在存儲(chǔ)器804中。盡管圖8A中未明示，但可以理解計(jì)算機(jī)系統(tǒng)800可經(jīng)由合適的網(wǎng)絡(luò)連接連接至一個(gè)或多個(gè)公共和/或私人網(wǎng)絡(luò)。還要理解，包含儲(chǔ)層模擬系統(tǒng)810的軟件指令可從⑶-ROM或其它合適的存儲(chǔ)介質(zhì)被加載到存儲(chǔ)器804。
[0036]在本發(fā)明的一個(gè)實(shí)施例中，使用本領(lǐng)域已知的儲(chǔ)層模擬軟件來(lái)實(shí)現(xiàn)儲(chǔ)層模擬系統(tǒng)810的一部分。這樣的儲(chǔ)層模擬軟件一般利用儲(chǔ)層的數(shù)值表示，因?yàn)樗粯?gòu)想為在任何井被鉆探之前并在任何油田開(kāi)發(fā)前就存在。儲(chǔ)層的這種表示結(jié)合與建議的或已存在的井和開(kāi)發(fā)策略有關(guān)的附加數(shù)據(jù)，允許該軟件預(yù)測(cè)儲(chǔ)層就流體噴射和生產(chǎn)方面如何執(zhí)行。在現(xiàn)有技術(shù)中，這樣的儲(chǔ)層模擬軟件單純將網(wǎng)格粗化成標(biāo)準(zhǔn)形狀單元并對(duì)標(biāo)準(zhǔn)形狀單元引入單元孔隙率、單元深度、和傳導(dǎo)能力數(shù)據(jù)。同等地看待具有和不具有間斷點(diǎn)的粗化單元，這導(dǎo)致不準(zhǔn)確地反映間斷點(diǎn)的傳導(dǎo)能力。本發(fā)明的一個(gè)目的是允許創(chuàng)建鋸齒形塊并更準(zhǔn)確地確定受間斷點(diǎn)約束的那些單元的傳導(dǎo)能力，將更準(zhǔn)確的傳導(dǎo)能力數(shù)據(jù)提供給儲(chǔ)層模擬器。
[0037]在一個(gè)實(shí)施例中，這樣的儲(chǔ)層模擬軟件有效地模擬非結(jié)構(gòu)化網(wǎng)格的能力可被影響(leverage),并且Mcoarse陣列的變型可結(jié)合本發(fā)明使用用于在完全非結(jié)構(gòu)化的網(wǎng)格上進(jìn)行粗糙化的塊(標(biāo)準(zhǔn)塊和非標(biāo)準(zhǔn)塊兩者)的網(wǎng)格塊映射或編組，如下文中更詳細(xì)描述的那樣。在下文中更詳細(xì)描述的另一實(shí)施例中，塊被粗化并且沒(méi)有鋸齒形塊的區(qū)域被映射到結(jié)構(gòu)化網(wǎng)格上，而鋸齒形塊被局部地實(shí)現(xiàn)為非結(jié)構(gòu)化區(qū)域。配合前述實(shí)施例的方法被涵蓋在系統(tǒng)810中。
[0038]轉(zhuǎn)向圖8B，示出由系統(tǒng)810采用的方法。在步驟812中，創(chuàng)建相對(duì)較高分辨率的地層模型，其中該模型由多種多樣的標(biāo)準(zhǔn)形狀的精細(xì)單元或塊表征。典型地，此時(shí)，標(biāo)準(zhǔn)形狀的精細(xì)塊被實(shí)現(xiàn)在結(jié)構(gòu)化網(wǎng)格上。通過(guò)定義的精細(xì)網(wǎng)格塊，標(biāo)識(shí)精細(xì)網(wǎng)格塊及其Btt鄰的精細(xì)網(wǎng)格塊之間的每個(gè)連接。對(duì)于一對(duì)精細(xì)網(wǎng)格塊之間的每個(gè)所標(biāo)識(shí)的連接，計(jì)算這對(duì)網(wǎng)格塊的傳導(dǎo)能力。本領(lǐng)域內(nèi)技術(shù)人員將了解，可用業(yè)內(nèi)已知的數(shù)量眾多的方式來(lái)確定這樣的標(biāo)準(zhǔn)形狀的精細(xì)單元對(duì)之間的傳導(dǎo)能力。所標(biāo)識(shí)的連接和所計(jì)算的傳導(dǎo)能力被保留存以在后面的步驟812中使用。
[0039]在步驟820，使用合適的粗化因子(例如2X2、4X4等)對(duì)步驟812中定義的精細(xì)單元網(wǎng)格進(jìn)行粗化。在一個(gè)實(shí)施例中，如下執(zhí)行該步驟。開(kāi)始于圖2A所示的網(wǎng)格200，并根據(jù)本文描述的實(shí)施例應(yīng)用2X2粗化，導(dǎo)致如圖9A所示的網(wǎng)格900。特定地，在實(shí)現(xiàn)方法中，應(yīng)用2X2粗化并隨后將包括多于一個(gè)結(jié)構(gòu)化單元的任何所得粗糙塊(例如在其內(nèi)部具有一個(gè)或多個(gè)間斷點(diǎn)的粗糙塊)進(jìn)一步分割以最終使粗化網(wǎng)格的每個(gè)塊僅包括一個(gè)結(jié)構(gòu)化單元。例如，參見(jiàn)圖9A和圖9B，不修正本發(fā)明地應(yīng)用2X2粗化將導(dǎo)致塊3和塊4 一起被組合在單個(gè)粗糙塊(被稱為組合塊3/4)中并導(dǎo)致塊7和塊8 一起被組合在單個(gè)粗糙塊(被稱為組合塊7/8)中。然而，根據(jù)本發(fā)明的實(shí)施例，在步驟820中，組合塊3/4沿貫穿其的間斷點(diǎn)被進(jìn)一步分割成單獨(dú)塊3、4，優(yōu)選地利用沿該間斷點(diǎn)存在的之前定義的精細(xì)塊。對(duì)于組合塊7/8執(zhí)行相同的進(jìn)程。該步驟導(dǎo)致諸如塊3和塊8之類的非標(biāo)準(zhǔn)塊或單元的形成。在圖9B示出的情形下，塊3和塊8是鋸齒形單元。另外，該結(jié)果是完全“非結(jié)構(gòu)化的”網(wǎng)格，因?yàn)樵谡麄€(gè)網(wǎng)格中，單元不具有統(tǒng)一的尺寸和形狀。
[0040]本領(lǐng)域內(nèi)技術(shù)人員將理解，盡管本發(fā)明為方便起見(jiàn)被描述為存在于一個(gè)平面內(nèi)的正方形或矩形的單元或塊，然而本發(fā)明不僅限于定義在一個(gè)平面內(nèi)的網(wǎng)格，而是可等同地適用于圖2B所示的網(wǎng)格。因此，這些塊不一定是具有直線邊并以90°相交的嚴(yán)格矩形。例如，圖2B邏輯上與圖2A相同。在圖2B的情形下，每個(gè)邊不是直線，塊具有不同的尺寸，它們不以90°角相交并且網(wǎng)格是旋轉(zhuǎn)的。然而，它仍然是其中貫穿有間斷點(diǎn)的8X4網(wǎng)格。每個(gè)塊在2d中具有4條邊(在3D中具有6條邊)并在2d中具有4個(gè)角(3d中8個(gè)角)。因此，本發(fā)明的方法仍然適用。
[0041]作為圖9B所示的實(shí)施例的可選實(shí)施例，如圖9C所示并表示步驟820的另一實(shí)施例，可定義一種基本結(jié)構(gòu)化的網(wǎng)格，這種網(wǎng)格僅具有局部的非結(jié)構(gòu)化例外。在圖9C中，執(zhí)行與圖9B所示相同的粗化方法，但它被表示為主要結(jié)構(gòu)化的網(wǎng)格，由結(jié)構(gòu)化1、J索引來(lái)表征。包含間斷點(diǎn)或結(jié)構(gòu)化單元的塊，即塊(1=3，J=I)和(1=2，J=2)，使用下文中描述的鋸齒方法將其對(duì)待為局部非結(jié)構(gòu)化例外。
[0042]已發(fā)現(xiàn)較老式的儲(chǔ)層模擬器無(wú)法非常有效地處理非結(jié)構(gòu)化模型。在其中利用較老式的儲(chǔ)層模擬器的情形下，優(yōu)選的是利用圖9C所示本發(fā)明的主要結(jié)構(gòu)化/局部非結(jié)構(gòu)化的實(shí)施例。相反，傾向于更有效地起作用以解決問(wèn)題的較新儲(chǔ)層模擬器，優(yōu)選地用于本發(fā)明的完全非結(jié)構(gòu)化實(shí)施例。
[0043]本領(lǐng)域內(nèi)技術(shù)人員將理解，例如前面描述的那些現(xiàn)有技術(shù)粗化方法，要求所有所得網(wǎng)格塊被表征為在任一平面內(nèi)恰好具有四條邊，稱為“標(biāo)準(zhǔn)”形狀單元，并且在非斷層的毗鄰單元之間具有單個(gè)連接。直接相反，并根據(jù)本文描述的實(shí)施例的各個(gè)特征，由系統(tǒng)810放寬假設(shè)，允許在任一平面內(nèi)可具有非四條邊的單元，即“鋸齒形塊”或“非標(biāo)準(zhǔn)”形狀單元，以及在這些單元之間的多個(gè)連接。結(jié)果，不管如前所述修正的粗化因子生成什么所得塊，不管是標(biāo)準(zhǔn)形狀的還是非標(biāo)準(zhǔn)形狀的，它們將被視為粗糙塊，如圖9B和/或圖9C所示。特定地，圖9B所示的網(wǎng)格900包括十個(gè)隨意編號(hào)的塊1-10。要注意塊3和塊8是L形的。在圖9B中，盡管不需要特定的編號(hào)方案，然而為方便起見(jiàn)，使用沿X方向然后沿y方向增加的連續(xù)號(hào)碼；然而，可使用任何合適的編號(hào)方案，例如與Mcoarse陣列一致的編號(hào)方案,只要結(jié)果是10個(gè)粗糙塊即可。同樣，在圖9C中，使用主要結(jié)構(gòu)化的實(shí)現(xiàn)仍然得到10個(gè)單元，但在這種情形下，可使用結(jié)構(gòu)化1、J索引來(lái)表征8個(gè)單元并作為非標(biāo)準(zhǔn)單元和連接添加另外兩個(gè)單元。
[0044]圖10示出被外推至更區(qū)域以更清楚地展示模擬模型中的間斷點(diǎn)的鋸齒形本質(zhì)的粗糙塊。該模擬問(wèn)題現(xiàn)在是“非結(jié)構(gòu)化的”，從這個(gè)意義上說(shuō)網(wǎng)格塊索引的標(biāo)準(zhǔn)模擬1、J表示不再是全局有效的，因?yàn)樵撃Ｐ驮诿啃泻?或列中具有不同數(shù)量的單元。使用圖9B中的完全非結(jié)構(gòu)化表示或圖9C中的局部非結(jié)構(gòu)化表示解決了這個(gè)問(wèn)題。本文描述的實(shí)施例——且特定地是系統(tǒng)810——的一個(gè)目的是自動(dòng)地計(jì)算每個(gè)網(wǎng)格塊(包括鋸齒形網(wǎng)格塊)的孔隙體積、連接、以及它們的傳導(dǎo)能力。
[0045]再次參見(jiàn)圖SB，在步驟822中，通過(guò)對(duì)在粗化以形成粗糙塊的期間編組在一起的精細(xì)塊的孔隙體積求和來(lái)確定每個(gè)粗糙塊的孔隙體積。對(duì)于鋸齒和非鋸齒形塊兩者都是如此。可使用許多已知方法容易地執(zhí)行精細(xì)塊的孔隙體積的計(jì)算，并因此不再予以詳細(xì)描述。
[0046]在步驟824中，計(jì)算全部粗糙塊的傳導(dǎo)能力(不管這些粗糙塊是由標(biāo)準(zhǔn)形狀單元還是非標(biāo)準(zhǔn)形狀單元構(gòu)成)。如本文所述，本發(fā)明的新穎方面之一是在儲(chǔ)層建模過(guò)程中使用非標(biāo)準(zhǔn)形狀單元。一旦已定義非標(biāo)準(zhǔn)形狀單元，期望該單元和其毗鄰單元之間的傳導(dǎo)能力的準(zhǔn)確計(jì)算。下面相對(duì)于圖13、圖14更詳細(xì)地討論用于計(jì)算粗化的非標(biāo)準(zhǔn)形狀單元的傳導(dǎo)能力的方法的非限制性示例。
[0047]在任何情況下，在步驟826中，計(jì)算每個(gè)粗糙塊的平均深度。在步驟828中，通過(guò)將井從精細(xì)參考坐標(biāo)系轉(zhuǎn)化至粗糙參考坐標(biāo)系，必要時(shí)維持所需的那些精細(xì)單元連接至粗糙單元計(jì)算要求精細(xì)單元對(duì)性質(zhì)的程度，來(lái)將井從精細(xì)網(wǎng)格映射到粗糙網(wǎng)格。已發(fā)現(xiàn)執(zhí)行這種映射的一種優(yōu)選方法是利用Mcoarse陣列。特定地，本領(lǐng)域內(nèi)技術(shù)人員將理解，通過(guò)首先標(biāo)識(shí)所有精細(xì)網(wǎng)格連接來(lái)完成映射。如所述，這些單元被粗化。在進(jìn)行此舉時(shí)，創(chuàng)建Mcoarse陣列，該Mcoarse陣列標(biāo)識(shí)每個(gè)精細(xì)單元將被編組到的粗糙單元。然后在逐個(gè)精細(xì)單元的基礎(chǔ)上，精細(xì)單元對(duì)被組合在一起而消除連接，即，令具有相同Mcoarse值的形成標(biāo)準(zhǔn)形狀的粗糙單元。此舉繼續(xù)，直到處理了所有精細(xì)單元連接對(duì)為止。剩下的是標(biāo)準(zhǔn)形狀粗糙單元之間的連接。使用Mcoarse編組陣列來(lái)完成精細(xì)單元編號(hào)和粗糙單元編號(hào)之間的映射。
[0048]最后，在步驟830中，利用合適的儲(chǔ)層模擬器來(lái)寫非結(jié)構(gòu)化模擬數(shù)據(jù)文件。迄今為止，非標(biāo)準(zhǔn)形狀單元尚未被用于如前所述和圖8B所示的儲(chǔ)層模型中。
[0049]下面將更詳細(xì)地描述粗化單元的傳導(dǎo)能力計(jì)算的示例性實(shí)施例。在圖9中，每個(gè)網(wǎng)格塊是由該網(wǎng)格塊和其毗鄰網(wǎng)格塊之間的連接來(lái)表征。這些“連接”由圖11中的箭頭表示。在圖11中，X和y表示連接強(qiáng)度或連接的傳導(dǎo)能力。同樣，由這樣一個(gè)事實(shí)，即粗糙塊3和粗糙塊2、2、4、4、5和9之間存在連接，來(lái)說(shuō)明非結(jié)構(gòu)化本質(zhì)。塊4在該列表中出現(xiàn)兩次的事實(shí)不是錯(cuò)誤，因?yàn)樵趬K3和4之間存在兩個(gè)連接。對(duì)于塊2和塊3之間的連接也是如此。如前面步驟812中提到的并為了促進(jìn)底層計(jì)算，檢測(cè)出精細(xì)塊之間的所有連接。如前面提到的，為了執(zhí)行儲(chǔ)層模擬，儲(chǔ)層模擬軟件所要求的主要數(shù)據(jù)關(guān)聯(lián)于單元(網(wǎng)格塊)、單元尺寸(孔隙體積)、單元之間的連接強(qiáng)度(傳導(dǎo)能力)、以及每個(gè)單元的平均深度。本發(fā)明提供了用于計(jì)算更準(zhǔn)確的傳導(dǎo)能力數(shù)據(jù)的方法，其中可對(duì)與一對(duì)單元相關(guān)聯(lián)的個(gè)別連接傳導(dǎo)能力求和來(lái)推導(dǎo)出共享至少一個(gè)連接的任何兩個(gè)單元之間的合計(jì)的傳導(dǎo)能力。
[0050]例如，在圖11中，粗化單元1、2是標(biāo)準(zhǔn)形狀的單元。因此，單元1、2之間僅存在單個(gè)連接并且單元1、2之間的連接僅存在單個(gè)傳導(dǎo)能力。因此可使用業(yè)內(nèi)已知的傳導(dǎo)能力計(jì)算的標(biāo)準(zhǔn)方法來(lái)計(jì)算出單元1、2的傳導(dǎo)能力。另一方面，粗化單元3(以及單元8)表示非標(biāo)準(zhǔn)形狀單元。因此，本發(fā)明的方法必須被應(yīng)用于粗化單元3以推導(dǎo)出粗化單元3和其每個(gè)毗鄰單元之間的合計(jì)的傳導(dǎo)能力，以供儲(chǔ)層模擬軟件使用。
[0051]在考慮單元3時(shí)，可看到單元3以單元2、5、4、和9為邊界。由于單元3是非標(biāo)準(zhǔn)形狀單元，為了確定單元3的合計(jì)的傳導(dǎo)能力，基于之前包括單元3的精細(xì)單元來(lái)操縱單元
3。由此，在本發(fā)明過(guò)程中的這點(diǎn)處，利用之前在步驟812確定和留存的精細(xì)單元連接和傳導(dǎo)能力數(shù)據(jù)。將包括粗糙單元3的精細(xì)單元考慮在內(nèi)，可以看到，單元3及其毗鄰單元之間總共存在六個(gè)連接。此外，每個(gè)連接具有精細(xì)單元傳導(dǎo)能力，從而xl_x4是沿X方向在單元3和毗鄰單元之間的連接的傳導(dǎo)能力，且yl_y2是沿y方向在單元3和毗鄰單元之間的連接的傳導(dǎo)能力。如通過(guò)單元2、3之間存在的兩個(gè)連接可看到的那樣，由于標(biāo)識(shí)非標(biāo)準(zhǔn)單元的連接和傳導(dǎo)能力目的而使用精細(xì)單元，鄰接標(biāo)準(zhǔn)單元和非標(biāo)準(zhǔn)單元的單邊之間可存在多于一個(gè)的連接。同樣，即便兩單元之間僅存在單邊，但是兩個(gè)傳導(dǎo)能力——即xl、x2——表征單元2、3之間的總連接。在其它情形下，兩個(gè)單元之間可存在多個(gè)邊，就像單元3和單元4之間的情形。在這些情形下，可由至少一個(gè)連接(也可能是多個(gè)連接)來(lái)表征每個(gè)邊。在圖11的示例中，可在單元3和單元4之間定義兩個(gè)連接。由于這些連接中的每一個(gè)具有與之關(guān)聯(lián)的傳導(dǎo)能力(如步驟812中確定的)，因此兩個(gè)傳導(dǎo)能力——即x4和y2——表征單元3、4之間的總連接。
[0052]一旦非標(biāo)準(zhǔn)單元及其毗鄰單元之間的每個(gè)連接和關(guān)聯(lián)的傳導(dǎo)能力被計(jì)算，可通過(guò)合計(jì)這對(duì)單元之間的傳導(dǎo)能力來(lái)確定非標(biāo)準(zhǔn)單元和每個(gè)毗鄰單元之間的總傳導(dǎo)能力。在該示例中，單元3、4之間的合計(jì)的傳導(dǎo)能力被定義為x4+y2。類似地，單元3、2之間的合計(jì)的傳導(dǎo)能力為xl+x2,單元3、5之間的合計(jì)的傳導(dǎo)能力為x3,而單元3、9之間的合計(jì)的傳導(dǎo)能力為yl。
[0053]維持在其間僅具有單邊的毗鄰單元之間的精細(xì)標(biāo)度連接的一個(gè)理由在于便于合計(jì)粗化的、非標(biāo)準(zhǔn)單元之間的傳導(dǎo)能力。在本發(fā)明的一個(gè)實(shí)施例中，可利用平行管道計(jì)算的求和來(lái)確定合計(jì)的傳導(dǎo)能力，其中可利用平行管道對(duì)粗糙塊進(jìn)行建模。本發(fā)明的一個(gè)新穎方面在于，利用該特定方法來(lái)確定非標(biāo)準(zhǔn)形狀單元的傳導(dǎo)能力，利用標(biāo)準(zhǔn)形狀(即矩形)的管道。每個(gè)管道可具有不同的長(zhǎng)度和性質(zhì)(傳導(dǎo)能力、滲透性等)，由此，必須分開(kāi)計(jì)算管道內(nèi)的每個(gè)傳導(dǎo)能力。
[0054]參見(jiàn)圖12，式子第一塊1201具有滲透性K1、長(zhǎng)度L1、高度H、寬度W、并對(duì)其施加壓力P1，而第二塊1202具有滲透性K2、長(zhǎng)度L2、高度H、寬度W、并對(duì)其施加壓力P2，第一塊1201和第二塊1202之間的傳導(dǎo)能力(“T12”)由下列式子定義:
[0055]T12=2HW/ ((Ll/Kl) + (L2/K2))
[0056]存在許多方法來(lái)計(jì)算本文描述的非標(biāo)準(zhǔn)單元或鋸齒形塊的傳導(dǎo)能力和滲透性。兩個(gè)這樣的方法包括:1)對(duì)傳導(dǎo)能力求和并從中計(jì)算有效滲透性(圖13所示)；以及(2)對(duì)每個(gè)塊中的滲透性求平均并隨后計(jì)算傳導(dǎo)能力(圖14所示)。
[0057]參見(jiàn)圖13，在步驟1300中，計(jì)算每個(gè)精細(xì)網(wǎng)格連接的傳導(dǎo)能力。在1302中，使用構(gòu)成粗糙網(wǎng)格塊的精細(xì)網(wǎng)格塊的精細(xì)網(wǎng)格傳導(dǎo)能力并考慮到連接的每個(gè)行和列的尺寸可不同，來(lái)計(jì)算每個(gè)粗糙網(wǎng)格連接的每邊的“半傳導(dǎo)能力”或“部分傳導(dǎo)能力”。本領(lǐng)域內(nèi)技術(shù)人員將理解，基于前面的式子和圖12，與兩單元之間的連接相關(guān)聯(lián)的傳導(dǎo)能力取決于所連接的兩單元的性質(zhì)。半傳導(dǎo)能力表示與單元關(guān)聯(lián)的傳導(dǎo)能力類屬性，而不是其本身的總連接。在該示例中，Txl和Tx2是分別被定義為單元1、2的半傳導(dǎo)能力，其中Txi= (KiHW)/Li (i=1或2)。在這種情形下，T12可被寫成Τ12=2/((1/Τχ1) + (1/Τχ2))。因此，連接屬性可被表示為個(gè)別塊的半傳導(dǎo)能力的調(diào)和平均數(shù)。
[0058]接著，在步驟1303中，計(jì)算每對(duì)標(biāo)準(zhǔn)形狀的粗化單元之間存在的每個(gè)粗糙網(wǎng)格連接的傳導(dǎo)能力。在該步驟中，本質(zhì)上說(shuō)，精細(xì)網(wǎng)格傳導(dǎo)能力被折疊成(collapsed into)這些單元的粗糙網(wǎng)格傳導(dǎo)能力。在步驟1304中，精細(xì)網(wǎng)格連接傳導(dǎo)能力，例如非標(biāo)準(zhǔn)形狀單元的半傳導(dǎo)能力，被求和以獲得粗糙網(wǎng)格連接的合計(jì)的傳導(dǎo)能力。盡管可利用對(duì)這些半傳導(dǎo)能力求和的各種方法，然而執(zhí)行這些求和的一個(gè)這樣的方法是一些管法(tubes method)。使用管方的求和，基于被合適地加權(quán)的精細(xì)標(biāo)度的半傳導(dǎo)能力的調(diào)和平均數(shù)，來(lái)計(jì)算粗糙傳導(dǎo)能力，不管是總傳導(dǎo)能力還是部分或半傳導(dǎo)能力。這非常類似于圖12所示的內(nèi)容，除了在每個(gè)粗糙塊(例如1201)可由多個(gè)精細(xì)塊構(gòu)成的情況外。然而，調(diào)和求平均是相同的。所以使用精細(xì)標(biāo)度的半傳導(dǎo)能力來(lái)計(jì)算粗糙標(biāo)度的半傳導(dǎo)能力。然后使用這些粗糙標(biāo)度的半傳導(dǎo)能力來(lái)計(jì)算總粗糙標(biāo)度傳導(dǎo)能力。如果在每個(gè)粗糙單元對(duì)之間存在多于一個(gè)的連接，則對(duì)傳導(dǎo)能力進(jìn)行求和。
[0059]最終，在步驟1306中，從合計(jì)的傳導(dǎo)能力來(lái)倒推計(jì)算有效滲透性。
[0060]圖13所不的方法的一個(gè)可選方法不出于圖14,其中不使用半傳導(dǎo)能力。相反,通過(guò)對(duì)每個(gè)粗糙塊中的滲透性求平均并使用前面示出的T12的式子，來(lái)計(jì)算粗糙標(biāo)度的傳導(dǎo)能力。
[0061]特定地，參見(jiàn)圖14，在步驟1400中，使用業(yè)內(nèi)已知的用戶指定的求平均技術(shù)，對(duì)精細(xì)網(wǎng)格滲透性求平均，諸如算術(shù)、調(diào)和、冪定律等。在步驟1402中，考慮尺寸可不同，計(jì)算每個(gè)塊內(nèi)的每個(gè)行和列的塊間傳導(dǎo)能力。最終，在步驟1404中，精細(xì)網(wǎng)格連接傳導(dǎo)能力被求和以獲得合計(jì)的塊間傳導(dǎo)能力。
[0062]本文描述的實(shí)施例，即便模型具有間斷點(diǎn)，允許真正可縮放的儲(chǔ)層模擬解決方案。在不改變底層模型的情況下，可自動(dòng)地改變模擬表示的保真度。高分辨率模型可完全應(yīng)用至(be taken all the way to)材料平衡模型或兩者間的任何,而不需要知道關(guān)于儲(chǔ)層模擬、粗化、或提升的任何信息。
[0063]本發(fā)明的一個(gè)實(shí)施例是粗化包括多個(gè)精細(xì)網(wǎng)格塊的精細(xì)網(wǎng)格的計(jì)算機(jī)實(shí)現(xiàn)的方法，該精細(xì)網(wǎng)格表示其中具有至少一個(gè)間斷點(diǎn)的地理模型。該方法包括:將若干精細(xì)網(wǎng)格塊編組在一起以形成粗糙網(wǎng)格塊，其中粗糙網(wǎng)格塊中的至少一個(gè)是非標(biāo)準(zhǔn)形狀網(wǎng)格塊；并計(jì)算每對(duì)毗鄰的粗糙網(wǎng)格塊的傳導(dǎo)能力。該計(jì)算包括:計(jì)算每對(duì)毗鄰精細(xì)網(wǎng)格塊的傳導(dǎo)能力；將每個(gè)精細(xì)網(wǎng)格塊對(duì)映射到粗糙網(wǎng)格塊對(duì)；并對(duì)于其中網(wǎng)格塊之一包括非標(biāo)準(zhǔn)形狀的網(wǎng)格塊的每個(gè)網(wǎng)格塊對(duì)，對(duì)該精細(xì)網(wǎng)格塊對(duì)的傳導(dǎo)能力求和。在一個(gè)實(shí)施例中，對(duì)精細(xì)網(wǎng)格塊值的求和涉及使用被映射到粗糙網(wǎng)格內(nèi)的所有精細(xì)塊的精細(xì)網(wǎng)格傳導(dǎo)能力和/或滲透性、以及網(wǎng)格塊尺寸和網(wǎng)格塊重疊性，進(jìn)行計(jì)算。
[0064]另一實(shí)施例包括計(jì)算機(jī)實(shí)現(xiàn)的儲(chǔ)層模擬系統(tǒng)。該系統(tǒng)包括:處理器；可由處理器訪問(wèn)的存儲(chǔ)介質(zhì)；以及存儲(chǔ)在存儲(chǔ)介質(zhì)上的軟件指令。軟件指令可由處理器執(zhí)行用于對(duì)包括多個(gè)精細(xì)網(wǎng)格塊的精細(xì)網(wǎng)格模型進(jìn)行粗化，其中粗化包括根據(jù)粗化因子將精細(xì)網(wǎng)格塊編組為包括粗糙網(wǎng)格模型的粗糙網(wǎng)格塊，并且粗糙網(wǎng)格塊中的至少一個(gè)具有非標(biāo)準(zhǔn)形狀；并計(jì)算每對(duì)毗鄰粗糙網(wǎng)格塊的傳導(dǎo)能力。
[0065]又一實(shí)施例是用于粗化包括多個(gè)精細(xì)網(wǎng)格塊的精細(xì)網(wǎng)格的計(jì)算機(jī)實(shí)現(xiàn)的系統(tǒng)，該精細(xì)網(wǎng)格表示其中具有至少一個(gè)間斷點(diǎn)的地理模型。系統(tǒng)包括處理器和可由處理器訪問(wèn)的存儲(chǔ)介質(zhì)。該系統(tǒng)進(jìn)一步包括存儲(chǔ)在存儲(chǔ)介質(zhì)上并可由處理器執(zhí)行用于將若干精細(xì)網(wǎng)格塊編組在一起以形成粗糙網(wǎng)格塊的指令，其中至少一個(gè)粗糙網(wǎng)格塊是非標(biāo)準(zhǔn)形狀的網(wǎng)格塊，該編組包括:根據(jù)用戶指定的粗化因子來(lái)粗化精細(xì)網(wǎng)格；以及計(jì)算其中粗糙網(wǎng)格塊對(duì)的至少一個(gè)網(wǎng)格塊是非標(biāo)準(zhǔn)形狀的網(wǎng)格塊的每對(duì)毗鄰粗糙網(wǎng)格塊的傳導(dǎo)能力。該計(jì)算包括:計(jì)算每對(duì)Btt鄰精細(xì)網(wǎng)格塊的傳導(dǎo)能力；將每個(gè)精細(xì)網(wǎng)格塊對(duì)映射到粗糖網(wǎng)格塊對(duì)；并對(duì)于其中網(wǎng)格塊之一是非標(biāo)準(zhǔn)形狀的每個(gè)網(wǎng)格塊對(duì)，對(duì)被映射至粗糙網(wǎng)格塊對(duì)的精細(xì)網(wǎng)格塊對(duì)的傳導(dǎo)能力進(jìn)行求和。
[0066]為了描述目的，術(shù)語(yǔ)“標(biāo)準(zhǔn)”網(wǎng)格塊和“非標(biāo)準(zhǔn)”網(wǎng)格塊的使用是指當(dāng)與網(wǎng)格中的多個(gè)塊的常見(jiàn)形狀比較時(shí)感興趣網(wǎng)格塊的形狀。例如，網(wǎng)格一般可由具有特定尺寸W、h和I的多個(gè)矩形網(wǎng)格塊來(lái)表征，藉此表示網(wǎng)格的塊的“標(biāo)準(zhǔn)”形狀。典型地，模擬單元是盒狀形狀的，其具有8個(gè)角和6個(gè)面。本發(fā)明的方法放寬這種假設(shè)以利用存在內(nèi)部特征或間斷點(diǎn)的非傳統(tǒng)形狀。由此，這種情形下的非標(biāo)準(zhǔn)網(wǎng)格塊將是非正方形的網(wǎng)格塊，諸如是鋸齒形網(wǎng)格塊、或具有與多個(gè)網(wǎng)格塊不同形狀的網(wǎng)格塊。本領(lǐng)域內(nèi)技術(shù)人員將理解，盡管傳統(tǒng)的矩形單元可相對(duì)于在本發(fā)明一個(gè)實(shí)施例中利用的非標(biāo)準(zhǔn)單元而被指定為基準(zhǔn)，然而在本發(fā)明的另一實(shí)施例中，“標(biāo)準(zhǔn)”單元可具有不同的形狀，只要總網(wǎng)格的多個(gè)單元也具有該形狀即可。例如，“標(biāo)準(zhǔn)”網(wǎng)格塊可由多種多樣的三角形單元構(gòu)成。在這種情形下，非標(biāo)準(zhǔn)單元將是具有非三角形的形狀的單元。
[0067]另外，在本發(fā)明的討論和說(shuō)明中，為了不過(guò)分使討論變得復(fù)雜，使用沿真實(shí)的X和y方向兩者使用相同程度粗化的恒定粗化因子。然而，本領(lǐng)域內(nèi)技術(shù)人員將理解，本發(fā)明適用于沿任何方向的任何粗化因子。它還支持在模型的不同部分中使用不同粗化因子。
[0068]盡管本文已詳細(xì)描述了本發(fā)明的某些特征和實(shí)施例，然而容易理解，本發(fā)明涵蓋落在下面權(quán)利要求書的范圍和精神內(nèi)的所有修正和改進(jìn)。此外，不旨在對(duì)本文所示的結(jié)構(gòu)或設(shè)計(jì)的細(xì)節(jié)作出限制，除了如下面權(quán)利要求書所述。另外，本領(lǐng)域內(nèi)技術(shù)人員將理解，各組件被表述為垂直取向或水平取向不旨在作為限制，而是為了方便描述本發(fā)明而設(shè)置。
[0069]因此，顯然，上述公開(kāi)的特定解說(shuō)性實(shí)施例可被修改或修正，并且所有這些變型都被認(rèn)為落在本發(fā)明的范圍和精神內(nèi)。另外，權(quán)利要求書中的術(shù)語(yǔ)具有其明白的、常見(jiàn)的含義，除非由專利權(quán)人另外明白且清楚地定義了其它。
【權(quán)利要求】
1.一種粗化包括多個(gè)精細(xì)網(wǎng)格塊的精細(xì)網(wǎng)格的計(jì)算機(jī)實(shí)現(xiàn)的方法，所述精細(xì)網(wǎng)格表示其中具有至少一個(gè)間斷點(diǎn)的地理模型，所述方法包括: 計(jì)算每對(duì)毗鄰精細(xì)網(wǎng)格塊的傳導(dǎo)能力； 將若干精細(xì)網(wǎng)格塊編組在一起以形成粗糙網(wǎng)格塊，其中所述粗糙網(wǎng)格塊中的至少一個(gè)是由多個(gè)精細(xì)網(wǎng)格塊表征的非標(biāo)準(zhǔn)形狀網(wǎng)格塊； 將非標(biāo)準(zhǔn)形狀網(wǎng)格塊的精細(xì)網(wǎng)格塊對(duì)中的每一對(duì)映射至粗糙網(wǎng)格塊對(duì)；以及計(jì)算每對(duì)毗鄰粗糙網(wǎng)格塊的傳導(dǎo)能力，其中所述粗糙網(wǎng)格塊對(duì)中的至少一個(gè)網(wǎng)格塊是非標(biāo)準(zhǔn)形狀網(wǎng)格塊。
2.如權(quán)利要求1所述的計(jì)算機(jī)實(shí)現(xiàn)的方法，其特征在于，所述編組是根據(jù)粗化因子執(zhí)行的。
3.如權(quán)利要求1所述的計(jì)算機(jī)實(shí)現(xiàn)的方法，其特征在于，所述編組還包括: 根據(jù)粗化因子來(lái)粗化所述精細(xì)網(wǎng)格以形成粗糙網(wǎng)格塊；以及 在所述粗化之后，將包括多于一個(gè)結(jié)構(gòu)化單元的粗糙網(wǎng)格塊中的那些細(xì)分以創(chuàng)建僅包括一個(gè)結(jié)構(gòu)化單元的粗糙網(wǎng)格塊。
4.如權(quán)利要求1所述的計(jì)算機(jī)實(shí)現(xiàn)的方法，其特征在于，還包括計(jì)算每對(duì)毗鄰粗糙網(wǎng)格塊的傳導(dǎo)能力，其中所述粗糙網(wǎng)格塊對(duì)中都不是非標(biāo)準(zhǔn)網(wǎng)格塊，所述計(jì)算包括: 計(jì)算每對(duì)毗鄰精細(xì)網(wǎng)格塊的傳導(dǎo)能力； 將所述精細(xì)網(wǎng)格塊對(duì)中的每一對(duì)映射至粗糙網(wǎng)格塊對(duì)；以及 對(duì)于僅具有標(biāo)準(zhǔn)形狀網(wǎng)格塊的每個(gè)網(wǎng)格塊對(duì)，計(jì)算粗糙標(biāo)度傳導(dǎo)能力。
5.如權(quán)利要求1所述的計(jì)算機(jī)實(shí)現(xiàn)的方法，其特征在于，所述計(jì)算粗糙網(wǎng)格塊傳導(dǎo)能力的步驟進(jìn)一步包括: 對(duì)包括非標(biāo)準(zhǔn)形狀網(wǎng)格塊的每個(gè)網(wǎng)格塊對(duì)，對(duì)非標(biāo)準(zhǔn)形狀網(wǎng)格塊的精細(xì)網(wǎng)格塊對(duì)的多個(gè)傳導(dǎo)能力求和，從而計(jì)算具有非標(biāo)準(zhǔn)形狀網(wǎng)格塊的粗糙網(wǎng)格塊對(duì)的總粗糙標(biāo)度。
6.如權(quán)利要求1所述的計(jì)算機(jī)實(shí)現(xiàn)的方法，其特征在于，所述非標(biāo)準(zhǔn)形狀網(wǎng)格塊是鋸齒形網(wǎng)格塊。
7.一種計(jì)算機(jī)實(shí)現(xiàn)的儲(chǔ)層模擬系統(tǒng)，所述系統(tǒng)包括: 處理器； 可由所述處理器訪問(wèn)的存儲(chǔ)介質(zhì)；以及 存儲(chǔ)在所述存儲(chǔ)介質(zhì)上并可由所述處理器執(zhí)行的軟件指令，所述處理器執(zhí)行所述軟件指令以: 粗化包括多個(gè)精細(xì)網(wǎng)格塊的精細(xì)網(wǎng)格模型，其中所述粗化包括根據(jù)粗化因子將所述精細(xì)網(wǎng)格塊編組成包括粗糙網(wǎng)格模型的粗糙網(wǎng)格塊，并且其中所述粗糙網(wǎng)格塊中的至少一個(gè)具有非標(biāo)準(zhǔn)形狀； 確定具有非標(biāo)準(zhǔn)形狀的所述粗糙網(wǎng)格塊中的每一個(gè)的孔隙體積；以及 計(jì)算每對(duì)毗鄰粗糙網(wǎng)格塊的總傳導(dǎo)能力。
8.如權(quán)利要求7所述的計(jì)算機(jī)實(shí)現(xiàn)的儲(chǔ)層模擬系統(tǒng)，其特征在于，還包括可由處理器執(zhí)行的用于計(jì)算每個(gè)粗糙塊的平均深度的軟件。
9.如權(quán)利要求7所述的計(jì)算機(jī)實(shí)現(xiàn)的儲(chǔ)層模擬系統(tǒng)，其特征在于，還包括可由所述處理器執(zhí)行用于將井從所述精細(xì)網(wǎng)格模型映射至所述粗糙網(wǎng)格模型的軟件。
10.如權(quán)利要求7所述的計(jì)算機(jī)實(shí)現(xiàn)的儲(chǔ)層模擬系統(tǒng)，其特征在于，還包括可由所述處理器執(zhí)行用于寫非結(jié)構(gòu)化數(shù)據(jù)文件的軟件。
11.如權(quán)利要求7所述的計(jì)算機(jī)實(shí)現(xiàn)的儲(chǔ)層模擬系統(tǒng)，其特征在于，所述計(jì)算傳導(dǎo)能力包括: 計(jì)算所述精細(xì)網(wǎng)格塊中每對(duì)毗鄰精細(xì)網(wǎng)格塊的傳導(dǎo)能力； 將所述精細(xì)網(wǎng)格塊對(duì)中的每一對(duì)映射至粗糙網(wǎng)格塊對(duì)；以及 對(duì)其中網(wǎng)格塊之一具有非標(biāo)準(zhǔn)形狀的每個(gè)網(wǎng)格塊對(duì)，對(duì)非標(biāo)準(zhǔn)形狀網(wǎng)格塊的精細(xì)網(wǎng)格塊對(duì)的多個(gè)傳導(dǎo)能力求和，從而計(jì)算所述粗糙網(wǎng)格塊對(duì)的總粗糙標(biāo)度的傳導(dǎo)能力。
12.如權(quán)利要求11所述的計(jì)算機(jī)實(shí)現(xiàn)的儲(chǔ)層模擬系統(tǒng)，其特征在于，所述計(jì)算傳導(dǎo)能力進(jìn)一步包括: 從所述計(jì)算出的傳導(dǎo)能力倒推計(jì)算有效滲透性。
13.如權(quán)利要求7所述的計(jì)算機(jī)實(shí)現(xiàn)的儲(chǔ)層模擬系統(tǒng)，其特征在于，所述計(jì)算傳導(dǎo)能力包括: 使用用戶指定的技術(shù)來(lái)對(duì)精細(xì)網(wǎng)格塊的滲透性求平均； 計(jì)算每對(duì)毗鄰精細(xì)網(wǎng)格塊的傳導(dǎo)能力； 將所述精細(xì)網(wǎng)格塊對(duì)中的每一對(duì)映射至粗糙網(wǎng)格塊對(duì)； 對(duì)于僅具有標(biāo)準(zhǔn)形狀網(wǎng)格塊的每個(gè)網(wǎng)格塊對(duì)，計(jì)算粗糙標(biāo)度的傳導(dǎo)能力；以及在其中所述網(wǎng)格塊對(duì)之一具有非標(biāo)準(zhǔn)形狀的情形下，對(duì)非標(biāo)準(zhǔn)形狀網(wǎng)格塊的精細(xì)網(wǎng)格塊對(duì)的多個(gè)傳導(dǎo)能力求和，從而計(jì)算所述粗糙網(wǎng)格塊對(duì)的總粗糙標(biāo)度的傳導(dǎo)能力。
14.如權(quán)利要求7所述的計(jì)算機(jī)實(shí)現(xiàn)的儲(chǔ)層模擬系統(tǒng)，其特征在于，確定每個(gè)粗糙網(wǎng)格塊的孔隙體積包括對(duì)包括所述粗糙網(wǎng)格塊的精細(xì)網(wǎng)格塊的孔隙體積求和。
15.如權(quán)利要求7所述的計(jì)算機(jī)實(shí)現(xiàn)的儲(chǔ)層模擬系統(tǒng)，其特征在于，所述非標(biāo)準(zhǔn)形狀是鋸齒形。
16.一種用于粗化包括多個(gè)精細(xì)網(wǎng)格塊的精細(xì)網(wǎng)格的計(jì)算機(jī)實(shí)現(xiàn)的系統(tǒng)，所述精細(xì)網(wǎng)格表示其中具有至少一個(gè)間斷點(diǎn)的地理模型，所述系統(tǒng)包括: 處理器； 可由所述處理器訪問(wèn)的存儲(chǔ)介質(zhì)；以及存儲(chǔ)在所述存儲(chǔ)介質(zhì)上并可由所述處理器執(zhí)行的指令，所述處理器執(zhí)行所述指令以:將若干精細(xì)網(wǎng)格塊編組在一起以形成粗糙網(wǎng)格塊，其中所述粗糙網(wǎng)格塊中的至少一個(gè)是非標(biāo)準(zhǔn)形狀網(wǎng)格塊，所述編組包括根據(jù)用戶指定的粗化因子來(lái)粗化所述精細(xì)網(wǎng)格；以及計(jì)算每對(duì)毗鄰粗糙網(wǎng)格塊的傳導(dǎo)能力，其中所述粗糙網(wǎng)格塊對(duì)中的至少一個(gè)網(wǎng)格塊是非標(biāo)準(zhǔn)形狀網(wǎng)格塊，所述計(jì)算包括: 計(jì)算每對(duì)毗鄰精細(xì)網(wǎng)格塊的傳導(dǎo)能力； 將所述精細(xì)網(wǎng)格塊對(duì)中的每一對(duì)映射至粗糙網(wǎng)格塊對(duì)；以及 對(duì)于包括非標(biāo)準(zhǔn)形狀網(wǎng)格塊的每個(gè)網(wǎng)格塊對(duì)，計(jì)算合適的粗糙標(biāo)度的傳導(dǎo)能力，并對(duì)對(duì)應(yīng)于映射至此的原始精細(xì)網(wǎng)格對(duì)的粗糙標(biāo)度的傳導(dǎo)能力求和。
17.如權(quán)利要求16所述的系統(tǒng)，其特征在于，所述編組還包括:在所述粗化之后，將粗糙網(wǎng)格塊中包括多于一個(gè)結(jié) 構(gòu)化單元的那些細(xì)分來(lái)創(chuàng)建僅包括一個(gè)結(jié)構(gòu)化單元的粗糙網(wǎng)格塊。
18.如權(quán)利要求16所述的系統(tǒng)，其特征在于，還包括存儲(chǔ)在所述存儲(chǔ)介質(zhì)上并可由所述處理器執(zhí)行用于計(jì)算每對(duì)毗鄰的粗糙網(wǎng)格塊的傳導(dǎo)能力的指令，其中所述粗糙網(wǎng)格塊對(duì)中沒(méi)有非標(biāo)準(zhǔn)形狀網(wǎng)格塊，所述計(jì)算包括: 計(jì)算每對(duì)毗鄰精細(xì)網(wǎng)格塊的傳導(dǎo)能力； 將每個(gè)所述精細(xì)網(wǎng)格塊對(duì)中的每一個(gè)映射至粗糙網(wǎng)格塊對(duì)；以及 對(duì)于僅具有標(biāo)準(zhǔn)形狀網(wǎng)格塊的每個(gè)網(wǎng)格塊對(duì)，計(jì)算所述粗糙標(biāo)度的傳導(dǎo)能力。
19.如權(quán)利要求16所述的系統(tǒng)，其特征在于，還包括存儲(chǔ)在所述存儲(chǔ)介質(zhì)上并可由所述處理器執(zhí)行用于計(jì)算每對(duì)毗鄰的粗糙網(wǎng)格塊的傳導(dǎo)能力的指令，其中所述粗糙網(wǎng)格塊對(duì)中沒(méi)有網(wǎng)格塊是非標(biāo)準(zhǔn)形狀網(wǎng)格塊。
20.如權(quán)利要求16所述的系統(tǒng)，其特征在于，所述非標(biāo)準(zhǔn)形狀網(wǎng)格塊是鋸齒形網(wǎng)格塊。
21.一種用于粗化包括多個(gè)精細(xì)網(wǎng)格塊的精細(xì)網(wǎng)格從而對(duì)儲(chǔ)層建模的計(jì)算機(jī)實(shí)現(xiàn)的方法，所述精細(xì)網(wǎng)格表示其中具有至少一個(gè)間斷點(diǎn)的地理模型，所述方法包括: 界定精細(xì)單元的網(wǎng)格； 標(biāo)識(shí)通過(guò)所述網(wǎng)格的間斷點(diǎn)； 粗化所述網(wǎng)格以產(chǎn)生粗糙單元； 對(duì)沒(méi)有間斷點(diǎn)的每對(duì)粗糙單元，確定所述粗糙單元對(duì)的粗糙的傳導(dǎo)能力； 對(duì)于其中在所述粗糙單元的至少一個(gè)中存在間斷點(diǎn)的每對(duì)粗糙單元，將所述精細(xì)單元對(duì)中的每一對(duì)映射至粗糙網(wǎng)格塊對(duì)； 對(duì)包括具有間斷點(diǎn)的粗糙單元的每個(gè)精細(xì)單元計(jì)算半傳導(dǎo)能力和/或平均傳導(dǎo)能力；通過(guò)采用(i)對(duì)每個(gè)精細(xì)單元計(jì)算出的半傳導(dǎo)能力和/或滲透性、(ii)所述精細(xì)單元對(duì)的映射位置，來(lái)計(jì)算其中在至少一個(gè)粗糙單元中存在間斷點(diǎn)的每對(duì)粗糙單元的粗糙標(biāo)度的傳導(dǎo)能力；以及 對(duì)在相同單元之間具有多個(gè)連接的任何粗糙單元的粗糙標(biāo)度的傳導(dǎo)能力求和；以及 生成建模的儲(chǔ)層的2D或3D表示。
22.如權(quán)利要求21所述的計(jì)算機(jī)實(shí)現(xiàn)的方法，其特征在于，還包括對(duì)具有非標(biāo)準(zhǔn)形狀的每個(gè)粗糙網(wǎng)格塊確定孔隙體積的步驟。
23.如權(quán)利要求21所述的計(jì)算機(jī)實(shí)現(xiàn)的方法，其特征在于，還包括計(jì)算每個(gè)粗糙塊的平均深度的步驟。
24.如權(quán)利要求21所述的計(jì)算機(jī)實(shí)現(xiàn)的方法，其特征在于，所得粗化網(wǎng)格是完全非結(jié)構(gòu)化的。
25.如權(quán)利要求21所述的計(jì)算機(jī)實(shí)現(xiàn)的方法，其特征在于，所得粗化網(wǎng)格是基本結(jié)構(gòu)化的。
26.如權(quán)利要求25所述的計(jì)算機(jī)實(shí)現(xiàn)的方法，其特征在于，所述粗化網(wǎng)格在其中至少一個(gè)粗糙單元中存在間斷點(diǎn)的所述粗糙單元周圍是非結(jié)構(gòu)化的。
27.如權(quán)利要求21所述的計(jì)算機(jī)實(shí)現(xiàn)的方法，其特征在于，還包括: 根據(jù)粗化因子來(lái)粗化所述精細(xì)網(wǎng)格來(lái)創(chuàng)建粗糙網(wǎng)格塊；以及 在所述粗化之后，將粗糙網(wǎng)格塊中包括多于一個(gè)間斷點(diǎn)的那些細(xì)分來(lái)創(chuàng)建僅包括一個(gè)間斷點(diǎn)的粗糙網(wǎng)格塊。
28.如權(quán)利要求28所述的計(jì)算機(jī)實(shí)現(xiàn)的方法，其特征在于，所述具有間斷點(diǎn)的網(wǎng)格塊是鋸齒形網(wǎng)格塊。
29.如權(quán)利要求28所述的計(jì)算機(jī)實(shí)現(xiàn)的方法，其特征在于，還包括通過(guò)利用平行管道對(duì)所述鋸齒形網(wǎng)格塊建模來(lái)對(duì)表征所述鋸齒形網(wǎng)格塊的半傳導(dǎo)能力進(jìn)行求和的步驟。
30.如權(quán)利要求1所述的計(jì)算機(jī)實(shí)現(xiàn)的方法，其特征在于，所述非標(biāo)準(zhǔn)形狀是三角形。
【文檔編號(hào)】G06G7/48GK103907118SQ201180073561
【公開(kāi)日】2014年7月2日 申請(qǐng)日期:2011年9月20日 優(yōu)先權(quán)日:2011年9月20日
【發(fā)明者】S·高瑞爾 申請(qǐng)人:界標(biāo)制圖有限公司