台差法提高低信噪比下走时变化探测精度研究.docVIP

台差法提高低信噪比下走时变化探测精度研究 　　摘要：为研究噪音对主动源信号走时变化提取精度的影响，采用对比数值实验的方法。按照常规流程，首先用反褶积去除震源激发时间的不确定性干扰，压缩子波以得到更加尖锐的震相，由于发现反褶积会降低信号的信噪比，于是提出了台差法来代替反褶积矫正气枪激发的不确定性，再结合干涉法提取走时变化，并定量对比了两者的走时变化提取精度，最后尝试用台差法提取主动源走时的日变化。结果表明，台差法可以有效矫正震源激发时刻的影响，且不降低信噪比，从而也能获得较高走时变化提取的精度。 　　关键词：主动源；走时变化；走时台差法；信噪比 　　中图分类号：P315.31文献标识码：A文章编号：1000-0666（2017）04-0595-10 　　0引言 　　地震波信号是照亮地球内部结构的一盏明灯。复杂的地震波信号是由断层的破裂过程和介质格林函数耦合而成。传统的地震学分为研究震源运动学和地球结构研究（Lay，Wallace，1995）。从地震信号研究震源S（t）或者介质结构G（t）就必须使两者解耦。因此，大多数震源或者介质的研究中，都必须假设另一者有先验的模型。 　　由于主动源具有精确的激发时刻、激发地点，为地下介质的研究提供了良好的可重复震源，从而使精细获取介质结构变化成为可能。Niu 等（2008）利用压电陶瓷作为主动源，发现在实验场附近2个小地震发生前，压电陶瓷产生地震波信号走时有显著的延迟。在物理上，由于地震发生之前都会有裂隙和孔隙压力的变化，这个变化在介质波速上有所体现（Simmons，1964；Yukutake，1988），当这种精细变化被主动源信号所捕获时，将有助于我们认识地震孕育的物理过程。压电陶瓷具有频率高、衰减快、传播距离十分有限的特点，不能在大区域的研究中作为广泛应用的主动源。因此，陈?J院士在其所主导建设的宾川地震信号发射台中，引入了4枪组合的气枪激发源作为主动源，从而克服了传播距离有限的问题（罗桂纯等，2006；王宝善等，2011；王伟涛等，2009；杨微等，2013）。 　　目前，宾川发射台气枪实验已顺利进行多年，积累了大量的实验数据。如何提取高精度的主动源信号走时变化是众多学者关心和研究的问题。笔者首先介绍目前提取走时变化常用的干涉法和处理流程，之后对提取走时变化流程中反褶积这一步骤提出可能存在放大噪音而影响走时变化提取精度的问题，并进行相应的数值实验来验证这一想法，进而提出一种可能替代反褶积的台差法来矫正震源的影响，并进行了相应的实证。 　　1干涉法提取走时变化简介 　　波速测量关键的问题在于走时的获取，对于相对走时变化在0.01%～1%量级的精细变化（Wang et al，2008；刘自凤等，2015）探测中，传统手动拾取震相到时的方法已经不再适用。由于主动源气枪激发的重复性，波形上也表现出一定的相似性，基于互相关提取相对走时变化的方法比较适用（刘自凤等，2015；王彬等，2012）。本文以Wang 等（2008）所提出的干涉法为基础，开展信号走时提取。 　　1.1互相关函数 　　对于两个能量有限、定义在[t1，t2]上的有限长地震实信号x（t）和y（t），相关函数定义为：Rxy（τ）=∫t2t1x（t）y（t-τ）dt（1）其中：τ为时间延迟。物理意义上，相关函数表征了将一个信号相对另一个信号移动之后的相似性，为了方便比较，引入无量纲归一化的相关系数ρxy：ρxy（τ）=Rxy（τ）∫t2t1x2（t）∫t2t1y2（t）（2）互相关系数在[-1，1]间变化，互相关系数值越大说明2个信号越相像。2个相同的信号在τ=0时，互相关系数取1。 　　1.2干涉法 　　这里的干涉法特指为气枪设计的一种基于互相关提取走时变化的方法。在气枪主动源的走时变化探测过程中，存在以下2个问题，必须设计与之相应、适合气枪主动源提取走时变化的方法： 　　（1）气枪信号由很多震相组成，不同震相代表了地震波在地下介质中传播的不同路径，因此形成了对地下不同的采样空间。不同路径上的走时变化可能不同，所以有必要截取出一定长度的时窗，使之只包含一个特定的震相，以保证窗内所提取的走时变化是稳定而一致的。 　　地震研究40卷第4期黄亦磊等：台差法提高低信噪比下走时变化探测精度的研究（2）地震仪的采样率为100 Hz，对应的时间采样间隙0.01 s，然而介质速度变化引起的走时变化可能小于时间采样间隙，因此我们有必要对记录的地震信号进行有效插值，将走时测量精度大幅提高至0.001 s以上。 　　干涉法的原理如图1所示。首先我们需要得到一个参考波形，将参考波形和待测信号分成若干个小时窗，每个小时窗可以左右移动一定的步长，将移动过程中最大的互相关系数（图1c）对应的时间移动定为该时窗的相对时移（图1d）