JASON随机反演的中文手册.doc

第五期用户通讯随机模拟与随机反演 地震 属性随机模拟、岩性指示随机模拟与随机反演：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、直方图与变异图的分析 并通过交会图确定两类数据的相关系数 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 (存入文件 usertypelist.txt) （二）利用 density 与 gamma 曲线建立岩性曲线 (深度域、lithology)： 例如：工区有六口井：1.wll 2.wll 3.wll 4.wll 5.wll 6.wll 经过分析：当 GR<85, density<2300kg/m3 时为 sand 当 GR>85, density.2300kg/m3 时为 shale 1、建立一个函数 （利用 density 与 gamma 曲线建立岩性曲线）： FunctionMod ® Edit ® Functions ® Newsand_shale_from_gamma_and_density Function name: 键入List Sand/shale Output 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 unit Ok¾Use domain range definitionsEdit Add Range: ® Gamma<85 and density<2300 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 ® sand_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 typeAdd 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 files ® 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 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.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 output ® 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/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 需要作如下直方图与变异图的分析： 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_ aitm_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 层段的岩性指示变差图，砂岩与泥岩要分别做，但要 存在同一个 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 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 ® 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_prob.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