李福元 韦成龙 邓桂林 张宝金 张 衡 杨 力

(①自然资源部海底矿产资源重点实验室，广东广州 510760；②自然资源部中国地质调查局广州海洋地质调查局，广东广州 510760)

0 引言

震源信号子波是反射地震勘探中重要的参数，在地震资料处理及后续很多环节中都起着关键作用，包括地震正演模拟、偏移成像、全波形反演、叠前叠后属性反演、多次波预测和多源数据合成处理等[1-8]。更为一般地说，震源信号是地震勘探信号系统的激励函数，也是地震波弹性波动方程中的源函数[9]。

研究震源信号的核心问题是获取震源信号子波。除了采集现场直接观测记录外，目前主要的方法有： ①利用近场记录拟合[10-12]； ②通过专用软件模拟[13-17]； ③从地震数据本身提取[18]。地震数据采集受现场环境、仪器装备等因素影响，从现场直接观测震源子波十分困难。利用近场记录拟合时，因近场记录中常伴有干扰，拟合结果往往存在误差； 在很多情况下地震数据中并没有近场记录，故该方法的使用受到限制。软件模拟方法受软件所使用理论模型、震源类型及环境参数等因素限制，模拟子波与实际子波存在差异。如果虑及检波器特性、仪器因素的影响，从地震数据中直接提取震源信号子波应更合理。随着数字信号处理理论、非线性理论、优化理论的发展，从反射地震数据中提取震源信号子波的方法从理论到实践都取得了巨大进展[19-20]，如何消除反射地震数据中包含的地层信息对提取结果的影响是这类方法成败的关键。……