jason地质统计学反演手册

1、StatMod MC入门手册Chapter 1. 工作流程Stage 1: 项目准备.5%Stage 2: 地质统计学参数分析.15%Stage 3: 模拟.5%百分数表示每个步骤所用时间占整个项目时间的百分比Stage 4: 反演.50%.10%Stage 5: 协模拟.15%Stage 6: 不确定性分析与风险评估Chapter 2.基本的输入输出数据输入数据输出数据岩性实现岩性概率体属性实现地质统计学参数岩石物理分析地层网格模型地震数据测井曲线Chapter 3. 详细操作步骤操作步骤以StatMod MC培训数据为例第一步.首先完成一个高质量的叠后CSSI反演这一步的目的是为地质统计学

2、提供一个好的研究基础, 这个“好”主要体现在：(1) 好的井震标定, 目标区的相关值达到0.85以上;(2) 好的叠后反演结果, 用来质控地质统计学模拟和反演结果, 是地质统计学反演结果横向预测准确度的参照物;(3) 利用叠后反演结果进行砂体雕刻, 对目标区的岩性展布、比例有一个总体上正确的把握, 这些认识都是地质统计学的初始输入。（说明：在提供的培训数据中已经为用户做了以上准备，用户可以从主界面中打开该培训数据所在工区, 然后用Map View看工区底图，用Section View查看地震数据、叠后CSSI反演数据、地质框架模型与层位数据以及井数据与子波 , 并用Well Editor检查井

3、震标定情况）第二步. 数据准备 井曲线重采样这一步将测井数据重采样至地质微层采样间隔，具体操作为：(1) JGW主界面 Analysis Processing toolkit；(2) Input Data selection Data type：选Well, 点击Input file(s)右边List选择任意井(可以选多井),然后在弹出的界面Select logs中选择任意井曲线(可以多选),点击OK退出;(3) Parameters Resample log, 在弹出界面Processing toolkit中填写重采样间隔(注意s 与ms单位), 点击OK退出;(4) Output Defi

4、ne process, 从Select from中选择Resample log, 点击输入到右边的Process里面;(5) Output Generate, 在弹出的界面中填写输出路径和输出文件名,然后点击Generate,开始计算重采样的曲线。 计算岩性曲线 ：这里说的岩性曲线和测井解释的岩性不同, 而是指能通过弹性属性(主要指纵波阻抗) 区分出来的岩性, 所以一般可分两到三种岩性。练习里通过交会图的方式教授如何划分岩性。 地震数据：线、道间距尽量一致，如果不一致需要进行跳道处理，避免往反演结果中引入各向异性。 地质框架模型：微层采样间隔与需要分辨的薄层厚度要匹配。第三步. 地质统计学参数

5、分析这里说的地质统计学参数主要指三个参数：概率密度函数(probability density function, 简称pdf，描述某一属性在空间的概率分布情况)、变差函数（描述某一属性随距离的变化，是距离的函数）、云变换（描述两个属性之间的相关关系）。下面先介绍前两个参数: 概率密度函数(1) JGW主界面 Modeling StatMod MC, 弹出StatMod MC主界面，主要填写内容如下：图 1 Solid Model： 点击Input，再点击Solid model，选择地质框架模型，如solidModel.smgLayers ： 点击Input，再点击Layers, 选择相应目标

6、层，如TopPad、Reservoir、BasePad, 利用输入至右边的界面中,并且修改Spacing值,一 般目标层厚度设为与需要分辨的薄层时间厚度值，而非目标层可以设置得稀一些，如与地震数据采样一致 Trace Gate ： 选择反演道门，如wholeProject QC Trace Gate ： 选择QC道门，如allWellsLine 填好后如图1所示。(2) StatMod MC主界面Mode, 共有四个选项：Generate Stats、Simulation、Inversion、Cosimulation，反映的其实就是整个流程。在这一步里主要用Generate Stats来分析p

