jason隨機(jī)反演的中文手冊(cè)

1、第五期用戶通訊隨機(jī)模擬與隨機(jī)反演 地震 屬性隨機(jī)模擬、巖性指示隨機(jī)模擬與隨機(jī)反演：1. 需要輸入的數(shù)據(jù)：1、wells (logs and tops)2、seismic.mod (D=4ms )3、horizons4、wavelets (D=4ms )5、aitm.mod (D=4ms )6、solid model (with TDC) 2. 主要工作流程： 1、建立巖性類型與單位 2、利用 density 與 gamma 曲線建立巖性曲線 (lithology) 3、建立巖性概率模型 (shale_prob.mod, sand_prob.mod) 4、直方圖與變異圖的分析 并通過交會(huì)圖確定兩

2、類數(shù)據(jù)的相關(guān)系數(shù) 5、隨機(jī)模擬 6、計(jì)算多次孔隙度隨機(jī)模擬結(jié)果的平均值 7、計(jì)算多次模擬的標(biāo)準(zhǔn)偏差 8、計(jì)算小層中的砂巖累計(jì)厚度 9、隨機(jī)反演 （一）建立巖性類型與單位： Utilities Project management Project parameters Edit Type and units Add Lothologic data typeEdit lithp-types sand/shale Data type : Unit: noneAdd Lithology type : shale 0Add Lithology type : sand 1 Ok (存入文件 usertyp

3、elist.txt) （二）利用 density 與 gamma 曲線建立巖性曲線 (深度域、lithology)： 例如：工區(qū)有六口井：1.wll 2.wll 3.wll 4.wll 5.wll 6.wll 經(jīng)過分析：當(dāng) GR85, density85, density.2300kg/m3 時(shí)為 shale 1、建立一個(gè)函數(shù) （利用 density 與 gamma 曲線建立巖性曲線）： FunctionMod Edit Functions Newsand_shale_from_gamma_and_density Function name: 鍵入List Sand/shale Output

4、type: none Output unit: Parameters :Add Single model,horizon or well-log Parameters kind 選gamma Parameters nameList gamma Parameters type 選Edit gAPI Parameters unit Ok Add Single model,horizon or well-log Parameters kind 選density Parameters nameList density Parameters type 選Edit Kg/m3 Parameters uni

5、t OkUse domain range definitionsEdit Add Range: Gamma85 and density2300 Case 鍵入Add Assignments : sand_shale_from_gamma_and_density Vriablename1 Expression to assign 鍵入 Ok OkEdit Default Range: Add Assignments sand_shale from_gamma_and_density Vriablename0 Expression to assign 鍵入 Ok Ok Ok save as san

6、d_shale_from_gamma_and_density.fun 2、 Time/Depth mode Depth Input Vertical gate of input where data avaiable Vertical gate of output where data avaiable Horiontal gate where data avaiable 3、Edit Assignments (1) File type WellSample interval : 0.125 mClear Output file settings :sand/shale Data typeAd

7、d List File name: 1.wll OkAdd 2.wll 3.wll 4.wll 5.wll 6.wll OkAdd （六口井全進(jìn)入窗口）Assign 點(diǎn)亮 1.wll sand_shale_from_gamma_and_density (2) Function: 選Assign parameters 點(diǎn)亮 gamma(sfile) 1.wll Well OkAssign parameters 點(diǎn)亮 density(sfile) 1.wll Well Ok OkCopy 點(diǎn)亮 1.wll 點(diǎn)亮 2.wll 3.wll 6.wll Ok Ok 4、Edit Existing fil

8、es Append and dont overwrite overlap 5、Output Generater Ok 6、Applications Display & edits Well editor (觀察巖性曲線)M Color: Attributes Plot attributes (調(diào)顏色) 7、EarthModel builder with TDC （建立時(shí)轉(zhuǎn)換的巖性曲線） tdc （三）建立巖性概率模型 (shale_prob.mod, sand_prob.mod) !應(yīng)針對(duì)某個(gè)小層(layer) 1、通過交會(huì)圖了解砂泥巖的阻抗分布范圍： Applications Display

