Geomagnetic transform function anomaly before the MS5.5 Eryuan earthquake on March 3, 2013
-
摘要: 利用有界影响估计地磁转换函数，分析了2013年3月3日云南洱源MS5.5地震前后滇西北地磁台阵各测点的地磁转换函数及地磁实感应矢量变化特征。结果显示，震前1—2个月云龙测点和洱源测点多个周期的转换函数实部绝对值存在变大—转折—恢复的过程，异常变化幅值超过1倍标准差，其余震中距较大或距离余震区较远的测点震前转换函数无明显异常变化。该结果表明洱源MS5.5地震前地磁转换函数异常范围小于50 km，可能与余震区空间分布有关。Abstract: The geomagnetic transform function represents the lateral inhomogeneity which could be used to study the changes of the underground electricity structure during the seismogenic process. In this study, the geomagnetic transform function and variation characteristics of geomagnetic induction vector for each measuring point of the geomagnetic array in northwestern Yunnan Province before and after the Eryuan MS5.5 earthquake on March 3, 2013 are analyzed by using the bounded influence estimation geomagnetic transform function. The results show that the absolute values of the real part of the transform function for Yunlong station and Eryuan station undergo a process of increase, turn and recover in more than one period during 1−2 months before the earthquake. The amplitude of the abnormal change is more than one time of the standard deviation. Before the earthquake, there is no obvious abnormal change in the transform function of the other measuring points which are far away from the epicenter or theaftershock area. The results show that the abnormal range of the geomagnetic transfer function before the Eryuan MS5.5 earthquake is less than 50 km, which may be related to the spatial distribution of the aftershock area.
-
引言
Benioff （1955）在1952年MW7.5科恩县地震的远场记录中发现破裂传播方向会引起地震动辐射能量随方位角变化的现象。在1957年怀尼米港的M4.7地震中，Housner和Hudson （1958）从距震中西北约8 km的美国海岸及地质调查局（United States Coast and Geodetic Survey）的一个强震仪的记录中观察到了第一个单脉冲强地震动，提醒工程师注意这类来自小震级地震的强烈地震作用所造成的工程结构的严重破坏。Housener和Trifunac （1967）在1966年帕克菲尔德M6.0地震的2号台站加速度记录中，首次认识并解释了近断层方向性效应引起的速度大脉冲（Aki，1967，1968），记录的速度脉冲幅值达到73 cm/s。在随后发生的1994年美国北岭市地震、1999年台湾集集地震和2011年日本3·11地震的近断层强震记录中多次观察到速度脉冲。方向性效应引起的速度脉冲型地震动以短持时、高能量为特征，工程结构，尤其是长周期结构，在其作用下会产生较大的位移和变形，故而引起了工程界的广泛关注。Anderson和Bertero （1987）的研究表明，如果脉冲周期与结构的自振周期相近，结构损伤更严重。脉冲型地震动的频谱形状与非脉冲型地震动相比有很大差别，对结构抗震设计有更严格的要求（Iervolino，Cornell，2008）。
一般来说，当观测台站位于破裂方向且破裂速度接近传播介质的剪切波速时，破裂方向性效应会引起速度大脉冲。对于走滑断层，滑冲效应引起的单侧速度大脉冲主要出现在断层平行向，常伴随着地表永久位移；破裂方向性效应引起的双侧速度大脉冲主要出现在断层垂直向，一般在强地震动的初始阶段出现，不一定造成永久位移。Somerville等（1997）和Somerville （2003）的研究表明，长周期脉冲型地震动不仅会造成地表永久位移，还会造成场地地震动幅值、持时等的空间差异。在高频部分，方向性效应令破裂方向上断层远端的加速度记录强震段持时较短且幅值较大，破裂起始点加速度记录的幅值较低、持时较长；低频部分，破裂方向性引起长周期地震波相长干涉和大幅值相干速度脉冲。刘启方（2005）在均匀弹性全空间以及基岩半空间内讨论了方向性产生的速度脉冲和震源参数之间的关系，表明断层的埋深、震级大小、初始破裂点和凹凸体的位置都会对脉冲的周期和幅值产生影响。Fayjaloun等（2017）的研究表明，断层的几何形状和破裂速度是影响脉冲持续时间的主要因素。Poiata等（2018）分析倾滑断层中脉冲型地震动的产生机制，发现下盘台站中，方向性效应对脉冲的产生起主导作用，而在上盘台站中，由Kagawa （2009）所提出的聚集效应是脉冲产生的原因。Scala等（2018）通过模拟2009年意大利拉奎拉地震的地震动，发现脉冲的出现和持续时间不完全取决于震源上升时间，与平均上升时间、台站位置、破裂速度和断层埋深关系较大，破裂面粗糙度对其影响较小。罗全波（2019）模拟了1999年台湾集集地震中的速度脉冲型地震动，发现震源模型中凹凸体埋深和上升时间对速度脉冲的波形、周期和幅值均有贡献。具体表现是，速度脉冲周期与凹凸体的上升时间成正比，速度脉冲的幅值与凹凸体上升时间呈反比。曹泽林（2020）借助频率-波数域格林函数（FK法）合成了1979年帝王谷MW6.5地震和1994年北岭市MW6.7地震的地震动场，基于等时线理论解释了速度脉冲的形成机理，认为震源参数、地壳结构的剪切波速、不连续界面的位置、界面两侧速度差都会影响脉冲区域，错动量分布对脉冲区域有显著影响。
本文以1979年帝王谷地震为例，分别用三种方法识别近断层强震观测记录的速度脉冲。随后，借助基于FK法合成的地震动，通过观察子源地震动叠加过程中的频率成分变化，速度脉冲出现的时间和凹凸体与台站的空间关系，初步探讨近断层脉冲型地震动的产生机制。
1. 强震记录的速度脉冲识别
1.1 研究数据
1979年帝王谷MW6.5地震发生在美国与墨西哥交界的圣安德烈斯断裂带附近，震中为（32.64°N，115.30°W）位于地震多发和重点观测地区，震源深度为8 km。这是一次走滑地震，破裂始于帝王断层并触发了布劳利断层。前者走向为323°、倾角为80°，在帝王谷断层产生了长达30.5 km的地表破裂；后者走向为360°、倾角为90°，产生总长13 km不连续的地表破裂（Archuleta，1984）。在近断层加速度时程中发现PGA的空间分布与破裂方向相关，呈现明显的方向性特征。本文从美国太平洋地震工程研究中心（Pacific Earthquake Engineering Research Center，缩写为PEER）收集断层距20 km内31个台站东西向和南北向两个水平分量的加速度记录，震中、断层地表投影和台站位置示于图1a，加速度时程见图1b。
![]()
图 1 帝王谷地震震中、断层投影和台站位置 （a）及31个台站的加速度时程 （b）Figure 1. The location of epicenter, fault projection and stations （a） and acceleration records at 31 stations （b） of the Imperial Valley earthquake1.2 速度脉冲识别方法及在帝王谷地震中的应用
近十余年出现了一些脉冲特征识别的方法，大体上可以分为两类：一类从地震学的角度出发，关注在地震记录早期出现的速度脉冲，认为其单纯由断层破裂引起；另一类面向工程抗震设计，不论脉冲是如何及何时产生，在整个地震动时程中识别脉冲，考虑其对工程结构的影响。第一类方法以Baker （2007）及Shahi和Baker （2014）提出的小波变换识别方向性效应产生的脉冲为代表。选取4阶Daubechies小波作为母波，对速度时程进行连续小波变换，绝对值最大的小波系数用来确定脉冲的周期和位置，第一次实现了通过脉冲指标量化识别速度脉冲型地震动。对美国西部地区下一代衰减关系（Next Generation Attenuation-West2，缩写为NGA-West2）的398条脉冲型地震动进行分类和统计，Baker （2007）给出了脉冲型地震动识别指标的表达式为
$$ {\rm{PI}} {\text{＝}} \frac{{\text{1}}}{{{\text{1}} {\text{＋}} {{\text{e}}^{ {\text{－}} 23.3 {\text{＋}} 14.6（{\rm{PGV\;ratio}}）{\text{＋}} 20.5（{\rm{enery\;ratio}}）}}}} {\text{，}} $$ (1) 式中，PGV ratio表示提取脉冲后残余地震动的PGV与原始地震动记录PGV之比；energy ratio表示提取脉冲后残余地震动与原始地震动记录能量之比。图2是借助该方法识别出脉冲型地震动的一个例子，识别标准为：① 脉冲指标PI≥0.85；② 前方向性效应引起的脉冲型地震动中，为排除晚到或非方向性效应产生的脉冲，要求速度时程记录前20%能量所在时刻迟于脉冲记录前10%能量所在时刻；③ 原始地震动的PGV＞30 cm/s。
![]()
图 2 小波变换法识别出的速度脉冲型地震动Figure 2. Identification of pulse-like ground motion by wavelet transform小波变换方法可以识别出单脉冲型地震动，但反复进行小波变换令计算效率不高（Mukhopadhyay，Gupta，2013），脉冲指标PI≥0.85这一标准可能会遗漏部分PI值处于0.15—0.85之间含有脉冲的速度记录（Mena，Mai，2011）。为克服这些局限性，Zhai等（2013）认为在短时间内能量骤增的脉冲型地震动对工程结构的危害更严重，从工程学角度提出了基于能量的识别和提取方法。该方法采用Dickinson和Gavin （2011）的脉冲模型，所用脉冲拟合函数能有效识别不同峰值、形状和宽度的脉冲，其表达式为
$$ {V_{\rm{p}}}（ {t；{V_{\rm{p}}}，{T_{\rm{p}}}，{N_{\rm{c}}}，{T_{{\rm{pk}}}}，\varphi } ） {\text{＝}} {V_{\rm{p}}}\exp \left[ { {\text{－}} \frac{{{\pi ^2}}}{4}{{\left( {\frac{{t {\text{－}} {T_{{\rm{pk}}}}}}{{{N_{\rm{c}}}{T_{\rm{p}}}}}} \right)}^2}} \right] \cos \left( {2\pi \frac{{t {\text{－}} {T_{{\rm{pk}}}}}}{{{T_{\rm{p}}}}} {\text{－}} \varphi } \right){\text{，}}$$ (2) 式中，Vp为拟合脉冲峰值，Tp为拟合脉冲周期，Nc为拟合脉冲个数，Tpk为脉冲峰值时刻，φ为脉冲相位角。
Zhai等（2013）认为保留Dickinson和Gavin （2011）模型中模拟脉冲的几个关键参数，如脉冲周期和脉冲峰值等，可让模型计算能力得到大幅度提升，故将模型进行了简化：由于只考虑单脉冲型地震动，令表示半循环脉冲个数的Nc＝1，半循环脉冲的相位为0。简化这两个参数后，利用式（2）与速度时程信号进行最小二乘拟合得到最匹配的脉冲峰值Vp和对应的周期Tp。脉冲周期的研究初期，Osaki （1994）考虑到地震动速度时程曲线很不规则，无法准确地定义脉冲的起点和终点，故而提出了峰点法（peak point method，缩写为PPM）。脉冲周期的值定义为与PGV相邻的波峰或者波谷的时间间隔。Zhai等（2013）根据PGV确定脉冲幅值Vp，以脉冲相对能量为判别指标，范围在0—1之间，超过阈值即为脉冲型地震动。为了避免将幅值较低或距离断层较远的地震动判别为脉冲型，与Baker （2007）相同，用PGV大于30 cm/s作为标准，从NGA数据库3 635条记录中筛选出357条脉冲型地震动记录。计算脉冲相对能量Ep，将Ep阈值设定为0.3，即Ep≥0.3的地震动被判定为脉冲型。图3为该方法识别的一例脉冲型地震动。
小波分析和能量法是识别单脉冲型地震动的两种方法。但方向性效应引起的脉冲一般包含主、次脉冲，主、次脉冲的周期一般相差较小。一个有效的脉冲识别方法应该能够识别出多脉冲型地震动，Zhao等（2016）所选用的零点法（zero velocity point method，缩写为ZVPM）提取速度脉冲的核心思想是通过定义速度时程中循环振动的区间长度作为脉冲周期，针对不同脉冲型记录定义不同的判别标准。通过识别地震动速度时程中不同循环振动区间的地震动参数获得等效速度脉冲。若识别的脉冲型地震动只有一个速度脉冲，则脉冲的能量比阈值为0.3；对包含多个脉冲的地震动，分别确定能量比为0.6，0.8和0.9作为双脉冲、三脉冲和四脉冲型地震动的阈值。综上，Baker （2007）及Shahi和Baker （2014）提出的小波变换法通过小波变换提取速度脉冲，第一次以客观指标识别脉冲型地震动；Zhai等（2013）提出的峰点法以Dickinson和Gavin （2011）模型为基础建立能量指标，在脉冲周期的识别上有明显进步；Zhao等（2016）用零点法确定脉冲周期的判别方法，给出识别多脉冲型地震动的阈值。图4为采用三种方法识别脉冲型地震动的一个例子。
![]()
图 4 三种方法提取RSN158台站记录的脉冲Figure 4. Identification of pulse-like ground motion by three methods以帝王谷地震为例，比较用这三种方法识别脉冲参数的差别。首先，将31个台站东西向和南北向两个水平分量的强震记录通过式（3）旋转为垂直断层方向FN和平行断层方向FP （Somerville et al，1997），即
$$ \left\{\begin{gathered} V_{\rm{FN}} {\text{＝}} V_{\rm{EW}} {\text{sin}}\phi {\text{＋}} V_{\rm{NS}}{\text{cos}} \phi {\text{，}} \hfill \\ V_{\rm{FP}} {\text{＝}} V_{\rm{EW}} {\text{cos}}\phi {\text{－}} V_{\rm{NS}}{\text{sin}} \phi {\text{．}} \end{gathered} \right. $$ (3) 小波变换法、峰点法、零点法三种方法分别识别出12条、12条、13条脉冲型地震动。三种方法共同识别出在FN方向上含有脉冲的台站所提取出的脉冲周期Tp和脉冲峰值Vp列于表1。
表 1 三种方法识别出的帝王谷地震中脉冲型地震动的脉冲参数Table 1. The pulse indicators of pulse-like ground motions in the Imperial Valley earthquake identified by three methodsRSN台站编号 脉冲周期Tp/s 脉冲峰值Vp/（cm·s−1） 小波变换法 峰点法 零点法 小波变换法 峰点法 零点法 158 − 1.9 1.9 − 59 60 159 2.3 1.9 1.9 53 50 53 161 4.3 4.2 3.3 51 55 52 170 4.4 3.9 3.5 79 63 71 171 3.4 3.2 3.1 121 111 115 173 4.5 3.3 4.9 69 71 70 178 4.5 4.4 4.4 69 61 60 179 4.8 4.4 3.8 81 89 82 180 4.1 3.9 2.1 130 133 136 181 3.8 3.5 4.8 155 156 160 182 4.4 3.1 4.4 98 100 109 185 2.8 3.5 3.0 68 71 72 三种方法识别出的脉冲峰值基本一致，大多数差别小于5%；脉冲周期因数学模型及算法的不同有一定差别，大多数差别大于10%，零点法的结果一般小于其它两种方法的结果。图5给出了三种方法识别出的含有脉冲的台站与断层地表投影的空间关系。
![]()
图 5 小波变换（a）、峰点法（b）和零点法（c）识别出的速度脉冲的台站Figure 5. Stations with velocity pulse identified by wavelet transform （a），PPM （b），and ZVPM （c）可见，三种方法识别出的脉冲记录台站主要分布在垂直于断层的方向上；小波变换法主要识别出的是近断层区域的速度脉冲，RSN160台站位于震中附近；峰点法和零点法可以识别出距断层较远的RSN175台站的脉冲型地震动。这是由前述两种方法的出发点不同所引起的。
2. 合成地震动的速度脉冲识别
2.1 基于FK法合成的宽频带地震动
本文采用的基于FK法的宽频带地震动场合成方法（曹泽林，2020）是在频率波数域求解格林函数结合有限断层震源模型的方法。该方法合成地震动的带宽主要取决于震源模型，高频分量主要由震源时间函数和上升时间控制，低频分量主要由错动分布控制，破裂速度的不均匀性和平均值也有影响。在有限断层震源模型中，断层多取为规则形状的破裂面，将其划分为NL×NW个沿走向和倾向相同大小的子源。从破裂起始点后破裂的能量以一定的速度传播，每个子源引起的地表场地上的地震动根据子源与场地的几何关系和区域地壳速度结构计算，各子源引起的地震动按一定的时滞进行叠加（式（4）），最后得到场地的地震动。
$$ a（t） {\text{＝}} \sum\limits_{i {\text{＝}} 1}^{{N_{\rm{L}}}} {\sum\limits_{j {\text{＝}} 1}^{{N_{\rm{W}}}} {{a_{ij}}［（t {\text{－}} （\Delta {t_{ij}'} {\text{＋}} \Delta {t_{ij}''}））］} } ， $$ (4) 式中，t为破裂时间，
$\Delta {t_{ij}'}$ 为第i个子源的破裂时间，$\Delta {t_{ij}''}$ 为第ij个子源到地表点的传播时间。本文以帝王谷地震为例，借助该方法合成的地震动，从时间和空间上分析震源破裂过程对脉冲型地震动的贡献。震源参数参考了Hartzell和Heaton （1983）的成果，详细信息列于表2。
表 2 帝王谷地震的震源参数Table 2. Source parameters of Imperial Valley earthquake震源
机制走向/° 倾角/° 断层参数 子源个数 凹凸体 埋深
/km长度
/km宽度
/km面积
/km2沿走向 沿倾向 个数 面积
/km2长度
/km宽度
/km中心沿走向、
倾向位置/km走滑型 323 90 8 42 10.4 436.8 16 8 1 102 13 7.8 22，6.4 断层面上的错动量和滑动角分布来自有限断层破裂模型数据库SRCMO （图6）。破裂面上仅有一个凹凸体，位于断层中部，破裂起始点如图中矩形所示，沿破裂面向右上方逐渐破裂，耗时13 s左右。
![]()
图 6 帝王谷地震断层面上的错动分布Figure 6. Slip distribution on the fault plane of the Imperial Valley earthquake2.2 借助子源叠加过程分析脉冲机制
现有的研究多关注于凹凸体的个数和位置、破裂速度、破裂方式以及上升时间对脉冲型地震动的影响（Somerville et al，1997；Somerville，2003）。曹泽林（2020）的研究表明，除了震源参数和地壳结构等影响因素外，其它因素也可能引起脉冲型地震动，例如，某些场地位于破裂面的凹凸体附近，或地震波穿过复杂地壳结构。Lin （2020）通过对台湾花莲地震中脉冲型地震动的前向建模的结果表明，速度脉冲的特征受米伦断层中的一个子断层影响。为探究凹凸体和台站空间位置等与脉冲产生机制的关系，本文以曹泽林（2020）基于FK方法合成的31个帝王谷台站的地震动为基础，在含有脉冲且位于破裂前方的10个台站中识别各个子源地震动叠加生成各个地表台站地震动时程过程中的速度脉冲，分析子源地震动叠加过程对形成速度脉冲的影响。由于篇幅所限，本文给出台站RSN158 （图7a），RSN179的例子（图7b），台站RSN158位于震中附近，凹凸体破裂的前方，台站RSN179位于凹凸体破裂后方且垂直于断层的方向。图中，阴影部分表示在叠加过程中识别出脉冲对应的子源，虚线圈出的长方形为凹凸体位置。在台站RSN158地震动的叠加过程中，在破裂起始点附近的子源就产生脉冲直至破裂完成均识别出脉冲；对于RSN179台站的地震动，叠加过程中第一次出现脉冲的子源位于凹凸体内部，随后破裂的所有子源叠加的地震动中均识别出脉冲。
![]()
图 7 台站RSN158 （a）和台站RSN179 （b）在破裂过程中出现脉冲的子源Figure 7. The velocity pulse appears during rupture process at station RSN158 （a） and station RSN179 （b）由图7可知，脉冲产生的早晚以及首次出现脉冲的子源位置和台站与凹凸体的空间关系有关。对所有台站均识别各个子源地震动叠加生成地震动时程过程中的速度脉冲，共有10个台站的实际地震记录和合成记录中都识别出了脉冲。从脉冲初次出现对应的子源、凹凸体和台站的空间位置关系，分析凹凸体对形成速度脉冲的影响（图8）。图中，红色矩形为首先出现脉冲对应的子源，虚线圈出的长方形为凹凸体位置。
![]()
图 8 最初产生脉冲的子源位置与凹凸体和错动量的关系Figure 8. The relationship between the sub-source position that the pulsefirst generated， the asperity and the slip distribution可见，在震中附近台站（RSN158，RSN159和RSN160），脉冲初次出现时对应的子源位置在破裂起始点周围。随着子源时程的不断叠加，当叠加到凹凸体内子源时，沿破裂方向、距震中较远台站（RSN171，RSN179，RSN180，RSN181，RSN182，RSN183，RSN184）的合成地震动中识别出脉冲，直至破裂完成，这些台站的合成地震动中始终能识别出脉冲。而脉冲的产生在子源地震动不断叠加的过程中得以体现，根据蒋连接和白国良（2016）的研究，脉冲型地震动的能量主要集中在低频部分。通过小波变换，给出台站RSN158和RSN179在叠加过程中不同子源地震动的时频分析图 （图9）。
从台站RSN158的地震动叠加过程可见，自破裂起始点开始，较高频率能量的释放转为低频能量集中，持续时间较短，通常单脉冲的脉冲周期大于0.5 s，在叠加到最初产生脉冲的子源位置时，低频能量的持时变长，低频脉冲能量明显，判定为速度脉冲。而RSN179台站，在叠加到凹凸体内部的子源时，伴随着长持时和显著的低频脉冲能量，速度脉冲初次产生。从位于凹凸体破裂前后方台站地震动的叠加过程来看，首次出现脉冲的子源以及相应的低频分量含量表明凹凸体与台站的空间关系对脉冲的产生有很大影响，这与上述曹泽林（2020）和Lin （2020）的结论一致。
3. 讨论与结论
本文以帝王谷地震为例，分别用三种方法识别脉冲型强震记录，借助FK法合成的地震动，通过逐步判别各子源地震动叠加过程中合成地震动的速度脉冲，探讨近断层速度脉冲的产生机制，得到以下几点初步结论：
1） 借助小波分析法、峰点法和零点法分别从31条近断层强震记录识别速度脉冲，三种方法提取的脉冲峰值差别较小，脉冲周期差别较大，零点法提取的脉冲周期整体上小于其它两种方法的结果。
2） 基于FK法的地震动合成方法能有效地表达速度脉冲且可以合成多脉冲地震动，而且合成结果包含震源、地壳速度结构的复杂影响，为进一步分析速度脉冲的产生机制奠定坚实基础。
3） 借助FK法合成的子源地震动记录进行叠加，结合凹凸体与台站的空间位置发现，位于凹凸体后方且垂直于断层的台站多数含有速度脉冲；在子源的叠加过程中，当脉冲初次出现对应的子源位于破裂起始点附近，在震中附近台站形成脉冲；位于凹凸体内部，在沿破裂方向上形成速度脉冲；凹凸体的位置对速度脉冲的形成有很大影响。
本文所用强震数据来自美国太平洋地震工程研究中心（The Pacific Earthquake Engineering Research Center），小波变换程序来自Baker （2007），西南交通大学讲师常志旺和哈尔滨工业大学讲师赵国臣提供了峰点法和零点法的计算程序，作者在此一并表示感谢。
-
![]()
图 1 滇西北地磁台阵测点分布及洱源两次MS≥5.0地震的震中和余震分布图
F1：怒江断裂；F2：澜沧江断裂；F3：维西—乔后断裂；F4：龙蟠—乔后断裂；F5：鹤庆—洱源断裂； F6：小金河—丽江断裂；F7：红河断裂；F8：程海—宾川断裂
Figure 1. Distribution of measuring points of geomagnetic array in northwest Yunnan Province and epicenters and aftershocks of two MS≥5.0 earthquakes in Eryuan
F1：Nujiang fault；F2：Lancangjiang fault；F3：Weixi-Qiaohou fault；F4：Longpan-Qiaohou fault；F5：Heqing-Eryuan fault； F6：Xiaojinhe-Lijiang fault；F7：Honghe fault；F8：Chenghai-Binchuan fault
![]()
图 2 不同求取方法所得地磁转换函数的稳定性对比
（a） 在水平分量中加入5%高斯白噪声的转换函数实部Ar；（b） 在垂直分量中加入5%高斯白噪声的转换函数实部Ar
Figure 2. Comparison of geomagnetic transfer function stability estimated by three different methods of BI estimator，least squares and weighted least squares
（a） Real component Ar of geomagnetic transfer function result with 5% Gaussian white noise in horizontal component；（b） Real component Ar of geomagnetic transfer function result with 5% Gaussian white noise in vertical component
![]()
图 3 云龙台不同周期T转换函数实部Ar的时间变化曲线
Figure 3. The time curves of the real component Ar of transfer function for different periods T at the Yunlong station
![]()
图 4 云龙台不同周期T转换函数实部Br的时间变化曲线
Figure 4. The time curves of the real component Br of transfer function at the Yunlong station for different periods T
![]()
图 5 周期T为18.3 min时各测点转换函数实部Br的时间曲线
Figure 5. The time curves of Br at the period of 18.3 minutes for five geomagnetic stations
![]()
图 6 不同周期实感应矢量随时间的变化
Figure 6. The real induction vector varies with time at different periods
（a） T＝18.3 min；（b） T＝25.6 min；（c） T＝36.6 min；（d） T ＝51.2 min
![]()
图 7 测点周期为18.3 min时实感应矢量长度L的时间曲线
Figure 7. The time curves of L at the period of 18.3 minutes for five geomagnetic stations
表 1 滇西北地磁台阵各测点与洱源两次MS≥5.0地震震中的距离
Table 1 The distances between the epicenters of two MS≥5.0 earthquakes in Eryuan and the measuring points of the geomagnetic array in northwest Yunnan Province
测点名称 与洱源MS≥5.0地震震中的距离/km 测点名称 与洱源MS≥5.0地震震中的距离/km MS5.5 MS5.0 MS5.5 MS5.0 云龙 37.15 42.85 牛街 47.12 50.29 剑川 52.73 52.55 鹤庆 76.53 77.34 炼铁 8.34 11.52 下关 64.50 55.88 洱源 22.94 27.38 丽江 116.05 118.71 
下载: 导出CSV
-
常祖峰,张艳凤,周青云,虎雄林,臧阳. 2014. 2013年洱源MS5.5地震烈度分布及震区活动构造背景研究[J]. 中国地震,30(4):560–570. doi: 10.3969/j.issn.1001-4683.2014.04.009 Chang Z F,Zhang Y F,Zhou Q Y,Hu X L,Zang Y. 2014. Intensity distribution characteristics and active tectonic background in area of the 2013 Eryuan MS5.5 earthquake[J]. Earthquake Research in China,30(4):560–570 (in Chinese).
车用太, 鱼金子. 2006. 地震地下流体学[M]. 北京: 气象出版社: 205–208. Che Y T, Yu J Z. 2006. Underground Fluids and Earthquake[M]. Beijing: China Meteorological Press: 205–208（in Chinese）.
陈伯舫. 1974. 渤海西岸的电导率异常[J]. 地球物理学报,17(3):169–172. Chen P F. 1974. Conductivity anomaly in west coast of Po Hai[J]. Acta Geophysica Sinica,17(3):169–172 (in Chinese).
丁鉴海,申旭辉,潘威炎,张晶,余素荣,李纲,关华平. 2006. 地震电磁前兆研究进展[J]. 电波科学学报,21(5):791–801. doi: 10.3969/j.issn.1005-0388.2006.05.031 Ding J H,Shen X H,Pan W Y,Zhang J,Yu S R,Li G,Guan H P. 2006. Seismo-electromagnetism precursor research progress[J]. Chinese Journal of Radio Science,21(5):791–801 (in Chinese).
高琼,金明培,邓嘉美. 2014. 2013年洱源5.5级地震前地震活动异常分析[J]. 地震地磁观测与研究,35(5/6):68–71. Gao Q,Jin M P,Deng J M. 2014. Analysis of the seismicity anomalies before Eryuan M5.5 earthquake in 2013[J]. Seismological and Geomagnetic Observation and Research,35(5/6):68–71 (in Chinese).
龚绍京,吴占峰. 1986. 唐山地震可能伴随的地电导率变化[J]. 地震学报,8(1):28–36. Gong S J,Wu Z F. 1986. Possible changes of earth conductivity accompanying the Tangshan earthquake of 1976[J]. Acta Seismologica Sinica,8(1):28–36 (in Chinese).
龚绍京,刘双庆,张明东. 2015. 2007年5月—2013年12月成都、西昌和重庆台地磁转换函数的时间变化[J]. 地震学报,37(1):144–159. doi: 10.11939/j.issn:0253-3782.2015.01.013 Gong S J,Liu S Q,Zhang M D. 2015. Time changes in the geomagnetic transfer functions at Chengdu,Xichang and Chongqing stations from May 2007 to December 2013[J]. Acta Seismologica Sinica,37(1):144–159 (in Chinese).
郝锦绮,黄平章,张天中,徐明发,周建国,马醒华. 1989. 岩石剩余磁化强度的应力效应[J]. 地震学报,11(4):381–391. Hao J Q,Huang P Z,Zhang T Z,Xu M F,Zhou J G,Ma X H. 1989. The stress effect on remanent magnetization of rocks[J]. Acta Seismologica Sinica,11(4):381–391 (in Chinese).
侯作中,史铁生. 1984. 云南省的地磁短周期变化异常[J]. 地震学报,6(3):287–293. Hou Z Z,Shi T S. 1984. Short-period geomagnetic variation anomaly in Yunnan Province[J]. Acta Seismologica Sinica,6(3):287–293 (in Chinese).
黄小龙,吴中海,蒋瑶,田婷婷,刘艳辉. 2015. 2013年3月3日云南大理洱源MS5.5地震烈度分布及发震构造[J]. 地质通报,34(1):135–145. doi: 10.3969/j.issn.1671-2552.2015.01.011 Huang X L,Wu Z H,Jiang Y,Tian T T,Liu Y H. 2015. Seismic intensity distribution and seismogenic structure analysis of the March 3,2013 Eryuan MS5.5 earthquake in Dali,Yunnan Province[J]. Geological Bulletin of China,34(1):135–145 (in Chinese).
钱复业,赵玉林,于谋明,王志贤,刘小伟,常思敏. 1982. 地震前地电阻率的异常变化[J]. 中国科学:B辑,33(9):831–839. Qian F Y,Zhao Y L,Yu M M,Wang Z X,Liu X W,Chang S M. 1982. Earth-resistivity anomalies before earthquake[J]. Scientia Sinica:Series B,33(9):831–839 (in Chinese).
孙洁,徐常芳,江钊,史书林,王继军,何明. 1989. 滇西地区地壳上地幔电性结构与地壳构造活动的关系[J]. 地震地质,11(1):35–45. Sun J,Xu C F,Jiang Z,Shi S L,Wang J J,He M. 1989. The electrical structure of the crust and upper mantle in the west part of Yunnan Province and its relation to crustal tectonics[J]. Seismology and Geology,11(1):35–45 (in Chinese).
王辉,许滔滔,罗景程,曾雅俊,魏文博. 2018. 基于地磁台网数据的长周期大地电磁数据远参考处理[J]. 地球物理学进展,33(6):2270–2277. doi: 10.6038/pg2018BB0446 Wang H,Xu T T,Luo J C,Zeng Y J,Wei W B. 2018. Remote reference processing of long period magnetotelluric data based on geomagnetic data[J]. Progress in Geophysics,33(6):2270–2277 (in Chinese).
王桥,黄清华. 2016. 华北地磁感应矢量时空特征分析[J]. 地球物理学报,59(1):215–228. doi: 10.6038/cjg20160118 Wang Q,Huang Q H. 2016. The spatio-temporal characteristics of geomagnetic induction vectors in North China[J]. Chinese Journal of Geophysics,59(1):215–228 (in Chinese).
吴乾蕃,祖金华,谢毅真,王都. 1988. 云南地区地热基本特征[J]. 地震地质,10(4):177–183. Wu Q F,Zu J H,Xie Y Z,Wang D. 1988. Characteristics of geothermal field in Yunnan region[J]. Seismology and Geology,10(4):177–183 (in Chinese).
邢西淳,刘春,毛娟,邵辉成. 2008. 关中地区地磁转换函数变化与地震关系初探[J]. 高原地震,20(2):11–16. doi: 10.3969/j.issn.1005-586X.2008.02.003 Xing X C,Liu C,Mao J,Shao H C. 2008. Discussion on the relationship between geomagnetic transfer characteristic variation and earthquakes in Guanzhong and its neighbouring area[J]. Plateau Earthquake Research,20(2):11–16 (in Chinese).
徐文耀,祁骙,王仕明. 1978. 甘肃省东部地区短周期地磁变化异常及其与地震的关系[J]. 地球物理学报,21(3):218–224. Xu W Y,Qi K,Wang S M. 1978. On the short period geomagnetic variation anomaly of the eastern Kansu Province[J]. Acta Geophysica Sinica,21(3):218–224 (in Chinese).
袁伊人,李琪,朱培育. 2015. 滇西北地区地磁短周期变化特征研究[J]. 地震地磁观测与研究,36(3):77–82. doi: 10.3969/j.issn.1003-3246.2015.03.015 Yuan Y R,Li Q,Zhu P Y. 2015. The geomagnetic short period variations in the northwestern Yunnan[J]. Seismological and Geomagnetic Observation and Research,36(3):77–82 (in Chinese).
曾小苹,林云芳. 1995. 地磁短周期变化异常对中国中强地震的响应[J]. 地震,(1):29–36. Zeng X P,Lin Y F. 1995. Response of anomalous short periodic geomagnetic variations to moderately strong earthquakes in China[J]. Earthquake,(1):29–36 (in Chinese).
Beamish D. 1982. A geomagnetic precursor to the 1979 Carlisle earthquake[J]. Geophys J Int,68(2):531–543. doi: 10.1111/j.1365-246X.1982.tb04913.x
Chave A D,Thomson D J. 2003. A bounded influence regression estimator based on the statistics of the hat matrix[J]. J R Stat Soc,52(3):307–322. doi: 10.1111/1467-9876.00406
Chave A D,Thomson D J. 2004. Bounded influence magnetotelluric response function estimation[J]. Geophys J Int,157(3):988–1006. doi: 10.1111/j.1365-246X.2004.02203.x
Egbert G D,Booker J R. 1986. Robust estimation of geomagnetic transfer functions[J]. Geophys J Int,87(1):173–194. doi: 10.1111/j.1365-246X.1986.tb04552.x
Eisel M,Egbert G D. 2001. On the stability of magnetotelluric transfer function estimates and the reliability of their variances[J]. Geophys J Int,144(1):65–82. doi: 10.1046/j.1365-246x.2001.00292.x
Gamble T D,Goubau W M,Clarke J. 1979. Error analysis for remote reference magnetotellurics[J]. Geophysics,44(5):959–968. doi: 10.1190/1.1440988
Nur A. 1972. Dilatancy,pore fluids,and premonitory variations of tS/tP travel times[J]. Bull Seismol Soc Am,62(5):1217–1222.
Parkinson W D. 1959. Directions of rapid geomagnetic fluctuations[J]. Geophys J Int,2(1):1–14. doi: 10.1111/j.1365-246X.1959.tb05776.x
Rikitake T. 1979. Changes in the direction of magnetic vector of short-period geomagnetic variations before the 1972 Sitka,Alaska,earthquake[J]. J Geomagn Geoelect,31(4):441–448. doi: 10.5636/jgg.31.441
Schmucker U. 1970. Anomalies of Geomagnetic Variations in the Southwestern United States[M]. Berkeley: University of California.
Scholz C H,Sykes L R,Aggarwal Y P. 1973. Earthquake prediction:A physical basis[J]. Science,181(4102):803–810. doi: 10.1126/science.181.4102.803
Sims W E,Bostick Jr F X,Smith H W. 1971. The estimation of magnetotelluric impedance tensor elements from measured data[J]. Geophysics,36(5):938–942. doi: 10.1190/1.1440225
Thomson D J. 1982. Spectrum estimation and harmonic analysis[J]. Proc IEEE,70(9):1055–1096. doi: 10.1109/PROC.1982.12433
Vargas J A,Ritter O. 2016. Source effects in mid-latitude geomagnetic transfer functions[J]. Geophys J Int,204(1):606–630.
Yanagihara K. 1972. Secular variation of the electrical conductivity anomaly in the central part of Japan[J]. Mem Kakioka Mag Obs,15(1):1–11.
-
期刊类型引用(15)
1. 刘文玉，程正璞，年秀清，陈闫，胡钰铃，覃祖建，邵明正. 基于三维剩余密度结构的松原地震成因. 地震地质. 2024(02): 462-476 . 
百度学术
2. 于晨，卢军，解滔，刘长生. 2017—2019年松原地区4次中强地震前绥化台地电阻率异常分析. 中国地震. 2023(02): 314-324 . 百度学术
3. 张志朋，李君，冯兵，王文青，柴旭超. 2021年青海玛多M_S7.4地震序列精定位与震源机制研究. 地震工程学报. 2022(01): 218-226 . 百度学术
4. 阮庆丰，刘俊清，田有，刘财，张宇，蔡宏雷. 松原地区地震发震机制与迁移特点研究. 地球物理学报. 2022(09): 3309-3321 . 百度学术
5. 李梦莹，张志宏，焦明若，方禹心. 2018-05-28松原M_S5.7地震地电场变化特征研究. 大地测量与地球动力学. 2022(12): 1276-1280+1299 . 百度学术
6. 李梦莹，杨士超，张志宏，孔祥瑞，孙庆山. 2018年松原M_S 5.7地震地电场变化特征. 地震地磁观测与研究. 2022(S1): 120-122 . 百度学术
7. Yu TANG，Aihua WENG，Yue YANG，Shiwen LI，Jianjun NIU，Yanhui ZHANG，Yabin LI，Jianping LI. Connection between earthquakes and deep fluids revealed by magnetotelluric imaging in Songyuan, China. Science China(Earth Sciences). 2021(01): 161-176 . 必应学术
8. 唐裕，翁爱华，杨悦，李世文，牛建军，张艳辉，李亚彬，李建平. 松原地震与流体作用联系的大地电磁证据. 中国科学:地球科学. 2021(01): 134-149 . 百度学术
9. 张志宏，郭安宁，李梦莹，黄明威，杨牧萍. 2018年吉林松原M_S5.7地震地磁异常分析. 科学技术与工程. 2021(29): 12406-12414 . 百度学术
10. 张洪艳，刘轶男，张帆，卢燕红，张宇. 吉林松原宁江地震序列精定位. 防灾科技学院学报. 2020(01): 32-37 . 百度学术
11. 李永生，赵谊，李继业，高峰，石伟. 2018年5月28日吉林松原M_S5.7地震发震构造分析. 地震学报. 2020(01): 12-23+120 . 本站查看
12. 阮庆丰，刘财，刘俊清，张宇，郑国栋. 2019年5月18日松原M5.1地震构造机制分析. 吉林大学学报(地球科学版). 2020(06): 1897-1904 . 百度学术
13. 李迎春，李锋，王俊菲，许田，王金艳. 2019年黄海M_L4.6地震序列的震源机制和发震构造. 地震学报. 2020(05): 543-551+507 . 本站查看
14. 王婷，延军平，李双双，万佳，张玉凤. 帕米尔高原Mw≥6.6级地震时间韵律特征. 高原地震. 2020(04): 6-16 . 百度学术
15. 李艳娥，邢成起，陈丽娟，郭祥云，康建红. 2017—2018年吉林松原地震序列研究. 地震学报. 2019(04): 435-444 . 本站查看
其他类型引用(5)
计量
- 文章访问数: 834
- HTML全文浏览量: 337
- PDF下载量: 43
- 被引次数: 20
