3.11地震基于主成分分析法震中方位角估算.docxVIP

基于主成分分析法的震中方位角估算分析马亮 卢建旗 李山有 马 强 (中国地震局工程力学研究所，中国地震局强震动观测研究室，哈尔滨 150080)摘要：地震定位的实现可以有多种方法，主要分为单台定位法、双台定位法、三台定位法和台网定位。单台定位法用于快速定位，优点是快速及时，能实现地震预警；缺点是精度小误差较大。本文使用主成分分析法对日本311地震多个台站测得的数据计算得出的方位角与实际方位角做比较误差分析，希望对实现快速定位有所裨益。关键词: 地震预警，震中方位角，主成分分析法，地震定位中图分类号：P315.61文献标识码：A文章编号:引言:1988年, 日本铁路技术研究所学者Nakamura在第九届世界地震工程大会上提出了UrEDAS单台地震预警系统，这个系统最初是为了让新干线列车遭遇地震时紧急制动不至于出轨。其中包含了地震的震中方位角确定方法，其原理是使用单个台站接收的最先到达的P波“EW分量到NS分量的振幅比值”导出的估计方位角。这个方法认为，在震源释放的能量的驱使下，地面上的质点M累积的运动轨迹和从震源→质点M的向量所在线段概率性重合，而且质点的累积位移向量与震源→质点向量一致，其的缺陷是它们的方向有可能相反，而且累计计算法得到的震中方位角离散型较大。1992年，吴兆熊、刘希强等的论文里提出了与“基于协方差矩阵的主成分分析法估算震中方位角”的方法有关的数学依据，但并未计算。2006年，刘建华、刘福田、婿颐在其论文里提出了“三分量地震记录的偏振分析”，其中的方法可以用于震中方位角估算。2008年，工程力学研究所学者马强在其博士论文里使用以上方法进行了大量的统计计算，显示了偏振分析法的良好价值。05-07年，李伟，李峰，唐子波等人的主持下完成了《兖州矿区数字遥测矿震台网建设及技术研究》的报告，其中应用了协方差主成分分析法。其主要目的是减轻地震造成的矿难引起的人员伤亡以及财产损失，并做了编程。2010年，山东省地震局学者周彦文、刘希强等再次在论文《基于单台P波记录的快速地震定位方法研究》中明确了协方差矩阵主成分分析法，具有很好的价值。学者发布的论文为震中方位角的估算提供了价值可观的论文与资料，在此基础上，本文对基于强震加速度三分量相关系数矩阵的主成分分析法做以下分析。相关系数矩阵的主成分分析：假如台站接收到的数据（可以使加速度、速度、位移）的垂直、南北、东西三分量存储在列向量（一维数组）D(t)、D2(t)、D3(t)中，令矩阵A=[D(t)，D2(t)，D3(t)]为三维数组。假如强震仪监测到的P波到时为Ka，而且在时间窗[Ka，Ka+F]内P波信息没有收到杂波污染，那么我们就可以对三维数组A1=[D(Ka,Ka+F)，D2(Ka,Ka+F)，D3(Ka,Ka+F)]计算而求得震中方位角，常取时间窗长F∈[0，0.6]秒。计算步骤：对三维数组A1做标准差标准化处理A2=zscore(A1)。并且令A2=(a2ij)3*3，（i，j=1，2，3）用zscore标准化的目的是：使得平均值为0，标准差为1，这样可以使不同量纲的数据放在同一个矩阵里。求数组A2的相关系数矩阵：R3=(rij)3*3rij（i，j=1，2，3）为原变量a2i与a2j的相关系数，其计算公式为：注意到rij=rji，R3是三阶对称方阵。计算R3的特征值与特征向量：用雅可比法（Jacobi）求出特征值，并使其按大小顺序排列λ1≥λ2≥λ3，再求出它们对应的特征向量e1、e2、e3，e1=(e11，e12，e13)，e2=(e21，e22，e23)，e3=(e31，e32，e33)，注意e1、e2、e3的排列顺序对应于λ，而不是从大到小。要求‖ei‖=1，即e2i1+e2i2+e2i3=1　，按这样要求，特征向量就是唯一的。计算震中方位角azi：[2]azi=f{arctan(|e13|/|e12|)}，具体地当e110，e120和e130时,方位角azi=arctan(|e13|/|e12|)当e110，e120和e130时，方位角azi=π-arctan(|e13|/|e12|)当e110，e120和e130时，方位角azi=π+arctan(|e13|/|e12|)当e110，e120和e130时，方位角azi=2π-arctan(|e13|/|e12|)当e110，e120和e130时，方位角azi=π+arctan(|e13|/|e12|)当e110，e120和e130时，方位角azi=2π-arctan(|e13|/|e12|)当e110，e120和e130时，方位角azi=arctan(|e13|/|e12|)当e110，e120和e130时，方位角azi=π-arctan(|e13|/|e12|)震中方位角是经过台站的