jason地質(zhì)統(tǒng)計(jì)學(xué)反演手冊

1、StatMod MC入門手冊Chapter 1. 工作流程Stage 1: 項(xiàng)目準(zhǔn)備.5%Stage 2: 地質(zhì)統(tǒng)計(jì)學(xué)參數(shù)分析.15%Stage 3: 模擬.5%百分?jǐn)?shù)表示每個(gè)步驟所用時(shí)間占整個(gè)項(xiàng)目時(shí)間的百分比Stage 4: 反演.50%.10%Stage 5: 協(xié)模擬.15%Stage 6: 不確定性分析與風(fēng)險(xiǎn)評估Chapter 2.基本的輸入輸出數(shù)據(jù)輸入數(shù)據(jù)輸出數(shù)據(jù)巖性實(shí)現(xiàn)巖性概率體屬性實(shí)現(xiàn)地質(zhì)統(tǒng)計(jì)學(xué)參數(shù)巖石物理分析地層網(wǎng)格模型地震數(shù)據(jù)測井曲線Chapter 3. 詳細(xì)操作步驟操作步驟以StatMod MC培訓(xùn)數(shù)據(jù)為例第一步.首先完成一個(gè)高質(zhì)量的疊后CSSI反演這一步的目的是為地質(zhì)統(tǒng)計(jì)學(xué)

2、提供一個(gè)好的研究基礎(chǔ), 這個(gè)“好”主要體現(xiàn)在：(1) 好的井震標(biāo)定, 目標(biāo)區(qū)的相關(guān)值達(dá)到0.85以上;(2) 好的疊后反演結(jié)果, 用來質(zhì)控地質(zhì)統(tǒng)計(jì)學(xué)模擬和反演結(jié)果, 是地質(zhì)統(tǒng)計(jì)學(xué)反演結(jié)果橫向預(yù)測準(zhǔn)確度的參照物;(3) 利用疊后反演結(jié)果進(jìn)行砂體雕刻, 對目標(biāo)區(qū)的巖性展布、比例有一個(gè)總體上正確的把握, 這些認(rèn)識都是地質(zhì)統(tǒng)計(jì)學(xué)的初始輸入。（說明：在提供的培訓(xùn)數(shù)據(jù)中已經(jīng)為用戶做了以上準(zhǔn)備，用戶可以從主界面中打開該培訓(xùn)數(shù)據(jù)所在工區(qū), 然后用Map View看工區(qū)底圖，用Section View查看地震數(shù)據(jù)、疊后CSSI反演數(shù)據(jù)、地質(zhì)框架模型與層位數(shù)據(jù)以及井?dāng)?shù)據(jù)與子波 , 并用Well Editor檢查井

3、震標(biāo)定情況）第二步. 數(shù)據(jù)準(zhǔn)備 井曲線重采樣這一步將測井?dāng)?shù)據(jù)重采樣至地質(zhì)微層采樣間隔，具體操作為：(1) JGW主界面 Analysis Processing toolkit；(2) Input Data selection Data type：選Well, 點(diǎn)擊Input file(s)右邊List選擇任意井(可以選多井),然后在彈出的界面Select logs中選擇任意井曲線(可以多選),點(diǎn)擊OK退出;(3) Parameters Resample log, 在彈出界面Processing toolkit中填寫重采樣間隔(注意s 與ms單位), 點(diǎn)擊OK退出;(4) Output Defi

4、ne process, 從Select from中選擇Resample log, 點(diǎn)擊輸入到右邊的Process里面;(5) Output Generate, 在彈出的界面中填寫輸出路徑和輸出文件名,然后點(diǎn)擊Generate,開始計(jì)算重采樣的曲線。 計(jì)算巖性曲線 ：這里說的巖性曲線和測井解釋的巖性不同, 而是指能通過彈性屬性(主要指縱波阻抗) 區(qū)分出來的巖性, 所以一般可分兩到三種巖性。練習(xí)里通過交會圖的方式教授如何劃分巖性。 地震數(shù)據(jù)：線、道間距盡量一致，如果不一致需要進(jìn)行跳道處理，避免往反演結(jié)果中引入各向異性。 地質(zhì)框架模型：微層采樣間隔與需要分辨的薄層厚度要匹配。第三步. 地質(zhì)統(tǒng)計(jì)學(xué)參數(shù)

5、分析這里說的地質(zhì)統(tǒng)計(jì)學(xué)參數(shù)主要指三個(gè)參數(shù)：概率密度函數(shù)(probability density function, 簡稱pdf，描述某一屬性在空間的概率分布情況)、變差函數(shù)（描述某一屬性隨距離的變化，是距離的函數(shù)）、云變換（描述兩個(gè)屬性之間的相關(guān)關(guān)系）。下面先介紹前兩個(gè)參數(shù): 概率密度函數(shù)(1) JGW主界面 Modeling StatMod MC, 彈出StatMod MC主界面，主要填寫內(nèi)容如下：圖 1 Solid Model： 點(diǎn)擊Input，再點(diǎn)擊Solid model，選擇地質(zhì)框架模型，如solidModel.smgLayers ： 點(diǎn)擊Input，再點(diǎn)擊Layers, 選擇相應(yīng)目標(biāo)

6、層，如TopPad、Reservoir、BasePad, 利用輸入至右邊的界面中,并且修改Spacing值,一 般目標(biāo)層厚度設(shè)為與需要分辨的薄層時(shí)間厚度值，而非目標(biāo)層可以設(shè)置得稀一些，如與地震數(shù)據(jù)采樣一致 Trace Gate ： 選擇反演道門，如wholeProject QC Trace Gate ： 選擇QC道門，如allWellsLine 填好后如圖1所示。(2) StatMod MC主界面Mode, 共有四個(gè)選項(xiàng)：Generate Stats、Simulation、Inversion、Cosimulation，反映的其實(shí)就是整個(gè)流程。在這一步里主要用Generate Stats來分析p

