最新jason地质统计学反演手册资料讲解

StatMod MC入门手册

Chapter 1.工作流程

Chapter 2.基本的输入输出数据

输入数据

输出数据

岩性实现岩性概率体

属性实现

地质统计学参数

岩石物理分析

地层网格模型地震数据

测井曲线

……………………...5%

...….………………..15%

..……………………..5%

……………………...50%

……………………...10%

….………………….15%

百分数表示每个步骤所用时

间占整个项目时间的百分比Stage 4:反演

Stage 2:地质统计学参数分析

Stage 3:模拟

Stage 5:协模拟

Stage 1: 项目准备

Stage 6: 不确定性分析与风险评估

Chapter 3.详细操作步骤

操作步骤以StatMod MC培训数据为例

第一步.首先完成一个高质量的叠后CSSI反演

这一步的目的是为地质统计学提供一个好的研究基础, 这个“好”主要体现在：

(1)好的井震标定, 目标区的相关值达到0.85以上;

(2)好的叠后反演结果, 用来质控地质统计学模拟和反演结果, 是地质统计学反演结果横向预测准确

度的参照物;

(3)利用叠后反演结果进行砂体雕刻, 对目标区的岩性展布、比例有一个总体上正确的把握, 这些认

识都是地质统计学的初始输入。

（说明：在提供的培训数据中已经为用户做了以上准备，用户可以从主界面中打开该培训数据所在工区, 然后用Map View看工区底图，用Section View查看地震数据、叠后CSSI反演数据、地质框架模型与

层位数据以及井数据与子波 , 并用Well Editor检查井震标定情况）

第二步. 数据准备

●●井曲线重采样

这一步将测井数据重采样至地质微层采样间隔，具体操作为：

(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→ Define process, 从Select from中选择Resample log, 点击››输入到右边的Process里面;

(5)Output→ Generate, 在弹出的界面中填写输出路径和输出文件名,然后点击Generate,开始计算重

采样的曲线。

●●计算岩性曲线：这里说的岩性曲线和测井解释的岩性不同, 而是指能通过弹性属性(主要指纵波

阻抗) 区分出来的岩性, 所以一般可分两到三种岩性。练习里通过交会图的方式教授如何划分岩性。

●●地震数据：线、道间距尽量一致，如果不一致需要进行跳道处理，避免往反演结果中引入各向

异性。

●●地质框架模型：微层采样间隔与需要分辨的薄层厚度要匹配。

第三步. 地质统计学参数分析

这里说的地质统计学参数主要指三个参数：概率密度函数(probability density function, 简称pdf，描述某一属性在空间的概率分布情况)、变差函数（描述某一属性随距离的变化，是距离的函数）、云变换（描述两个属性之间的相关关系）。下面先介绍前两个参数:

●●概率密度函数

(1)JGW主界面→ Modeling→ StatMod MC, 弹出StatMod MC主界面，主要填写内容如下：

图 1

Solid Model：点击Input，再点击Solid model，选择地质框架模型，如solidModel.smg Layers ：点击Input，再点击Layers, 选择相应目标层，如TopPad、Reservoir、

BasePad, 利用→输入至右边的界面中,并且修改Spacing值,一般目标层厚度

设为与需要分辨的薄层时间厚度值，而非目标层可以设置得稀一些，如与

地震数据采样一致

Trace Gate ：选择反演道门，如wholeProject

QC Trace Gate ：选择QC道门，如allWellsLine

填好后如图1所示。

(2)StatMod MC主界面→Mode, 共有四个选项：Generate Stats、Simulation、Inversion、

Cosimulation，反映的其实就是整个流程。在这一步里主要用Generate Stats来分析pdf函数

（即概率密度函数）。点击Generate Stats，弹出界面（如图2所示）。

图2 界面上部由空白栏和右边的控键组成。根据分析的需要，选择不同的控键往空白栏处充填内容。在这一步只需用到最上边的一个控键Base/Joint，表示分析某一属性的pdf函数；

界面中部选择Output directory填写输出路径；

界面下方选择Generate and run产生脚本并运行；

现在开始分层位分岩性分属性地填写空白栏内容：

①TopPad层：点击Base/Joint, 弹出小界面，如图3所示。

图3

点击Layer对应空白栏右边的三角下拉菜单，选择相应层位，如TopPad；

点击Discrete Property对应空白栏右边的三角下拉菜单,选择相应离散属性,一般指岩性属

性（对于非目标层，这一项可以不选）；

点击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选定义好的岩性曲线,

点击Continuous Property右边的Edit,选择P-Impedance。选井。填写完后Generate Stats界

面空白栏显示如图4所示。

图4

还可以通过选择界面中的Copy，选中所有的行，最后点击Paste，选择Paste wells，可以将井复制到选中的位置；

最后在Generate Stats界面下半部分的Output directory 中填写相应的路径，然后点击Generate and run 产生脚本并运行。

(3)查看分析结果。在Linux系统中打开网络流览器Konqueror, 并在Location中输入相应的脚本地

址。在该路径下除了脚本还有很多运行脚本产生的重要文件：

➢Index.html文件：里面是分层位分岩性分属性分析得到的所有pdf图件。以Reservoir层nopay岩性中P-Impedance的分析为例，如图5所示。