地震波阻抗反演方法综述

地震波阻抗反演方法综述、地震反演技术研究现状地震反演方法是一门综合运用数学、物理、计算机科学等学科发展起来的新技术新方法，每当数学方法、物理理论有了新的认识和发展时，就会有新的地震反演技术、方法的提出。随着计算机技术的不断发展、硬件设施的不断升级，这些方法技术得到了实践验证和提升，反过来地震反演技术运用中出现的新问题、新思路又不断促使数学方法、地球物理学理论的再次发展。时至今日，地震反演技术仍然是一个不断发展、不断成熟、不断丰富着的领域。反演是正演的逆过程，在地震勘探中正演是已知地下的地质构造情况、岩性物性分布情况，根据地震波传播规律和适当的数学计算方法模拟地震波在地下传播以及接收地震波传输到地表信息的过程。地球物理反演就是使用已知的地震波传播规律和计算方法，将地表接收到的地震数据通过逆向运算，预测地下构造情况、岩性物性分布情况的过程。地震波阻抗正演是对反演的理论基础和实现手段。1959年美国人EdwinLaurentineDrake在宾夕法尼亚州开凿的第一口钻井揭开了世界石油工业的序幕。从刚开始的查看地质露头、寻找构造高点寻找石油，到通过地震剖面的亮点技术寻找石油，再到现在运用多种科学技术手段进行油气资源的预测，石油勘探经历了一个飞速的发展历程。声波阻抗（AI）是介质密度和波在介质中传播速度的乘积，它能够反映地下地质的岩性信息。声波阻抗反演技术是20世纪70年代加拿大RoyLindseth博士提出的，通过反演能够将反映地层界面信息的地震数据变为反映岩性变化的波阻抗（或速度）信息。由于波阻抗与地下岩石的密度、速度等信息紧密联系，又可以直接与已知地质、钻井测井信息对比，因此广泛应用于储层的预测和油藏描述中，深受石油工作者的喜爱。70年代后期，从地震道提取声波资料的合成声波技术得到了快速发展，以此为基础发展的基于模型的一维有井波阻抗反演技术，提高了反演结果的可靠性。进入80年代，Cooke等人将数学中的广义线性方法运用于地震资料反演，提出了广义线性地震反演。此后Seymour等人又提出了测井声波资料和地震数据正反演相结合求取地下声波阻抗的测井约束反演，大大拓宽了反演结果的纵向分辨能力。90年代，在基于前人对地质统计学研究的基础上Bortoli和Haas提出了地质统计学反演，Dubrule等人对该方法进行了改进和推广。 在国内随着油田对地震反演技术的广泛应用，以周竹生为主提出的地震、地质和测井资料联合反演方法，将地质信息引入地震反演中，提高的反演结果与地质认识的联系，克服了线性反演存在的缺陷。 1996年，李宏兵等人将宽频带约束方法应用于递推反演并对其进行改进，减弱了噪音对反演结果的影响。1999年，任职于英国石油公司的Connolly在《弹性波阻抗》一文中介绍了弹性波阻抗（EI）的概念和计算方法，阐述了不同入射角度（偏移距）地震道集部分叠加反演波阻抗随入射角之间的关系，但是该方法求取的弹性阻抗随入射角变化很大，无法与常规叠后反演波阻抗直接比较，因此推广应用较为困难。2002年，Whitcombe通过修正PatrickConnolly的计算公式，得到了弹性波阻抗的归一化求取方法，消除了弹性阻抗随入射角变化大的难题。2003年，西北大学马劲风教授从Zoeppritz方程简化出发提出了广义弹性波阻抗的概念，克服了以往波阻抗反演要求地震波垂直入射到地表的假设条件， 推导出了任意入射角下纵波反射系数的递推公式，提高了中等入射角度下弹性波阻抗反演的精度。

