jason地質(zhì)統(tǒng)計學(xué)反演基礎(chǔ)手冊

StatModMC入門手冊Chapter1.工作流程5%15%5%50%10%15%百分數(shù)表達每個環(huán)節(jié)所用時間占整個項目時間比例Chapter2.基本輸入輸出數(shù)據(jù)輸入數(shù)據(jù)地層網(wǎng)格模型巖石物理分析地質(zhì)記錄學(xué)參數(shù)測井曲線1 —?地層網(wǎng)格模型巖石物理分析地質(zhì)記錄學(xué)參數(shù)測井曲線1 —?輸出數(shù)據(jù)屬性實現(xiàn)巖性概率體屬性實現(xiàn)巖性概率體Chapter3.詳細操作環(huán)節(jié)操作環(huán)節(jié)以StatModMC培訓(xùn)數(shù)據(jù)為例第一步.一方面完畢一種高質(zhì)量疊后CSSI反演這一步目是為地質(zhì)記錄學(xué)提供一種好研究基本，這個“好”重要體當前：好井震標定，目的區(qū)有關(guān)值達到0.85以上；好疊后反演成果，用來質(zhì)控地質(zhì)記錄學(xué)模仿和反演成果，是地質(zhì)記錄學(xué)反演成果橫向預(yù)測精確度參照物；運用疊后反演成果進行砂體雕刻，對目的區(qū)巖性展布、比例有一種總體上對的把握，這些結(jié)識都是地質(zhì)記錄學(xué)初始輸入。(闡明：在提供培訓(xùn)數(shù)據(jù)中已經(jīng)為顧客做了以上準備，顧客可以從主界面中打開該培訓(xùn)數(shù)據(jù)所在工區(qū)，然后用MapView看工區(qū)底圖，用SectionView查看地震數(shù)據(jù)、疊后CSSI反演數(shù)據(jù)、地質(zhì)框架模型與層位數(shù)據(jù)以及井數(shù)據(jù)與子波，并用WellEditor檢查井震標定狀況)第二步.數(shù)據(jù)準備?井曲線重采樣這一步將測井數(shù)據(jù)重采樣至地質(zhì)微層采樣間隔，詳細操作為：JGW主界面一Analysis一Processingtoolkit；Input一Dataselection一Datatype：選Well，點擊Inputfile(s)右邊List選取任意井(可以選多井),然后在彈出界面Selectlogs中選取任意井曲線(可以多選)，點擊OK退出；Parameters—Resamplelog，在彈出界面Processingtoolkit中填寫重采樣間隔(注意s與ms單位)，點擊OK退出；Output一Defineprocess，從Selectfrom中選取Resamplelog，點擊＞＞輸入到右邊Process里面；Output—Generate，在彈出界面中填寫輸出途徑和輸出文獻名，然后點擊Generate，開始計算重采樣曲線。計算巖性曲線：這里說巖性曲線和測井解釋巖性不同，而是指能通過彈性屬性(重要指縱波阻抗)區(qū)別出來巖性，因此普通可分兩到三種巖性。練習里通過交會圖方式專家如何劃分巖性。地震數(shù)據(jù)：線、道間距盡量一致，如果不一致需要進行跳道解決，避免往反演成果中引入各向異性。?地質(zhì)框架模型：微層采樣間隔與需要辨別薄層厚度要匹配。第三步.地質(zhì)記錄學(xué)參數(shù)分析這里說地質(zhì)記錄學(xué)參數(shù)重要指三個參數(shù)：概率密度函數(shù)(probabilitydensityfunction，簡稱pdf，描述某一屬性在空間概率分布狀況)、變差函數(shù)(描述某一屬性隨距離變化，是距離函數(shù))、云變換(描述兩個屬性之間有關(guān)關(guān)系)。下面先簡介前兩個參數(shù)：概率密度函數(shù)(1)JGW主界面—Modeling—StatModMC,彈出StatModMC主界面，重要填寫內(nèi)容如下：SolidModel：點擊Input，再點擊Solidmodel,選取地質(zhì)框架模型，如solidModel.smgLayers點擊InputSolidModel：點擊Input，再點擊Solidmodel,選取地質(zhì)框架模型，如solidModel.smgLayers點擊Input，再點擊Layers，選取相應(yīng)目的層，如TopPad、Reservoir、BasePad，運用一輸入至右邊界面中，并且修改Spacing值，一般目的層厚度設(shè)為與需要辨別薄層時間厚度值，而非目的層可以設(shè)立得稀某些，如與地震數(shù)據(jù)采樣一致TraceGate:選取反演道門，如wholeProjectQCTraceGate:選取QC道門，如allWellsLine填好后如圖1所示。(2)StatModMC主界面一Mode，共有四個選項：GenerateStats、Simulation>Inversion>Cosimulation，反映其實就是整個流程。在這一步里重要用GenerateStats來分析pdf函數(shù)（即概率密度函數(shù)）。點擊GenerateStats，彈出界面（如圖2所示）。Mrr^pirtaiiaIp*Ul+yirsHartqrwMrr^pirtaiiaIp*Ul+yirsHartqrwWvUi,imimimimimi■■■■ihiiWd ―j-ihfitrJoint...AiiEhMtU-.-

hiuf11*...

CdMUJE..-.Edit..Mils-CHthAHMlhnRv*Hi機Mfh11CMt*l界面上部由空白欄和右邊控鍵構(gòu)成。依照分析需要，選取不同控鍵往空白欄處充填內(nèi)容。在這一步只需用到最上邊一種控鍵Base/Joint，表達分析某一屬性pdf函數(shù)；界面中部選取Outputdirectory填寫輸出途徑；界面下方選取Generateandrun產(chǎn)生腳本并運營；當前開始分層位分巖性分屬性地填寫空白欄內(nèi)容：①TopPad層：點擊Base/Joint，彈出小界面，如圖3所示。點擊Layer相應(yīng)空白欄右邊三角下拉菜單，選取相應(yīng)層位，如TopPad；點擊DiscreteProperty相應(yīng)空白欄右邊三角下拉菜單，選取相應(yīng)離散屬性，普通指巖性屬性（對于非目的層，這一項可以不選）；點擊ContinuousProperty右邊Edit，選取P-Impedance；在Wells一行點擊Edit選取重采樣井途徑，然后在Wellsincurrentdirectory中選取相應(yīng)井，輸入至Selectedwells；點擊OK退出編輯狀態(tài)。這時在GenerateStats主界面中就會浮現(xiàn)相應(yīng)內(nèi)容。②BasePad層：同TopPad層。③Reservoir層：點擊Base/Joint,Layer選Reservoir,DiscreteProperty選定義好巖性曲線，點擊ContinuousProperty右邊Edit,選取P-Impedance。選井。填寫完后GenerateStats界面空白欄顯示如圖4所示。圖4還可以通過選取界面中Copy，選中所有行，最后點擊Paste，選取Pastewells，可以將井復(fù)制到選中位置；最后在GenerateStats界面下半某些Outputdirectory中填寫相應(yīng)途徑，然后點擊Generateandrun產(chǎn)生腳本并運營。(3)查看分析成果。在Linux系統(tǒng)中打開網(wǎng)絡(luò)流覽器Konqueror，并在Location中輸入相應(yīng)腳本地址。在該途徑下除了腳本尚有諸多運營腳本產(chǎn)生重要文獻：?Index.html文獻：里面是分層位分巖性分屬性分析得到所有pdf圖件。以Reservoir層nopay巖性中P-Impedance分析為例，如圖5所示。