9、 and Editors Well edit . 選井，選曲線：Porosity 、 impedancen Cross Plot在交會(huì)圖上了解砂泥巖在阻抗域的大致分布范圍。例如： 5e+6 7e+6 砂 砂、泥 泥 1 1 0 0 2、通過交會(huì)圖來建立巖性與阻抗轉(zhuǎn)換的概率表： Utilities create / Edit Tabls (1) 建立由阻抗轉(zhuǎn)換成砂巖概率的表（ai_to_sand_pro.tab）：X_name Y_nameAIsand_probX_type Y_typeprobabilityimpedance5e+6 17e+6 0Save as ai_to_sand_prob

10、.tab Ok(2) 建立由阻抗轉(zhuǎn)換成泥巖概率的表（ai_to_shale_prob.tab):X_name Y_nameshale_probAIX_type Y_typeprobabilityimpedance5e+6 07e+6 1Save as ai_to_shale_prob.tab Ok 3、生成巖性概率模型，作為隨機(jī)指示模擬的趨勢(shì)約束： 在 FunctionMod 模塊中完成 (1) Input Time/Depth mode Time Vertical gate of input Solid model Select solid b,c,d Vertical gate of ou

11、tput Solid model tdc Horizontal gate Trace gate Set trace gate (2) Edit Assignments ClearSeismic/PropertyFile type : Sample interval : 0.001 s Data type : probabilitysand_prob.mod File name : 鍵入Add 抹亮 sand_prob.modAssign b,c,dAssign 抹亮 c ai_to_sand_prob.tab OkAssign parameters . 抹亮 x(sfile ) Seismic

12、/Property aitm_ai.modCopy Ok 抹亮 c b,d Ok Ok Ok (3) Output Generate Ok (4) 同理產(chǎn)生： shale_prob.mod(5) 剖面顯示:砂巖概率剖面（sand_pro.mod）、 泥巖概率剖面（shale_prob.mod）， 作為后邊隨機(jī)指示模擬的趨勢(shì)約束。 （四）直方圖與變差圖的分析（針對(duì)一個(gè)小層layer）: 例如：對(duì) c 層作分析：(存入文件：analysis.var) D=0.0005 s 如果作： collocated co-simulationb indicate simulation with trend 需

13、要作如下直方圖與變異圖的分析： 1、井?dāng)?shù)據(jù) c_sand_porosity_hist (1535) 井曲線上 c 層段、砂巖的孔隙度直方圖，數(shù)據(jù)范圍限制在 c 層段 砂巖孔隙度的范圍內(nèi)（例如：1535） 2、井?dāng)?shù)據(jù) c_shale_porosity_hist (120) 井曲線上 c 層段、泥巖的孔隙度直方圖，數(shù)據(jù)范圍限制在 c 層段 泥巖孔隙度的范圍內(nèi)（例如：120） 3、數(shù)據(jù)體 c_sand_aitm_hist (1e+068e+06) 阻抗數(shù)據(jù)體上 c 層段、砂巖的孔隙度直方圖，數(shù)據(jù)范圍限制在 c 層段砂巖孔隙度的范圍內(nèi)（例如：1e+068e+06） 4、數(shù)據(jù)體 c_ shale_ ai

14、tm_hist (4e+061e+07) 阻抗數(shù)據(jù)體上 c 層段、泥巖的孔隙度直方圖，數(shù)據(jù)范圍限制在 c 層段泥巖孔隙度的范圍內(nèi)（例如：4e+061e+07） 5、井?dāng)?shù)據(jù) c_sand_porosity_var 井曲線上 c 層段、砂巖的孔隙度變差圖，數(shù)據(jù)范圍限制在 c 層段 砂巖孔隙度的范圍內(nèi)（例如：1535） 6、井?dāng)?shù)據(jù) c_shale_porosity_var 井曲線上 c 層段、泥巖的孔隙度變差圖，數(shù)據(jù)范圍限制在 c 層段 泥巖孔隙度的范圍內(nèi)（例如：120） 7、 井?dāng)?shù)據(jù) litho_var sand shale 井曲線上 c 層段的巖性指示變差圖，砂巖與泥巖要分別做，但要 存在同一個(gè)