7、df函數(shù)（即概率密度函數(shù)）。點(diǎn)擊Generate Stats，彈出界面（如圖2所示）。 圖2 界面上部由空白欄和右邊的控鍵組成。根據(jù)分析的需要，選擇不同的控鍵往空白欄處充填內(nèi) 容。在這一步只需用到最上邊的一個(gè)控鍵Base/Joint，表示分析某一屬性的pdf函數(shù)；界面中部選擇Output directory填寫輸出路徑；界面下方選擇Generate and run產(chǎn)生腳本并運(yùn)行；現(xiàn)在開始分層位分巖性分屬性地填寫空白欄內(nèi)容： TopPad層：點(diǎn)擊Base/Joint, 彈出小界面，如圖3所示。圖3點(diǎn)擊Layer對應(yīng)空白欄右邊的三角下拉菜單，選擇相應(yīng)層位，如TopPad；點(diǎn)擊Discrete Pr

8、operty對應(yīng)空白欄右邊的三角下拉菜單,選擇相應(yīng)離散屬性,一般指巖性屬性（對于非目標(biāo)層，這一項(xiàng)可以不選）；點(diǎn)擊Continuous Property右邊的Edit,選擇P-Impedance；在Wells一行點(diǎn)擊Edit選擇重采樣的井路徑，然后在Wells in current directory中選擇相應(yīng)井，輸入至Selected wells；點(diǎn)擊OK退出編輯狀態(tài)。這時(shí)在Generate Stats主界面中就會出現(xiàn)相應(yīng)的內(nèi)容。 BasePad層：同TopPad層。 Reservoir層：點(diǎn)擊Base/Joint，Layer選Reservoir, Discrete Property選定義好的

9、巖性曲線, 點(diǎn)擊Continuous Property右邊的Edit,選擇P-Impedance。選井。填寫完后Generate Stats界面空白欄顯示如圖4所示。 圖4還可以通過選擇界面中的Copy，選中所有的行，最后點(diǎn)擊Paste，選擇Paste wells，可以將井復(fù)制到選中的位置；最后在Generate Stats界面下半部分的Output directory 中填寫相應(yīng)的路徑，然后點(diǎn)擊Generate and run 產(chǎn)生腳本并運(yùn)行。(3) 查看分析結(jié)果。在Linux系統(tǒng)中打開網(wǎng)絡(luò)流覽器Konqueror, 并在Location中輸入相應(yīng)的腳本地址。在該路徑下除了腳本還有很多運(yùn)行腳

10、本產(chǎn)生的重要文件：Ø Index.html文件：里面是分層位分巖性分屬性分析得到的所有pdf圖件。以Reservoir層nopay巖性中P-Impedance的分析為例，如圖5所示。分析井段穿過Reservoir層nopay巖性的樣本點(diǎn)統(tǒng)計(jì)結(jié)果Reservoir層中nopay巖性中的P-Impedance的pdf函數(shù)分析結(jié)果圖5圖5中顯示了Reservoir層中nopay巖性中P-Impedance的pdf函數(shù)分布。其它分析結(jié)果與之相似，為了直觀地了解分析結(jié)果，用戶應(yīng)該仔細(xì)檢查每一層每一巖性中各種屬性的pdf函數(shù)分析結(jié)果圖件。Ø *.png文件：一系列*.png文件, 實(shí)際

11、是把Index.html文件里的圖件分開來顯示, 一個(gè)*.png文件就是一張圖, 這些文件方便用戶在做報(bào)告或文檔時(shí)直接插入感興趣的圖件, 而不用再去抓圖。Ø scriptCopy.py文件：就是根據(jù)在Generate Stats界面中填寫的參數(shù)產(chǎn)生的腳本文件，有點(diǎn)類似程序語言，但是相對簡單多了。用戶應(yīng)該逐漸熟悉腳本文件，先嘗試閱讀腳本，領(lǐng)會每部分的意思。在逐漸熟悉后可以嘗試修改腳本或自己編寫。足夠熟練后甚至可以跨過填寫界面這一步，直接編寫腳本并運(yùn)行，來完成地質(zhì)統(tǒng)計(jì)學(xué)參數(shù)分析、模擬、反演等工作。圖6是該工區(qū)地質(zhì)統(tǒng)計(jì)學(xué)參數(shù)分析腳本例子。輸入地層框架模型輸入分析道門準(zhǔn)備信息輸入QC道門輸入

12、分析層段以及設(shè)定相應(yīng)的采樣間隔BasePad層為非目標(biāo)層，對其中的P-Impedance不用分巖性分析開始逐層分析Reservoir層為目標(biāo)層，對其中的P-Impedance應(yīng)該分巖性進(jìn)行pdf函數(shù)的分析paynopay 圖6 Ø stats.py文件：以上分析結(jié)果的具體數(shù)據(jù)都存放在該文件內(nèi)，可以說是最終的結(jié)果文件，以后在做模擬與反演的過程中調(diào)用的就是該文件。以Reservoir層的nopay巖性為例顯示該文件的部分內(nèi)容，如圖7所示。Reservoir層nopay巖性的P-Impedance的pdf函數(shù)的分析結(jié)果，以函數(shù)（Q,X）表示，含義為P-Impedance小于X的概率為Q圖7

