jason地质统计学反演手册.docx

《jason地质统计学反演手册.docx》由会员分享，可在线阅读，更多相关《jason地质统计学反演手册.docx（34页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

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

2、, 这个“好”主要表达在：(1) 好的井震标定, 目标区的相关值到达以上;(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) Outp

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

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

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

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

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

9、t，Layer 选 Reservoir, Discrete Property 选定义好的岩性曲线, 点击 Continuous Property 右边的 Edit,选择 P-Impedance。选井。填写完后Generate Stats 界面空白栏显示如图 4 所示。图 4还可以通过选择界面中的 Copy，选中全部的行，最终点击 Paste，选择 Paste wells，可以将井复制到选中的位置；最终在 Generate Stats 界面下半局部的 Output directory 中填写相应的路径，然后点击Generate and run 产生脚本并运行。(3) 查看分析结果。在 Linux

10、 系统中翻开网络流览器 Konqueror, 并在 Location 中输入相应的脚本地址。在该路径下除了脚本还有很多运行脚本产生的重要文件：文件：里面是分层位分岩性分属性分析得到的全部 pdf 图件。以 Reservoir 层 nopay 岩性中 P-Impedance 的分析为例，如图 5 所示。分析井段穿过Reservoir 层nopay 岩性的样本点统计结果Reservoir 层中 nopay 岩性中的P-Impedance 的 pdf 函数分析结果图 5图 5 中显示了 Reservoir 层中 nopay 岩性中 P-Impedance 的 pdf 函数分布。其它分析结果与之相像，

11、为了直观地了解分析结果，用户应当认真检查每一层每一岩性中各种属性的 pdf 函数分析结果图件。预备信息输入地层框架模型输入分析道门输入 QC 道门输入分析层段以及设定相应的采样间隔BasePad 层为非目标层，对其中的 P-Impedance 不用分岩性分析开头逐层分析*.png 文件：一系列*.png 文件, 实际是把文件里的图件分开来显示, 一个*.png 文件就是一张图, 这些文件便利用户在做报告或文档时直接插入感兴趣的图件, 而不用再去抓图。文件：就是依据在 Generate Stats 界面中填写的参数产生的脚本文件，有点类似程序语言，但是相对简洁多了。用户应当渐渐生疏脚本文件，先尝

12、试阅读脚本，领悟每局部的意思。在渐渐生疏后可以尝试修改脚本或自己编写。足够娴熟后甚至可以跨过填写界面这一步，直接编写脚本并运行，来完成地质统计学参数分析、模拟、反演等工作。图 6 是该工区地质统计学参数分析脚本例子。Reservoir 层为目标层，对其中的P-Impedance 应当分岩性进展 pdf 函数的分析paynopay图 6ReservoirP-Impedance 的 pdf 函数的分析结果，以函数Q,X表示，含义为 P-Impedance 小于 X 的概率为 Q文件：以上分析结果的具体数据都存放在该文件内，可以说是最终的结果文件，以后在做模拟与反演的过程中调用的就是该文件。以 Re

13、servoir 层的 nopay 岩性为例显示该文件的局部内容，如图 7 所示。图 7变差函数的求取(1) JGW 主界面 Modeling Specials StatMod StatMod Analysis, 弹出 StatMod Analysis 界面；(2) Input Time/Depth mode: TimeSolid model： SOLIDLayers：点击 Select layers 选择目标层位， 如 Reservoir 层，在 Selected interval 设置采样间隔为，点击 Assign current interval to selected layers 应用

14、修改的采样间隔Lithology masks ：先把 Enable lithology masking 左边的小方框点亮，再点击Lithology 右边的 List，在弹出的界面中选择事先划分好的岩性曲线，如 Lithology, 这时在 Select lithotypes touse 中选择某一岩性，如 pay其它岩性稍后也需要分析，只是一次只分析一种岩性中的属性变差Data for histograms and transforms ： 在弹出界面的上半部Data 的下方的 Well file(s)右边 List 中在./WELLS/RESAMPLED 选择相应的分析井曲线，如 P-Imp

