龚绍京 马 骥 刘双庆 栗连弟

（天津市地震局，天津300201）

对《转换函数与汶川大地震关系的初步研究》一文的分析*

龚绍京 马 骥 刘双庆 栗连弟

（天津市地震局，天津300201）

我们不能苟同该文作者的一些概念，认为该文用的计算公式以及所得的结果存在较多问题，其所得转换函数的量级以及与地震对应关系的分析叫人难以理解。

复转换函数；帕金森矢量；电性结构；复数最小二乘法

1 关于计算方法

关于复转换函数的计算公式，我们发现一些作者的描述概念不清。例如该文［1］及引文［2-4］的作者仅求出了功率谱，利用功率谱求出的振幅谱和（1）式求出了复转换函数或其模|A|、|B|或帕金森矢量系数a，b。仅用功率谱是无法算出复转换函数的，因此，他们求出的不是真正的复转换函数及其模，也不是真正的帕金森矢量系数。

1.1 复转换函数的计算方法

求复转换函数可有两种计算方法：（a）对某一事件三个分量的时间序列，求它们的自谱和互谱，然后利用下面公式求复转换函数：

式中，Phh、Pdd是自谱，即功率谱。Pzh、Pzd、Phd等是互谱。A＝Au＋i Av，B＝Bu＋iBv。下标u代表实部，v代表虚部。功率谱互谱方法由于只用一个事件即可求转换函数，所以要求样本长度较长。一般为24小时，也有用12小时的。即使如此，由于连续的长时间扰动多发生在磁暴期间，源场准均匀的假设很难满足，所以计算结果涨落较大。有人用5个事件的转换函数求平均。我们从对青光台短周期事件的时间序列分析［6］中发现，源场的扰动有一定的“惯性”，因此，认为5个事件太少，结果很难理想。……