7、df函数（即概率密度函数）。点击Generate Stats，弹出界面（如图2所示）。 图2 界面上部由空白栏和右边的控键组成。根据分析的需要，选择不同的控键往空白栏处充填内 容。在这一步只需用到最上边的一个控键Base/Joint，表示分析某一属性的pdf函数；界面中部选择Output directory填写输出路径；界面下方选择Generate and run产生脚本并运行；现在开始分层位分岩性分属性地填写空白栏内容： TopPad层：点击Base/Joint, 弹出小界面，如图3所示。图3点击Layer对应空白栏右边的三角下拉菜单，选择相应层位，如TopPad；点击Discrete Pr

8、operty对应空白栏右边的三角下拉菜单,选择相应离散属性,一般指岩性属性（对于非目标层，这一项可以不选）；点击Continuous Property右边的Edit,选择P-Impedance；在Wells一行点击Edit选择重采样的井路径，然后在Wells in current directory中选择相应井，输入至Selected wells；点击OK退出编辑状态。这时在Generate Stats主界面中就会出现相应的内容。 BasePad层：同TopPad层。 Reservoir层：点击Base/Joint，Layer选Reservoir, Discrete Property选定义好的

9、岩性曲线, 点击Continuous Property右边的Edit,选择P-Impedance。选井。填写完后Generate Stats界面空白栏显示如图4所示。 图4还可以通过选择界面中的Copy，选中所有的行，最后点击Paste，选择Paste wells，可以将井复制到选中的位置；最后在Generate Stats界面下半部分的Output directory 中填写相应的路径，然后点击Generate and run 产生脚本并运行。(3) 查看分析结果。在Linux系统中打开网络流览器Konqueror, 并在Location中输入相应的脚本地址。在该路径下除了脚本还有很多运行脚

10、本产生的重要文件：Ø Index.html文件：里面是分层位分岩性分属性分析得到的所有pdf图件。以Reservoir层nopay岩性中P-Impedance的分析为例，如图5所示。分析井段穿过Reservoir层nopay岩性的样本点统计结果Reservoir层中nopay岩性中的P-Impedance的pdf函数分析结果图5图5中显示了Reservoir层中nopay岩性中P-Impedance的pdf函数分布。其它分析结果与之相似，为了直观地了解分析结果，用户应该仔细检查每一层每一岩性中各种属性的pdf函数分析结果图件。Ø *.png文件：一系列*.png文件, 实际

11、是把Index.html文件里的图件分开来显示, 一个*.png文件就是一张图, 这些文件方便用户在做报告或文档时直接插入感兴趣的图件, 而不用再去抓图。Ø scriptCopy.py文件：就是根据在Generate Stats界面中填写的参数产生的脚本文件，有点类似程序语言，但是相对简单多了。用户应该逐渐熟悉脚本文件，先尝试阅读脚本，领会每部分的意思。在逐渐熟悉后可以尝试修改脚本或自己编写。足够熟练后甚至可以跨过填写界面这一步，直接编写脚本并运行，来完成地质统计学参数分析、模拟、反演等工作。图6是该工区地质统计学参数分析脚本例子。输入地层框架模型输入分析道门准备信息输入QC道门输入

12、分析层段以及设定相应的采样间隔BasePad层为非目标层，对其中的P-Impedance不用分岩性分析开始逐层分析Reservoir层为目标层，对其中的P-Impedance应该分岩性进行pdf函数的分析paynopay 图6 Ø stats.py文件：以上分析结果的具体数据都存放在该文件内，可以说是最终的结果文件，以后在做模拟与反演的过程中调用的就是该文件。以Reservoir层的nopay岩性为例显示该文件的部分内容，如图7所示。Reservoir层nopay岩性的P-Impedance的pdf函数的分析结果，以函数（Q,X）表示，含义为P-Impedance小于X的概率为Q图7

13、l 变差函数的求取(1) JGW主界面 Modeling Specials StatMod StatMod Analysis, 弹出StatMod Analysis界面；(2) Input Time/Depth mode: Time Solid model ： SOLID Layers ：点击Select layers选择目标层位，如Reservoir层，在Selected interval设置采样间隔为0.001s，点击Assign current interval to selected layers应用修改的采样间隔 Lithology masks ：先把Enable lithology