13、l 變差函數(shù)的求取(1) JGW主界面 Modeling Specials StatMod StatMod Analysis, 彈出StatMod Analysis界面；(2) Input Time/Depth mode: Time Solid model ： SOLID Layers ：點(diǎn)擊Select layers選擇目標(biāo)層位，如Reservoir層，在Selected interval設(shè)置采樣間隔為0.001s，點(diǎn)擊Assign current interval to selected layers應(yīng)用修改的采樣間隔 Lithology masks ：先把Enable lithology

14、 masking左邊的小方框點(diǎn)亮，再點(diǎn)擊Lithology右邊的List，在彈出的界面中選擇事先劃分好的巖性曲線，如Lithology, 這時(shí)在Select lithotypes to use中選擇某一巖性，如pay（其它巖性稍后也需要分析，只是一次只分析一種巖性中的屬性變差） Data for histograms and transforms： 在彈出界面的上半部Data的下方的Well file(s) 右邊List中在./WELLS/RESAMPLED選擇相應(yīng)的分析井曲線，如P-Impedance。 Trace gate : wholeProjectEdit Transforms, 彈出

15、兩個(gè)界面，一個(gè)是主界面StatMod Transform，一個(gè)是輔界面Transform, 點(diǎn)擊輔界面Transform最上方的Histogram, 彈出一個(gè)小界面StatMod Histogram 并帶著一個(gè)不停閃爍的紅色提示Build the histogram, 點(diǎn)擊該紅色提示，則主界面出現(xiàn)相應(yīng)的樣本點(diǎn)分析，調(diào)試小界面中的Nr of intervals的值, 盡量使pdf函數(shù)樣本點(diǎn)分布規(guī)則連續(xù)。調(diào)試好后點(diǎn)擊Dismiss退出。點(diǎn)擊輔界面Transform上方第二個(gè)選項(xiàng)Transform, 彈出界面如圖8，同時(shí)在主界面中出現(xiàn)一條藍(lán)色擬合線，如圖9。用鼠標(biāo)左鍵調(diào)整藍(lán)色擬合線，使得該曲線形態(tài)能

16、大致反映樣本點(diǎn)的分布形態(tài)，最后用ctrl+鼠標(biāo)左鍵停止編輯藍(lán)色擬合線。然后在圖8顯示的界面下方選擇Save as first, 在彈出界面中給出名字，如pay_p-imp_histogram，再連續(xù)點(diǎn)三次OK回到StatMod Transform主界面, 點(diǎn)擊Dismiss退出界面。 輔界面主界面 圖8 圖9 (3) Input Data for variogram sampling and modeling Primary data Well log data, 在 彈出界面中的上部Data部分的Well file(s)右邊的List中選擇相應(yīng)的分析井曲線，如P-Impedance, 然后點(diǎn)

17、OK, 這時(shí)彈出兩個(gè)界面，一個(gè)StatMod Variogram主界面，一個(gè)Variogram輔界面。點(diǎn)擊 Variogram 輔界面上方第一個(gè)選項(xiàng)Sample Vario, 彈出界面StatMod Sample Variogram Parameters, 如圖10所示，先點(diǎn)擊在該界面下方不停閃爍的紅色提示Build the Sample Variogram, 則主界面中出現(xiàn)變差函數(shù)的樣本點(diǎn)分析曲線，修改圖10界面中的Select Sampling Parameters部分的Intervalm下方空白處的值（即變程的單位距離），則主界面中的變差函數(shù)的樣本點(diǎn)分析曲線形態(tài)發(fā)生改變，不斷調(diào)試Inte

18、rvalm值，使得主界面中的變差函數(shù)的樣本點(diǎn)分析曲線形態(tài)盡量光滑，如圖11所示，主要看上半部的紅色樣本曲線，最后點(diǎn)OK結(jié)束編輯。 輔界面主界面圖 10 圖11 （說明：圖11的主界面顯示的是進(jìn)行變差函數(shù)分析的樣本點(diǎn)曲線，變差函數(shù)是三維的，有縱向變差，也有橫向上的變差，圖11的輔界面下方表示的是將三個(gè)方向上的變差函數(shù)曲線用三種顏色顯示。但是因?yàn)檫@里是利用井上的樣本點(diǎn)進(jìn)行變差函數(shù)的分析，一方面井曲線在縱向上的樣本點(diǎn)的個(gè)數(shù)一般能滿足統(tǒng)計(jì)需要（一般50個(gè)樣本點(diǎn)就可以得到比較理想的pdf函數(shù)和變差函數(shù)），另一方面井曲線在橫向上的分布是遠(yuǎn)不能達(dá)到分析需要的。所以在這里只相信縱向上的變差函數(shù)分析結(jié)果，即只需

19、分析紅色的曲線）點(diǎn)擊Variogram輔界面的選項(xiàng)Model Vario，彈出界面如圖12所示，同時(shí)把鼠標(biāo)箭頭放在圖13的上部，按住左鍵可以得到一條擬合曲線，可以任意拖動(dòng)該擬合曲線，使得它能最好地?cái)M合縱向變差樣本點(diǎn)，特別是前幾個(gè)樣本點(diǎn)（思考為什么？），ctrl+鼠標(biāo)左鍵停止編輯。 圖12 圖13這時(shí)圖12顯示界面的上部空白欄處會出現(xiàn)變差函數(shù)表達(dá)式（Type, Sill, X, Y, Z, Power）對應(yīng)的具體參數(shù)。在本例中表示變差函數(shù)類型是指數(shù)型的，基臺值（Sill,即擬合曲線的拐點(diǎn)值）為3.07, 三個(gè)方向(x,y,z)的變程均為26.9m（當(dāng)然就像之前所說的只相信縱向上的變程分析結(jié)果）。

