1、利用petrel詳細建模操作方法利用petrel詳細建模操作方法 編輯整理：尊敬的讀者朋友們：這里是精品文檔編輯中心，本文檔內容是由我和我的同事精心編輯整理后發布的，發布之前我們對文中內容進行仔細校對，但是難免會有疏漏的地方，但是任然希望（利用petrel詳細建模操作方法）的內容能夠給您的工作和學習帶來便利。同時也真誠的希望收到您的建議和反饋，這將是我們進步的源泉，前進的動力。本文可編輯可修改，如果覺得對您有幫助請收藏以便隨時查閱，最后祝您生活愉快 業績進步，以下為利用petrel詳細建模操作方法的全部內容。petrel操作流程1. 前期數據準備地震數據體，斷層線fault lins or 斷

2、層棍fault sticks，fault polygons，數字化的等值線。工區內各井的坐標，頂深，海拔，底深（完鉆井深)，東西偏移,方位角，傾角，砂巖分層數據,砂層等厚圖，測井曲線（公制單位)，單井相，各層沉積相圖，砂巖頂面構造圖，單井巖性劃分，測井解釋成果表，含油面積圖。（在編輯數據的過程中，命名文件時最好數據文件名都和井名一致)2. 數據加載加載井口數據（well headers）well_name x y kb top bottom symbol井名 x坐標 y坐標 海拔 頂深 底深(完鉆） 井的類型加載井斜數據(well path） 第一種數據格式md tvd dx dy azim

3、incl 斜深 垂深 東西偏移 南北偏移 方位角 傾角第二種數據格式md incl azim第三種數據格式tvd dx dy(單井用well logs，多井加井斜可用 production logs）加載分層數據（well tops）（包括斷點數據)md wellpoint 層名 well name-1500 horizon nm31 np11600 fault nm32 np1 以well tops加載之后刪除系統的缺省項,新建4項,對應輸入數據的列，名稱進行編輯,sub-sea z values must be negative!(低于海平面的z值都為負），該選項在編輯時不要選中加載測井曲

4、線（well logs） las格式文件md resis ac sp gr 曲線采用0.125m的點數據（1m8個點數據），注意有的曲線單位要由英制轉換為公制，如：ac 英制單位 s/in要換成工制單位s/m，再用轉換程序轉換為las格式文件進行輸入，以提高數據的加載速度. 如果有孔滲飽數據,按相同格式依次排列即可。在/input data中設置數據的排列順序,曲線內容較多，系統缺省項只有md，所以要用specify to be loaded 定義新的曲線，對應加載數據的列數,名稱和屬性進行編輯。如:dept 1 measured depth dept sp 4 spontaneous pot

5、en sp在/units中可設置輸入數據的單位在/setting 中設置參考高度:矯正海拔or平均海平面or其他加載單井相（well logs) prn格式文件md(斜深) face巖性代碼在曲線上劃分不同層所屬的相，劃相時讀取每層的頂值，不同的相取不同的代碼，如:河道 1，天然堤 2，決口扇 3 .以well logs加載之后，使用默認的md,使用類似曲線的添加方式，屬性取巖性lithologies。備注 點沙壩 bar point 天然堤 natural levee 河道 channel 決口扇 crevasse splay 泛濫平原 background floodplain 河道沙壩

6、channel bar在加載完數據之后,global well logs 中會出現face項，在/setting中設置每個相代碼所代表的名稱，顏色和圖案。加載測井解釋成果表（well logs） prn格式文件top bottom explain well name頂深 底深 解釋結果 井名加載方式和加相一樣，將油層,差油層,氣層，水層等進行編碼，在加載完數據之后，global well logs 中會出現explain項，在/setting中設置每個代碼所代表的名稱，顏色和圖案.使用甲方給定的最后一次測井解釋成果表加載構造圖x坐標 y坐標 類型 z值加載之后可選擇以點或線輸入，然后設置列數加

7、載斷層polygon（閉合）數據 x y 類型 （z)加載地震數據3. 構造模型 stractural modelingmake/edit surface 成面選擇input 數據中加載進來的每個層的等值線或等厚線或polygons，用戶自己根據對應的層名命名，如果輸入的有邊界，就用邊界，沒輸入邊界的就點擊選擇自動生成邊界。在geometry中，幾何網格模型選擇自動模式，設置grid increment，如果輸入的有邊界，就不選boundary。在algorithm中選擇成面所用的算法，不同的面用不同的方法，常用的有構造面,等厚圖,相圖，在對所研究區域的地質情況有充分了解的時候可以在不同的算法

8、里進行不同的設置。在做相圖的時候可選擇協克里金方法或者模擬算法。define model新建一個模型，用戶自己命名.fault modeling斷層模型在3d窗口顯示斷層polygon，用戶輸入的斷層沒有z值，要給斷層賦z值：雙擊斷層polygon，在/calculation下用公式a=surface，surface選擇該斷層所對應的面.編輯斷層polygon數據,例如，在3d窗口顯示編輯好的 1號斷層polygon，該斷層在每個層位的polygon都顯示出來，然后create new fault from polygons，然后根據斷層的掉向和組合方式對新建的斷層進行編輯。新生成的斷層可根據

9、用戶需要，自己定義編號,以便管理。pillar griding生成三維框架編輯完斷層之后，系統自動彈出2d窗口，顯示出等值線和斷層線，用戶自己繪制一個邊界，(如要設定斷層走向要加趨勢線），點擊create new 創建新3d grid，設定網格的平面步長increment，并且選擇斷層處的網格編輯方式。其他設置按用需要進行設置。appiy,觀察網格分布，不符合要求再次調整，調整好后,按ok完成。make horizons 在horizon中根據用戶分層添加層數，添加對應層位的分層，和surface，在fault中設置斷層的距離,一般情況下，不同的層有不同的斷層，在列表中選擇斷層，某一層中沒有的

10、斷層，就不激活.設置完所有的層位，在 /setting 中完成相關設置，點擊完成.make zones 對每個層進行編輯,在隔層存在的情況下，每層有頂底兩個surface，再加上該層的分層,如有其他需要 ，在 setting和well adjustment 中再做設置,如不整合設置和光滑設置。 layering在layering中，設置每個層在縱向上的分層，根據甲方要求劃分縱向分層，隔層用1層就行了 。 4. 相建模scale up well logs 離散化測井曲線選取工區內有單井相數據的井，對單井相數據進行離散化。data analysis 數據分析相數據分析proportion（比例分析

11、)原井數據會形成直方圖，點擊fit active/all curve（s) to histograme，觀察在比例分布里是否有異常值,如有，則要進行更改thickness（厚度分析）觀察各相的厚度比例。probability（概率分析）variograms（變差函數分析）根據數據分布的相關性，調節主，次和垂向變程,調節物源方向，選擇分布模型.河道中,曲流河的主變程一般范圍在12001500m，辮狀河主變程一般在900-1200m左右。（經驗值）一般情況下，曲流河的變程比辨狀河的稍大，次變程一般去主變程的1/3左右垂向上的變程一般取河道沉積的單砂體厚度？數據分析中的地震屬性體是起約束作用？變差函

