jason反演 第三讲 子波及其极性 100403
地震数据的准备--地震数据的极性认识
地震数据的极性判断： １、SEGY格式规定，初至波起跳向下，记录数值是负值，对应正反 射系数界面，此称“正常极性地震记录”。即波谷对应正反射系数 (Normal Format)； ２、地震子波是混合相位的； ３、不同的处理方法可以得到不同相位的子波，……地震数据的“视极 性”也就各异； ４、通常处理技术还很难解决子波的零相位化； ５、有一些资料处理模块是可以改变相位的。 （据李庆忠院士《走向精确勘探的道路》P134） 所以：子波是未知的！判断地震数据的极性有一定的困难。 平常我们拿到手的地震数据，其极性是什么？没有答案。因为野外 是保证极性下跳；处理中心数据入站，极性没有变化；但处理流程中， 每一个处理环节没有进行相位监控，所以到数据出站时，地震数据的相 位被改造成什么样？这是一笔糊涂帐。所以很难说清地震数据的极性是 什么，也就是子波相位是说不清楚的。这给我们的反演带来困惑。（解 决办法见合成记录）。
从上面两个例子来看，如果不注意子波的极性（即 相位），那么可能把阻抗计算反了，在九十年代初期， 曾发生过这种事情。也就是从那个时候起，才引起我们 的注意，可是目前还有许多解释员，只简单的用零相位 雷克子波做合成记录，然后用这个子波进行反演，这种 做法是不对的，有可能把高速与低速层做反了，那么地 震反演的结果就全错了，因此，子波的相位是至关重要 的！
《地震勘探数据处理规范》中关于极性的规定：
记录初至下跳（负值），称正常极性：反之，称为反正常极性。数据 处理中通常采用正常极性。 对于没有提供极性信息的测线，应放大显示少许单炮记录的初至，以 鉴定其极性。对反正常极性记录应进行反极性处理。在成果剖面中应 标注清楚记录道的极性。
地震子波的极性（SEGY格式规定：normal格式，正反射系数，负子波）
极性和层位标定
层位标定是储层研究工作最重要基础之一； 搞清基础资料的极性是正确标定的前提； 搞清子波相位是标定的基础。
有关极性的约定
地震勘探的正反射和负反射 --上覆为疏介质、下覆为密介质反射系数为正； --上覆为密介质、下覆为疏介质反射系数为负。
有关极性的约定
SEG极性标准：
1、对炸药震源，当爆炸压缩开始时初至下跳，产生一负数，图形显示下 跳（即初至相当正反射系数，处理剖面正常极性显示为波谷）； 2、对零相位子波，正反射系数由波峰中的规定
《地震勘探采集技术规范》中关于极性的规定：
年度开工前或每次改变大线、检波器和仪器输入电路时，必须采用炸 药激发，检查地震数据采集系统的极性（包括仪器、地震有关标准对 极性的规定电缆、检波器串），并将极性记录记带，磁带上的能量初 至应为负值。 经计算机显示，要求初至下跳，宽行打印初至样值是一负值。
子波时窗选择： 在估算子波时，应选择一个适当的时间范围（时 窗），地震数据中的子波在纵向上其振幅和相位都会 发生变化（子波是空变的），因此，如果子波参与反 演处理，则应选择短时窗（1000ms或更短）；时窗长 度应为子波长度的3-4倍。
子波长度:
子波应该有足够的长度，这样可以使子波的能量削弱到零，我们不希望在子波的边 缘处有尖锐的切除现象，因为可能有褶积人工痕迹，也可在子波谱中产生边界效应。频 率低的数据需要长的子波，一般情况下，我们更想让子波的长度越短越好，典型的长度 范围在60-128ms之间。 子波简单化：因为两个主要原因我们追求一个简单的子波： 1、处理的结果导致子波趋向简单化，除震源子波外，所有处理过程最少做过一次反褶 积处理，反复多次反褶积，我们希望简单些。 2、我们追求一个简单、可信的地质模型，该模型能用于油气勘探，在子波估算中，有 一些不确定因素，特别是第二个、三个旁瓣的存在，用户不喜欢在子波导致的第二个、 三个旁瓣获得的地质目标上钻探。 子波带宽：频率成分应与地震数据的频率完全一致，如果子波的频率超过高信噪比数据 的频率范围，则子波看起来很杂乱。 子波相位：因为所有因素，我们希望相位随频率有一个简单变化，这些扰乱相位的因素 （检波器、记录、处理）有平滑频率响应的能力，当子波太长时，相位就会反复无常， 与噪声相位响应相似，当子波比较长时，杂乱的相位就会在高频部分产生，再向低频部 分传播，因此，在有用的频宽内，不允许有反复无常相位出现。
下跳子波 0 空气 大地
反射系数
V1，ρ1 V2，ρ2
RC = ( ρ2v2 - ρ1v1) / (ρ2v2 + ρ1v1)
V1 < V2 ，ρ1 < ρ2
反射系数
正常极性子波（normal）
地震响应
V1 1
+
低速层
波谷
V2 2
波峰
高速层
V3 3 V1 < V2 > V3 ρ1 < ρ2 > ρ3
激发
接收
子波－由脉冲震源产生的信号（可有几个相位）。较厚的地层 的反射子波为零相位 ； 实际地震数据很难见到零相位子波，这是因为陆相沉积很难见 到产生零相位反射的又平又厚的常速地层 ； 地震数据是由不同形状子波和反射系数的“褶积”形式组成的；
子波形状的变化是由波阻抗的垂直变化或地质体边界的横向变 化产生的，每一种形状都和振幅变化有关 。
目前在子波估算中存在的问题： 1、简单地应用零相位雷克子波； 2、没有注意子波的相位（极性）； 3、空变子波要慎用! 以上因素都会影响地震反演的结果!
子波估算中应注意的问题：是子波估算,而不是子波提取。 1、子波的长度(保证一个主峰,两个旁瓣)； 2、估算子波的井旁道(最好是走向方向)； 3、估算子波的时窗(至少是子波长度的三倍以上，边界不要卡在强轴上)； 4、兼顾子波的波形及频谱(振幅谱,相位谱)，特别要注意在地震主频带内,相 位接近常相位。
几种判断子波极性的方法（也是判断地震数据的极性）：
a) 海上资料 ； b)如果工区地层有煤层沉积 ； c)如果工区地层有碳酸盐岩沉积； d)制作合成记录时，要兼顾浅、中深层的反射同相轴，然后用零相位子波 和负180度子波分别做两个合成记录，经过仔细对比；
e)需要通过多口井的合成记录，来判别子波的相位。
从地震数据中估算子波的三种方法：振幅谱子波、 相位谱子波和振幅谱相位谱子波。
三种子波估算方法说明
振幅谱子波：子波的振幅谱来自于子波自相关的付氏变换，这一计算是通过地 震自相关的付氏变换除以反射系数来实现的，反射系数来自声阻抗曲线；通过 子波长度和能量衰减来减少旁瓣，将估算的子波的自相关截短。
相位谱子波：一个不随频率变化的常相位输入子波，是通过输入的地震数据和 声阻抗曲线来估算的，振幅谱与输入的子波相同。子波估算过程由相关步骤完 成：首先，输入的反射系数从声阻抗曲线中导出，然后，最初始的子波的关键 参数，如相位、时移和比例因子等，是通过对反射系数和地震数据做最小平方 运算得到的，这些系数用于产生子波，接下来，在井周围的数据道上对子波进 行稀疏脉冲反演处理，得到的结果是一个改进了的声阻抗模型和一个改进的反 射系数估算，这一过程可以反复迭代，直到子波达到要求为止。 振幅和相位谱子波：输入的地震数据中的子波，其振幅谱和随频率变化的相位 谱被同时估算，提供一个优先的子波用于分析。在连续的处理步骤中，声阻抗 曲线用于导出反射系数，然后从这些反射系数和地震数据中得出一个最好的最 小平方子波，接下来用到White(1980)方法，然后用得到的子波做稀疏脉冲反 演，这一结果是一个新的改进了的声阻抗模型和一个改进的反射系数估算，第 一步可以按要求反复迭代，其它在子波估算中的关键输入参数是子波长度、子 波时移和数据中最大的有效频率，子波长度一定要足够长，与子波的能量相适 应。当所选的输入地震与井曲线没有优化结合时，子波开始时间就代表子波的 时移时间，依据Jason习惯，如果一个子波的开始时间与一半的子波长度的负 值相同的话，子波不做时移处理。
低速层
波谷对应正反射系数（+），代表地震波由低速层进入高速层；波峰对应负 反射系数（-），代表地震波由高速层进入低速层。
反射系数
反极性子波（reverse）
地震响应
V1 1
+ 低速层
波峰
V2 2
高速层
波谷
V3 3
V1 < V2 > V3 ρ1 < ρ2 > ρ3
低速层
波峰对应正反射系数(+)，代表地震波从低速层进入高速层;波谷对应负 反射系数(-)，代表地震波从高速层进入低速层。
影响地震子波相位（极性）的因素：
（1）大地吸收的影响； （2）野外激发接收的影响；
（3）不同年份施工、不同仪器施工、不同处理参数，也会造成子 波相位的混乱；
（4）有一些处理模块会改变子波的极性，但目前在处理中，没有 全程对子波相位进行监控，所以最后出站的结果，不知道子波相 位被改造成什么样； （5）目前的处理技术还很难解决子波的零相位化，这些因素都给 子波相位的估算带来极大的困难！
多井子波提取和优化
子波应用的目的
子波应用的目的:
1、利用多种技术估算地震子波:单独利用地震数据、 利用一口井或多口井反射、利用部分叠加或多个地震 数据部分叠加;
2、最优化的井震联合； 3、在横向非均匀介质中内插子波； 4、产生地质建模所用的合成子波；
5、计算零相位地震数据的反褶积算子；
6、估算Q值和产生基于Q值的子波。 使用一个合适的子波，是获得满意的油藏描述结果的关键！
子波的合理估算--子波的相位
不同的地震反演软件，对子波相位的设计是不同的： 1、Jason 软件中,对子波的振幅谱和相位谱是同时估算，软件通过子波的 相位谱会自动识别是波峰代表正反射系数，还是波谷代表正反射系数，不 需要人为干预。 2、Strata软件中，需要人工指定是波峰代表正反射系数,还是波谷代表正 反射系数。一般情况下，该软件缺省的极性设置是：波峰代表正反射系数； 因此需要首先调查地震数据的极性，把地震数据处理为接近零相位；然后 进行反演，以保证反演结果的正确性。 有人说：“进行地震反演，必须首先把地震数据乘以（-1），然后再进行 反演，否则结果就是错误的。”这个说法不确切。 1、首先要确认地震数据的极性，不是所有的地震数据都要乘以（-1）； 2、要看所应用的软件是如何处理相位这个问题的； 3、要认真学习有关的软件，了解其功能的具体设置，不要以偏概全。
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*、递归 等
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地震反演
关于波阻抗反演中的低频分量有两个基本的问题：（1）在波阻抗反演中，为什么要构造低频信息？地震资料里为什么没有低频信息？（2）为什么井资料有低频信息？如何用它构造用于波阻抗反演的低频信息大家在实际勘探开发过程中都在讲我们的地下储层薄，地震资料不能满足要求，讲提高分辨率的问题，楼主提了一个很有意思的问题，低频分量。
在这里我想尝试对此进行一个简单分析：1、什么是低频？学过信号处理的人都知道，我们可以把一个波（各种波场都可以）变换到频率域，在频率域频率较低的叫低频.如果用一个通俗的说法，低频就是“趋势”，如果你炒股的话，每天的成交量属于高频，而30日均线上看到的算低频，如果说到天气，下一个小时的气温算高频，但北京今年秋天的平均气温变化算低频。
对于我们搞地学的来讲，地下的压实规律是低频，某套地层的砂地比是多少算低频，某套地层的平均阻抗也算低频。
2、低频有什么用？通过上面的分析可以看出，低频反映的是某种事物或事物的某种特征的变化规律，是更本质的东西。
如果你把握住了股市的低频，你就可以低价买，高价卖。
如果你把握了房价的低频，你也可以发大财。
如果你把握了地下储层变化的低频，你就知道那里砂多，那里泥多，你的勘探成功率也会提高。
3、地震资料有没有低频？通过讲地震资料变换到频率域，你很容易知道，地震资料中没有低频，这主要是由于地震检波器的工作原理决定的。
对模拟检波器来说，<10Hz的低频很难收到。
对现在逐渐应用的数字检波器，小于3Hz的信号也很难被收到。
4、井资料中有没有低频？一般来说，大家在使用井资料的时候，主要用她的高频，看某个砂层的情况，内部岩性、物性、含油气性变化等，但如果测量井段较长，你就会看到低频：比如对渤海湾上第三系地层，馆陶组总体砂多泥少，而明化镇泥多砂少，具体到某口井，明下段又有什么规律，这就是低频。
5、反演中为什么需要低频？反演是我们把得到的地震信号反推回地下地质结构的过程，由于地震资料中缺失低频，自然就需要通过各种方法进行补偿。
Jason反演技术在建南地区储层预测中的应用
地边缘沉积了生物礁、 生屑滩 ， 低部位则沉积了硅泥 质， 其沉积相 出现 了明显 的分异 （ １ 。建南研究 图 ）
区处 于 开 阔台地 一 台地边 缘相 带 。
震体的约束稀疏脉冲反演， 得到的反演平面图中， 可以清晰刻画出建南构造南部长兴组有利储层的 展布范围， 其预测效果与钻井资料吻合较好 ， 对下步勘探部署有一定的指导作用。
关 键词 ．ａｏ 演技 术 ； ＂ ｓｎ反 Ｊ 约束稀 疏脉 冲反 演 ； 长兴组 ； 层预 测 储
１ 研究 区地质背景
南部进行波阻抗反演计算。通过对井的波阻抗反演剖面分析后认为， 阻抗反演质量较好 ， 波 与井的
符合程度较 高， 分辨率达到要求， 可以清楚反映长兴组储层 的横 向展布变化 ， 并通过未参加反 演计 算的 Ｊ３井检验， Ｘ 其反演结果与实钻情况相符 。根据试验线的成功反演 ， 进而在全区进行 了三维地
建南地 区位于湖北省利川市建南镇一带， 建南 构造处于川东褶皱带石柱复向斜中部。晚二叠世长
南
兴期 ， 自东南方 向的古太平洋和西南方 向的古特 来 提斯海水广泛入侵四川海域 ， 除西侧的康滇古陆有 较轻微的活动外 ， 北边 的马尔康、 摩天岭、 巴山等 大 古隆均没于水下 ， 使本区形成广阔的碳酸盐台地。
第２ １卷
第 ３期
江
汉
石
油
科
技
Ｖ０ ． Ｎｏ． １ ２ｌ ３
Ｓ ｐ． ０１ ｅ ２ １
பைடு நூலகம்
２１ ０ １年 ９月
Ｊ ＧＨＡＮ Ｐ Ｔ Ｌ ＵＭ Ｃ ＥＮＣ ＮＤ ’ ＨＮＯＬ Ｌ Ｅ Ｒ０ Ｅ ＳＩ ＥＡ Ｉ ＥＣ ＯＧＹ
JASON需要注意以下几点
需要注意以下几点
1. 对井数据进行进行地质统计学参数（概率密度函数和变差函数）分析之前，需要对井资料进行重新采样，采样间隔应该与最终地质统计学反演结果的采样率一致（如1ms、0.5ms）；
2. 在对重新采样井资料进行地质统计学参数（概率密度函数和变差函数）分析之前，要确定Tracegate，确保地震数据的线、道间距要一致（如25×25m），该Tracegate应与地质统计学反演时用的Tracegate一致；
3. 针对反演目的层段进行分层位、分岩性、分属性（Pimp、Por、Sw等）进行统计学参数（概率密度函数和变差函数）分析；
4. 在设置变差函数时，X和Y方向的变程一致，因为在MC算法中是假设XY是相等的，只有满足这一假设算法才会稳定、收敛；
5. 子波左右样本点必需对称，即要求子波零点左右的样点数一致；
6. 要在目地层上下设置边界层，边界层厚度要在半个子波长度以上，以消除反演时子波的旁瓣效应，
为了节约计算时间可以将上下界面的Layer的微层间隔设置大一些，但最大与地震的采样率一致；
7. 影响计算时间的几个因素：
1). 微层间隔，即layer的spacing，在完成solid model以后要检查模型的闭合情况，因为模型局部的变形会耗费大量的计算时间；
2). 地震道数的多少，数据量越大，时间越长，需要注意的是在做地质统计学反演时跳道和不跳道，以及垂向上采样率的不同，统计学参数（概率密度函数和变差函数）是不同的；
3). 地震信噪比S/N的大小，信噪比设置越大，运算时间越长；
4). 变差函数X和Y方向的变程的大小，变程越大，时间越长；
5). 参与反演井的数量不会对运算时间有太大的影响。
jason地震反演软件培训课件
jason地震反演软件培训
31
做好质量控制
jason地震反演软件培训
32
合理使用井的约束作用
反演过程中井的约束控制
F L p ( r ) L q ( s d ) 1 L1 ( Ztrend )
约束方法： 无 约 束 （ disable any constraint ） 软 约 束 （ soft constraint） ： α －1， 硬 约 束 （ hard constraint） ： 趋 势 ＋ 边 界 条 件
约束稀疏脉冲反演是基于道的反演，它的实质就是在阻抗趋势的约束 下，用最少数目的反射系数脉冲达到合成记录与地震道的最佳匹配。
jason地震反演软件培训
26
数据加载
分析和预处理、极性分析
井旁子波
相关性分析 校正
初标定 子波优选 综合标定
层位解释
断层解释及趋势面分析 地层结构分析 框架模型
初始模型 选择反演井、范围、子波等
4、反演结果特点
体现地震资料振幅、频率、相位特征， 适用范围较广 。但反演结果分辨率不太高。
（一） 稀疏脉冲地震波阻抗约束反演方法
Y8
Y118-
4
1
Y8
Y118-
4
1
10m reservoir
delineation
稀疏脉冲反演波阻抗剖面与对应的地震剖面对比
The Menu of Constrained sparse spike inversion
jason地震反演软件培训
7
5、近几年，从叠后反演发展到叠前反演 。 叠前同步反演能够反演出纵横波阻抗等，直接进行岩性 和流体识别。 Zoeppritz方程组是叠前反演的理论基础，但关系复杂 难以直接求解。Shuey和Aki、Richards等在二十世纪 八十年代对其进行了简化研究。Ali Tura（1999）、 David M Dolberg（2000）等AVO反演、AVA反演方 法。JASON公司等推出了实用的叠前和叠后地震反演 软件系统。
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反演第二讲 曲线标准化 100404
红3410井声波-密度曲线归一化前的分布形态(红色为HS4井曲线) 3410井声波-密度曲线归一化前的分布形态(红色为HS4井曲线) 井声波 HS4井曲线
1、地震反演所需资料准备（地震资料：叠后地震数据体、地震解释层位；井资料： 测井曲线、钻井分层数据）； 2、地震资料与井资料的数据加载； 3、使用直方图统计法进行曲线标准化处理； 4、子波估算与合成记录制做； 5、初始阻抗模型的建立（来自井信息的阻抗模型）； 6、约束稀疏脉冲反演的实现过程（相对波阻抗反演）； 7、反演地震道的合并（补充低频与高频成分）； 8、反演剖面上目的层的砂体识别（确定砂体的阻抗门槛值）与砂体平面预测（求取 砂体样点的时间累计厚度图）； 9、如何将反演结果输入到LandMark解释系统中？ 10、 LandMark解释系统中的砂体预测（自动求取砂体样点的时间累计厚度图及手工 在剖面上进行砂体追踪解释）。
9口井归一化处理后的声波曲线直方图统计（呈近似高斯正态分布） 口井归一化处理后的声波曲线直方图统计（呈近似高斯正态分布）
9口井归一化处理后的密度曲线直方图统计（呈近似高斯正态分布） 口井归一化处理后的密度曲线直方图统计（呈近似高斯正态分布）
Jason反演软件主界面 Jason反演软件主界面
红3417井声波-密度曲线归一化后的分布形态(红色为HS4井曲线) 3417井声波-密度曲线归一化后的分布形态(红色为HS4井曲线) 井声波 HS4井曲线
红74井声波-密度曲线归一化后的分布形态(红色为HS4井曲线) 74井声波-密度曲线归一化后的分布形态(红色为HS4井曲线) 井声波 HS4井曲线