20、另外，還可以改變變差函數(shù)類型，在圖12界面中上部分，F(xiàn)unc Type下方，點(diǎn)擊小控鍵，彈出幾種類型，一般用到的為Gaussian、Expo兩種類型。在原點(diǎn)附近，Gaussian型變化較慢，Expo型變化較快，用戶可以分別嘗試。如果不滿意分析結(jié)果，可以在圖12中選中上方空白欄里的函數(shù)，然后點(diǎn)擊下面的Remove, 就可以刪掉該函數(shù)，然后重新編輯。最后記住測試好的變差函數(shù)類型與縱向上的變程（單位為深度域單位m）(4) 依次對Reservoir層的其它巖性類型中的屬性（如P-Impedance）進(jìn)行變差函數(shù)分析，并記住測試好的變差函數(shù)類型與縱向上的變程。再對其它目的層進(jìn)行分巖性的屬性變差函數(shù)測試。

21、過程同上。(5) 在分層確定了不同巖性的變差函數(shù)類型與縱向變程后, 橫向變程怎么確定呢? 這里提供兩中思路: 一是與地質(zhì)概念相結(jié)合, 與地質(zhì)專家或客戶進(jìn)行交流, 大體掌握工區(qū)的沉積模式與巖性展布, 以此來確定一個(gè)大概的橫向變差范圍(只要求一個(gè)大概的輸入)；一是依據(jù)疊后CSSI反演結(jié)果進(jìn)行巖性的地質(zhì)規(guī)模和分布等地分析, 來得到一個(gè)大概的橫向變差范圍。(6) 不管是橫向還是縱向上的變差，經(jīng)過分析得到的都是一個(gè)大概的值，對反演的結(jié)果起到的是軟約束的作用，而起到硬約束影響作用的是地震。第四步.地質(zhì)統(tǒng)計(jì)學(xué)模擬-目的是測試第三步中確定的地質(zhì)統(tǒng)計(jì)學(xué)參數(shù)是否合理(1) StatMod MC主界面 Mode S

22、imulation, 彈出Simulation主界面，該界面上部和Generate Stats界面基本一致。以模擬目標(biāo)層Reservoir中不同巖性的P-Impedance為例來介紹需要填寫的內(nèi)容。點(diǎn)擊Base/Joint，在彈出的小界面中Layer選Reservoir層，Discrete Property選巖性曲線, Continuous中選P-Impedance, OK退出。這時(shí)Simulation主界面中出現(xiàn)相關(guān)信息。如圖14所示。對界面中的內(nèi)容進(jìn)行編輯：包括巖性比例、巖性的變差函數(shù)、不同巖性中P-Impedance的pdf函數(shù)與變差函數(shù)，也就是前面幾步分析得到的結(jié)果在這里都要用到。選井

23、。圖14相關(guān)參數(shù):Ø 巖性比例：選中Reservoir層第一行，如圖14所示，然后點(diǎn)擊Edit，彈出界面如 圖15所示。在Proportions部分編輯巖性比例（ 說明： 巖性比例缺省值為1.0，本練習(xí)工區(qū)巖性比例設(shè)置可以參考pay ：nopay =13 ：87, 可以從井上目標(biāo)層段不同巖性的樣本點(diǎn)統(tǒng)計(jì)得到一個(gè)粗略的估計(jì)，但仍需測試）Ø 巖性的變差函數(shù)：在圖15中的Variograms部分, 點(diǎn)擊Single右邊的Edit，在彈出菜單中依 次 將Type、Vertical range、Lateral range設(shè)定為gaussian、0.010、600×600。(

24、 說明： Type、Vertical range均為在第三步中測試到的值, Lateral range可以由第一步得到,即利用疊后反演結(jié)果進(jìn)行砂體雕刻, 以求對目標(biāo)區(qū)的巖性展布、比例有一個(gè)總體上正確的把握, 同時(shí)結(jié)合工區(qū)的地質(zhì)沉積模式等地質(zhì)方面的信息來指導(dǎo)Lateral range的輸入。用戶可以放心的是Lateral range的輸入可以只是一個(gè)大概的值，它對真實(shí)的各種巖性體的分布起到的是軟約束的作用，在后面的反演部分要講到，對巖性體的分布起硬約束作用的是地震數(shù)據(jù)。并且Type、Vertical range、Lateral range的設(shè)定也不是一蹴而就的，需要調(diào)試，尋找合理值的范圍)編輯巖

25、性比例編輯屬性的變差函數(shù)編輯巖性變差函數(shù)圖15 圖16Ø 不同巖性中P-Impedance的pdf函數(shù)與變差函數(shù)：在Simulation主界面中選中Reservoir目 標(biāo)層某一巖性對應(yīng)的行，如選擇Lithology.nopay所在行，然后點(diǎn)擊右邊控鍵欄Property下方的Edit，彈出界面如圖16所示，只需對紅框圈出的部分進(jìn)行編輯。點(diǎn)擊Histogram對應(yīng)的Edit，選擇相應(yīng)的存放pdf函數(shù)分析結(jié)果的文件。大家還記得這個(gè)文件在哪嗎？就在第三步做地質(zhì)統(tǒng)計(jì)學(xué)參數(shù)分析時(shí)，在Generate Stats界面運(yùn)行腳本產(chǎn)生的stats.py文件，里面存放所有目標(biāo)層中不同巖性中的P-Impe

