MP算法在地震波阻抗反演中的应用.pdfVIP

第4l卷 第 1期 成都理工大学学报 (自然科学版) Vo1．41 No．1 2014年 2月 JOURNALOFCHENGDUUNIVERSITYOFTECHNOLOGY(ScienceTechnologyEdition) Feb．2014 [文章编号]1671—9727(2014)O1—0087—07 MP算法在地震波阻抗反演中的应用 周东勇 ，文晓涛 ，贺振华 ，蒲 勇 ，黄德济 ，胡军辉 (1．油气藏地质及开发工程国家重点实验室(成都理工大学)，成都 610059； 2．中国石化股份有限公司 勘探南方分公司，成都 610041) [摘要]当孔隙中充填油气后 ，波阻抗相对较低 ，准确反演地下地质体波阻抗信息可在一定程 度上为油气检测提供依据。利用基于雷克子波原子库 的匹配追踪 (MP)算法对地震信号进行 稀疏分解 ，进而可实现对地下地质体 的高精度波阻抗反演 。模型试算表 明，MP算法可 以用于 地震信号的稀疏分解，且两层 间厚度 2／4调谐时，仍可大体反映层位及系数的相对大小关系。 通过对实际地震资料的计算表 明，该算法也适用于实际资料 ，且与实钻情况吻合。 [关键词]匹配追踪算法；波阻抗反演；稀疏分解；雷克子波 [分类号]P631．443 [文献标志码]A 波阻抗反演开始出现于 20世纪 70年代初 ， 函数集 (函数集中的任意一个函数称为基 函数)上 是一种 以研究岩性油气藏为主 的地震勘探技术 。 进行分解 ，然后在变换域上表达原始信号或原始 目前 ，波阻抗反演 的方法有很多 ，如模拟退火法 、 信号的主要成分 。若在变换域上用尽量少 的基 函 遗传算法等等。本文研究采用匹配追踪算法对地 数来表示原始信号或原始信号 的主要成分，就是 震数据进行分解 ，得相对反射系数 ，进而依据波阻 信号的稀疏表示，得到信号稀疏表示的过程就是 抗公式进行波阻抗反演。匹配追踪算法 (Matc— 信号的稀疏分解_6]。若用这种分解过程分解地震 hingPursuits，本文统一简写为MP算法)由Mal— 信号，则是地震信号的稀疏分解。MP算法就是 lat和 Zhang于 1993年首先提出口]，之后很快被 实现地震信号稀疏分解的算法之一。 应用于地震信号处理领域。2003年 Castagna等 设研究的地震信号为 -厂，信号长度为 N，且 厂 将 MP算法用于地震谱分解_2；2007年 Y．H_ ∈H，H为希尔伯特空间。MP算法将原始地震 Wang将该算法应用于时频分析口；2012年张繁 信号 ．厂稀疏表示为 昌等将该算法用于瞬时谱识别来确定三角洲砂体 n 1 尖灭线 等等 。虽然应用 MP算法对信号稀疏分 f一 cig ， ( 《N ) (1) 解时，计算量十分大，但该算法有很多优 良的性 0 式 中：g 为函数集 中的基函数 (在 MP算法 中函 质 ，如对信号 自适应 的灵活表达 ，较高 的时间、频 数集被称为过完备原子库 ，基 函数被称为原子)； 率分辨率，暂态结构的局部 自适应性，信号结构的 参数表示更加灵活等_l5]。