【正文】 點(diǎn)距離通過求解克里格方程組來確定 u 點(diǎn)的概率，并沒有考慮 dn的幾何形狀和各條件數(shù)據(jù)的配位關(guān)系。 多點(diǎn)地質(zhì)統(tǒng)計(jì)學(xué)與兩點(diǎn)地質(zhì)統(tǒng)計(jì)學(xué)的主要區(qū)別在于上面的概率的確定方法不同，它首先引入一訓(xùn)練圖像，通過在訓(xùn)練圖像中尋找與待估點(diǎn)內(nèi)條件數(shù)據(jù)分布完全相同的事件的個(gè)數(shù)來確定概率分布，因此它可以反映出多個(gè)位置的聯(lián)合變異性。 圖 62 變差函數(shù)與多點(diǎn)統(tǒng)計(jì)地質(zhì)學(xué)對比 如上圖，彎曲河道的 3種不同的空間結(jié)構(gòu) (圖 a， b， c)在橫向上 (東西方向，圖 d)和縱向上 (南北方向，圖 e)的變差函數(shù)十分相似，這說明應(yīng)用變差函數(shù)不能區(qū)分這 3種不同的空間結(jié)構(gòu)及幾何形態(tài)，因此，基于變差函數(shù)的傳統(tǒng)地質(zhì)統(tǒng)計(jì)學(xué)插值和模擬方法難于精確表征具有復(fù)雜空間結(jié)構(gòu)和幾何形態(tài)的地質(zhì)體。 舉例來說，當(dāng) Petrel的 Geobodies技術(shù)從地震數(shù)據(jù)中提取了一些特定的地質(zhì)體，并且想把這個(gè)從地震中局部提取的地質(zhì)體做到模型中去，用戶怎樣在不改變模型其它部分（模擬好的相模型）的情況下而實(shí)現(xiàn)這一目的呢？或者模型新打了一口井，由于儲層分布的變化，需要更新這個(gè)相模型，很顯然，我們建模 人員在具體實(shí)現(xiàn)的操作過程中顯得非常困難。 因此，從目前的發(fā)展來看，建模發(fā)展具有以下兩個(gè)優(yōu)勢： （ 1）多點(diǎn)地質(zhì)統(tǒng)計(jì)學(xué)是今后地質(zhì)統(tǒng)計(jì)學(xué)發(fā)展的主要方向，它可以聯(lián)合反映空間多個(gè)位置點(diǎn)的幾何形狀和相 互配位關(guān)系，在模擬具有復(fù)雜形狀地質(zhì)體分布時(shí)，它比序貫指示模擬算法具有更大的優(yōu)勢。但是，建立在兩點(diǎn)統(tǒng)計(jì)關(guān)系上的變差函數(shù)本身在描述儲層非均質(zhì)性上有很大的不足， 它反映的僅僅是空間兩點(diǎn)之間的相關(guān)性，不能充分描述復(fù)雜幾何形狀砂體的空間連續(xù)性和變異性（如河道砂體和沖積扇砂體）。 總之，通過這兩年的創(chuàng)新及研發(fā)， Petrel軟件的相建模方法不斷得到了改進(jìn)，為用戶提供了更多的思路和好的實(shí)現(xiàn)方法。新的算法運(yùn)用了截?cái)喔咚鼓M，充分考慮了所有井資料從而生成每條河道，這種方法的最大特征是能夠以較快的速度完成模擬井密度較大的區(qū)塊。在這個(gè)版本中，提供了一個(gè)估算以后用到相模擬中的概率趨勢的方法。采用沿海平面設(shè)定的插值算法在某些狀態(tài)下更加合理，例如對飽和度的插值模擬（屬性插值）。 2） 基于快速傅立葉變換的高斯模擬算法 一種新的新的高斯模擬算法在 Petrel ，這種算法與 GSLIB的序貫隨機(jī)模擬方式不同。同時(shí)用戶可以導(dǎo)入自己的算法和人工賦值的方法，建立沉積相模型。 Petrel2020 在原有版本基礎(chǔ)上對相建模方法做了較大的改進(jìn)，主要體現(xiàn) 在以下四個(gè)方面： 1）全新的 MPS多點(diǎn)統(tǒng)計(jì)相模擬算法 在 Petrel2020 版本中，引進(jìn)了多點(diǎn)地質(zhì)統(tǒng)計(jì)學(xué)相模擬方法，該方法的引進(jìn)改變了過去傳統(tǒng)的兩點(diǎn)統(tǒng)計(jì)地質(zhì)學(xué)方法，而發(fā)展為多點(diǎn)地質(zhì)學(xué)，解決了過去 兩點(diǎn)統(tǒng)計(jì)關(guān)系上變差函數(shù)的不足，特別是對儲層非均質(zhì)性描 述上的不足，多點(diǎn)統(tǒng)計(jì)地質(zhì)學(xué)能夠充分描述復(fù)雜幾何形狀砂體的空間連續(xù)性和變異性。 A 它比 SGS運(yùn)算速度提高了很多 B 它不是序貫算法 C 它可以并行運(yùn)算 D 它可以進(jìn)行快速的協(xié)同模擬設(shè)定 如同上面提到的，這種高斯算法不同于序貫?zāi)M的序貫算法，允許并行計(jì)算，采用的算法是傅立葉變換算法，這種算法具有快速、并行、在大的范圍變程內(nèi)優(yōu)選最合理的變程等優(yōu)點(diǎn)，這種算法的界面與序貫高斯模擬算法有些類似。 4） 基于目標(biāo)模型的條件模擬 在河道模擬流程中，不論是河流相模擬還是適應(yīng)性模擬，在 ，都新增加了 3D屬性模擬約束條件。這可以被認(rèn)為是垂向比例曲線的標(biāo)準(zhǔn)化。 另外，可以在多口特定的井中通過指定同一條河道連同體的方式，對河道目標(biāo)體的模擬加以限制，增強(qiáng)模擬的確定性信息。 但 由于相建模涉及因素較多，要想建立較準(zhǔn)確的相模型， 還 需要大量的練習(xí)，和豐富的工作經(jīng)驗(yàn)，本手冊只給出常規(guī)的操作步驟，還需要用戶在以后 的學(xué)習(xí)和工作中不斷實(shí)踐。 當(dāng)井資料較少時(shí)，用于計(jì)算變差函數(shù)的點(diǎn)對很少，它也就不能正確反應(yīng)空 3 間兩點(diǎn)之間的相關(guān)性。 （ 2）利用多點(diǎn)地質(zhì)統(tǒng)計(jì)學(xué)模擬算法具有快速靈活性，模擬的巖相展布圖具有一定的真實(shí)性，它為儲層參數(shù)的模擬奠定了基礎(chǔ)。 多點(diǎn)統(tǒng)計(jì)地質(zhì)學(xué)相模擬方法對這類問題的解決提供了可能，對于該方法的更多說明，可以參閱多點(diǎn)統(tǒng)計(jì)地質(zhì)學(xué)相模擬用戶手冊，在第一次采用這種方法之前建議最好查閱手冊說明。 在 常規(guī)兩點(diǎn)地質(zhì)統(tǒng)計(jì)學(xué)中 ，變程起到了一個(gè)顯著地作用，主要是在模擬時(shí)把變程用來作為控制條件模擬屬性，而對于多點(diǎn)統(tǒng)計(jì)地質(zhì)學(xué)相模擬 ，訓(xùn)練圖像也 起到了相應(yīng)的顯著作用，一個(gè)訓(xùn) 練圖像就是三維模型的基礎(chǔ)，用來表明我們想要模擬模型的典型幾何分布形態(tài)，在 ，訓(xùn)練圖像可以通過以下幾個(gè)方式來得到，一是在用布爾模型在規(guī)則的網(wǎng)格上建立一個(gè)簡單的無條件 3D模型，二是通過手工畫相，三是從 Petrel軟件外部導(dǎo)入模型（如 3D數(shù)字露頭模型）。 例如，計(jì)算圖 1 (a) 中 u 點(diǎn)的概率時(shí)，相應(yīng)的條件數(shù)據(jù)場為 dn = { z ( u1) , z ( u2) , z ( u3) , z ( u4) }，其基本方法是首先要在訓(xùn)練圖像 (b)中尋找與圖 (a)中數(shù)據(jù)分布完全相同的事件的個(gè)數(shù)，即要在訓(xùn)練圖像中找出與圖 (a)幾何完全相同的區(qū)域，同時(shí)在該區(qū)域中相同的位置處 z ( u1) , z ( u2) , z ( u3) ,z ( u4) 的值完全相同。 基于上述原理， Sebastien St rebelle 提出了 snesim 模擬算法，利用該算法可以快速、靈活地模擬巖相分布。 (4) 重復(fù)上一步模擬，直到所有的點(diǎn)全被模擬。這樣的信息內(nèi)容被保存在一種叫做多點(diǎn)統(tǒng)計(jì)地質(zhì)學(xué)模式樹中，盡管這種典型的模式樹比較大，而且用了相對簡單的數(shù)據(jù)格式，并且這種模式是不能看見的，內(nèi)置到軟件中的，但是非常重的一點(diǎn)是這種模式樹能夠非常簡單的被產(chǎn)生和使用（這個(gè)圖僅僅來演示模式，但具體的模式內(nèi)容，用戶是看不到的）（圖 74）。 對不同的模擬區(qū)域單元可以此采用不同的訓(xùn)練圖像模式（圖 78），并且允許用戶改變各相之間的相關(guān)系數(shù)，只有激活的區(qū)域才被模擬，各模擬區(qū)域單元之間可以單獨(dú)或者按照順序來模擬，后來模擬的區(qū)域單元會依賴于早先的 區(qū)域單元模擬結(jié)果，其他系統(tǒng)的參數(shù)控制著模擬時(shí)相的分布比例。 圖 610 多點(diǎn)相模擬時(shí)圖像旋轉(zhuǎn)及調(diào)整 10 5）查看模擬結(jié)果 自定義目標(biāo)體 用戶定義目標(biāo)體的生成 用戶定義目 標(biāo)體的生成可以生成并使用包含一個(gè) 或多個(gè)相代碼的非參數(shù)目標(biāo)體。兩種類型 之間可以進(jìn)行轉(zhuǎn)換，盡管需要根據(jù)目標(biāo)體的復(fù)雜情況而定，結(jié)果的質(zhì)量卻 是千差萬別。得到訓(xùn)練圖像 的方法有很多種，許多公司應(yīng)用露頭 來建立訓(xùn)練圖像（以 3D 屬性 輸入 Petrel）。操作不需要離開這個(gè)界面，訓(xùn)練圖像生成后還可以自動(dòng)進(jìn)行檢測。這種類型網(wǎng)格唯一的用途就是存儲訓(xùn)練圖像。用戶定義的目標(biāo)體中插入了一個(gè)按鈕，如果相建模流程打開訓(xùn)練圖像的網(wǎng)格是激活的話，這個(gè)按鈕只用于目標(biāo)體建模（“ Model”面板下黑體顯示）。 總之，現(xiàn)在通過用訓(xùn)練圖像得到的概念模型來產(chǎn)生一個(gè)地質(zhì)實(shí)現(xiàn)的結(jié)果會更容易。這是一個(gè)非常有用的選項(xiàng)，可以用來再屬性中復(fù)制地質(zhì)體，編輯它們并且用它們來生成訓(xùn)練圖像。在這個(gè)版本中，這個(gè)功能只能用在訓(xùn)練圖像網(wǎng)格中的無條件模擬。作為輸入，算法在井處得到粗化的相數(shù)據(jù) ，然后用塊克里金算法估算每種相的平均概率。 練習(xí)一：如何建立相的概率趨勢模型？ 1） Upscaled well log data，相數(shù)據(jù)離散化； 2）定義變差函數(shù)，也可以使用數(shù)據(jù)分析功能進(jìn)行變差函數(shù)分析； 3）定義克里金均值，設(shè)置 2D 或 3D 約束（如下圖）。 在“ Expert”一欄中，可以選擇簡單或普通克里金算法 。我們可以看到河道和決口壩估算的概率體。 3) 在建模設(shè)置標(biāo)簽下， 在 Existing Property Status下拉菜單中選擇 Fluvial Facies Modpel。點(diǎn)擊蘭色箭頭允許用戶從 條件數(shù)據(jù) 評估相 比例 (圖 711)。 10) 最后 點(diǎn)擊 Apply生成模型 。 圖 612 河道連同體 2） 激活模型 New Model/3D Grid (DC)。 6） 從下拉菜單的方法選擇的相體標(biāo)簽選擇目標(biāo)體建模（隨機(jī)）并點(diǎn)擊 添加一個(gè) adaptive channel， 如圖 613。 10）模擬結(jié)果如圖 714. 圖 614 自適應(yīng)河道模擬結(jié)果 C 新增加功能介紹 在河道模擬流程中，不論是河流相模擬還是適應(yīng)性模擬，在 Petrel2020版本 中，都新增加了 3D屬性模擬約束 條件（圖 615） 。 采用這種河道遷移算法，就會更加具有適應(yīng)性，解決了簡單的嵌入式點(diǎn)算法只能簡單考慮河道初始點(diǎn)分布的問題，而能夠使這種概率約束到河道的生長方向和偏移方向了。 序貫指示模擬 是 一種基于象元的隨機(jī)算法，針對每種相分別設(shè)置各自的變差函數(shù)和所占比例，在沉積概念模型和地震屬性模型約束下建立沉積相模型 A 練習(xí)步驟 1) 在處理進(jìn)程表里雙擊 Facies Modeling打開設(shè)置對話框（如圖 717） 20 圖 617 序貫指示模擬設(shè)置 2) 先 選擇離散化屬性數(shù)據(jù) Facies（測井曲線 scale up或者用戶自己解釋的單井相）。如果數(shù)據(jù)分析都完整作過，而且結(jié)果也已經(jīng)保存，可以使用 這一系列按鈕來保證每個(gè) Zone的分析結(jié)果被調(diào)用。 21 圖 618 序貫指示模擬協(xié)同變量設(shè)置 多級相控建模 前面已經(jīng)提到， Petrel2020的相建模增加了一個(gè)相控的功能（從 2020版本中引進(jìn)），相控建模不再是專指相控的物性建模了，對于離散化數(shù)據(jù)類型同 樣可以進(jìn)行相控，這種模擬方法源于分級建模的思想，對于具有較多的相類型，尤其像大型三角洲等沉積環(huán)境，我們就可以采用從相－亞相－微相這種分布建模的方式，充分刻畫出沉積相模型。 圖 620 未采用分級相控的設(shè)置 23 圖 621 多級相控設(shè)置 截?cái)喔咚冠厔?相建模 和截?cái)喔咚鼓M 對于儲層連續(xù)性性好的沉積以及碳酸鹽巖的儲層沉積模擬， 可以考慮截?cái)喔咚冠厔菹嘟?和截?cái)喔咚鼓M。這個(gè)進(jìn)程是不可返回的，因此在您開始畫相之前，做一個(gè)模型備份。 5) 點(diǎn)擊刷子 按鈕。 8) 改變刷子類型并保持繪制（如圖 625）。當(dāng)新的局部數(shù)據(jù)需要加入一個(gè)大的模型（比如一口新鉆的井），剩下的區(qū)域已經(jīng)進(jìn)行了歷史匹配并且不想對以前的工作進(jìn)行重做時(shí)，進(jìn)行局部更新就是一個(gè)特別好的辦法。多邊形是任意的并且可以把多個(gè)封閉的多邊形組成一個(gè)多邊形體。系統(tǒng)建議的缺省值是 3 個(gè)網(wǎng)格。有一口新鉆的井，需要把它的相數(shù)據(jù)和巖石物理屬性數(shù)據(jù)加入模型。 3. 打開“ Facies modeling”流程更新相屬性。 5. 對所有要更新的層設(shè)置模擬的參數(shù)并運(yùn)行，所有多邊形內(nèi)的網(wǎng)格會得到一個(gè)新的值。 Petrel2020裂縫模塊主要功能 Petrel2020 裂縫模塊主要有兩部分功能，即建立離散型裂縫幾何模型建立及和模型粗化裂縫等效物性模型建 立。 可以通過右擊裂縫圖標(biāo)選擇插入屬 性來添加屬性。因此對不同性質(zhì)裂縫類型歸類非常重要。 累積曲線：用來將油藏細(xì)分為力學(xué)區(qū)域的一種客觀方法。對于密度測井構(gòu)造窗的長度對于曲線的光滑度很重要，例如：排出個(gè)別的裂縫的作用，用平均的值畫點(diǎn)，裂縫密度用累積對數(shù)的導(dǎo)數(shù)計(jì)算，與使用者定義的權(quán)一致或有偏差。 ? 右擊測井曲線文件夾下的裂縫數(shù)據(jù)，在下拉菜單選擇插入全井密度曲線。 圖 75 螞蟻體形成的裂縫片（已經(jīng)刪除斷層及錯(cuò)誤信息）及確定性建立 DFN 2）隨機(jī)建模方法 隨機(jī)性建?？梢杂脭?shù) 學(xué) 或者屬性體計(jì)算 ，同時(shí)可以 地震數(shù)據(jù) 或其它二 維 成果 圖件 加以 約束 模擬 。 密度屬性 體 定義每個(gè)單元體積內(nèi)的裂縫數(shù)。 裂縫幾何形態(tài)：在 Petrel中裂縫被簡化為面。 需要指出的是指數(shù)和冪函數(shù)分布模型要求輸入最大值和最小值，控制裂縫數(shù)字和裂縫尺度。 這些值可以在空間中以實(shí)際方位來衡量，或者在 3D模型或者平面模型中用坐標(biāo)來計(jì)量。 裂縫發(fā)育 密度通常與斷層距離 、 油藏沉積類型和地層曲率有關(guān)，其中一些可以通過建立幾何模型得到。 ? 在裂縫網(wǎng)絡(luò)模型流程下打開 ‘Create discrete fracture work’ ? 選擇創(chuàng)建新的裂縫模型 ‘Create new fracture work’ ? 在輸入數(shù)據(jù)窗口中，選擇自動(dòng)斷裂系統(tǒng)解釋結(jié)果輸入。 aperture / transmissivity 關(guān)系確定主要以來于裂縫滲流公式。 Petrel 2020 中包含有兩種滲透 率離散方法，兩種方法 — Oda方法和基于流體的離散化方法，兩種方法都是全張量模擬方法。 基于流體的離散化是一個(gè)嚴(yán)格的方法，它為每各網(wǎng)格都進(jìn)行特別的限定，并在壓力梯度下進(jìn)行流動(dòng)模擬，以計(jì)算每個(gè)方向的滲透率。 常壓：選擇界于壓力最高值和最低值之間的