双谱信号重构技术在重力数据去噪中应用(1).pdfVIP

第48卷第3期 石油物探 V01．48，No．3 2009年5月 GEOPHYSICALPROSPECTINGFORPETROLEUM May，2009 文章编号：1000一1441(2009)03一0307一07 双谱信号重构技术在重力数据去噪中的应用 俞建宝1’2，陈 高3，李正文2 (1．中国石油化工股份有限公司石油勘探开发研究院南京石油物探研究所，江苏南京210014；2．成 都理工大学，四川成都610059；3．中国石油化工股份有限公司勘探南方分公司，四川成都610041) 摘要：一个实信号可以用一个谐波信号序列来拟合，计算实信号的双谱，通过Fourier变换推算出谐波信号的系 数和相位角与实信号双谱振幅谱和相位谱之间的关系，可以实现信号重构。谐波基频的准确估算是信号重构的 关键，经过试验分析，取双谱振幅谱的第一个峰值对应的频率作为估算基频。构成实信号的谐波项数是未知的， 通过设立系数因子控制谐波项数，可以满足不同频率(或波数)段信号的拟合，实现噪声压制。数值模拟和实际 应用结果表明，双谱信号重构技术不仅能有效压制随机噪声，还可以进行异常分离。 关键词：双谱分析；信号重构；重力异常；去噪；异常分离；归一化总梯度 中图分类号：P631．4 文献标识码：A 双谱分析技术形成于20世纪60年代[1]，之后1．1双谱信号重构技术原理 被应用于很多领域[2~6]。在应用地球物理领域，双 双谱信号重构技术是指基于双谱分析技术的 谱分析技术在地震勘探上，被用于子波估计[7’8]、 谐波信号重构‘28|。设z(￡)是均值为。的实函数， 自动层位追踪[9]、小断层检测[10’11]、裂缝检测[12|、 其Fourier变换表示为 砂体预测[13]、油气预测[14]、剩余静校正[15]、动校叠 (1) X(厂)一A(厂)exp[j声(厂)] 加与速度分析[1617]、波阻抗反演[18|、属性参数提 取[11’14’19|、多次波压制[20|、数值模拟[21]等；在重磁 电测量中，被用于大地电磁信号噪声压制[22~24|、高 位。z(￡)的3阶累积量函数表达式为 精度航磁异常油气特征信息提取等[2印；在遥感测 r斗。o C(￡1，￡2)一I 量中，被用于航空遥感图像的质量评定[26]等。 J——。。 基于前人的研究结果，对双谱分析技术的基础 式中：￡。和￡z分别表示2个不同的时间变量。 理论及其计算方法进行了讨论，对双谱分析技术重 构信号的方法进行了改进，并开发编制了方法实现 r。。r。。 软件。通过数学模型检验了双谱分析技术压制噪 B(厂1，厂2)一I C(￡l，￡2)· √—。。J——∞ 声的有效性和可靠性，并将其应用于实际高精度重 力测量数据的处理和解释。 (3) =X(厂，)X(厂2)X。(厂1+^) 1 双谱信号重构技术 ^是与时间变量￡-，￡z相对应的频率变量。进一 功率谱分析实际上是将二阶统计量作为数学 步，可得z(￡)的双谱表达式为 分析工具，在时间域为相关函数，在频率域为功率 谱。相关函数和功率谱都存在一些缺陷，如具有等