￡dd技i r-whQ_^utUO C* %Ldi果L#L￡K』i>m■=NaLaV?nru=HiJfCLi>3i^jiin?Ciwi?-|diriiai'iLihi|iJEJi.Iihil■卜日StatsfarlayerRcscirwoirwherelithotestIsrsapaiyP-Imp&dliarkec'分析井段穿過Reservoir'分析井段穿過Reservoir層nopay巖性樣本點記錄成果區(qū)地質(zhì)本例早cc-lidn-Ddelp?clidM-^del3jtq1rttAreo-ste區(qū)地質(zhì)本例早cc-lidn-Ddelp?clidM-^del3jtq1rttAreo-ste8■1whclePrsjertrrqcrracegat?:■’allNtlliBlil]匕「rmdd■HodeLldepthunit-dutine)準備信息ncde-ljpEa-s*F-s-J,'|■nii:iclayer5pa!：i?gn0.C'04mdd|"I'F-Iriped(jk日],Stype-1P-InF?dar<epr'■niic^jran■c^onenti^iL10.DO]r77"577*^BrplcfleLlo-11wlL/ieE!anFlei,,l■wl111?'vell/iei3aiiFLe7J.wll?r我簽"出 「e^^pRJ’E.vll'rll.''r<Mrple/-3vll1f1w^ell/reSBjple-/4STncde-ll 11imL-riolayeijF^TiD-g■=<'hi"1ncdelL r11|1Lit:fw_tKt11-DL-ccreteFicifertyI早type-1Lithc^te-rt'.ihi3tc;<iran■|'pay,■L1ncip^ay'■]J-r■.■luiG-qrausezn-oiieDti-al?■.".ljlj1輸入地層框架模型輸入分析道門輸入QC道門輸入分析層段以及設(shè)定相應(yīng)采樣間隔BasePad層為非目的層，對其中P-Impedance不用分巖性分析“Reservoir層中nopay巖性中P-Impedancepdf函數(shù)分析成果圖5中顯示了Reservoir層中nopay巖性中P-Impedancepdf函數(shù)分布。其他分析成果與之相似，為了直觀地理解分析成果，顧客應(yīng)當仔細檢查每一層每一巖性中各種屬性pdf函數(shù)分析成果圖件。>*.png文獻：一系列*.png文獻，實際是把Index.html文獻里圖件分開來顯示，一種*.png文獻就是一張圖，這些文獻以便顧客在做報告或文檔時直接插入感興趣圖件，而不用再去抓圖。>scriptCopy.py文獻：就是依照在GenerateStats界面中填寫參數(shù)產(chǎn)生腳本文獻，有點類似程序語言，但是相對簡樸多了。顧客應(yīng)當逐漸熟悉腳本文獻，先嘗試閱讀腳本，領(lǐng)略每某些意思。在逐漸熟悉后可以嘗試修改腳本或自己編寫。足夠純熟后甚至可以跨過填寫界面這一步，直接編寫腳本并運營，來完畢地質(zhì)記錄學(xué)參數(shù)分析、模仿、反演等工作。圖6是該工記錄學(xué)參數(shù)分析腳子。開始逐級分析開始逐級分析r.cdill￡ti!：-：iz'..'l.:.i:l.:i'..'f-l- ；-[jp±J±nc.■Btype-'F-Inpedflij-re^r■.■mic-granezpc-M-ntiaLli.)i]F、1)ZBiplirtfeLla■|'wlL,1jezanFl*-''!.mL1!,1veIl/zesaTpLe/S.arilr1Krll/ilEkaifilEHET.V11\^nll/rEsaifilE^.vlL'j1wllfzEJuplklZE-vll1]]iKid^llF恥XTE二jr'II'Litb:_tent'll'mpsy11I'P-InpedarFre「]?E-sseFrcperty(3typs='P-InFGdar>rfl?Pvaiicgrjn=Dxpco^ntidL|G_lit!]r[丁1oaMuplWgLlz=IH.sic^anpL<;L.hLI",'-Jil'.■i>；2aipLfl.il3-.M,llprggJ'irsll/rfljajpla/3ST.vllprrma11/rasaifilB{4.?1L11wll^Ea^aiplD^SI_vll1]|

