6，所述流程被配置成用來(lái)確定或接收用于執(zhí)行模擬的調(diào)節(jié)數(shù)據(jù)。這個(gè)信息可通過(guò)獲取地質(zhì)信息的地震勘查、測(cè)井或其他技術(shù)采集。例如，圖5為圖示根據(jù)本公開(kāi)實(shí)施方案用于執(zhí)行模擬的調(diào)節(jié)數(shù)據(jù)的實(shí)例的示意圖。在描繪的實(shí)施方案中，黑點(diǎn)(例如，點(diǎn)502和504)表不油氣井遇到砂巖的位置，而白點(diǎn)(例如，點(diǎn)506和508)則表不遇到頁(yè)巖的位置。
[0026]—旦確定了調(diào)節(jié)數(shù)據(jù)，流程400就將開(kāi)始在步驟408中為每條儲(chǔ)層河道選擇種子點(diǎn)。例如，圖6為描繪與一組河道對(duì)應(yīng)的選定種子點(diǎn)的示意圖。在某些實(shí)施方案中，所述流程被配置成用來(lái)確定河道種子點(diǎn)的位置，因此當(dāng)向選定種子點(diǎn)添加了模擬寬度之后，河道將覆蓋已知會(huì)出現(xiàn)砂巖的位置，而不覆蓋已知會(huì)出現(xiàn)頁(yè)巖的任何位置。
[0027]在步驟410，所述流程被配置成用來(lái)為每個(gè)種子點(diǎn)分配傳播方向，然后在步驟412迭代生成每條儲(chǔ)層河道。例如，在一個(gè)實(shí)施方案中，所述流程為每個(gè)點(diǎn)分配描述最大連續(xù)性局部方向的向量，所述方向?yàn)槭荜P(guān)注體積中的方向，當(dāng)?shù)刭|(zhì)性質(zhì)從種子點(diǎn)位置移至最遠(yuǎn)距離時(shí)，所述地質(zhì)性質(zhì)沿所述方向的大小保持大致相同。如本文所提及，所謂的“大致相同”在不同的情況下可能會(huì)發(fā)生變化。例如，如果在受關(guān)注體積很大一部分之上地質(zhì)性質(zhì)高度可變，則“大致相同”的程度可能更高(例如，15%)。然而，如果在受關(guān)注體積很大一部分之上地質(zhì)性質(zhì)穩(wěn)定，則“大致相同”的程度可能更低(例如，5%)。在某些實(shí)施方案中，用戶可調(diào)節(jié)確定“大致相同”的涵義的變量。在一個(gè)實(shí)施方案中，所述流程從地震數(shù)據(jù)或地層傾角測(cè)井?dāng)?shù)據(jù)中推導(dǎo)出每個(gè)種子點(diǎn)的傳播方向。
[0028]在某些實(shí)施方案中，種子點(diǎn)可在任一個(gè)方向被傳播。例如，在一個(gè)實(shí)施方案中，所述流程在步驟418使用地質(zhì)統(tǒng)計(jì)方法(如具有趨勢(shì)模型416的克里格方法414)來(lái)為方位角計(jì)算或構(gòu)建概率分布。作為所屬領(lǐng)域公知的方法，克里格方法為一種從一組有限點(diǎn)處的一組函數(shù)估值建立函數(shù)近似的方法。在一個(gè)實(shí)施方案中，在生成河道時(shí)，所述流程在步驟422被配置成用來(lái)對(duì)方位角的概率分布進(jìn)行取樣，以確定傳播的可信方向。在一些實(shí)施方案中，所述流程被配置成使用附近方位角作為調(diào)節(jié)數(shù)據(jù)，用來(lái)估算每條生長(zhǎng)中心線尖端的方位角分布O
[0029]此外，在某些實(shí)施方案中，克里格方法414可利用變異函數(shù)模型420來(lái)控制河道的曲折度。例如，在一個(gè)實(shí)施方案中，變異函數(shù)模型的范圍控制河道彎曲的距離。舉個(gè)實(shí)例，在一個(gè)實(shí)施方案中，高斯變異函數(shù)模型可用于創(chuàng)建方位角的高度連續(xù)性。舉另一實(shí)例，在曲流河環(huán)境中，所述流程可被配置成使用孔穴效應(yīng)變異函數(shù)模型，使河道呈周期性的方式來(lái)回?cái)[動(dòng)。
[0030]在一個(gè)實(shí)施方案中，在步驟422作為對(duì)方位角概率分布進(jìn)行取樣以確定可信傳播方向的流程的一部分，所述流程被配置成用來(lái)確保中心線不會(huì)太接近已知出現(xiàn)頁(yè)巖的位置。例如，在一個(gè)實(shí)施方案中，當(dāng)中心線被傳播至新點(diǎn)時(shí)，所述流程模擬所述點(diǎn)的河道寬度(例如，使用具有趨勢(shì)模型的克里格方法來(lái)捕獲河道在遠(yuǎn)端方向變窄的趨勢(shì))，并檢查附近(例如，在等于模擬寬度的半徑范圍內(nèi))的任何井。在一個(gè)實(shí)施方案中，如果任何附近井為頁(yè)巖，模擬的方位角將被拒絕，并且將對(duì)方位角概率分布重新取樣。所述流程重復(fù)上述步驟，直到發(fā)現(xiàn)新點(diǎn)與新寬度的可接受組合不與已發(fā)現(xiàn)頁(yè)巖的任何位置相矛盾。
[0031]舉個(gè)實(shí)例，圖7為圖示根據(jù)本公開(kāi)實(shí)施方案的中心線模擬演進(jìn)的三個(gè)截圖的示意圖。從圖7可以看出，中心線為迭代生成，并可從原始種子點(diǎn)的任何方向傳播。
[0032]圖8為圖示根據(jù)本公開(kāi)實(shí)施方案的可能最終實(shí)現(xiàn)方式的三個(gè)實(shí)例的示意圖。最終實(shí)現(xiàn)方式在不使用基礎(chǔ)網(wǎng)格的情況下捕獲了具有一組中心線與河道寬度的砂巖-頁(yè)巖幾何形狀。如圖所描繪，基于中心線的生成方法及模擬寬度，可能的最終實(shí)現(xiàn)方式可能會(huì)有所不同。
[0033]本公開(kāi)實(shí)施方案的一個(gè)優(yōu)點(diǎn)在于，所述流程被配置成只在必要之處生成中心線位置和河道寬度的模擬值(即，只在必要之處執(zhí)行克里格方法和概率分布取樣)，因此，這種風(fēng)格的模擬方法的計(jì)算速度更快。例如，使用本公開(kāi)實(shí)施方案生成如圖8所示的實(shí)現(xiàn)方式只涉及幾千次克里格/取樣計(jì)算。相比而言，如果要使用精細(xì)得足以分辨最窄河道的網(wǎng)格對(duì)所述區(qū)域進(jìn)行柵格化，將會(huì)有成千上萬(wàn)個(gè)像素，每個(gè)像素均要求克里格計(jì)算和概率分布取樣。類似地，使用多點(diǎn)統(tǒng)計(jì)方法時(shí)，用采集所需統(tǒng)計(jì)數(shù)據(jù)替換克里格法進(jìn)一步降低了計(jì)算速度。
[0034]此外，另一個(gè)優(yōu)點(diǎn)在于，本公開(kāi)實(shí)施方案以便于后處理的方式表示河道結(jié)構(gòu)。例如，通過(guò)將河道結(jié)構(gòu)表現(xiàn)為一組中心線和河道寬度，就很容易對(duì)實(shí)現(xiàn)過(guò)程進(jìn)行后處理，以提取河道形態(tài)信息。例如，使用本公開(kāi)實(shí)施方案，就很容易識(shí)別彎曲河道的內(nèi)緣與外緣。然而，要從像素化圖像中提取同樣的信息卻要困難得多。在將模擬孔隙度-滲透率值疊加到儲(chǔ)層巖相模擬上時(shí)，或者在完善巖相模擬以便包含如越岸沉積等次要性質(zhì)時(shí)，此類信息通常都非常有用。
[0035]另外，本公開(kāi)實(shí)施方案的另一個(gè)優(yōu)點(diǎn)在于，如圖8所示的所述最終實(shí)現(xiàn)方式比如圖3所示的、采用SNES頂算法的最終實(shí)現(xiàn)方式“更清楚”(例如，河道的邊緣更清楚、更清晰)。雖然很容易把視覺(jué)清晰度只看作為美學(xué)優(yōu)點(diǎn)，但是有生產(chǎn)規(guī)劃研究表明，砂巖中連續(xù)通路的存在與否非常重要。例如，某些研究發(fā)現(xiàn)短尺度變異過(guò)多的像素化圖像可能會(huì)低估儲(chǔ)層流體與生產(chǎn)模擬流體長(zhǎng)距離有效移動(dòng)的能力。
[0036]此外，與其他地質(zhì)統(tǒng)計(jì)模擬程序要求的主要參數(shù)相比，上述模擬方法要求的主要參數(shù)更便于大多數(shù)地質(zhì)學(xué)家使用。例如，其中兩個(gè)關(guān)鍵參數(shù)為河道走向與寬度，這兩個(gè)參數(shù)可從訓(xùn)練圖像(如果有的話)中提取，或者可由儲(chǔ)層地質(zhì)學(xué)家直接提供，因?yàn)榕c兩點(diǎn)統(tǒng)計(jì)法(變異函數(shù))和多點(diǎn)統(tǒng)計(jì)法相比，這些幾何性質(zhì)更加易于理解。唯一其他參數(shù)，即控制河道曲折度的變異函數(shù)模型，也可從訓(xùn)練圖像中提取，或者如果沒(méi)有訓(xùn)練圖像，也可通過(guò)試錯(cuò)法進(jìn)行校準(zhǔn)，因?yàn)樗龀绦蚩煽焖僭龠\(yùn)行，以便目視檢查選擇不同變異函數(shù)模型所產(chǎn)生的影響。
[0037]圖9為圖示用于實(shí)現(xiàn)本公開(kāi)實(shí)施方案的特性與功能的系統(tǒng)900的一個(gè)實(shí)施方案的方框圖。通常，在一個(gè)實(shí)施方案中，除了其他組件以外，系統(tǒng)900包括處理器900、主存儲(chǔ)器902、二級(jí)存儲(chǔ)單元904、輸入/輸出接口模塊906和通訊接口模塊908。處理器900可為能夠執(zhí)行用于實(shí)現(xiàn)本公開(kāi)實(shí)施方案的特征和功能之指令的任何類型或任何數(shù)量的單核或多核處理器。
[0038]輸入/輸出接口模塊906使系統(tǒng)900能夠接收用戶輸入(例如，從鍵盤和鼠標(biāo))并將信息輸出至一個(gè)或多個(gè)設(shè)備，如(但不限于)打印機(jī)、外部數(shù)據(jù)存儲(chǔ)設(shè)備和音箱。系統(tǒng)900可視情況包括獨(dú)立顯示模塊910，從而使信息可被顯示在集成或外部顯示設(shè)備上。例如，顯示模塊910可包括用于提供與一個(gè)或多個(gè)顯示設(shè)備相關(guān)的增強(qiáng)型圖形、觸屏和/或多點(diǎn)觸控功能的指令或硬件(例如，顯卡或芯片)。
[0039]主存儲(chǔ)器902是存儲(chǔ)正在執(zhí)行的指令/數(shù)據(jù)或用于執(zhí)行的預(yù)取指令/數(shù)據(jù)的易失性存儲(chǔ)器。二級(jí)存儲(chǔ)單元904是用于存儲(chǔ)持久性數(shù)據(jù)的非易失性存儲(chǔ)器。二級(jí)存儲(chǔ)單元904可為或可包括任何類型的數(shù)據(jù)存儲(chǔ)組件，如硬盤驅(qū)動(dòng)器、閃盤驅(qū)動(dòng)器、或存儲(chǔ)卡。在一個(gè)實(shí)施方案中，二級(jí)存儲(chǔ)單元904存儲(chǔ)計(jì)算機(jī)可執(zhí)行代碼/指令以及使用戶能夠?qū)崿F(xiàn)本公開(kāi)實(shí)施方案的特性與功能的其他相關(guān)數(shù)據(jù)。
[0040]例如，根據(jù)本公開(kāi)實(shí)施方案，二級(jí)存儲(chǔ)單元904可永久性地存儲(chǔ)用于實(shí)現(xiàn)上述河流三角洲環(huán)境流程的無(wú)網(wǎng)格模擬的可執(zhí)行代碼/指令920。然后可執(zhí)行代碼/指令920可由處理器900從二級(jí)存儲(chǔ)單元904載入主存儲(chǔ)器902中用于實(shí)現(xiàn)本公開(kāi)實(shí)施方案。
[0041]此外，在一些實(shí)施方案中，系統(tǒng)900使用通訊接口模塊908與通訊網(wǎng)絡(luò)930進(jìn)行通訊。例如，網(wǎng)絡(luò)接口模塊908可包括網(wǎng)絡(luò)接口卡和/或無(wú)線收發(fā)機(jī)，使系統(tǒng)900能夠通過(guò)通訊網(wǎng)絡(luò)930和/或直接利用其他設(shè)備發(fā)送和接收數(shù)據(jù)。通訊網(wǎng)絡(luò)930可為任何類型的網(wǎng)絡(luò)，包括下列一個(gè)或多個(gè)網(wǎng)絡(luò)的組合:廣域網(wǎng)、局域網(wǎng)、一個(gè)或多個(gè)專用網(wǎng)、因特網(wǎng)、電話網(wǎng)(如公用交換電話網(wǎng)(PSTN))、一個(gè)或多個(gè)蜂窩網(wǎng)