14、 masking左边的小方框点亮，再点击Lithology右边的List，在弹出的界面中选择事先划分好的岩性曲线，如Lithology, 这时在Select lithotypes to use中选择某一岩性，如pay（其它岩性稍后也需要分析，只是一次只分析一种岩性中的属性变差） Data for histograms and transforms： 在弹出界面的上半部Data的下方的Well file(s) 右边List中在./WELLS/RESAMPLED选择相应的分析井曲线，如P-Impedance。 Trace gate : wholeProjectEdit Transforms, 弹出

15、两个界面，一个是主界面StatMod Transform，一个是辅界面Transform, 点击辅界面Transform最上方的Histogram, 弹出一个小界面StatMod Histogram 并带着一个不停闪烁的红色提示Build the histogram, 点击该红色提示，则主界面出现相应的样本点分析，调试小界面中的Nr of intervals的值, 尽量使pdf函数样本点分布规则连续。调试好后点击Dismiss退出。点击辅界面Transform上方第二个选项Transform, 弹出界面如图8，同时在主界面中出现一条蓝色拟合线，如图9。用鼠标左键调整蓝色拟合线，使得该曲线形态能

16、大致反映样本点的分布形态，最后用ctrl+鼠标左键停止编辑蓝色拟合线。然后在图8显示的界面下方选择Save as first, 在弹出界面中给出名字，如pay_p-imp_histogram，再连续点三次OK回到StatMod Transform主界面, 点击Dismiss退出界面。 辅界面主界面 图8 图9 (3) Input Data for variogram sampling and modeling Primary data Well log data, 在 弹出界面中的上部Data部分的Well file(s)右边的List中选择相应的分析井曲线，如P-Impedance, 然后点

17、OK, 这时弹出两个界面，一个StatMod Variogram主界面，一个Variogram辅界面。点击 Variogram 辅界面上方第一个选项Sample Vario, 弹出界面StatMod Sample Variogram Parameters, 如图10所示，先点击在该界面下方不停闪烁的红色提示Build the Sample Variogram, 则主界面中出现变差函数的样本点分析曲线，修改图10界面中的Select Sampling Parameters部分的Intervalm下方空白处的值（即变程的单位距离），则主界面中的变差函数的样本点分析曲线形态发生改变，不断调试Inte

18、rvalm值，使得主界面中的变差函数的样本点分析曲线形态尽量光滑，如图11所示，主要看上半部的红色样本曲线，最后点OK结束编辑。 辅界面主界面图 10 图11 （说明：图11的主界面显示的是进行变差函数分析的样本点曲线，变差函数是三维的，有纵向变差，也有横向上的变差，图11的辅界面下方表示的是将三个方向上的变差函数曲线用三种颜色显示。但是因为这里是利用井上的样本点进行变差函数的分析，一方面井曲线在纵向上的样本点的个数一般能满足统计需要（一般50个样本点就可以得到比较理想的pdf函数和变差函数），另一方面井曲线在横向上的分布是远不能达到分析需要的。所以在这里只相信纵向上的变差函数分析结果，即只需

19、分析红色的曲线）点击Variogram辅界面的选项Model Vario，弹出界面如图12所示，同时把鼠标箭头放在图13的上部，按住左键可以得到一条拟合曲线，可以任意拖动该拟合曲线，使得它能最好地拟合纵向变差样本点，特别是前几个样本点（思考为什么？），ctrl+鼠标左键停止编辑。 图12 图13这时图12显示界面的上部空白栏处会出现变差函数表达式（Type, Sill, X, Y, Z, Power）对应的具体参数。在本例中表示变差函数类型是指数型的，基台值（Sill,即拟合曲线的拐点值）为3.07, 三个方向(x,y,z)的变程均为26.9m（当然就像之前所说的只相信纵向上的变程分析结果）。

