【正文】 壓力為常壓，適用于出口和入口界面。 1. 加入的反映裂縫分布的點線面 圖 811 圖 712 ，可以將加載的數(shù)據(jù)應用到裂縫建模中去。 Moving average 滑動平均模擬 PETREL包含幾種確定性建模方法：克里格方法、函數(shù)法、滑動平均法、近點距離等方法， 確定性建模會產(chǎn)生光滑的效果 . 克里格方法包含變差的信息，因此產(chǎn)生的非均質(zhì)模型包含了 3D空間各點的相互地質(zhì)統(tǒng)計的相關(guān)性 . 確定性方法不會產(chǎn)生局部變量，如果選擇 100次實現(xiàn)，各實現(xiàn)都相同，下面是滑動平均法的實現(xiàn)練習。 序貫高斯模擬 Petrel 可以根據(jù)序貫高斯模擬建立隨機巖石物理模型 . 這是用于產(chǎn)生多變量高斯域?qū)崿F(xiàn)的直接算法。 圖 86 46 9） 按 OK 。 圖 811用于約束的趨勢面 圖 812約束后的孔隙度模型 高斯隨機模擬 ，這種算法與 GSLIB的序貫隨機模擬方式不同。 圖 814 49 第 九 章 模型粗化與輸出輸入 模型粗化 模型粗化是使細網(wǎng)格的精細地質(zhì)模型 “ 轉(zhuǎn)化 ” 為粗網(wǎng)格模型的過程。 2. Scale Up Properties: 應用適當?shù)乃惴ò丫毦W(wǎng)格中的屬性逐個重新采樣到粗化網(wǎng)格。可能會報錯，需要重新定義趨勢線和斷層方向使 Pillar Gridding完成， 4) Ok。 3) 確定勾選了 Keep locked nodes unchanged 。 圖 92 Scale Up Properties 1) 雙擊 Process diagram/ Upscaling / Scale Up Properties。 新的局部網(wǎng)格加密 ( LGR )方法： ? 基于層面生成 ； ? 在井射孔段周圍生成 ； ? 漸變的局部網(wǎng)格加密 (井 , 層面 )。如圖（ 107） 圖 97 模型輸出 Petrel 可以輸出多種類型的網(wǎng)格數(shù)據(jù)，涵蓋了絕大多數(shù)主流數(shù)值模擬軟件的輸入格輸入上面步驟建立的Polygong 設置 x， y， z 要細化的網(wǎng)格數(shù) 選擇要加密的 zone 選擇加密區(qū)域所在的 Segment 53 式。 2) 文件名給 ，文件類型選 ECLIPSE Grid keywords (grdecl)(ASCII)(*.*)。 該文件的例子如下： 5) 輸出 Well Completion Data，在 3D Grid (SIMU)上按右鍵選擇 Export? 。MainFaul39。MainFaul39。MainFaul39。INCOMP 39。 該文件的例子如下： COMPLETION SPECIFICATION DATA WELL LOCATION NAME I J K1 K2 COMPDAT 39。 30 46 1 1/ 39。 。 31 46 1 1/ ... / Well connection data for Schedule TRAJECTORY_COLUMN_ORDER MD_ENTRY GRID_I GRID_J GRID_K WELL_ENTRY ENTRY_FACE MD_EXIT WELL_EXIT EXIT_FACE PROP1 PROP2 WELLNAME 34/10A21 WELLHEAD_I 30 WELLHEAD_J 45 TRAJECTORY 30 45 1 Z X 1222 30 46 1 X+ Y 1223 ... END_TRAJECTORY 55 模型輸入 操作步驟 1) 激活 MODEL標簽，單擊 INSERT菜單中 NEW MODEL選項，加入一個新的模型，在 MODEL標簽下會產(chǎn)生一個名為 NEW MODEL 的模型。 30 45 1 1/ 39。彈出文件選擇對話框。X+39。Y+39。Y+39。 6) 文件名 ，文件類型選 ECLIPSE Well Completion Data (ASCII)(*.*)。彈出文件選擇對話框。 1) 首先輸出角點網(wǎng)格，在 3D Grid (SIMU)上按右鍵選擇 Export? 。如圖（ 105） 圖 95 3） 雙擊 ProcessDiagram/structuralmodeling/make local grids 4） 彈出如圖（ 106）的參數(shù)設置對話框， 進行參數(shù)設置。 圖 9－ 3 51 局部網(wǎng)格加密 Petrel 中油藏地質(zhì)模型粗化模塊增加模擬網(wǎng)格設置技術(shù)，允許用戶對完成 的油藏模型進行局部網(wǎng)格加密。 5) OK。 圖 91 50 Scale Up Structure 1) 激活練習 61創(chuàng)建的 3D Grid (SIMU)。 2) 定義粗化網(wǎng)格的網(wǎng)格大小，勾選 Make zigzag type faults。 Petrel2020 引入了新的 IJK fault 粗化算法 ，使得 逆斷層附近不再有洞。 新的協(xié)同高斯隨機函數(shù)模擬通常比 類似的序貫高斯模擬優(yōu)越，主要原因如下： 1） 通常該算法比序貫模擬運算快 2） 快速改變主次變量之間的相關(guān)性參數(shù)設定 3） 在模擬的變量或者主次變量之間的相關(guān)性之間沒有系統(tǒng)誤差。 圖 87沉積相模型 圖 88孔隙度模型 10） 在 Properties 文件夾下 , 右鍵單擊，打開 Calculator 標簽 . 選擇所有的孔隙度實現(xiàn)，并計算他們的平均值 . Mean_por=(por01+por02+por03+por04+por05+por06+por07+por08+por09)/9 11） 顯示平均 Porosity 模型 ， 和 Fluvial Facies 模型進行對比 . 使用 Zone filter 只顯示 Zone1 層 . 有趨勢約束的模擬 1） 運行 Petrophysical Modeling （ 雙擊 Petrophysical Modeling） 2） 選擇 Porosity 屬性進行建模 . 3） 打開 Common settings 設置實現(xiàn)的數(shù)目 10. 4） 打開 Zone Settings ， 選擇 Zone1 進行建模 ， 點擊當前層的 Reset settings ， 按 設置為缺省值 . 5） 選擇 Sequential Gaussian Simulation 方法 . 6） 打開 Petrophysical Modeling 處理流程的設置標簽 ， 定義變差的變程 ， 方向 ， 在Output data range 下點擊 Estimate 定義 minimum 和 maximum 孔隙度值 .對其他相作相同的操作 . 圖 89 7） 建立孔滲飽模型時，我們除了可以用相約束，同時也可以利用其他的趨勢面進行約束，具體的操作步驟為： a). 打開 Trend標簽 (如圖 86) 47 b). 把描述該層的空隙度的平面趨勢圖輸進來 c). 按 APPLY 鍵觀察模型的變化 圖 810 8） 在 Properties 文件夾下 , 右鍵單擊，打開 Calculator 標簽 . 選擇所有的孔隙度實現(xiàn)，并計算他們的平均 值 . Mean_por=(por01+por02+por03+por04+por05+por06+por07+por08+por09)/9 9） 觀察利用趨勢圖約束的效果和相約束效果的區(qū)別。 操作步驟 1． 產(chǎn)生一個模型 : 1) 激活 Geomodel項目下一個模型 ) 2) 打開 Petrophysical Modeling 處理流程 , 選擇 Use Existing Property 和 Porosity 模型 . 3) 激活 Zone層的 Lock icon 按鈕取消 Zone1層的鎖定狀態(tài) . 4) 選擇 Sequential Gaussian Simulation方法 . 5) 在 Variogram選項卡 , 選擇 Exponential Variogram類型 , 2020 =Major Range, 1000 =Minor Range, 3 = Vertical Range 和 － 44=》 OK(如圖93) 圖 83 2． 使用其他方法，觀察其變化： 1) 從新打開 Petrophysical Modeling 處理流程的 Distribution選項卡 2) 在 Output data range 右邊點擊 Estimate 按鈕 (將估計孔隙度的范圍 ). 設定輸出范圍 (例如 . ). 圖 84 45 3) 按 OK，觀察模型的變化 相控模擬 一般情況下巖石物理模型要受沉積相的約束， 這將約束巖石物理數(shù)值， 很明顯在泥巖背景下的地質(zhì)沉積環(huán)境，如河流相沉積，主要是代表河道的網(wǎng)格擁有高的孔隙度值，并不是在整個層中平均化孔隙度值。 8) 在 Properties 文件夾下雙擊 Porosity 屬性模型，在 Statistics 標簽下檢查統(tǒng)計結(jié)果（如圖 92）。隨機建模可以采用巖相模型、地震屬性模型等作為屬性模擬的約束條件。 周期：如果滲透率沿 I方向計算， J+和 J方向?qū)⒕哂型瑯拥闹?，所有的流體從 J+方向流出（或流進）將出 J流進（或流出）， W+/W相同。 邊界條件應用于每個網(wǎng)格單元。 Oda方法是數(shù)據(jù)統(tǒng)計的計算方法，它以單個網(wǎng)格內(nèi)裂縫的總面積及裂縫的不同參數(shù)為基準，估算每個單元滲透率和其它參數(shù)。 39 圖 710 ? 粗化裂縫網(wǎng)格屬性 計算裂縫網(wǎng)格屬性：將裂縫模型轉(zhuǎn)換為雙孔、雙滲模型中的裂縫屬性。 ? 右擊 fracture，選擇 View Calculator 。 37 圖 76 裂縫密度及分布范圍定義 圖 77 裂縫的幾何形態(tài)定義 圖 78 裂縫的方位定義 38 裂縫幾何模型建立 在以上工作基礎上，可以直接生成裂縫幾何模型。 d. 裂縫發(fā)育長度 ： 裂縫發(fā)育長度 可以 通過露頭觀察得到， 可以根據(jù)經(jīng)驗賦值或 由螞蟻追蹤得到。 c. 裂縫方向：由網(wǎng) 格或者趨勢面定義裂縫方向，裂縫傾角和方位角的平均數(shù)和集中度被輸入到 Fisher模型中。切斷需要常態(tài)和冪數(shù)法分配（最大和最小的復選框變灰），并且定義在裂縫網(wǎng)絡明確模擬的并在第一標簽中定義的裂縫強度。 裂縫面積 / 體積 =裂縫表面積 /模型體積，通常應用密度測井曲線 裂縫長度 / 體積 =裂縫長度 /模型體積 裂縫 /體積 = 裂縫數(shù)量 /裂縫體積 使用密度曲線模擬時，允許用戶通過定義裂縫長度，選擇是否模擬最小的裂縫，因為對滲透率貢獻小，因此，模擬結(jié)果可能與 最初的密度屬性值略有不同。 ：定義模型哪個部分參與計算， 垂向上是 某個 zone 或者 某個區(qū)域， 還是整個網(wǎng)格 。 圖 73 裂縫密度曲線生成 35 圖 74 裂縫 累 計 密度 曲線 密度曲線 裂縫參數(shù)確定及預測 Petrel 2020 裂縫建模提供了確定性及隨機模擬兩種預測方法。當產(chǎn)生裂縫強度時，用戶需要提供窗口長度，這個窗口長度是用來平均平滑數(shù)值。 密度測井：對于模型密度是重要的元素，是裂縫網(wǎng)格的基本數(shù)據(jù)。 Petrel2020 裂縫模塊 具有立體圖分析功能，即將所有統(tǒng)計結(jié)果在同一方位坐標系統(tǒng)內(nèi)顯示，根據(jù)其本身性質(zhì)加以濾波和定量歸類。azimuth’ 作為算法 ? 選擇要計算的層面 ? 檢查 ‘Relative’ 檢驗欄 ? 點擊 Apply，新的地層傾角和地層傾角方位被存儲為裂縫點數(shù)據(jù)屬性 裂縫性質(zhì)分析 對于不同地區(qū)，不同的構(gòu)造運動也會產(chǎn)生不同性質(zhì)的裂縫類型（一些可能與構(gòu)造變動有關(guān) ，一些與沉積因素有關(guān)），裂縫網(wǎng)絡是代表裂縫的平面組。數(shù)據(jù)有效格式為 ‘Point well data (ASCII)’ ，即（ x， y， z，裂縫參數(shù)），數(shù)據(jù)以點數(shù)據(jù)類型存儲于 Global Well Logs 文件夾下。 32 第 七 章 裂縫建模 裂縫建模需要多步實現(xiàn)，涉及到油藏描述和模擬的方方