Reservoir層為目的層，對其中P-payImpedance應(yīng)當分巖性進行pdf函數(shù)分析nopay圖6>stats.py文獻：以上分析成果詳細數(shù)據(jù)都存儲在該文獻內(nèi)，可以說是最后成果文獻，后來在做模仿與反演過程中調(diào)用就是該文獻。以Reservoir層nopay巖性為例顯示該文獻某些內(nèi)容，如圖7所示。#5n:itc-f5 e壯Miwdpre-p-ertic-5-L*P-Imp^danM11]s#5n:itc-f5 e壯Miwdpre-p-ertic-5-L*P-Imp^danM11]sfMst-D.jra.Tis?drrespontfvicftrjhi^tograizs-[i1qHintngr^m|quantiles-[(O.OOIh弟DBQ50r4W叩(0.12^也馳”DLQBZ50r4^73*345jj?口.3"『"典GH(D.SDDji5142502：-^對43戲口：、[CCi.?50,3駐“站L^D.B75P由荷S-細％b474097^(QB999r5661155).DietdiMe*齡枷商￡redtlannjptwoWhereIl-sye-r1F--5trrcir1rlithG-lcgy~1litho_t?si：1r1abel "CoEE^liAtiariPdt(Reservoir層nopay巖性P-Impedancepdf函數(shù)分析成果，以函數(shù)（Q,X）表達，含義為P-Impedance不大于X概率為Q?變差函數(shù)求取(1)JGW主界面一Modeling一Specials一StatMod一StatModAnalysis，彈出StatModAnalysis界面;(2)Input一Time/Depthmode：TimeSolidmodelSOLIDLayers :點擊Selectlayers選取目的層位，如Reservoir層，在Selectedinterval設(shè)立采樣間隔為0.001s，點擊Assigncurrentintervaltoselectedlayers應(yīng)用修改采樣間隔Lithologymasks:先把Enablelithologymasking左邊小方框點亮，再點擊Lithology右邊List，在彈出界面中選取事先劃分好巖性曲線，如Lithology，這時在Selectlithotypestouse中選取某一巖性，如pay(其他巖性稍后也需要分析，只是一次只分析一種巖性中屬性變差)Dataforhistogramsandtransforms:在彈出界面上半部Data下方Wellfile(s)右邊List中在../WELLS/RESAMPLED選取相應(yīng)分析井曲線，如P-Impedance。Tracegate:wholeProjectEdit一Transforms，彈出兩個界面，一種是主界面StatModTransform，一種是輔界面Transform，點擊輔界面Transform最上方Histogram，彈出一種小界面StatModHistogram并帶著一種不斷閃爍紅色提示Buildthehistogram，點擊該紅色提示，則主界面浮現(xiàn)相應(yīng)樣本點分析，調(diào)試小界面中Nrofintervals值，盡量使pdf函數(shù)樣本點分布規(guī)則持續(xù)。調(diào)試好后點擊Dismiss退出。點擊輔界面Transform上方第二個選項Transform，彈出界面如圖8,同步在主界面中浮現(xiàn)一條藍色擬合線，如圖9。用鼠標左鍵調(diào)節(jié)藍色擬合線，使得該曲線形態(tài)能大體反映樣本點分布形態(tài)，最后用ctrl+鼠標左鍵停止編輯藍色擬合線。然后在圖8顯示界面下方選取Saveasfirst，在彈出界面中給出名字，如pay_p-imp_histogram，再持續(xù)點三次OK回到StatModTransform主界面，點擊Dismiss退出界面。-1 jik*wiMuturetKmmcwdlfiid心"? UM…5Wn alVdurn-MitUi-1 jik*wiMuturetKmmcwdlfiid心"? UM…5Wn alVdurn-MitUir±L*[pm囑h >■■!.?AmH■CM'Wlxwina+i?iuHakmuiUEfilUtCU-jCaUs,輔界面主界面圖8 圖9⑶Input一Dataforvariogramsamplingandmodeling一Primarydata一Welllogdata,在彈出界面中上部Data某些Wellfile（s）右邊List中選取相應(yīng)分析井曲線，如P-Impedance，然后點OK,這時彈出兩個界面，一種StatModVariogram主界面，一種Variogram輔界面。點擊Variogram輔界面上方第一種選項SampleVario,彈出界面StatModSampleVariogramParameters，如圖10所示，先點擊在該界面下方不斷閃爍紅色提示BuildtheSampleVariogram，則主界面中浮現(xiàn)變差函數(shù)樣本點分析曲線，修改圖10界面中SelectSamplingParameters某些Interval[m]下方空白處值（即變程單位距離），則主界面中變差函數(shù)樣本點分析曲線形態(tài)發(fā)生變化，不斷調(diào)試Interval[m]值，使得主界面中變差函數(shù)樣本點分析曲線形態(tài)盡量光滑，如圖11所示，重要看上半部紅色樣本曲線，最后點OK結(jié)束編輯。FLritliia…]?'EIKri-HC.*r輔界面FLritliia…]?'EIKri-HC.*r輔界面圖10圖11（闡明：圖11主界面顯示是進行變差函數(shù)分析樣本點曲線，變差函數(shù)是三維，有縱向變差，也有橫向上變差,圖11輔界面下方表達是將三個方向上變差函數(shù)曲線用三種顏色顯示。但是由于這里是運用井上樣本點進行變差函數(shù)分析，一方面井曲線在縱向上樣本點個數(shù)普通能滿足記錄需要（普通50個樣本點就可以得到比較抱負pdf

函數(shù)和變差函數(shù)），另一方面井曲線在橫向上分布是遠不能達到分析需要。因此在這里只相信縱向上變差函數(shù)分析成果，即只需分析紅色曲線）點擊Variogram輔界面選項ModelVario，彈出界面如圖12所示，同步把鼠標箭頭放在圖13上部，按住左鍵可以得到一條擬合曲線，可以任意拖動該擬合曲線，使得它能最佳地擬合縱向變差樣本點，特別是前幾種樣本點（思考為什么？），ctrl+鼠標左鍵停止編輯。.ii MlMdliYrLafl--*Rir>:dMAiww:HfimriiiikalirmftvIiKnS##*uc-M-n-d.....hcrtailwf…talfcL?_rr皿.ii MlMdliYrLafl--*Rir>:dMAiww:HfimriiiikalirmftvIiKnS##*uc-M-n-d.....hcrtailwf…talfcL?_rr皿'Lfrh曰心|pribV爾H蟬微說liHAfig前的帕做- ?*Ifd-kriUm編間1喩皿1>-1-齢Xct:utl?T5?]att■HEt?dFirmSuii>ri?3iifc■cLIdfmHMEk*v?itcaLirllsrija.SB.如ATid-PltZCMMk.圖12 圖13這時圖12顯示界面上部空白欄處會浮現(xiàn)變差函數(shù)表達式（Type，Sill，X，Y，Z，Power）相應(yīng)詳細參數(shù)。在本例中表達變差函數(shù)類型是指數(shù)型，基臺值（Sill，即擬合曲線拐點值）為3.07，三個方向（x,y,z）變程均為26.9m（固然就像之前所說只相信縱向上變程分析成果）。此外，還可以變化變差函數(shù)類型，在圖12界面中上某些，F(xiàn)uncType下方，點擊小控鍵，彈出幾種類型，普通用到為Gaussian、Expo兩種類型。在原點附近，Gaussian型變化較慢，Expo型變化較快，顧客可以分別嘗試。如果不滿意分析成果，可以在圖12中選中上方空白欄里函數(shù)，然后點擊下面Remove，就可以刪掉該函數(shù)，然后重新編輯。最后記住測試好變差函數(shù)類型與縱向上變程（單位為深度域單位m）依次對Reservoir層其他巖性類型中屬性(如P-Impedance)進行變差函數(shù)分析，并記住測試好變差函數(shù)類型與縱向上變程。再對其他目層進行分巖性屬性變差函數(shù)測試。過程同上。在分層擬定了不同巖性變差函數(shù)類型與縱向變程后，橫向變程怎么擬定呢?這里提供兩中思路：一是與地質(zhì)概念相結(jié)合，與地質(zhì)專家或客戶進行交流，大體掌握工區(qū)沉積模式與巖性展布，以此來擬定一種大概橫向變差范疇(只規(guī)定一種大概輸入)；一是根據(jù)疊后CSSI反演成果進行巖性地質(zhì)規(guī)模和分布等地分析，來得到一種大概橫向變差范疇。不論是橫向還是縱向上變差，通過度析得到都是一種大概值，對反演成果起到是軟約束作用，而起到硬約束影響作用是地震。第四步.地質(zhì)記錄學(xué)模仿 目是測試第三步中擬定地質(zhì)記錄學(xué)參數(shù)與否合理(1)StatModMC主界面一Mode一Simulation，彈出Simulation主界面，該界面上部和GenerateStats界面基本一致。以模仿目的層Reservoir中不同巖性P-Impedance為例來簡介需要填寫內(nèi)容。點擊Base/Joint，在彈出小界面中Layer選Reservoir層，DiscreteProperty選巖性曲線，Continuous中選P-Impedance，OK退出。這時Simulation主界面中浮現(xiàn)有關(guān)信息。如圖14所示。對界面中內(nèi)容進行編輯：涉及巖性比例、巖性變差函數(shù)、不同巖性中P-Impedancepdf函數(shù)與變差函數(shù)，也就是前面幾步分析得到成果在這里都要用到。選井。