20、另外，还可以改变变差函数类型，在图12界面中上部分，Func Type下方，点击小控键，弹出几种类型，一般用到的为Gaussian、Expo两种类型。在原点附近，Gaussian型变化较慢，Expo型变化较快，用户可以分别尝试。如果不满意分析结果，可以在图12中选中上方空白栏里的函数，然后点击下面的Remove, 就可以删掉该函数，然后重新编辑。最后记住测试好的变差函数类型与纵向上的变程（单位为深度域单位m）(4) 依次对Reservoir层的其它岩性类型中的属性（如P-Impedance）进行变差函数分析，并记住测试好的变差函数类型与纵向上的变程。再对其它目的层进行分岩性的属性变差函数测试。

21、过程同上。(5) 在分层确定了不同岩性的变差函数类型与纵向变程后, 横向变程怎么确定呢? 这里提供两中思路: 一是与地质概念相结合, 与地质专家或客户进行交流, 大体掌握工区的沉积模式与岩性展布, 以此来确定一个大概的横向变差范围(只要求一个大概的输入)；一是依据叠后CSSI反演结果进行岩性的地质规模和分布等地分析, 来得到一个大概的横向变差范围。(6) 不管是横向还是纵向上的变差，经过分析得到的都是一个大概的值，对反演的结果起到的是软约束的作用，而起到硬约束影响作用的是地震。第四步.地质统计学模拟-目的是测试第三步中确定的地质统计学参数是否合理(1) StatMod MC主界面 Mode S

22、imulation, 弹出Simulation主界面，该界面上部和Generate Stats界面基本一致。以模拟目标层Reservoir中不同岩性的P-Impedance为例来介绍需要填写的内容。点击Base/Joint，在弹出的小界面中Layer选Reservoir层，Discrete Property选岩性曲线, Continuous中选P-Impedance, OK退出。这时Simulation主界面中出现相关信息。如图14所示。对界面中的内容进行编辑：包括岩性比例、岩性的变差函数、不同岩性中P-Impedance的pdf函数与变差函数，也就是前面几步分析得到的结果在这里都要用到。选井

23、。图14相关参数:Ø 岩性比例：选中Reservoir层第一行，如图14所示，然后点击Edit，弹出界面如 图15所示。在Proportions部分编辑岩性比例（ 说明： 岩性比例缺省值为1.0，本练习工区岩性比例设置可以参考pay ：nopay =13 ：87, 可以从井上目标层段不同岩性的样本点统计得到一个粗略的估计，但仍需测试）Ø 岩性的变差函数：在图15中的Variograms部分, 点击Single右边的Edit，在弹出菜单中依 次 将Type、Vertical range、Lateral range设定为gaussian、0.010、600×600。(

24、 说明： Type、Vertical range均为在第三步中测试到的值, Lateral range可以由第一步得到,即利用叠后反演结果进行砂体雕刻, 以求对目标区的岩性展布、比例有一个总体上正确的把握, 同时结合工区的地质沉积模式等地质方面的信息来指导Lateral range的输入。用户可以放心的是Lateral range的输入可以只是一个大概的值，它对真实的各种岩性体的分布起到的是软约束的作用，在后面的反演部分要讲到，对岩性体的分布起硬约束作用的是地震数据。并且Type、Vertical range、Lateral range的设定也不是一蹴而就的，需要调试，寻找合理值的范围)编辑岩

25、性比例编辑属性的变差函数编辑岩性变差函数图15 图16Ø 不同岩性中P-Impedance的pdf函数与变差函数：在Simulation主界面中选中Reservoir目 标层某一岩性对应的行，如选择Lithology.nopay所在行，然后点击右边控键栏Property下方的Edit，弹出界面如图16所示，只需对红框圈出的部分进行编辑。点击Histogram对应的Edit，选择相应的存放pdf函数分析结果的文件。大家还记得这个文件在哪吗？就在第三步做地质统计学参数分析时，在Generate Stats界面运行脚本产生的stats.py文件，里面存放所有目标层中不同岩性中的P-Impe

