jason地质统计学反演基础手册

Stage1：项目准备……Stage2：地质记录学参数分析..Stage3：模仿..……Stage4：反演……Stage5：协模仿……Stage6：不拟定性分析与风险评估…StatModMC入门手册Chapter1.工作流程………………...5%.….………………..15%………………..5%百分数表达每个环节所用时………………...50%间占整个项目时间比例………………...10%.………………….15%Chapter2.基本输入输出数据Chapter3.详细操作环节操作环节以StatModMC培训数据为例第一步.一方面完毕一种高质量叠后CSSI反演这一步目是为地质记录学提供一种好研究基本，这个“好”重要体当前：(1)好井震标定，目的区有关值达到0.85以上;(2)好叠后反演成果，用来质控地质记录学模仿和反演成果，是地质记录学反演成果横向预测精确度参照物;(3)运用叠后反演成果进行砂体雕刻，对目的区岩性展布、比例有一种总体上对的把握，这些结识都是地质记录学初始输入。（阐明：在提供培训数据中已经为顾客做了以上准备，顾客可以从主界面中打开该培训数据所在工区，然后用MapView看工区底图，用SectionView查看地震数据、叠后CSSI反演数据、地质框架模型与层位数据以及井数据与子波，并用WellEditor检查井震标定状况）输入数据地层网格模型岩石物理分析地质记录学参数地震数据测井曲线输出数据岩性实现属性实现岩性概率体第二步.数据准备●●井曲线重采样这一步将测井数据重采样至地质微层采样间隔，详细操作为：(1)JGW主界面→Analysis→Processingtoolkit；(2)Input→Dataselection→Datatype：选Well，点击Inputfile(s)右边List选取任意井(可以选多井),然后在弹出界面Selectlogs中选取任意井曲线(可以多选),点击OK退出;(3)Parameters→Resamplelog，在弹出界面Processingtoolkit中填写重采样间隔(注意s与ms单位)，点击OK退出;(4)Output→Defineprocess，从Selectfrom中选取Resamplelog，点击››输入到右边Process里面;(5)Output→Generate，在弹出界面中填写输出途径和输出文献名,然后点击Generate,开始计算重采样曲线。●●计算岩性曲线：这里说岩性曲线和测井解释岩性不同，而是指能通过弹性属性(重要指纵波阻抗)区别出来岩性，因此普通可分两到三种岩性。练习里通过交会图方式专家如何划分岩性。●●地震数据：线、道间距尽量一致，如果不一致需要进行跳道解决，避免往反演成果中引入各向异性。●●地质框架模型：微层采样间隔与需要辨别薄层厚度要匹配。第三步.地质记录学参数分析这里说地质记录学参数重要指三个参数：概率密度函数(probabilitydensityfunction，简称pdf，描述某一属性在空间概率分布状况)、变差函数（描述某一属性随距离变化，是距离函数）、云变换（描述两个属性之间有关关系）。下面先简介前两个参数:●●概率密度函数(1)JGW主界面→Modeling→StatModMC，弹出StatModMC主界面，重要填写内容如下：图1SolidModel：点击Input，再点击Solidmodel，选取地质框架模型，如solidModel.smgLayers：点击Input，再点击Layers，选取相应目的层，如TopPad、Reservoir、BasePad，运用→输入至右边界面中,并且修改Spacing值,一般目的层厚度设为与需要辨别薄层时间厚度值，而非目的层可以设立得稀某些，如与地震数据采样一致TraceGate：选取反演道门，如wholeProjectQCTraceGate：选取QC道门，如allWellsLine填好后如图1所示。(2)StatModMC主界面→Mode，共有四个选项：GenerateStats、Simulation、Inversion、Cosimulation，反映其实就是整个流程。在这一步里重要用GenerateStats来分析pdf函数（即概率密度函数）。点击GenerateStats，弹出界面（如图2所示）。图2界面上部由空白栏和右边控键构成。依照分析需要，选取不同控键往空白栏处充填内容。在这一步只需用到最上边一种控键Base/Joint，表达分析某一属性pdf函数；界面中部选取Outputdirectory填写输出途径；界面下方选取Generateandrun产生脚本并运营；当前开始分层位分岩性分属性地填写空白栏内容：①TopPad层：点击Base/Joint，弹出小界面，如图3所示。图3点击Layer相应空白栏右边三角下拉菜单，选取相应层位，如TopPad；点击DiscreteProperty相应空白栏右边三角下拉菜单,选取相应离散属性,普通指岩性属性（对于非目的层，这一项可以不选）；点击ContinuousProperty右边Edit,选取P-Impedance；在Wells一行点击Edit选取重采样井途径，然后在Wellsincurrentdirectory中选取相应井，输入至Selectedwells；点击OK退出编辑状态。这时在GenerateStats主界面中就会浮现相应内容。②BasePad层：同TopPad层。③Reservoir层：点击Base/Joint，Layer选Reservoir，DiscreteProperty选定义好岩性曲线，点击ContinuousProperty右边Edit,选取P-Impedance。选井。填写完后GenerateStats界面空白栏显示如图4所示。图4还可以通过选取界面中Copy，选中所有行，最后点击Paste，选取Pastewells，可以将井复制到选中位置；最后在GenerateStats界面下半某些Outputdirectory中填写相应途径，然后点击Generateandrun产生脚本并运营。(3)查看分析成果。在Linux系统中打开网络流览器Konqueror，并在Location中输入相应脚本地址。在该途径下除了脚本尚有诸多运营脚本产生重要文献：Index.html文献：里面是分层位分岩性分属性分析得到所有pdf图件。以Reservoir层nopay岩性中P-Impedance分析为例，如图5所示。图5图5中显示了Reservoir层中nopay岩性中P-Impedancepdf函数分布。其他分析成果与之相似，为了直观地理解分析成果，顾客应当仔细检查每一层每一岩性中各种属性pdf函数分析成果图件。*.png文献：一系列*.png文献，实际是把Index.html文献里图件分开来显示，一种*.png文献就是一张图，这些文献以便顾客在做报告或文档时直接插入感兴趣图件，而不用再去抓图。scriptCopy.py文献：就是依照在GenerateStats界面中填写参数产生脚本文献，有点类似程序语言，但是相对简朴多了。顾客应当逐渐熟悉脚本文献，先尝试阅读脚本，领略每某些意思。在逐渐熟悉后可以尝试修改脚本或自己编写。足够纯熟后甚至可以跨过填写界面这一步，直接编写脚本并运营，来完毕地质记录学参数分析、模仿、反演等工作。图6是该工分析井段穿过Reservoir层nopay岩性样本点记录成果Reservoir层中nopay岩性中P-Impedancepdf函数分析成果区地质记录学参数分析脚本例子。准备信息输入地层框架模型输入分析道门输入QC道门输入分析层段以及设定相应采样间隔BasePad层为非目的层，对其中P-Impedance不用分岩性分析开始逐级分析Reservoir层为目的层，对其中P-payImpedance应当分岩性进行pdf函数分析nopay图6stats.py文献：以上分析成果详细数据都存储在该文献内，可以说是最后成果文献，后来在做模仿与反演过程中调用就是该文献。以Reservoir层nopay岩性为例显示该文献某些内容，如图7所示。图7变差函数求取(1)JGW主界面→Modeling→Specials→StatMod→StatModAnalysis，弹出StatModAnalysis界面；(2)Input→Time/Depthmode：TimeSolidmodel：SOLIDReservoir层nopay岩性P-Impedancepdf函数分析成果，以函数（Q,X）表达，含义为P-Impedance不大于X概率为QLayers：点击Selectlayers选取目的层位，如Reservoir层，在Selectedinterval设立采样间隔为0.001s，点击Assigncurrentintervaltoselectedlayers应用修改采样间隔Lithologymasks：先把Enablelithologymasking左边小方框点亮，再点击Lithology右边List，在弹出界面中选取事先划分好岩性曲线，如Lithology，这时在Selectlithotypestouse中选取某一岩性，如pay（其他岩性稍后也需要分析，只是一次只分析一种岩性中属性变差）Dataforhistogramsandtransforms：在弹出界面上半部Data下方Wellfile(s)右边List中在../WELLS/RESAMPLED选取相应分析井曲线，如P-Impedance。Tracegate：wholeProjectEdit→Transforms，弹出两个界面，一种是主界面StatModTransform，一种是辅界面Transform，点击辅界面Transform最上方Histogram，弹出一种小界面StatModHistogram并带着一种不断闪烁红色提示Buildthehistogram，点击该红色提示，则主界面浮现相应样本点分析，调试小界面中Nrofintervals值，尽量使pdf函数样本点分布规则持续。调试好后点击Dismiss退出。点击辅界面Transform上方第二个选项Transform，弹出界面如图8，同步在主界面中浮现一条蓝色拟合线，如图9。用鼠标左键调节蓝色拟合线，使得该曲线形态能大体反映样本点分布形态，最后用ctrl+鼠标左键停止编辑蓝色拟合线。然后在图8显示界面下方选取Saveasfirst，在弹出界面中给出名字，如pay_p-imp_histogram，再持续点三次OK回到StatModTransform主界面，点击Dismiss退出界面。辅界面图8图9(3)Input→Dataforvariogramsamplingandmodeling→Primarydata→dWaetla，l在log弹出界面中上部Data某些Wellfile(s)右边List中选取相应分析井曲线，如P-Impedance，然后点OK，这时弹出两个界面，一种StatModVariogram主界面，一种Variogram辅界面。点击Variogram辅界面上方第一种选项SampleVario，弹出界面StatModSampleVariogramParameters，如图10所示，先点击在该界面下方不断闪烁红色提示BuildtheSampleVariogram，则主界面中浮现变差函数样本点分析曲线，修改图10界面中SelectSamplingParameters某些Interval[m]下方空白处值（即变程单位距离），则主界面中变差函数样本点分析曲线形态发生变化，不断调试Interval[m]值，使得主界面中变差函数样本点分析曲线形态尽量光滑，如图11所示，重要看上半部红色样本曲线，最后点OK结束编辑。图11（阐明：图11主界面显示是进行变差函数分析样本点曲线，变差函数是三维，有纵向变差，也有横向上变差，图11辅界面下方表达是将三个方向上变差函数曲线用三种颜色显示。但是由于这里是运用井上样本点进行变差函数分析，一方面井曲线在纵向上