、地震反演基本原理正问题和反问题就像物理学中的作用力和反作用力一样是同时存在的。 地球物理学中的正问题一般是已知我们容易直接测量或者获取的物理量， 通过两者内在关系求解未知的不容易获取物理量问题，反问题是已知我们不容易获取的物理量求解问题应该满足的条件或物理量。求解正问题的过程称为正演，求解反问题的过程称为反演。一般情况下正演相对简单，求解结果唯一，而反演的求解大多是欠定的，求解过程较为复杂，其结果也可能出现多解性，因此减少反演结果多解性是反演问题的一个重要研究内容。在油气地震勘探中，当地下地质体的岩性、物性发生变化时，经常会引起岩石密度和地震波传播速度的变化，导致波阻抗产生变化，波阻抗在岩性、物性界面处的变化会生成一个反射系数脉冲，当地震波传播到波阻抗界面时，就会在其界面处产生反射波。因此，若已知地下介质的波阻抗，通过求取反射系数与并地震子波褶积运算得到地震波通过地下介质后观测记录的过程就是地震波阻抗正演。 当我们根据观测到的地震记录， 通过地震波反射的褶积理论等求取反映地下介质的岩性、 物性信息的波阻抗等信息的过程就是地震波阻抗反演。 广义的波阻抗反演包括地震资料的处理、 解释以及波阻抗求取，一般讨论的波阻抗反演是在地震资料处理和解释的基础上通过反演方法对波阻抗的求取。可以直观的进地震资料是反射界面的表现，是一种界面型信息，只能间接表达地下地层的地质特征，不能直接进行储层的描述。而作为储层预测重要工具的地震反演技术可以将地震数据转换成波阻抗数据，波阻抗数据由于是地层型剖面， 可以和测井地质信息直接对比，可以直观的进行储层的识别和描述。按反演结果分

表2-1按反演结果分

表2-1是四种分类、按测井在地震反演中的作用大小来分和按反演方法的实现方式来分等。类方法对应的主要反演方法。表2-1地震反演方法分类表反演分类方法反演方法按使用地震资料分类叠后反演叠前反演按反演实现思路上分类地震道直接反演模型反演随机反演按反演结果分类构造反演波阻抗反演岩性反演物性参数反演按测井资料在反演中起到的作用分类无井地震反演井控地震反演测井-地震联合反演地震控制下的测井反演

