《JASON随机反演的中文手册.doc》由会员分享，可在线阅读，更多相关《JASON随机反演的中文手册.doc（12页珍藏版）》请在三一办公上搜索。
1、第五期用户通讯随机模拟与随机反演 地震 属性随机模拟、岩性指示随机模拟与随机反演：1. 需要输入的数据：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、直方图与变异图的分析 并通过交会图确定两
2、类数据的相关系数 5、随机模拟 6、计算多次孔隙度随机模拟结果的平均值 7、计算多次模拟的标准偏差 8、计算小层中的砂岩累计厚度 9、随机反演 （一）建立岩性类型与单位： 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)： 例如：工区有六口井：1.wll 2.wll 3.wll 4.wll 5.wll 6.wll 经过分析：当 GR85, density85, density.2300kg/m3 时为 shale 1、建立一个函数 （利用 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 （六口井全进入窗口）Assign 点亮 1.wll sand_shale_from_gamma_and_density (2) Function: 选Assign parameters 点亮 gamma(sfile) 1.wll Well OkAssign parameters 点亮 density(sfile) 1.wll Well Ok OkCopy 点亮 1.wll 点亮 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 (调颜色) 7、EarthModel builder with TDC （建立时转换的岩性曲线） tdc （三）建立岩性概率模型 (shale_prob.mod, sand_prob.mod) !应针对某个小层(layer) 1、通过交会图了解砂泥岩的阻抗分布范围： Applications Display
9、 and Editors Well edit . 选井，选曲线：Porosity 、 impedancen Cross Plot在交会图上了解砂泥岩在阻抗域的大致分布范围。例如： 5e+6 7e+6 砂 砂、泥 泥 1 1 0 0 2、通过交会图来建立岩性与阻抗转换的概率表： Utilities create / Edit Tabls (1) 建立由阻抗转换成砂岩概率的表（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) 建立由阻抗转换成泥岩概率的表（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、生成岩性概率模型，作为随机指示模拟的趋势约束： 在 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) 同理产生： shale_prob.mod(5) 剖面显示:砂岩概率剖面（sand_pro.mod）、 泥岩概率剖面（shale_prob.mod）， 作为后边随机指示模拟的趋势约束。 （四）直方图与变差图的分析（针对一个小层layer）: 例如：对 c 层作分析：(存入文件：analysis.var) D=0.0005 s 如果作： collocated co-simulationb indicate simulation with trend 需
13、要作如下直方图与变异图的分析： 1、井数据 c_sand_porosity_hist (1535) 井曲线上 c 层段、砂岩的孔隙度直方图，数据范围限制在 c 层段 砂岩孔隙度的范围内（例如：1535） 2、井数据 c_shale_porosity_hist (120) 井曲线上 c 层段、泥岩的孔隙度直方图，数据范围限制在 c 层段 泥岩孔隙度的范围内（例如：120） 3、数据体 c_sand_aitm_hist (1e+068e+06) 阻抗数据体上 c 层段、砂岩的孔隙度直方图，数据范围限制在 c 层段砂岩孔隙度的范围内（例如：1e+068e+06） 4、数据体 c_ shale_ ai
14、tm_hist (4e+061e+07) 阻抗数据体上 c 层段、泥岩的孔隙度直方图，数据范围限制在 c 层段泥岩孔隙度的范围内（例如：4e+061e+07） 5、井数据 c_sand_porosity_var 井曲线上 c 层段、砂岩的孔隙度变差图，数据范围限制在 c 层段 砂岩孔隙度的范围内（例如：1535） 6、井数据 c_shale_porosity_var 井曲线上 c 层段、泥岩的孔隙度变差图，数据范围限制在 c 层段 泥岩孔隙度的范围内（例如：120） 7、 井数据 litho_var sand shale 井曲线上 c 层段的岩性指示变差图，砂岩与泥岩要分别做，但要 存在同一个
15、 litho_var 中 8、两类数据的相关系数 通过交会图确定两类数据的相关系数。利用井曲线做 porosity与 impedance 的交会图，用该图来确定两类数据的相关系数。 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 负相关 正相关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 相关性不好 相关性好 （也可用井上的孔隙度与阻抗数据做交会图） （五）随机模拟： 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_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