15、edance。Trace gate:wholeProjectEdit Transforms, 弹出两个界面，一个是主界面 StatMod Transform，一个是辅界面Transform, 点击辅界面 Transform 最上方的 Histogram, 弹出一个小界面 StatMod Histogram 并带着一个不停闪耀的红色提示 Build the histogram , 点击该红色提示，则主界面消灭相应的样本点分析，调试小界面中的 Nr of intervals 的值, 尽量使 pdf 函数样本点分布规章连续。调试好后点击 Dismiss 退出。辅界面点击辅界面 Transform 上

16、方其次个选项 Transform, 弹出界面如图 8，同时在主界面中消灭一条蓝色拟合线，如图 9。用鼠标左键调整蓝色拟合线，使得该曲线形态能大致反映样本点的分布形态，最终用 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

17、 log data, 在弹出界面中的上部 Data 局部的 Well file(s)右边的 List 中选择相应的分析井曲线，如P-辅界面主界面Impedance, 然后点 OK, 这时弹出两个界面，一个 StatMod Variogram 主界面，一个Variogram 辅界面。点击 Variogram 辅界面上方第一个选项 Sample Vario, 弹出界面StatMod Sample Variogram Parameters, 如图 10 所示，先点击在该界面下方不停闪耀的红色提示 Build the Sample Variogram , 则主界面中消灭变差函数的样本点分析曲线，修改图

18、 10 界面中的 Select Sampling Parameters 局部的 Intervalm下方空白处的值即变程的单位距离，则主界面中的变差函数的样本点分析曲线形态发生转变，不断调试Intervalm值，使得主界面中的变差函数的样本点分析曲线形态尽量光滑，如图 11 所 示，主要看上半部的红色样本曲线，最终点 OK 完毕编辑。图 10图 11说明：图 11 的主界面显示的是进展变差函数分析的样本点曲线，变差函数是三维的，有纵向变差，也有横向上的变差，图 11 的辅界面下方表示的是将三个方向上的变差函数曲线用三种颜色显示。但是由于这里是利用井上的样本点进展变差函数的分析，一方面井曲线在纵向

19、上的样本点的个数一般能满足统计需要一般50 个样本点就可以得到比较抱负的pdf 函数和变差函数，另一方面井曲线在横向上的分布是远不能到达分析需要的。所以在这里只信任纵向上的变差函数分析结果，即只需分析红色的曲线点击 Variogram 辅界面的选项 Model Vario，弹出界面如图 12 所示，同时把鼠标箭头放在图 13 的上部，按住左键可以得到一条拟合曲线，可以任意拖动该拟合曲线，使得它能最好地拟合纵向变差样本点，特别是前几个样本点思考为什么， ctrl+鼠标左键停顿编辑。图 12图 13这时图 12 显示界面的上部空白栏处会消灭变差函数表达式 Type, Sill, X, Y, Z,

20、Power对应的具体参数。在本例中表示变差函数类型是指数型的，基台值 Sill,即拟合曲线的拐点值为, 三个方向(x,y,z)的变程均为固然就像之前所说的只信任纵向上的变程分析结果。另外，还可以转变变差函数类型，在图 12 界面中上局部，Func Type 下方，点击小控键，弹出几种类型，一般用到的为 Gaussian、Expo 两种类型。在原点四周，Gaussian 型变化较慢，Expo 型变化较快，用户可以分别尝试。假设不满足分析结果，可以在图 12 中选中上方空白栏里的函数，然后点击下面的 Remove, 就可以删掉该函数，然后重编辑。最终记住测试好的变差函数类型与纵向上的变程单位为深度

