稀疏脉冲反演方法及Jason软件实现
Jason稀疏脉冲波阻抗反演流程
地震反演是储层预测和油藏描述中非常关键的技术。
地震反演技术将地震数据、地质构造模型、井点测井资料有机地结合为一体；把地震资料的常规界面反射型剖面，反演成能与测井资料直接对比的岩层型波阻抗剖面,以及岩性剖面和孔隙度等地层物性剖面，反演结果分辨储层的能力也有了明显提高。
Jason稀疏脉冲波阻抗反演流程(附流程图)1 子波提取子波长度的选取要适宜，一般为80ms左右，在浅层，地震频带较宽，子波可短些；深层，地震频带窄，子波可略长些；对提取的子波，要从子波的波形、振幅谱、相位谱等方面进行子波质量的判断，要求子波的波形稳定、频域内单峰平滑、有效频带内相位稳定。
2 合成地震记录标定反演过程中，可以利用反演波阻抗数据，对地震地质标定进行精细微调，提高标定的质量。
子波提取与合成地震记录标定有着密切的关系，两者是循环反复进行的。
先用初始理论子波进行初步标定，标定的时深关系基本确定后，再用井旁地震道和测井反射系数重新估算子波。
完成最终的地震地质标定之后，再重新提取更精确合理的子波。
3地质层位解释解释的层位要尽量靠近反演目的层，要十分精细，反映出目的储层的空间变化特点来。
同时，使用时还要对层位中的缺失、奇异点进行修补、剔除，保证反演所用层位的正确性。
4 地震、地质、测井约束模型的建立建立尽可能接近实际地层沉积情况的波阻抗约束模型，是减少反演最终结果多解性的十分重要的一环。
建立波阻抗模型的过程实际上就是把地震界面信息与测井波阻抗信息正确结合起来的过程。
5 反演参数选择选择合理的反演参数是做好反演处理的前提，稀疏脉冲反演中的主要参数有稀疏约束项、拟合误差约束项以及模型约束项等等。
针对不同的地质情况，有目的性地选取合适的反演参数，同时利用相应的质量监控手段，加强对反演信息的反馈处理，可以有效地改进反演结果的质量，起到事半功倍的效果。
6反演结果检查处理过程中的反演数据体反复检查、各种资料对比分析，发现异常、寻找原因、反复反演，知道获得满意的结果为止。
jason
RockTrace同时AVA 约束稀疏脉冲反演自从2000 年秋季推出以来，RockTrace 就对行业应用和集成PSTM 地震数据的方法产生了重大影响。
它是唯一的定量集成测井曲线弹性岩石物性和AVA 地震数据的技术，可以生成标定的定量岩石物性三维数据体。
RockTrace 以InverTrace Plus技术为基础，并将该技术扩展到了AVO 域。
在InverTrace Plus中，应用的约束以波阻抗(Z p) 为依据。
在RockTrace 模块中，目标是在波阻抗之外求解出横波阻抗(Z s) 和密度，因此，对这三个参数分别地设置约束条件。
在并行处理过程中，能够生成的弹性参数类型组合为：∙纵波阻抗、横波阻抗和密度∙纵波阻抗、纵/横波速度比和密度∙纵波时差、横波时差和密度∙纵波速度、横波速度和密度和InverTrace Plus一样，应用全局模式时，一个空间控制项被加到目标函数上，同时，一个相当大的地震道数据体被整体转换。
RockTrace反演应用多个不同角度叠加道集的地震数据体，并且生成三个弹性参数数据体。
该算法是InverTrace Plus模块整体多道反演算法的扩展。
它是以一个地震数据体作为输入信息，并且只生成一个阻抗数据体（上述数据体之一）。
通常，RockTrace 算法可以产生三个弹性参数数据体，并且保留了许多其它模块的原理和约束算法。
这是业界一项独特的技术，具有以下技术优势：∙反演得到的弹性参数是岩石的真实属性，与储层属性有关。
∙当采用Knott-Zoeppritz 方程式选项时，无任何物理近似假设。
∙允许振幅和相位随偏移距变化。
通过对每隔输入部分叠加数据体，计算唯一的子波来实现反演。
∙在反演过程中，弹性参数可以直接进行各自约束。
∙岩石物理关系可用于约束弹性参数对。
∙由于所有输入数据必须和单一的输出模型相一致，降低了噪音的影响。
∙最终的弹性参数模型可重构输入地震数据，这也是反演优化算法的一部分。
Jason约束稀疏脉冲反演流程
12
JGW_8v2_workflow_SYB
【ASCII Horizon Import-(Untitled Session)】 ParametersExisting files，选择 Append and overwrite overlap Existing XY coordinates，选择 Do not read new XY coordinates。
在外部打开层位文件，查看内容。 第二列为 Inline，第四列为 Xline，第五列为 X 坐标，第六列为 Y 坐标，第七列为 Z 值。 【Edit ascii horizon generic colum format definition】 选中 Column 2， 点击 Edit， 弹出窗口 【Edit column definition】 ， Column definition 设置为 Line number。
11
JGW_8v2_workflow_SYB
【Import ASCII horizon files】 点击 Add，选择要加载的层位和断层文件，如果有深度域的层位，需要点击 Edit 将其设置 为 TVD 格式； 在 JGW horizon file 填写内部文件名称，如 demo。OK。
【ASCII Horizon Import-(Untitled Session)】 设置完成。
弹出窗口【Configure project “Jason_demo”】 选择米制单位：Metric units，勾选或不勾选 Add extended data type list，OK。
2
JGW_8v2_workflow_SYB
新建工区完成，主窗口变成【Jason-Jason_demo】
3
jason地震反演软件培训.pptx
2、技术关键
（1）硬约束和软约束条件建立。 （2）λ值等关键参数试验。 （3）频带补偿、拓宽。
3、实际资料使用特点
使用地震资料品质较好。可用测井速度 或地震处理速度进行约束。
反演方法分类
1、基于原始地震资料 叠前反演类、叠后反演类
叠后地震反演使用叠后地震资料，反演 的纵波阻抗能够在一定程度上反映储层的变 化规律，但不能给出横波波阻抗。
叠前地震反演使用叠前地震道集资料， 具有良好的保真性和多信息性。能同时反演 纵、横波阻抗等反映地层岩性、物性和流体 特征的多种参数。 2、基于数学算法
线性反演类、非线性反演类
反演方法分类
3、基于实现方式
基于反褶积的反演方法：包括基于地层反褶积的道积分、 递归、广义线性等反演方法；基于最大似然反褶积的稀疏脉冲 反演方法;基于最大后验概率准则算法反褶积的模拟退火等。
基于波动方程的反演：Born反散射。 基于随机过程的反演：随机反演、随机模拟等。 基于特征分析的反演：特征反演、神经网络反演等。 其它：混沌反演等
第二章 叠后约束地震反演原理及常用反演方法
第一节 稀疏脉冲地震波阻抗反演
第二节 随机模拟地震波阻抗反演 第三节 协模拟岩性及孔隙度反演 第四节 常用反演流程
Largo
RockTrace
EarthModel Wavelets
JASON Environment
InverTrace StatMod
InverModel
反演方法分类
4、基于地震和测井的相对作用
带限反演：道积分、Vlog、Seislog和Verilog*、递归 等
迭后约束稀疏脉冲反演中文培训教程
Data_loading 工区
该数据只是用于数据加载练习，在后面的反演和解释练习中并没有用到该数据。
Wedge2d 工区
该数据是合成的二维数据工区，该工区让您充分理解利用波阻抗数据来开展工作的好 处。
Big_gom 工区
该工区是更大范围的 GOM 工区，您将通过它学习如何解释 InverTracePlus 生成的反演 结果，在该练习中我们会用到工区内的 6 口井，在工区中有一个 InverTracePlus 反演数据体 aitm_ai.mod 文件。
第五章：建立低频模型............................................................................42 5.1 构造格架建立 ............................................................................42 5.2 低频模型的产生.........................................................................44
启动 Map View 模块
在 JGW 主菜单点击 File > Select project 并选择 GOM 工区，点击 Analysis > MAP VIEW 启动该模块。
数据的选择显示
通常有五种方法可以选择您需要加载的数据： －在 MAP VIEW 窗口顶部通过 Input 菜单。 －通过菜单下面的分类快捷按钮 , , , 进行数据加载。 －通过 Plot List，点击显示窗口右边的快捷按钮 可以进入，进行数据加载。
第二章：工区的建立与数据的加载..........................................................24 2.1 工区的建立 ................................................................................24 2.2 数据的加载 ................................................................................25
Jason 软件详细介绍与基本操作手册
Jason模块功能简单描述Jason地学平台是一个将各学科有用信息综合在一起，即为用户提供丰富的、从油气勘探、开发到生产阶段将地震、测井和地质信息有机结合，进行反演、模拟和预测不同阶段储层岩石物性和流体变化的综合研究工具。
您根据可用信息和经济条件，您可以明智的选择应用哪个“工作流程”更好。
1、Enveriment 运行环境（工区管理地震解释和接口）与其它解释系统紧密相连、数据格式灵活、方便。
2、EarthMod 地质框架模型，包括模型建造器，模型生成器和测井曲线生成器。
把地质、地震、测井曲线（斜井轨迹）、沉积模式（整合，不整合，礁体及河道等）建立以地震道为坐标的模型上，成为储层预测的基础，并赋予地质意义。
3、VelMod 速度模型建立工具VelMod 除具有建立速度模型的基本功能外，且有其独特特征：•地质引导：结合地质信息进行插值运算，能把稀疏的速度数据拓展为三维速度场。
•基于地质概念的速度调整•真三维插值4 、Wavelets 二维和三维地震子波分析●无井、单井、多井提子波和空变子波分析●从斜井中进行真三维分析●估算信噪比谱●为使结果稳定利用模型驱动5、 Modtrace 完整的地震反演系统( 包括对大倾角地层)●所用约束稀疏脉冲反演是基于道的递归反演，产生地震带宽内的波阻抗数据。
●其算法的唯一特征应用解释控制，为获得可靠的低频信息提供工具。
在EarthMod 的基础上建立空变约束条件，其结果突出各向异性，提高了反演分辨率。
●多种方式提取子波，包括无井、单井多道、多井多道、直井、斜井提子波和在EarthMod 的基础上提空变子波。
●由于大量地震信息的反演是基于EarthMod 的基础，所以适用与多条复杂断层的地质情况，使其反演结果更接近真实的地质模型。
●给出产层有效厚度图和孔隙度分布图６、InverMod 精细储层描述（多井）•采用一欧洲石油公司的专利，专门致力于薄层预测和精细描述。
对油田滚动勘探开发十分有用。
JASON软件介绍
JASON软件介绍RockTrace同时AVA 约束稀疏脉冲反演⾃从2000 年秋季推出以来，RockTrace 就对⾏业应⽤和集成PSTM 地震数据的⽅法产⽣了重⼤影响。
它是唯⼀的定量集成测井曲线弹性岩⽯物性和AVA 地震数据的技术，可以⽣成标定的定量岩⽯物性三维数据体。
RockTrace 以InverTrace Plus技术为基础，并将该技术扩展到了AVO 域。
在InverTrace Plus中，应⽤的约束以波阻抗(Z p) 为依据。
在RockTrace 模块中，⽬标是在波阻抗之外求解出横波阻抗(Z s) 和密度，因此，对这三个参数分别地设置约束条件。
在并⾏处理过程中，能够⽣成的弹性参数类型组合为：纵波阻抗、横波阻抗和密度纵波阻抗、纵/横波速度⽐和密度纵波时差、横波时差和密度纵波速度、横波速度和密度和InverTrace Plus⼀样，应⽤全局模式时，⼀个空间控制项被加到⽬标函数上，同时，⼀个相当⼤的地震道数据体被整体转换。
RockTrace反演应⽤多个不同⾓度叠加道集的地震数据体，并且⽣成三个弹性参数数据体。
该算法是InverTrace Plus模块整体多道反演算法的扩展。
它是以⼀个地震数据体作为输⼊信息，并且只⽣成⼀个阻抗数据体（上述数据体之⼀）。
通常，RockTrace 算法可以产⽣三个弹性参数数据体，并且保留了许多其它模块的原理和约束算法。
这是业界⼀项独特的技术，具有以下技术优势：反演得到的弹性参数是岩⽯的真实属性，与储层属性有关。
当采⽤Knott-Zoeppritz ⽅程式选项时，⽆任何物理近似假设。
允许振幅和相位随偏移距变化。
通过对每隔输⼊部分叠加数据体，计算唯⼀的⼦波来实现反演。
在反演过程中，弹性参数可以直接进⾏各⾃约束。
岩⽯物理关系可⽤于约束弹性参数对。
由于所有输⼊数据必须和单⼀的输出模型相⼀致，降低了噪⾳的影响。
最终的弹性参数模型可重构输⼊地震数据，这也是反演优化算法的⼀部分。
稀疏脉冲反演方法及Jason软件实现
第四篇 约束稀疏脉冲波阻抗反演1、概述：约束稀疏脉冲波阻抗反演方法（简称CSSI ）是目前应用比较广泛的一种方法，它是以趋势约束稀疏脉冲反褶积算法为基础。
对地下的反射吸收序列作如下假设：它的分布是稀疏分布的，而不是像其它的反演方法那样假设其是连续分布的。
在反演过程中，它是通过寻找一个使目标函数最小的反射系数脉冲数目来得到波阻抗数据的，此时合成的地震记录可以和实际地震记录匹配的最佳。
2、基本假设和方法原理约束稀疏脉冲反演（CSSI ）是一种递推地震波阻抗反演方法，基于范数反褶积和最大似然反褶积。
它假设地层的反射系数序列是由一系列稀疏且不连续的大反射系数和高斯分布的小反射系数叠加而成。
从地质角都解释，大反射系数代表地下不连续界面、岩性分界面。
首先由约束稀疏脉冲的基本假设，得到目标函数：)1ln()(2)ln(2)(1)(1122122λλ----+=∑∑==M L M K n NK r RJ LKLK（1）2R 为反射系数均方值，2N 为噪音均方值；第K 个采样点的反射系数为r(K),n(K)为噪音；M 为反射层数；L 为总采样数；λ为既定反射系数似然值；可经迭代计算出反射系数。
其次用最大似然反褶积反演得到反射系数。
采用递推算法结合初始波阻抗趋势，得到波阻抗模型。
)()(1)1()(i R i R i Z i Z +-= （2）Z(i)为第i 层的波阻抗值；R(i)为第i 层的反射系数。
最后约束稀疏脉冲反演依据目标函数逐一对地震道计算的初始波阻抗进行调整，并调整反射系数，得到目标函数：Z L a d s L r L F q p ∆+-+=-11)()(λ （3）r 是反射系数，Z ∆是阻抗趋势差；d 是地震道；s 是合成地震道；λ是残差权重因子；a 是趋势权重因子；p,q 是L 模因子。
式（3）右边分别是反射系数的绝对值和，合成声波记录与原始地震数据的差值，趋势约束项。
约束稀疏脉冲反演基于地震数据，用波阻抗趋势进行约束，以期用最少的反射系数脉冲，对合成记录与地震道进行最佳匹配。
0-叠后反演-Jason-讲
四、Jason6.2平台反演模块的构成
Invermod
Invertrace
Velmod Environment Plus s
Statmod
（地震反演）
Functionmod
• Recursive trace inversion (RTI 递归反演) • Constrained Sparse Spike Inversion (CSSI 约束稀疏脉冲反 演) • Trace merging (AITM 道合并) • Netpay and porosity estimation (NPPE 有效厚度和孔隙度估 算) 依据绝对阻抗与储集岩和孔隙度的关系计算有效厚度和孔隙度。
三、反演方法分类
1、基于原始地震资料 叠前反演类、叠后反演类 2、基于数学算法 线性反演类、非线性反演类 3、基于实现方式 基于反褶积的反演方法：包括基于地层反褶积的道积分、 递归、广义线性等反演方法； 基于最大似然反褶积的稀疏脉冲反演方法; 基于最大后验概率准则算法反褶积的模拟退火等。 基于波动方程的反演：Born反散射。 基于随机过程的反演：随机反演、随机模拟等。 基于特征分析的反演：特征反演、神经网络反演等。 其它：混沌反演等
• Principal component analysis （主因子分析） • Model estimation （模型估算） • Model generator （模型生成器）
四、Jason6.2平台反演模块的构成
Invermod
Rocktrace
Velmod Environment Plus Earthmodel
Velmod
Environment Plus Earthmodel
Statmod
Functionmod
Jason反演软件适用性分析
Jason反演软件适用性分析一）Jason软件适用性分析1. 概述波阻抗反演是目前应用最广泛的储层预测和油藏描述技术之一，其发展也日新月异，出现了各种各样的反演计算方法，这些方法由于其基于的数学物理模型的不同，各有优缺点，要求的基础资料和能够达到的反演效果也是不同的。
地质院于1999年首次引进Jason测井约束反演软件（部分模块），版本为5.1；2002年对该软件进行升级并扩充新模块，版本为6.0。
该软件于2002年11月进行了软件安装，并随之举行了为期两周的现场培训。
随后该软件正式投入生产，目前运行情况良好。
Jason软件自引进以来，在济阳坳陷隐蔽油藏勘探中发挥了非常重要的作用，目前已成为国内应用最普遍的波阻抗反演软件之一。
2. 软件功能和技术特点Jason软件是一个综合性、多方法的平台，它可以实现由递归反演－稀疏脉冲反演－随机模拟、随机反演－地震地质特征反演－弹性反演的处理过程，用户可根据需要由Jason 的模块构建自己的研究流程，完成储层预测和油藏描述。
其重要特点是随着越来越多的非地震信息（测井，测试，地质）的引入，由地震数据推演的油气藏参数模型的分辨率和细节会得到不断的改善。
主要包括Environment Plus、Wavelets、Velmod、Invertrace、Invermod、Rocktrace、Statmod、Largo和Functionmod等9个模块。
1)、vironment PlusEnvironment-Plus 集总了Jason的运行环境、各种分析工具和各种辅助性模块。
实现数据输入、输出（地震Segy格式，测井和层位ASCII格式）；各种数据，各种方式的显示(井、层位、地震等的2D/3D显示)；合成记录标定、井曲线编辑（方便、灵活）；绘图/直方图分析（确定属性间的关系）；三维（地质/储集）体自动解释（用门槛值/多边形定义解释标准；解释结果（层顶、底，层厚度，层间属性）可保存；报告解释体的叠置及井钻遇等情况；沿层属性提取（沿层面，层间等提取属性最大/最小/平均值等）；层位类数据的各种处理计算（平滑，加/减等）；地层异常检测（用属性的梯度异常确定地层/异常体的边界）；处理工具包（重采样，滤波，互相关等）；等值线（构造、属性等等值线）。
稀疏脉冲反演原理
稀疏脉冲反演原理稀疏脉冲反演原理是一种在信号处理领域中常用的技术，它可以通过对传感器接收到的信号进行处理，从而获取到目标物体的信息。
本文将介绍稀疏脉冲反演原理的基本概念、原理和应用。
稀疏脉冲反演原理是一种通过发送稀疏脉冲信号并对接收到的信号进行反演处理来获得目标信息的方法。
在传感器接收到信号后，通过对信号进行处理，可以对目标物体的位置、形状、大小等信息进行推测和分析。
这种方法在地震勘探、医学成像、雷达探测等领域有着广泛的应用。
稀疏脉冲反演原理的基本原理是利用信号的稀疏性来提取目标信息。
稀疏性是指信号在某个域中具有较少的非零元素，这意味着信号在这个域中可以通过较少的参数进行表示。
通过对信号进行稀疏变换，可以将信号从时间域转换到稀疏域，然后利用稀疏性对信号进行反演处理，从而获得目标的信息。
稀疏脉冲反演原理的应用非常广泛。
在地震勘探中，可以通过发送地震波信号并对接收到的反射波信号进行反演处理，从而获取到地下地层的信息。
在医学成像中，可以通过发送超声波信号并对接收到的回波信号进行反演处理，从而获得人体组织的信息。
在雷达探测中，可以通过发送雷达波信号并对接收到的回波信号进行反演处理，从而获取到目标的位置和形状等信息。
稀疏脉冲反演原理的优点是可以在较短的时间内获得目标信息，并且对信号的要求较低。
它可以通过对信号进行稀疏变换和反演处理，从而提取出目标的关键信息，避免了对整个信号进行复杂处理的过程。
同时，稀疏脉冲反演原理还可以通过优化算法对稀疏性进行进一步提取，从而提高目标信息的准确性和可靠性。
然而，稀疏脉冲反演原理也存在一些局限性。
首先，稀疏脉冲反演需要对信号进行稀疏变换和反演处理，这些处理过程需要消耗大量的计算资源和时间。
其次，稀疏脉冲反演原理对信号的要求较高，需要信号具有较好的稀疏性才能实现较好的反演效果。
最后，稀疏脉冲反演原理对噪声和干扰的抗干扰能力较弱，容易受到外界环境的影响。
稀疏脉冲反演原理是一种通过对信号进行稀疏变换和反演处理来获取目标信息的方法。
jason软件培训-zap
1、软件简介引进时间：2001.3 美国APEX 公司软件概况：6.0(h) 10大模块拥有6大模块开发程度：6大模块主要功能开发并应用推广程度：4大模块应用××××一、软件的技术原理和方法2、反演原理概念：利用地震资料,恢复地下地质结构和岩石性质的方法从有限频带宽度的地震数据中恢复出宽带波阻抗。
原理:反演的基础:褶积模型假设前提:1）地震模型2）反射系数序列3）地震子波4）噪音分量3、反演方法①递归反演：基于反射系数的计算得到波阻抗Z i+1=Z i((1+r i)/(1-r i))优点1、完全使用地震资料进行波阻抗反演2、计算速度快缺点：1、反射系数的带限严重、低频损失严重2、误差累积问题严重3、反演方法②稀疏脉冲反演：它假设地震反射系数是一系列大的反射系数叠加在高斯分布的小反射系数，大的反射系数相当于主要的不整合面或岩性界面，其目的是寻找一个使目标函数最小的脉冲数目，然后得到波阻抗数据。
（dj-sj）q+α2∑(tj-zj)2OBJF=∑│rj│p+λq∑j优点：1、使用地震资料进行波阻抗反演2、产生一个地质特征的结果3、反演结果中包含了低频信息缺点：1、子波选取不当，准确度降低2、只能达到地震资料的分辨率③基于模型的反演：在地质解释的基础上，同时结合测井资料构造反演的初始模型，采用模型优选迭代算法，不断修改更新模型，使模型正演合成地震资料与实际地震资料数据最佳吻合，最终的模型数据便是反演结果优点：1、分辨率、反演效果明显提高2、对薄互层反演有明显的优势缺点：1、初始模型要准确2、反演结果具有多解性④随机模拟地质统计学是以区域化变量理论为基础，以变差（或变异）函数为基本工具来研究那些展布于空间并且呈现出一定结构性和随机性的自然现象。
储层参数的空间分布问题正好符合这一自然现象。
而随机反演技术是随机模拟与地震反演相结合的产物，是将随机模拟的思想引入到反演中，用地震做约束，用随机模拟算法去实现储层预测。
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。
反演及Jason软件介绍
通过这一部分，用户可以建立自己的工区， 并向工区内加载钻井、测井、地震、地质 数据，数据的加载方式有两种，一种是其 他应用软件（ProMax LandMak GeoQuest 等）直接生成的文件，通过Jason相应的接 口程序可直接加载到Jason应用软件中；而 另一种则是通过Geoshare数据接口将其他 应用软件（如IES/IESX,Geoframe等）中有 关数据传到Jason数据库中。
基于模型反演方法思路
Wavelet
Seismic Data
பைடு நூலகம்
Initial Model
Best fit No
Synthietic Section
Yes
Model Perturbation
Final Model
测井约束反演
测井约束反演：是一种基于模型的波阻 抗反演技术，这种方法利用测井资料， 以地震解释的层位为控制，从井点出发 进行外推内插形成初始波阻抗模型，然 后利用共轭梯度法，对初始波阻抗模型 不断进行更新，使得模型的合成记录最 佳逼近于实际地震纪录，此时的波阻抗 模型便是反演结果。
预处理 输入地震、测井、地质等资料
建立地质模型（EarthMod）
子波估算 进行主成分分析PCA 进行模型估算Model Estination 生成反演数据体 Model Generator (Al、KTH、Rt、Sp…..) 用Function Mod生成相应的平面图
InverMod反演流程图
七、地质统计与随机模拟系统
叠前反演 地震反演 叠后反演
旅行时 振幅 旅行时 振幅 递推法
CT成像 AVO成像 构造反演 波阻抗反演
直接反演
道积分 无井
波阻抗反演 广义线性反演
JASON软件模块操作说明
JASON软件操作使用说明编写人：钮学民、徐美娥慎国强、孙振涛刘会芹特殊处理所JASON (5.0) 一．FileSelect project …Create project…Copy project …Exit二．ModelingLargo …正演及井校正EarthModel …地质框架模型Wavelets …子波VelMod …速度模型InverTrace …地震道反演ModTrace …RockTrace …弹性反演InverMod …储层精细描述StatMod …随机反演FunctionMod …Specials三．AnalysisVolume view and body checking …Section view and interpretation …Map view and calculation …Well log editing and seismic tie …Cross plots and histograms …Attribute extraction …Auto snapper …Graph view …Well log view …Stratigraphic anomaly checkerProcessing toolkit …Cross correlation …四．DatalinksLandmarkGeoQuestGeoshare half_link …Seismic/property dataHorizonsWellsWaveletsVelocity data …TablesCaltural data …Export for upscaling …Specials五．UtilitiesProject managementDowngrade (64 to 32 bit) …Downgrade (5.0 to 4.x) …Well upgrade (4.x to 5.0) …Project copy …Project parameters …Project coordinates …Project description …File manager …Horizon managementMerge data files …Inquire data files …Downscale utility …Create/Edit Tables …Data attributes editor …File Manager …Specials→JGW←→ASCII 六．SystemSubmit a software change request …Show system load/disk usage …Unix terminal window …Inspect current license features …Show batch queue …Show postscript printer queue …Screen snapping tools …Screendump …XV …Hints and tips …Netscale navigator一、File菜单：Select project:选择工区，必须是已经建好的工区，或重新恢复的工区。
稀疏脉冲波阻抗反演参数对反演效果的影响研究_穆星
2 反演参数
根据目标函数 , 可以看出控制反演的参数有 p 、q 、λ(目标函数控制参数), 再加上 反演的迭代 次数参数以及子波的波形参数等 。 2. 1 子波的波形参数
关键词 :储层预测 ;波阻抗反演 ;稀疏脉冲反演 ;反 演参数
中图分类号 :P631. 5
文献标识码 :A
收稿日期 :2004 — 10 — 18
EFFECT OF PARAMETERS CONTROL ON SPARSE SPIKE IMPEDANCE INVERSION
M U Xing
( The College of the Earth’ s Resources and I n formation , University of Petroleum , Dongy ing 257000, China)
根据最大似然反褶积计算得到的反射系数 ,
结合初始阻抗模型 , 采用递推算法 , 反演得到初始
的波阻抗模型 :
Z(i)=Z(i
-
1)11
+R(i) - R(i)
(2)
式中 , Z(i)为第 i 层的波阻抗值 ;R(i)为第 i 层的
反射系数 。
3)结合井的约束条件进行波阻抗反演
约束稀疏脉冲反演对每一道依据目标函数对
稀疏脉冲波阻抗反演参数对反演效果 的影响研究
穆 星
(石油大学(华东)信息与资源学院 , 东营 257000)
摘 要 :稀疏脉冲波阻抗反演是目前最普遍应用的波阻抗反演方法之一 。 该方法在计算过程中涉及各 种类
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约束稀疏脉冲反演在苏里格气田下石盒子组砂岩储层预测中的应用
Jason约束稀疏脉冲反演在苏里格气田下石盒子组砂岩储层预测中的应用X刘盼喜,时义刚(中国石化中原油田分公司物探研究院,河南濮阳 457001) 摘 要:苏里格气田为鄂尔多斯盆地岩性油气藏,充分利用高分辨率地震资料,在经过测井曲线标准化后,在储层精细标定基础上进行约束稀疏脉冲反演(C SSI)对下盒子组砂岩进行储层预测,优选滚动探井,钻遇32.8m/L 储层并获高产气流,取得了滚动勘探的成功,为同类岩性油气藏的勘探开发提供了参考和借鉴。
关键词:苏里格气田;岩性气藏;滚动勘探;约束稀疏脉冲反演 中图分类号:P 631.4 文献标识码:A 文章编号:1006—7981(2012)01—0135—02 苏里格气田位于长庆靖边气田西北侧的苏里格庙地区。
鄂尔多斯盆地陕北斜坡的西北侧.目前已发现上古生界二叠系下石盒子组盒7、盒8、山西组山1、山2及下古生界奥陶系马家沟组马五5、马五6等多套含气层段。
苏里格气田是一个低孔、低渗、低压、低丰度、大面积分布的岩性气藏(图1)。
图1 研究区位置图1 资料解释1.1 合成记录及层位标定层位标定是储层横向预测的关键。
它将地震资料与钻井资料相互沟通,使二者之间建立一个准确的对应关系。
为后续的地震地质层位追踪解释、储层预测等提供可靠目标。
利用Welledit 模块作合成记录,针对目的层段从波形、频率、相位等特征与地震资料对比,对合成记录进行微小的压缩、拉伸,使其与地震相关性最佳(图2)。
1.2 层位解释利用Discovery 解释系统,以储层标定为基础,对目的层,进行对比追踪精细解释,获得所需的地震层位数据,使其作为目标反演中保持地震数据能够在同一低频背景下运算的约束条件,从而提高初始模型精度。
分别编制了盒8底与山西组底的深度构造图(图3)。
2 约束稀疏脉冲波阻抗反演(C SSI)2.1 测井响应特征盒8砂岩储层主要电性特征(表1)为“三低一大”即低自然伽玛、低密度、低声波时差;大幅度自然电位异常。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
第四篇 约束稀疏脉冲波阻抗反演
1、概述：
约束稀疏脉冲波阻抗反演方法（简称CSSI ）是目前应用比较广泛的一种方法，它是以趋势约束稀疏脉冲反褶积算法为基础。
对地下的反射吸收序列作如下假设：它的分布是稀疏分布的，而不是像其它的反演方法那样假设其是连续分布的。
在反演过程中，它是通过寻找一个使目标函数最小的反射系数脉冲数目来得到波阻抗数据的，此时合成的地震记录可以和实际地震记录匹配的最佳。
2、基本假设和方法原理
约束稀疏脉冲反演（CSSI ）是一种递推地震波阻抗反演方法，基于范数反褶积和最大似然反褶积。
它假设
地层的反射系数序列是由一系列稀疏且不连续的大反射系数和高斯分布的小反射系数叠加而成。
从地质角都解释，大反射系数代表地下不连续界面、岩性分界面。
首先由约束稀疏脉冲的基本假设，得到目标函数：
)1ln()(2)ln(2)(1
)(1
1
22
1
22
λλ----+=
∑∑==M L M K n N
K r R
J L
K
L
K
（1）
2R 为反射系数均方值，2N 为噪音均方值；第K 个采样点的反射系数为r(K),n(K)为噪
音；M 为反射层数；L 为总采样数；λ为既定反射系数似然值；可经迭代计算出反射系数。
其次用最大似然反褶积反演得到反射系数。
采用递推算法结合初始波阻抗趋势，得到波阻抗模型。
)
()
(1)
1()(i R i R i Z i Z +-= （2）
Z(i)为第i 层的波阻抗值；R(i)为第i 层的反射系数。
最后约束稀疏脉冲反演依据目标函数逐一对地震道计算的初始波阻抗进行调整，
并调整反射系数，得到目标函数：
Z L a d s L r L F q p ∆+-+=-11)()(λ （3）
r 是反射系数，Z ∆是阻抗趋势差；d 是地震道；s 是合成地震道；λ是残差权
重因子；a 是趋势权重因子；p,q 是L 模因子。
式（3）右边分别是反射系数的绝对值和，合
成声波记录与原始地震数据的差值，趋势约束项。
约束稀疏脉冲反演基于地震数据，用波阻抗趋势进行约束，以期用最少的反射系数脉冲，对合成记录与地震道进行最佳匹配。
3 、CSSI反演流程
先对地震资料细致处理并提取子波；对地震层位精细解释，建立合理地质地层模型；关键是对测井资料做标准化分析，结合地震子波做合成记录，以保证时间深度域准确对应，结合地层框架模型和合成记录，以三维叠后地震数据为基础，利用地震子波以及初始模型趋势模型，进行波阻抗反演，最终得到精细反演结果。
反演流程如下：
4、二维模型测试
选用二维mamous模型，利用Jason软件进行稀疏脉冲反演，目标油气层大约在0.35s处。
其基本流程如下：
图4-1：原始地震道
: 图4-2：解释层位
图4-3：井震标定：相关系数为0.82
图4-4：建立低频模型
图4-5：反演结果
5、实际资料反演
选用某工区的实际叠后资料来进行叠后的稀疏约束脉冲反演，本文采用的是Jason软件来进行反演。
其基本流程如下：
（1）确定反演的范围，解释层位：
5-1：原始地震剖面及解释的层位
5-2：原始地震的连井剖面
5-2：原始地震剖面在频率域的显示（2）对反演范围内的井进行井震标定：
5-2：对工区CB32井标定，相关系数0.41，
5-3：对工区CB804井标定，相关系数0.76 （3）构建低频模型
5-4：固体模型
5-5：低频模型
（4）进行反演
调节好各个参数（特别是稀疏约束因子、地震信噪比、合并频率、子波刻度的选择）进行反演，结果如下：
5-6波阻抗反演剖面
:
5-7：连井波阻抗反演剖面。