12、數是區域化變量空間變異性的一種度量,反映了空間變異程度隨距離而變化的特征。強調三維空間上的數據構形，從而可定量描述區域化變量的空間相關性，使克里金技術以及隨機模擬的一個重要工具。變程（range) ：指區域化變量在空間上具有相關性的范圍.在變程范圍之內，數據具有相關性；而在變程之外，數據之間互不相關，即在變程以外的觀測值不對估計結果產生影響.塊金值(nugget） ：變差函數如果在原點間斷，在地質統計學中稱為“塊金效應”，表現為在很短的距離內有較大的空間變異性，無論h多小,兩個隨機變量都不相關 .它可以由測量誤差引起，也可以來自礦化現象的微觀變異性。在數學上，塊金值c0相當于變量純隨機性的部分

13、.基臺值(sill):代表變量在空間上的總變異性大小.即為變差函數在h大于變程時的值，為塊金值c0和拱高cc之和。 拱高為在取得有效數據的尺度上，可觀測得到的變異性幅度大小。當塊金值等于0時,基臺值即為拱高。 已知井位數據，為了計算某個方向的實驗變差函數，通常計算該方向上的若干不同距離的實驗變差函數值。此時選取某個距離為計算實驗變差函數的基本距離l，稱之為步長，分別計算l， 2l,3l，.。ml距離的實驗變差函數值，此時就可以得到m個實驗變差函數點，m稱之為步長個數.單個步長可選為指定方位上的平均井距,步長數目15-25個。 實際情況下是不可能在精確的某個方向上,或精確的某個步長上能獲得需要的

14、點對數目來計算實驗變差函數值,此時。一個變通的方法就是給方向一個容許的范圍，稱之為角度容限，同時，可以給每個步長一個容許范圍,稱之為步長容限，只要點對記錄落入到該容限內就可以認為該點對參與計算。這樣就可以計算出某個方向相應特定步長的實驗變差函數.步長容限可選為1/2該方位上的井距，角度容限可選為/8, 值得注意的是，隨著步長增加，雖然有時點對符合方向容限與步長容限,但是偏差將會增大，為此需要用偏離方向主線的一個固定寬度的帶子來限制，使超出該范圍的點對不參與計算，這樣有利于得到合理實驗變差函數計算值，帶寬可選為2倍井距。 在擬合過程中，對于短距離的實驗變差函數值應予以重視，是擬合模型形狀盡量照顧

15、到這些點。因為這些點所提供的信息對空間變異型及預測結果有重要影響。另外注意大多數情況下,由于單個步長一般不會太小，從而對選取塊金常數不能提供精確的信息。通常的做法是，若研究目標為區域上的物性參數變化情況，那么小的塊金常數會告訴我們該物性參數具有很好的連續性。因此多數情況下，可置塊金常數為0。 變程的選取依據是適當步長大到何值時，實驗變差函數值基本穩定在某個值范圍鄰近擺動，這個值就是拱高，而相應的步長就是變程.數據分析中提供3種分布模型球狀模型接近原點處,變差函數呈線性形狀，在變程處達到基臺值。原點處變差函數的切 線在變程的2/3處與基臺值相交。指數模型變差函數漸近地逼近基臺值。在實際變程處，變

16、差函數為0.95c。模型在原點處為直線高斯模型變差函數漸近地逼近基臺值。在實際變程處，變差函數為0。95c.模型在原點處為拋物線。facies modeling 相建模選擇數據分析的相,分zones,選擇算法,一般用序貫指示模擬算法,利用數據分析的結果，并加上每層的沉積相分布圖加以約束，相關系數取值？幾種不同的算法：克里金方法是一種實用的、有效的插值方法,它優于三角剖分法、距離反比加權法，它不僅考慮到被估點位置與已知數據位置的相互關系，而且還考慮到已知點位置之間的相互關系，因此估值精度相對較高. 指示克里金插值-是對原始數據的指示變換值進行克里金估計.所謂指示變換是指按照一定的原則將原始數據變

17、換為0或1的過程。屬于非參數統計方法,不同于其他克里金方法（通過參數的均值和方差來估值），它是以概率形式考慮特異值得存在，在不舍棄特異值數據的條件下進行有效的空間估計。 序貫模擬的總體思路是沿著隨機路徑序貫求取各節點的累積條件分布函數ccdf,并從ccdf中提取模擬值。其中用于求取ccdf的條件數據不僅包括原始的樣品點，還包括已模擬好的點。概算法的主要目的是充分利用更多的條件數據來恢復變量的空間相互性。 序貫指示模擬-無需假設原始樣本服從正態分布,而是通過給出一系列的門檻值,估計某一類型變量低于某一門檻值的概率，以此確定隨機變量的分布。實際上是應用指示克里金求取ccdf的序貫模擬方法,主要特點

18、是變量的指示變換、指示克里金和序貫模擬算法.5. 屬性建模 scale up well logs離散化對屬性體數據進行離散化，一般建模過程是poreperm-soil-ntg data analysis 數據分析選擇要分析的屬性體數據，分層，分相分析transformations （數據轉換分析）選擇屬性體后,會出現數據的分布圖,觀察數據的分布，有異常值，在輸入削截和輸出削截中調節數值范圍控制數據的分布 輸入削截的數據不參與運算，輸出削截的數據參與運算，而不向外輸出？其他參數的調節variograms（變差函數分析）和相的變差函數分析一樣,調節主次，垂向變程,調節物源方向，選擇分布模型，擬合分

19、布曲線，便于下步操作。 petrophysical modeling 屬性建模選擇要進行建立模型的屬性體數據，選擇zones，在單zone里分相處理,例如對于ng3-1，存在泥巖，砂巖和差砂巖三個相,用序貫指示模擬的方法計算，采用數據分析的結果。在協-克里金方法下設置相關屬性體或相關平面.其他設置6. 儲量計算make contacts 建立油水界面,油氣界面。劃分界面時,要根據不同的斷塊，不同的層進行劃分，每個層的界面應該都不一樣,每個斷塊的也是一樣.volume calculation 儲量計算操作界面 選擇所在油藏的含有相，以南堡為例，有油氣水三相，因此在前期建立的就應有油水和油氣界面，

20、gen。props中選擇相關屬性體ntg/por,并設置相關系數。oil.props中設置含油飽和度，含氣飽和度（三相時設置），并且設置油的體積系數和溶解氣油比（solution gas_oil radio）,并且設置采收率系數？gas。props 在有氣相存在時,設置氣的體積系數和溶解油氣比(vapourized oil_gas radio)，設置采收率系數？result/output 中進行輸出內容設置，在report setting中可以設置輸出報告的內容以及輸出結果的單位。 /facies是怎么應用的？/boundaries 中可設置工區邊界？井邊界?幾個相關公式:bulkvolume

21、vn net=bulkvolume n/gv por=netvolume * porosityhcpvo= porvolume * sohcpvg= porvolume sgstoilp=hcpvo/bo + (hcpvg/bg） * rvgilp=hcpvg/bg + （hcpvo/bo) rsrecoverable oil=stoilp recfo(oil recoverable factor)recoverable gas = gilp recfg(gas recoverable factor)經過儲量計算之后,得出一分報告，在報告中，會給出軟件名，作者，制作時間，項目名稱，模型名稱以及

22、所含相;作者的bo,bg,rs，rv，recfo，recfg,so，sw，sg等參數設置都會列表顯示出來。列表顯示上面每個公式所得到的值，最后得到可采儲量。而后，分層zones，分斷塊segments，顯示結果.7. 井設計及相關預測well design畫一條polygan。新成一個剖面，顯示pore模型,add new points，在預計井所過的位置加點設計軌跡。在完成設計的井軌跡之后，在input/wells中多出propsed wells ,打開，雙擊選擇make logs tab，就會生成預測曲線；選擇reporttab/iconnize point as :horizon，選擇網格體,按execute就可以整理出分層數據，同時出一份報告.8. 出圖及相關報告所需資料繪制平面圖map window繪制截面圖intersection windowfunct