安庆地震前地磁台垂直分量的变化特征.docxVIP

安庆地震前地磁台垂直分量的变化特征 1 通过物理和数学的分析,突出的震磁行为 现在，从地磁场变化中提取地震异常的一般方法通常是在掌握地磁场正常变化规律的前提下，根据地磁场的时空变化规律，比较和分析多站、多测点，消除地磁场自身的正常背景和外空场干扰等非地震因素的影响，寻找接近正常变化的因素。通过物理和数学分析，选择有用的地震磁体信号。例如，提取异常的简单值法、复杂值法、时空参考场法、空间相关法和加权差分法。从变化磁体的变化中提取地磁前身信息的振幅相位法、地磁垂直分量的日变形状、地磁性能和释放的响应比法、低位移法和地磁谱分析法。利用地磁场短周期变化的地磁转换函数法和地磁声波谱法进行了比较。 2011年1月19日, 安徽安庆发生了M4.8地震, 同年, 9月10日, 江西瑞昌与湖北阳新交界发生M4.6地震, 两个地震时间间隔不到8个月, 震级相近, 震中距相差190 km。本文利用此两次地震震中周边主要地磁台站绝对地磁观测资料, 采用地磁转换函数、地磁低点位移、垂直分量日变形态三种方法探讨地震前后地震震中周边地磁台地磁场变化特征, 探索地磁场变化与中强地震的孕育和发生的相关性。 2 异常地震数据与地震之间的关系 2.1 静海地磁台转换函数参数的确定 地磁转换函数法可了解和监视震源及邻区的地下电性结构, 如果在地震发生前捕捉到这种变化, 就可进行地震预报。早在上个世纪90年代, 汪江田等就有过利用转换函数法对我国上海及邻近地区中强地震进行预报的尝试, 并确认该方法对中强地震中期预报的有效性。2008年, 李琪等分析了2006年7月4日河北文安地震前后, 距离震中分别为52 km和298 km的静海地磁台和泰安地磁台不同周期的转换函数参数随时间的变化过程, 发现静海地磁台转换函数参数的数据异常点会随着时间的向后推移从短周期向长周期迁移。近年来, 地磁转换函数在我国不同区域、不同省份区域地磁台网均有不同程度的应用研究, 积累了部分地磁转换函数预报地震的宝贵经验和实例。 2.2 与地震活动的关系 地磁低点位移是指地磁垂直分量日变极小值出现时间因受到外来因素的影响偏离正常时间的现象, 研究结果表明, 全国较大范围的地磁低点位移异常与我国中强地震活动有一定的关系。目前研究人员利用地磁日变低点位移的两条以上分界线交汇来研究其与地震之间的关系, 以此缩小预测范围。由于应力的积累, 孕震介质的结构等性质发生变化而产生的各种前兆也呈多样性, 地磁低点位移异常在我国不同地区有着不同的表现形式, 反应不同地区、不同强度的地震信息。 2.3 强震前孕震体整体电导率测量 变化磁场内场强弱既与外场有关, 又与地下物质电导率的分布、大小有关。强震前孕震体环境由于应力增强引起介质电导率等电磁性质的改变, 可观测到与地震有关的局部或区域性磁异常。国内多次震例显示, 中强地震发生前, 区域地磁场垂直分量日变化会出现的相位和幅度变化。 3 安庆地震前后地震资料的特征分析 3.1 异常现象及分析 地磁短周期变化的垂直分量ΔZ、水平分量ΔH和磁偏角ΔD存在如下关系: 式中的A和B为转换函数。转换函数A和B的时空变化特征为地下电性构造的函数, 因而可用于研究局部和部分地区地壳和上地幔电性结构及其变化, 转换函数A和B是地点 (λ, ue788) 、深度 (h) 、频率 (ω) 和时间 (t) 的函数。由于每一个地磁短周期变化都可以对应一个方程 (如式 (1) ) , 于是可以利用多个相同类型的地磁短周期变化建立一个矩阵方程, 用最小二乘法得出转换函数A和B、方位角αp以及总方差σz等。 本方法计算选用湖北、安徽、湖南等震中相邻省份7套FHD质子矢量磁力仪数字化分钟值数据, 资料选用时间段为2009年1月—2012年5月。数据计算前先采用粗差处理技术对数据进行高通滤波, 去除分钟值中数据突跳等不可用数据, 然后, 计算预处理后地磁各要素分钟值虚实富氏谱, 在此基础之上, 计算并选取不同周期的转换函数参数, 最后, 采用滑动平均方法依次去除转换函数中的月以及年变趋势。安庆地震前后, 地磁转换函数出现异常的台站主要有湖北九峰地磁台、安徽泾县地磁台和安徽金寨地磁台, 研究区域内异常地磁台站及显著地震事件震中分布如图1所示, 三个地磁台站泾县台、金寨台、九峰台三个台站距离安徽安庆地震震中依次为120、160、150 km, 距离湖北阳新地震震中依次为310、230、120 km。经台站所属台网人员核实, 异常出现时间段观测仪器的运行率较高, 数据资料稳定可靠。 在安徽安庆4.8级地震前, 2010年3月初开始, 安徽泾县台转换函数A不同周期曲线同步出现异常低值, 异常时间持续6月余, 10月开始, 不同周期异常曲线逐步恢复, 12月底恢复到正常水平 (图2) 。在此期间, 湖北九峰台,