jason随机模拟与随机反演模块中文教程

第五期用户通讯随机模拟与随机反演地震属性随机模拟、岩性指示随机模拟与随机反演：1.需要输入的数据：1、wells(logsandtops)2、seismic.mod(=4ms…)3、horizons4、wavelets(=4ms…)5、aitm.mod(=4ms…)6、solidmodel(withTDC)2.主要工作流程：1、建立岩性类型与单位2、利用density与gamma曲线建立岩性曲线(lithology)3、建立岩性概率模型(shale_prob.mod,sand_prob.mod)4、直方图与变异图的分析并通过交会图确定两类数据的相关系数5、随机模拟6、计算多次孔隙度随机模拟结果的平均值7、计算多次模拟的标准偏差8、计算小层中的砂岩累计厚度9、随机反演（一）建立岩性类型与单位：UtilitiesProjectmanagementProjectparametersEditTypeandunits…LothologicdatatypeDatatype:Unit:noneLithologytype:shale0Lithologytype:sand1sand/shaleAdd…Editlithp-types…Add…Add…sand_shale_from_gamma_and_densitySinglemodel,horizonorwell-loggammaList…Singlemodel,horizonorwell-logdensityList…Add…Edit…Ok(存入文件usertypelist.txt)（二）利用density与gamma曲线建立岩性曲线(深度域、lithology)：例如：工区有六口井：1.wll2.wll3.wll4.wll5.wll6.wll经过分析：当GR85,density2300kg/m^3时为sand当GR85,density.2300kg/m^3时为shale1、建立一个函数（利用density与gamma曲线建立岩性曲线）：FunctionModEditFunctionsNewFunctionname:键入Outputtype:Outputunit:Parameters:Parameterskind选ParametersnameParameterstype选ParametersunitOkParameterskind选ParametersnameParameterstype选ParametersunitOkUsedomainrangedefinitionsRange:Add…Add…Sand/shaleList…nonegammagAPIEdit…densityKg/m^3Edit…1Default0sand/shaleAdd…Add…Assignsand_shale_from_gamma_and_densityGamma85anddensity2300Case键入Assignments:VriablenameExpressiontoassign键入OkOkRange:AssignmentsVriablenameExpressiontoassign键入OkOkOksaveas…sand_shale_from_gamma_and_density.fun2、Time/DepthmodeDepthInputVerticalgateofinputwheredataavaiableVerticalgateofoutputwheredataavaiableHoriontalgatewheredataavaiable3、EditAssignments…(1)FiletypeWellSampleinterval:0.125mOutputfilesettings:DatatypeFilename:1.wllOk2.wll3.wll4.wll5.wll6.wllOk（六口井全进入窗口）点亮1.wll(2)Function:选点亮gamma(sfile)…Well1.wllOkClearAdd…sand_shale_from_gamma_and_densityEdit…sand_shalefrom_gamma_and_densityList…Add…Assignparameters…Add1.wllCopyM点亮density(sfile)…WellOkOk点亮1.wll点亮2.wll3.wll…6.wllOkOk4、EditExistingfiles…Appendanddon’toverwriteoverlap5、OutputGenerater…Ok6、ApplicationsDisplay&editsWelleditor(观察岩性曲线)Color:AttributesPlotattributes…(调颜色)7、EarthModelbuilderwithTDC（建立时转换的岩性曲线）tdc（三）建立岩性概率模型(shale_prob.mod,sand_prob.mod)!!应针对某个小层(layer)1、通过交会图了解砂泥岩的阻抗分布范围：ApplicationsDisplayandEditors…Welledit...选井，选曲线：Porosity、impedanceCrossPlot在交会图上了解砂泥岩在阻抗域的大致分布范围。例如：5e+67e+6砂砂、泥泥11002、通过交会图来建立岩性与阻抗转换的概率表：Utilitiescreate/EditTabls…(1)建立由阻抗转换成砂岩概率的表（ai_to_sand_pro.tab）：X_nameY_nameX_typeY_typeSaveas…ai_to_sand_prob.tabOkAIsand_probimpedanceprobability5e+67e+610Assignparameters…AIshale_probsand_prob.modAssign(2)建立由阻抗转换成泥岩概率的表（ai_to_shale_prob.tab):X_nameY_nameX_typeY_typeSaveas…ai_to_shale_prob.tabOk3、生成岩性概率模型，作为随机指示模拟的趋势约束：在FunctionMod模块中完成(1)InputTime/Depthmode…TimeVerticalgateofinputSolidmodelSelectsolidb,c,dVerticalgateofoutputSolidmodeltdcHorizontalgateTracegate…Settracegate…(2)EditAssignments…Filetype:Sampleinterval:0.001sDatatype:probabilityFilename:键入抹亮sand_prob.modb,c,d抹亮cai_to_sand_prob.tabOk抹亮x(sfile…)Seismic/Propertyaitm_ai.modOk抹亮cOkOkb,dOk(3)OutputGenerate…Ok(4)同理产生：shale_prob.mod(5)剖面显示:砂岩概率剖面（sand_pro.mod）、泥岩概率剖面（shale_prob.mod），impedanceprobability5e+67e+601ClearSeismic/PropertyAdd…AssignAssignparameters..Copy作为后边随机指示模拟的趋势约束。（四）直方图与变差图的分析（针对一个小层layer）:例如：对c层作分析：(存入文件：analysis.var)=0.0005s如果作：collocatedco-simulationbindicatesimulationwithtrend需要作如下直方图与变异图的分析：1、井数据c_sand_porosity_hist(15~35)井曲线上c层段、砂岩的孔隙度直方图，数据范围限制在c层段砂岩孔隙度的范围内（例如：15~35）2、井数据c_shale_porosity_hist(1~20)井曲线上c层段、泥岩的孔隙度直方图，数据范围限制在c层段泥岩孔隙度的范围内（例如：1~20）3、数据体c_sand_aitm_hist(1e+06~8e+06)阻抗数据体上c层段、砂岩的孔隙度直方图，数据范围限制在c层段砂岩孔隙度的范围内（例如：1e+06~8e+06）4、数据体c_shale_aitm_hist(4e+06~1e+07)阻抗数据体上c层段、泥岩的孔隙度直方图，数据范围限制在c层段泥岩孔隙度的范围内（例如：4e+06~1e+07）5、井数据c_sand_porosity_var井曲线上c层段、砂岩的孔隙度变差图，数据范围限制在c层段砂岩孔隙度的范围内（例如：15~35）6、井数据c_shale_porosity_var井曲线上c层段、泥岩的孔隙度变差图，数据范围限制在c层段泥岩孔隙度的范围内（例如：1~20）7、井数据litho_varsandshale井曲线上c层段的岩性指示变差图，砂岩与泥岩要分别做，但要存在同一个litho_var中8、两类数据的相关系数通过交会图确定两类数据的相关系数。利用井曲线做porosity与impedance的交会图，用该图来确定两类数据的相关系数。ApplicationsDisplay&EditorsWelledit..InputWell…选井EditDisplay&Edit选曲线(porosity、impedance)CrossPlotporositysandshaleimpedance负相关正相关相关性不好相关性好（也可用井上的孔隙度与阻抗数据做交会图）（五）随机模拟：StatModStatModModelingInput1.Simulationmode…RunpropertysimulationMode:Gaussiancollocatedco-simulationRunlithologysimulationMode:Indicatorsimulationwithtrend2.Time/Depthmode…timeSolidmodel…tdc3.Layerc0.0005s4.PrimaryinoutdataWelllogdatatdc*.wllimpedance5.SecondaryinputdataModeldataaitm_ai.mod6.Lithologyinputdatawelllogdatatdc*.wllsand/shale7.Variograms/transformsfiles…analysis.var8.Lithologyprobabilitytrends…Trendtype:ModelSelecttrendmodels:Edit9.Tracegate1、Realiationparameters:10Outputsampleinterval(s):0.0005PerformkrigingalsoUseordinarykrigingUsecokrigingsecondnon-biascondition(克里金方法可选）2、Searchparameters:1000（视井距而定）3、Transforms…cshalec_shale_porosity_histcsandc_sand_porosity_histcshalec_shale_aitm_histcsandc_sand_aitm_hist4、V