文档介绍：该【地震反演方法综述 】是由【1781111****】上传分享，文档一共【9】页，该文档可以免费在线阅读，需要了解更多关于【地震反演方法综述 】的内容，可以使用淘豆网的站内搜索功能，选择自己适合的文档，以下文字是截取该文章内的部分文字，如需要获得完整电子版，请下载此文档到您的设备，方便您编辑和打印。在上世纪70~80年代,地震反演作为地球物理学的一个重要进展得到了广泛的赞扬,获得广泛应用;地震反演技术能够帮助解释人员确定地层单元而不仅仅是通过反射波确定地层单元的边界,而且能直接进行深度域成图。在一个竞争的市场环境中,开发出了很多不同的反演算法,在基本递归反演方法的基础上不断取得进进展,一下简要介绍几种基本的地震反演方法。主要分三大类:1、基于地震数据的声波阻抗反演:其结果有两种:相对阻抗反演(常说的道积分)与绝对阻抗反演。主要算法有:递归反演(早期的地震反演算法)与约束稀疏脉冲反演(优化的地震反演算法)。这种反演受初始模型的影响小,忠实于地震数据,反映储层的横向变化可靠;但分辨率有限,无法识别10米以下的薄砂层。2、基于模型的测井属性反演:此种反演可以得到多种测井属性的反演结果,分辨率较高(可识别2-6米的薄层砂岩);但受初始模型的影响严重,存在多解性,只有井数多(工区内至少有10口以上的井,分布合理,且要求反演的属性与阻抗相关),才能得到较好的结果。3、基于地质统计的随机模拟与随机反演:此种算法可以进行各种测井属性的模拟与岩性模拟,分辨率高(可识别2-6米的薄层砂岩),能较好的反映储层的非均质性,受初始模型的影响小,在井点处忠实于井数据,在井间忠实于地震数据的横向变化,最终得到多个等概率的随机模拟结果;但要求工区内至少有6-7口井,且分布较合理,才能得到好的模拟结果。道积分道积分技术出现,为广大少井无井地区岩性及油气预测提供了新的途径,它能得到类似于虚速度测井的新方法,其结果对应于地层的波阻抗,它最大优点是不像虚速度测井那样依赖于井的资料和地球物理学家的经验。尽管道积分剖面不能像GLOG波阻抗剖面那样反映地层绝对速度,而只能反映其相对速度大小,但是它反映出的层位与GLOG剖面是一样的,甚至在反映的细节上还比它多,对薄层识别也非常有利,因此道积分剖面能用于岩性和油气层解释。道积分的基本原理:式中V为地震波速度,K为常数,*表示褶积,Fn为零相位带通滤波因子,Xi为地震道。该式反映了相对速度的变化,在处理上,它是将叠偏后的地震剖面的何一道分别进行积分(振幅累加),边积分边输出,便得到与地震剖面相对应的相对速度剖面,即道积分剖面。道积分剖面,对应于地层的相对波阻抗剖面,当地层进入一定深度后,不同岩性密度()大致相等,也可视为地层的相对速度剖面。因此它非常有利于地层岩性解释。在正极性道积分剖面中,“正波”(正相位)代表相对高波阻抗层,即相对高速度、高密度地层,“负波”(负相位)代表相对低波阻抗层,即相对低速度、低密度地层,正、负波的极值点(峰值),分别对应地层中点,零相位则指示了波阻抗界面或近似的地层分界面,因此当我们知道了地区岩相和油气的波阻抗特征后便可用道积分剖面进行油气和岩层解释。行道积分剖面的标定解释。分析可知,堤北地区馆陶组上段砂岩速度高于泥岩速度,而且从现有井资料看,该地层为砂泥岩互层结构。因此,在道积分剖面上波峰同相轴应当反映砂层,波谷同相轴应反映泥岩层或气层。但这并不是绝对的。统计资料表明,上馆陶组地层在道积分剖面上有90%以上的波峰同相轴与井中砂层相对应。但也有例外,如当砂层很厚时,道积分剖面不是表现为一个宽波峰,而是表现出顶部为一条强同相轴,其下为弱波峰波谷的组合形式(即空白带)。这里以CB14井及通过该井的645。1剖面为例,来说明道积分剖面的解释方法。1将声波曲线进行积分,转换成等时间序列的层速度曲线(即CVL转换曲线)。图1是CB14井的CVL转换曲线,图中左边的数字为深度值(间隔不等),右边的数字为时间(等间隔)。图1CB14井CVL转换曲线式中,R2和N2分别表示反射系数和噪音的均方值;r(K)和n(K)分别表示第个采样点的反射系数和噪音;M表示反射层数,L表示采样总数,表示给定反射系数的似然值。通过多次迭代,求取反射系数。。如果从最大似然反褶积中求得的反射系数是R(t),则波阻抗在上述过程中为了得到可靠的反射系数估计值,可以单独输入波阻抗信息作为约束条件,从而求得最合理的波阻抗模型。,一是针对单口井存在的问题而进行校正:①在储层井段,测井资料受泥饼的影响会出现测量误差;②在井壁不规则时,测量结果包含泥浆浸泡因素导致测量误差。二是测井资料通常由不同仪器、不同年代和不同处理流程得到,虽然其声波和密度值的相对关系明确,但有时其基值和相对幅值相差很大,不适合于区域的波阻抗建模,经过匹配处理和资料校正从而获得比较符合地质特征的波阻抗模型。地震资料处理主要是提高原纯波地震资料的分辨率。反演的分辨率主要取决于所用地震资料的分辨率。影响地震资料分辨率的因素,一是地震资料的采样率,二是主频高低,三是频带宽窄。提高地震资料分辨率的措施是提高采样率和主频,拓展频带。,子波提取是关键。首先是初始子波选择。一般选用零相位Ricker子波作为初始子波,在子波频率选择上,通常是选择与地震资料主频相近的子波频率。其次,通过循环反馈不断逼近最终得到最佳子波。从每口标定井提取子波,并利用该子波重新对合成记录进行微调,再重新提取子波,将每口精细标定过的井提取子波进行比较,子波的相位、振幅比较一致,最终得到最佳子波。同时,在反演面积较大的情况下,子波特征会有相应的微小差异,与地层物性特征或速度特征变化有关。针对这种情况,适当采用随空间特征变化的空变子波,有助于对地震资料及地下地层特征的正确反演处理,得到的反演结果也会更加接近地下岩层的真实地质特征。图2表示了空变子波体的形成过程。。对地震来说,关键是对正确解释起控制作用的波阻抗界面;对测井而言,主要是对地震波阻抗界面赋予合适的波阻抗信息。建立合理的初始模型,要注意区分地质模型和地球物理模型的差别。地质模型表达的是地质体的空间存在,用精细的点描述地质体的空间分布,点的精度可以用分米记;地球物理模型表达的是地质体的空间实在,用测量的结果寻找地质体的空间分布,一般地,整体精度只能达到十米以上。波阻抗初始模型本身是地球物理模型,其精度上限是地震资料分辨率,而测井资料和层位、断层等约束条件构建的是地质模型。建立初始模型过程中,两者要做到统一。具体方法是:在精细制作合成记录的基础上,精细解释层位和断层,作多井过井剖面,分析区域的沉积体系,从而建立地层构造模型,再在此基础上,形成合理的初始波阻抗模型。再利用校正后的测井曲线作为约束,从而形成约束波阻抗模型。这样就避免了强行将测井曲线内插、注意地质模型的精细性而忽略地质背景的低频特征以及地球物理模型的实际分辨率限制。:min∑ri+(di-si),式中:ri为反射系数;di为地震数据;si为合成道数据;di-si为地震残差;λ为控制参数。因此,反射系数的稀疏、合成记录与原始地震道的残差最小是相互矛盾的。反演参数λ控制反演脉冲的稀疏程度,即对反演起着协调低频一致性和高频约束性的作用。通过调整反演参数λ,可以使反演剖面在保持细节和不损失低频背景中取得平衡,兼顾纵向和横向分辨率,使其取得最佳效果。具体调试过程,主要是通过信噪比、测井曲线相关性、地震道不匹配度、反射系数不匹配度几个指标控制来实现。,因此反演结果缺少地震数据所缺失的低频信息,需要进行低频补偿。实现方法是:将参与反演的测井数据进行内插,得到一个阻抗体;然后设计一个低通滤波器对测井曲线进行滤波,将滤波后的结果与约束稀最终得到一个全频带的绝对波阻抗数据体。实际操作中,波阻抗低频分量使用了多井加层位控制的低频分量数据体,实现了低频分量的空变,使波阻抗数据体更加准确,与井资料有良好的对比关系。基于模型的测井属性反演测井约束反演或宽带约束反演(Broad-bandConstrainedInversion,简称为BCI),就是一种基于模型的反演。其基本思想是:测井资料有很高的垂向分辨率,但只是点上的一孔之见;地震勘探的垂向分辨率虽不高,却具有线上和面上的详细资料。将两者结合起来,取长补短,在垂向上充分利用井的高分辨率信息,在横向上充分利用地震资料的可对比性作为控制,建立起较为可靠、分辨率较高的初始地质模型;对初始模型进行正演,计算出合成地震剖面,再与实际地震剖面相比较,求取速度、密度及深度等模型修改参数(摄动量);反复修改更新初始模型,使合成地震剖面与实际地震剖面在最小平方意义下最为接近,最终得出高分辨率的波阻坑反演剖面。其处理流程如图3所示。反演所用的公式为图3测井约束反演处理技术流程式中,M为更新的模型;M为初始模型;G为灵敏度矩阵,或称雅可比算子,它是由0一系列偏导数组的矩阵;G为噪音协方差矩阵;G为模型协方差矩阵;I为单位矩阵;Snm是地震数据;D是计算的地震数据;S一D称剩余偏差或残差;M-M为模修改量或称摄动0量。每次模型修改后,再重复以上计算,直至残差小到可以接受的程度,中止迭代,即可得到相应的波阻抗模型。基于地质统计的随机模拟与随机反演地质统计学反演(GeostatisticalInversion)又称随机反演(StochasticInversion),它结合了地震反演储层预测方法和随机建模储层预测方法的优势,充分利用随机建模技术综合不同尺度的数据(如模拟过程可综合构造和地层信息),结果的多个实现可用于不确定性评价。目前该方法在国外已广泛应用于油气田开发阶段的储层预测,取得了很好效果。1地质统计学反演总体思路地质统计学反演算法由两个部分组成:①序贯随机模拟;②优化随机模拟结果。通过地质统计学反演使随机模拟结果符合地震数据,使最终反演结果同时符合井数据、地震数据以及已知信息的地质统计学特征,从而建立与三维地震数据CDP面元大小一致的三维网格,井间受地震数据约束,同时以地震反射层解释数据文件和高分辨层序地层学分析结果建立精细构造和地层模型,决定数据内插的空间位置。地质统计学反演中用到克里金技术,而克里金技术要求利用测井资料计算垂向变差函数,利用地震数据计算水平方向变差函数。地质统计学反演算法的基本思路是:首先随机选取一个点,并用普通克里金技术计算该点的局部概率密度函数,通过序贯高斯随机模拟建立井间波阻抗,然后将波阻抗转换成反射系数,并利用由确定性反演方法求得的子波与反射系数褶积产生地震道,通过反复迭代直至合成地震道与原始地震数据达到一定程度的匹配为止;按照上述方式对逐个网格点进行模拟和优化,反演结果是多个等概率的波阻抗数据体。2技术关键和基本流程地质统计学反演过程中的关键问题是如何使下一个模拟值(从pdf中抽取值)比前一个模拟值更快地达到规定的与地震数据的匹配程度,以避免陷入局部极小或由于大量的迭代次数而耗费机时。目前多采用模拟退火算法解决这个问题。地质统计学反演的基本流程:(1)测井数据预处理首先对声波和密度测井数据进行环境校正,以消除泥浆浸泡和井眼因素对测井数据的影响;在此基础上进行测井数据标准化处理(即在目的层附近选取一个或多个标准层利用趋势面分析方法进行校正);最后将预处理后的声波和密度曲线合成波阻抗曲线。(2)测井波阻抗数据体的统计分析由于序贯高斯模拟要求数据服从高斯分布,因此模拟前应对数据进行分析,若不服从高斯分布需要进行数据转换。(3)建立初始模型由于模拟退火算法并不依赖初始模型,因此可利用确定性反演结果作为初始模型。(4)拟合变差函数求取空间变差函数是地质统计学反演的关键步骤之一,在常规随机建模中一般使用井点数据,但井点数据只在垂向上数据点密集,求取变差函数精度高,而在水平方向可供用于拟合的数据点稀少,点距大,拟合出的变差函数精度低(图1)。在地质统计学反演中可利用井数据求取垂向变差函数,而水平方向的变差函数可利用稀疏脉冲约束反演波阻抗数据体计算(图2),以充分发挥地震波阻抗数据体在横向上密集的优势。通过对各个方向进行扫描计算,求取各个方向的变差函数,充分表征储层性质空间分布的各向异性。由图1、图2可以看出:①在垂向上波阻抗变程很小,这说明波阻抗在垂向上变化快,在水平方向上,不同方向变差函数变程不同,反映了不同方向数据相关性的差异;②图2的水平方向变差函数的精度明显高于图1垂向变差函数。因此可充分发挥井点数据和稀疏脉冲约束反演波阻抗数据体在垂向和水平方向的分辨率优势,通过实验选取合适的变差函数模型拟合三维变差函数。变差函数计算可在地层模型控制下分层进行,然后分别赋予相应的层。由于地质模型是在地震资料精细解释和高分辨层序地层学分析的基础上建立的,从而保证了层位对比方案的精确性。1根据井资料求取变差函数图2根据地震波阻抗求取变差函数(5)确定控制参数对于模拟退火算法而言,冷却进度表包括4个参数:①起始温度一般要求足够大,以保证能跳出局部极小,但过大的初始温度会无端增大计算量;②温度下降参数该参数主要是保证让温度能缓慢下降,以防止淬火现象。温度下降函数多种多样,常用的是一种直观的线性下降函数Tk+1=Tk,0<α<1(α越接近于1,温度下降越慢,温度下降太慢将增加计算量;下降太快可能会错过最优解);③终止准则确定给定一个接近于0的很小值作为计算终止条件,以保证计算能收敛到最优解;④循环迭代次数从理论上讲,迭代次数越多系统越能趋于热平稳,合成地震数据愈能和原始地震数据达到最佳匹配,但这样将增加计算工作量。