21、域单位 m(4) 依次对 Reservoir 层的其它岩性类型中的属性如 P-Impedance进展变差函数分析，并记住测试好的变差函数类型与纵向上的变程。再对其它目的层进展分岩性的属性变差函数 测试。过程同上。(5) 在分层确定了不同岩性的变差函数类型与纵向变程后, 横向变程怎么确定呢 这里供给两中思路: 一是与地质概念相结合, 与地质专家或客户进展沟通, 大体把握工区的沉积模式与岩性展布, 以此来确定一个或许的横向变差范围(只要求一个或许的输入)；一是依据叠后 CSSI 反演结果进展岩性的地质规模和分布等地分析, 来得到一个或许的横向变差范围。(6) 不管是横向还是纵向上的变差，经过分析得

22、到的都是一个或许的值，对反演的结果起到的 是软约束的作用，而起到硬约束影响作用的是地震。第四步.地质统计学模拟目的是测试第三步中确定的地质统计学参数是否合理(1) StatMod MC 主界面 Mode Simulation, 弹出 Simulation 主界面，该界面上部和Generate Stats 界面根本全都。以模拟目标层 Reservoir 中不同岩性的 P-Impedance 为例来介绍需要填写的内容。点击 Base/Joint，在弹出的小界面中 Layer 选 Reservoir 层， Discrete Property 选岩性曲线, Continuous 中选 P-Impeda

23、nce, OK 退出。这时Simulation 主界面中消灭相关信息。如图 14 所示。对界面中的内容进展编辑：包括岩性比例、岩性的变差函数、不同岩性中 P-Impedance 的 pdf 函数与变差函数，也就是前面几步分析得到的结果在这里都要用到。选井。图 14相关参数:岩性比例：选中 Reservoir 层第一行，如图 14 所示，然后点击 Edit，弹出界面如 图 15 所示。在 Proportions 局部编辑岩性比例 说明： 岩性比例缺省值为，本练习工区岩性比例设置可以参考pay ：nopay =13 ：87, 可以从井上目标层段不同岩性的样本点统计得到一个粗略的估量，但仍需测试岩性

24、的变差函数：在图 15 中的 Variograms 局部, 点击 Single 右边的 Edit，在弹出菜单中依次 将 Type、Vertical range、Lateral range 设定为 gaussian、600600。( 说明： Type、Vertical range 均为在第三步中测试到的值, Lateral range 可以由第一步得到,即利用叠后反演结果进展砂体雕刻, 以求对目标区的岩性展布、比例有一个总体上正确的把握, 同时结合工区的地质沉积模式等地质方面的信息来指导Lateral range 的输入。用户可以放心的是Lateral range 的输入可以只是一个或许的值，它

25、对真实的各种岩性体的分布起到的是软约束的作用，在后面的反演局部要讲编辑岩性比例编辑属性的变差函数编辑岩性变差函数到，对岩性体的分布起硬约束作用的是地震数据。并且Type、Vertical range、Lateral range 的设定也不是一蹴而就的，需要调试，查找合理值的范围)图 15图 16不同岩性中 P-Impedance 的 pdf 函数与变差函数：在 Simulation 主界面中选中 Reservoir目标层某一岩性对应的行，如选择所在行，然后点击右边控键栏 Property 下方的 Edit，弹出界面如图 16 所示，只需对红框圈出的局部进展编辑。点击 Histogram 对应的

26、 Edit，选择相应的存放 pdf 函数分析结果的文件。大家还记得这个文件在哪吗就在第三步做地质统计学参数分析时，在 Generate Stats 界面运行脚本产生的文件，里面存放全部目标层中不同岩性中的 P-Impedance 等属性的 pdf 函数分析结果。图 7 呈现文件的局部内容。大家需要好好回忆一下，整理好头绪。然后点击 Variogram 对应的 Edit, 在弹出的小界面中依此设定 Type、Vertical range、Lateral range: exponential、400400。这个结果也是由第三步的分析得到的。连续两次点击 OK 退出编辑状态。编辑其它岩性中属性的 p

