利用电磁台站观测数据推算汶川Ms8.0地震等效电偶矩.pdfVIP

利用电磁台站观测数据推算汶川Ms8．0地震等效电偶矩+ 李美卜2谭捍东I 1j中国地覆大学 j℃哀 2j中营地震台两中心 摘要 本文用单位电流源置于地震震源处产生的电磁场模型，推算2008年5月12日汶川胍8．0 地震前在河北高碑店台观测到电场信号等效电偶矩。计算在观测频率为10Hz时，距震中1000 公里处观测到EW1．3mVm-1电场强度，震源处等效电偶矩为』4l-4．2×1013A．m。研究表明， 地层电导率越大，电场强度随距离变化衰减越快。在现有仪器观测精度0．5mV／mm 合电场强度 0．02mVm-1 下，地层电导率ol lO之S／m时探测范围仅为20km；在6l低于10’4S／m时，能在较 大范围记录到异常信息，例如1500km；另外，在观测平面上，电场分布极不均匀，具有明显方 向性特征，在电流源方向最强，其它方向快速衰减，因此，若观测台站位于非轴向地区，即使在 探测范围内，也记录不到电磁信号。 引言 最近几十年，地震相关电磁辐射现象逐渐被认识和探索，并且已经在地震前兆异常研究中起 et 到重要作用。Li 个超低频 0．1—10Hz 电磁观测台站有3个记录到了电场信息，最大记录幅度达到70mm，合电 场强度为1．3mVm～，这是本台网自上世纪80年代运行观测以来持续时间最长、幅度最大的一次 异常。这期间，中国标准时间2008年5月12日14：28：01分，一次历时Ms8．0地震袭击四川I 汶川地区，震中位于东经103．4。、北纬31．0。，震源深度19kmli】。震源处多大电偶矩才能在距 离震中1300km外的高碑店台产生1．3mVm一1场强?带着这个疑问，本文用一个简单的物理模型 计算不同介质参数下电磁信息衰减和分布特征。 模型描述 为了研究地电流产生的电磁波场，我们延用经Kingetal． 1981 12垸善的置于界面附近一根天 线理论模型，原理见图1，一个X方向单位电流源放置在X 0、Y 0、z d处，z轴正方向向下。 图中，我们假定，一个无限短水平偶极子放置在深度为d的地下半无限空间，填充为介质1 地 球 ，其电性参数为磁导率^、介电常数气、电导率oi：介质2 空气 相应参数为鸬，毛和吒。 如果我们定义复介电常数喜 昌+灯，母，则介质j中的Maxwell方程为： Vx弓 im略 Vx％ 玛一嘞+J 这里，西 2j矿为偶极子圆频率。 经2D空间Fourier变换，我们可以计算在各介质中产生的6个电磁场分量 乓、毛、巨、 墨、当J和芝1。 y 图1地下单位电流源产生电磁场示意图 以汶川地震为例 et 根据Li a1． 2013 描述，考虑主破裂长度AI 240km，设电流源位于震源深度，即 d 19km，观测台站距离震中1000km，所以观测点位于 Ax 1000km，AY 0，△Z O 点。与其 它参数相比，地球电导率吒在计算中起主要作用。对汶川地震，观测频率 厂 10Itz， I 其阈值听由趋肤深度公式占 a厅V力2求出。经计算，在观测点处产生EWl．3mVm-1电场强 度所需电偶矩为 埘Z24．2X1013A．m。 探测能力和电场分布 均Z24．2X 我们固定电流源为常数值 1013A．m，计算不同电性参数 q 下电场强 度随AX变化 图2 。 图2电流源 1．．4／24．2x1013A．m在不同电性参数下产生地表电场强度值IEx 2 由图2，地下电导率越大，电场信号强度衰减越快。如果把探测最小强度设定为0．02mVm一1 这个值接近观测仪器分辨率 ，在Gl 10。3S／m和ol 10。2S／m时，信号探测范围分别为450km 和20km。若GI小于10。4S／m，得到较大的探测范围，如1500km。 芒xVm ■‘m B x日 疆 帅 口 卸 100 ■ o 图3在深度d 19ion电流源』盈22．5x101aA．m产生电场空间分布 a 、 b 、 c 和 d 分别代表z--0、3km、6kin和9km电场分布图 km、 为了研究汶川地震电流源产生的电场在上方沿z轴向的变化，在区域laxI 100lAyl 100 km、0 Az 9 Elx km范围内计算电场分布Elx和Ely。研究发现，在z--0地表面，IElyI比I km 小十几个数量级，z逐渐靠近震源深度，因此，0sAz 9 Elxl起控制作用。图3 a ～ d 为 分别为z Okm、3km、6km和9km电场强度I Elxl分布图。由图3可见，对比不同深度电场变化， 强度随观测点接近电流源而逐渐增强。对每一个xy电场平面，X轴向最强，Y轴向次之，其余非 主轴方向由中心点衰减迅速。这可以部分地解释，尽管位于龙门山断裂带延伸方向的高碑店台距 离汶JII地震震中距离为1300