Petrel建模流程

1、contents一、 數(shù)據(jù)準(zhǔn)備二、 數(shù)據(jù)輸入三、 pillar gridding四、 make horizon五、 laying六、 建立幾何模型七、 離散化測(cè)井曲線八、 對(duì)vsh數(shù)據(jù)進(jìn)行分析九、 相建模十、 對(duì)連續(xù)數(shù)據(jù)進(jìn)行分析十一、 屬性建模十二、 網(wǎng)格粗化及屬性粗化的操作十三、 儲(chǔ)量計(jì)算十四、 產(chǎn)生stoiip （烴體積密度分布圖）十五、 輸出數(shù)模所需要的文件主要模塊介紹一、 數(shù)據(jù)準(zhǔn)備本實(shí)例中的數(shù)據(jù)整理如下：wellhead井位坐標(biāo)文件jinghaoxykbtopdepthbottomdepthx21-233973816364714261433.0821502195x21-24397407

2、0364716291433.082156.12193.1x21-253974257364718491433.082154.42190.4x21-263974480364720961436.52154.82189.8x22-193972535364705161407.562120.32152.3x22-203972803364707951417.462139.12165.1x22-213973010364710401379.72102.62135.6welltop分層文件xyhbwellpointsurfacejinghao397381636471426-716.92horizon c811x2

3、1-23397381636471426-724.92horizonc8121x21-23397381636471426-735.92horizon c8122x21-23397381636471426-755.92horizon c813x21-23397381636471426-761.92horizon c821x21-23397407036471629-723.02horizon c811x21-24397407036471629-731.02horizonc8121x21-24397407036471629-742.02horizon c8122x21-2439740703647162

4、9-754.02horizon c813x21-24397407036471629-760.02horizon c821x21-24測(cè)井文件準(zhǔn)備depthperm_kpor_ksw_kvsh_kntg2140.1250.005901002140.250.005901012140.3750.005901002140.50.00590010二、 數(shù)據(jù)輸入1 輸入well header（井位坐標(biāo)文件）右鍵點(diǎn)擊輸入well header:文件類型里選：well heads (*.*)2 輸入well tops（分層文件）: 右鍵點(diǎn)擊well tops文件夾并選擇import (on selection

5、)；文件類型里選：petrel well tops (ascii)3 輸入well logs右鍵點(diǎn)擊wells文件夾，選擇import (on selection)；文件類型：well logs (ascii)input data logs specify logs to be load 加載per, perm, sw, vash, ntg 等數(shù)據(jù)。設(shè)置template settings -9999全選左擊ok然后下圖設(shè)置左擊ok for all 。 在流程窗口左鍵雙擊 出現(xiàn)下圖，在窗口中輸入建模名稱點(diǎn)apply,再點(diǎn)cancel三、pillar gridding建一個(gè)網(wǎng)格邊界（工區(qū)范圍）邊界

6、標(biāo)定了3d網(wǎng)格的側(cè)向延伸。僅在邊界內(nèi)形成3d網(wǎng)格，因此在邊界外不會(huì)進(jìn)行儲(chǔ)量計(jì)算，也不存在構(gòu)造層面和屬性單元?？捎脛?chuàng)建邊界工具，在2d窗口中創(chuàng)建一個(gè)邊界。同時(shí)在這里設(shè)置最終模型的橫向網(wǎng)格大小。步驟：左擊下拉菜單中的并把資源管理器中的wells打勾出現(xiàn)建模的2d圖形：把pillar gridding 單擊點(diǎn)亮：左擊在2d里瀏覽全圖，出現(xiàn)下圖：左擊，然后用鼠標(biāo)左鍵畫邊界：邊界將要封閉前雙擊左鍵，就封閉了.然后左鍵雙擊左擊apply，出現(xiàn)下圖左擊是（y）。pillar網(wǎng)格化的過(guò)程就是一個(gè)空間結(jié)構(gòu)生成的過(guò)程。在i、j方向上定義網(wǎng)格單元的大小。生成的骨架網(wǎng)格（也叫做pillar網(wǎng)格）定義出了空間結(jié)構(gòu)。創(chuàng)建