26、dance等属性的pdf函数分析结果。图7展示stats.py文件的部分内容。大家需要好好回忆一下，整理好头绪。然后点击Variogram对应的Edit, 在弹出的小界面中依此设定Type、Vertical range、Lateral range: exponential、0.006、400×400。这个结果也是由第三步的分析得到的。连续两次点击OK退出编辑状态。编辑其它岩性中属性的pdf函数与变差函数（其实pdf函数调用的都是stats.py文件，只需编辑变差函数即可）。变差函数值可参考以下表1（例子）： VariogramLithologyTypeVertical rangeLa

27、teral rangepayexponential0.0161200×1200nopayexponential0.0161200×1200 表 1在编辑参数时可以使用主界面中的Copy、Paste功能，实现Histogram与Variogram的快速编辑。界面中下方的Run Parameters部分也有一些选项，在Nr of realizations可以填写想要产生实现的个数，Random seed中的值随之改变。如Nr of realizations填1，表示只产生一个实现，这时Random seed自动将种子点设为0；如Nr of realizations填2， 表示产

28、生两个实现，Random seed自动将种子点设为0,1。 在实际的工作中用户可能需要产生更多的实现。用户可以自己熟悉其它选项。最后别忘了给出输出路径与输出文件名。点击Generate script产生脚本文件，可以从网络流览器中调出脚本文件查看，如图17展示脚本中部分内容。最后点击Run运行该脚本文件。 （留住Simulation主界面） nopayReservoir层中每一岩性中P-Impedance的变差函数 pay图17 (2) 查看结果并质控： 打开网络流览器，查看运行脚本产生Index.html文件，选择该文件，然后在Run report下面选择QC histograms(deta

29、iled), 里面是所有反映模拟信息的图件和输入输出参数统计表,图18展示其中部分内容。nopay中输出P-Impedance的pdf函数的主要参数岩性比例标准偏差平均值nopay中输入P-Impedance的pdf函数的主要参数标准偏差平均值nopay中P-Impedance输入输出的pdf函数曲线pay中P-Impedance输入输出的pdf函数曲线图 18 在该图件文件里面主要检查两方面内容：目标层中每一种岩性内相应属性，如P-Impedance 输入输出的pdf函数形态应该大体一致，主要比较入输出pdf函数的平均值与标准偏差（图18中蓝色曲线表示输入，红色曲线表示输出）；目标层输入输出

30、的岩性比例应该大体一致（偏差在5%左右）。 JGW主界面 Analysis Section View, 先选好tracegate (最好选择在StatMod MC主界面中填入的QC Trace Gate：本例中是MyQCGate), 然后从刚才填写的模拟结果输出路径中调出结果查看。打开这个结果的方式主要有两种：Section View界面 File Open session 进入模拟结果输出路径 选择QC.jvis，这时显示模拟结果，包括岩性体与P-Impedance体；Section View界面 Input Stratigraphy data 进入模拟结果输出路径 选择lithology.

31、hor/P-Impedance查看。质控方法与标准：和叠后反演的结果进行比较，如图19所示。与叠后CSSI反演的结果（右图），在反映砂体的展布、地质体规模、岩性比例等方面大体比较一致，这时的模拟结果是比较合理的。记住相应的模拟参数。 图19 比较输入输出的井曲线。在运行脚本的过程中，自动在相应的目录下产生在井点处抽取的“伪井曲线”（Lithology与 P-Impedance曲线）。检查这些伪井曲线，Jason主界面AnalysisWell Manage中选择相应路径下的井，会发现对应Lithology与P-Impedance，分别有Lithology_in, Lithology_out与P-

32、Impedance_in, P-Impedance_out。*_in表示输入的重采样曲线，而*_out表示输出的从模拟结果中的井点处抽取出来的岩性或属性曲线，如图20所示，蓝色表示输入，红色表示输出，两者偏差很小，这是比较理想的模拟结果。图20 建议多测试几个参数，如改变目标层岩性比例、变差函数的类型和纵横向的变程范围来产生实现，并比较差别，尤其注意在Section View / Map View中与叠后反演结果进行比较，选择令实现与叠后反演结果最相似的参数组合；在比较剖面的同时注意检查生成Index.html文件，比较标准与上面所说的一致。建议做一个列表来纪录试验的参数：Proportion