15、 litho_var 中 8、兩類數(shù)據(jù)的相關(guān)系數(shù) 通過交會(huì)圖確定兩類數(shù)據(jù)的相關(guān)系數(shù)。利用井曲線做 porosity與 impedance 的交會(huì)圖，用該圖來確定兩類數(shù)據(jù)的相關(guān)系數(shù)。 Applications Display & Editors Well edit . Input Well 選井Edit Display & Edit 選曲線(porosity、impedance)n Cross Plotp p p p p p p p l porosity p sand shale impedance 負(fù)相關(guān) 正相關(guān)p p p p p p p p p p p p p p p p p p p p p

16、 p p l l l lpp p p p p p p pp l 相關(guān)性不好 相關(guān)性好 （也可用井上的孔隙度與阻抗數(shù)據(jù)做交會(huì)圖） （五）隨機(jī)模擬： StatMod StatMod Modeling Input 1. Simulation mode Run property simulation Mode: Gaussian collocated co-simulation Run lithology simulation Mode: Indicator simulation with trend 2. Time/Depth mode time Solid model tdc 3. Layer c

17、 0.0005 s 4. Primary inout data Well log data tdc*.wll impedance 5. Secondary input data Model data aitm_ai.mod 6. Lithology input data well log data tdc*.wll sand/shale 7. Variograms/transforms files analysis.var 8. Lithology probability trends Trend type : Model Select trend models:shale shale_pro

18、b.modsand sand_prob.mod 9. Trace gate Edit 1、Realiation parameters : 10 Output sample interval (s): 0.0005 Perform kriging also Use ordinary kriging Use cokriging second non-bias condition (克里金方法可選） 2、Search parameters :1000 （視井距而定） 3、Transforms Transforms c shale c_shale_porosity_hist c sand c_sand

19、_porosity_histSecond Transforms c shale c_shale_aitm_hist c sand c_sand_aitm_hist 4、Variograms Variograms c shale c_shale_porosity_var c sand c_sand_porosity_varIndicator variograms c litho_var (包括shale、 sand) 5、Super codes Layer name Super code c c 6、 Correlation coefficients c shale -0.7 c sand -0

20、.7 Output Generate Primary filename: Primary file name : ssm_first_c_poro.mod Generate primary kiging error: ssm_first_c_poro_error.mod Secondary filename : Generate secondary data: ssm_second.mod Generate secondary kriging error: ssm_second_error.mod Lithology filename Lithology file name : ssm_sis

21、.modOk （六）計(jì)算多次孔隙度隨機(jī)模擬結(jié)果的平均值 (mean_poro.mod): FunctionMod Input Time/Depth mode Time Vertical gate of input Solid Vertical gate of output Solid Horiontal gate Tracee gate Edit Assignments File type: Seismic/Property Sample interval: 0.0005 sclear Output file settings : Data type : porosity File name:

22、 mean_poro.mod Add 點(diǎn)亮 mean_poro.mod AssignMean_predefined_mfile Function: 點(diǎn)亮 v(mfile) Assignments parameters Seismic/PropertyFile: (選多次孔隙度模擬的結(jié)果） ssm_c_first_poro_01.mod ssm_c_first_poro_02.mod ssm_c_first_poro_10.mod Add Ok Ok Ok Output Generate Ok （七）計(jì)算多次模擬的標(biāo)準(zhǔn)偏差： FunctionMod Input Time/Depth mode Time Vertical gate of input Solid Vertical gate of output Solid Horiontal gate Trace gate Edit Assignments File type: Seismic/Property Sample interval: 0.0005 sclear Ou