圖14有關(guān)參數(shù)：>巖性比例：選中Reservoir層第一行，如圖14所示，然后點擊Edit，彈出界面如圖15所示。在Proportions某些編輯巖性比例（闡明：巖性比例缺省值為1.0,本練習工區(qū)巖性比例設(shè)立可以參照pay:nopay=13:87，可以從井上目的層段不同巖性樣本點記錄得到一種粗略預(yù)計，但仍需測試）>巖性變差函數(shù)：在圖15中Variograms某些，點擊Single右邊Edit，在彈出菜單中依次將Type、Verticalrange>Lateralrange設(shè)定為gaussian、0.010、600x600。（闡明：Type、Verticalrange均為在第三步中測試到值，Lateralrange可以由第一步得到，即運用疊后反演成果進行砂體雕刻，以求對目的區(qū)巖性展布、比例有一種總體上對的把握，同步結(jié)合工區(qū)地質(zhì)沉積模式等地質(zhì)方面信息來指引Lateralrange輸入。顧客可以放心是Lateralrange輸入可以只是一種大概值，它對真實各種巖性體分布起到是軟約束作用，在背面反演某些要講到U，對巖性體分布起硬約束作用是地震數(shù)據(jù)。并且Type、Verticalrange>Lateralrange設(shè)定也不是一蹴而就，需要調(diào)試，尋找合理值范疇）iw Mi.rhiw Mi.rh> -忡土…Ml"呷imMLt..編輯屬性變差函數(shù)ETIiUltkbiT [-|E^e…MiLLs1hcwrrartdirK-forpMl：!圖15圖16圖15>不同巖性中P-Impedancepdf函數(shù)與變差函數(shù)：在Simulation主界面中選中Reservoir目標層某一巖性相應(yīng)行，如選取Lithology.nopay所在行，然后點擊右邊控鍵欄Property下方Edit，彈出界面如圖16所示，只需對紅框圈出某些進行編輯。點擊Histogram相應(yīng)Edit，選取相應(yīng)存儲pdf函數(shù)分析成果文獻。人們還記得這個文獻在哪嗎？就在第三步做地質(zhì)記錄學(xué)參數(shù)分析時，在GenerateStats界面運營腳本產(chǎn)生stats.py文獻，里面存儲所有目的層中不同巖性中P-Impedance等屬性pdf函數(shù)分析成果。圖7展示stats.py文獻某些內(nèi)容。人們需要好好回憶一下，整頓好頭緒。然后點擊Variogram相應(yīng)Edit，在彈出小界面中依此設(shè)定Type、Verticalrange>Lateralrange：exponential、0.006、400x400。這個成果也是由第三步分析得到。持續(xù)兩次點擊OK退出編輯狀態(tài)。編輯其他巖性中屬性pdf函數(shù)與變差函數(shù)（其實pdf函數(shù)調(diào)用都是stats.py文獻，只需編輯變差函數(shù)即可）。變差函數(shù)值可參照如下表1（例子）：?--VariogramLithologyTypeVerticalrangeLateralrangepayexponential0.0161200x1200nopayexponential0.0161200x1200表1在編輯參數(shù)時可以使用主界面中Copy、Paste功能，實現(xiàn)Histogram與Variogram迅速編輯。界面中下方RunParameters某些也有某些選項，在Nrofrealizations可以填寫想要產(chǎn)生實現(xiàn)個數(shù)，Randomseed中值隨之變化。如Nrofrealizations填1，表達只產(chǎn)生一種實現(xiàn)，這時Randomseed自動將種子點設(shè)為0；如Nrofrealizations填2,表達產(chǎn)生兩個實現(xiàn)，Randomseed自動將種子點設(shè)為0,1。在實際工作中顧客也許需要產(chǎn)生更多實現(xiàn)。顧客可以自己熟悉其他選項。最后別忘了給出輸出途徑與輸出文獻名。點擊Generatescript產(chǎn)生腳本文獻，可以從網(wǎng)絡(luò)流覽器中調(diào)出腳本文獻查看，如圖17展示腳本中某些內(nèi)容。最后點擊Run運營該腳本文獻。（留住Simulation主界面）rc-d.elHcde1?!■adLlnx-de1■!a<-lidHcdel>.芬口qr』■traceq'a.te-'MticlePrcject■qcTiaceijatr-'^llHeLLcLLnEr)ntde1.EtsGrTC-le?.nicrclaycr^pacLng .J'Jlnc-dcJ.E?￡￡-r7,c-lr|'1.11hc-_r??5tr|■DteergtePrepsrt},(type-■r.1.11hc-_tcst-.hLStc-qran-| rpa^'■15.0!■,■-?aiLeirans-o.oi*nuggett?gussLantChdLL12ULUQS0.QL^EKp&nEntlairo.OLli,L2nnHlZO0hKtLLs-\ L.wll\'ijfHjJf'[tsanpLe-/3..wll\'WtLL/iesanple/SST.ijsJLB|1nc-del.Reaerrc-Lc|rlLl:i__LeL：lrr1|'7Z:v.e' -BaEeProp*erLyi*早.type=rP-Irpedanc±*i ■:■:.caT2.r-J.JZ ■:cr.:l1 .1‘.JJ,：[二二丨】 nOpAVvellc-|rwEll.rc3-ar7ple.!1Luwll'we11rc^arpleJJ11'weLL/rc^anpl iReservoir層中每一巖性中P-Impedanceredel.Ejcscrrc-LE|'1ithc_Lcstf||Fpij-']|1P-Inp^danccF|-Ba^Epre^rty；S : 變差函數(shù)b^jrLc.qra.n-U.□]anuggst* !gKpcncntlaL(U.11PL2Udr]jiaa|, \?H..尸-mi- c■'1.Al1'r'-J-：]■■rt^^rp.- .w..■, pay圖17(2)查當作果并質(zhì)控：①打開網(wǎng)絡(luò)流覽器，查看運營腳本產(chǎn)生Index.html文獻，選取該文獻，然后在Runreport下面選取 ，里面是所有反映模仿信息圖件和輸入輸出參數(shù)登記表，圖18展示其中某些內(nèi)容。平均值原則偏差nopay中輸出P-nopay中P-Impedance輸入平均值原則偏差輸出pdf函數(shù)曲線Impedancepdf函數(shù)重要參pay中P-Impedance輸pay中P-Impedance輸■尹入輸出pdf函數(shù)曲線圖18在該圖件文獻里面重要檢查兩方面內(nèi)容：目的層中每一種巖性內(nèi)相應(yīng)屬性，如P-Impedance輸入輸出pdf函數(shù)形態(tài)應(yīng)當大體一致，重要比較入輸出pdf函數(shù)平均值與原則偏差（圖18中藍色曲線表達輸入，紅色曲線表達輸出）；目的層輸入輸出巖性比例應(yīng)當大體一致（偏差在5%左右）。②JGW主界面一Analysis一SectionView，先選好tracegate（最佳選取在StatModMC主界面中填入QCTraceGate：本例中是MyQCGate），然后從剛才填寫模仿成果輸出途徑中調(diào)出成果查看。打開這個成果方式重要有兩種：SectionView界面一File—Opensession一進入模仿成果輸出途徑一選取QC.jvis,這時顯示模仿成果，涉及巖性體與P-Impedance體；SectionView界面一Input一Stratigraphydata一進入模仿成果輸出途徑一選取lithology.hor/P-Impedance查看。質(zhì)控辦法與原則：和疊后反演成果進行比較，如圖19所示。與疊后CSSI反演成果（右圖），在反映砂體展布、地質(zhì)體規(guī)模、巖性比例等方面大體比較一致，這時模仿成果是比較合理。記住相應(yīng)模仿參數(shù)。圖19③比較輸入輸出井曲線。在運營腳本過程中，自動在相應(yīng)目錄下產(chǎn)生在井點處抽取“偽井曲線”（Lithology與P-Impedance曲線）。檢查這些偽井曲線，Jason主界面-Analysis WellManage中選取相應(yīng)途徑下井，會發(fā)現(xiàn)相應(yīng)Lithology與P-Impedance,分別有Lithology_in，Lithology_out與P-Impedance_in,P-Impedance_out。*_in表達輸入重采樣曲線，而*_out表達輸出從模仿成果中井點處抽取出來巖性或?qū)傩郧€，如圖20所示，藍色表達輸入，紅色表達輸出，兩者偏差很小，這是比較抱負模仿成果。圖20④建議多測試幾種參數(shù)，如變化目的層巖性比例、變差函數(shù)類型和縱橫向變程范疇來產(chǎn)生實現(xiàn)，并比較差別，特別注旨在SectionView/MapView中與疊后反演成果進行比較，選取令實現(xiàn)與疊后反演成果最相似參數(shù)組合；在比較剖面同步注意檢查生成Index.html文獻，比較原則與上面所說一致。建議做一種列表來紀錄實驗參數(shù):ProportionValueDPVarioValueP-impedanceValueVarioOutputdirectorywell-basedGaussian-shortExponential-longwell-basedGaussian-longExponential-longLowerPayGaussian-shortExponential-longLowerPayGaussian-longExponential-longHigerPayGaussian-shortExponential-longHigerPayGaussian-longExponential-longwell-basedGaussianlongPay：Exp-longNonpay：Exp-short第五步.地質(zhì)記錄學(xué)反演? 無井約束反演——重要目是測試地震權(quán)重(1)StatModMC主界面一Mode一Inversion，彈出Inversion主界面。顧客這時也許已經(jīng)發(fā)現(xiàn)Inversion主界面與GenerateStats界面、Simulation主界面基本一致。區(qū)別大地方重要在界面左下方SeismicBelief某些。在該某些需要填寫地震數(shù)據(jù)權(quán)重。此外，由于是反演，因此需要在目的層上下加個“蓋子”，即加入TopPad層和BasePad層作為邊界層(邊界層厚度要在半個子波長度以上，以消除反演時子波旁瓣效應(yīng))。這一步得回到StatModMC主界面選取，把TopPad層和BasePad層Spacing設(shè)為0.004s。打開Inversion主界面，對于目的層參數(shù)填寫與在Simulation主界面填寫一致，固然也可以用測試效果更好參數(shù)代替。對于非目的層，不用分巖性反演屬性。不要選井。(2)當前來考慮地震權(quán)重。在Inversion主界面左下方SeismicBelief某些，有Seismic和Wavelet兩個控鍵，點擊Seismic，則彈出選項如圖21所示。Seisidcbeli&f $點5?2kavolvtS觸址kNvissdB▼DafMilth&lsaI□Coirtiiit圖21在StatModMC反演模塊里，地震信息是通過設(shè)定一定噪音引入，或者更好理解一些，可以通過設(shè)定信噪比引入。信噪比計算公式如下所示：S/NRatio血dB)=20*lg(std(seismic)/std(residual))其中，lg表達10為底對數(shù)；std表達原則偏差；在該公式中，可以把殘差residual理解為噪音(即地震數(shù)據(jù)與合成記錄差值，或者說是合成記錄與地震數(shù)據(jù)匹配不上某些)，地震數(shù)據(jù)是一定，信噪比越高，表達殘差越小(既噪音越弱，地震數(shù)據(jù)與合成記錄越相似)；反之亦然。因此可以通過控制信噪比大小來控制引入噪音。點擊Seismic相應(yīng)Edit，選取全疊加地震數(shù)據(jù)，然后在Noise一欄，一方面點擊Noise右邊下拉菜單，可以看見三個選項value/dB/file，選取dB，表達引入信噪比來控制合成記錄與地震有關(guān)性。然后填入詳細信噪比值。那么如何設(shè)立這個值呢？可以參照疊后反演時生成一種文獻(當前體會到在做地質(zhì)記錄學(xué)反演前需要做疊后CSSI反演重要性了吧？)。一方面從JGW主界面下Analysis中打開MapView,然后選取Input一Horizon,找到途徑../INVER_TRACE_PLUS,在這一目錄中有一質(zhì)控層位文獻qc.hor,選取該文獻，然后在彈出各種層位中選取invertedsignaltonoiseratio（Signal-to-noiseratio）,然后OK,這時MapView中顯示目的區(qū)疊后反演信噪比成果，可以打開色標對照看，如圖22所示。圖22由圖22可看出本工區(qū)信噪比范疇大概在5?25dB之內(nèi)，且工區(qū)中大某些信噪比在10dB以上。因此咱們可以依照疊后反演信噪比來設(shè)立地質(zhì)記錄學(xué)反演信噪比，它不是某一確定值，而是要在某一數(shù)值范疇內(nèi)進行測試。對于本工區(qū)建議測試5dB、10dB、15dB、20dB幾組值。也就是要產(chǎn)生四個腳本并運營。固然，如果顧客感興趣可以測試更多值，并感受隨之帶來變化。填入地震信息后，一定別忘了還需填入子波信息。點擊圖24中Wavelet選項，在彈出界面中缺省是用單一子波反演，即Single,點擊Single相應(yīng)Edit,選取相應(yīng)子波（在疊后CSSI反演時用最后子波）。在RunParameters某些，顧客可以自己設(shè)立需要計算實現(xiàn)，由于這步重要測試地震權(quán)重，因此不建議產(chǎn)生太多實現(xiàn)（要花時間）。最后給出輸出途徑與腳本名，再產(chǎn)生腳本并且運營即可。（3）檢查并質(zhì)控成果。?檢查反演成果并進行抽井檢查：從JGW主界面下Analysis中打開SectionView,先在InputTracegate中選取適當?shù)篱T，最佳是聯(lián)井線，這里選QCTraceGateo然后選反演數(shù)據(jù)，點擊Input,選取Stratigraphicdata,或者是在界面上找到圖標吳并點擊，即可查找反演體文獻*.hor（到剛才輸出途徑中查找）并顯示。一種實現(xiàn)相應(yīng)兩個反演成果：P-

