STRATA地震反演培训-理论和练习.ppt

《STRATA地震反演培训-理论和练习.ppt》由会员分享，可在线阅读，更多相关《STRATA地震反演培训-理论和练习.ppt（286页珍藏版）》请在金锄头文库上搜索。

1、STRATA 地震反演培训地震反演培训理论和练习理论和练习2STRATA 培训STRATA 课程大纲课程大纲地震反演基础地震反演基础- 3 - 3 练习练习1 1：楔型模型简单练习（递推、模型）：楔型模型简单练习（递推、模型）- 6161练习练习2 2：Erskine 3DErskine 3D 数据加载数据加载- 99- 99测井曲线时深对比-118练习3：Erskine 3D 测井曲线时深对比-131测井曲线层间插值建模-156子波提取子波提取-174-174练习练习4 4：BlackfootBlackfoot 子波提取子波提取-198-198基于模型反演参数基于模型反演参数-225-225

2、练习练习5 5：BlackfootBlackfoot 基于模型反演基于模型反演-241-241其它反演参数其它反演参数-256-256练习练习6 6：BlackfootBlackfoot 其它反演方法其它反演方法-271-2713STRATA 培训地震反演基础地震反演基础反演反演：从地震数据中提取它所包含的潜在地质信息的过程.传统上，反演是在叠后地震数据上进行,目的是提取声波阻抗体声波阻抗体.近年来，反演已经被扩展到叠前数据体,目的是既提取声波阻抗又声波阻抗又提取横波阻抗体提取横波阻抗体，从而可以推得孔隙流体。另外，最新的发展认为可以利用反演结果直接预测岩性参数如孔隙度和含水饱和度。4STRA

3、TA 培训几种不同反演方法几种不同反演方法递归法递归法: 传统带限反演传统带限反演颜色反演颜色反演: 递归反演的改进递归反演的改进稀疏脉冲法稀疏脉冲法(两个不同算法两个不同算法): 约束以获得尽可能少的同相轴约束以获得尽可能少的同相轴基于模型法基于模型法: 迭代更新层状初始模型迭代更新层状初始模型弹性波阻抗弹性波阻抗: 对对AVO数据的改进数据的改进LMR: 对对AVO数据的改进数据的改进联合反演：联合反演： 对对AVO数据的改进数据的改进在在STRATA软件中提供以下这些反演方法软件中提供以下这些反演方法(六种方法七个模块六种方法七个模块):5STRATA 培训 阻阻 抗抗纵波纵波横波横波弹

4、性波弹性波反射系数反射系数子波子波地震道地震道所有反演方法中共同的正演模型所有反演方法中共同的正演模型:反演中的一般正演模型反演中的一般正演模型6STRATA 培训波阻抗波阻抗反射系数反射系数反演反演子波子波地震道地震道反演试图反推得到这个正演模型反演试图反推得到这个正演模型:反演模型反演模型纵波、横波、弹性阻抗纵波、横波、弹性阻抗7STRATA 培训纵波阻抗纵波阻抗或或横波阻抗横波阻抗或或 弹性阻抗弹性阻抗阻抗阻抗 反射系数反射系数纵纵 波波 阻阻 抗抗 =横横 波波 阻阻 抗抗 =弹弹 性性 阻阻 抗抗 = 复杂公式复杂公式 (后面介绍后面介绍)8STRATA 培训反射系数反射系数 地震道

5、地震道地震道地震道 = 子波与反射系数的褶积子波与反射系数的褶积+噪音噪音.注释注释(1)不存在模型化的多次波不存在模型化的多次波.(2)没考虑传播损失和几何扩散没考虑传播损失和几何扩散.(3)没考虑频率吸收没考虑频率吸收.(4)子波可能是时变的子波可能是时变的.9STRATA 培训反射系数与子波褶积反射系数与子波褶积的结果是移除了大量的结果是移除了大量高频细节高频细节:10STRATA 培训在时间域中的褶积就是频率在时间域中的褶积就是频率域中的乘积域中的乘积.从这些图中可以看出从这些图中可以看出,子波的子波的作用是将地震频谱中高频和作用是将地震频谱中高频和低频都消除了低频都消除了.理论上讲理

6、论上讲,反演就是试图将这反演就是试图将这些失去的频率区域进行恢复些失去的频率区域进行恢复.11STRATA 培训反演中反演中“多解性多解性”问题问题所有反演算法都有多解性问题所有反演算法都有多解性问题.存在多于一种地质模型可以与地震数据存在多于一种地质模型可以与地震数据相一致，要在这些可能的模型中决定一相一致，要在这些可能的模型中决定一个模型，就需要去选择地震数据以外的个模型，就需要去选择地震数据以外的的一些信息。的一些信息。通常用以下两种方法来使用地震以外的通常用以下两种方法来使用地震以外的信息：信息： 初始猜测模型的建立初始猜测模型的建立 对最终结果与初始猜测模型背离对最终结果与初始猜测模

7、型背离幅度进行约束幅度进行约束最后反演结果既依赖于最后反演结果既依赖于 “其它信息其它信息”也取也取决于地震数据决于地震数据.12STRATA 培训带限带限 (递归递归) 反演反演递归反演递归反演, 也称为带限反演是最简单和最早一种反演方法也称为带限反演是最简单和最早一种反演方法.从第一层开始从第一层开始, 可以用这个公式递推后续可以用这个公式递推后续每层的阻抗每层的阻抗. :ith +1 层的阻抗可以从上一层层的阻抗可以从上一层 ith 确定确定:从定义反射系数开始从定义反射系数开始:13STRATA 培训在这个简单的实例中可知在这个简单的实例中可知:(a) 告诉我们如果给出一个脉冲，告诉我

8、们如果给出一个脉冲，那么我们就可以恢复阻抗的真实值那么我们就可以恢复阻抗的真实值, 但是但是(b) 中看到如果我们用一个子波与中看到如果我们用一个子波与脉冲进行褶积，那就无法恢复阻抗脉冲进行褶积，那就无法恢复阻抗的真实值。的真实值。Z2 = 818Z2 = 1500Z1 = 1000Z3 = 1227Z4 = 1004Z1 = 1000带限反演带限反演14STRATA 培训输入地震数据输入地震数据递归反演数据递归反演数据递归反演递归反演获得获得的结果与输入的结果与输入的地震数据频的地震数据频带一样带一样.从图中可以看从图中可以看出出,与测井比较与测井比较反演结果损失反演结果损失了高频细节了高频

9、细节.15STRATA 培训步骤步骤1:递归反演递归反演 的初始背景模型是通过对井阻抗滤的初始背景模型是通过对井阻抗滤波获得波获得:10-Hz 高截滤波高截滤波16STRATA 培训步骤步骤 2:对地震道运用递归运算对地震道运用递归运算. (注释注释: 这几乎等同于在相位域进这几乎等同于在相位域进行行-90 度转换度转换):17STRATA 培训步骤步骤 3:滤波后的模型加上滤波后的模型加上比例化比例化反演道集就得到了最终结果反演道集就得到了最终结果:+=18STRATA 培训递归反演中几个要点递归反演中几个要点:(1)不考虑子波不考虑子波. 这意味输入的地震数据必须是零相位这意味输入的地震数

10、据必须是零相位. 如果有如果有提取子波可行提取子波可行Strata 软件可以自动地对地震数据软件可以自动地对地震数据 “去相位去相位”.(2)即使地震数据是零相位即使地震数据是零相位, 实际子波的旁瓣实际子波的旁瓣在递归算法中也会在递归算法中也会被解释为岩性的变化被解释为岩性的变化.(3)反演的结果其频带与地震数据反演的结果其频带与地震数据频带一样频带一样.(4)对对地震数据的加权值地震数据的加权值(与反射系数匹配与反射系数匹配)决定了递归反演的结决定了递归反演的结果阻抗值范围是否合适果阻抗值范围是否合适.19STRATA 培训基于模型的反演基于模型的反演基于模型反演基于模型反演 是从褶积模型

11、方程式开始的是从褶积模型方程式开始的:假定地震道假定地震道S,和子波和子波 W, 是已知的是已知的.假定噪音是随机的并与地震信号不相关假定噪音是随机的并与地震信号不相关.求解反射系数求解反射系数, R, 以满足这个等式以满足这个等式. 这是一个非线性问题这是一个非线性问题,所所以求解过程是迭代进行的以求解过程是迭代进行的.20STRATA 培训步骤步骤 1:基于模型反演基于模型反演 的初始背景模型是通过对井阻抗分段形的初始背景模型是通过对井阻抗分段形成的成的:用户用毫秒用户用毫秒(MS)来定义层分块大小来定义层分块大小.所有的层分段开始时都被设成一样所有的层分段开始时都被设成一样大小大小(用毫

12、秒用毫秒).21STRATA 培训步骤步骤 2:利用层状模型和利用层状模型和已知已知的子波计算人工合成地震道的子波计算人工合成地震道.这是与实际地震道的比较这是与实际地震道的比较.通过分析人工合成地震记录通过分析人工合成地震记录与实际地震道的误差，每层与实际地震道的误差，每层(分段分段)的厚度和震幅值都进的厚度和震幅值都进行修改以减少误差行修改以减少误差.这个步骤通过一系列递归不这个步骤通过一系列递归不断重复断重复.SyntheticSeismic22STRATA 培训输入地震数据输入地震数据基于模型反演基于模型反演基于模型反演基于模型反演 获得了宽频的结获得了宽频的结果果.问题是高频成分问题

13、是高频成分可能来自于初始可能来自于初始猜测模型而不是猜测模型而不是来自地震数据来自地震数据.23STRATA 培训递归反演递归反演基于模型反演基于模型反演这是递归反演和基这是递归反演和基于模型反演的比较于模型反演的比较.通常基于模型反演通常基于模型反演可以获得更细致信可以获得更细致信息息,但这个结果实际但这个结果实际上相当类似上相当类似.24STRATA 培训基于模型反演的几个要点基于模型反演的几个要点:(1)由于已经知道子波由于已经知道子波,在计算过程中它的影响从地震数据中排除了在计算过程中它的影响从地震数据中排除了.也就是说也就是说地震数据没有必要是零相位地震数据没有必要是零相位,而只要而

14、只要子波与地震数据等相位子波与地震数据等相位就可以就可以.(2)估算的估算的子波若有误差子波若有误差,将导致反演结果出差错将导致反演结果出差错.(3)地震有效地震有效分辨率分辨率得到提高得到提高.(4)反演结果可能在很大程度上依赖原始猜测初始模型反演结果可能在很大程度上依赖原始猜测初始模型.解决的方法是对初始解决的方法是对初始模型模型进行滤波进行滤波.(5)与其它反演方法一样与其它反演方法一样,存在存在多解性多解性问题问题.25STRATA 培训稀疏脉冲反演稀疏脉冲反演稀疏脉冲反演稀疏脉冲反演 假定实际假定实际反射可以认为是由一系列反射可以认为是由一系列大脉冲里夹杂有小脉冲背大脉冲里夹杂有小脉

15、冲背景景.稀疏脉冲反演稀疏脉冲反演假定假定 只有只有大脉冲有意义大脉冲有意义.该方法通该方法通过检查地震道来寻找大脉过检查地震道来寻找大脉冲的位置冲的位置.26STRATA 培训稀疏脉冲反演稀疏脉冲反演 每次建每次建立一个脉冲反射序列，立一个脉冲反射序列， 通过增加脉冲直到地震通过增加脉冲直到地震道被足够准确地反演完道被足够准确地反演完成。成。阻抗块的振幅值是由阻抗块的振幅值是由基基于模型反演于模型反演算法来确定算法来确定27STRATA 培训输入地震输入地震稀疏脉冲反演稀疏脉冲反演稀疏脉冲反演稀疏脉冲反演生成一个宽带生成一个宽带高频反演结果高频反演结果.28STRATA 培训基于模型反演基于

16、模型反演稀疏脉冲反演稀疏脉冲反演稀疏脉冲反演稀疏脉冲反演结果与结果与基于模基于模型反演型反演结果类结果类似似.主要区别是缺主要区别是缺少非常薄层的少非常薄层的细节细节.29STRATA 培训稀疏脉冲反演几个要点稀疏脉冲反演几个要点:(1)只有当地震数据只有当地震数据有脉冲有脉冲，稀疏脉冲反演才得到同相轴；，稀疏脉冲反演才得到同相轴；(2)它试图利用地震数据获得它试图利用地震数据获得最简单的可能模型最简单的可能模型；(3)通常得到的结果比地质本身实际的通常得到的结果比地质本身实际的同相轴少同相轴少；(4)它比它比基于模型反演基于模型反演更更少的依赖于初始猜测模型少的依赖于初始猜测模型。30STR