7、出的骨架網(wǎng)格不代表任何表面，而是代表了pillar頂部、中部和底部的位置。在下一個(gè)進(jìn)程中（創(chuàng)建地層層面）地層層面會(huì)被插入（make horizon），并連接到pillar上，z方向上的網(wǎng)格單元也將被定義。pillar網(wǎng)格化進(jìn)程完成后，首先會(huì)生成一個(gè)3d grid網(wǎng)格。網(wǎng)格化的目的就是要?jiǎng)?chuàng)建均勻分布的矩形網(wǎng)格單元。對(duì)生成的網(wǎng)格結(jié)果感覺(jué)滿意后，點(diǎn)擊ok以開(kāi)始構(gòu)建頂部和底部的骨架網(wǎng)格。在彈出窗口（詢問(wèn)是否將開(kāi)始構(gòu)建頂部和底部骨架網(wǎng)格）中點(diǎn)擊"yes"。四、make horizon在3d骨架網(wǎng)格中加入層面左鍵雙擊：窗口中出現(xiàn)下圖：左擊5次，建立5個(gè)層面。出現(xiàn)下圖：左擊1,然后依次左

8、擊2 、3 ，4把層面加上去。1432出現(xiàn)下圖。并把光滑窗口數(shù)字0改為2。雙擊窗口中出現(xiàn)下圖：直接左擊。再cancel.五、laying細(xì)分層僅僅是網(wǎng)格精細(xì)化的過(guò)程，不是所有輸入資料都用在了這個(gè)進(jìn)程中。用戶可通過(guò)設(shè)置單元的厚度、單元的個(gè)數(shù)或用比例數(shù)，來(lái)定義網(wǎng)格垂向的分辨率。給定單元厚度時(shí)，zone的劃分既可以跟隨頂部也可以按底部。小層本應(yīng)該根據(jù)將要建立的屬性模型來(lái)定義。通常，小層的厚度應(yīng)該是模擬的最薄相的厚度（最薄層處一米一個(gè)layer）。但是，有很重要的一點(diǎn)應(yīng)當(dāng)記住，小層厚度減少時(shí)，單元數(shù)目會(huì)增加，所以不應(yīng)該插入太多的細(xì)節(jié)。步驟：下面創(chuàng)建細(xì)分層：左鍵雙擊出現(xiàn)下圖:把分層數(shù)改成如下圖，兩個(gè)主力

9、層各分成10個(gè)小層。左擊apply,再cancel.六、建立幾何模型（可以不做）幾何屬性是通過(guò)先前定義好的方法，如網(wǎng)格高度（cell height）、總體積（bulk volume）、創(chuàng)建出的3d屬性。每個(gè)網(wǎng)格都將賦予一個(gè)與所選方法相對(duì)應(yīng)的數(shù)值。在進(jìn)行儲(chǔ)量計(jì)算和巖石物理屬性間的數(shù)學(xué)運(yùn)算（如生成含水飽和度屬性）時(shí)可能會(huì)用到。幾何屬性建模進(jìn)程允許用戶建立幾何屬性模型，另外還可進(jìn)行簡(jiǎn)單的建模操作，如接觸面之上的計(jì)算，它是計(jì)算用戶定義的接觸面之上的網(wǎng)格單元的高度。選擇“單元體積”方法；用總體積作為屬性模板，然后點(diǎn)擊應(yīng)用即可生成；左鍵雙擊出現(xiàn)下圖，左擊apply,再左擊cancel.卡泥質(zhì)含量： 打開(kāi)i

10、nput資源管理器（出現(xiàn)下圖）對(duì)右擊出現(xiàn)下圖：左擊設(shè)置如下圖：把vsh測(cè)井曲線處理成離散的值0（代表砂巖段）和1（代表泥質(zhì)段）七、離散化測(cè)井曲線離散化進(jìn)程就是給井曲線穿過(guò)的網(wǎng)格單元賦值。因?yàn)槊總€(gè)網(wǎng)格單元僅能得到一個(gè)值，那就要求測(cè)井曲線要均勻分布，即離散化。其目的就是要在屬性建模時(shí)能把井的信息作為輸入，即控制井間的屬性分布。有一點(diǎn)要明確，離散化之后得到的網(wǎng)格單元將作為屬性的一部分，而不是獨(dú)立出的一項(xiàng)。沿井軌跡的網(wǎng)格單元內(nèi)分布的值與整個(gè)3d離散化之后得到的屬性分布是一致的。左鍵雙擊定義離散化設(shè)置。算法選平均法，以線數(shù)據(jù)處理測(cè)井曲線，使用neighbor cell方法。注意：*對(duì)vshale設(shè)置成如