27、df 函数与变差函数其实 pdf 函数调用的都是文件，只需编辑变差函数即可。变差函数值可参考以下表 1例子：VariogramTypeVertical rangeLateral rangeLithologypayexponential12001200nopayexponential12001200表 1在编辑参数时可以使用主界面中的 Copy、Paste 功能，实现 Histogram 与 Variogram 的快速编辑。界面中下方的 Run Parameters 局部也有一些选项，在 Nr of realizations 可以填写想要产生实现的个数，Random seed 中的值随之转变。如

28、 Nr of realizations 填 1，表示只产生一个实现，这时 Random seed 自动将种子点设为 0；如 Nr of realizations 填 2， 表示产生两个实现，Random seed 自动将种子点设为 0,1。 在实际的工作中用户可能需要产生更多的实现。用户可以自己生疏其它选项。最终别忘了给出输出路径与输出文件名。点击Generate script 产生脚本文件，可以从网络流览器中调出脚本文件查看，如图 17 呈现脚本中局部内容。最终点击 Run 运行该脚本文件。 留住 Simulation 主界面nopayReservoir 层中每一岩性中 P-Impedanc

29、e的变差函数pay图 17(2) 查看结果并质控：翻开网络流览器，查看运行脚本产生文件，选择该文件，然后在 Run report 下面选择QC histograms(detailed), 里面是全部反映模拟信息的图件和输入输出参数统计表,图18 呈现其中局部内容。的 pdf 函数的主要参数nopay 中 P-Impedance 输入输出的 pdf 函数曲线nopay 中输入 P-Impedance的 pdf 函数的主要参数pay 中 P-Impedance 输入输出的 pdf 函数曲线nopay 中输出 P-Impedance图 18在该图件文件里面主要检查两方面内容：目标层中每一种岩性内相应

30、属性，如P-Impedance 输入输出的 pdf 函数形态应当大体全都，主要比较入输出 pdf 函数的平均值与标准偏差图 18 中蓝色曲线表示输入，红色曲线表示输出；目标层输入输出的岩性比例应当大体全都偏差在 5%左右。JGW 主界面 Analysis Section View, 先选好 tracegate (最好选择在 StatMod MC 主界面中填入的 QC Trace Gate：本例中是 MyQCGate), 然后从刚刚填写的模拟结果输出路径中调出结果查看。翻开这个结果的方式主要有两种： Section View 界面 File Open session 进入模拟结果输出路径 选择，

31、这时显示模拟结果，包括岩性体与 P- Impedance 体；Section View 界面 Input Stratigraphy data 进入模拟结果输出路径 选择 P-Impedance 查看。质控方法与标准：和叠后反演的结果进展比较，如图19 所示。与叠后 CSSI 反演的结果右图，在反映砂体的展布、地质体规模、岩性比例等方面大体比较全都，这时的模拟结果是比较合理的。记住相应的模拟参数。图 19比较输入输出的井曲线。在运行脚本的过程中，自动在相应的名目下产生在井点处抽取 的“伪井曲线”Lithology 与 P-Impedance 曲线。检查这些伪井曲线，Jason 主界面Analys

32、isWell Manage 中选择相应路径下的井，会觉察对应 Lithology 与 P-很小，这是比较抱负的模拟结果。Impedance，分别有 Lithology_in, Lithology_out 与 P-Impedance_in, P- Impedance_out。*_in 表示输入的重采样曲线，而*_out 表示输出的从模拟结果中的井点处抽取出来的岩性或属性曲线，如图 20 所示，蓝色表示输入，红色表示输出，两者偏差图 20建议多测试几个参数，如转变目标层岩性比例、变差函数的类型和纵横向的变程范围来 产生实现，并比较差异，尤其留意在 Section View / Map View中与

33、叠后反演结果进展比较，选择令实现与叠后反演结果最相像的参数组合；在比较剖面的同时留意检查生成文件，比较标准与上面所说的全都。建议做一个列表来纪录试验的参数：Proportion ValueDP Vario ValueP-impedance Value Vario Output directorywell-basedGaussian-shortExponential-longwell-basedGaussian-longExponential-longLowerPayGaussian-shortExponential-longLowerPayGaussian-longExponential-lo