26、dance等屬性的pdf函數(shù)分析結(jié)果。圖7展示stats.py文件的部分內(nèi)容。大家需要好好回憶一下，整理好頭緒。然后點(diǎn)擊Variogram對應(yīng)的Edit, 在彈出的小界面中依此設(shè)定Type、Vertical range、Lateral range: exponential、0.006、400×400。這個(gè)結(jié)果也是由第三步的分析得到的。連續(xù)兩次點(diǎn)擊OK退出編輯狀態(tài)。編輯其它巖性中屬性的pdf函數(shù)與變差函數(shù)（其實(shí)pdf函數(shù)調(diào)用的都是stats.py文件，只需編輯變差函數(shù)即可）。變差函數(shù)值可參考以下表1（例子）： VariogramLithologyTypeVertical rangeLa

27、teral rangepayexponential0.0161200×1200nopayexponential0.0161200×1200 表 1在編輯參數(shù)時(shí)可以使用主界面中的Copy、Paste功能，實(shí)現(xiàn)Histogram與Variogram的快速編輯。界面中下方的Run Parameters部分也有一些選項(xiàng)，在Nr of realizations可以填寫想要產(chǎn)生實(shí)現(xiàn)的個(gè)數(shù)，Random seed中的值隨之改變。如Nr of realizations填1，表示只產(chǎn)生一個(gè)實(shí)現(xiàn)，這時(shí)Random seed自動(dòng)將種子點(diǎn)設(shè)為0；如Nr of realizations填2， 表示產(chǎn)

28、生兩個(gè)實(shí)現(xiàn)，Random seed自動(dòng)將種子點(diǎn)設(shè)為0,1。 在實(shí)際的工作中用戶可能需要產(chǎn)生更多的實(shí)現(xiàn)。用戶可以自己熟悉其它選項(xiàng)。最后別忘了給出輸出路徑與輸出文件名。點(diǎn)擊Generate script產(chǎn)生腳本文件，可以從網(wǎng)絡(luò)流覽器中調(diào)出腳本文件查看，如圖17展示腳本中部分內(nèi)容。最后點(diǎn)擊Run運(yùn)行該腳本文件。 （留住Simulation主界面） nopayReservoir層中每一巖性中P-Impedance的變差函數(shù) pay圖17 (2) 查看結(jié)果并質(zhì)控： 打開網(wǎng)絡(luò)流覽器，查看運(yùn)行腳本產(chǎn)生Index.html文件，選擇該文件，然后在Run report下面選擇QC histograms(deta

29、iled), 里面是所有反映模擬信息的圖件和輸入輸出參數(shù)統(tǒng)計(jì)表,圖18展示其中部分內(nèi)容。nopay中輸出P-Impedance的pdf函數(shù)的主要參數(shù)巖性比例標(biāo)準(zhǔn)偏差平均值nopay中輸入P-Impedance的pdf函數(shù)的主要參數(shù)標(biāo)準(zhǔn)偏差平均值nopay中P-Impedance輸入輸出的pdf函數(shù)曲線pay中P-Impedance輸入輸出的pdf函數(shù)曲線圖 18 在該圖件文件里面主要檢查兩方面內(nèi)容：目標(biāo)層中每一種巖性內(nèi)相應(yīng)屬性，如P-Impedance 輸入輸出的pdf函數(shù)形態(tài)應(yīng)該大體一致，主要比較入輸出pdf函數(shù)的平均值與標(biāo)準(zhǔn)偏差（圖18中藍(lán)色曲線表示輸入，紅色曲線表示輸出）；目標(biāo)層輸入輸出

30、的巖性比例應(yīng)該大體一致（偏差在5%左右）。 JGW主界面 Analysis Section View, 先選好tracegate (最好選擇在StatMod MC主界面中填入的QC Trace Gate：本例中是MyQCGate), 然后從剛才填寫的模擬結(jié)果輸出路徑中調(diào)出結(jié)果查看。打開這個(gè)結(jié)果的方式主要有兩種：Section View界面 File Open session 進(jìn)入模擬結(jié)果輸出路徑 選擇QC.jvis，這時(shí)顯示模擬結(jié)果，包括巖性體與P-Impedance體；Section View界面 Input Stratigraphy data 進(jìn)入模擬結(jié)果輸出路徑 選擇lithology.

31、hor/P-Impedance查看。質(zhì)控方法與標(biāo)準(zhǔn)：和疊后反演的結(jié)果進(jìn)行比較，如圖19所示。與疊后CSSI反演的結(jié)果（右圖），在反映砂體的展布、地質(zhì)體規(guī)模、巖性比例等方面大體比較一致，這時(shí)的模擬結(jié)果是比較合理的。記住相應(yīng)的模擬參數(shù)。 圖19 比較輸入輸出的井曲線。在運(yùn)行腳本的過程中，自動(dòng)在相應(yīng)的目錄下產(chǎn)生在井點(diǎn)處抽取的“偽井曲線”（Lithology與 P-Impedance曲線）。檢查這些偽井曲線，Jason主界面AnalysisWell Manage中選擇相應(yīng)路徑下的井，會發(fā)現(xiàn)對應(yīng)Lithology與P-Impedance，分別有Lithology_in, Lithology_out與P-

32、Impedance_in, P-Impedance_out。*_in表示輸入的重采樣曲線，而*_out表示輸出的從模擬結(jié)果中的井點(diǎn)處抽取出來的巖性或?qū)傩郧€，如圖20所示，藍(lán)色表示輸入，紅色表示輸出，兩者偏差很小，這是比較理想的模擬結(jié)果。圖20 建議多測試幾個(gè)參數(shù)，如改變目標(biāo)層巖性比例、變差函數(shù)的類型和縱橫向的變程范圍來產(chǎn)生實(shí)現(xiàn)，并比較差別，尤其注意在Section View / Map View中與疊后反演結(jié)果進(jìn)行比較，選擇令實(shí)現(xiàn)與疊后反演結(jié)果最相似的參數(shù)組合；在比較剖面的同時(shí)注意檢查生成Index.html文件，比較標(biāo)準(zhǔn)與上面所說的一致。建議做一個(gè)列表來紀(jì)錄試驗(yàn)的參數(shù)：Proportion