17、ATA 培训颜色反演颜色反演颜色反演颜色反演 是对是对递归反演递归反演的修改的修改, 最早是由最早是由BP公司的公司的 Lancaster and Whitcombe 在在2000 年年SEG年会上提出来的年会上提出来的.这种反演过程这种反演过程, 寻找一个简单操作因子寻找一个简单操作因子 O, 对地震道进行直接转换求得反演对地震道进行直接转换求得反演结果结果:作者在频率域确定简单操作因子作者在频率域确定简单操作因子, O . 通过比较地震数据和实际反演结果，认为简单操作因子的相位是通过比较地震数据和实际反演结果，认为简单操作因子的相位是-90度。度。31STRATA 培训简单操作因子震幅谱用

18、这种方式来求取简单操作因子震幅谱用这种方式来求取:利用工区的一组井利用工区的一组井, 所有所有井的声阻抗的振幅谱做井的声阻抗的振幅谱做交汇图交汇图.通过理论预测通过理论预测,我们可以我们可以拟合一条直线来拟合一条直线来 代表代表“理想理想” 输出阻抗谱输出阻抗谱.Log(Frequency)Log(Impedance)声阻抗震幅谱声阻抗震幅谱32STRATA 培训然后然后, 利用井旁一组地震道利用井旁一组地震道, 求求取取平均地震谱平均地震谱.地震谱地震谱Frequency (Hz)操作因子谱操作因子谱从以上两个频谱从以上两个频谱,求取操作因子谱求取操作因子谱. 使得地震谱形状在地震频带内转换

19、使得地震谱形状在地震频带内转换为为阻抗谱阻抗谱.33STRATA 培训颜色反演操作因子颜色反演操作因子Time (ms)将获得的震幅谱做将获得的震幅谱做-90 度相移就生成了度相移就生成了颜色反颜色反演操作因子演操作因子.将它与所有地震道褶积将它与所有地震道褶积.34STRATA 培训输入地震输入地震颜色反演颜色反演颜色反演颜色反演 结果非常类似结果非常类似于递归反演于递归反演.其中一个区别是颜色反其中一个区别是颜色反演的结果是相对声阻抗演的结果是相对声阻抗,既有正值也有负值既有正值也有负值.-30000+300035STRATA 培训递归反演递归反演颜色反演颜色反演-30000+300046

20、00830012000相对相对 AI绝对绝对 AI36STRATA 培训颜色反演几点小结颜色反演几点小结(1)除了定义整体阻抗趋势外除了定义整体阻抗趋势外,与初始模型基本无关与初始模型基本无关.(2)运算非常快运算非常快.(3)用户所需要定义的参数非常少用户所需要定义的参数非常少.(4)假定地震数据是零相位的假定地震数据是零相位的.(5)在最初实施中，该方法可以获得在最初实施中，该方法可以获得相对相对阻抗体，虽然我们可阻抗体，虽然我们可以选择添加低频阻抗背景趋势。以选择添加低频阻抗背景趋势。37STRATA 培训AVO 反演反演基本褶积模型反演假定的是零偏移距地震数据基本褶积模型反演假定的是零