11、下方式，然后左擊*對(duì)于連續(xù)數(shù)據(jù)perm離散化時(shí)應(yīng)該選用harmonic方法，其它設(shè)置完后左擊*每次設(shè)置完一道測(cè)線參數(shù)并apply后，都要把點(diǎn)上，以免被沖掉（那就白做了）。 *對(duì)于por, sw, ntg選用的方法一樣，都用arithmetic 八、對(duì)vsh數(shù)據(jù)分析數(shù)據(jù)分析的結(jié)果可以直接被相建模和屬性建模的模塊調(diào)用，數(shù)據(jù)分析分為兩類：對(duì)離散數(shù)據(jù)的分析和連續(xù)數(shù)據(jù)的分析。vsh離散數(shù)據(jù)的分析:1.打開(kāi)分析窗口,界面如下: 2.選擇分析的對(duì)象,是經(jīng)過(guò)離散化的井點(diǎn),還是未離散化的測(cè)井曲線,或整個(gè)模型（默認(rèn)）.3.是否使用濾波功能，很少使用（默認(rèn)不用）。4.左擊標(biāo)簽,進(jìn)行相應(yīng)的分析：選中vshu的zone

12、 1依次左擊1 、2 、3 、4 24135、vsh變差函數(shù)的分析對(duì)vshale的4個(gè)層（zone1,zone2,zone3,zone4）分別進(jìn)行變差函數(shù)的分析：以zone1層的泥shale為例,其它2、3、4層的操作相同左擊標(biāo)簽,進(jìn)行相應(yīng)的變差分析：出現(xiàn)下圖：首先設(shè)置主方向的分析參數(shù),包括帶寬,搜索半徑,步長(zhǎng)等,然后再設(shè)置次方向和垂向上的參數(shù),這些參數(shù)的設(shè)置需要用戶對(duì)本地區(qū)數(shù)據(jù)的大概了解的基礎(chǔ)上,否則分析的結(jié)果的可信度大大降低.在該例中的分析參數(shù)和結(jié)果如下圖.在分析變差之前,首先大概了解數(shù)據(jù)的分布情況,然后再調(diào)整這些分析參數(shù),這樣才能達(dá)到比較好的分析效果。1） 先把major directi

13、on （主方向值）左擊，（出現(xiàn)下圖）然后設(shè)置其它參數(shù)，手動(dòng)調(diào)參完畢后左擊apply, 手動(dòng)調(diào)整下圖的1和2的幅度，使得藍(lán)線盡量與小黑點(diǎn)重合，調(diào)好后塊金值nugget設(shè)為0。 12注意：major range 值比minor range 值一般要大些。type值均設(shè)置成spherical(球形的)，no lags 一般根據(jù)經(jīng)驗(yàn)設(shè)置，,本例設(shè)為30 ,這樣條柱較多，容易調(diào)參。2） major direction 手動(dòng)調(diào)參完畢后，再左擊次方向（出現(xiàn)下圖）,開(kāi)始手動(dòng)調(diào)參，調(diào)完后，左擊apply 3） major direction調(diào)參完畢后，左擊vertical direction （出現(xiàn)下圖）垂向厚

14、度一般小于10米，但有時(shí)主力層厚要超過(guò)10米。 塊金值nugget設(shè)為0至此，vshale的zone 1層的泥質(zhì)操作完畢。然后是zone 1層的砂層是同樣操作。本例中調(diào)好的圖形如下：主方向調(diào)試如下： 次方向調(diào)試如下：垂向調(diào)試如下：然后調(diào)試vshu的zone 2 ,zone 3,zone 4 ，重復(fù)前面操作過(guò)程，4個(gè)層都依次變差分析完成后，對(duì)vsh的數(shù)據(jù)分析結(jié)束。九、相建模petrel中有幾種方法可以用來(lái)生成相模型：在這里常使用序貫指示模擬法的相模型隨機(jī)計(jì)算（隨機(jī)建模）。用序貫指示模擬法(sis)，建立一個(gè)基本的相模型的進(jìn)程:1、雙擊下的。 2、打開(kāi)相建模流程3、選擇use existing p

15、roperty,屬性選擇4、選擇zone settings，選擇zone 1 (頂層) ，左擊，把打勾，并設(shè)為2，并把leave zone unchanged 按鈕按下。選中zone1打勾以zone1層為例：5 、選擇zone 1。取消選擇leave zone unchanged按鈕以改變?cè)O(shè)置并選擇序貫指示模擬方法。6、左擊1、2、3、4、各一次，并左擊5、兩次、出現(xiàn)下圖：45321左擊apply運(yùn)行，然后依次對(duì)zone 2,zone 3,zone 4 層做和上述相同步驟。（本例做出的2、3、4層圖形如下）zone 2層相建模如下：zone 3層相建模如下：zone 4層相建模如下：這樣相模型

