最新jason地质统计学反演手册资料讲解
jason 反演说明1
第一章数据的加载jason是目前最常用的反演软件，它操作上的特点是它需要什么数据或参数就给它什么数据或参数。
下面是它的主窗口（图1）。
图1因为作反演之前已经将坨163区块进行了构造解释，所以可以直接从lanmark中将地震、测井、层位数据导入jason中，操作比较简单。
步骤如下：1、选择工区（即一个文件夹）主窗口——File——Select Project（图2），弹出图3。
选择一个工区，ok。
图2图32、数据的导入主窗口——Datalinks——Landmark——Landmark Link（2003）（图3），弹出图4。
图3图4A 工区的选择File——Seisworks project：选地震工区t163，ok。
（图5）图5File——Openworks project——选SHNEGCAI, 选井列表t163，ok。
（图5）此时，图5 窗口的状态栏将会发生变化，以上选择的工区将会显示。
（图6）图6B 地震数据的导入Select——Import——Seismic/property data（图7）,弹出图8。
选cb 3dv（纯波数据，作反演时一定要用纯波数据），ok。
图7图8C 层位数据的导入Select——Import——Horizons，选择反演时需要的层位和断层（图9）。
图9D：井数据的导入Select（图7）——wells，弹出图10。
选择需要的井，ok。
图10E：数据的传输Transport——Import，以上所选的landmark中的数据将传入jason中。
图11。
Jason反演Wavelets中文手册
Jason反演Wavelets中文手册13．WaveletsWavelets应用程序的目标是：利用各种技术评估地震子波：单独从地震数据利用来自于一口井或多口井的反射系数利用地震数据的部分（partial）叠加或multiple partial 叠加定位井轨迹以优化井和地震数据的关系横向变化时，进行子波插值为建模生成合成子波为得到零相位地震数据计算反褶积评估Q和产生基于子波的Q13.1 快速开始子波估算作为科学，带有艺术性和经验性。
应用一个正确的子波是获取满意的储层描述结果的关键。
为了运行Wavelets，必须下面的数据：◆地震数据◆阻抗井或者选择的位置（Impedance well(s) or selected locaitons）◆选择一个时窗，在时窗内来执行估算◆每口井选择一定数量的地震道或者位置用于估算13.2 菜单导航13.2.1 File 菜单13.2.2 Input 菜单Seismic mode子波能够通过Zero-offset（ZO）地震数据和A V A非Zero-offset地震数据估算而来。
选择你喜欢工作的模式。
Preferred Types选择喜欢的查询输入测井曲线的数据类型。
对于Zero offset模式，它是P-Impedance类型。
对于A VA模式，可以是P-Sonic，S-Sonic和密度类型。
SeismicSeismic 菜单有两部分。
在第一部分中，指定了地震数据文件。
你必须提供一个包含2D或3D数据的地震数据（.mod）。
地震数据必须代表真正振幅偏移反射信息。
在Zero-offset模式中，合成地震记录是通过子波与零偏移反射系数曲线褶积而来的，其中反射系数曲线来源于输入的P-Impedance曲线。
在AV A模式中，角度依赖的反射系数来源于P-Sonic，S-Sonic和密度井曲线，所使用的角度是是在输入的A VA地震数据文件中指定的。
在第二部分中，需要设置地震数据的A V A参数。
Jason 软件详细介绍与基本操作手册
Jason模块功能简单描述Jason地学平台是一个将各学科有用信息综合在一起，即为用户提供丰富的、从油气勘探、开发到生产阶段将地震、测井和地质信息有机结合，进行反演、模拟和预测不同阶段储层岩石物性和流体变化的综合研究工具。
您根据可用信息和经济条件，您可以明智的选择应用哪个“工作流程”更好。
1、Enveriment 运行环境（工区管理地震解释和接口）与其它解释系统紧密相连、数据格式灵活、方便。
2、EarthMod 地质框架模型，包括模型建造器，模型生成器和测井曲线生成器。
把地质、地震、测井曲线（斜井轨迹）、沉积模式（整合，不整合，礁体及河道等）建立以地震道为坐标的模型上，成为储层预测的基础，并赋予地质意义。
3、VelMod 速度模型建立工具VelMod 除具有建立速度模型的基本功能外，且有其独特特征：•地质引导：结合地质信息进行插值运算，能把稀疏的速度数据拓展为三维速度场。
•基于地质概念的速度调整•真三维插值4 、Wavelets 二维和三维地震子波分析●无井、单井、多井提子波和空变子波分析●从斜井中进行真三维分析●估算信噪比谱●为使结果稳定利用模型驱动5、 Modtrace 完整的地震反演系统( 包括对大倾角地层)●所用约束稀疏脉冲反演是基于道的递归反演，产生地震带宽内的波阻抗数据。
●其算法的唯一特征应用解释控制，为获得可靠的低频信息提供工具。
在EarthMod 的基础上建立空变约束条件，其结果突出各向异性，提高了反演分辨率。
●多种方式提取子波，包括无井、单井多道、多井多道、直井、斜井提子波和在EarthMod 的基础上提空变子波。
●由于大量地震信息的反演是基于EarthMod 的基础，所以适用与多条复杂断层的地质情况，使其反演结果更接近真实的地质模型。
●给出产层有效厚度图和孔隙度分布图６、InverMod 精细储层描述（多井）•采用一欧洲石油公司的专利，专门致力于薄层预测和精细描述。
对油田滚动勘探开发十分有用。
地质统计学反演技术浅析
长江大学学报 （ 自然 科 学 版 ） 理 工 ２１ 年 ６ 第 ９ 第 ６ ０２ 月 卷 期 Ｊ ｕｎ ｌ ｆＹ ｎ ｔｅ ｎｖｒｉ （ ａ Ｓｉ ｄｔ Ｓｉ ｎ Ｊｎ ２ １ ，Ｖｏ． ． ｏｒ ａ ｏ ａ ｇｚ ｉ ｓｔ Ｎ ｔ ｃ Ｅ ｉ ｃ＆Ｅ ｇ ｕ ． ０ ２ Ｕ ｅ ｙ ） １９Ｎｏ ６
演 技 术 的基 本 原 理 和 实现 过 程进 行 了阐 述 ， 并 交 待 了反 演 过 程 中 的计 学 反 演 ； 随 机模 拟 ； 储 层 预 测 ；薄 储 层 识 别 关
［ 图分 类 号 ］ Ｐ ２ ． 中 ６８２
［ 献标识码］Ａ 文
者 波 阻抗 无法 描述 的储 层或 油藏 的地 质建模 。
２ 地 质 统 计学 反 演过 程
２ １ 地 质 统 计 学 分 析 ．
地 质统 计学 分析 是 指对 目的层 段所 要模 拟 的属 性进 行概 率 分 布统 计 。首 先 建 立 属性 概 率 密 度 函数 ， 其 次进 行 空间 变差 函数 分析 以确 立空 间上 的结 构关 系 。在 地层 研究 的最 小单 元格 架 内针 对不 同层 段 、不 同岩性 或沉 积微 相及 各种 岩性 内的属性 值进 行统 计分 析 ，得 到具有 地质 意义 的不 同层 段不 同岩性 的变 差 函数及 所模 拟 属性 的变差 函数 。在地质 统计 部分 要 引人尽 量 多 的地 质信 息 ，包 括储层 展 布 的横 向和纵 向 非 均 质性研 究 成果 以及 区域沉 积 的研究 成果 ，使 井 上及地 震 信息统 计得 到 的空 间结构 关 系与前 期 地质 研 究 成果 统一 起来 。只有全 面正 确地 应用 地 质信息 ，所 建立 的地 质统 计模 型才 能更 真实 地 反映储 层 的空 间
jason地震反演
关于波阻抗反演中的低频分量有两个基本的问题：（1）在波阻抗反演中，为什么要构造低频信息？地震资料里为什么没有低频信息？（2）为什么井资料有低频信息？如何用它构造用于波阻抗反演的低频信息大家在实际勘探开发过程中都在讲我们的地下储层薄，地震资料不能满足要求，讲提高分辨率的问题，楼主提了一个很有意思的问题，低频分量。
在这里我想尝试对此进行一个简单分析：1、什么是低频？学过信号处理的人都知道，我们可以把一个波（各种波场都可以）变换到频率域，在频率域频率较低的叫低频.如果用一个通俗的说法，低频就是“趋势”，如果你炒股的话，每天的成交量属于高频，而30日均线上看到的算低频，如果说到天气，下一个小时的气温算高频，但北京今年秋天的平均气温变化算低频。
对于我们搞地学的来讲，地下的压实规律是低频，某套地层的砂地比是多少算低频，某套地层的平均阻抗也算低频。
2、低频有什么用？通过上面的分析可以看出，低频反映的是某种事物或事物的某种特征的变化规律，是更本质的东西。
如果你把握住了股市的低频，你就可以低价买，高价卖。
如果你把握了房价的低频，你也可以发大财。
如果你把握了地下储层变化的低频，你就知道那里砂多，那里泥多，你的勘探成功率也会提高。
3、地震资料有没有低频？通过讲地震资料变换到频率域，你很容易知道，地震资料中没有低频，这主要是由于地震检波器的工作原理决定的。
对模拟检波器来说，<10Hz的低频很难收到。
对现在逐渐应用的数字检波器，小于3Hz的信号也很难被收到。
4、井资料中有没有低频？一般来说，大家在使用井资料的时候，主要用她的高频，看某个砂层的情况，内部岩性、物性、含油气性变化等，但如果测量井段较长，你就会看到低频：比如对渤海湾上第三系地层，馆陶组总体砂多泥少，而明化镇泥多砂少，具体到某口井，明下段又有什么规律，这就是低频。
5、反演中为什么需要低频？反演是我们把得到的地震信号反推回地下地质结构的过程，由于地震资料中缺失低频，自然就需要通过各种方法进行补偿。
Jason反演入门培训教程
第一章数据的加载jason是目前最常用的反演软件，它操作上的特点是它需要什么数据或参数就给它什么数据或参数。
下面是它的主窗口（图1）。
图1因为作反演之前已经将坨163区块进行了构造解释，所以可以直接从lanmark中将地震、测井、层位数据导入jason中，操作比较简单。
步骤如下：1、选择工区（即一个文件夹）主窗口——File——Select Project（图2），弹出图3。
选择一个工区，ok。
图2图32、数据的导入主窗口——Datalinks——Landmark——Landmark Link（2003）（图3），弹出图4。
图3图4A 工区的选择File——Seisworks project：选地震工区t163，ok。
（图5）图5File——Openworks project——选SHNEGCAI, 选井列表t163，ok。
（图5）此时，图5 窗口的状态栏将会发生变化，以上选择的工区将会显示。
（图6）图6B 地震数据的导入Select——Import——Seismic/property data（图7）,弹出图8。
选cb 3dv（纯波数据，作反演时一定要用纯波数据），ok。
图7图8C 层位数据的导入Select——Import——Horizons，选择反演时需要的层位和断层（图9）。
图9D：井数据的导入Select（图7）——wells，弹出图10。
选择需要的井，ok。
图10E：数据的传输Transport——Import，以上所选的landmark中的数据将传入jason中。
图11第二章合成记录的建立在jason上建合成记录的特点是精度高，但随意性大。
建立合成记录的步骤是:井曲线、地震数据、子波的加载，子波的编辑和评价，合成记录的生成和编辑。
1、井曲线、地震数据、子波的输入主窗口——Analysis——Well log editing and seismic tie（图1），弹出图2。
图1图2Input——well，选择要输入的井，例如：T714_e。
地质统计学反演
地质统计学反演（StatMod）一、方法原理JASON的StatMod是一个集多种随机模拟技术的软件包，是以概率论为其理论基础的。
其目的是提供一个或多个在某种概率条件下的，既满足数据的地质统计学特征又满足地质、测井和地震信息的三维储层参数概率模型。
数据的地质统计学特征由数据的概率分布图和变异函数描述。
由于地质统计模拟是基于概率意义上的随机模拟。
为满足概率条件必须有足够多的井资料。
软件要求的已知井数不少于6口。
地质统计学主要的算法是岩性指示模拟和序贯高斯模拟技术（SGS）。
序贯高斯模拟方法是一种产生来自高斯场模型实现的方法。
它基于序贯模拟思想。
该方法首先是将研究区域离散为网格系统，然后序贯地处理每一个网格节点。
由于每个节点处随机变量是服从条件化的正态分布，因此，网格节点值完全由均值和方差两个参数确定。
通过求解克里金方程组就可给出该网格节点处的均值和方差，从而将节点处的正态分布确定下来，并采用相应的抽样方法得到该网格节点处的一个样本。
直至全部网格节点计算完毕。
需要指出的是：求解克里金方程组时的条件数据包括原始数据，先前已模拟的、落在模拟邻域内所有被模拟的网格节点处的值。
序贯高斯模拟方法是一种条件模拟，它保证原始数据和直方图及变异函数都被条件化。
在地震储层预测中，每一道就是一个网格节点。
在模拟过程中，需要求取的最典型的属性是波阻抗和孔隙度。
地质统计学考虑了模拟过程中结果的不唯一因素，故而增加了结果的误差分析。
通过用户定义方式在三维地质模型的每个网格节点上计算出的概率密度函数，可以计算出结果数据体。
概率密度函数是数据体中能够对不确定性进行正常估算的参数分布规律。
地质统计学反演对测井曲线的应用方式，与Jason其它的反演方法截然不同。
Jason其它反演方法再队测井曲线进行应用过程中，除子波估算外。
只有在地震数据中没有低频信息时，在最终的反演结果中才会反映测井曲线的信息。
然而地质统计学将重采样后的测井曲线沿井轨迹复制到三维网格点中。
[讲解]地质统计学反演
![[讲解]地质统计学反演](https://img.taocdn.com/s3/m/cb2ed9ecfbb069dc5022aaea998fcc22bcd14307.png)
地质统计学反演（StatMod）一、方法原理JASON的StatMod是一个集多种随机模拟技术的软件包，是以概率论为其理论基础的。
其目的是提供一个或多个在某种概率条件下的，既满足数据的地质统计学特征又满足地质、测井和地震信息的三维储层参数概率模型。
数据的地质统计学特征由数据的概率分布图和变异函数描述。
由于地质统计模拟是基于概率意义上的随机模拟。
为满足概率条件必须有足够多的井资料。
软件要求的已知井数不少于6口。
地质统计学主要的算法是岩性指示模拟和序贯高斯模拟技术（SGS）。
序贯高斯模拟方法是一种产生来自高斯场模型实现的方法。
它基于序贯模拟思想。
该方法首先是将研究区域离散为网格系统，然后序贯地处理每一个网格节点。
由于每个节点处随机变量是服从条件化的正态分布，因此，网格节点值完全由均值和方差两个参数确定。
通过求解克里金方程组就可给出该网格节点处的均值和方差，从而将节点处的正态分布确定下来，并采用相应的抽样方法得到该网格节点处的一个样本。
直至全部网格节点计算完毕。
需要指出的是：求解克里金方程组时的条件数据包括原始数据，先前已模拟的、落在模拟邻域内所有被模拟的网格节点处的值。
序贯高斯模拟方法是一种条件模拟，它保证原始数据和直方图及变异函数都被条件化。
在地震储层预测中，每一道就是一个网格节点。
在模拟过程中，需要求取的最典型的属性是波阻抗和孔隙度。
地质统计学考虑了模拟过程中结果的不唯一因素，故而增加了结果的误差分析。
通过用户定义方式在三维地质模型的每个网格节点上计算出的概率密度函数，可以计算出结果数据体。
概率密度函数是数据体中能够对不确定性进行正常估算的参数分布规律。
地质统计学反演对测井曲线的应用方式，与Jason其它的反演方法截然不同。
Jason其它反演方法再队测井曲线进行应用过程中，除子波估算外。
只有在地震数据中没有低频信息时，在最终的反演结果中才会反映测井曲线的信息。
然而地质统计学将重采样后的测井曲线沿井轨迹复制到三维网格点中。
jason地震反演软件培训PPT参考课件
线性反演类、非线性反演类
9
反演方法分类
3、基于实现方式
基于反褶积的反演方法：包括基于地层反褶积的道积分、 递归、广义线性等反演方法；基于最大似然反褶积的稀疏脉冲 反演方法;基于最大后验概率准则算法反褶积的模拟退火等。
基于波动方程的反演：Born反散射。 基于随机过程的反演：随机反演、随机模拟等。 基于特征分析的反演：特征反演、神经网络反演等。 其它：混沌反演等
2
国内外研究现状
3
概述
隐蔽性油气藏已成为勘探开发重要目标， 岩性油气藏占有其中很大的比重。对河道砂 体、三角洲浊积砂体、砂砾岩体等复杂岩性 储层进行高精度预测和描述，成为重要的研 究问题。
将具有高横向分辨能力的地震资料、高 垂向分辨能力的测井资料与地质结合起来进
行地震反演，把界面反射型地震资料反演成 岩层单元型波阻抗等地层物性参数，以岩层 为目标进行储层解释。
地震反演软件培训
1
目录
第一章 国内外研究现状 第二章 叠后约束地震反演原理及常用反演方法 第三章 波阻抗反演影响因素 第四章 复杂岩性储层反演关键技术环节 第五章 多种类型复杂岩性储层反演技术对策及应用 第六章 Jason 软件结构 第七章 Geoframe 系统与JASON 系统数据交换 第八章 反演数据的综合研究 第九章 叠前地震纵、横波阻抗（速度）约束反演
融合多种资料和技术所长的约束地震反演处理 与解释应用技术是研究复杂岩性储层变化规律的十 分有效的储层描述新技术。
该项技术通过科学的数学算法和计算机软件， 综合应用地震、测井、地质等资料，并将它们有机 地结合为一体，把地震资料中包含的丰富岩性、物 性、流体信息反演成岩层单元型波阻抗数据或地层 层速度数据以及岩性和孔隙度等地层物性资料。使 其能与钻井、测井资料直接对比分析，以岩层为目 标进行地质、储层解释，在研究储层的空间特征和 厚度变化规律方面具有独特的优势。
jason 反演说明3
第三章模型的构造反演模型的构造是进行反演的基础。
模型的构造可分为两步：模型的建立和模型的产生。
(一)、模型的建立图1Modelling（图1）——Earthmodel——Model builder ( without TDC )，弹出图2。
1、Input :（1）Time/Depth Model:Time,ok.（2）Horizons …→选层文件t163.hor（包括断层）（3）Framework … →第一次做，可不选；如果已有*.frw 文件，可选一个参考,此处选frw.frw。
（这是一个格架）(4) Wells …→选已经编辑调整好的深时转换的井文件（*T*.wll），ctrl +左键（选多井），此处选t712_ed.wll, t714_ed.wll, t717-1_ed.wll, t710_ed.wll(5)Logs …→选多种井曲线(只选三种: p_sonic,density,p_impedance)(6)Trace gate …→选道范围???2、Edit：(1) Edit framewo rk … -------关键( *.frw )注意：a)制表应注意的事项：b)从底层开始向上逐层编辑先建断层下盘的地层，后建断层上盘的地层被断层切割的层不能作为datum有断层时最好加顶、底层注意顶底的接触关系c) Save as …(*.frw)d) Ok最后格架编辑的结果如图3。
3、Output → Generate …(1) 点亮 Run now(2) 选输出内容: 全选(3)Output directory:(4) Interface & layers :产生文件： ../ Test163gz /Tinterface.hor, Tthickness.hor, Ttopbase.hor4、File → Save & Exit5、剖面显示，检查时间域的 井曲线和分层数据（二）模型的产生(二)EarthModel →Model generator …1、Input :(1) Solid model …→ zmod (选Model builder without TDC 的结果)(2)Data for seimic modeling——seimic data&Wavelets(选纯波cb 3dv和子波Creat_Wavelets25.mtr)(3) Trace gate …选整个三维数据体2、Edit :(1)Model parameters …Time sample interval : 0.001 s （0.0005）（采样间隔：用于InverMod 和StatMod 时，采样间隔可小些。
Jason综合地震反演和储层预测
一、Jason综合地震反演和储层预测地质框架模型（Earthmodel）子波提取和分析（wavelets）约束稀疏脉冲反演（invertrace）多参数岩性特征反演（invermod）地质统计随机模型与随机反演（statmod）(一)地质框架模型（Earthmodel）为储层和油气藏定量描述创建一个由地震坐标描述的地质框架模型。
这个模型融合了构造（层位、断层）、地质沉积模式、测井资料和初始权重分布等信息。
地质框架模型是整个Jason地震反演和储层、油气藏定量描述的基础。
1.目的1)构造以层为基础的、用于Jason其它模块的参数描述模型，即地质框架模型。
2)生成以地质框架模型为基础的测井曲线内插模型。
3)提供用于地震反演的低频模型。
4)生成平滑、封闭的层位顶、底、厚度平面图。
2.模块及功能（1）Model builder with/without TDC(模型建造器)用构造、地质、沉积、测井等资料形成一个参数化的时间、深度域的三维封闭模型。
这个参数化的模型包括了层位、断层、地层接触关系（整合、不整合、河道、礁等），测井曲线和基于层位的权重系数。
（2）Model generator（模型生成器）由模型建造器形成的参数化三维封闭模型创建不同测井曲线的三维属性模型（如波阻抗、声波速度、孔隙度、SP等）（3）Model interpolator（模型内插器）在参数化的三维封闭模型控制下生成不同网格密度的三维属性体。
（4）Well curve generator（测井曲线生成器）从三维属性体中抽取任意位置上的测井曲线。
3.特点1)地质框架模型含有生成属性模型所有参数。
2)地质框架模型包括地震、地质、沉积、测井等资料和信息。
3)提供多种内插算法，供用户选择。
(二)子波提取和分析（wavelets）提供用于合成记录与反演的地震子波估算或理论子波计算工具。
在单井、多井或无井的情况下，都可以由单道或多道地震信息估算出最佳地震子波。
地质统计学反演
地质统计学反演一、引言地质学中的反演是指根据地质数据反演出与之相适应的模型或参数，从而探究地球内部地质结构和地球演化历史。
地质统计学反演是一种基于数学方法的地质反演技术，它可以利用地球物理测量数据反演地质模型，是研究地球结构和演化过程的重要手段之一。
本文将介绍地质统计学反演的基本概念、方法和应用。
二、基本概念地质统计学反演是一种基于概率论与统计学的数学方法，它可以将地球内部的物理场量（如重力、磁场、地震波等）与地质结构联系起来，从而推导出地下地质结构的空间分布和特征。
其基本假设是地下地质结构满足一定的随机性和空间相关性，因此可以用统计学的方法进行描述和分析。
三、反演方法地质统计学反演的主要方法包括：贝叶斯反演、最小二乘反演、滤波反演、参数估计反演、基于神经网络的反演等。
（一）贝叶斯反演贝叶斯反演是一种基于概率论的反演方法，它可以通过对先验分布、误差分布和观测数据分布的综合分析，计算出后验分布，并据此推导出反演模型。
贝叶斯反演可以有效地处理噪声信号，适用于高维数据反演和模型参数不确定性较大的情况。
（二）最小二乘反演最小二乘反演是一种基于最小化误差平方和的反演方法，它可以用于处理线性或非线性反演问题，并可用于求解一些简单的地球物理模型。
（三）滤波反演滤波反演是一种基于频谱分析的反演方法，它可以将地球物理数据进行频域分析，通过滤波器对频域数据进行处理，并反演得到地下物质的空间分布。
（四）参数估计反演参数估计反演是一种基于参数估计的反演方法，它将地球物理数据与地下模型之间的关系表示为一组参数化的方程式，并通过估计参数，得到最优的地下模型。
（五）基于神经网络的反演基于神经网络的反演是一种新兴的反演方法，它使用深度学习神经网络模型来反演地球物理数据，可以有效解决高维数据反演与非线性反演问题。
四、应用地质统计学反演已经广泛应用于地球物理勘探、地震监测、矿产勘探、环境监测等领域。
（一）地球物理勘探地质统计学反演可用于地球物理勘探，明确地下地质结构、分层、岩性、溶洞、断裂带等地质信息，从而指导勘探过程。
地质统计学反演流程
地质统计学反演流程一、准备工作。
咱得先收集数据呀，这数据就像是盖房子的砖头一样重要。
包括地震数据、测井数据等等。
地震数据呢，就像是大地给我们的一种特殊信号，它能告诉我们地下大概是个啥情况。
测井数据则更像是从地下各个小角落直接得来的信息，比如地层的岩性呀、孔隙度之类的。
这些数据可得收集全乎了，要是少了点啥，就好像做饭少了调料，做出来的东西味道就不对啦。
而且收集数据的时候要特别小心，可不能弄错了数值，不然后面的工作就全乱套喽。
二、建立地质模型。
有了数据之后呢，就开始建立地质模型啦。
这就像是搭积木一样，不过这个积木可是根据我们收集的数据来搭的哦。
我们要根据对地下地质结构的理解，把不同的地层、岩石类型等在模型里表示出来。
这个过程中要充分发挥我们的想象力，但是又不能天马行空，得符合地质规律。
比如说，砂岩一般会在某个特定的深度范围出现，我们就得按照这个规律来构建模型。
这就好像是给我们要探索的地下世界画一个草图，虽然还不是很精确，但是有个大概的框架了。
三、确定变差函数。
变差函数可重要啦。
它有点像描述地下地质属性变化规律的一把尺子。
我们要根据数据去确定这个变差函数。
比如说，某个地区的地层厚度变化是比较平缓的呢，还是变化得很剧烈的呢？变差函数就能告诉我们这些信息。
这就好比我们知道了一个人的性格特点，就能更好地预测他在不同情况下的反应。
确定变差函数需要我们对数据进行分析和计算，有时候可能得试好几次才能找到最适合的那个变差函数，就像试衣服一样，得找到最合身的那件。
四、进行反演计算。
接下来就是进行反演计算啦。
这就像是一场魔法表演，我们把前面收集的数据、建立的模型还有确定的变差函数都放进一个特殊的算法里，然后这个算法就像一个小魔法师一样，开始计算出地下地质属性的分布情况。
这个过程可能会有点漫长，就像等待面包发酵一样，我们得有耐心。
在计算的过程中，计算机可是在飞速运转，做着各种复杂的计算，而我们就像在旁边焦急等待结果的小粉丝，心里既期待又有点小紧张呢。
地震资料反演及Jason软件学习共55页文档
31、园日涉以成趣，门虽设而常关。 32、鼓腹无所思。朝起暮归眠。 33、倾壶绝余沥，窥灶不见烟。
34、春秋满四泽，夏云多奇峰，秋月 扬明辉 ，冬岭 秀孤松 。 35、丈夫志四海，我愿不知老。
谢谢你的阅读
❖ 知识就是财富 ❖ 丰富你的人生
71、既然我已经踏上这条道路，那么，任何东西都不应妨碍我沿着这条路走下去。——康德 72、家庭成为快乐的种子在外也不致成为障碍物但在旅行之际却是夜间的伴侣。——西塞罗 73、坚持意志伟大的事业需要始终不渝的精神。——伏尔泰 74、路漫漫其修道远，吾将上下而
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
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主界面，主要填写内容如下：图 1Solid Model：点击Input，再点击Solid model，选择地质框架模型，如solidModel.smg Layers ：点击Input，再点击Layers, 选择相应目标层，如TopPad、Reservoir、BasePad, 利用→输入至右边的界面中,并且修改Spacing值,一般目标层厚度设为与需要分辨的薄层时间厚度值，而非目标层可以设置得稀一些，如与地震数据采样一致Trace Gate ：选择反演道门，如wholeProjectQC 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所示。
图5图5中显示了Reservoir 层中nopay 岩性中P-Impedance 的pdf 函数分布。
其它分析结果与之相似，为了直观地了解分析结果，用户应该仔细检查每一层每一岩性中各种属性的pdf 函数分析结果图件。
➢ *.png 文件：一系列*.png 文件, 实际是把Index.html 文件里的图件分开来显示, 一个*.png 文件就是一张图, 这些文件方便用户在做报告或文档时直接插入感兴趣的图件, 而不用再去抓图。
➢ scriptCopy.py 文件：就是根据在Generate Stats 界面中填写的参数产生的脚本文件，有点类似程序语言，但是相对简单多了。
用户应该逐渐熟悉脚本文件，先尝试阅读脚本，领会每部分的意思。
在逐渐熟悉后可以尝试修改脚本或自己编写。
足够熟练后甚至可以跨过填写界面这一步，直接编写脚本并运行，来完成地质统计学参数分析、模拟、反演等工作。
图6是该工区地质统计学参数分析脚本例子。
图6Reservoir 层中nopay 岩性中的 P-Impedance 的pdf 函数分析结果分析井段穿过Reservoir 层nopay 岩性的样本点统计结果准备信息开始逐层分析输入分析道门 输入地层框架模型 输入QC 道门输入分析层段以及设定相应的采样间隔 BasePad 层为非目标层，对其中的P-Impedance 不用分岩性分析 Reservoir 层为目标层，对其中的P-Impedance 应该分岩性进行pdf 函数的分析pay nopay➢ stats.py 文件：以上分析结果的具体数据都存放在该文件内，可以说是最终的结果文件，以后在做模拟与反演的过程中调用的就是该文件。
以Reservoir 层的nopay 岩性为例显示该文件的部分内容，如图7所示。
图7●● 变差函数的求取(1) JGW 主界面→ Modeling → Specials→ StatMod→ StatMod Analysis , 弹出StatMod Analysis 界面；(2) Input → Time/Depth mode: TimeSolid model ： SOLIDLayers ：点击Select layers 选择目标层位，如Reservoir 层，在Selectedinterval 设置采样间隔为0.001s ，点击Assign current interval to selected layers 应用修改的采样间隔Lithology masks ：先把Enable lithology 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, 弹出两个界面，一个是主界面StatMod Transform ，一个是辅界面Transform, 点击辅界面Transform 最上方的Histogram, 弹出一个小界面StatMod Histogram 并带着一个不停闪烁的红色提示Build the histogram , 点击该红色提示，则主界面出现相应的样本点分析，调试小界面中的Nr of intervals 的值, 尽量使pdf 函数样本点分布规则连续。