Impedance.hor和lithology.hor。調(diào)出剖面后會發(fā)現(xiàn)地質(zhì)記錄學(xué)反演成果具備很高辨別率，但是這一成果與否就是顧客想要對的成果呢？固然要進行質(zhì)控。這一步最重要是把井選進來，進行抽井檢測（這已經(jīng)是Jason在顧客中普及得最深刻質(zhì)控手段了，不論是疊前還是疊后反演，要判斷反演成果好壞最基本就是進行抽井檢查，地質(zhì)記錄學(xué)反演也不例外）。例如你比較反演巖性體與井上巖性曲線，那么一方面你要看到在井周邊反演巖性體與井與否是比較吻合，然后是觀測井間巖性與否有比較好連通性，好反演成果應(yīng)當可以預(yù)測大某些規(guī)模較大地質(zhì)體并且對于大多數(shù)細節(jié)某些（如薄層）也能有較好反映。如果在這一步得不到好成果，那么就算把井加上作為約束，成果也還是不可靠。如果成果不滿意怎么辦？需要繼續(xù)測試輸入?yún)?shù)，如記錄學(xué)參數(shù)，甚至井震標定，直到得到滿意成果為止。如圖23為抽井檢查反演巖性體實例。圖23＞輸入輸出pdf函數(shù)、巖性比例質(zhì)控：不論在記錄參數(shù)測試環(huán)節(jié)、模仿環(huán)節(jié)還是在反演步驟，都能在相應(yīng)途徑找到index.html文獻。打開網(wǎng)絡(luò)流覽器，在相應(yīng)途徑下選取該文獻，然后在Runreport下面選取QChistograms（detailed），能看到反映模仿信息圖件和輸入輸出參數(shù)登記表，如圖24顯示是其中某些內(nèi)容——reservoir層中nopay巖性P-Impedancepdf函數(shù)、巖性比例輸入輸出狀況。LayerReservoirwherelhho_testisnopayfiMeFwair-PUmpedanicenuirPoints1W71S阪m5.07^4-06StDevcoeFVw3.歐）'-Ii5Be+06MasS取76Prcportian0.93Ri-iwrvcjr.P-imiwdLiirw*nop畔inputMranSIDev'2.97e^-D5coefV#@0579圖24圖中藍色線與紅色線分別為輸入輸出reservoir層中nopayP-Impedancepdf函數(shù)，兩者相似度較高，是比較抱負輸入輸出狀況。圖中右邊是描述輸入輸出pdf函數(shù)參數(shù)表，其中決定pdf函數(shù)形態(tài)參數(shù)是平均值Mean、原則偏差StDev，顧客應(yīng)當比較輸入輸出平均值與原則偏差；此外還應(yīng)比較輸入輸出巖性比例，即圖表中Proportion，兩個比較原則都是輸入輸出應(yīng)當大體一致。這某些其實與模仿環(huán)節(jié)比較內(nèi)容一致。＞反演成果與疊后反演成果比較：在比較過程中重要注意兩方面內(nèi)容一是巖性比例二是巖性形態(tài)與連通性。由于使用相似正演褶積算子來模仿地震數(shù)據(jù)，因此地質(zhì)記錄學(xué)反演成果應(yīng)當與疊后反演成果比較一致（這個“一致”重要指大套巖性規(guī)模分布等，而不強調(diào)細節(jié)，如薄互層），如圖25所示，左邊是地質(zhì)記錄學(xué)反演P-Impedance成果，右邊是疊后反演成果，兩者對大套巖性體規(guī)模尺度及分布反映都比較一致，區(qū)別僅是地質(zhì)記錄學(xué)反演成果能反映出更豐富細節(jié)。圖25＞巖性體沿層切片檢查：比較地質(zhì)記錄學(xué)反演巖性體沿層切片與地震數(shù)據(jù)沿層提取均方根振幅屬性。地質(zhì)記錄學(xué)反演巖性體沿層切片上地質(zhì)體形態(tài)應(yīng)當與地震比較一致。如果不一致，得考慮如下兩個問題：一是設(shè)立信噪比過低了；二是選取變差函數(shù)與地震信息不符，如變程不不大于重要地質(zhì)體尺度與規(guī)模。打開地質(zhì)記錄學(xué)反演巖性體沿層切片方式為：打開MapView一Input一Stratalslicer，在彈出界面中，相應(yīng)StratigraphicData選取相應(yīng)反演巖性體，在Layer中選取相應(yīng)層，在Microlayer一行用先后翻鍵查看微層中巖性分布狀況。?無井約束反演——重要目是學(xué)會產(chǎn)生各種實現(xiàn)并進行記錄分析運用界面或腳本產(chǎn)生各種實現(xiàn)：如果對腳本不熟悉，那么運用界面產(chǎn)生各種實現(xiàn)也是非常簡樸，做法其實在前面也提到過。在上一步反演界面中擬定最佳地質(zhì)記錄學(xué)參數(shù)與信噪比并填寫好，然后在界面右下方RunParameters某些第一種選項Nrofrealizations填寫想要產(chǎn)生實現(xiàn)個數(shù)，如10(上一步中建議為了節(jié)約時間，只產(chǎn)生一種實現(xiàn))，然后把鼠標放在第二個選項Randomseeds相應(yīng)選項內(nèi)單擊，就會浮現(xiàn)0，1，2，3 9,這時準備工作就做好了，再點擊Generatescript產(chǎn)生腳本，感興趣顧客可以打開該腳本看看腳本中用來產(chǎn)生各種實現(xiàn)語句是什么？最后點擊RUN運營反演。這個過程也許會花較長時間，因此白天做測試，下午離開辦公室之前運營腳本是最佳選取。運營完畢后，在相應(yīng)目錄下產(chǎn)生10個途徑：seed_0、seed_l、seed_2 、seed_9，每個途徑下都是一組完整反演成果，涉及P-Impedance.hor和lithology.hor，以及index.html文獻等，顧客應(yīng)當按照上一步講述質(zhì)控環(huán)節(jié)檢查成果，并注意每個實現(xiàn)之間差別。對各種實現(xiàn)進行記錄分析：當前產(chǎn)生了各種實現(xiàn)，并且這些實現(xiàn)都是等概率，那么哪個實現(xiàn)更合理呢？StatModMC提供了對各種實現(xiàn)進行記錄性分析辦法。在反演界面中點擊OK退出，StatModMC主界面一Utilities一Summarystatistics，彈出界面，在Directory相應(yīng)Edit中選取相應(yīng)目錄，這時在中部左邊方框中浮現(xiàn)seed_0、seed_l、seed_2、 、seed_9十個實現(xiàn)，全選然后用一輸入到右邊方框中，最后在界面下方Outputdirectory中填寫輸出途徑(缺省為summary)，點擊RUN運營，運營速度較快。運營結(jié)束后在相應(yīng)途徑下查當作果：對于巖性體來說，分析成果涉及每一種巖性概率體(例如四種巖性就有四種巖性概率體)以及一種極大似然巖性概率體；對于P-Impedance體來說，分析成果重要有四個分別是平均值P-Impedance_mean.hor、最小值P-Impedance_min.hor、最大值P-Impedance_max.hor、原則偏差體P-Impedance_stddev.hor。顧客可以從相應(yīng)途徑下打開這些成果來查看。檢查并質(zhì)控成果：＞比較各個實現(xiàn)與疊后反演成果：重要比較這些實現(xiàn)(P-Impedance.hor和lithology.hor)與疊后CSSI反演成果對大套巖性體規(guī)模尺度反映上與否基本一致，而不用太摳細節(jié)。＞對巖性概率體進行抽井檢查：在SectionView中打開巖性概率體，如此說來圖26所示，并投上井上巖性曲線進行抽井檢查。左圖為泥巖概率體，右圖為氣砂概率體。紅色代表概率為1,即10個實現(xiàn)中普通以上都反映有該種巖性；黑色表達概率為0,即所有實現(xiàn)都沒有這種巖性。比較時注意兩點：一是在井點處巖性概率體與否與井上巖性相符，如果不符，那么不擬定性有多大？二是井間或遠離井點處巖性不擬定性有多大？如果在井點處吻合得較好，而井間不具備連通性，即巖性概率減少，那么要考慮這種變化與否與實際狀況相符，如果不相符，那么要考慮加入信噪比與否適當？選取巖性橫向變程與否過?。繄D26＞平均P-Impedance體與疊后反演成果比較：準則是兩者對大套巖性規(guī)模尺度反映大體一致。＞查看極大似然巖性概率體：調(diào)出極大巖性概率體并且投上井巖性曲線進行抽井檢查，在井周比較極大似然巖性概率體與井上巖性曲線吻合限度，在井間檢查巖性連通性等等。＞在井點處比較井與從平均P-Impedance體中抽出P-Impedance曲線：從平均P-Impedance體中抽出P-Impedance曲線要可以反映井上P-Impedance大體趨勢。＞輸入輸出pdf函數(shù)、巖性比例質(zhì)控：打開網(wǎng)絡(luò)流覽器，在各個實現(xiàn)對應(yīng)途徑下選取index.html文獻文獻，然后在Runreport下面選取 ,查看每層輸入輸出巖性比例、每層每種巖性中屬性輸入輸出直方圖，并且比較至少三個實現(xiàn)這些信息，注意每一實現(xiàn)以上信息之間差別不應(yīng)當太大。無井約束反演——重要目是把握地質(zhì)記錄學(xué)反演不擬定性地質(zhì)記錄學(xué)反演具備不擬定性和顧客對工區(qū)地質(zhì)、地震結(jié)識關(guān)于，這些結(jié)識不擬定性導(dǎo)致地質(zhì)記錄學(xué)反演輸入不是絕對精確，在StatModMC中這些對反演成果產(chǎn)生重要影響輸入為：.巖性比例.巖性變差函數(shù)類型與變程范疇/屬性變差函數(shù)類型與變程范疇圖圖27■地震噪音（輸入信噪比）顧客對以上輸入應(yīng)當重復(fù)測試，辦法是先找到大體合理參數(shù)值范疇，然后在該值附近范疇對以上參數(shù)可以進行組合測試，并且相應(yīng)每個參數(shù)組合產(chǎn)生3個實現(xiàn)。顧客可以依照實際工區(qū)進行測試參數(shù)范疇設(shè)立。最后產(chǎn)生腳本并且運營，對反演各種實現(xiàn)進行上一步列出一系列質(zhì)控工作。?加井約束反演 在擬定最佳輸入后最后反演成果在通過上一步進一步測試后擬定最抱負輸入?yún)?shù)（可以是一組也可以是幾組），然后把井加進來進行約束，在界面上加井操作如前面幾步所述，然后產(chǎn)生腳本并運營，最后對最后反演成果進行質(zhì)控。質(zhì)控如前所述，重要涉及：＞輸入輸出pdf函數(shù)、巖性比例質(zhì)控：打開網(wǎng)絡(luò)流覽器，在各個實現(xiàn)相應(yīng)途徑下選取index.html文獻文獻，詳細檢查比較內(nèi)容如前所述。＞地質(zhì)記錄學(xué)反演成果與疊后CSSI反演成果比較：比較原則如前所述。＞巖性體/屬性體沿層切片檢查：檢查辦法如前所述。＞查看巖性概率體：檢查辦法如前所述。＞平均P-Impedance體與疊后反演成果比較：比較辦法如前所述。＞P-Impedance原則偏差體：在井附近應(yīng)當能看到明顯“牛眼”，如圖27所示，表達井周邊P-Impedance原則偏差明顯減?。芟朊靼讍?？由于加入了井做約束，井的確加入了地震之外信息）。反而如果沒有“牛眼”效應(yīng)，也許暗示存在問題有：一是輸入信噪比過低；二是反演成果并沒有受井約束。＞檢查極大似然巖性概率體：檢查過程如上一步說述。 ?第六步.地質(zhì)記錄學(xué)協(xié)模仿 在講操作之間，先來考慮兩個問題： （1）在地質(zhì)記錄學(xué)協(xié)模仿時引入了云變換這個概念，那么為什么要用云變換來描述P-Impedance與孔隙度（飽和度 .....）關(guān)系？這一步是在地質(zhì)記錄學(xué)反演之后進行。在得到多組等概率合 理巖性體與P-Impedance體實現(xiàn)基本上，如果要進一步得到儲層屬性體（如孔隙度體、飽和 度體等）,可以通過協(xié)模仿得到。運用已得到巖性體與P-Impedance體，并結(jié)合地質(zhì)記錄學(xué) 參數(shù)中P-Impedance與孔隙度（飽和度 ）兩兩之間關(guān)系（在StatModMC中專業(yè)述語叫云 變換），在兩者基本上來計算儲層屬性體?？倎碚f，運用云變換計算孔隙度體（飽和度體 ），比通過單純用線性關(guān)系擬合P-Impedance與孔隙度（飽和度 ）關(guān)系來計算孔隙度體（飽和度體 .....）要精準諸多，為什么？由于真實地質(zhì)環(huán)境中P-Impedance與孔隙度 （飽和度…）關(guān)系普通都是非線性，而是如圖28所示像云團同樣分布關(guān)系，而用線性關(guān)系 對這種分布關(guān)系進行擬合，只是一種近似關(guān)系（這個大多數(shù)顧客都應(yīng)當很清晰），例如在圖 中當P-Impedance值為5e+06時，相應(yīng)孔隙度值有紅色雙箭頭線涉及范疇那么多，但是用線 性或曲線關(guān)系擬合時，普通取是紅色雙箭頭線中值，其他也許值就不考慮了。而協(xié)模仿會將 紅色雙箭頭線涉及整個孔隙度值范疇都考慮進來，將這些也許值作為一種概率分布來考慮，這種帶著地質(zhì)記錄學(xué)觀點思路顯然更科學(xué)，固然也更復(fù)雜。StatModMC就為顧客提供這種辦法（并且實現(xiàn)起來很輕松），從而使得更客觀地結(jié)識真實儲層屬性成為也許。圖28 （2）為什么不直接反演儲層屬性體？如果要得到儲層參數(shù)體，例如說孔隙度體，需要分析井段中有孔隙度曲線。那么有顧客也許會問：為什么不象反演P-Impedance體那樣，直接運用第三