33、 ValueDP Vario ValueP-impedance Value VarioOutput directorywell-basedGaussian-shortExponential-longwell-basedGaussian-longExponential-longLowerPayGaussian-shortExponential-longLowerPayGaussian-longExponential-longHigerPayGaussian-shortExponential-longHigerPayGaussian-longExponential-longwell-basedGa

34、ussian-longPay: Exp-longNonpay: Exp-short 表 2 第五步.地質(zhì)統(tǒng)計(jì)學(xué)反演l 無井約束反演-主要目的是測試地震的權(quán)重(1) StatMod MC主界面 Mode Inversion, 彈出Inversion主界面。用戶這時(shí)可能已經(jīng)發(fā)現(xiàn)Inversion主界面與Generate Stats界面、Simulation主界面基本一致。區(qū)別大的地方主要在界面的左下方的Seismic Belief部分。在該部分需要填寫地震數(shù)據(jù)的權(quán)重。另外，由于是反演，所以需要在目標(biāo)層上下加個(gè)“蓋子”，即加入TopPad層和BasePad層作為邊界層（邊界層厚度要在半個(gè)子波長度以上

35、，以消除反演時(shí)子波的旁瓣效應(yīng)）。這一步得回到StatMod MC主界面選擇，把TopPad層和BasePad層的Spacing設(shè)為0.004s。打開Inversion主界面，對于目標(biāo)層的參數(shù)填寫與在Simulation主界面的填寫一致，當(dāng)然也可以用測試效果更好的參數(shù)代替。對于非目標(biāo)層，不用分巖性反演屬性。不要選井。(2) 現(xiàn)在來考慮地震權(quán)重。在Inversion主界面左下方的Seismic Belief部分，有Seismic和Wavelet兩個(gè)控鍵, 點(diǎn)擊Seismic，則彈出選項(xiàng)如圖21所示。 圖21 在StatMod MC反演模塊里，地震的信息是通過設(shè)定一定的噪音引入的，或者更好理解一些，

36、可以通過設(shè)定信噪比引入。信噪比計(jì)算公式如下所示：S/N Ratio (in dB) = 20 * lg (std (seismic)/ std (residual)其中，lg表示10為底的對數(shù)；Std表示標(biāo)準(zhǔn)偏差；在該公式中，可以把殘差residual理解為噪音（即地震數(shù)據(jù)與合成記錄的差值，或者說是合成記錄與地震數(shù)據(jù)匹配不上的部分），地震數(shù)據(jù)是一定的，信噪比越高，表示殘差越?。仍胍粼饺?，地震數(shù)據(jù)與合成記錄越相似）；反之亦然。所以可以通過控制信噪比的大小來控制引入的噪音。點(diǎn)擊Seismic對應(yīng)的Edit，選擇全疊加地震數(shù)據(jù)，然后在Noise一欄，首先點(diǎn)擊Noise右邊下拉菜單，可以看見三個(gè)選項(xiàng)

37、value/dB/file，選擇dB，表示引入信噪比來控制合成記錄與地震的相關(guān)性。然后填入具體的信噪比值。那么如何設(shè)置這個(gè)值呢？可以參考疊后反演時(shí)生成的一個(gè)文件（現(xiàn)在體會到在做地質(zhì)統(tǒng)計(jì)學(xué)反演前需要做疊后CSSI反演的重要性了吧？）。首先從JGW主界面下的Analysis中打開Map View，然后選擇Input Horizon,找到路徑./INVER_TRACE_PLUS, 在這一目錄中有一質(zhì)控的層位文件qc.hor，選擇該文件，然后在彈出的多個(gè)層位中選擇inverted signal to noise ratio(Signal-to-noise ratio), 然后OK，這時(shí)Map View

38、中顯示目標(biāo)區(qū)的疊后反演信噪比結(jié)果，可以打開色標(biāo)對照看，如圖22所示。 圖 22由圖 22可看出本工區(qū)的信噪比范圍大約在525dB之內(nèi)，且工區(qū)中大部分的信噪比在10dB以上。所以我們可以根據(jù)疊后反演的信噪比來設(shè)置地質(zhì)統(tǒng)計(jì)學(xué)反演的信噪比，它不是某一確定的值，而是要在某一數(shù)值范圍內(nèi)進(jìn)行測試。對于本工區(qū)建議測試5dB、10 dB、15 dB、20 dB幾組值。也就是要產(chǎn)生四個(gè)腳本并運(yùn)行。當(dāng)然，如果用戶感興趣可以測試更多值，并感受隨之帶來的變化。填入地震信息后，一定別忘了還需填入子波信息。點(diǎn)擊圖24中的Wavelet選項(xiàng)，在彈出界面中缺省是用單一子波反演，即Single，點(diǎn)擊Single對應(yīng)的Edit，

39、選擇相應(yīng)子波（在疊后CSSI反演時(shí)用的最終子波）。在Run Parameters部分，用戶可以自己設(shè)置需要計(jì)算的實(shí)現(xiàn)，因?yàn)檫@步主要測試地震的權(quán)重，所以不建議產(chǎn)生太多實(shí)現(xiàn)（要花時(shí)間）。最后給出輸出路徑與腳本名，再產(chǎn)生腳本并且運(yùn)行即可。(3) 檢查并質(zhì)控結(jié)果。Ø 檢查反演結(jié)果并進(jìn)行抽井檢查：從JGW主界面下的Analysis中打開Section View, 先在Input 的Trace gate中選擇合適的道門，最好是聯(lián)井線，這里選QC Trace Gate。然后選反演數(shù)據(jù)，點(diǎn)擊Input，選擇Stratigraphic data, 或者是在界面上找到圖標(biāo)并點(diǎn)擊,即可查找反演體文件*.h