34、ngHigerPayGaussian-shortExponential-longHigerPayGaussian-longExponential-longPay: Exp-longwell-basedGaussian-longNonpay: Exp-short表 2第五步.地质统计学反演无井约束反演主要目的是测试地震的权重(1) StatMod MC 主界面 Mode Inversion, 弹出 Inversion 主界面。用户这时可能已经觉察Inversion 主界面与 Generate Stats 界面、Simulation 主界面根本全都。区分大的地方主要在界面的左下方的 Seismic

35、 Belief 局部。在该局部需要填写地震数据的权重。另外，由于是反演，所以需要在目标层上下加个“盖子”，即参与 TopPad 层和 BasePad 层作为边界层边界层厚度要在半个子波长度以上，以消退反演时子波的旁瓣效应。这一步得回到StatMod MC 主界面选择，把 TopPad 层和 BasePad 层的 Spacing 设为。翻开 Inversion 主界面，对于目标层的参数填写与在 Simulation 主界面的填写全都，固然也可以用测试效果更好的参数代替。对于非目标层，不用分岩性反演属性。不要选井。(2) 现在来考虑地震权重。在 Inversion 主界面左下方的 Seismic

36、Belief 局部，有 Seismic 和Wavelet 两个控键, 点击 Seismic，则弹出选项如图 21 所示。图 21在 StatMod MC 反演模块里，地震的信息是通过设定确定的噪音引入的，或者更好理解一些，可以通过设定信噪比引入。信噪比计算公式如下所示：S/N Ratio (in dB) = 20 * lg (std (seismic)/ std (residual)其中，lg 表示 10 为底的对数；Std 表示标准偏差；在该公式中，可以把残差 residual 理解为噪音即地震数据与合成记录的差值，或者说是合成记录与地震数据匹配不上的局部，地震数据是确定的，信噪比越高，表示

37、残差越小既 噪音越弱，地震数据与合成记录越相像；反之亦然。所以可以通过把握信噪比的大小来控 制引入的噪音。点击 Seismic 对应的 Edit，选择全叠加地震数据，然后在 Noise 一栏，首先点击 Noise 右边下拉菜单，可以观看三个选项 value/dB/file，选择 dB，表示引入信噪比来把握合成记录与地震的相关性。然后填入具体的信噪比值。那么如何设置这个值呢可以参考叠后反演 时生成的一个文件现在体会到在做地质统计学反演前需要做叠后 CSSI 反演的重要性了吧。首先从 JGW 主界面下的 Analysis 中翻开 Map View，然后选择 Input Horizon,找到路径./

38、INVER_TRACE_PLUS, 在这一名目中有一质控的层位文件，选择该文件，然后在弹出的多个层位中选择 inverted signal to noise ratio(Signal-to-noise ratio), 然后 OK， 这时 Map View 中显示目标区的叠后反演信噪比结果，可以翻开色标比照看，如图 22 所示。图 22由图 22 可看出本工区的信噪比范围大约在 525dB 之内，且工区中大局部的信噪比在 10dB 以上。所以我们可以依据叠后反演的信噪比来设置地质统计学反演的信噪比，它不是某一确 定的值，而是要在某一数值范围内进展测试。对于本工区建议测试 5dB、10 dB、15

39、 dB、20 dB 几组值。也就是要产生四个脚本并运行。固然，假设用户感兴趣可以测试更多值，并感受随之带来的变化。填入地震信息后，确定别忘了还需填入子波信息。点击图 24 中的 Wavelet 选项，在弹出界面中缺省是用单一子波反演，即 Single，点击 Single 对应的 Edit，选择相应子波在叠后CSSI 反演时用的最终子波。在 Run Parameters 局部，用户可以自己设置需要计算的实现，由于这步主要测试地震的权重，所以不建议产生太多实现要花时间。最终给出输出路径与脚本名，再产生脚本并且 运行即可。(3) 检查并质控结果。检查反演结果并进展抽井检查： 从 JGW 主界面下的

40、Analysis 中翻开 Section View, 先在Input 的 Trace gate 中选择适宜的道门，最好是联井线，这里选 QC Trace Gate。然后选反演数据，点击 Input，选择 Stratigraphic data, 或者是在界面上找到图标 并点击,即可查找反演体文件*.hor(到刚刚输出的路径中查找)并显示。一个实现对应两个反演结果：P-和 。调出剖面后会觉察地质统计学的反演的结果具有很高的区分率，但是这一结果是否就是用户想要的正确的结果呢固然要进展质控。这一步最重要的是把井选进来，进展抽井检测这已经是 Jason 在用户中普及得最深刻的质控手段了，不管是叠前还是叠

41、后反演，要推断反演结果的好坏最根本的就是进展抽井检查，地质统计学反演也不例外。比方你比较反演的岩性体与井上的岩性曲线，那么首先你要看到在井四周反演的岩性体与井是否是比较吻合的，然后是观看井间岩性是否有比较好的连通性，好的反演结果应当能够推想大局部规模较大的地质体并且对于大多数细节局部如薄层也能有很好的反映。假设在这一步得不到好的结果，那么就算把井加上作为约束，结果也还是不行靠的。假设结果不满足怎么办需要连续测试输入参数，如统计学参数，甚至井震标定，直到得到满足的结果为止。如图 23 为抽井检查反演岩性体的实例。图 23输入输出 pdf 函数、岩性比例质控：不管在统计参数测试步骤、模拟步骤还是在

42、反演步骤，都能在相应路径找到文件。翻开网络流览器，在相应路径下选择该文件，然后在 Run report 下面选择 QC histograms(detailed)，能看到反映模拟信息的图件和输入输出参数统计表，如图 24 显示的是其中局部内容-reservoir 层中 nopay 岩性的 P-Impedance 的 pdf 函数、岩性比例输入输出的状况。图 24图中蓝色线与红色线分别为输入输出的 reservoir 层中 nopay 的 P-Impedance 的 pdf 函数，两者的相像度较高，是比较抱负的输入输出状况。图中右边是描述输入输出 pdf 函数的参数表，其中打算 pdf 函数形态的

43、参数是平均值 Mean、标准偏差 StDev，用户应当比较输入输出的平均值与标准偏差；此外还应比较输入输出的岩性比例，即图表中Proportion，两个比较的原则都是输入输出应当大体全都。这局部其实与模拟步骤比较 的内容全都。反演结果与叠后反演结果的比较: 在比较的过程中主要留意两方面内容一是岩性比例二是岩性的形态与连通性。由于使用一样的正演褶积算子来模拟地震数据，所以地质统计学反 演的结果应当与叠后反演结果比较全都这个“全都”主要指大套岩性的规模分布等，而 不强调细节，如薄互层，如图 25 所示，左边是地质统计学反演 P-Impedance 结果，右边是叠后反演结果，两者对大套岩性体的规模尺

44、度及分布反映都比较全都，区分仅是地质 统计学反演结果能反映出更丰富的细节。图 25岩性体沿层切片检查： 比较地质统计学反演岩性体的沿层切片与地震数据的沿层提取的均方根振幅属性。地质统计学反演岩性体的沿层切片上的地质体形态应当与地震的比较全都。假设不全都，得考虑以下两个问题：一是设置的信噪比过低了；二是选择的变差函数 与地震信息不符，如变程大于主要地质体的尺度与规模。翻开地质统计学反演岩性体的沿层切片的方式为：翻开 Map View Input Stratal slicer, 在弹出界面中，对应StratigraphicData 选择相应的反演岩性体，在 Layer 中选择相应的层，在 Microlayer 一行用前后翻