16、就做完了。十、對(duì)連續(xù)數(shù)據(jù)進(jìn)行分析雙擊流程圖中的出現(xiàn)下圖：分別對(duì)perm,por,sw等連續(xù)數(shù)據(jù)進(jìn)行如下操作：1用鼠標(biāo)選中中的perm，對(duì)perm進(jìn)行操作。2操作：分別對(duì)此窗口中的、左擊，然后依次點(diǎn)擊加入到右側(cè)窗口。出現(xiàn)下圖：先用鼠標(biāo)選中 （反白）進(jìn)行如下操作：依次左擊1、2、3、4、5、6、77432561對(duì)不做任何操作（圖形如下）再用鼠標(biāo)選中 （反白）進(jìn)行如下操作：依次左擊1、2、3、4、5、6431562再用鼠標(biāo)左擊variograms,出現(xiàn)下圖：進(jìn)行perm變差函數(shù)的分析：perm的主方向調(diào)變程圖如下（前面已經(jīng)對(duì)vsh進(jìn)行了分析，perm變差分析與vsh變差分析的步驟相似），分別調(diào)好主方

17、向，次方向，和垂向的變差。（下面選主方向變差圖，次方向和垂向省略）再選中perm zone1的砂體如下圖進(jìn)行操作（步驟重復(fù)以上）：解釋1.選擇分析對(duì)象,是經(jīng)過(guò)離散化的井點(diǎn),還是未離散化的測(cè)井曲線,或整個(gè)模型（默認(rèn)即可）2.是否使用濾波功能;是否用相約束（打勾）3.分別按以下2個(gè)標(biāo)簽,進(jìn)行相應(yīng)的分析4.打開(kāi)每個(gè)標(biāo)簽后,按鍵,刷新顯示5.在進(jìn)行transformation分析時(shí)，可以能夠多種數(shù)據(jù)的處理，包括輸入截?cái)?，?duì)數(shù)變換，正態(tài)分布變換等。6.例如要進(jìn)行輸入截?cái)嗪驼龖B(tài)分布轉(zhuǎn)換,處理的對(duì)話框如下:十一、屬性建模在建立油藏模型時(shí)，petrel提供了確定性模擬和隨機(jī)模擬兩種算法：輸入數(shù)據(jù)為粗化的井模型

18、和趨勢(shì)數(shù)據(jù)以及各種設(shè)置參數(shù)。當(dāng)建立屬性模型時(shí)，所有的網(wǎng)格都賦于數(shù)值。井?dāng)?shù)據(jù)和趨勢(shì)數(shù)據(jù)分布在3d網(wǎng)格中。在建立模型之前，用戶必須進(jìn)行詳細(xì)的數(shù)據(jù)分析，定義變差函數(shù)。隨機(jī)屬性建模petrel 可以根據(jù)序貫高斯模擬建立隨機(jī)巖石物理模型。這是用于產(chǎn)生多變量高斯域?qū)崿F(xiàn)的直接算法。該算法可以產(chǎn)生局部變量。意思是如果基于相同的輸入定義100 個(gè)實(shí)現(xiàn)(使用不同的種子點(diǎn))，可以得到100個(gè)不同的結(jié)果，他們都能和輸入相匹配，但既然通過(guò)分布來(lái)輸入，每個(gè)網(wǎng)格的數(shù)值的值會(huì)根據(jù)這個(gè)分布賦值。如果運(yùn)行50-100個(gè)實(shí)現(xiàn)，各模型之間的差別反映了模型的不確定性。步驟：分別對(duì)perm,por,sw ,ntg等屬性進(jìn)行如下操作。雙擊

19、流程圖里。選中perm，設(shè)置如下：左擊1 、2 、3231然后如下設(shè)置：左擊1、2、3、4、5、6164352perm的zone 1的泥質(zhì)做完，然后做perm的zone1的砂體如下：然后如下設(shè)置：左擊1、2、3、44321依次對(duì)zone 2、 zone 3、zone4操作（ 與zone1層操作一致）。por,sw屬性建模與perm 的操作一致。（選幾幅圖如下） ntg的屬性建模如下：注意：對(duì)于ntg離散數(shù)據(jù)由于沒(méi)有進(jìn)行變程分析，所以對(duì)其操作如下：如下圖所設(shè)置：選中ntg，四個(gè)層同時(shí)做 ，facies先選泥質(zhì)，左擊進(jìn)行變差設(shè)置，塊金值設(shè)為0，變差類型設(shè)為spherical,【spherical