40、or(到剛才輸出的路徑中查找)并顯示。一個(gè)實(shí)現(xiàn)對應(yīng)兩個(gè)反演結(jié)果：P-Impedance.hor和 lithology.hor。調(diào)出剖面后會發(fā)現(xiàn)地質(zhì)統(tǒng)計(jì)學(xué)的反演的結(jié)果具有很高的分辨率，但是這一結(jié)果是否就是用戶想要的正確的結(jié)果呢？當(dāng)然要進(jìn)行質(zhì)控。這一步最重要的是把井選進(jìn)來，進(jìn)行抽井檢測（這已經(jīng)是Jason在用戶中普及得最深刻的質(zhì)控手段了，不管是疊前還是疊后反演，要判斷反演結(jié)果的好壞最基本的就是進(jìn)行抽井檢查，地質(zhì)統(tǒng)計(jì)學(xué)反演也不例外）。比如你比較反演的巖性體與井上的巖性曲線，那么首先你要看到在井周圍反演的巖性體與井是否是比較吻合的，然后是觀察井間巖性是否有比較好的連通性，好的反演結(jié)果應(yīng)該能夠預(yù)測大部分

41、規(guī)模較大的地質(zhì)體并且對于大多數(shù)細(xì)節(jié)部分（如薄層）也能有很好的反映。如果在這一步得不到好的結(jié)果，那么就算把井加上作為約束，結(jié)果也還是不可靠的。如果結(jié)果不滿意怎么辦？需要繼續(xù)測試輸入?yún)?shù)，如統(tǒng)計(jì)學(xué)參數(shù)，甚至井震標(biāo)定，直到得到滿意的結(jié)果為止。如圖23為抽井檢查反演巖性體的實(shí)例。圖23 Ø 輸入輸出pdf函數(shù)、巖性比例質(zhì)控：不管在統(tǒng)計(jì)參數(shù)測試步驟、模擬步驟還是在反演步 驟，都能在相應(yīng)路徑找到index.html文件。打開網(wǎng)絡(luò)流覽器，在相應(yīng)路徑下選擇該文件，然后在Run report下面選擇QC histograms(detailed)，能看到反映模擬信息的圖件和輸入輸出參數(shù)統(tǒng)計(jì)表，如圖24顯

42、示的是其中部分內(nèi)容-reservoir層中nopay巖性的P-Impedance的pdf函數(shù)、巖性比例輸入輸出的情況。 圖24圖中藍(lán)色線與紅色線分別為輸入輸出的reservoir層中nopay的P-Impedance的pdf函數(shù)，兩者的相似度較高，是比較理想的輸入輸出情況。圖中右邊是描述輸入輸出pdf函數(shù)的參數(shù)表，其中決定pdf函數(shù)形態(tài)的參數(shù)是平均值Mean、標(biāo)準(zhǔn)偏差StDev，用戶應(yīng)該比較輸入輸出的平均值與標(biāo)準(zhǔn)偏差；此外還應(yīng)比較輸入輸出的巖性比例，即圖表中Proportion，兩個(gè)比較的原則都是輸入輸出應(yīng)該大體一致。這部分其實(shí)與模擬步驟比較的內(nèi)容一致。Ø 反演結(jié)果與疊后反演結(jié)果的比

43、較: 在比較的過程中主要注意兩方面內(nèi)容一是巖性比例二是巖性的形態(tài)與連通性。由于使用相同的正演褶積算子來模擬地震數(shù)據(jù)，所以地質(zhì)統(tǒng)計(jì)學(xué)反演的結(jié)果應(yīng)該與疊后反演結(jié)果比較一致（這個(gè)“一致”主要指大套巖性的規(guī)模分布等，而不強(qiáng)調(diào)細(xì)節(jié)，如薄互層），如圖25所示，左邊是地質(zhì)統(tǒng)計(jì)學(xué)反演P-Impedance結(jié)果，右邊是疊后反演結(jié)果，兩者對大套巖性體的規(guī)模尺度及分布反映都比較一致，區(qū)別僅是地質(zhì)統(tǒng)計(jì)學(xué)反演結(jié)果能反映出更豐富的細(xì)節(jié)。 圖25Ø 巖性體沿層切片檢查：比較地質(zhì)統(tǒng)計(jì)學(xué)反演巖性體的沿層切片與地震數(shù)據(jù)的沿層提取的均方根振幅屬性。地質(zhì)統(tǒng)計(jì)學(xué)反演巖性體的沿層切片上的地質(zhì)體形態(tài)應(yīng)該與地震的比較一致。如果不一

44、致，得考慮以下兩個(gè)問題：一是設(shè)置的信噪比過低了；二是選擇的變差函數(shù)與地震信息不符，如變程大于主要地質(zhì)體的尺度與規(guī)模。打開地質(zhì)統(tǒng)計(jì)學(xué)反演巖性體的沿層切片的方式為：打開Map View Input Stratal slicer, 在彈出界面中，對應(yīng)Stratigraphic Data選擇相應(yīng)的反演巖性體，在Layer中選擇相應(yīng)的層，在Microlayer一行用前后翻鍵查看微層中的巖性分布情況。l 無井約束反演-主要目的是學(xué)會產(chǎn)生多個(gè)實(shí)現(xiàn)并進(jìn)行統(tǒng)計(jì)分析(1) 利用界面或腳本產(chǎn)生多個(gè)實(shí)現(xiàn)：如果對腳本不熟悉，那么利用界面產(chǎn)生多個(gè)實(shí)現(xiàn)也是非常簡單的，做法其實(shí)在前面也提到過。在上一步的反演界面中確定最佳地質(zhì)

