高斯牛顿迭代法在地震计传递函数计算中应用.pdfVIP

第 期 （总第 期） 四 川 地 震 ! !! 234! 年 月 #$$% ’()*+,-(.’ )’/’()0+ 12 /10+-(2 5678#$$% !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 高斯牛顿迭代法在地震计传递函数计算中的应用 邵玉平 （四川省地震局水库地震研究所，四川 成都 ） %!$$=! 摘要：介绍了利用地震计对阶跃电流标定或稳态正弦电流标定的响应波形获取地震计传递函数的方法。对阶跃 标定电流的响应波形进行傅立叶变换 （ ）和扣除阶跃信号频谱的运算获得地震计的幅频特性；从稳态正弦 * 标定电流的响应波形可直接提取幅频特性，然后用高斯牛顿迭代法拟合出地震计传递函数的各个系数。该方法 的最大优点是只需要已知幅频特性或其有限离散点的数值就可以拟合出地震计的传递函数，较目前常用的需要 同时已知幅频和相频特性的拟合方法更为方便和实用。 关键词：阶跃标定；高斯牛顿迭代；传递函数；地震计 中图分类号： 文献标识码： 文章编号： （ ） *+%#4# ( !$$!:!!9 #$$% $!:$$#;:$9 地震计的工作特性是用传递函数来表达的。如果知道了地震计的传递函数，就知道了地震计线性工 作条件下的特性，如仪器灵敏度、稳态输入响应中的幅频特性、相位特性和群时延特性、瞬态输入响应 中的冲击响应、阶跃响应等。至于地震计内部如何实现该传递函数，对于地震数据使用者来说是不太关 心的。目前测定地震计传递函数主要有两种方法：一是用瞬态信号激励在时域直接用响应波形测定，如 常用的标定线圈阶跃电流标定；二是对上述激励的响应通过变换在频域测定。时域测定法中最常用的是 利用地震计对阶跃脉冲的响应波形，直接测量响应波形的峰值幅度和间隔并计算得到相应的固有周期和 阻尼，从而得到传递函数表达式。这种方法受干扰影响较大，测量阻尼大于 的二阶系统和高于二阶的系 ! 统均有困难。而且测定阻尼周期受采样率和噪声干扰的影响存在较大的误差，所以最后的计算结果误差 也就偏大。 变换测定法仍然是利用地震计对阶跃标定的响应波形，通过傅立叶变换获得阶跃信号的幅频和相频 特性，利用?@A 拟合法 （根据在不同频率下的系统幅频和相频特性进行回归）拟合出传递函数表达式， B 但是?@A 法需要同时已知幅频特性和相频特性，二者不一定都能通过标定方便地得到。例如，在地震部 B 门常用不同频率的稳态正弦激励来标定系统，这样就只能获得幅频特性中数个离散频点上的数值，而得 不到相频特性，难于满足?@A 法同时已知幅频和相频特性的要求。为了解决这个问题，我们采用了高斯 B 牛顿迭代法，只使用幅频特性来计算传递函数。高斯牛顿迭代法的实质是把非线性回归问题逐次化为一 系列线性回归问题来迭代求解。 ! 速度传感地震计的标定 !4! 标定电流与地动速度之间的等效关系 速度传感地震计的标定信号传递框图如图 所示。图中： 为工作线圈的电动常数； 为标定线圈 ! ! ! ! # 的电动常数； 为折合摆长； 为转动惯量； 为拉普拉斯算子。从框图可以得到下面计算公式： # C