步中分析得到孔隙度概率密度函數(shù)與變差函數(shù)并結(jié)合地震與井信息來反演孔隙度體。這個觀點很容易從對P-Impedance體反演中引用過來啊。為什么還要加入云變化作為求取孔隙度體橋梁呢？對此，人們可以考慮一種問題，如圖29所示，相應(yīng)含油井段，P-Impedance值明顯減少，而孔隙度明顯增大，相應(yīng)這種關(guān)系如果單純從兩個屬性概率密度函數(shù)上來看是看不出這種反映好儲層屬性之間關(guān)系，而通過云變換，卻能較好反映出當P-Impedance值減小時孔隙度值增大關(guān)系，從而提示這段也許是儲層，如圖30所示。圖29圖30圖29（3）操作:>StatModMC主界面一Mode一GenerateStats，

彈出界面，該過程與第三步中環(huán)節(jié)類似，只是在選取Reservoir層時，由于在記錄P-Impedance同步還要記錄其與孔隙度（飽和度 ）云變換關(guān)系，因此選取了Base/Joint選項后在彈出界面中Continuousproperties中除了P-Impedance外還要多選，如孔隙度（飽和度 .....）。填好后界面如圖31所示。圖31產(chǎn)生腳本并運營。并且在相應(yīng)途徑檢查運營成果:①Index.html文獻：里面是分層位分巖性分屬性分析得到所有pdf圖件。以Reservoir層nopay中P-Impedance、Porositydensity>WaterSaturation分析為例，如圖32所示。