45、統(tǒng)計(jì)學(xué)參數(shù)與信噪比并填寫好，然后在界面的右下方Run Parameters部分的第一個(gè)選項(xiàng)Nr of realizations填寫想要產(chǎn)生的實(shí)現(xiàn)個(gè)數(shù)，如10（上一步中建議為了節(jié)省時(shí)間，只產(chǎn)生一個(gè)實(shí)現(xiàn)），然后把鼠標(biāo)放在第二個(gè)選項(xiàng)Random seeds對應(yīng)的選項(xiàng)內(nèi)單擊，就會出現(xiàn)0，1，2，3.9，這時(shí)準(zhǔn)備工作就做好了，再點(diǎn)擊Generate script產(chǎn)生腳本，感興趣的用戶可以打開該腳本看看腳本中用來產(chǎn)生多個(gè)實(shí)現(xiàn)的語句是什么？最后點(diǎn)擊RUN運(yùn)行反演。這個(gè)過程可能會花較長時(shí)間，所以白天做測試，下午離開辦公室之前運(yùn)行腳本是最好的選擇。運(yùn)行完畢后，在相應(yīng)目錄下產(chǎn)生10個(gè)路徑：seed_0、seed_

46、1、seed_2. 、seed_9, 每個(gè)路徑下都是一組完整的反演結(jié)果，包括P-Impedance.hor和 lithology.hor，以及index.html文件等，用戶應(yīng)該按照上一步講述的質(zhì)控步驟檢查結(jié)果，并注意每個(gè)實(shí)現(xiàn)之間的差別。(2) 對多個(gè)實(shí)現(xiàn)進(jìn)行統(tǒng)計(jì)分析：現(xiàn)在產(chǎn)生了多個(gè)實(shí)現(xiàn)，并且這些實(shí)現(xiàn)都是等概率的，那么哪個(gè)實(shí)現(xiàn)更合理呢？StatMod MC提供了對多個(gè)實(shí)現(xiàn)進(jìn)行統(tǒng)計(jì)性分析的方法。在反演界面中點(diǎn)擊OK退出，StatMod MC主界面 Utilities Summary statistics, 彈出界面，在Directory對應(yīng)的Edit中選擇相應(yīng)目錄，這時(shí)在中部左邊的方框中出現(xiàn)se

47、ed_0、seed_1、 seed_2 、. 、seed_9十個(gè)實(shí)現(xiàn)，全選然后用輸入到右邊方框中，最后在界面下方的Output directory中填寫輸出路徑（缺省為summary），點(diǎn)擊RUN運(yùn)行，運(yùn)行速度較快。運(yùn)行結(jié)束后在相應(yīng)路徑下查看結(jié)果：對于巖性體來說，分析結(jié)果包括每一種巖性的概率體（比如四種巖性就有四種巖性的概率體）以及一個(gè)極大似然巖性概率體；對于P-Impedance體來說，分析結(jié)果主要有四個(gè)分別是平均值P-Impedance_mean.hor、最小值P-Impedance_min.hor 、最大值P-Impedance_max.hor 、標(biāo)準(zhǔn)偏差體P-Impedance_std

48、dev.hor。用戶可以從相應(yīng)路徑下打開這些結(jié)果來查看。(3) 檢查并質(zhì)控結(jié)果：Ø 比較各個(gè)實(shí)現(xiàn)與疊后反演結(jié)果：主要比較這些實(shí)現(xiàn)（P-Impedance.hor和 lithology.hor）與疊后CSSI反演結(jié)果對大套巖性體的規(guī)模尺度的反映上是否基本一致，而不用太摳細(xì)節(jié)。Ø 對巖性概率體進(jìn)行抽井檢查：在Section View中打開巖性概率體，如此說來圖26所示，并投上井上的巖性曲線進(jìn)行抽井檢查。左圖為泥巖概率體，右圖為氣砂概率體。紅色代表概率為1，即10個(gè)實(shí)現(xiàn)中一般以上都反映有該種巖性；黑色表示概率為0，即所有實(shí)現(xiàn)都沒有這種巖性。比較時(shí)注意兩點(diǎn)：一是在井點(diǎn)處巖性概率體是

49、否與井上的巖性相符，如果不符，那么不確定性有多大？二是井間或遠(yuǎn)離井點(diǎn)處巖性的不確定性有多大？如果在井點(diǎn)處吻合得很好，而井間不具有連通性，即巖性概率降低，那么要考慮這種變化是否與實(shí)際情況相符，如果不相符，那么要考慮加入的信噪比是否合適？選擇的巖性橫向變程是否過??？ 圖 26Ø 平均P-Impedance體與疊后反演結(jié)果的比較：準(zhǔn)則是兩者對大套巖性的規(guī)模尺度反應(yīng)大體一致。Ø 查看極大似然巖性概率體：調(diào)出極大巖性概率體并且投上井的巖性曲線進(jìn)行抽井檢查, 在井周比較極大似然巖性概率體與井上巖性曲線的吻合程度, 在井間檢查巖性的連通性等等。Ø 在井點(diǎn)處比較井與從平均P-Impedance體中抽出的P-Impedance曲線: 從平均P-Impedance體中抽出的P-Impedance曲線要能夠反映井上的P-Impedance的大致趨勢。Ø 輸入輸出pdf函數(shù)、巖性比例質(zhì)控：打開網(wǎng)絡(luò)流覽器，在各個(gè)實(shí)現(xiàn)對 應(yīng)的路徑下選