33、 ValueDP Vario ValueP-impedance Value VarioOutput directorywell-basedGaussian-shortExponential-longwell-basedGaussian-longExponential-longLowerPayGaussian-shortExponential-longLowerPayGaussian-longExponential-longHigerPayGaussian-shortExponential-longHigerPayGaussian-longExponential-longwell-basedGa

34、ussian-longPay: Exp-longNonpay: Exp-short 表 2 第五步.地质统计学反演l 无井约束反演-主要目的是测试地震的权重(1) StatMod MC主界面 Mode Inversion, 弹出Inversion主界面。用户这时可能已经发现Inversion主界面与Generate Stats界面、Simulation主界面基本一致。区别大的地方主要在界面的左下方的Seismic Belief部分。在该部分需要填写地震数据的权重。另外，由于是反演，所以需要在目标层上下加个“盖子”，即加入TopPad层和BasePad层作为边界层（边界层厚度要在半个子波长度以上

35、，以消除反演时子波的旁瓣效应）。这一步得回到StatMod MC主界面选择，把TopPad层和BasePad层的Spacing设为0.004s。打开Inversion主界面，对于目标层的参数填写与在Simulation主界面的填写一致，当然也可以用测试效果更好的参数代替。对于非目标层，不用分岩性反演属性。不要选井。(2) 现在来考虑地震权重。在Inversion主界面左下方的Seismic Belief部分，有Seismic和Wavelet两个控键, 点击Seismic，则弹出选项如图21所示。 图21 在StatMod MC反演模块里，地震的信息是通过设定一定的噪音引入的，或者更好理解一些，