F-limpadanciaiPesterstlydensity^WaterSoturatlonPorositydensitypdf分布P-Impedance寫Porositydensity云變換關(guān)系P-Impedancepdf分布圖32②stats.py文獻：以上分析成果詳細數(shù)據(jù)都存儲在該文獻內(nèi)，在協(xié)模仿過程中調(diào)用就是該文獻中存儲云變換關(guān)系。以Reservoir層nopay巖性為例顯示該文獻某些內(nèi)容，如圖33所示。I4h*x-d't -' 匸Porositydensitypdf分布P-Impedance寫Porositydensity云變換關(guān)系P-Impedancepdf分布圖32②stats.py文獻：以上分析成果詳細數(shù)據(jù)都存儲在該文獻內(nèi)，在協(xié)模仿過程中調(diào)用就是該文獻中存儲云變換關(guān)系。以Reservoir層nopay巖性為例顯示該文獻某些內(nèi)容，如圖33所示。I4h*x-d't -' 匸,*J.1thc-lc-g^—111fho_t<uc1r]ab41—1ncpdbp"*i: ntJ.l*CaL：x*?l^tlonPDr(ak?xL^P-Tmp^Jjbn(rPcTCtdbyd*ncdbp*j,""MhE^v euvatLan.ur?n?n11avaL1Jhb]quancilec■I■I'-^r；■;'.>Eiirnqi?mih-s-fF-Inpidanc-±!-j2u-s-jzA-.n-e、：t■:<i.fJ-JJfJTbSJij：>飛 . 、以八- , Reservoir層nopay巖性P-Impedancepdf函數(shù)分E.LMXdBEw.、 ='：tl.J50,45j：Jb|-i :zg.agM、. 析成果，以函數(shù)(Q,X)表達，含義為P-■v.A^14Jt-QJ■■■:■ :Impedance不大于X概率為Qj0。泓Cd.L^br(d-2bU^<a-37sB(□.boo,9.日丁丄州956州9泡mf沿r:H性W、1Si.、丄!>Ohl4iri0.324^123jr:0.134^2L-4ir:fl.2^2590-3I「Q.ST&O223I?l.JHI:^inl/ t丄gLQ?<2“4臼里Ir2￡d.O!>ar口玷I,V,.L±J3r0.Irf0ajm.stJO2j-1IP:I：a.37!>?a.J2H2Sj\<a-hOOa0.g"ll73J,.J?：a-lIrReservoir層nopay巖性Porositydensitypdf函數(shù)分析成果，以函數(shù)（Q,X）表達，含義為Porositydensity不大于X概率為QReservoir層nopay巖性WaterSaturationpdf函數(shù)分析成果，以函數(shù)（Q,X）表達，含義為WaterSaturation不大于X概率為Q以矩陣表達P-Impedance、Porositydensity>WaterSaturation兩兩之間有關(guān)關(guān)系，即云變換關(guān)系圖33>StatModMC主界面一Mode一Cosimulation，彈出Cosimulation主界面。這時顧客也許會驚喜地發(fā)現(xiàn)這個界面與GenerateStats>Simulation>Inversion界面基本相似，因此它操作也和以上三個界面那么簡樸。在界面右邊點擊Cosimulation選項，彈出小界面AddCosimulationProperty，如圖34所示。在第一行Layer選取目的層，如Reservoir層；在第二行ContinuousProperties中運用右邊Edit選取將要通過協(xié)模仿得到儲層屬性，如Porositydensity>WaterSaturation；然后在第三行與第四行DiscreteProperty與ContinuousProperties中分別選取巖性曲線與P-Impedance，表達要從反演得到巖性體與P-Impedance體來求得Porositydensity>WaterSaturation體。在Histogram一行應(yīng)選取上一步分析產(chǎn)生stats.py文獻。在Variogram一行填寫Porositydensity>WaterSaturation變差函數(shù)，詳細值獲取方式與P-Impedance基本一致，不再贅述。在中間Inputs某些，輸入巖性體Lithology.hor與P-Impedance.hor，表達要從這兩個體中計算Porositydensity與WaterSaturation，而選取巖性體與P-Impedance體必要是上一步中反演得到比較抱負成果。選中其中一行，然后運用旁邊Edit選取相應(yīng)文獻。最后選井。點擊OK退出。這時Cosimulation主界面上半部如圖35所示。Edit.LayerEdit.CoiTtinuouspr-DpertiesFo-rosity WattrSaturationf「MOdtc「砒eprcperiji *Edit,Dat-aHistograaLoad %~tarti7istarit：g.py*'■Edit...1i 1V&rLogiw§fiussian(B.01,AM,8M)EditInputs-millOicWry14 Edit.-.|HellsincurrentSelected rillsML ，well/「《相■卩L*/357*11.all圖34圖35填寫好界面后，在RunParameters某些填寫想要產(chǎn)生實現(xiàn)個數(shù)，并給出輸出途徑與腳本名。最后點擊GenerateScript產(chǎn)生腳本，并運營腳本，得到相應(yīng)儲層屬性體。接著可以變換巖性體與P-Impedance體來求儲層屬性體，在計算過程中通過測試變差函數(shù)與設(shè)置計算實現(xiàn)個數(shù)來得到各種儲層屬性體實現(xiàn)。最后對這些實現(xiàn)進行質(zhì)控。（4）檢查并質(zhì)控成果。＞比較輸入輸出云變換。原則是輸出云變換（既屬性之間交會關(guān)系）應(yīng)當與輸入云變換大體形態(tài)（在第三步中分析得到）比較一致，如圖36所示，橫坐標為P-Impedance，縱坐標為porosity，顏色反映巖性。輸出（1）中云變換要與輸入云變換形態(tài)基本一致，但分布范疇大些，這正反映了協(xié)模仿結(jié)果不擬定性。而輸出（2）中云變換則與輸入差別較大，屬于不太抱負協(xié)模仿成果。＞抽井檢查儲層屬性體，在井點處井上儲層屬性曲線與儲層屬性體匹配得好嗎？如果不好，要考慮與否是之前計算巖性體與P-Impedance體就與井對不好，或者是使用變差函數(shù)有問題。＞觀測儲層屬性體與相應(yīng)巖性體，在某一巖性中儲層屬性變化平緩嗎？如果變化過于平緩，則需要再進一步測試變差函數(shù)。如圖37顯示是將孔隙度體與相應(yīng)巖性體與P-Impedance體進行比較質(zhì)控。PorosityLithologyP-ImpedancePorosityLitho