地震反演原理及其应用
地球物理反演方法的综述
地球物理反演方法的综述
地球物理反演是一种利用地球物理方法来推断地下构造和物质分布的技术。通过观测和测量地球物理场，如重力、地磁、电磁、地震等，结合数理统计和计算机模拟方法，可以对地下的地质构造、岩石性质和地下水资源等进行精确的推断。本文将综述地球物理反演方法的原理、分类及应用。
一、地球物理反演方法的原理
地球物理反演方法的原理在于根据地球物理场的观测数据，通过数学模型和计算方法，将地球物理场与地下介质属性之间的关系联系起来。根据电磁波传播、物质密度、电阻率、磁化率等反演参数的变化规律，推断地下介质的结构和成分。
其中常用的地球物理反演方法包括重力法、磁法、电磁法、地电法和地震法等。不同的反演方法适用于不同的地质介质和研究目标，各有其优势和限制。
二、地球物理反演方法的分类
1. 重力反演法：利用重力场观测数据，通过计算物质的密度分布，来推断地下构造的方法。重力反演法在石油勘探、地质灾害分析、水资源评价等领域具有广泛应用。
2. 磁法反演法：通过磁场观测数据，推断地下磁化率和磁性物质的空间分布。磁法反演在矿产勘探、地震预测等方面发挥重要作用。
3. 电磁法反演法：通过电磁场观测数据，推断地下电阻率分布，来研究地下水资源、矿产和工程勘探。电磁法反演在地下水资源评价、油气勘探、环境地球物理和岩土工程等方面有广泛应用。
4. 地电法反演法：通过电场和电位观测数据，推断地下电阻率分布，用于研究地下水位、地下水性质、污染监测和地下工程等。地电法反演在工程地球物理勘探和水文地球物理领域具有广泛应用。 5. 地震法反演法：通过地震波在地下的传播与变化，推断地下介质的速度和密度分布，用于研究地质构造、地震预测和石油勘探等。地震法反演是地球物理反演方法中应用最广泛的方法之一。
三、地球物理反演方法的应用
地球物理反演方法广泛应用于地质探测、资源勘探、环境监测和工程勘察等领域。以下是几个常见的应用领域：
1. 石油勘探：地震反演方法可用于确定油气藏的位置、大小和分布，辅助油田开发和管理。
反演常用方法
地球物理反演就是利用观测数据恢复地下地
质结构和岩石性质的方法 ,狭义地说 ,反演就是从
有限频带宽度的地震数据中恢复出宽带波阻抗。
地震反演从方法上大致可分为基于波动理论
的波动方程反演和基于褶积模型的反演两大类。
以波动方程为基础的地震反演由于算法结构复杂 ,
对地震噪音敏感等因素 ,要得到一个稳定的解是比
较困难的 ,因此 ,这类方法还未得到广泛的实际应
用。以褶积模型为基础的地震反演 ,由于算法简
单 ,对地震噪音敏感性小 ,一般情况下都能得到一
个稳定的解 ,在生产中得到了广泛地应用。
稀疏脉冲法
包括最大似然反褶积、L1范数反褶积、最小熵反褶积、最大熵反褶积、同态反褶积等，稀疏脉冲反演是基于脉冲反褶积基础上的递推反演方法,其基本假设是地层的强反射系数是稀疏分布的。从地震道中根据稀疏的原则提取反射系数,与子波褶积后生成合成地震记录;利用合成地震记录与原始地震道残差的大小修改参与褶积的反射系数个数,再作合成地震记录;如此迭代,最终得到一个能最佳逼近原始地震道的反射系数序列。该方法适用于井数较少的地区,其主要优点是能够获得宽频带的反射系数,较好地解决地震反演的多解性问题,从而使反演结果更趋于真实。
约束稀疏脉冲反演采用一个快速的趋势约束脉冲反演算法,用解释层位和井约束控制波阻抗的趋势和幅值范围,脉冲算法产生了宽带结果,恢复了缺失的低频和高频成分;同时,再加入根据井的波阻抗的趋势约束。约束稀疏脉冲反演最小误差函数是J=∑(ri)p+λq∑(di-si)q++α2∑(ti-Zi)2(1)
式中:ri为样点的反射系数;zi为样点的波阻抗;di是原始地震道;si是合成地震道;Zi介于井约束的最大和最小波阻抗之间;ti是用户提供的波阻抗趋势;α为趋势最小匹配加权因子;p,q为L模因子;i是地震道样点序号;λ为数据不匹配加权因子。
如果从最大似然反褶积中求反射系数r(t),则在上述过程中为了得到可靠的反射系数估计值,可以单独输入波阻抗信息作为约束条件,从而求得最合理的波阻抗模型
全波形反演的工作原理
全波形反演的工作原理
全波形反演的工作原理
简介
全波形反演是一种用于地球物理勘探的前处理技术，旨在通过地震波的传播和反射信息来重建地下介质的物理属性。它是一项复杂而强大的技术，被广泛应用于石油勘探、地震监测以及其他地球物理领域。
传统反演方法的局限性
在传统的地震勘探中，通常只使用地震波的到达时间和振幅信息来判断地下结构。然而，这种方法无法提供足够详细的地下信息，并且很难区分不同类型的地质体。
全波形反演的原理
全波形反演通过利用地震波在地下介质中的传播和反射规律，以及地震波在地下界面上的反射和透射现象，来反推地下介质的物理属性。其工作原理可概括为以下几个步骤：
1. 正演模拟：根据已知的地下介质模型和波源信息，使用数值模拟方法计算地震波的传播过程和在地表或接收器上的观测结果。 2. 误差函数定义：将地震波的观测数据与模拟结果进行比对，计算它们之间的差异，常用的误差函数包括最小二乘法、互相关等。
3. 反演算法：通过不断调整地下介质的模型参数，以使误差函数最小化，来逐步重建地下模型。常用的反演算法有梯度下降法、共轭梯度法等。
4. 迭代计算：反演过程是一个迭代的过程，需要不断更新地下模型并重新进行正演模拟和误差函数计算，直到满足反演收敛的条件。
虚拟实验与验证
全波形反演是一项理论和实践相结合的技术，虚拟实验和验证是其关键步骤之一。虚拟实验通过利用已知的地下介质模型，产生合成的地震波观测数据，并将其输入到反演算法中进行反演。通过与实际观测数据进行比对，可以验证反演结果的准确性和可靠性。
应用领域
全波形反演技术在地球物理勘探领域具有广泛的应用，包括但不限于以下几个方面：
• 油气勘探：全波形反演可帮助解决油气勘探中的各种问题，如油藏形态预测、岩性识别、储层描述等。
• 地震监测：全波形反演可用于地震活动的监测和分析，包括地震源研究、震源机制解析等。 • 地下水勘查：全波形反演可用于地下水资源勘查，帮助确定地下水的分布、流动性质等。
波阻抗反演
波阻抗反演通常指利用叠后地震资料进行反演的一种技术，它将地震资料、测 井数据、地质解释相结合，充分利用测井资料具有较高的垂向分辨率和地震剖 面有较好的横向连续性的特点，将地震剖面“转换”成波阻抗剖面，不仅便于 解释人员将地震资料与测井资料连接对比，而且能有效地对地层物性参数的变 化进行研究，从而得到物性参数在空间的分布规律，指导油气的勘探开发， 地震反演的方法主要有两种，一种是叠前反演，一种是叠后反演，叠前反演 主要有旅行时反演和振幅反演，叠后反演主要分为振幅反演和波场反演 。而我 们这里所说到的波阻抗反演属叠后振幅反演， 主要有 递归反演、稀疏脉冲反演和 基于模型的反演这三种方法 。
二、 波阻抗反演方法介绍
1、波阻抗反演的基本假设前提
1、波阻抗反演的基本假设前提 目前我们常用的波阻抗反演软件所用方法基本都是基于褶积模型的基础 上建立的， 因此要求资料都要满足褶积模型的假设前提， 基本可概括为下面的四 个方面：
(1)、地震模型 假设地层是水平层状介质，地震波为平面波法向入射，其地震剖面为正入 射剖面，并且假设地震道为地震子波与地层反射系数的褶积。
(2) 、反射系数序列 在普通递归反演中，假设反射系数为完全随机的序列，而在稀疏脉冲反演 中，假设反射系数为由一系列大的反射系数叠加在高斯分布的小反射系数的背景 上构成的。
(3)、地震子波 假设反射系数剖面中的每一道都可以看作是地下反射率与一个零相位子波 的褶积。实际情况下往往需要对地震剖面进行相位校正处理
(4)、噪音分量 通常假设波阻抗反演输入的地震数据其振幅信息反映了地下波阻抗变化情 况，地震剖面没有多次波和绕射波的噪音分量。 因此，在资料处理时可以考虑的 处理流程是反褶积、 噪音剔除，尤其是多次波， 处理的最终目标是得到真振幅剖 面。类似二维滤波和多道混波这样可以改变地震振幅和特征的处理模块应当避免
使用。
有许多反演技术都存在两个问题：一是多解性，即存在多个反演结果与地 震数据相吻合；另一个问题是地震数据的带限问题。
地震波阻抗反演方法综述
地震波阻抗反演方法综述
一、地震反演技术研究现状
地震反演方法是一门综合运用数学、物理、计算机科学等学科发展起来的新技术新方法，每当数学方法、物理理论有了新的认识和发展时，就会有新的地震反演技术、方法的提出。随着计算机技术的不断发展、硬件设施的不断升级，这些方法技术得到了实践验证和提升，反过来地震反演技术运用中出现的新问题、新思路又不断促使数学方法、地球物理学理论的再次发展。时至今日，地震反演技术仍然是一个不断发展、不断成熟、不断丰富着的领域。
反演是正演的逆过程，在地震勘探中正演是已知地下的地质构造情况、岩性物性分布情况，根据地震波传播规律和适当的数学计算方法模拟地震波在地下传播以及接收地震波传输到地表信息的过程。地球物理反演就是使用已知的地震波传播规律和计算方法，将地表接收到的地震数据通过逆向运算，预测地下构造情况、岩性物性分布情况的过程。地震波阻抗正演是对反演的理论基础和实现手段。
1959年美国人Edwin Laurentine Drake在宾夕法尼亚州开凿的第一口钻井揭开了世界石油工业的序幕。从刚开始的查看地质露头、寻找构造高点寻找石油，到通过地震剖面的亮点技术寻找石油，再到现在运用多种科学技术手段进行油气资源的预测，石油勘探经历了一个飞速的发展历程。
声波阻抗（AI）是介质密度和波在介质中传播速度的乘积，它能够反映地下地质的岩性信息。声波阻抗反演技术是20世纪70年代加拿大Roy Lindseth博士提出的，通过反演能够将反映地层界面信息的地震数据变为反映岩性变化的波阻抗（或速度）信息。由于波阻抗与地下岩石的密度、速度等信息紧密联系，又可以直接与已知地质、钻井测井信息对比，因此广泛应用于储层的预测和油藏描述中，深受石油工作者的喜爱。70年代后期，从地震道提取声波资料的合成声波技术得到了快速发展，以此为基础发展的基于模型的一维有井波阻抗反演技术，提高了反演结果的可靠性。进入80年代，Cooke等人将数学中的广义线性方法运用于地震资料反演，提出了广义线性地震反演。此后Seymour等人又提出了测井声波资料和地震数据正反演相结合求取地下声波阻抗的测井约束反演，大大拓宽了反演结果的纵向分辨能力。
地震反演
1 第一章 反演理论
第一节 基本概念
一． 反演和正演
1． 反演
反演是一个很广的概念，根据地震波场、地球自由振荡、交变电磁场、重力场以及热学等地球物理观测数据去推测地球内部的结构形态及物质成分，来定量计算各种有关的物理参数，这些都可以归结为反演问题。在地震勘探中，反演的一个重要应用就是由地震记录得到波阻抗。
有反演，还有正演。要正确理解反演问题，还要知道正演的概念。
2． 正演
正演和反演相反，它是对一个假设的地质模型，给定某些参数（如速度、层数、厚度）用理论关系式（数学模型）推导出某种可测量的量（如地震波）。在地震勘探中，正演的一个重要应用就是制作合成地震记录。
3． 例子
考虑地球内部的温度分布，假定地球内部的温度随深度线性增加，其关系式可表示成：
T(z)=a+bz
正演：给定a和b，求不同深度z的对应温度T(z)
反演：已经在不同点z测得T(z)，求a和b。
二．反演问题描述和公式表达的几个重要问题
1．应用哪种参数化方式——离散的还是连续的？
2．地球物理数据的性质是什么？观测中的误差是什么？
3．问题能不能作为数学问题提出，如果能够，它是不是适定的？
4．对问题有无物理约束？
5．能获得什么类型的解，达到什么精度？要求得到近似解、解的范围、还是精确解？
6．问题是线性的还是非线性的？
7．问题是欠定的、超定的、还是适定的？
8．什么是问题的最好解法？
9．解的置信界限是什么？能否用其它方法来评价？
第二节 反演的数学基础
2 一．解超定线性反问题
1．简单线性回归
可利用最小平方法确定参数a、b使误差的平方和最小。
22)()(xxnyxxynbxbynxbya （1-2-1）
拟合公式为：
bxayˆ （1-2-2）
该方法的公式原来只适用于解超定问题，但同样适用于欠定问题，当我们有多个参数时，称为多元回归，在地球物理领域广泛采用这种方法。此过程用矩阵形式表示，则称为广义最小平方法矩阵方演。
地震反演 2
地震反演
姓名：李雪松 班级：油气田s101 学号：201071059
一、地震反演的基本定义
通俗的讲就是由地震为基础加上其他条件为约束推测出地层岩性构造的过程叫地震反演。
把常规的界面型反射剖面转换成岩层型的测井剖面，将地震资料变成可与测井资料直接对比的形式，实现这种转换的处理过程叫地震反演。
地震反演：地震反演是利用地表观测地震资料，以已知地质规律和钻井、测井资料为约束，对地下岩层空间结构和物理性质进行成像(求解)的过程，广义的地震反演包含了地震处理解释的整个内容。
地震多解性和粗略性
地震反演多解性是指同一地震资料可对应用不同的岩层结构，粗略性是指推断的参数少，分辨率低，前者可能导致地下模型的错误，后者影响模型的精度。
理论基础：
鲁宾逊褶积模型基础。
其实现过程是：
（1）子波提取；
（2）逐道修改波阻抗模型道，与子波褶积合成地震记录，使之与相对应地震记录相似度最大化（相关系数最大或绝对误差最小），逐道外推，直到完成全剖面的波阻抗反演。
叠偏地震记录X(t)可表示为：
X(t)=K(t)*W(t) (1)
式中W(t)为地震子波，K(t)为反射系数。
从井出发优化波阻抗模型，并正演合成地震记录，使之与相对应地震道相关度最大化，形成反演成果道，选择反演成果道中相关系数达到标准的阻抗模型，以此为基础点建立下一道的初始波阻抗模型，并进行上述优化，直至完成全剖面的反演。
技术特点：
1.采用逐道相关外推建立（优选）初始阻抗模型。 2.反演成果纵向分辨薄层的能力较强。
前提条件：要有地震偏移资料，构造沉积解释层位，标准化后的声波和密度测井曲线，如有其它相关资料更好。
优点：逐道外推波阻反演对井的依赖较小，单井时通常也有较高的精度，整体建模反演，适应于岩性剧烈变化的地带，井多时反演精度较高。
反演原理——精选推荐
3 第二章 地震反演技术 地震反射波法的基础是由于地下不同的地层存在着波阻抗的差异，从而形成了反射波法。所以从本质上来讲，地震反演的目标就是根据已经获得的地震反射波形，以已知地质规律和钻井、测井资料为约束，对地下岩层空间结构和物理性质所进行的成像（求解），广义的地震反演包含了地震资料处理解释的整个内容。 波阻抗反演是利用地震资料反演地层波阻抗（或）速度的地震特殊处理解释技术。与地震模式识别预测油气、神经网络预测地层参数、振幅拟合预测储层厚度等统计性方法相比，波阻抗反演具有明确的物理意义，是储层岩性预测、油藏特征描述的确定性方法，在实际应用中取得了明显的地质效果，因此地震反演通常特指波阻抗反演。李庆忠院士指出：“波阻抗反演是高分辨率地震资料处理的最终表达方式”，说明了波阻抗反演在地震技术中的特殊地位。 地震数据的反演可以用图2.1所示的框图来概括（R.Brain，1993）。最完善的反演方法是叠前反演，它分为振幅反演（如AVO分析，叠前波动方程波形反演等）和旅行时反演（常称层析法）。然而，由于叠前反演信噪比低，稳定性差，分辨率低，正演模拟困难以及计算量大等原因，所以，叠后反演仍然是当前最常用的方法。 图2.1 地震反演方法概况示意图 叠后地震道反演方法可分为：①递推反演，如D.W.Oldenoburg（1983）和C.Walker（1983）的最大熵（MED）及自回归（AR）方法，B.Ursin（1985）的最大似然反褶积（MLD）方法，vielle（1991）的贝叶斯估计反褶积（MED）；②广义线性反演（GLI），如D.A.Cooke（1983）；③非线性约束反演（BCI），如B.Cornish （1988）的宽带约束反演（BCI）、S.Gluck（1989）的随机反演（ROVIN）、振 幅 法 旅 行 时 法 稀疏脉冲法 模 型 法 带 限 法 地 震 反 演 叠 前 反 演 叠 后 反 演 宽 带 法 4 R.D.Martinez（1988）的多参数约束反演（包括BCI、WLI和LCI）。根据反演结果的频带特性又可分为：①带限法（如GLI、合成声波测井）；②稀疏脉冲法（如MED、AR、MLD、BED方法等）；③模型法（如BCI、ROVIN、LCI、Strata和Jason方法等）。但不管怎样，拓宽频带，提高分辨率和稳定性以及提高反演结果的精度，永远是叠后反演方法所追求的目标。基于这一原则，R.Brain（1993）根据数值模拟和实际资料的计算结果说明了模型法是最好的。根据测井资料在其中所起作用的大小又可分成：①地震直接反演；②测井控制下的地震反演；③测井—地震联合反演；④地震控制下的测井内插外推；分别用于油气勘探的不同阶段。从实现方法上可分为：①递推反演；②基于模型反演；③地震属性反演。下面就几种反演方法进行描述。 第一节 递推反演 递推反演是早期的波阻抗反演方法，如Vlog、Seislog等反演软件。假定地震记录经噪音消除，振幅恢复和反褶积以后，可以得到一条反射率或反射系数剖面。这个反射系数剖面实际就是反射系数序列与一理想零相位子波的褶积，并称这一子波为剩余子波。也就是说我们得到的不能是脉冲化的反射率，而是脉冲反射率经剩余子波带通滤波的结果。 将反射率函数转换为波阻抗，其目的在于极薄地层由于子波的滤波作用而不可恢复，常规反射系数剖面上显示不出来，比极限稍厚一点的地层也因常规剖面显示动态范围不足而看不出来。在有利条件下，波阻抗反演后，这部分地层可以鉴别出来。 递推反演的关键在于从地震记录估算地层反射系数，得到能与已知钻井最佳吻合的波阻抗信息。递推反演方法中测井资料主要起标定和控制的作用，因而递推反演又称之为直接反演后测井控制下的地震反演。 2.1.1 方法原理 无噪偏移地震记录的理论模型为： S(t) = r(t)W(t) (2.1.1) 其中：S(t)为地震记录，r(t)为地层反射系数，W(t)为地震子波。通过子波反褶积处理，可由地震记录求得反射系数，进而递推计算出地层波阻抗或层速度。 5 波阻抗是重要的岩石物理参数，可直接与钻井对比进行储层岩性解释和物性分析。 2.1.2 方法实现 递推反演是对地震资料的转化处理过程，其结果的分辨率、信噪比以及可靠程度完全依赖于地震资料本身的品质，因此用于反演的地震资料应具有较宽的频带、较低的噪声、相对振幅保持和准确成像。测井资料，尤其是声波测井和密度测井资料，是地震横向预测的对比标准和解释依据，在反演处理之前应进行仔细的编辑和校正，使其能够正确反映岩石的物理特征。 递推反演的技术核心在于地震资料正确估算地层反射系数（或）消除地震子波的影响），比较典型的实现方法有：基于地层反褶积方法、稀疏脉冲反演和测井控制地震反演等。地层反褶积方法是根据已有测井资料（声波和密度）与井旁道地震记录，利用最小平方估算数学意义上的“最佳”子波或反射系数。这种方法的优点是把子波求解的“欠定”问题变成了确定问题，在井点已有测井井段范围内可获得与测井最吻合的反演结果。局限性主要有：①该方法完全忽略了测井误差和地震噪音，这些因素尤其是前者的客观存在使“子波”确定更加困难；②地层反褶积因子的估算是在计算时窗内数学意义上的最佳逼近，实际处理范围与该时窗的不同超过了该方法的适用范围，即便是在井点位置，得到的反演结果已不可能是“误差最小”。不难看出，影响基于地层反褶积递推反演效果的主要因素是测井资料的质量和地震资料的信噪比以及地震噪声的一致性。 稀疏脉冲反演（Sparse—spike Inversion）是基于稀疏脉冲基础上的递推反演方法，主要包括最大似然反褶积（MLD）、L1模反褶积和最小熵反褶积（MED）。这类方法针对地震记录的欠定问题，提出了地层反射系数由一系列迭加于高斯背景上的强轴组成的基本假设，在此条件下以不同方法估算地下“强”反射系数和地震子波，这种方法的优点是无需测井资料，直接由地震记录计算反射系数，实现递推反演，其缺点为很难得到与测井曲线相吻合的最终结果。 基于频率域反褶积和相位校正的递推反演方法从方法实现上回避了计算子波或反射系数的欠定问题，以井旁反演结果与实际测井曲线的吻合程度作为参数优选的基本数据，从而保证了反演资料的可信度可解释性，是递推反演的主导技术。 图2.2是其处理流程，其主要技术关键有：恢复地层反射系数振幅谱
地球物理学中的地震波传播与反演
地球物理学中的地震波传播与反演
地震波是地震发生时产生的波动，是研究地震学的基础。地震学家借助地震波的传播与反演，可以了解地下构造的情况，进而研究地震活动与岩石物理性质等问题。本文将从地震波的传播机制、地震波反演理论及方法等方面探讨地球物理学中的地震波传播与反演。
一、地震波的传播机制
地震波的传播引起了地壳中的微小变形和位移，导致地震波在地球上传播。地震波主要分为纵波和横波两种，纵波又叫P波（Primary wave），横波又叫S波（Secondary wave）。P波是一种纵波，具有直线传播、传播速度快、能穿透岩石和液态物质的特点；而S波是一种横波，具有像水波一样的传播方式、传播速度慢、只能穿透固体岩石等性质。
地震波在地壳中传播的速度与介质的密度、压缩模量以及剪切模量等因素密切关联。另外，地震波的传播速度受到地壳中不均匀性的影响，地壳中有不同密度的层次，地震波通过不同密度层次时会出现反射、折射等现象，使得地震波路径发生曲折，从而研究地壳结构时要对这些影响因素进行较为精细的考虑。
二、地震波反演理论与方法
能否将地震波数据反演成有关介质结构的有用信息，是地震勘探、地球物理勘探中常常需要考虑的问题。
地震波反演的基本思想是借助地震波在地下介质中传播的情况来推断地下介质的物理参数。通常情况下，为了研究介质的速度、密度、弹性模量、剪切模量等参数，需要通过处理地震波在地下的传播路径和传播时间，从而反演得到介质的物理结构。
地震波反演的方法有很多种，主要包括正演法、反演法和拟合法。正演法指利用已知参数的介质来计算地震波在介质中的传播规律。反演法是利用地震波在介质中所传递的信息，探索出地下介质的物理参数。拟合法主要是利用地震波在介质中的传播速度随深度分布变化的规律来拟合地下介质的物理结构。
在地震波反演中，数据处理也是非常重要的一环。地震波的反演可以通过复杂的图形工具和数学模型来完成。比如从地震带上提取的地震记录中得到横波和纵波，分别对横波和纵波进行分析、处理，再分别反演有关介质信息。
测井约束地震反演技术分析及其应用
第１３卷第５期 断 块 油 气 田 ＦＡＵＬＴ．ＢＩ ＯＣＫ ０ＩＬ＆ＧＡＳ ＦＩＥＩ Ｄ ２００６年９月
测井约束地震反演技术分析及其应用
王 波
（中国地质大学能源学院，北京１０００８３） 摘 要 首先从测井约束地震反演技术的基本原理出发，然后从储层地球物理特征、地震 子波提取、初始波阻抗模型的建立３个重要技术环节进行分析，最后根据文明寨油田Ｍ１５块 实际资料情况和地质特点，应用测井约束反演技术对其进行了储层反演，反演结果与钻井情况 吻合良好，为下一步对该区的开发钻探指明了方向。 关键词测井约束地震反演地球物理地震子波波阻抗
随着东濮凹陷勘探程度的不断提高，勘探工 作已由原来的寻找构造油气藏逐渐向寻找隐蔽岩 性油气藏过渡，测井约束反演作为一种储集层横 向预测技术，利用地震资料在井的约束下进行目 标储集层反演，目前广泛应用于各油田的勘探开 发过程中。在对文明寨地区Ｍ１５块的实际资料 情况和地质特点进行综合分析的基础上，应用 Ｊａｓｏｎ软件采用约束稀疏脉冲反演方法对其进行 储层预测，取得了很好的效果。 １ 方法原理 稀疏脉冲反演的思路是假定地下地层的反射 系数是稀疏分布的，是在高斯分布的一系列小反 射系数的背景上，叠加了一些大的反射系数。通 过对每一道逐点增加脉冲，反射系数被逐渐改变， 直到找到合适的位置使该反射系数序列合成的地 震道和实际地震道的误差达到最小。这时的反射 系数序列模型就是最终的反射系数序列。当反射 系数求出后，就可以求出相对波阻抗。加上井的 低频约束就可以得到绝对波阻抗…。约束稀疏 脉冲反演的核心是优化如下目标函数： Ｅ＝∑（ ） ＋Ａ ∑（ｄ —ｓｉ） ＋ Ｅ（ｔｉ一 ） 式中： 为反射系数序列； 为声阻抗序列；ｄ 为 地震道序列；ｓ 为合成地震道序列；ｔｉ为用户提供 的趋势序列；Ａ为残差权重因子；ｃｔ为趋势权重因 子；ｐ、ｑ为Ｌ模因子。在质量控制下可以找到一 个合适的Ａ值，在一定的可信度下提高剖面的分 辨率。 ２主要技术环节