2.1直接反演直接反演方法是在反褶积理论的基础上， 通过对地震道进行数学计算求取地下波阻抗体的反演方法。假设地震记录褶积模型为(2-1)S(t)=W(t)*R(t)+N(t)(2-1)式中，S(t)为地震记录，W(t)为地下介质中传播的地震子波， R(t)为地下波阻抗界面的反射系数，N(t)射系数，N(t)为采集等原因产生的噪音,般认为是白化的。这就是地震记录的褶积模型,由模型可以看出，地震记录由地下波阻抗界面反射系数与地震子波的褶积加上一定的白化噪音构成。地震道直接反演是根据地震波传播理论，用数学计算工具消除地震记录中的子波W(t)、弱化白化噪音N(t)，反射界面的反射系数序列 R(t)，再通过积分或递推等方法求取地下波阻抗。道积分反演和递推反演是最常见的直接反演方法。2.1.1道积分该方法在地下地层波阻抗是深度变量连续且可再对反射系数进行深度积该方法在地下地层波阻抗是深度变量连续且可再对反射系数进行深度积反演波阻抗数据体是微函数的假设条件下，通过地震道自相关方法求取地层反射系数，分计算出地层波阻抗。由于道积分反演实际上只是对反射系数的积分,相对值。地震波垂直入射到波阻抗界面时，反射系数的表达式为R_E十R_E十Vi十~EViPm十+Pvi Aid+AliiVi =Aid—AIi(2-2)V十分别是第i、i+1层的式中Ri为第i、i+1V十分别是第i、i+1层的密度和速度，AIi、Alj十为第i、i+1层的波阻抗。根据地层波阻抗是深度连续可微函数的假设，可以将地层看成是很多很薄的地层，因此相邻地层的波阻抗差异不会太大AIi十止AIi，所以AIh^-AI^Ml、Ali++Ali俺2AI，公式(2-2)可以改写成2AI(2-3)2AI对反射系数在深度t上进行积分，可得地层相对波阻抗 AI(t)tJ。RdttJ。Rdt屯MIdt=丄[1nAl(t)-lnAI(0)]2AI2(2-4)2Rdt(2-5)Al(t)=AI(O)eb(2-5)式中R为地层反射系数，t为地层深度，AI(t)为地层波阻抗。反演过程具有积分误差地震处理等原因，缺失低频和道积分反演方法无反演过程具有积分误差地震处理等原因，缺失低频和道积分反演方法无仅适用于地震勘探初期所以道积分反演不能高频信息，频带较窄，所以该方法求取的波阻抗纵向分辨率低、精度小。法使用测井等资料约束，因此反演结果是相对值、纵向分辨率不高，以及没有井数据的研究区。由于反演结果不能反映地层的真实波阻抗，定量计算地层的岩性、物性参数。2.1.2递推反演首先利用实际地递推反演是在地下波阻抗界面的反射系数为稀疏分布的假设条件下，首先利用实际地震记录估算波阻抗界面的反射系数序列，然后用递推方法求取地层波阻抗的反演方法。地下介质波阻抗的递推公式推导如下：公式(2-2)中当第0层反射系数&和波阻抗AI(0)已知时，可以导出第1层的波阻抗AI(1)为1+R1+RAI(1)=AI(0) 01-Ro(2-6)当第0层反射系数Ro、第1层波阻抗Al(1)已知时，可推导出第 0层AI(0)为1—RoAI(0)=AI1—RoAI(0)=AI⑴1+R0(2-7)当已知第m(0<mvn)层波阻抗AI(m)和m到n层间波阻抗界面的反射系数序列R(i=0,1,2，…,n-1R(i=0,1,2，…,n-1时，第n层地层波阻抗AI(n)为1+RAI(n)=AI(n—1)• =AI(n-2)•1-Rn_l1+也1—Rz1+也1—也=...=AI(m)•^Rn^1-Rn41+R=...=AI(m)•^Rn^1-Rn41+Rm41一RmJ1+Rm1—Rm=AI(m心(0兰men)(2-8)当已知第n层波阻抗Al(n)和m到n间波阻抗界面的反射系数R(i=0,1,2，…，R(i=0,1,2，…，n-1)时，第m(0<mcn)层地层波阻抗AI(m)为AI(m)=AI(m+1)•^^Rm=Al(mAI(m)=AI(m+1)•^^Rm=Al(m+2)•1+Rm1+Rm 1+Rm半=...=AI(n)•^l^m1+Rm1-Km*1+1-尺_21+Riz1~Rn^1+备21-R=AI(叫齢(0<mcn)(2-9)通过公式（2-8）和（2-9）可知，当已知第k层地层波阻抗Al（k）和该层与待求层第 n层的间各层界面的反射系数序列时， 可以递推出k和n层间各层的波阻抗值。实际应用中一般选取大套泥岩层作为标准层求取其波阻抗，然后根据反射系数序列递推出所有地层的波阻抗值。递推反演中最重要的部分是反射系数序列的求取， 反射系数的稀疏程度对反演结果影响很大，求取合适的反射系数序列关乎反演的成败。稀疏脉冲反褶积是实际生产中常用的求取递推反演反射系数的方法。 图2-1是稀疏脉冲反褶积求取反射系数序列流程图。首先采用最大似然反褶积等方法估算波阻抗界面的反射系数序列，并经褶积求取合成地震记录与实际地震资料的差异反馈给反射系数序列求取过程，根据反馈结果适当修正反射系数序列， 当合成地震记录与实际地震资料的残差满足预设条件时，反射系数序列既是波阻抗界面的稀疏反射系数，然后利用公式（ 2-8）和（2-9）就可以计算出波阻抗。递推反演结果能够反映地下波阻抗的分布规律， 而且递推算法使得反演结果不存在多解现象；递推反演采用简单的递推计算方法，计算速度快、 对计算机要求较低；递推反演纵向分辨率高于道积分反演，反演结果是绝对波阻抗， 能够用于储层预测与和砂体雕刻。但是递推反演采用的递推计算可以导致波阻抗误差的层层累计， 反演结果受地震资料固有频带宽度的影响，对薄储层的识别能力有限，不能满足薄储层的识别和分析， 因此主要用于勘探开发前期。

2.2模型反演基于模型的反演是充分利用地震、地质、测井等资料建立研究区的宽频带波阻抗模型，根据地震波传播原理在该模型上进行正演模拟求取合成地震记录， 对比其与实际地震资料差异并将两者残差反馈到模型上， 根据残差不断修正波阻抗模型再次正演， 直到正演模拟结果与实际地震资料匹配较好、残差满足预定要求时，正演模型就是我们要求的反演结果。 基于模型的反演是一个不断正演、不断修改模型再次正演的过程，图 2-2是模型反演的流程图。I地震资料II地震资料I［测井资料II地质资料II地震子波I地震子波］I波阻抗模型J正演模拟I合成地震记录I合成地震记录I修改模型是否满条件/?（反演结果图2-2模型反演流程图根据模型反演的思路，适当的初始波阻抗模型和子波以及合理的残差是模型反演的关键，它们决定着反演结果的模型化程度和模型修改次数。 因此，基于模型的反演中要注意以下几个方面：（1） 波阻抗模型的建立，模型反演是在初始波阻抗模型的基础上进行正演模拟计算，模型对反演结果影响很大， 因此建立一个好的波阻抗模型是进行模型反演的基础。 在建立波阻抗模型之前，一般需要对测井数据进行标准化处理，对地质信息 （构造发育史、地层断层接触关系、地层岩性信息）有一个充分的理解，搭建合适的地层框架， 在建模时需要选择合适的内插函数。（2） 地震子波的求取，进行高精度的地震记录标定是求取子波的关键，只有选取合适的地震子波（包括频率和相位），正演结果和实际地震资料才能匹配。（3） 残差的大小，合成记录与地震资料的残差决定了模型修改次数，只有设定一个合理的残差范围才能使反演结果真实反映地下波阻抗情况，又不会进行过多的迭代次数。基于模型的反演相比于地震道直接反演方法有很多优点： 首先，它不需要假设波阻抗反射截面的反射系数为深度的连续可微函数或者稀疏分布， 使得反演结果上更能反映地下真实信息；其次，地震子波可以通过精细层位标定求取，使得子波更接近真实情况；该反演方法充分利用了地震、地质和测井等资料建立地质模型， 一定程度上约束了反演结果， 使反演误差不随深度累计；最后，模型反演使用测井资料建模，一定程度上拓宽了反演结果频带宽度，提高了纵向分辨率。模型反演方法也有其致命的缺点： 由于反演是在模型基础上进行的正演计算，反演结果的模型化严重；不同的地下地质情况可能产生相似的地震记录， 反演过程对地质模型有很高的依赖性，使得反演结果具有多解性。不管是地震道直接反演还是模型反演， 反演过程都依赖于实际地震资料， 反演结果受地不管是地震道直接反演还是模型反演， 反演过程都依赖于实际地震资料， 反演结果受地2.3约束稀疏脉冲反演约束稀疏脉冲反演是波阻抗模型约束下的地震道递推反演方法，它与稀疏脉冲反演的

区别就是初始波阻抗模型和测井资料的约束， 约束稀疏脉冲反演既能像递推反演一样快速的求解出反演结果，又不会使递推的误差不随深度累积。约束稀疏脉冲反演的实现方法是将测井资料中声波和密度曲线生成的波阻抗通过地质框架内插得到一个初始波阻抗模型，然后用波阻抗模型以及测井曲线的趋势来约束反演结果，求取地下地层反射系数序列。 约束稀疏脉冲反演得到的波阻抗数据体既与地震资料匹配又与测井数据吻合，更能反映地下地质体的真实分布。在满足合成地震记录与实际地震资约束稀疏脉冲反演的关键是求取稀疏反射系数序列，

料精度的前提下，最少脉冲数目反射系数序列的目标函数为在满足合成地震记录与实际地震资J=2rJ=2riP+汀送(di-Si)q2(ti-z)(2-10)式中J式中J是目标函数，ri是估算的反射系数序列,di为原始地震资料，Si为合成地震记录，A是合成地震记录与实际地震资料残差的权重系数,t是合成地震记录与实际地震资料残差的权重系数,ti是测井波阻抗变化趋势， z是反演选取a为测井波阻抗趋势对反演结Ja为测井波阻抗趋势对反演结J由三部分组成，第一项2ri在勘探期可指导地质统计学反演的以井上波阻抗变化趋势为中心的硬约束范围内的反演值,果的约束程度。一般情况下取ot=1、p=1、q=2。通过式（2-10）可知，约束稀疏脉冲反演最优化目标函数是反射系数序列r的P次方之和，第二项汕无（dj-sjq是反演合成地震记录s与实际地质资料di残差q次方之和，第三项a2送（ti-乙）2是反演波阻抗与测井趋势均方差之和。式中前两项是互相制约、互相影响的，两者对反演结果的影响依靠 A值大小调控：A值太大时，过分强调合成地震记录与实际地震资料的残差小， 反演结果会引入地震噪声、忽略了反射系数序列的稀疏性；兀值太小时，又过分强调反射系数序列的稀疏性、 反射系数个数太少，使得反演结果缺失细节、降低反演结果分辨率，此时合成地震记录与地震资料的吻合性太低，残差太大。因此为了使反演结果分辨率较高，又不引入太多噪音，需要正确选择几值的大小。测井曲线以及初始波阻抗模型的约束使得约束稀疏脉冲反演结果的频带宽度较地震资料有所提高，但是地震资料在反演中还是起着主要作用， 反演结果的频带没有太大提高， 而且容易产生窗帘效应。为了提高反演结果的可靠性及分辨率， 需要将地质、钻井测井和地震这样才能获得相对宽资料建立的波阻抗模型的低频、高频成分通过滤波补偿到反演结果中，频带的反演波阻抗体。这样才能获得相对宽约束稀疏脉冲反演结合了直接反演和模型反演的优点，有很广泛的应用，用于储层预测，指导井位部署；在开发期可预测储层的横向变化规律，参数的选择、反演结果的优选。2.4地质统计学反演

-几十赫兹，因此反演结果不能直接应用到精细储层识别的工作-几十赫兹，因此反演结果不能直接应用到精细储层识别的工作将研究区井上资料分析的地质统计分析合成地震记录与实际地震资料模拟结果就是反因此反将研究区井上资料分析的地质统计分析合成地震记录与实际地震资料模拟结果就是反因此反使得反演结果具有不根据前期地质认识演结果。地质统计学反演的优点是利用了地震数据连续性和井上资料纵向高分辨率,演结果能用于薄储层的储层预测。使得反演结果具有不根据前期地质认识地质统计学反演由于地质统计学的统计规律和随机模拟的随机性，确定性和多解性，因此一般地质统计学反演会求出多个等概率的数据体，和确定性反演结果分析，优选最符合实际情况的结果。2.4.1地质统计学地质统计学(Geostatistics)是法国G.Matheron教授通过大量理论实验研究提出的一种统计学方法。地质统计学通过分析样本数据的变差函数和概率密度函数， 研究区域化变量的空间分布规律，地质统计学考虑了传统统计学中样本值，同时又考虑了样本间的空间关系。241.1区域化变量以往我们认为空间中点 x是坐标(xu,Xv,Xw)的函数(xu,Xv,Xw)=z(x)以往我们认为空间中点 x是坐标(xu,Xv,Xw)的函数(xu,Xv,Xw)=z(x)，是一个确定的值，只是空间坐标(Xu，Xv,Xw)的函数,但是在现在我们认为区域化变量是随机变化的,只有进行一次观测我们才能得到它的一个实现z(x)，在进行观测之前我们不能通过空间函数来确定。因此，区域化变量既是随机函数又是空间场的函数。 区域化变量的随机性在资源勘查中能够有效的反映地质体的空间局限性、连续性和各向异性。区域化变量的特征可以用方差函数、协方差函数以及变差函数等进行描述。2.4.1.2协方差函数在地质统计学中，可以用协方差来描述区域化变量的差异。随机向量为X与X的协方差Cov(X,Y)=E[(X—EX)(Y—EY)](2-11)区域化变量z(x)NXu’Xv'Xw)的(自)协方差可是变量 z(x)在空间X和X+h两点处z(x)和z(x+h)的二阶混合中心矩，即Cov(z(x),z(x+h))=E[z(x)z(x+h)]-E[z(x)]E[z(x+h)](2-12)由公式可知，区域化变量z(x)的协方差是空间点X和偏移量h的函数。假设z(xi)为z(x)在空间位置xi处的一个观测值(实现)，z(Xi+h)是z(x)在xi+h的观测值(实现)[i=12…，N(h)]，由式(2-12)可得:(2-13)c^h)=丘F[z(x)—z(Xi)][z(x+h)_z(Xi+h)]N(h)y(2-13)式中Z(Xi)、Z(Xi+h)］为空间x和x+h两点的平均期望，假设 z(Xi^z(x+h^m式中Z(Xi)、Z(Xi+h)］为空间为常数)，c(h)可简化为c(h)=N^F)[z(x)z(Xc(h)=N^F)[z(x)z(Xi+h)-m2]N(h)y(2-14)2.4.1.3变差函数地质统计学分析中，以直方图和变差函数的形式来反映地质变量空间分布的随机性和相关性。它是指地质变变差函数是一种能定量表示空间相关关系或地质变量连续性的数学方法。它是指地质变量z(x)在X和x+h两点处的增量的方差之半， 即区域变化量在相距为 h的任意两点处的平方差值的一半。其数学表达式为1Yg^EIzg—z(x+h)]1Yg^EIzg—z(x+h)]2(2-15)对于z(x)的观测值(实现)，z(x)是离散的，表达式(2-15)可改写为N(h)2N(h)N(h)2N(h)泸Xi)—z(Xi+h)]2(2-16)CCo、基台值c和变程a。图在地震储层反演中，变差函数