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

Click to edit Master title style,Click to edit Master text styles,Second level,Third level,Fourth level,Fifth level,*,Last Updated:October 2004,*,STRATA 培训,STRATA,地震反演培训,理论和练习,2,STRATA 培训,STRATA,课程大纲,地震反演基础,-3,练习,1,：楔型模型简单练习（递推、模型）,-,61,练习,2,：,Erskine 3D,数据加载,-99,测井曲线时深对比,-118,练习,3,：,Erskine,3D,测井曲线时深对比,-131,测井曲线层间插值建模,-156,子波提取,-174,练习,4,：,Blackfoot,子波提取,-198,基于模型反演参数,-225,练习,5,：,Blackfoot,基于模型反演,-241,其它反演参数,-256,练习,6,：,Blackfoot,其它反演方法,-271,3,STRATA 培训,地震反演基础,反演,：从地震数据中提取它所包含的潜在地质信息的过程,.,传统上，反演是在叠后地震数据上进行,目的是提取,声波阻抗体,.,近年来，反演已经被扩展到叠前数据体,目的是既提取,声波阻抗又提取横波阻抗体,，从而可以推得孔隙流体。,另外，最新的发展认为可以利用反演结果直接预测岩性参数如,孔隙度和含水饱和度。,4,STRATA 培训,几种不同反演方法,递归法,:,传统带限反演,颜色反演,:,递归反演的改进,稀疏脉冲法,(,两个不同算法,):,约束以获得尽可能少的同相轴,基于模型法,:,迭代更新层状初始模型,弹性波阻抗,:,对,AVO,数据的改进,LMR:,对,AVO,数据的改进,联合反演：对,AVO,数据的改进,在,STRATA,软件中提供以下这些反演方法,(,六种方法七个模块,):,5,STRATA 培训,阻 抗,纵波,横波,弹性波,反射系数,子波,地震道,所有反演方法中共同的正演模型,:,反演中的一般正演模型,6,STRATA 培训,波阻抗,反射系数,反演,子波地震道,反演试图反推得到这个正演模型,:,反演模型,纵波、横波、弹性阻抗,7,STRATA 培训,纵波阻抗,或,横波阻抗,或,弹性阻抗,阻抗 反射系数,纵 波 阻 抗,=,横 波 阻 抗,=,弹 性 阻 抗,=,复杂公式,(,后面介绍,),8,STRATA 培训,反射系数 地震道,地震道,=,子波与反射系数的褶积,+,噪音,.,注释,不存在模型化的多次波,.,没考虑传播损失和几何扩散,.,没考虑频率吸收,.,子波可能是时变的,.,9,STRATA 培训,反射系数与子波褶积的结果是移除了大量高频细节,:,10,STRATA 培训,在时间域中的褶积就是频率域中的乘积,.,从这些图中可以看出,子波的作用是将地震频谱中高频和低频都消除了,.,理论上讲,反演就是试图将这些失去的频率区域进行恢复,.,11,STRATA 培训,反演中,“,多解性,”,问题,所有反演算法都有多解性问题,.,存在多于一种地质模型可以与地震数据相一致，要在这些可能的模型中决定一个模型，就需要去选择地震数据以外的的一些信息。,通常用以下两种方法来使用地震以外的信息：,初始猜测模型的建立,对最终结果与初始猜测模型背离幅度进行约束,最后反演结果既依赖于,“,其它信息,”,也取决于地震数据,.,12,STRATA 培训,带限,(,递归,),反演,递归反演,也称为带限反演是最简单和最早一种反演方法,.,从第一层开始,可以用这个公式递推后续每层的阻抗,.,:,i,th,+1,层的阻抗可以从上一层,i,th,确定,:,从定义反射系数开始,:,13,STRATA 培训,在这个简单的实例中可知,:,(a),告诉我们如果给出一个脉冲，那么我们就可以恢复阻抗的真实值,但是,(b),中看到如果我们用一个子波与脉冲进行褶积，那就无法恢复阻抗的真实值。,Z,2,=818,Z,2,=1500,Z,1,=1000,Z,3,=1227,Z,4,=1004,Z,1,=1000,带限反演,14,STRATA 培训,输入地震数据,递归反演数据,递归反演,获得的结果与输入的地震数据频带一样,.,从图中可以看出,与测井比较反演结果损失了高频细节,15,STRATA 培训,步骤,1:,递归反演,的初始背景模型是通过对井阻抗滤波获得,:,10-Hz,高截滤波,16,STRATA 培训,步骤,2:,对地震道运用递归运算,.(,注释,:,这几乎等同于在相位域进行,-90,度转换,):,17,STRATA 培训,步骤,3:,滤波后的模型加上,比例化,反演道集就得到了最终结果,:,+,=,18,STRATA 培训,递归反演中几个要点,:,不考虑子波,.,这意味输入的地震数据必须是零相位,.,如果有提取子波可行,Strata,软件可以自动地对地震数据,“,去相位,”.,即使地震数据是零相位,实际子波的旁瓣,在递归算法中也会被解释为岩性的变化,.反演的结果其频带与地震数据,频带一样,.,对,地震数据的加权值,(,与反射系数匹配,),决定了递归反演的结果阻抗值范围是否合适,.,19,STRATA 培训,基于模型的反演,基于模型反演,是从褶积模型方程式开始的,:,假定地震道,S,和子波,W,是已知的,.,假定噪音是随机的并与地震信号不相关,.,求解反射系数,R,以满足这个等式,.,这是一个非线性问题,所以求解过程是迭代进行的,.,20,STRATA 培训,步骤,1:,基于模型反演,的初始背景模型是通过对井阻抗分段形成的,:,用户用毫秒,(MS),来定义层分块大小,.,所有的层分段开始时都被设成一样大小,(,用毫秒,).,21,STRATA 培训,步骤,2:,利用层状模型和,已知,的子波计算人工合成地震道,.,这是与实际地震道的比较,.,通过分析人工合成地震记录与实际地震道的误差，每层,(,分段,),的厚度和震幅值都进行修改以减少误差,.,这个步骤通过一系列递归不断重复,.,Synthetic,Seismic,22,STRATA 培训,输入地震数据,基于模型反演,基于模型反演,获得了宽频的结果,.,问题是高频成分可能来自于初始猜测模型而不是来自地震数据,.,23,STRATA 培训,递归反演,基于模型反演,这是递归反演和基于模型反演的比较,.,通常基于模型反演可以获得更细致信息,但这个结果实际上相当类似,.,24,STRATA 培训,基于模型反演的几个要点,:,由于已经知道子波,在计算过程中它的影响从地震数据中排除了,.,也就是说地震数据没有必要是零相位,而只要,子波与地震数据等相位,就可以,.,估算的,子波若有误差,将导致反演结果出差错,.,地震有效,分辨率,得到提高,.,反演结果可能在很大程度上依赖原始猜测初始模型,.,解决的方法是对初始模型,进行滤波,.,与其它反演方法一样,存在,多解性,问题,.,25,STRATA 培训,稀疏脉冲反演,稀疏脉冲反演,假定实际反射可以认为是由一系列大脉冲里夹杂有小脉冲背景,.,稀疏脉冲反演,假定,只有大脉冲有意义,.,该方法通过检查地震道来寻找大脉冲的位置,.,26,STRATA 培训,稀疏脉冲反演,每次建立一个脉冲反射序列，,通过增加脉冲直到地震道被足够准确地反演完成。,阻抗块的振幅值是由,基于模型反演,算法来确定,27,STRATA 培训,输入地震,稀疏脉冲反演,稀疏脉冲反演,生成一个宽带高频反演结果,.,28,STRATA 培训,基于模型反演,稀疏脉冲反演,稀疏脉冲反演,结果与,基于模型反演,结果类似,.,主要区别是缺少非常薄层的细节,.,29,STRATA 培训,稀疏脉冲反演几个要点,:,只有当地震数据,有脉冲,，稀疏脉冲反演才得到同相轴；,它试图利用地震数据获得,最简单的可能模型,；,通常得到的结果比地质本身实际的,同相轴少,；,它比,基于模型反演,更,少的依赖于初始猜测模型,。,30,STRATA 培训,颜色反演,颜色反演,是对,递归反演,的修改,最早是由,BP,公司的,Lancaster and Whitcombe,在,2000,年,SEG,年会上提出来的,.,这种反演过程,寻找一个简单操作因子,O,对地震道进行直接转换求得反演结果,:,作者在频率域确定简单操作因子,O,.,通过比较地震数据和实际反演结果，认为简单操作因子的相位是,-90,度。,31,STRATA 培训,简单操作因子震幅谱用这种方式来求取,:,利用工区的一组井,所有井的声阻抗的振幅谱做交汇图,.,通过理论预测,我们可以拟合一条直线来 代表,“,理想,”,输出阻抗谱,.,Log(Frequency),Log(Impedance),声阻抗震幅谱,32,STRATA 培训,然后,利用井旁一组地震道,求取,平均地震谱,.,地震谱,Frequency(Hz),操作因子谱,从以上两个频谱,求取操作因子谱,.,使得地震谱形状在地震频带内转换为,阻抗谱,.,33,STRATA 培训,颜色反演操作因子,Time(ms),将获得的震幅谱做,-90,度相移就生成了,颜色反演操作因子,.,将它与所有地震道褶积,.,34,STRATA 培训,输入地震,颜色反演,颜色反演,结果非常类似于递归反演,.,其中一个区别是颜色反演的结果是相对声阻抗,既有正值也有负值,.,-3000,0,+3000,35,STRATA 培训,递归反演,颜色反演,-3000,0,+3000,4600,8300,12000,相对,AI,绝对,AI,36,STRATA 培训,颜色反演几点小结,除了定义整体阻抗趋势外,与初始模型基本无关,.,运算非常快,.,用户所需要定义的参数非常少,.,假定地震数据是零相位的,.,在最初实施中，该方法可以获得,相对,阻抗体，虽然我们可以选择添加低频阻抗背景趋势。,37,STRATA 培训,AVO,反演,基本褶积模型反演假定的是零偏移距地震数据,.,褶积反演法不能运用到带有,AVO,响应的地震数据,因为褶积反演法没有明确地考虑,V,P,/V,S,的变化,.,将反演方法延伸到叠前,AVO,地震数据的处理,目前这些算法主要运用于,:,弹性阻抗反演,Lambda-Mu-Rho(LMR),反演,联合反演,38,STRATA 培训,弹性阻抗,弹性阻抗概念最初由,Connolly,先生提出,(The Leading Edge,18,no.4,438-452(1999).,他利用,Aki-Richards,方程,将反射震幅与入射角联系起来,:,注释,:,在常规反演理论中假定,B=C=0,而没有考虑,V,P,/V,S,的变化,.,39,STRATA 培训,注意,对于零偏移距,:,通过类推,Connolly,定义了一种新阻抗类型,:,通过数学处理,最终,:,40,STRATA 培训这张图重叠显示了同一口井位置弹性阻抗与声波阻抗,.,在碳氢聚集区域弹性阻抗值显示异常低值,.,41,STRATA 培训,道集,AVO,分析,近角度叠加,stack at,q,1,远角度叠加,stack at,q,2,反演成弹性阻抗,EI(,q,1,),反演成弹性阻抗,EI(,q,2,),这种反演的工作流程是从迭前道集开始，首先产生两组角度叠加剖面,分别对其进行反演,.,42,STRATA 培训,这就获得两种反演结果,:,远角度反演,近角度反演,43,STRATA 培训,近角度反演与远角度反演进行交绘显示,.,远角度反演异常低值区可被解释为潜在含有碳氢聚合物,.,近角度反演,远角度反演,44,STRATA 培训,将交绘图中的异常低值点从交绘图中成像到原始地震数据剖面上以进行解释,.,45,STRATA 培训,Lambda-Mu-Rho(LMR),LMR,方法运用了,V,P,V,S,r,之间的关系式及拉梅常数,l,和,m,:,注意：最终结果表达式是,lr,和,mr,，可以通过纵波、横波阻抗,Z,P,和,Z,S,的关系式来表示。,LMR,方法最初是由,Goodway et al,(SEG Expanded Abstracts,1997),提出来的。类似于弹性阻抗法，这种方法对处理含,AVO,效应的数据比常规反演更有利。,46,STRATA 培训,Gathers,AVO AnalysisR,P,Estimate,R,S,Estimate,Cross-plot,Invert to Z,P,Invert to Z,S,Transform to,lr,and,mr,LMR,反演流程首先从叠前数据分析得到纵横波反射系数体,R,P,和,R,S,对,R,P,和,R,S,分别反演会生成纵横波阻抗体,Z,P,和,Z,S,.,这些数据体会通过类似,Goodway,等关系式或交汇图进行转换。,LMR,反演流程,47,STRATA 培训,如图展示了对,Albert,地区典型含气砂岩运用,LMR,法得到的结果，图中椭圆内指示为含气砂岩。,上剖面为,体,下剖面为,体,LMR,反演例子,48,STRATA 培训,左图为,与,交汇图，红色区域指示气层,(,低,),蓝色区域指示泥岩和水砂岩。,这些区域投到剖面上如下图所示，可以指示含气砂岩区域（红色部分）,lambda-rho,mu-rho,LMR,反演实例,49,STRATA 培训,联合反演,联合反演,可以对输入的叠前角道集数据反演得到,Z,P,Z,S,或者可能密度体。这种方法的优点在于它可以将这些参数进行约束，这样就可以使结果更稳定并且尽量避免多解性问题。,where:,我们从,Fattis version of the Aki-Richards,方程开始。这个公式用角度作为变量来表示反射系数：,50,STRATA 培训,联合反演,Ln(Zs),Ln(Zp),Ln(,),Ln(Zp),假定水砂岩可以用一个线性关系式来表示。,联合反演,寻找背景趋势的偏差,:,51,STRATA 培训,联合反演,联合反演,生成,Zp,、,Zs,、,Density,体,以及推导得到的各种组合。,52,STRATA 培训,联合反演,联合反演的解释与其它,AVO,反演结果类似：,53,STRATA 培训,它是地质统计反演的一种形式，可以通过对输入地震数据体运算产生大量的反演结果以此来解决多解性的问题。每一个运算结果要与地震数据一致，并且保持如方差图中所包含的所期望的连续条件。通过分析这些结果以估算得到一个非确定性结果，以及最大可能结果。,随机反演,Stochastic Inversion,This is a form of geostatistical inversion which explicitly addresses 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,uncertainty,in the result,along with the,most probable,result.,54,STRATA 培训,一般地震反演流程,(1),建模型,选择测井,每口井进行相关,提取子波,读取,/,拾取地震反演,(2),完成反演,选择反演方法和参数,对反演结果进行,QC,(3),解释反演结果,做数据切片,做交绘图,输出到,EMERGE,软件,55,STRATA 培训,反演结果,QC,两种方法,:,误差图,交互验证,我们如何知道反演结果是否有效,?,输入地震,反演结果,56,STRATA 培训,从阻抗道利用已知的子波可以计算人工合成道,.,理想情况,这个人工合成道应该与输入地震道非常相似,.,输入地震,人工合成道,57,STRATA 培训,从输入地震道中减去人工合成道,就得到,反演误差图,.,如果反演效果不错,误差图上应该只有非常小的震幅,并没有聚集在某处,.,由于存在多解性，即使微小的误差图也不能保证反演结果就是正确的,.,输入地震,反演误差,58,STRATA 培训,第二种质量控制方法是,交互验证,.,在这个验证中,我们在初始模型中删除一口井,在其位置做反演,然后将反演结果与隐藏井进行比较,.,隐藏井,反演结果,反演误差,59,STRATA 培训,每口井分析误差,以便识别有问题的井,.,60,STRATA 培训,反演作为,EMERGE,一种属性,最新使用反演结果是将其做为,EMERGE,软件的输入,EMERGE,可以直接预测孔隙度和其它岩性体,.,反演,EMERGE,孔隙度体,61,STRATA 培训,练习,1:,楔形模型练习,第一个练习对楔型模型做反演,.,目的是在简单例子上学习反演的基本步骤,.,在,WINDOWS(UNIX),上键入,geoview,或从,Windows,启动,GEOVIEW,软件,.,我们将创建一个新的,GEOVIEW,数据库,.,选择,New,如显示,然后在显示的第一个菜单上点击,OK,（如上图）,点击,Well Explorer,出现如下图所示界面。,62,STRATA 培训,现在我们为这个工区创建一个新的数据库,如果,Geoview,工区已打开可以点,Database/New,创建,如果这是,GEOVIEW,第一次在这台机器上运行，你可以看到如右图窗口，选择,New,并点击,OK.,63,STRATA 培训,将数据库命名,wedge_database,Ok,:,现在，可以看到,Well Explorer,中没有井显示：,64,STRATA 培训,现在，为这个数据库添加井曲线，如下步骤：,点,Import Data,:,选择,Wedge_well_logs.las,文件并点,Next,65,STRATA 培训,在接下来的菜单中，选择默认的即可：,Note that we are creating a new well called,Wedge_well,:,最后，我们可以读取声波曲线等如下：,在这个菜单,选择默认的参数即可，如,XY,坐标：,66,STRATA 培训,接下来点,OK,，出现如下提示窗口，问你选择哪种单位，,点击,Ok,即可：,67,STRATA 培训,当曲线加载完成后，出现如下窗口：,点“,1”,选择这口井,然后点,Display well,.,曲线就会显示出来。,68,STRATA 培训,接下来启动,Strata,，选择,GEOVIEW,主窗口上的,STRATA,按钮,.,选择,Start New Project,:,给新工区命名：,“,wedge_project”,:,69,STRATA 培训,在,STRATA,窗口上,通过,Data Manager/Import Data/Open Seismic/From SEG-Y File,加载地震数据,:,在,File Selection,菜单,选择文件,“wedge.sgy”,.,点击,Next,:,70,STRATA 培训,这菜单有几页来定义文件格式和数据的几何形状,(,观测系统,).,下一页,尽管数据是一条,2D,测线,我们可以建立,3D,观测系统,.,接下来,需要告诉程序在道头中没有,Inline&Xline,数,也没有,X&Y,坐标,.,71,STRATA 培训,点击,Next,两次,.,地震数据扫描后,出现这样界面,:,电击,Ok,利用这些缺省值,.,后面的练习将更仔细地观测分析地震数据的加载,.,72,STRATA 培训,现在可以看见这个楔型模型。,另外，还会有如下图所示窗口显示出来，定义井在地震上的位置，可以在,Xline,处填,45,，然后点击,OK,即可。,73,STRATA 培训,这口,井就会显示在这个模型上,然而，井曲线并不一定对应准确，所以需要做校正，点,Well/Edit/Correlate Well,:,让你选择井来进行校正，选择这口井，然后点,Edit,.,74,STRATA 培训,eLog,窗口就会出现，右侧显示的是模型数据,75,STRATA 培训,eLog,软件窗口弹出,(,显示井曲线,wedge log).,完成井曲线对比,在井口位置求一个平均地震复合道,.,缺省值是在井口,+/-1,周围道求平均地震道,.,点击,Ok,使用缺省值,.,进行井曲线对比，点击,Correlate,按钮,:,76,STRATA 培训,eLog,软件窗口显示如下,:,蓝色显示人工合成记录道,.,红色显示复合地震道,.,原始地震道用黑色显示,.,在对井曲线进行对比前需要一个更准确的子波,.,77,STRATA 培训,用于计算合成记录的子波是缺省的称为,“current wavelet”,.,查看缺省子波,点击,Wavelet/Display Current Wavelet,:,缺省子波是带通滤波,(,频率,:5/10-50/60.),点击,Frequency,标签可在频率域显示子波,.,78,STRATA 培训,求取,“,统计,”,子波（不使用井信息），生成零相位子波,其频率域与地震数据的响应一致，点击,Wavelet/Extract Wavelet/Statistical,:,在子波提取的界面,第一页是设置分析窗口,.,使用缺省值,.,点击,Next,两次,.,79,STRATA 培训,使用缺省值,.,这些参数后面再讨论,.,点击,OK,提取统计子波,.,新子波如左图所示,.,检查它的频率响应,点击,Frequency,标签,.,80,STRATA 培训,在窗口的底部子波名字也变了,.,井对比的过程是将蓝色的合成记录道上的点对应于红色地震复合道上的点,.,用鼠标左键点中，将四个同相轴连接,.,如右图所示点击,Stretch,按钮,:,新子波自动成为,“current wavelet”.,这意味它被用于在井对比窗口中合成记录的制作,.,注意到蓝色合成记录道有微小的变化,.,81,STRATA 培训,点击,Ok,(,在,Check Shot Analysis,窗口,),.,井对比类似于,check shot,校正,.,故出现两个新窗口,允许设置一些参数,.,后面更仔细谈论这些参数,.,82,STRATA 培训,点击,Wavelet/Extract Wavelet/Use Well,:,选择,wedge_well,点击,Next,使用缺省值点击,OK,(,在最后一页,).,这时你会发现合成记录拉伸后 与地震数据拟合的更好,.,现在可以利用井与地震之间的相关求取更好的子波,.,83,STRATA 培训,接下来都选择默认即可,这一页定义合成道在井位置处进行合成,选择时窗：,84,STRATA 培训,这一页定义子波提取方法以及基本的子波参数：,接下来的课程讨论一下这些参数,85,STRATA 培训,新子波看起来与上一个很相象,:,但是点击,Frequency,标签显示有一个非常小的相位校正,.,86,STRATA 培训,求取了子波，点击,Ok,查看新的对比结果,.,菜单告知新对比后的声波曲线将用名字,P-wave_corr,.,点击,Ok,:,点击,File/Exit Window,(,在,eLog,窗口,),:,87,STRATA 培训,现在加载地震层位用于初始猜测模型的内插。,在,Horizon Selection,菜单,使用缺省值名字点击,Ok,.,点击,No,(,在,dialog,窗口询问是否显示层位的平面图,).,点击,Horizon/Pick Horizons88,STRATA 培训,在,STRATA,窗口底部显示层位拾取的一些选项,:,89,STRATA 培训,如下图所示：,注意：一个选择参数“,Mode”,，我们选择,“Rubber Band”.,如果我们选择这项后，必须将剖面上目的层点亮，点下左键，释放鼠标后，层位会被拾取出来。,也可以选择“,Right and Left Repeat”.,我们选择用这种方式。,90,STRATA 培训,拾取第一层位,置鼠标在波峰,150 ms,周围,:,点击左键,.,同相轴自动左右拾取,:,T,开始拾取第二层位,点击,(,在,Horizon,下拉菜单,).,设置拾取模式,“Left&Right Repeat”,.,91,STRATA 培训,使用缺省名字,Horizon 2,.,然后改变,抓获,Snap,选项为波谷,Trough,:,点击楔型顶部的波谷,.,用类似方式,拾取楔型底部以及在,450 ms,周围的,波谷如图示,.,完成后点击,Ok,存储拾取的层位,.,92,STRATA 培训,现在可以建立初始猜测模型,利用井曲线和四个拾取的层位,.,点击,Model/Build/Rebuild a Model,(,在,STRATA,窗口,),:,第一页设置新模型的名字,.,利用缺省参数,参数含义后面课程再详细解释,.,连续点击,Next,最后点击,OK,建立模型,.,93,STRATA 培训,初始猜测模型生成，。,在,wedge,地震数据后面初始猜测阻抗模型用彩色显示,.,94,STRATA 培训,完成第一种反演,点击,Invert/Bandlimited,:,在反演菜单上,将输出名字,Output Volume File,改为,“,bandlimited,其它用缺省值,:,95,STRATA 培训,点击,OK,得到反演结果,:,96,STRATA 培训,现在进行第二种反演,(,在初始模型窗口,).,点击,Invert/Model Based,:,改变输出名字,Output Volume File,为,“,model_based,”:,使用其它缺省获得基于模型的反演结果,Model Based,:,97,STRATA 培训,鉴别反演结果好不好的一个重要标准是,“error display”,点击眼睛可以从反演窗口查看此项,如右图,在参数菜单选择,“model_based derived Synthetic Error”,如右图所示,点击,Ok,.,如右图,可以看到原始楔形模型与反演结果生成合成记录之间的差异,这个值接近于,0,指示反演结果是成功的。,98,STRATA 培训,在,Open Window Data Menu,显示了刚生成的所有数据体,:,该练习结束,.,点击,File/Exit Project,(,在,STRATA,软件任何窗口,),点击,Yes,(,在警告诊断窗口,),使,STRATA,关闭所有窗口并保存工区,.,99,STRATA 培训,练习,2,对加拿大西部的一块,3D,数据做反演,.,第一步是读一口井曲线到已经定义好的,GEOVIEW,数据库中,.,在做这个之前，先创建一个数据库,点击,Database/New,.,定义这个数据库名称,erskine_database,.,练习,2:Erskine 3D,初始化,100,STRATA 培训,在,GEOVIEW Well Explorer,窗口,点击,Import Data/Logs,Check Shots,Tops,Deviated Geometry from Files,:,101,STRATA 培训,改变,Destination Well Name,域为,“,erskine_well,”,点击,Next,:,在,File Import,页面,选择文件,e,“erskine_log.las”,点击,Next,.,102,STRATA 培训,显示有一口井声波曲线,.,点击,Ok,读取这根曲线,.,点击,Next,(,利用缺省位置信息,),选取默认单位,103,STRATA 培训,在井曲线读进,GEOVIEW,后,点击,erskine_well,然后点击,Display Well,:,erskine,声波显示如右图,.,104,STRATA 培训,在,GEOVIEW,主窗口上点击STRATA,启动,STRATA,软件,.,定义工区为“,erskine_project,”,选择,Start New Project,:,在,STRATA,建立新工区对,erskine,地震数据进行反演,.,105,STRATA 培训,第一步将,3D,地震数据读进,STRATA,软件,.,点击,Data Manager/Import Data/Open Seismic,/From SEG-Y File,:,选择文件“,ersk3d.sgy,”,点击,Next,:,106,STRATA 培训,同样做为三维地震数据体加载地震数据,.,这个文件在道头中没有,Inline&Xline,号或,X&Y坐标,.,设置,primary key to“,CDP(Rectangular),”.,这是指数据体是由一系列测线组成,每条线有相同的道数,.,107,STRATA 培训,改变观测系统键入,“,155,”,如图显示,.INLINE,线将由程序计算,.,点击,Ok,加载地震数据体,.,点击,Next,两次选择缺省直到最后一页,.,软件最初假定只有单一一条,INLINE,线,.,108,STRATA 培训,在,Well To Seismic Map,菜单,如下所示将井设置在地震数据上,.,在,Plot,box,中,wedge_well,井没有被选用因为它不属于这个数据体,:,在地震显示中,键入,24,”,作为,Inline,号显示的结果如右图所示,:,109,STRATA 培训,对这个数据体要建立初始猜测模型,需要一组层位文件,.,首先拾取一个层位然后读进以前拾取的一组层位,.,点击,Horizon/Pick Horizons,:,使用缺省名字,“,Horizon 1,”.,点击,Ok,.,点击,“,Yes,”(,诊断窗口询问是否显示层位平面图,).,110,STRATA 培训,向以前一样,将鼠标点在靠近层位处拾取单个层位,:,平面图窗口显示这条,INLINE,线上拾取的时间,.,111,STRATA 培训,软件利用这个人工拾取层作为引导可以自动拾取整个,3D,数据体层位,.,点击,Options/Automatic Picking,:,点击,Ok,(,在,Automatic Picking,菜单,),整个,3D,层位就自动拾取,.,在层位的平面图上观测,最左边上的几条,INLINE,线上可能拾取的层位位置有问题,.,112,STRATA 培训,显示,inline,2,测线,(,在,STRATA,窗口,).,用手动方式在同相轴上周围重新拾取层位,:,我们发现自动拾取的结果在一个小区域上层位跳出一段,.,113,STRATA 培训,现在再重新自动拾取,.,在自动拾取菜单上发现缺省参数是清除以前自动拾取的层位,而保留人工拾取的层位作为新引导,:,点击,Ok,得到新层位结果人工拾取第二条,INLINE,线后模型得到了改进去引导自动拾取,:,114,STRATA 培训,下面我们读进层位,.,首先删除刚拾取的层位,:,层位,Ok,删除拾取的层位,.,点击,Horizon/Import Horizons/From File,:,选择五个文件从,erskine1.pik,到,erskine5.pik,.,115,STRATA 培训,文件类型是缺省,Geoquest,:,点击,Next,选择缺省参数包括名字和层位颜色,.,点击,OK,读进层位,.,116,STRATA 培训,层位读进后,再显示,inline,24,STRATA,窗口显示为,:,117,STRATA 培训,点击,Model/Build/Rebuild a Model,:,在建模型菜单,我们修改其中两个参数,其它用缺省参数,:,现在建立反演所需的初始模型,.,118,STRATA 培训,每个地震道初始猜测模型是由,阻抗曲线组成,的，阻抗曲线通常是将实际声波曲线乘以实际的密度曲线。阻抗曲线模型是以双程旅行时来计算的，原始的测井曲线是按深度测量的，关键的一步是作,深,-,时转换,：,测井曲线时深对比,119,STRATA 培训,深,-,时转换是用深,-,时表来完成，转换过程是将每个深度点转到双程旅行时,(,从基准面,(,地表,),到那个深度和从那个深度回到基准面的时间,),：,120,STRATA 培训,到达一个同相轴的时间取决于这个层以上的所有速度,包括到地表的第一层的速度,V,1,，这个速度是未知的,它通常是由测量的第一个速度近似地外推到地表：,深,-,时表用这个等式从声波曲线速度中计算得来,:,其中,:t,i,=,到达第,i,层的时间,d,j,=j,层的厚度,V,j,=j,层的速度,121,STRATA 培训,如果井是斜井,需要校正,(,投影,),成直井并从,KB,校到基准面,:,D,M,=,从,KB,计算的测量深度,D,V,=,从,KB,计算的垂直深度,D,S,=,从基准面计算的垂直深度,T,=,从基准面计算的双程旅行时,D,M,D,V,D,S,122,STRATA 培训,从声波曲线计算的深,-,时表很少能获得理想的阻抗模型使之能与地震数据适当的匹配，原因是：,地震,基准面,与井曲线基准面可能不一样,.,第一层的平均速度是未知的,.,声波速度的误差会造成计算的传播旅行时有,累积误差,.,地震上的同相轴由于地震,偏移处理,误差会造成,时间偏移,.,地震数据由于频率吸收和短周期多次波等影响可能带有,时间拉伸,123,STRATA 培训,改进深,-,时表可使用两种途径,:,运用,check shot,校正,.,应用手动,(,交互图形化,),法将井曲线与地震数据进行对比,Check Shot,校正,一个,check shot,表是由一组深度点计算的实际双程旅行时组成的,:124,STRATA 培训,从声波曲线计算的深,-,时表需要修改以获得理想的,check shot,时间,:,原始深,/,时,曲线,理想深,/,时曲线,125,STRATA 培训,线性拟合,:,用直线段内插,样条拟合,:,用平滑曲线内插,多项式拟合,:,利用最小平方拟合一条光滑曲线,漂移曲线两点间的内插方法有三种选项,:,126,STRATA 培训,改变深,-,时表意味着原始声波曲线速度可能发生变化,.,在,STRATA,软件有三种选项,:,(1),井曲线的,所有速度都改变,使新井曲线可获得理想对应的双程旅行时,.,注释,:,这意味着在第一个测量深度前用,斜波速度,来处理整体时移以最小化合成记录中可能的伪反射。,在,STRATA,称此选项为,“,Apply All Changes,”,。,127,STRATA 培训,(2),只改变,check shot,中起始深度和最后深度之间的层速度。,这是指在第一个测量深度前不加斜坡速度,.,获得新井曲线可得理想对应时间但不考虑整体时移,在,STRATA,称为,“,Apply Relative Changes,”.,128,STRATA 培训,(3),不改变声波曲线的速度,.,新井曲线不能对应理想的双程旅行时间，但,GEOVIEW,和,STRATA,将利用这个新的深,-,时表,.,这个选项可以保留原始反射系数进行人工合成记录的计算,.,在,STRATA,称为,“,Change Depth-Time Table Only,”.,129,STRATA 培训,注释,:,例子中时间的拉伸大的不现实,.,不同的内插选项可能使声波曲线的变化相当大,:,130,STRATA 培训,对比就是对深,-,时曲线进行,手工校正,以优化初始模型和地震数据间的相关,.,对比应在,check shot,校正之后进行,理想情况时深校正对井的,改变很小,的,.,测井曲线对比改变深,-,时曲线的方式与,check shot,校正一样,.,测井曲线对比是在合成记录上,选择几个同相轴,使之对应于实际地震道上,.,由于使用合成记录,子波的选择,是至关重要的,.,测井曲线