压缩感知理论重构地震数据研究.pdfVIP

·672· 中国地球物理2011 (18)储层地球物理 压缩感知理论重构地震数据研究 on onseismicdata based sense Study recovery compressivetheory 曹思远’路交通 Cao Siyuan+LuJiaotong 中国石油大学(北g)CNPC物探重点实验室北京昌平102249 1．研究意义 伴随着国内外油气勘探对象由简单到复杂、由构造向岩性转变，时移地震油气藏监测等技术对地震数 据的规则性以及完整性提出更高要求，然而由于各种因素导致地震数据不规则、不完整。缺失的地震数据 不进行处理对后期多道处理将产生较大影响，如将导致三维面元覆盖的不均匀，偏移成像时产生构造假象 等等。常规处理方法难以满足实际生产需要，将最近提出的压缩感知理论用于地震数据重构，通过模拟的 地震记录以及Marmousi2模型，验证该方法的可行性与有效性。 2．研究方法 压缩感知理论认为当一个信号在某个正交空间具有稀疏特性或者该信号是可压缩时，可通过一个测量 矩阵(与稀疏变换正交基不相关)将稀疏变换得到的高维信号投影到低维空间，然后通过一定重构优化算 法可较好恢复原始信号。它包括三个方面：稀疏变换、测量矩阵、重构算法fI2】。依据压缩感知理论体系， 首先将缺失地震数据变换到Curvelet稀疏域，然后通过测量矩阵M进行测量，同时引入一种约束矩阵R， 使缺失的地震数据满足随机缺失的特性，这样得到的缺失地震数据在稀疏域表现为随机噪声，再通过一种 地震信号自适应阈值迭代算法【31就可以实现缺失地震数据恢复，最后通过逆Curvelet变换得到重构后地震 数据。选取第二代离散Curvelet变换作为稀疏变换，因其良好的方向性、局部性以及各向异性，使得缺失 的地震数据表现更优稀疏特性，便于后面重构处理；引入约束矩阵R使得地震数据不需满足等间距缺失， 克服常规f-x与F．K插值算法的局限性：同时采用地震信号自适应阈值迭代重构算法，是因其能克服常规 软阈值与硬阈值迭代算法的缺点，取得较好重构效果。 3．模型测试 首先通过高阶有限差分方法模拟一个四层水平均匀介质，通过中间放炮两端接收得到236道地震记录， 通过约束矩阵R使得模拟得到的地震数据满足随机缺失的特性，然后通过在Curvelet稀疏域的自适应阈值 迭代重构算法，经过100迭代，缺失的地震数据得到较好重构。同时通过对比重构前后地震数据的F．K谱 也可明显看出重构效果较好。 为进一步验证压缩感知理论重构缺失地震数据的可行性与有效性，选取国际通用的Marmousi2模型用于测 理论体系重构后，模型数据信噪比达到74．96db，可发现信噪比提高10倍左右，地震数据得到良好恢复。 4．结论 根据压缩感知理论的三个主要部分，分别选取Curvelet变换作为稀疏变换，提出一种约束矩阵R以及 采用地震信号自适应阈值迭代重构算法，使得地震数据信噪比有明显提高，较好实现合成地震记录与 Marmosi2模型的重构，证实该方法有效性。同时该重构体系克服常规插值算法对地震数据等间距缺失的要 求：实际处理过程中也可根据需要选取不同的稀疏变换以及重构算法， 以获得不同运算效率或者不同重 构精度的数据，具有较强的灵活性。 参考文献 【l】HerrmannJ，HennenfentG seismicdata withcurvelet Non—parametric recovery frames[J]．Geophysical 2008，173(1)：A1～A5． 【2】Donoho Trans．On D．Compressedsensing[J]．IEEEInformationTheory,2006，52(4)：1289-1306． [3】吴爱弟，赵秀玲．基于曲波变换的地震信号去噪方法【J】．中国石油大学学报，2010，34(3)：30一-3