20、aberration（球面像差）】主、次軸、垂向和物源方向設(shè)置如右表所示：設(shè)置完畢后，左擊applyfacies再選砂體操作如下： 至此，屬性建模完畢。十二、網(wǎng)格粗化及屬性粗化的操作為了防止粗化后的網(wǎng)格不符合數(shù)模要求，在做此步驟之前最好copy 3d如下圖：左擊1，再左擊2（復(fù)制）21左擊（粘貼） 出現(xiàn)下圖： 左擊1、2 3，再雙擊3，出現(xiàn)下圖，設(shè)置4、5（網(wǎng)格粗化：本例粗化為80*80網(wǎng)格），左擊6 apply645231 左擊，出現(xiàn)下圖： 左擊，出現(xiàn)下圖： 左擊1、2、3 出現(xiàn)下圖：321 該步驟需要粗化垂向網(wǎng)格數(shù)，將圖中3，10，10，2改為原始分層數(shù)，也可不改（因數(shù)模人員要求而定）本例

21、分層數(shù)據(jù)改成下圖：左擊,再左擊，出現(xiàn)下圖，左擊 在如下的流程圖里，雙擊出現(xiàn)下圖： 左擊，再粗化屬性模型：步驟： 雙擊 出現(xiàn)下圖：這里需要粗化的數(shù)模屬性體：perm、por、sw、ntg，分別選中perm、por、sw、ntg，下面四幅圖如下所示：perm的方法選用harmonic方法，而por、sw、ntg均選用arithmetic【與離散化測(cè)井曲線時(shí)一致！】perm的粗化如下：por的粗化如下：sw的粗化如下： ntg的粗化：十三、儲(chǔ)量計(jì)算儲(chǔ)量計(jì)算是常用模塊之一?？梢跃_計(jì)算每個(gè)層位，每個(gè)斷塊的儲(chǔ)量， 在油水界面的定義(make contact)中形成的油水界面作為儲(chǔ)量計(jì)算的輸入，儲(chǔ)量計(jì)算的

22、設(shè)置直接明了，但還須用戶仔細(xì)檢查所有標(biāo)簽。所有的設(shè)置參數(shù)可以在volume calculation 窗口中進(jìn)行設(shè)置，用戶須定義輸出的格式以產(chǎn)生3d 屬性,報(bào)告,分布函數(shù)) ，以及使用哪些輸入數(shù)據(jù)。用戶定義的報(bào)告將在運(yùn)行后產(chǎn)生，報(bào)告中將詳細(xì)列出每個(gè)層位、每個(gè)斷塊或每個(gè)相的儲(chǔ)量，用戶也可以生成體積密度分布圖(典型的是 hcpv 圖或stoiip 圖).（一）、計(jì)算油水界面上的總儲(chǔ)量要計(jì)算模型的儲(chǔ)量，首先激活正確的模型步驟:1、左擊2、雙擊儲(chǔ)量計(jì)算流程打開(kāi)對(duì)話框：3、 定義一般設(shè)置：a 打開(kāi)要產(chǎn)生的屬性以及定義報(bào)告中包括哪些參數(shù)。b 凈毛比和孔隙度標(biāo)簽，如果沒(méi)有n/g 屬性，使用常數(shù)值（經(jīng)驗(yàn)定）.c 體積系數(shù)bo (無(wú)單位) 和氣油比rs （單位m3/m3）bo和rs的設(shè)置要查報(bào)告得出，（本例西峰白馬中區(qū)數(shù)值分別為1.313和78.4）d 油水界面標(biāo)簽 本例由于沒(méi)有提供油水界面，所以如下設(shè)置：如果有油水界面值，比如-745.8，則可進(jìn)行如下操作:加入油水界面：在資源管理器中左擊，出現(xiàn)下圖： 在流程圖中把下拉菜單打開(kāi)，雙擊，出現(xiàn)下圖：選中，左擊移除（因?yàn)楸纠龥](méi)有油氣接觸）， 出現(xiàn)下圖：在空格處輸入油水界面（查看報(bào)告得出，或從地質(zhì)那里搜集），比如輸入-74