36、可以通过设定信噪比引入。信噪比计算公式如下所示：S/N Ratio (in dB) = 20 * lg (std (seismic)/ std (residual)其中，lg表示10为底的对数；Std表示标准偏差；在该公式中，可以把残差residual理解为噪音（即地震数据与合成记录的差值，或者说是合成记录与地震数据匹配不上的部分），地震数据是一定的，信噪比越高，表示残差越小（既噪音越弱，地震数据与合成记录越相似）；反之亦然。所以可以通过控制信噪比的大小来控制引入的噪音。点击Seismic对应的Edit，选择全叠加地震数据，然后在Noise一栏，首先点击Noise右边下拉菜单，可以看见三个选项

37、value/dB/file，选择dB，表示引入信噪比来控制合成记录与地震的相关性。然后填入具体的信噪比值。那么如何设置这个值呢？可以参考叠后反演时生成的一个文件（现在体会到在做地质统计学反演前需要做叠后CSSI反演的重要性了吧？）。首先从JGW主界面下的Analysis中打开Map View，然后选择Input Horizon,找到路径./INVER_TRACE_PLUS, 在这一目录中有一质控的层位文件qc.hor，选择该文件，然后在弹出的多个层位中选择inverted signal to noise ratio(Signal-to-noise ratio), 然后OK，这时Map View

38、中显示目标区的叠后反演信噪比结果，可以打开色标对照看，如图22所示。 图 22由图 22可看出本工区的信噪比范围大约在525dB之内，且工区中大部分的信噪比在10dB以上。所以我们可以根据叠后反演的信噪比来设置地质统计学反演的信噪比，它不是某一确定的值，而是要在某一数值范围内进行测试。对于本工区建议测试5dB、10 dB、15 dB、20 dB几组值。也就是要产生四个脚本并运行。当然，如果用户感兴趣可以测试更多值，并感受随之带来的变化。填入地震信息后，一定别忘了还需填入子波信息。点击图24中的Wavelet选项，在弹出界面中缺省是用单一子波反演，即Single，点击Single对应的Edit，

39、选择相应子波（在叠后CSSI反演时用的最终子波）。在Run Parameters部分，用户可以自己设置需要计算的实现，因为这步主要测试地震的权重，所以不建议产生太多实现（要花时间）。最后给出输出路径与脚本名，再产生脚本并且运行即可。(3) 检查并质控结果。Ø 检查反演结果并进行抽井检查：从JGW主界面下的Analysis中打开Section View, 先在Input 的Trace gate中选择合适的道门，最好是联井线，这里选QC Trace Gate。然后选反演数据，点击Input，选择Stratigraphic data, 或者是在界面上找到图标并点击,即可查找反演体文件*.h

40、or(到刚才输出的路径中查找)并显示。一个实现对应两个反演结果：P-Impedance.hor和 lithology.hor。调出剖面后会发现地质统计学的反演的结果具有很高的分辨率，但是这一结果是否就是用户想要的正确的结果呢？当然要进行质控。这一步最重要的是把井选进来，进行抽井检测（这已经是Jason在用户中普及得最深刻的质控手段了，不管是叠前还是叠后反演，要判断反演结果的好坏最基本的就是进行抽井检查，地质统计学反演也不例外）。比如你比较反演的岩性体与井上的岩性曲线，那么首先你要看到在井周围反演的岩性体与井是否是比较吻合的，然后是观察井间岩性是否有比较好的连通性，好的反演结果应该能够预测大部分

41、规模较大的地质体并且对于大多数细节部分（如薄层）也能有很好的反映。如果在这一步得不到好的结果，那么就算把井加上作为约束，结果也还是不可靠的。如果结果不满意怎么办？需要继续测试输入参数，如统计学参数，甚至井震标定，直到得到满意的结果为止。如图23为抽井检查反演岩性体的实例。图23 Ø 输入输出pdf函数、岩性比例质控：不管在统计参数测试步骤、模拟步骤还是在反演步 骤，都能在相应路径找到index.html文件。打开网络流览器，在相应路径下选择该文件，然后在Run report下面选择QC histograms(detailed)，能看到反映模拟信息的图件和输入输出参数统计表，如图24显

42、示的是其中部分内容-reservoir层中nopay岩性的P-Impedance的pdf函数、岩性比例输入输出的情况。 图24图中蓝色线与红色线分别为输入输出的reservoir层中nopay的P-Impedance的pdf函数，两者的相似度较高，是比较理想的输入输出情况。图中右边是描述输入输出pdf函数的参数表，其中决定pdf函数形态的参数是平均值Mean、标准偏差StDev，用户应该比较输入输出的平均值与标准偏差；此外还应比较输入输出的岩性比例，即图表中Proportion，两个比较的原则都是输入输出应该大体一致。这部分其实与模拟步骤比较的内容一致。Ø 反演结果与叠后反演结果的比

43、较: 在比较的过程中主要注意两方面内容一是岩性比例二是岩性的形态与连通性。由于使用相同的正演褶积算子来模拟地震数据，所以地质统计学反演的结果应该与叠后反演结果比较一致（这个“一致”主要指大套岩性的规模分布等，而不强调细节，如薄互层），如图25所示，左边是地质统计学反演P-Impedance结果，右边是叠后反演结果，两者对大套岩性体的规模尺度及分布反映都比较一致，区别仅是地质统计学反演结果能反映出更丰富的细节。 图25Ø 岩性体沿层切片检查：比较地质统计学反演岩性体的沿层切片与地震数据的沿层提取的均方根振幅属性。地质统计学反演岩性体的沿层切片上的地质体形态应该与地震的比较一致。如果不一

44、致，得考虑以下两个问题：一是设置的信噪比过低了；二是选择的变差函数与地震信息不符，如变程大于主要地质体的尺度与规模。打开地质统计学反演岩性体的沿层切片的方式为：打开Map View Input Stratal slicer, 在弹出界面中，对应Stratigraphic Data选择相应的反演岩性体，在Layer中选择相应的层，在Microlayer一行用前后翻键查看微层中的岩性分布情况。l 无井约束反演-主要目的是学会产生多个实现并进行统计分析(1) 利用界面或脚本产生多个实现：如果对脚本不熟悉，那么利用界面产生多个实现也是非常简单的，做法其实在前面也提到过。在上一步的反演界面中确定最佳地质

45、统计学参数与信噪比并填写好，然后在界面的右下方Run Parameters部分的第一个选项Nr of realizations填写想要产生的实现个数，如10（上一步中建议为了节省时间，只产生一个实现），然后把鼠标放在第二个选项Random seeds对应的选项内单击，就会出现0，1，2，3.9，这时准备工作就做好了，再点击Generate script产生脚本，感兴趣的用户可以打开该脚本看看脚本中用来产生多个实现的语句是什么？最后点击RUN运行反演。这个过程可能会花较长时间，所以白天做测试，下午离开办公室之前运行脚本是最好的选择。运行完毕后，在相应目录下产生10个路径：seed_0、seed_

46、1、seed_2. 、seed_9, 每个路径下都是一组完整的反演结果，包括P-Impedance.hor和 lithology.hor，以及index.html文件等，用户应该按照上一步讲述的质控步骤检查结果，并注意每个实现之间的差别。(2) 对多个实现进行统计分析：现在产生了多个实现，并且这些实现都是等概率的，那么哪个实现更合理呢？StatMod MC提供了对多个实现进行统计性分析的方法。在反演界面中点击OK退出，StatMod MC主界面 Utilities Summary statistics, 弹出界面，在Directory对应的Edit中选择相应目录，这时在中部左边的方框中出现se

47、ed_0、seed_1、 seed_2 、. 、seed_9十个实现，全选然后用输入到右边方框中，最后在界面下方的Output directory中填写输出路径（缺省为summary），点击RUN运行，运行速度较快。运行结束后在相应路径下查看结果：对于岩性体来说，分析结果包括每一种岩性的概率体（比如四种岩性就有四种岩性的概率体）以及一个极大似然岩性概率体；对于P-Impedance体来说，分析结果主要有四个分别是平均值P-Impedance_mean.hor、最小值P-Impedance_min.hor 、最大值P-Impedance_max.hor 、标准偏差体P-Impedance_std

48、dev.hor。用户可以从相应路径下打开这些结果来查看。(3) 检查并质控结果：Ø 比较各个实现与叠后反演结果：主要比较这些实现（P-Impedance.hor和 lithology.hor）与叠后CSSI反演结果对大套岩性体的规模尺度的反映上是否基本一致，而不用太抠细节。Ø 对岩性概率体进行抽井检查：在Section View中打开岩性概率体，如此说来图26所示，并投上井上的岩性曲线进行抽井检查。左图为泥岩概率体，右图为气砂概率体。红色代表概率为1，即10个实现中一般以上都反映有该种岩性；黑色表示概率为0，即所有实现都没有这种岩性。比较时注意两点：一是在井点处岩性概率体是

49、否与井上的岩性相符，如果不符，那么不确定性有多大？二是井间或远离井点处岩性的不确定性有多大？如果在井点处吻合得很好，而井间不具有连通性，即岩性概率降低，那么要考虑这种变化是否与实际情况相符，如果不相符，那么要考虑加入的信噪比是否合适？选择的岩性横向变程是否过小？ 图 26Ø 平均P-Impedance体与叠后反演结果的比较：准则是两者对大套岩性的规模尺度反应大体一致。Ø 查看极大似然岩性概率体：调出极大岩性概率体并且投上井的岩性曲线进行抽井检查, 在井周比较极大似然岩性概率体与井上岩性曲线的吻合程度, 在井间检查岩性的连通性等等。Ø 在井点处比较井与从平均P-Impedance体中抽出的P-Impedance曲线: 从平均P-Impedance体中抽出的P-Impedance曲线要能够反映井上的P-Impedance的大致趋势。Ø 输入输出pdf函数、岩性比例质控：打开网络流览器，在各个实现对 应的路径下选