21、偏移距地震数据.褶积反演法不能运用到带有褶积反演法不能运用到带有AVO响应的地震数据响应的地震数据, 因为褶积反演法没有明因为褶积反演法没有明确地考虑确地考虑 VP/VS 的变化的变化.将反演方法延伸到叠前将反演方法延伸到叠前AVO 地震数据的处理地震数据的处理, 目前这些算法主要运用于目前这些算法主要运用于:(1)弹性阻抗反演弹性阻抗反演(2)Lambda-Mu-Rho (LMR)反演反演(3)联合反演联合反演38STRATA 培训弹性阻抗弹性阻抗弹性阻抗概念最初由弹性阻抗概念最初由Connolly 先生提出先生提出(The Leading Edge, 18, no. 4, 438-452

22、(1999).他利用他利用Aki-Richards 方程方程,将反射震幅与入射角联系起来将反射震幅与入射角联系起来:注释注释:在常规反演理论中假定在常规反演理论中假定 B=C=0, 而没而没有考虑有考虑VP/VS 的变化的变化.39STRATA 培训注意注意, 对于零偏移距对于零偏移距:通过类推通过类推, Connolly 定义了一种新阻抗类型定义了一种新阻抗类型:通过数学处理通过数学处理,最终最终:40STRATA 培训这张图重叠显示了同一口井位置弹性阻抗与声波阻抗这张图重叠显示了同一口井位置弹性阻抗与声波阻抗.在碳在碳氢聚集区域弹性阻抗值显示异常低值氢聚集区域弹性阻抗值显示异常低值.41S

23、TRATA 培训道集AVO 分析近角度叠加 stack at q1远角度叠加 stack at q2反演成弹性阻抗 EI(q1)反演成弹性阻抗 EI(q2)这种反演的工作流程是从迭前道集开始，首先这种反演的工作流程是从迭前道集开始，首先产生两组角度叠加剖面产生两组角度叠加剖面, 分别对其进行反演分别对其进行反演.42STRATA 培训这就获得两种反演结果这就获得两种反演结果:远角度反演远角度反演近角度反演近角度反演43STRATA 培训近角度反演与远角度反近角度反演与远角度反演进行交绘显示演进行交绘显示.远角度反演异常低值区远角度反演异常低值区可被解释为潜在含有碳可被解释为潜在含有碳氢聚合物氢

24、聚合物.近角度反演近角度反演远角度反演远角度反演44STRATA 培训将交绘图中的异常低将交绘图中的异常低值点从交绘图中成像值点从交绘图中成像到原始地震数据剖面到原始地震数据剖面上以进行解释上以进行解释.45STRATA 培训Lambda-Mu-Rho (LMR)LMR方法运用了方法运用了VP, VS, r r之间的关系式及拉梅常数之间的关系式及拉梅常数l l和和 m m:注意：最终结果表达式是注意：最终结果表达式是 lrlr 和和mrmr ，可以通过纵波、，可以通过纵波、横波阻抗横波阻抗 ZP和和ZS的关系式的关系式来表示。来表示。LMR 方法最初是由方法最初是由 Goodway et al

25、 (SEG Expanded Abstracts, 1997)提出来的。提出来的。类似于弹性阻抗法，这种方法对处理含类似于弹性阻抗法，这种方法对处理含AVO效应的数据比常规反演更有利。效应的数据比常规反演更有利。46STRATA 培训GathersAVO AnalysisRP Estimate RS EstimateCross-plotInvert to ZPInvert to ZSTransform to lr and mrLMR反演流程首先从叠前数据分反演流程首先从叠前数据分析得到纵横波反射系数体析得到纵横波反射系数体RP 和和 RS对对RP 和和 RS分别反演会生成纵横波分别反演会生成纵

26、横波阻抗体阻抗体ZP 和和 ZS .这些数据体会通过类似这些数据体会通过类似Goodway等关系式或交汇图进行转换。等关系式或交汇图进行转换。LMR 反演流程反演流程47STRATA 培训如图展示了对如图展示了对Albert地区地区典型含气砂岩运用典型含气砂岩运用LMR法法得到的结果，图中椭圆内得到的结果，图中椭圆内指示为含气砂岩。指示为含气砂岩。 上剖面为上剖面为体体, 下剖面为下剖面为体体 LMR 反演例子反演例子48STRATA 培训左图为左图为与与交汇图，红色区域指示气层交汇图，红色区域指示气层 (低低) 蓝色区域指示泥岩和水砂岩。蓝色区域指示泥岩和水砂岩。这些区域投到剖面上如下图所示

27、，可以指示含这些区域投到剖面上如下图所示，可以指示含气砂岩区域（红色部分）气砂岩区域（红色部分）lambda-rhomu-rhoLMR 反演实例反演实例49STRATA 培训联合反演联合反演联合反演联合反演可以对输入的叠前角道集数据反演得到可以对输入的叠前角道集数据反演得到ZP, ZS, 或者可能密度体。这种方法或者可能密度体。这种方法的优点在于它可以将这些参数进行约束，这样就可以使结果更稳定并且尽量避免多的优点在于它可以将这些参数进行约束，这样就可以使结果更稳定并且尽量避免多解性问题。解性问题。where:我们从 Fattis version of the Aki-Richards 方程开始

28、。这个公式用角度作为变量来表示反射系数：50STRATA 培训联合反演联合反演Ln(Zs)Ln(Zp)Ln()Ln(Zp)假定水砂岩可以用一个线性关系式来表示。假定水砂岩可以用一个线性关系式来表示。联合反演联合反演寻找背景趋势的偏差寻找背景趋势的偏差 :51STRATA 培训联合反演联合反演联合反演联合反演 生成生成Zp、Zs、Density体体,以及推导得到以及推导得到的各种组合。的各种组合。52STRATA 培训联合反演联合反演联合反演的解释与其它联合反演的解释与其它AVO反演结果类似：反演结果类似：53STRATA 培训它是地质统计反演的一种形式，可以通过对输入地震数据体运算产生大量的反

29、它是地质统计反演的一种形式，可以通过对输入地震数据体运算产生大量的反演结果以此来解决多解性的问题。每一个运算结果要与地震数据一致，并且保演结果以此来解决多解性的问题。每一个运算结果要与地震数据一致，并且保持如方差图中所包含的所期望的连续条件。通过分析这些结果以估算得到一个持如方差图中所包含的所期望的连续条件。通过分析这些结果以估算得到一个非确定性结果，以及最大可能结果。非确定性结果，以及最大可能结果。随机反演随机反演Stochastic InversionThis is a form of geostatistical inversion which explicitly addresses

30、the non-uniqueness problem by producing a large range of inversion results for a given input seismic volume. Each of the results is consistent with the seismic data, and honors the expected continuity conditions, as contained in the variograms. These results are analyzed to give an estimate of the u

31、ncertainty in the result, along with the most probable result.54STRATA 培训一般地震反演流程一般地震反演流程(1) 建模型建模型选择测井选择测井每口井进行相关每口井进行相关提取子波提取子波读取读取 / 拾取地震反演拾取地震反演(2) 完成反演完成反演选择反演方法和参数选择反演方法和参数对反演结果进行对反演结果进行QC(3) 解释反演结果解释反演结果做数据切片做数据切片做交绘图做交绘图输出到输出到EMERGE 软件软件55STRATA 培训反演结果反演结果QC两种方法两种方法:(1)误差图误差图(2)交互验证交互验证我们如何知

32、道反演结果是我们如何知道反演结果是否有效否有效?输入地震输入地震反演结果反演结果56STRATA 培训从阻抗道利用已知的子波从阻抗道利用已知的子波可以计算人工合成道可以计算人工合成道.理想情况理想情况, 这个人工合成这个人工合成道应该与输入地震道非常道应该与输入地震道非常相似相似.输入地震输入地震人工合成道人工合成道57STRATA 培训从输入地震道中减去人工合从输入地震道中减去人工合成道成道,就得到就得到反演误差图反演误差图.如果反演效果不错如果反演效果不错,误差图上误差图上应该只有非常小的震幅应该只有非常小的震幅,并没并没有聚集在某处有聚集在某处.由于存在多解性，即使微小由于存在多解性，即

33、使微小的误差图也不能保证反演结的误差图也不能保证反演结果就是正确的果就是正确的.输入地震输入地震反演误差反演误差58STRATA 培训第二种质量控制方法是第二种质量控制方法是交互验证交互验证. 在这个验证中在这个验证中,我们在初始我们在初始模型中删除一口井模型中删除一口井,在其位置做反演在其位置做反演,然后将反演结果与隐藏然后将反演结果与隐藏井进行比较井进行比较.隐藏井隐藏井反演结果反演结果反演误差反演误差59STRATA 培训每口井分析误差每口井分析误差,以便识别有问题的井以便识别有问题的井.60STRATA 培训反演作为反演作为 EMERGE一种属性一种属性最新使用反演结果是将其做为最新使

34、用反演结果是将其做为EMERGE软件的软件的输入输入, EMERGE可以直接预测孔隙度和其它岩性可以直接预测孔隙度和其它岩性体体.反演反演EMERGE孔隙度体孔隙度体61STRATA 培训练习练习 1: 楔形模型练习楔形模型练习第一个练习对楔型模型做反演第一个练习对楔型模型做反演. 目的是在简单例子上学习反演的基本步骤目的是在简单例子上学习反演的基本步骤.在在WINDOWS(UNIX)上键入上键入geoview或或从从Windows启动启动GEOVIEW 软件软件.我们将创建一个新的我们将创建一个新的GEOVIEW 数据库数据库. 选择选择New, 如显示如显示, 然后在显示的第一个然后在显示

35、的第一个菜单上点击菜单上点击 OK（如上图）（如上图）点击点击Well Explorer出现如下图所示界面。出现如下图所示界面。62STRATA 培训现在我们为这个工区创建一个新的数据库现在我们为这个工区创建一个新的数据库如果如果Geoview工区已打开可工区已打开可以点以点 Database / New创建创建如果这是如果这是GEOVIEW第一次在第一次在这台机器上运行这台机器上运行 ，你可以看，你可以看到如右图窗口，选择到如右图窗口，选择New并并点击点击OK.63STRATA 培训将数据库命名将数据库命名 wedge_database Ok:现在，可以看到现在，可以看到 Well Exp

36、lorer 中没有井显示：中没有井显示：64STRATA 培训现在，为这个数据库添加井曲线，如下步骤：现在，为这个数据库添加井曲线，如下步骤：点点 Import Data:选择选择 Wedge_well_logs.las 文文件并点件并点Next65STRATA 培训在接下来的菜单中，选择默认的即可：在接下来的菜单中，选择默认的即可：Note that we are creating a new well called Wedge_well:最后，我们可以读取声波曲线等如下：最后，我们可以读取声波曲线等如下：在这个菜单在这个菜单, 选择默认的参数即可，如选择默认的参数即可，如XY坐标：坐标：6

37、6STRATA 培训接下来点接下来点OK，出现如下提示窗口，问你选择哪种单位，出现如下提示窗口，问你选择哪种单位， 点击点击 Ok 即可：即可：67STRATA 培训当曲线加载完成后，出现如下窗口：当曲线加载完成后，出现如下窗口：点点 “1” 选择这口井选择这口井, 然后点然后点Display well. 曲线就会曲线就会显示出来。显示出来。68STRATA 培训接下来启动接下来启动Strata，选择，选择 GEOVIEW 主窗口上的主窗口上的 STRATA 按钮按钮.选择选择 Start New Project:给新工区命名：给新工区命名： “wedge_project”: 69STRATA

38、 培训在在STRATA 窗口上窗口上,通过通过 Data Manager / Import Data / Open Seismic / From SEG-Y File 加载加载地震数据地震数据:在在File Selection 菜单菜单, 选择文件选择文件 “wedge.sgy”. 点击点击 Next :70STRATA 培训这菜单有几页来定义文件格式和数据的几何形状这菜单有几页来定义文件格式和数据的几何形状(观测系统观测系统). 下一页下一页,尽管数据是一条尽管数据是一条2D测线测线,我们可以建立我们可以建立3D观测系统观测系统.接下来接下来, 需要告诉程序在道头中没有需要告诉程序在道头中没

39、有Inline & Xline 数数 也没有也没有 X & Y坐标坐标.71STRATA 培训点击点击Next 两次两次. 地地震数据扫描后震数据扫描后,出现这出现这样界面样界面:电击电击Ok 利用这些缺省利用这些缺省值值. 后面的练习将更后面的练习将更仔细地观测分析地震仔细地观测分析地震数据的加载数据的加载.72STRATA 培训现在可以看见这个楔型模型。现在可以看见这个楔型模型。另外，还会有如下图所示窗口另外，还会有如下图所示窗口显示出来，定义井在地震上的显示出来，定义井在地震上的位置，可以在位置，可以在Xline处填处填45，然后点击然后点击OK即可。即可。73STRATA 培训这口这口

40、井就会显示在这个模型上井就会显示在这个模型上 然而，然而，井曲线并不一定对应准确，所以需要井曲线并不一定对应准确，所以需要做校正，点做校正，点Well / Edit/Correlate Well:让你选择井来进行校正，选择这口让你选择井来进行校正，选择这口井，然后点井，然后点Edit.74STRATA 培训eLog 窗口就会出现，窗口就会出现，右侧显示的是模型右侧显示的是模型数据数据75STRATA 培训eLog 软件窗口弹出软件窗口弹出(显示井曲线显示井曲线 wedge log).完成井曲线对比完成井曲线对比,在井口位置求一个平均在井口位置求一个平均地震复合道地震复合道. 缺省值是在井口缺省

41、值是在井口+/-1周围道周围道求平均地震道求平均地震道.点击点击Ok 使用缺省值使用缺省值.进行井曲线对比，点击进行井曲线对比，点击Correlate 按钮按钮:76STRATA 培训eLog 软件窗口显示如下软件窗口显示如下: 蓝色显示人工合成记录蓝色显示人工合成记录道道. 红色显示复合地震道红色显示复合地震道. 原始地震道用黑色显示原始地震道用黑色显示.在对井曲线进行对比前需在对井曲线进行对比前需要一个更准确的子波要一个更准确的子波.77STRATA 培训用于计算合成记录的子波是缺省的称为用于计算合成记录的子波是缺省的称为“current wavelet”. 查看缺省子波查看缺省子波, 点

42、击点击 Wavelet / Display Current Wavelet:缺省子波是带通滤波缺省子波是带通滤波(频率频率: 5/10 - 50/60. ) 点击点击Frequency 标签可在频标签可在频率域显示子波率域显示子波.78STRATA 培训求取求取 “统计统计” 子波（不使用井信息），生成零相位子波子波（不使用井信息），生成零相位子波,其其频率域与地震数据的响应一致，点击频率域与地震数据的响应一致，点击 Wavelet / Extract Wavelet / Statistical:在子波提取的界面在子波提取的界面,第一页是设置分第一页是设置分析窗口析窗口. 使用缺省值使用缺省值

43、. 点击点击Next 两次两次.79STRATA 培训使用缺省值使用缺省值. 这些参数后面再讨论这些参数后面再讨论.点击点击OK 提取统计子波提取统计子波.新子波如左图所示新子波如左图所示. 检查它的频率响应检查它的频率响应,点点击击Frequency 标签标签.80STRATA 培训在窗口的底部子波名字也变了在窗口的底部子波名字也变了.井对比的过程是将蓝色的合成记录道上的井对比的过程是将蓝色的合成记录道上的点对应于红色地震复合道上的点点对应于红色地震复合道上的点. 用鼠标用鼠标左键点中，将四个同相轴连接左键点中，将四个同相轴连接.如右图所示点击如右图所示点击Stretch 按钮按钮:新子波自

44、动成为新子波自动成为 “current wavelet”.这意味它被用于在井对这意味它被用于在井对比窗口中合成记录的制作比窗口中合成记录的制作. 注意到蓝色合成记录道有微小的注意到蓝色合成记录道有微小的变化变化.81STRATA 培训点击点击 Ok (在在Check Shot Analysis窗口窗口).井对比类似于井对比类似于check shot 校正校正. 故出现两个新窗口故出现两个新窗口, 允许设允许设置一些参数置一些参数. 后面更仔细谈论这些参数后面更仔细谈论这些参数.82STRATA 培训点击点击Wavelet / Extract Wavelet / Use Well:选择选择 we

45、dge_well 点击点击Next 使用缺省使用缺省值点击值点击 OK (在最后一页在最后一页).这时你会发现合成记录拉伸后这时你会发现合成记录拉伸后 与地震数据拟合的更好与地震数据拟合的更好. 现在可以利用井与地震之间的相关求取更好的子波现在可以利用井与地震之间的相关求取更好的子波.83STRATA 培训接下来都选择默认即可接下来都选择默认即可这一页定义合成道在井这一页定义合成道在井位置处进行合成位置处进行合成选择时窗：选择时窗：84STRATA 培训这一页定义子波提取方法这一页定义子波提取方法以及基本的子波参数：以及基本的子波参数：接下来的课程讨论一下这接下来的课程讨论一下这些参数些参数8

46、5STRATA 培训新子波看起来与上一个很相象新子波看起来与上一个很相象:但是点击但是点击Frequency 标签显示有标签显示有一个非常小的相位校正一个非常小的相位校正.86STRATA 培训求取了子波，点击求取了子波，点击Ok 查看新查看新的对比结果的对比结果.菜单告知新对比后的声波曲线将用名字菜单告知新对比后的声波曲线将用名字P-wave_corr. 点击点击Ok :点击点击File / Exit Window ( (在在在在 eLog 窗口窗口):87STRATA 培训现在加载地震层位用于初始猜测现在加载地震层位用于初始猜测模型的内插。模型的内插。在在Horizon Selection

47、 菜单菜单菜单菜单, 使用缺使用缺省值名字点击省值名字点击Ok. 点击点击No (在在 dialog 窗口询问是否显示层位的平窗口询问是否显示层位的平面图面图).点击点击 Horizon / Pick Horizons:88STRATA 培训在在 STRATA 窗口底部显示层位拾取的一些选项窗口底部显示层位拾取的一些选项:89STRATA 培训如下图所示：如下图所示：注意：一个选择参数注意：一个选择参数“Mode”，我们选择，我们选择 “Rubber Band”. 如果我们选择这项后，必须将如果我们选择这项后，必须将剖面上目的层点亮，点下左键，剖面上目的层点亮，点下左键，释放鼠标后，释放鼠标后

48、， 层位会被拾取出层位会被拾取出来。来。也可以选择也可以选择 “Right and Left Repeat”. 我们选择用这种方我们选择用这种方式。式。90STRATA 培训拾取第一层位拾取第一层位, 置鼠标在波置鼠标在波峰峰150 ms周围周围:点击左键点击左键. 同相轴自动左右同相轴自动左右拾取拾取:T开始拾取第二层位开始拾取第二层位, 点击点击 (在在 Horizon 下拉菜单下拉菜单).设置拾取模式设置拾取模式“Left & Right Repeat”.91STRATA 培训使用缺省名字使用缺省名字, Horizon 2. 然后改变然后改变 抓抓获获Snap 选项为波谷选项为波谷 Tr

49、ough:点击楔型顶部的波谷点击楔型顶部的波谷.用类似方式用类似方式, 拾取楔型底部拾取楔型底部以及在以及在450 ms 周围的波谷波谷如图示如图示.完成后点击完成后点击Ok存储拾取的存储拾取的层位层位.92STRATA 培训现在可以建立初始猜测模型现在可以建立初始猜测模型, 利用井曲线和四个拾取的层位利用井曲线和四个拾取的层位. 点击点击 Model / Build/Rebuild a Model (在在STRATA窗口窗口):第一页设置新模型的名字第一页设置新模型的名字. 利用缺利用缺省参数省参数,参数含义后面课程再详细参数含义后面课程再详细解释解释. 连续点击连续点击Next 最后点击最

50、后点击OK 建立模型建立模型.93STRATA 培训初始猜测模型生成，。初始猜测模型生成，。在在wedge 地震数据地震数据后面初始猜测阻抗模后面初始猜测阻抗模型用彩色显示型用彩色显示.94STRATA 培训完成第一种反演完成第一种反演, 点击点击 Invert / Bandlimited:在反演菜单上在反演菜单上, 将输出名字将输出名字Output Volume File改为改为 “bandlimited ,其它用缺省其它用缺省值值:95STRATA 培训点击点击OK 得到反演结果得到反演结果:96STRATA 培训现在进行第二种反演现在进行第二种反演(在初始模型窗口在初始模型窗口). 点击

51、点击Invert / Model Based:改变输出名字改变输出名字 Output Volume File 为为为为“model_based”:使用其它缺省获得基于模型的使用其它缺省获得基于模型的反演结果反演结果Model Based:97STRATA 培训鉴别反演结果好不好的一个重鉴别反演结果好不好的一个重要标准是要标准是 “error display”,点击点击眼睛可以从反演窗口查看此项眼睛可以从反演窗口查看此项,如右图如右图在参数菜单选择在参数菜单选择“model_based derived Synthetic Error”, 如右图所示如右图所示,点击点击Ok.如右图如右图,可以看到

52、原始楔形模型与可以看到原始楔形模型与反演结果生成合成记录之间的差反演结果生成合成记录之间的差异异, 这个值接近于这个值接近于0,指示反演结果指示反演结果是成功的。是成功的。98STRATA 培训在在Open Window Data Menu 显示了刚生成的所有数据体显示了刚生成的所有数据体:该练习结束该练习结束. 点击点击File / Exit Project (在在 STRATA 软件任何窗口软件任何窗口)点击点击Yes (在警告诊断窗口在警告诊断窗口)使使STRATA 关关闭所有窗口并保存工区闭所有窗口并保存工区.99STRATA 培训练习练习 2, 对加拿大西部的一块对加拿大西部的一块3

53、D数据做反演数据做反演. 第一步是读一口井曲线到已经定义好第一步是读一口井曲线到已经定义好的的GEOVIEW 数据库中数据库中.在做这个之前，先创建一个数据库在做这个之前，先创建一个数据库点击点击 Database / New.定义这个数据库名称定义这个数据库名称 erskine_database.练习练习 2: Erskine 3D 初始化初始化100STRATA 培训在在 GEOVIEW Well Explorer 窗口窗口, 点击点击Import Data / Logs, Check Shots, Tops, Deviated Geometry from Files:101STRATA

54、培训改变改变 Destination Well Name 域为域为 “erskine_well”点击点击Next :在在 File Import 页面页面, 选择文件选择文件e“erskine_log.las” 点击点击Next .102STRATA 培训显示有一口井声波曲线显示有一口井声波曲线. 点击点击Ok 读取这根曲线读取这根曲线.点击点击Next (利用缺省位置信息利用缺省位置信息)选取默认单位选取默认单位103STRATA 培训在井曲线读进在井曲线读进GEOVIEW后后, 点击点击erskine_well 然后点击然后点击Display Well:erskine 声波显示如右图声波显

55、示如右图.104STRATA 培训在在GEOVIEW 主窗口上点击主窗口上点击STRATA启动启动 STRATA 软软件件.定义工区为定义工区为 “erskine_project”选择选择 Start New Project:在在STRATA 建立新工区对建立新工区对erskine 地震数据进行反演地震数据进行反演.105STRATA 培训第一步将第一步将3D 地震数据读进地震数据读进STRATA软件软件.点击点击Data Manager / Import Data / Open Seismic / From SEG-Y File:选择文件选择文件 “ersk3d.sgy”点击点击Next :

56、106STRATA 培训同样做为三维地震数据体加载地震数据同样做为三维地震数据体加载地震数据.这个文件在道头中没有这个文件在道头中没有Inline & Xline 号或号或 X & Y 坐标坐标. 设置设置 primary key to “CDP (Rectangular)”. 这是指数据体是由一系列测线组成这是指数据体是由一系列测线组成,每条线有相同的每条线有相同的道数道数.107STRATA 培训改变观测系统键入改变观测系统键入 “155” 如图显示如图显示. INLINE线将由程序计算线将由程序计算.点击点击Ok 加载地震数据体加载地震数据体.点击点击Next 两次选择缺省直到最后一页两

57、次选择缺省直到最后一页. 软件最初假定只有软件最初假定只有单一一条单一一条INLINE线线.108STRATA 培训在在 Well To Seismic Map 菜单菜单, 如下所示将井设置在地震如下所示将井设置在地震数据上数据上. 在在Plot box 中中wedge_well 井没有被选用因为它井没有被选用因为它不属于这个数据体不属于这个数据体:在地震显示中在地震显示中,键入键入 “24” 作为作为Inline 号显示的号显示的结果如右图所示结果如右图所示:109STRATA 培训对这个数据体要建立初始猜测模型对这个数据体要建立初始猜测模型,需要一组层位文件需要一组层位文件.首首先拾取一个

58、层位然后读进以前拾取的一组层位先拾取一个层位然后读进以前拾取的一组层位.点击点击Horizon / Pick Horizons:使用缺省名字使用缺省名字 “Horizon 1”.点击点击Ok.点击点击“Yes” (诊断窗口询问是否诊断窗口询问是否显示层位平面图显示层位平面图).110STRATA 培训向以前一样向以前一样, 将鼠标点在靠近层位处拾取单个层位将鼠标点在靠近层位处拾取单个层位:平面图窗口显示这条平面图窗口显示这条INLINE线上拾取的时线上拾取的时间间.111STRATA 培训软件利用这个人工拾取层作为引导可以自动拾取整个软件利用这个人工拾取层作为引导可以自动拾取整个3D数数据体层

59、位据体层位.点击点击Options / Automatic Picking:点击点击Ok (在在Automatic Picking 菜单菜单菜单菜单) )整个整个整个整个3D3D层位就自动拾取层位就自动拾取层位就自动拾取层位就自动拾取.在层位的平面图上观测在层位的平面图上观测,最左边上的几条最左边上的几条INLINE线上可能拾取的层位位置有问题线上可能拾取的层位位置有问题.112STRATA 培训显示显示inline 2 测线测线(在在STRATA 窗口窗口).用手动方式在同相轴上周围重新拾取层位用手动方式在同相轴上周围重新拾取层位:我们发现自动拾取的结果我们发现自动拾取的结果在一个小区域上层

60、位跳出在一个小区域上层位跳出一段一段.113STRATA 培训现在再重新自动拾取现在再重新自动拾取.在自动拾取菜单上发现缺省参数是清除以前自动拾在自动拾取菜单上发现缺省参数是清除以前自动拾取的层位取的层位,而保留人工拾取的层位作为新引导而保留人工拾取的层位作为新引导:点击点击Ok 得到新层位结果得到新层位结果.人工拾取第二条人工拾取第二条INLINE线后模型得到了改进去引导自动拾取线后模型得到了改进去引导自动拾取:114STRATA 培训下面我们读进层位下面我们读进层位.首先首先删除刚拾取的层位删除刚拾取的层位:层位层位Ok 删除拾取的层位删除拾取的层位.点击点击Horizon / Impor

61、t Horizons / From File:选择五个文件从选择五个文件从 erskine1.pik 到到 erskine5.pik.115STRATA 培训文件类型是缺省文件类型是缺省Geoquest:点击点击Next 选择缺省参数包括名字和层选择缺省参数包括名字和层位颜色位颜色 .点击点击OK 读进层位读进层位.116STRATA 培训层位读进后层位读进后, 再显示再显示inline 24 ,STRATA 窗口显示为窗口显示为:117STRATA 培训点击点击Model / Build/Rebuild a Model:在建模型菜单在建模型菜单,我们修改其中两个参数我们修改其中两个参数,其它

62、用其它用缺省参数缺省参数:现在建立反演所需的初始模型现在建立反演所需的初始模型.118STRATA 培训每个地震道初始猜测模型是由阻抗曲线组成的，阻抗曲线通常是将实际声波曲线乘以实际的密度曲线。阻抗曲线模型是以双程旅行时来计算的，原始的测井曲线是按深度测量的，关键的一步是作深-时转换：测测 井井 曲曲 线线 时时 深深 对对 比比119STRATA 培训深-时转换是用深-时表来完成，转换过程是将每个深度点转到双程旅行时(从基准面(地表)到那个深度和从那个深度回到基准面的时间)：120STRATA 培训到达一个同相轴的时间取决于这个层以上的所有速度,包括到地表的第一层的速度V1，这个速度是未知的

63、,它通常是由测量的第一个速度近似地外推到地表：深深-时表用这个等式从声波曲线速度中计算得来时表用这个等式从声波曲线速度中计算得来:其中其中: ti = 到达第到达第i层的时间层的时间 dj = j层的厚度层的厚度 V j = j层的速度层的速度121STRATA 培训如果井是斜井如果井是斜井,需要校正需要校正(投影投影)成直井并从成直井并从KB校到基准面校到基准面:DM = 从从KB计算的测量深度计算的测量深度DV = 从从KB计算的垂直深度计算的垂直深度DS = 从基准面计算的垂直深度从基准面计算的垂直深度T = 从基准面计算的双程旅行时从基准面计算的双程旅行时DMDVDS122STRATA

64、 培训从声波曲线计算的深-时表很少能获得理想的阻抗模型使之能与地震数据适当的匹配，原因是：n地震基准面与井曲线基准面可能不一样.n第一层的平均速度是未知的.n声波速度的误差会造成计算的传播旅行时有累积误差.n地震上的同相轴由于地震偏移处理误差会造成时间偏移.n地震数据由于频率吸收和短周期多次波等影响可能带有时间拉伸123STRATA 培训改进深改进深-时表可使用两种途径时表可使用两种途径:n运用运用 check shot 校正校正.n应用手动应用手动(交互图形化交互图形化)法将井曲线与地震数据进行对比法将井曲线与地震数据进行对比Check Shot 校正校正一个一个check shot 表是由

65、一组深度表是由一组深度点计算的实际双程旅行时组成的点计算的实际双程旅行时组成的:124STRATA 培训从声波曲线计算的从声波曲线计算的深深-时表需要修改以时表需要修改以获得理想的获得理想的check shot时间时间:原始深原始深/时时曲线曲线理想深理想深/时时曲线曲线125STRATA 培训线性拟合线性拟合: 用直线段内插用直线段内插样条拟合样条拟合: 用平滑曲线内插用平滑曲线内插多项式拟合多项式拟合: 利用最小平方拟合一条光滑曲线利用最小平方拟合一条光滑曲线漂移曲线两点漂移曲线两点间的内插方法间的内插方法有三种选项有三种选项:126STRATA 培训改变深改变深-时表意味着原始声波曲线速

66、度可能发生时表意味着原始声波曲线速度可能发生变化变化. 在在STRATA软件有三种选项软件有三种选项:(1) 井曲线的井曲线的所有速度都改变所有速度都改变使新井曲线可获得使新井曲线可获得理想对应的双程旅行时理想对应的双程旅行时.注释注释: 这意味着在第一个测量深度前用这意味着在第一个测量深度前用斜波速度斜波速度来处理整体时移以最小化合成记录中可能的伪反来处理整体时移以最小化合成记录中可能的伪反射。射。在在STRATA称此选项为称此选项为“Apply All Changes” 。127STRATA 培训(2) 只改变只改变check shot中起始深度和最后深度之中起始深度和最后深度之间的层速度

67、。间的层速度。这是指在第一个测量深度前不加斜坡速度这是指在第一个测量深度前不加斜坡速度.获得新井曲线可得理想对应时间但不考虑整体获得新井曲线可得理想对应时间但不考虑整体时移时移在在STRATA称为称为 “Apply Relative Changes” .128STRATA 培训(3) 不改变声波曲线的速度不改变声波曲线的速度.新井曲线不能对应理想的双程旅行时间，但新井曲线不能对应理想的双程旅行时间，但GEOVIEW 和和STRATA 将利用这个新的深将利用这个新的深-时表时表.这个选项可以保留原始反射系数进行人工合成记这个选项可以保留原始反射系数进行人工合成记录的计算录的计算.在在STRATA

68、称为称为 “Change Depth-Time Table Only” .129STRATA 培训注释注释: 例子中时间的例子中时间的拉伸大的不现实拉伸大的不现实.不同的内插选项可不同的内插选项可能使声波曲线的变能使声波曲线的变化相当大化相当大:130STRATA 培训对比就是对深对比就是对深-时曲线进行时曲线进行手工校正手工校正以优化初始模型和地震数据间的相关以优化初始模型和地震数据间的相关.对比应在对比应在check shot 校正之后进行校正之后进行, 理想情况时深校正对井的理想情况时深校正对井的改变很小改变很小的的.测井曲线对比改变深测井曲线对比改变深-时曲线的方式与时曲线的方式与ch

69、eck shot 校正一样校正一样.测井曲线对比是在合成记录上测井曲线对比是在合成记录上选择几个同相轴选择几个同相轴使之对应于实际地震道上使之对应于实际地震道上.由于使用合成记录由于使用合成记录,子波的选择子波的选择是至关重要的是至关重要的.测井曲线时深对比测井曲线时深对比131STRATA 培训练习练习 3: Erskine 3D 测井曲线时深对比测井曲线时深对比现在在现在在erskine井井上作井对比上作井对比.点击点击 Well / Edit/Correlate Well:选择选择 “erskine_well”作相关作相关:132STRATA 培训在在 eLog 窗口窗口, 点击点击Co

70、rrelate:在提取地震复合道在提取地震复合道Extract Composite Trace 菜单上菜单上,使用缺省参数使用缺省参数, 即使用即使用INLINE和和XLINE旁边旁边(+/- 1)道提取复合道道提取复合道:133STRATA 培训井对比窗口井对比窗口:首先首先, 提取一个新子波提取一个新子波. 由于井曲线尚未由于井曲线尚未对比对比,利用统计法利用统计法Statistical 求零相位求零相位子波而使震幅谱与地震数据的一样子波而使震幅谱与地震数据的一样.134STRATA 培训设置分析窗口使用较小时设置分析窗口使用较小时间间Time 窗口窗口,选择在井旁小选择在井旁小范围范围I

71、NLINE和和XLINE地震地震道道:使用缺省参数使用缺省参数(在第三页在第三页) Statistical Wavelet Extraction 菜单菜单:135STRATA 培训提取的子波为提取的子波为:136STRATA 培训井对比窗口如下井对比窗口如下:可以看到合成记录同可以看到合成记录同相轴与地震同相轴之相轴与地震同相轴之间有一个闭合差。间有一个闭合差。程序也会给出提示建程序也会给出提示建议时移议时移14ms，相关，相关系数最高。系数最高。为了得到更好的效果，为了得到更好的效果，也可以点也可以点Parameters 按钮按钮137STRATA 培训井对比窗口会显示合成记录与地井对比窗口

72、会显示合成记录与地震之间的相关系数。震之间的相关系数。注意：当时移注意：当时移14ms时会得到最大时会得到最大相关系数。相关系数。另外还要注意，通过优化道计算另外还要注意，通过优化道计算时窗也可以使运算结果会得到提时窗也可以使运算结果会得到提高高 。遇到这种情况，先放着。遇到这种情况，先放着。138STRATA 培训现在回到井对比窗口应用建议的时移，点现在回到井对比窗口应用建议的时移，点Apply Shift:发生两种改变发生两种改变 1、 曲线会时移曲线会时移; 2、 相关系数提高相关系数提高139STRATA 培训这个时移是我们时移得到的最好的这个时移是我们时移得到的最好的效果。效果。要想

73、进一步提高相关系数，我们需要想进一步提高相关系数，我们需要手工应用时变时移。要手工应用时变时移。要做这个，可以选择合成道要做这个，可以选择合成道(蓝色蓝色)上上的轴与实际地震道的轴与实际地震道(red)同相轴点中，同相轴点中，然后再点然后再点Stretch 按钮。按钮。140STRATA 培训缺省选项为缺省选项为Spline 漂移曲线间两点进行样条内插漂移曲线间两点进行样条内插.141STRATA 培训改变内插方式为线性改变内插方式为线性 Linear ,点击点击Apply. 注意到漂移曲线注意到漂移曲线形状的变化形状的变化.142STRATA 培训改变菜单选项如下，点改变菜单选项如下，点Ap

74、ply。注意选项。注意选项 Apply all changes 后后后后测井曲线上面多了测井曲线上面多了个斜坡速度，测井曲线值也变了。个斜坡速度，测井曲线值也变了。143STRATA 培训最后最后, 改变菜单如图示点击改变菜单如图示点击 Apply. 然后在然后在Check Shot窗口窗口点击点击 Ok 接接受这些参数受这些参数.144STRATA 培训井对比窗口现在显示井对比窗口现在显示如右图如右图. 注意到相关注意到相关值达到值达到85% .要查看要查看计算相关值使用的参计算相关值使用的参数数,点击点击Parameters 按钮按钮.145STRATA 培训从相关系数窗口可以看到在时间从

75、相关系数窗口可以看到在时间0值处值处出现一个峰值。出现一个峰值。实际上，它建议时移实际上，它建议时移-1ms，再点击，再点击Apply Shift 同样从合成记录相关系数上可以看到不同样从合成记录相关系数上可以看到不需要再做相位调整了。需要再做相位调整了。146STRATA 培训井对比结束井对比结束, 在在 eLog 窗口底部点击窗口底部点击OK.下一个菜单允许给新生成的声波下一个菜单允许给新生成的声波曲线命名曲线命名.点击点击Ok 选择缺省名选择缺省名(P-wave_corr).最后点击最后点击File / Exit Window 退出退出eLog 窗口窗口.147STRATA 培训现在运行

76、基于模型反演建立初始模型：现在运行基于模型反演建立初始模型：需要分两步完成：第一步，在井位处应需要分两步完成：第一步，在井位处应用反演决定反演参数并且让程序决定最用反演决定反演参数并且让程序决定最优比例因子。点优比例因子。点Analysis / Post-stack Analysis / Model Based:在第一个菜单项，选择在第一个菜单项，选择ersk3d 作为反演输入数据，再点作为反演输入数据，再点 Next148STRATA 培训第二项，我们确定应用了正确的子波，点第二项，我们确定应用了正确的子波，点 Set Current Wavelet 来设置来设置图示显示之前提取的子波，已经

77、设图示显示之前提取的子波，已经设置对了，点置对了，点Cancel 退出窗口。退出窗口。然后点然后点 Next and Ok 来接受默来接受默认参数和产生反演分析窗口。认参数和产生反演分析窗口。149STRATA 培训反演分析窗口显示了一系反演分析窗口显示了一系列有用的曲线，这些曲线列有用的曲线，这些曲线帮助我们确认反演运行正帮助我们确认反演运行正确。确。SyntheticReal DataErrorInversion TraceInitial ModelReal Log150STRATA 培训许多显示参数也可以通过点眼睛来设置：许多显示参数也可以通过点眼睛来设置：选择选择 Curves 标签标

78、签选择选项应用滤波如右选择选项应用滤波如右图，然后点击图，然后点击Ok151STRATA 培训从这个显示可以看到，应用从这个显示可以看到，应用2 QC标准标准 得得到一个好的反演结果。反演道（红色）到一个好的反演结果。反演道（红色）与地震道（蓝色、滤波后的）相关系数与地震道（蓝色、滤波后的）相关系数很高，同时，合成记录与地震道之间的很高，同时，合成记录与地震道之间的误差也接近于误差也接近于0。如果喜欢的话，我们可如果喜欢的话，我们可以修改反演参数，点击以修改反演参数，点击 Apply 查看新的结果。查看新的结果。然而，这个反演体已经然而，这个反演体已经足够好了，继续运行。足够好了，继续运行。点

79、点 File / Exit 退出分析退出分析窗口。窗口。152STRATA 培训在结果菜单在结果菜单, 所有的默认参数是正所有的默认参数是正确的。鉴于我们确的。鉴于我们 在分析时给出的在分析时给出的参数。点击参数。点击 Next 直到达到比例直到达到比例因子选项而，可以看到井点处产生因子选项而，可以看到井点处产生的比例会应用于整个工区。点击的比例会应用于整个工区。点击 Ok 对整个工区进行反演。对整个工区进行反演。现在，对整个工区应用反演。点击现在，对整个工区应用反演。点击 Inversion / Post-stack Inversion / Model Based Inversion:153

80、STRATA 培训单条线反演结束后得到以下剖面单条线反演结束后得到以下剖面:154STRATA 培训检查反演结果的质量一种方法是绘制误差图检查反演结果的质量一种方法是绘制误差图Error Plot. 这这是利用反演结果计算的合成记录和实际地震道之间的差别是利用反演结果计算的合成记录和实际地震道之间的差别. 点击点击 眼球眼球“eyeball” 图标图标.设设 Trace Data Volume 为为 “model_based derived Synthetic Error”.点击点击 Apply 观测误差图观测误差图.155STRATA 培训误差图的加权与实际地震数据一样误差图的加权与实际地震

81、数据一样. 结果是几乎没有相关误结果是几乎没有相关误差存在差存在,表明求得阻抗模型可以忠实地反映了地震数据表明求得阻抗模型可以忠实地反映了地震数据.(练习练习 3 结束结束)点点 File / Exit Project 退出工退出工区。区。156STRATA 培训模型中使用一口井产生均衡的水平模型模型中使用一口井产生均衡的水平模型:测井曲线层间插值建模157STRATA 培训注释注释: 对每道而言对每道而言, 单个拾取的同相轴只能在每道进行简单井曲线的块状移动单个拾取的同相轴只能在每道进行简单井曲线的块状移动. 这这个等同于用单点进行个等同于用单点进行check shot 校正校正.拾取一个同

82、相轴指导井的内插拾取一个同相轴指导井的内插:158STRATA 培训拾取多个同相轴相当于在每道运用可变的拾取多个同相轴相当于在每道运用可变的check-shot校正校正. 在两个层位之间，其阻抗进行相应的拉伸在两个层位之间，其阻抗进行相应的拉伸/挤压挤压.159STRATA 培训层位间内插由层位间内插由STRATA中的模型选项来控制中的模型选项来控制:160STRATA 培训如右图有三种内插方法如右图有三种内插方法.缺省为所有层位间除第一层和最后一层缺省为所有层位间除第一层和最后一层都用整合面来处理都用整合面来处理.顶超顶超底超底超整合整合161STRATA 培训尖灭点可以通过使两个层位聚合来

83、处理尖灭点可以通过使两个层位聚合来处理:162STRATA 培训如果层位交叉如果层位交叉, 每个层位被赋予优先权值来处理每个层位被赋予优先权值来处理.这种情形这种情形, H1 优先级高优先级高,故故H2层位被删层位被删这是相反的例子这是相反的例子.163STRATA 培训STRATA 目前不处理断层目前不处理断层目前不处理断层目前不处理断层. 但是可以在断层的两侧但是可以在断层的两侧,断层面拾断层面拾取同样的层位来模拟取同样的层位来模拟 :164STRATA 培训模型中有多口井时模型中有多口井时,井之间可用反距离加权法内插井之间可用反距离加权法内插:165STRATA 培训其中其中:一般公式一

84、般公式:假设有两口输入井曲线假设有两口输入井曲线 L1 和和 L2. 我们希望计算输出我们希望计算输出 Lout.这是两个输入曲线的线性组合这是两个输入曲线的线性组合: Lout = w1*L1 + w2*L2从目标点到每个输入曲线从目标点到每个输入曲线,加权变化与距离成加权变化与距离成反比反比:166STRATA 培训井间的内插选项如下井间的内插选项如下:167STRATA 培训较少使用的选项是三角内插法较少使用的选项是三角内插法, 在相邻井拟在相邻井拟合一系列平面段合一系列平面段 和克里金法和克里金法, 其要求有变量图其要求有变量图输入输入:168STRATA 培训多口井多个层位使用反距离

85、加权法内插由层位指导多口井多个层位使用反距离加权法内插由层位指导:169STRATA 培训加层位的一般准则加层位的一般准则:(1) 拾取的层位要贯穿拾取的层位要贯穿整个工区整个工区.没有拾取的点程序将进行内插没有拾取的点程序将进行内插.(2) 只拾取可以确定的层位只拾取可以确定的层位.(3) 拾取大的构造层位拾取大的构造层位, 不拾取细节层位不拾取细节层位.(4) 通常解释系统拾取的层位正好可以是通常解释系统拾取的层位正好可以是STRATA软件所需要的层位软件所需要的层位.有层位和没有层位内插的有层位和没有层位内插的比较比较:170STRATA 培训STRATA软件缺省自动在内插后对模型通过软

86、件缺省自动在内插后对模型通过高截滤波进行平高截滤波进行平滑滑,使模型只含有低频成分使模型只含有低频成分.这样可以避免模型中的高频成分干这样可以避免模型中的高频成分干扰反演结果扰反演结果.171STRATA 培训你可以检查通过内插你可以检查通过内插模型道得到的高频模模型道得到的高频模型型,但不进行任何平滑但不进行任何平滑处理处理.还有一种选项是对相邻还有一种选项是对相邻层位进行模型分块化层位进行模型分块化.172STRATA 培训High Frequency ModelSmooth ModelInversion ResultInversion Result用高频模型反演与平滑初始模型反演后的对比

87、：用高频模型反演与平滑初始模型反演后的对比：173STRATA 培训用用高频模型反演高频模型反演的结果的结果平滑模型反演平滑模型反演的结果的结果两个结果非常相似，但第二两个结果非常相似，但第二个更可靠个更可靠 。174STRATA 培训褶积模型是反演基础：褶积模型是反演基础： 地震道地震道 = 子波子波 * 反射系数反射系数 + 噪音噪音在频率域，褶积对应相乘：在频率域，褶积对应相乘：反演结果可以通过除子波得到：反演结果可以通过除子波得到： 反射系数反射系数 = 地震道地震道 / 子波子波通过窄带子波可以限制频率域有用信息的范围。通过窄带子波可以限制频率域有用信息的范围。子波提取子波提取175

88、STRATA 培训从频率域，相位谱很可能从频率域，相位谱很可能近似表达为一条直线。近似表达为一条直线。The intercept of the line is the constant phase rotation which best characterizes this wavelet.The slope of the line measures the time-shift of the wavelet.这条线的截距对应子波常这条线的截距对应子波常相位旋转，可以很好的描相位旋转，可以很好的描述这个子波。述这个子波。这条线的梯度表示子波时这条线的梯度表示子波时移程度移程度子波可以完全通过震

89、幅谱和相位谱来描述：子波可以完全通过震幅谱和相位谱来描述：176STRATA 培训0o45o90o180o-90o这些子波有同样的振幅谱，但是相位谱却存在差这些子波有同样的振幅谱，但是相位谱却存在差异。异。177STRATA 培训子波时移子波时移可以改变可以改变相位谱的梯度相位谱的梯度，但，但截距保持不变截距保持不变。178STRATA 培训子波相位和时移曲线之子波相位和时移曲线之间存在一个基本的间存在一个基本的双解双解问题。问题。考虑到这点，如果波阻考虑到这点，如果波阻抗没有对比好，那么存抗没有对比好，那么存在一个闭合差。在一个闭合差。最好的办法是时移这条最好的办法是时移这条曲线曲线 (改变

90、时深曲线改变时深曲线). 然后再提取子波。然后再提取子波。ShiftExtract179STRATA 培训理论上说理论上说, 闭合差也可以闭合差也可以通过时移子波来解决。通过时移子波来解决。然而这样做是然而这样做是 不理想的不理想的 因为：因为：(1)如果有多口井存在这如果有多口井存在这个问题个问题, 一个子波的一个子波的时移不能解决全部井时移不能解决全部井的问题。的问题。(2)用时移子波做反演会用时移子波做反演会错误定位阻抗值到相错误定位阻抗值到相应的时间上应的时间上Dont ShiftExtract180STRATA 培训当用当用STRATA 提取子波提取子波时，它可以检测一口井时，它可以

91、检测一口井上需不需要做时移并且上需不需要做时移并且给出提示：给出提示：用户可以通过点击用户可以通过点击 Apply Shift 来来应用这个时移就可以了应用这个时移就可以了181STRATA 培训极性约定可以在极性约定可以在 Synthetic Polarity Convention 菜单进行设置：菜单进行设置：而而 alternate convention 的零相位地震数据默认阻抗增加对应波谷响应；的零相位地震数据默认阻抗增加对应波谷响应；一个特殊的子波相位问题就是极性约定（一个特殊的子波相位问题就是极性约定（ Polarity Convention）常规的零相位地震数据默认阻抗增加对应波峰

92、响应：常规的零相位地震数据默认阻抗增加对应波峰响应：182STRATA 培训子波在地球内部传播可能会空变，也可能会时变，原因有下：子波在地球内部传播可能会空变，也可能会时变，原因有下：n地表因素影响地表因素影响 (空变空变)n频率域吸收衰减频率域吸收衰减 (空变和时变空变和时变)n层间多次波层间多次波 (空变和时变空变和时变)n动校拉伸动校拉伸n处理上的人为因素处理上的人为因素STRATA假定子波是常相位的假定子波是常相位的 （不管时变还是空间）（不管时变还是空间）:n时变时变: 意味着只能在限定的时窗内进行反演。意味着只能在限定的时窗内进行反演。n空变空变: 这意味着对这个数据要选择性的进行

93、反演以消除子波的空变带来的这意味着对这个数据要选择性的进行反演以消除子波的空变带来的影响。影响。183STRATA 培训子波提取有两种方法：子波提取有两种方法：(1)用用 井和地震井和地震共同来提取子共同来提取子波及其对应的振幅谱和相位波及其对应的振幅谱和相位谱谱(2)只用只用地震地震 来提取子波的振来提取子波的振幅谱，这时假定子波的相位幅谱，这时假定子波的相位是是零相位零相位的。的。184STRATA 培训为什么不只用从井上为什么不只用从井上提取子波提取子波?因为因为井相关误差井相关误差 (拉拉伸伸) 可能导致非常大可能导致非常大的相位问题的相位问题。解决方案解决方案: 用井和地用井和地震提

94、子波前先从地震震提子波前先从地震上提子波并做好井的上提子波并做好井的对比工作。对比工作。ExtractExtract185STRATA 培训这一步只用自动相关，相位假定是已知的：这一步只用自动相关，相位假定是已知的：主要参数：主要参数：道范围道范围 (通常设定这个大一些来提高统计效果）通常设定这个大一些来提高统计效果） 时窗时窗 (应该至少应该至少2倍子波长度倍子波长度)子波长度子波长度统计子波提取统计子波提取 （不用井）（不用井）186STRATA 培训统计子波提取实例统计子波提取实例: 注意子波是零相位的，因为这是已设定好的。注意子波是零相位的，因为这是已设定好的。187STRATA 培训

95、用井进行子波提取用井进行子波提取:188STRATA 培训这一步用井曲线来提取子波的振幅谱和相这一步用井曲线来提取子波的振幅谱和相位谱。它的好坏很大程度上信赖于之前井位谱。它的好坏很大程度上信赖于之前井和地震对比的质量。和地震对比的质量。主要参数主要参数:n选择井选择井 n时窗时窗n子波长度子波长度n提取类型（提取类型（常相位常相位）用井进行统计子波提取用井进行统计子波提取:189STRATA 培训它可以通过解决时域算子（描述地震合成道对应的井上的反射系数）问题来精确地它可以通过解决时域算子（描述地震合成道对应的井上的反射系数）问题来精确地提取振幅谱和相位谱，这项只有在提取振幅谱和相位谱，这项

96、只有在对比非常好的情况下使用对比非常好的情况下使用：Full Wavelet 选项选项:190STRATA 培训通过与地震道之间的互相关，用统计方法准确地计算子波的振幅谱和相位谱，相位谱通过与地震道之间的互相关，用统计方法准确地计算子波的振幅谱和相位谱，相位谱可以粗略的记作一个常值。这个方法可以粗略的记作一个常值。这个方法比比Full Wavelet要稳定要稳定，尤其是井震对比不是太，尤其是井震对比不是太好的情况下，一般把它作为默认选项。好的情况下，一般把它作为默认选项。Constant Phase 选项选项:191STRATA 培训计算相位谱的流程计算相位谱的流程:(1) 用统计子波提取方法

97、用统计子波提取方法提取子波提取子波 (不用井不用井).(2) 对提取的子波用一系列的对提取的子波用一系列的常相位旋转常相位旋转。(3) 对每一个相位旋转后的子波，生成对每一个相位旋转后的子波，生成合成记录道与地震道做对比合成记录道与地震道做对比。 (4) 选择得到选择得到最大相关系数的相位最大相关系数的相位。192STRATA 培训这一项应用了这一项应用了Roy White (White, R.E., and Simm, R., 2003, Tutorial:Good practice in well ties: First Break 21, 75-83）发明的一个算法。）发明的一个算法。提

98、取方法提取方法与与Full Wavelet方法相似方法相似，另外，还会生成一个鉴定结果允许用户估算子波，另外，还会生成一个鉴定结果允许用户估算子波提取中的误差。提取中的误差。Roy White 选项选项:193STRATA 培训Full Wavelet如果井对比很好，各方法提取的子波会如果井对比很好，各方法提取的子波会比较相似：比较相似：Constant PhaseRoy White194STRATA 培训子波提取中的一个基本问题子波提取中的一个基本问题:n用井提取子波前，必须用井提取子波前，必须优化对比结果优化对比结果，把它做好。，把它做好。n要合理地应用对比，要合理地应用对比，子波必须已知

99、子波必须已知。子波提取流程子波提取流程:(1) 用统计子波提取用统计子波提取 一个一个初始的子波初始的子波，这意味着子波的相位近似已知，这意味着子波的相位近似已知 .(2) 对曲线进行对曲线进行拉伸拉伸与地震对齐与地震对齐.(3) 用井数据提取一个用井数据提取一个新的子波新的子波.(4) 再再重复第重复第 (2) 和和 (3)步步多次得最终结果多次得最终结果195STRATA 培训对子波判别非常有用的相关窗口对子波判别非常有用的相关窗口当前井当前井时移后时移后的最大相关的最大相关建议时移。只有在没做建议时移。只有在没做拉伸前使用拉伸前使用两侧是两侧是平行对称平行对称的说明没有剩余相的说明没有剩

100、余相位误差位误差196STRATA 培训如果有多口井，一个很好的判别方法是如果有多口井，一个很好的判别方法是 Multi-well Analysis:197STRATA 培训相关图相关图 展示了每一口井的相关系数，展示了每一口井的相关系数， 它可以用来它可以用来标记不好的井，在接下来子波提取中进行删除。标记不好的井，在接下来子波提取中进行删除。Multi-well analysis198STRATA 培训练习练习 4: Blackfoot 子波提取练习子波提取练习这个练习这个练习, 对一个新数据体作反演对一个新数据体作反演.这个数据来自加拿大西部的这个数据来自加拿大西部的Blackfoot(印

101、地安人种族印地安人种族) 地区地区, 整块三维有整块三维有13口井口井. 这这些井已经加载到些井已经加载到GEOVIEW 数据库数据库. 第一步就是打开那个数据库第一步就是打开那个数据库.选择数据库选择数据库 “blackfoot.wdb”点击点击Database / Open (在在GEOVIEW 窗口窗口):199STRATA 培训GEOVIEW井数据库井数据库( Well Explorer)窗口如下窗口如下:查看任一口井，查看任一口井，发现都有声波和发现都有声波和密度曲线密度曲线.我们已经修改了外部数据库我们已经修改了外部数据库, 故需要完整退出故需要完整退出STRATA软件软件. 然后

102、重新启动然后重新启动STRATA软件软件(在在 GEOVIEW 窗口点击窗口点击STRATA 按钮按钮).200STRATA 培训启动新工区启动新工区:命名为命名为 “blackfoot_project”:点击点击Data Manager / Import Data / Open Seismic / From SEG-Y 打开打开blackfoot 3D数据数据:文件为文件为 “blackfoot_seismic.sgy”:201STRATA 培训这次利用道头中这次利用道头中Inline & Xline 读地震数据读地震数据 :第四页第四页, 检查道头以确认检查道头以确认Inline and

103、Xline 数是正确的数是正确的. 如下点击如下点击Detail Specification :202STRATA 培训新菜单出现带有数页菜单新菜单出现带有数页菜单. 第一页包第一页包含测线道头的一般信息含测线道头的一般信息 .内容如果有内容如果有误可以修改误可以修改.点击点击Next 两次得到第两次得到第三页三页.203STRATA 培训这页显示这页显示Inline 和和 Xline 字节位置字节位置.第四页允许观测实际道头数值第四页允许观测实际道头数值. 各道头各道头正确点击正确点击OK 按钮采用标准的按钮采用标准的SEGY 格格式式.204STRATA 培训扫描后扫描后, 观测系统几何观

104、测系统几何Geometry 页显示如下页显示如下: 点击点击 Ok 加载数据。加载数据。205STRATA 培训井到地震成像是正确的因为井到地震成像是正确的因为在在GEOVIEW软件中井的软件中井的X&Y坐标已经设置好坐标已经设置好.点击点击 Ok (在在Well To Seismic Map菜单菜单).STRATA窗口中的地震显示窗口中的地震显示.206STRATA 培训观测井的位置在观测井的位置在STRATA窗口点击窗口点击View / Base Map :底图为底图为:207STRATA 培训修改参数如下显示修改参数如下显示Xline 42 测线测线:208STRATA 培训在显示中插入

105、不同的井曲线在显示中插入不同的井曲线. 点击眼球点击眼球 “eyeball” 图标出现菜单图标出现菜单.修改参数菜单中修改参数菜单中Insert 页如右图页如右图. 然然后点击后点击 Apply 可看到下图可看到下图.在加载密度曲线后在加载密度曲线后, 点击点击Cancel(参数参数菜单菜单)恢复显示声波曲线恢复显示声波曲线.209STRATA 培训现在读进层位现在读进层位.点击点击 Horizon / Import Horizons / From File:选择文件选择文件 “blackfoot_horizons.txt”.这文件有两个层位这文件有两个层位(以自由以自由格式格式) Free

106、Format:210STRATA 培训下一页定义有两个层位下一页定义有两个层位:最后一页最后一页, 点击层位名点击层位名,然后点击然后点击 “Display selected file” 按钮便按钮便可以确定格式参数可以确定格式参数: 211STRATA 培训层位文件格式层位文件格式:如左图插入参数如左图插入参数. 改变层位名字改变层位名字.212STRATA 培训层位读进后下一步将对每口井进行对比层位读进后下一步将对每口井进行对比.实际上所有井已经做过对比实际上所有井已经做过对比. 检查其检查其中一口井中一口井.点击点击 井井well “08-08” 然后点击然后点击 Edit 按钮按钮.点

107、击点击 Well / Edit/Correlate Well:213STRATA 培训在在eLog 窗口窗口,点击点击 Correlate button.点击点击 Ok (在在Extract Composite Trace 菜单菜单)利用缺省值提取复合地震道利用缺省值提取复合地震道.214STRATA 培训先看一下，对比比较差先看一下，对比比较差.点击点击 Wavelet / Display Current Wavelet 看缺省的子波看缺省的子波:但是在提取更好子波前不应做决定但是在提取更好子波前不应做决定.215STRATA 培训提取统计子波点击提取统计子波点击 Wavelet / Ext

108、ract Wavelet / Statistical:设时间窗口为设时间窗口为800 to 1200 ms, 而而Inline范围范围设为一条设为一条Inline线线.子波参数子波参数(Wavelet Parameters)页使用缺省参数页使用缺省参数:216STRATA 培训提取的子波为提取的子波为:注意高频成分注意高频成分.217STRATA 培训注意到如果子波有误而进行井对比注意到如果子波有误而进行井对比, 其结果将会使误差变得不可估量。其结果将会使误差变得不可估量。测量的相关系数是只有测量的相关系数是只有0.46, 但这可能是计算窗口影响但这可能是计算窗口影响. 点击点击 Parame

109、ters 按钮按钮:用新子波相关已经变得更好用新子波相关已经变得更好.218STRATA 培训时移时移6ms后显示相关系数变好后显示相关系数变好（%70） ，点击点击 Apply Shift 应用：应用：219STRATA 培训现在，微调相关系数。实际上，现在，微调相关系数。实际上，已经非常好了。选择右侧的同相已经非常好了。选择右侧的同相轴并点击轴并点击Stretch.然后点然后点Ok 应用应用Stretch220STRATA 培训相关系数显示了剩余相位误差，我们可以试相关系数显示了剩余相位误差，我们可以试图校正这口井，但是实际上我们需要一个满图校正这口井，但是实际上我们需要一个满足所有井的最

110、优子波。足所有井的最优子波。要做到这点，我们需要用井来提取子波。要做到这点，我们需要用井来提取子波。点击点击 Wavelet / Extract Wavelet / Use Well.221STRATA 培训在第一个菜单而，选择在第一个菜单而，选择Add All 来应用所来应用所有的井。有的井。设置时窗从设置时窗从 800 到到 1200 ms.222STRATA 培训在最后一页上，对这个子波应用所有默认参在最后一页上，对这个子波应用所有默认参数，点击数，点击 Next and Ok 来提取子波。来提取子波。223STRATA 培训现在显示的相关系数图形上显示对称的，同现在显示的相关系数图形上

111、显示对称的，同样建议再一次时移，所以点击样建议再一次时移，所以点击Apply Shift 按钮再应用一次。按钮再应用一次。我们已经完成了这口井上的对比工作，点击我们已经完成了这口井上的对比工作，点击Ok 接受新的纵波时差名称：接受新的纵波时差名称：最后，点击最后，点击File / Exit Window 退出窗口。退出窗口。224STRATA 培训我们现在完成了我们现在完成了13口井中的一口井的对比工作，并且为所有的井提取了一口井中的一口井的对比工作，并且为所有的井提取了一个子波。个子波。在一个真实工区中，我们必须为接下来的在一个真实工区中，我们必须为接下来的12口井也完成同样的工作，为了口井

112、也完成同样的工作，为了节省时间，我们已经完成了这部分工作，假定其它井已对比好了。节省时间，我们已经完成了这部分工作，假定其它井已对比好了。一个完整的多井工区基本的工作流程如下：一个完整的多井工区基本的工作流程如下：(1)提取一个提取一个统计子波统计子波；(2)对对每一口井每一口井完成同一上述工作，做好对比；完成同一上述工作，做好对比；(3)提取所有井的提取所有井的井旁道子波井旁道子波；(4)每一口井每一口井做一遍，精细微调对比，通常这一步只需要应用一个时移就做一遍，精细微调对比，通常这一步只需要应用一个时移就可以了。可以了。(练习练习 4结束结束)225STRATA 培训基于模型反演的参数基于

113、模型反演的参数主要参数有主要参数有: 平均块大小 反复迭代次数 比例化类型次要参数：次要参数： 反演选项 最大阻抗变化基于模型反演菜单基于模型反演菜单:226STRATA 培训由于由于STRATA 通过一系列递推聚合收敛，这个参数决定收敛的程度。实际应用中，通过一系列递推聚合收敛，这个参数决定收敛的程度。实际应用中，用用3次递推次递推就可完成大部分工作就可完成大部分工作.递推次数多递推次数多没有任何坏处没有任何坏处,它只影响运行时间它只影响运行时间.递推多少次才使结果递推多少次才使结果收敛取决于反演中的块大小，收敛取决于反演中的块大小，块状小块状小需要需要更多的递推次数更多的递推次数.检查递归

114、次数是否足够可通过检查递归次数是否足够可通过误差图进行分析误差图进行分析.反复递推次数反复递推次数:227STRATA 培训最后反演结果可能会改变块的大小最后反演结果可能会改变块的大小, 但是块的总数没有变化但是块的总数没有变化. 这意味着有些块变大这意味着有些块变大而有些块变小而有些块变小,而平均保持常数不变而平均保持常数不变.利用利用小块状小块状如如 (2 ms) 将将提高分辨率提高分辨率, 但增加的细节可能来自初始猜测阻抗模型但增加的细节可能来自初始猜测阻抗模型.利用利用小块状小块状总会总会改善改善最后反演道与输入地震的拟合最后反演道与输入地震的拟合.此参数控制最后结果的分辨率此参数控制

115、最后结果的分辨率. 初始猜测阻抗模型被分块化成一系列均衡的小块初始猜测阻抗模型被分块化成一系列均衡的小块:平均块状大小平均块状大小228STRATA 培训平均块大小平均块大小反演用反演用6ms块大小块大小:用用2ms块大小块大小建议：设置块大小与采样间隔一致，但是要用平滑的初始模型。建议：设置块大小与采样间隔一致，但是要用平滑的初始模型。229STRATA 培训主要反演菜单设置好后，接下来显示数据的比例因子：主要反演菜单设置好后，接下来显示数据的比例因子：比例参数比例参数230STRATA 培训反演的基本理论褶积模型：反演的基本理论褶积模型：地震道地震道 = 子波子波 * 反射系数反射系数 +

116、 噪音噪音频率域频率域:反射系数反射系数=地震道地震道 / 子波子波求反射系数，子波需要已知。求反射系数，子波需要已知。正常情况下，子波提取后，它的形态是已知的，但绝对振幅强度是未知的。反演需要正常情况下，子波提取后，它的形态是已知的，但绝对振幅强度是未知的。反演需要绝对振幅也已知。绝对振幅也已知。从上面的等式可知，如果子波乘以从上面的等式可知，如果子波乘以2，结果反射系数会是原来的，结果反射系数会是原来的1/2。STRATA 通过初始模型合成道（滤波前）的均方根振幅与实际地震道的均方根振幅通过初始模型合成道（滤波前）的均方根振幅与实际地震道的均方根振幅相同来自动定义子波的比例大小。相同来自动

117、定义子波的比例大小。为什么要设置比例为什么要设置比例?231STRATA 培训定义比例因子唯一的方法是看最终反演道与初始猜测曲线匹配的好不好。定义比例因子唯一的方法是看最终反演道与初始猜测曲线匹配的好不好。求取整个数据体的单比例因子（求取整个数据体的单比例因子（single scalar for dataset）对每一道求取一个比例因子（对每一道求取一个比例因子（separate scalar for each trace）第一项，单比例因子第一项，单比例因子原理上更理想原理上更理想，因为它假定子波用一个比例用于整个数据体，这，因为它假定子波用一个比例用于整个数据体，这样可以保持道与道之间的振

118、幅变化。样可以保持道与道之间的振幅变化。 第二项，选用不同比例因子，第二项，选用不同比例因子，更稳定更稳定。它有效地假定道与道之间的变化可以通过重新。它有效地假定道与道之间的变化可以通过重新比例化来消除，这并不是基于岩性数据。比例化来消除，这并不是基于岩性数据。不管是哪个选项，在实际反演中选用不管是哪个选项，在实际反演中选用不同时窗对应不同的比例因子不同时窗对应不同的比例因子。对于许多数据设。对于许多数据设置，尤其稀疏模型，自动比例可能并不理想。这种情况下，可能需要手工调整，给出置，尤其稀疏模型，自动比例可能并不理想。这种情况下，可能需要手工调整，给出多个自动比例因子。多个自动比例因子。比例因

119、子比例因子232STRATA 培训比例对比比例对比蓝色蓝色 = 初始阻抗初始阻抗曲线曲线红色红色 = 反演阻抗反演阻抗结果结果Scaling too lowScaling too highScaling just right233STRATA 培训比例选项比例选项建议建议: 用反演分析在井旁道处求得一个全局比例因子。用反演分析在井旁道处求得一个全局比例因子。234STRATA 培训这个参数控制如何使用约束这个参数控制如何使用约束.基于模型反演基于模型反演 用以下形式使目标函数最小化用以下形式使目标函数最小化:J = weight1 x (T - W*r) + weight2 x (M - H*

120、r)其中其中:T = 地震道地震道W= 子波子波 r = 最后反射系数最后反射系数M= 初始猜测阻抗模型初始猜测阻抗模型H = 综合因子与最后反射系数褶积得到最后阻抗综合因子与最后反射系数褶积得到最后阻抗* = 褶积褶积反演选项反演选项235STRATA 培训在模型约束在模型约束( Model Constraint)( Model Constraint)中的参数是目标函数中的加权量中的参数是目标函数中的加权量 weight weight2 2这个参数设为这个参数设为 0 0 使地震数据起决定作用；使地震数据起决定作用；参数设为参数设为 1 1 使初始猜测阻抗起决定作用；使初始猜测阻抗起决定作用

121、；这称为这称为软约束软约束因为最终模型可与初始猜测阻抗背离可大可小甚至完全背离。因为最终模型可与初始猜测阻抗背离可大可小甚至完全背离。最小化第一部分最小化第一部分 (T - W*r)，使模型与地震道逼近，使模型与地震道逼近.最小化第二部分最小化第二部分 (M - H*r)，是通过利用定义的块大小逼近初始猜测阻抗，是通过利用定义的块大小逼近初始猜测阻抗.这两种方法其条件通常不兼容这两种方法其条件通常不兼容. 加权：加权：weight1 和和weight2决定两部分如何均衡决定两部分如何均衡. 在在软约束软约束(Soft Constraint inversion) 选项选项, 目标函数正如以上所示

122、目标函数正如以上所示. 加权是由这加权是由这个参数确定个参数确定:目标函数有两部分目标函数有两部分: 236STRATA 培训在实际应用中在实际应用中, 通常优先选择通常优先选择硬约束硬约束因为因为最大阻抗变化最大阻抗变化比比模型约束模型约束更现实合理更现实合理.在最大阻抗变化在最大阻抗变化(Maximum Impedance Change) 是井曲线平均误差的百分比是井曲线平均误差的百分比. 注意对于这个模型有效的范围为注意对于这个模型有效的范围为:在在硬约束硬约束(Hard Constraint inversion), 目标函数中目标函数中只有第只有第1项项。但是。但是, 算法被约束为其的

123、最后阻抗值将被约束到以下参数算法被约束为其的最后阻抗值将被约束到以下参数所定义的范围内所定义的范围内:这称为这称为硬约束硬约束(hard constraint), 因为因为最后阻抗值不允许变化超过给定的边界最后阻抗值不允许变化超过给定的边界范围范围.237STRATA 培训误差图误差图误差图误差图 (Error Plot )显示利用反演阻抗结果计算的合成记录和实际地震的区别显示利用反演阻抗结果计算的合成记录和实际地震的区别:理想状态理想状态,误差图不应有相干能量而应该是一些总体很低的震幅误差图不应有相干能量而应该是一些总体很低的震幅.238STRATA 培训误差中的低频成分误差中的低频成分-可

124、能是子波有误可能是子波有误:239STRATA 培训误差出现在某一边区域化误差出现在某一边区域化-可能是使用的层位可能不够可能是使用的层位可能不够:240STRATA 培训整个数据体都有相干误差整个数据体都有相干误差 可能是可能是: n反演块状太大反演块状太大(SIZE) n递推次数不够递推次数不够n约束的太紧（硬约束）约束的太紧（硬约束）241STRATA 培训练习练习 5: Blackfoot 基于模型反演基于模型反演现在给现在给Blackfoot数据体建初始模型数据体建初始模型.点击点击 Model / Build/Rebuild a Model.称为称为 Model_1.使用缺省参数但

125、使用所有使用缺省参数但使用所有井井.242STRATA 培训点击点击 Model / Set Model Parameters:现在分析建模型中的各个参数现在分析建模型中的各个参数.243STRATA 培训每页都看每页都看.点击点击 Advanced 标签标签, 如右图改变如右图改变.点击点击 Apply .看到看到STRATA窗口利用窗口利用这些新参数进行刷新这些新参数进行刷新.菜单由一系列分页组成菜单由一系列分页组成,控制模型建立各方面的参数控制模型建立各方面的参数.244STRATA 培训这个高频模型并不适合去反演，所以点击这个高频模型并不适合去反演，所以点击 Cancel 恢恢复平滑模

126、型。复平滑模型。模型显示了高频细节，用没滤波的井插值的结果。模型显示了高频细节，用没滤波的井插值的结果。245STRATA 培训我们需要用默认模型参数完成反演。正如之前练习，需要用两步完成反演。我们需要用默认模型参数完成反演。正如之前练习，需要用两步完成反演。选择体选择体 blackfoot_seismic:点击点击 Analysis / Post-stack Analysis / Model Based:采用所有默认参数，点击采用所有默认参数，点击 Next and Ok.246STRATA 培训当出现反演分析窗口。注当出现反演分析窗口。注意，可以通过点击箭头滚意，可以通过点击箭头滚动查看每

127、一口井。动查看每一口井。247STRATA 培训可以通过对波阻抗滤波提高显示效果。这时，可以通过对波阻抗滤波提高显示效果。这时，我们需要查看地震数据的频带范围，在我们需要查看地震数据的频带范围，在 STRATA 主菜单，点击主菜单，点击 Process / Utility / Amplitude Spectrum 并且填入如下参数：并且填入如下参数：248STRATA 培训现在，通过点眼睛修改如下参现在，通过点眼睛修改如下参数来修改反演分析显示：数来修改反演分析显示：249STRATA 培训通过放大目标区域，我通过放大目标区域，我们可以看到反演结果有们可以看到反演结果有点高，如果这个确实是点高

128、，如果这个确实是研究目标区域，并且在研究目标区域，并且在其它井处是准确的。我其它井处是准确的。我们可能希望修改它，应们可能希望修改它，应用比例因子来调整。用比例因子来调整。然而，对于这个练习，我们采然而，对于这个练习，我们采用默认参数。用默认参数。250STRATA 培训另一个有用的显示是另一个有用的显示是Error Profile. 在在 Error Analysis / Logs / Error Profile:这个显示了每口井的误差这个显示了每口井的误差 (实际实际阻抗与预测阻抗之差）。阻抗与预测阻抗之差）。如果有很大的误差我们需要检查如果有很大的误差我们需要检查那口井，看可能存在的问题

129、。那口井，看可能存在的问题。现在，关掉反演分析窗口现在，关掉反演分析窗口File / Exit.251STRATA 培训既然在井位处检查了参数，并且定义了既然在井位处检查了参数，并且定义了优化比例因子，现在对整个工区应用反优化比例因子，现在对整个工区应用反演。点击演。点击 Inversion / Post-stack Inversion / Model Based Inversion:既然我们已经得到了参数，在这个菜单既然我们已经得到了参数，在这个菜单除了输出体的名字外不需要再做任何事除了输出体的名字外不需要再做任何事情了。情了。当你已经改变了这项，点当你已经改变了这项，点Next and O

130、k 应用整个反演。大概需要应用整个反演。大概需要10分钟。分钟。252STRATA 培训当反演完成后，如下图，注意你可以通过点击箭头来查看其它线。当反演完成后，如下图，注意你可以通过点击箭头来查看其它线。253STRATA 培训插入井上的阻抗曲线插入井上的阻抗曲线很有用，可以通过点很有用，可以通过点击击“眼睛眼睛” 改变菜单改变菜单如下，应用如下，应用Apply.放大后的显示如下：放大后的显示如下：254STRATA 培训另一个需要提高的是改变色标范另一个需要提高的是改变色标范围围Color Key ，点，点 Data Range.改变值域范围如下图改变值域范围如下图.255STRATA 培训

131、最终反演结果如下显示：最终反演结果如下显示：256STRATA 培训约束稀疏脉冲反演参数约束稀疏脉冲反演参数:约束稀疏脉冲反演和基于模型反演用同样的参数约束稀疏脉冲反演和基于模型反演用同样的参数.这些额外的参数决定每道上将会检测到多少脉冲：这些额外的参数决定每道上将会检测到多少脉冲：Maximum Number of Spikes 最大脉冲数最大脉冲数Spike Detection Threshold 脉冲检测门槛脉冲检测门槛其其 它它 反反 演演 参参 数数最大似然稀疏脉冲反演参数最大似然稀疏脉冲反演参数257STRATA 培训最大脉冲数最大脉冲数 这个参数设置每道最多允许的脉冲数。缺省与时

132、间窗口内的总采样点数一样这个参数设置每道最多允许的脉冲数。缺省与时间窗口内的总采样点数一样. 实实际正常情况下际正常情况下,该参数不操作该参数不操作.脉冲检测门坎脉冲检测门坎Spike Detection Threshold每增加一个脉冲每增加一个脉冲, 它的振幅与它的振幅与已经检测到的所有脉冲的平均振幅进行比较已经检测到的所有脉冲的平均振幅进行比较. 当新振当新振幅比平均震幅给定的百分数小幅比平均震幅给定的百分数小, 软件将停止增加脉冲软件将停止增加脉冲.258STRATA 培训反演道功率谱反演道功率谱功率谱包含功率谱包含3个频带。个频带。Seismic Band是输入地震数据的频带，该信息

133、涉及到的低频和是输入地震数据的频带，该信息涉及到的低频和高频信息必须从初始模型或其它结合反演过程中预设信息中获取高频信息必须从初始模型或其它结合反演过程中预设信息中获取 。在反演菜单中的两个参数，在反演菜单中的两个参数，稀疏性稀疏性和和约束频率约束频率，用来控制这两个频带。，用来控制这两个频带。线性稀疏脉冲反演参数线性稀疏脉冲反演参数259STRATA 培训Constraint frequency = 10Constraint frequency = 20Constraint frequency = 1约束频带决定从初始模型中采用约束频带决定从初始模型中采用多少低频初始信多少低频初始信息息 。

134、 如果应用太多模型信息如果应用太多模型信息 (约束频带约束频带太大太大)，将出现，将出现一个与地震一个与地震拟合非常差拟合非常差的结果。的结果。 约束约束太小太小会导致会导致横向连续性变差横向连续性变差。约束频带约束频带260STRATA 培训Sparseness = 100Sparseness = 20Sparseness = 0.1稀疏性决定稀疏性决定多少高频信息多少高频信息要添加到模型中。要添加到模型中。 100% 就是添加所有高频信息。就是添加所有高频信息。 0% 就是不添加高频信息就是不添加高频信息,模型是带通的模型是带通的 。 稀疏性对与地震数据的拟合没影响。稀疏性对与地震数据的拟

135、合没影响。稀疏性稀疏性261STRATA 培训Window length = 128Window length = 256Window length时窗长度时窗长度线性稀疏脉冲反演是在一个时窗道内完线性稀疏脉冲反演是在一个时窗道内完成的。输入地震数据被分成多个重叠窗成的。输入地震数据被分成多个重叠窗口。每个窗口内完成反演并最终全部垒口。每个窗口内完成反演并最终全部垒加起来。时窗越小，程序运行越快，这加起来。时窗越小，程序运行越快，这是影响运行速度的主要参数。是影响运行速度的主要参数。262STRATA 培训带通反演菜单带通反演菜单:带通反演唯一的参数：带通反演唯一的参数：高截约束：这个参数控制

136、对初始模型的滤波以得到一个低频初始模型。初始模型中高截约束：这个参数控制对初始模型的滤波以得到一个低频初始模型。初始模型中大于这个值的所有频带被删除大于这个值的所有频带被删除 ，递推反演道中低于这个值的所有频带被删除。然，递推反演道中低于这个值的所有频带被删除。然后将这个低频部分与反演高频部分相加得到最终结果后将这个低频部分与反演高频部分相加得到最终结果 。带通反演参数带通反演参数263STRATA 培训不同的约束高截频设置对应的不同的最终初始模型：不同的约束高截频设置对应的不同的最终初始模型：264STRATA 培训首先，必须做频谱分析创首先，必须做频谱分析创建操作因子，结果显示在建操作因子

137、，结果显示在另一侧：另一侧：颜色反演一个关键参数是颜色反演一个关键参数是阻抗输出选项阻抗输出选项，这里，这里高频高频残差残差或者相对阻抗是默认或者相对阻抗是默认的。的。另一个选项是通过另一个选项是通过从模型从模型中定义频率范围创建全频中定义频率范围创建全频带带.颜色反演参数菜单颜色反演参数菜单265STRATA 培训频谱分析和操作因子创建结果有两部分，顶部图片显示了频谱分析和操作因子结果。频谱分析和操作因子创建结果有两部分，顶部图片显示了频谱分析和操作因子结果。颜色反演操作因子求取颜色反演操作因子求取井频率井频率井阻抗井阻抗地震频率地震频率能量能量操作因子频率操作因子频率O=I/S266STR

138、ATA 培训操作因子是用最小平方拟合法自动创建的。下面显示的部分允许你改变自动参数。操作因子是用最小平方拟合法自动创建的。下面显示的部分允许你改变自动参数。另一部分描述参数如菜单所示。另一部分描述参数如菜单所示。267STRATA 培训第一个菜单描述的是一系列井的阻第一个菜单描述的是一系列井的阻抗的不同频段的振幅谱。抗的不同频段的振幅谱。红色线是回归曲线，它代表颜色反红色线是回归曲线，它代表颜色反演的理想输出。演的理想输出。这些参数是允许你替换自动运算得到的这些参数是允许你替换自动运算得到的回归曲线。回归曲线。268STRATA 培训第二个菜单描述了两条曲线，一条第二个菜单描述了两条曲线，一条

139、是输入地震数据的振幅谱，另一条是输入地震数据的振幅谱，另一条是从先前的菜单中得到的理想输出。是从先前的菜单中得到的理想输出。注意这是一条新曲线，因为我们有注意这是一条新曲线，因为我们有一个频率域的线性比例。一个频率域的线性比例。这个参数允许你对地震频谱应用平滑。这个参数允许你对地震频谱应用平滑。Seismic SpectrumDesired Spectrum269STRATA 培训平行的红色是频谱门槛，它设置一个值平行的红色是频谱门槛，它设置一个值使反演算子在这个频率域进行，只有当使反演算子在这个频率域进行，只有当地震谱大于那个门槛的频段才参与运算。地震谱大于那个门槛的频段才参与运算。这个门槛

140、可以预防除这个门槛可以预防除0或小噪音值干扰。或小噪音值干扰。这个参数允许你改变门槛值这个参数允许你改变门槛值:Spectrum Threshold这些可选参数允许你手工设置频率范围这些可选参数允许你手工设置频率范围 (如果设置成非如果设置成非0值值) :270STRATA 培训最终，这些菜单显示了将要运算的时间域和频率域操作因子。最终，这些菜单显示了将要运算的时间域和频率域操作因子。菜单上其它参数发生变化，这些因子会自动更新来反映这些变化。菜单上其它参数发生变化，这些因子会自动更新来反映这些变化。这些参数直接反映了操作因子本身：这些参数直接反映了操作因子本身：271STRATA 培训来到这个

141、窗口（不是先来到这个窗口（不是先前的反演结果）点击前的反演结果）点击 Inversion / Post-stack Inversion / Bandlimited既然我们已经为既然我们已经为Blackfoot创建了初始模型，我们可以轻松的对数据应用这些反演创建了初始模型，我们可以轻松的对数据应用这些反演方法并将其与基于模型反演的结果进行对比。方法并将其与基于模型反演的结果进行对比。 首先从带通反演开始。首先从带通反演开始。练习练习 6: Blackfoot 其它反演方法其它反演方法272STRATA 培训Bandlimited Inversion在第一个菜单，改变输出在第一个菜单，改变输出体的

142、名字为体的名字为bandlimited, 如图所示如图所示.并且，对一条联络线进行并且，对一条联络线进行支演。如右图支演。如右图点击点击 Next 273STRATA 培训选择默认参数，主要参数如右图所示，选择默认参数，主要参数如右图所示，点击点击OK 得到最终结果得到最终结果:注意带通反演比基于注意带通反演比基于模型反演运算要快。模型反演运算要快。然而，结果少了些细然而，结果少了些细节信息。节信息。274STRATA 培训来到有地震数据的窗口，点击来到有地震数据的窗口，点击Inversion/Post-stack Inversion / Linear Programming Sparse S

143、pike在第一个菜单，改变输出体的名字在第一个菜单，改变输出体的名字 lp_inversion, 设置设置 Xline 范围范围, 如如右：右：Click on Next 线性稀疏脉冲反演线性稀疏脉冲反演275STRATA 培训应用默认选项，主要参应用默认选项，主要参数如右图所示：数如右图所示：276STRATA 培训这是最终结果，注意比较其与基于模型反演的细节部分类似。这是最终结果，注意比较其与基于模型反演的细节部分类似。277STRATA 培训回到含有地震数据模型的窗口，回到含有地震数据模型的窗口，点击点击Inversion / Post-stack Inversion / Colored

144、 Inversion在第一个菜单，改变输出体在第一个菜单，改变输出体的名字及设置联络线如右图的名字及设置联络线如右图所示。所示。点击点击 Next 颜色反演颜色反演278STRATA 培训应用默认参数。但是，我们需应用默认参数。但是，我们需要为反演创建操作因子，点击要为反演创建操作因子，点击Run Spectral Analysis ，如，如右图所示：右图所示：将会生成一系列图，在下页显将会生成一系列图，在下页显示。示。279STRATA 培训颜色反演操作因子颜色反演操作因子280STRATA 培训点击点击 OK，将会要求给出操作因子的名字，如左图所示。点，将会要求给出操作因子的名字，如左图所

145、示。点击击OK 接受默认选项，在接下来的菜单中采用默认选项，得接受默认选项，在接下来的菜单中采用默认选项，得到最终结果如下图所示。到最终结果如下图所示。颜色反演颜色反演281STRATA 培训将会生成一个相对阻抗体。现将会生成一个相对阻抗体。现在回到含有地震数据和模型的在回到含有地震数据和模型的界面重新运行界面重新运行 Colored Inversion, 改变名字改变名字colored_inversion_abs 并且改变输出阻抗体如右图所并且改变输出阻抗体如右图所示。示。282STRATA 培训应用同一个操作因子，采用默认参数，最终结果如下所示：应用同一个操作因子，采用默认参数，最终结果如

146、下所示：283STRATA 培训Model basedBandlimited这一页和接下来的这一页和接下来的2页会展示各种方法之页会展示各种方法之间效果的对比：间效果的对比：284STRATA 培训Sparse SpikeAbsolute colored285STRATA 培训Model basedBandlimitedSparse SpikeAbsolute colored(End of Exercise 6)286STRATA 培训总结总结反演反演 是从地震数据中提取隐含的地质信息。是从地震数据中提取隐含的地质信息。反演可以是多解的，低频模型尤其重要。反演可以是多解的，低频模型尤其重要。成功的反演信赖于每口井精细的对比和好的子波提取。成功的反演信赖于每口井精细的对比和好的子波提取。在这节课程中我们研究了以下类型的反演在这节课程中我们研究了以下类型的反演:基于模型反演基于模型反演带通反演带通反演颜色反演颜色反演稀疏脉冲反演稀疏脉冲反演AVO反演（只讨论）反演（只讨论）对大部分反演工作，基于模型反演是我们推荐的反演选择。对大部分反演工作，基于模型反演是我们推荐的反演选择。