Spatio-temporal variation characteristic of the ultra-low frequency magnetic field prior to strong earthquakes of western Chinese mainland
-
摘要: 对2015年至2019年期间中国大陆西部的磁通门磁力仪秒采样观测资料开展了5—100 s频段的地磁垂直强度极化分析，并运用一些数学方法对分析结果进行了处理。结果显示，极化高值在经向和纬向均无明显的形态和幅值变化，且极化高值与地磁外源场扰动无关。在此基础上筛选出18次极化高值异常事件，利用插值方法得到了极化高值异常的空间分布图，并分析了18次高值异常与中国大陆西部及周边15次强震的时空关系。分析结果表明：地磁场出现多台同步极化高值现象后的半年内，高值区可能发生M6.0以上强震；高值现象出现后，多个高值区均有可能发生强震；后续强震的震级与高值区面积呈正相关。Abstract: Based on the one second sampling data from fluxgate magnetometers in western Chinese mainland from 2015 to 2019, this paper carried out the vertical component polarization analysis of the frequency band between 5 s to 100 s, and then processed the analysis results by some mathematical methods. The results show that the high polarization value has no obvious shape and amplitude change in both meridional and zonal directions, and the high polarization value has nothing to do with geomagnetic field disturbance. On this basis, 18 high-value anomalies were screened out, and their spatial distribution map was obtained by interpolation method. The results show that within half a year after the synchronous appearance of multiple high-value anomalies of the vertical component polarization of geomagnetic field, the high-value region may have a strong earthquake with magnitude over M6.0. After the appearance of high value phenomenon, several high value regions are likely to have strong earthquakes; furthermore the magnitude of subsequent strong earthquake seems to be positively correlated with the area of the high-value zone.
-
Keywords:
- geomagnetic field /
- polarization /
- strong earthquake /
- western Chinese mainland
-
引言
由于城市交通的立体化发展，地下轨道交通和城市立交桥交叉修建。地下结构和邻近地面桥梁相互作用，在地震过程中，高烈度地区的地下和地面结构之间存在复杂的动力相互作用，这种地下和地面结构的耦联破坏作用及其破坏机理已成为城市岩土地震工程及防震减灾领域重要的研究课题之一。
针对地下结构与地面结构的地震相互作用已有相关研究。杨书燕等（2007）对矩形隧道下穿地表框架结构和紧邻地表框架结构两个二维有限元模型进行了数值模拟，结果表明：地表建筑紧邻地下结构时地震响应最大；傅玉勇等（2009）对并行圆形隧道与地表框架结构动力相互作用体系进行了二维有限元数值模拟表明：输入水平向地震动时，地下隧道会使邻近地表建筑顶层水平位移放大30%，但输入竖向地震动时影响较小；何伟和陈健云（2012）通过采用二维有限元数值模拟方法研究了不同场地条件下矩形地下结构对地表建筑地震响应的影响，认为地下结构会使得上部区域的地震动减小而使远处区域增大，并且场地土体刚度越小地下结构对地表建筑地震响应的影响越大；郭靖等将地表建筑简化为弯剪梁模型并采用层间位移谱分析方法计算地表建筑的地震响应，分别采用二维有限元数值模型计算矩形地下结构和圆形隧道存在时的地面运动，并以此作为地表建筑弯剪梁模型底部的地震激励分析了矩形地下结构对地表建筑地震响应的影响，其结果表明，地下结构对地表建筑地震响应的影响与地表建筑的自振周期相关（郭靖，2013；郭靖，陈健云，2013；郭靖等，2017）；Wang等（2013）通过对ANSYS二次开发使其能在频域内求解，利用该方法对地下结构与地表建筑的三维有限元模型进行了求解，分析出矩形地下结构对邻近低矮结构地震响应的影响更加明显。
上述研究采用的数值模型大多未考虑或仅考虑了部分非线性特征，难以给出实际的结构地震破坏过程，因此尚无法合理给出地下结构与地面结构耦联灾变机理。近些年来，由于数值计算技术的发展，杜修力等（2017，2018）首次利用三维动力有限元数值模型模拟了大开车站塌毁过程。董瑞（2020）采用非线性数值模型精细地研究了矩形地下结构的地震灾变过程，得到了矩形地下结构的渐进式破坏模式，并指出结构形式（顶板和侧墙的线刚度比）对其地震响应有较大的影响。Xu等（2020）模拟了地下框架结构的破坏过程，认为中柱由于受到高轴压作用而发生完全脆性破坏。
精细的非线性有限元数值模型可以实现对结构地震破坏过程的高度仿真，但针对地下结构与地面结构的耦联灾变过程的模拟鲜有研究。为了研究地铁车站与横跨桥梁耦联破坏机理，本文拟采用非线性动力时程数值模拟方法，对软土场地的地铁车站与上部交叉桥梁模型进行数值模拟分析。模拟出一种潜在的地铁车站−桥梁耦联灾变模式；并对比弹性工况和弹塑性工况下单一车站模型和车站−桥梁模型的模拟结果，以期分析出地铁车站−桥梁耦联灾变机理，为地下结构与地面结构的选址及抗震设计提供科学依据。
1. 计算模型及模拟工况
1.1 计算模型
本文以两层三跨形式地铁车站及其上部横跨的高架桥结构-土-结构相互作用体系为研究对象，其横断面尺寸如图1所示。地铁车站的横断面（宽×高）为21.00 m×12.5 m；侧墙厚度为0.70 m （纵筋配筋率为0.8%），顶板厚度为0.70 m （纵筋配筋率为1.0%），底板厚度为0.8 m （纵筋配筋率为1.0%），中柱直径为0.8 m，中心距为5.0 m，（纵筋配筋率为6.0%，箍筋为直径9 mm；钢筋间距为350 mm）。上部横跨桥梁的桥面宽度为16 m；桥柱高度为8 m，截面尺寸为0.8 m×0.6 m；基础形式为桩基础，每个桥柱下面一根直径为0.6 m的混凝土圆形桩，桩长为11.5 m，承台厚度为0.8 m。
![]()
图 1 两层三跨地铁车站（a）及高架桥（b）示意图Figure 1. Diagram of two-story and three-span subway station （a） and viaduct （b）根据广州某工程场地勘探资料选取场地土层，土层参数列于表1。土层从上至下分别为素填土层、粉质黏土层、粉细砂层、粉质黏性土层、全风化混合花岗岩层和中风化混合花岗岩层，厚度分别为4.4，3.5，3.4，8.7，8.1和5.9 m。地铁车站与桥梁相互作用模型剖面如图2所示。车站埋深为5.1 m，邻近车站的桩基础距离车站侧边净距为1.2 m。
表 1 场地土层参数Table 1. Soil parameters of site土层序号 土层类别 密度/（103kg·m−3） 剪切波速/（m·s−1） 弹性模量/kPa 泊松比 黏聚力/kPa 摩擦角/° 厚度/m 1 素填土 1.80 24.39 3000 0.4 10 10 4.4 2 粉质黏土 1.87 27.59 3700 0.3 18 14 3.5 3 粉细砂 1.85 32.12 5000 0.31 0 28 3.4 4 粉质黏性土 1.91 34.37 6000 0.33 25 21 8.7 5 全风化混合花岗岩 1.93 38.37 7500 0.32 30 23 8.1 6 中风化混合花岗岩 2.50 742.42 3500000 0.27 300 35 5.9 在Abaqus中建立两层三跨车站与桥梁交叉的有限元模型，如图3所示。土体和混凝土均采用8结点6面体缩减积分单元进行离散；钢筋采用2结点杆单元进行离散；采用Abaqus中的嵌固模型模拟钢筋和混凝土间的协同工作。根据廖振鹏（2002）给出的离散网格中的波传播条件，选择土体单元尺寸（长×宽×深）为2 m×2 m×2 m，并在地下结构附近进行适当加密处理。
根据董瑞（2020）研究结果，选取模拟人工边界条件为：底部采用黏性边界，侧向边界采用自由度绑定边界（tied degrees of freedom，缩写为TDOF）（Doltsinis，1989），即分别在x和y向两组侧边界处设置位移约束条件使对应节点的位移协调；为保证计算精度，车站-桥梁模型土体计算范围为147 m×30 m×34 m。在Abaqus/Explicit模块中，黏性边界采用无限元（CIN3D8），绑定边界采用多点约束（multi-point constraints，缩写为MPC）方式实现。混凝土结构和周围土体间的摩擦滑移通过Abaqus中的接触面模型模拟，法向采用“硬接触”，切向采用摩擦接触（
$\; \mu {\text{＝}}0.4 $ ）。1.2 本构模型
土体本构模型采用基于Montáns （2000）模型提出的改进边界面本构模型（Dong et al，2020），这一模型已经被二次开发并可以用于通用有限元软件Abaqus。该模型在归一化偏平面内建立屈服函数；并采用土体材料模型与围压的经验公式描述围压改变引起土体模量变化；根据偏应力不变量建立应力等效关系，将一维滞回本构（Pyke）模型拓展到三维应力空间并推导塑性模量。归一化偏平面如图4所示，边界面（即Mises强度面）在归一化偏平面中的投影为一个圆，其方程可以表示为：
![]()
图 4 归一化偏平面内边界面及投影法则Figure 4. Bounding surface and mapping law in normalized deviatoric plane$$ F{\text{＝}}\frac{3}{2}{\boldsymbol{m}}_{{\rm{p}}{\text{，}}ij}{\boldsymbol{m}}_{{\rm{p}}{\text{，}}ij}{\text{－}}{m}_{\mathrm{c}}^{2}{\text{＝}}0 {\text{，}} $$ (1) 式中，
$ {\boldsymbol{m}}_{{\rm{p}}{\text{，}}ij} $ 表示投影点在一般应力空间内的i，j两个方向上的应力张量（i，j＝1，2，3表示空间内x，y，z方向），${m}_{\mathrm{c}}{\text{＝}}6\mathrm{s}\mathrm{i}\mathrm{n}\phi /（3{\text{－}}\mathrm{s}\mathrm{i}\mathrm{n}\phi ）$ 为边界面的半径。图4中的mi为与mij对应的主应力空间内的归一化应力张量，其中${\boldsymbol{m}}_{i}{\text{＝}}{\boldsymbol{s}}_{i}/\bar{{p}}$ ，$\bar{{p}}{\text{＝}}{p}{\text{＋}}c/\mathrm{t}\mathrm{a}\mathrm{n}\phi$ 为等效围压，p为围压，$ c $ 为黏聚力，$ \phi $ 为内摩擦角。基于第二偏应力不变量的应力状态等效关系如图5所示。根据应力状态参量
$ \eta $ 可将一维滞回本构模型拓展到三维应力空间。增量形式的一维滞回本构模型可以表示为![]()
图 5 应力等效关系。图中q为偏应力，$\bar {q}$ 为等效偏应力Figure 5. Equivalent relation of stress。 q is deviatoric stress， and$\bar {q}$ is equivalent deviatoric stress$$ {\rm{d}}\tau {\text{＝}}G{\Bigg(1{\text{－}}\frac{{\tau }_{{\rm{r}}}{\text{－}}\tau }{{\tau }_{{\rm{p}}}{\text{＋}}{\tau }_{{\rm{r}}}}\Bigg)}^{2}{\rm{d}}\gamma {\text{，}} $$ (2) 式中：
$ \tau $ 表示剪应力，$ {\tau }_{{\rm{p}}} $ 和$ {\tau }_{{\rm{r}}} $ 分别为投影点和反向点的剪应力；$ \mathrm{d}\gamma $ 为工程剪应变增量；初始剪切模量$ G $ 采用经验公式计算$$ G{\text{＝}}k{\mathrm{P}}_{\mathrm{a}}{\left(\frac{p}{{\mathrm{P}}_{\mathrm{a}}}\right)^{n}} {\text{，}} $$ (3) 式中，
$ G $ 表示最大剪切模量，$ {\mathrm{P}}_{{\rm{a}}} $ 表示一个标准大气压强，k和$ n $ 为与土性相关的材料常数。根据塑性模量和偏应力增量可以表示为（Dong et al，2020）：
$$ {K}_{\mathrm{p}}{\text{＝}}12G{（{m}_{{\rm{c}}}{\bar{{p}}}_{{\rm{p}}}）}^{2}\frac{{（1{\text{－}}\eta ）}^{2}}{1{\text{－}}{（1{\text{－}}\eta ）}^{2}} {\text{，}} $$ (4) $$ \mathrm{d}{\boldsymbol{s}}_{{\rm{c}}，ij}{\text{＝}}2G\mathrm{d}{\boldsymbol{e}}_{ij}{\text{－}}\frac{3G\left[1{\text{－}}{（1{\text{－}}\eta ）}^{2}\right]}{{m}_{{\rm{c}}}^{2}}（{\boldsymbol{m}}_{{\rm{p}}，{{kl}}}{\rm{d}}{\boldsymbol{e}}_{kl}）{\boldsymbol{m}}_{{\rm{p}}，{{ij}}} {\text{.}} $$ (5) 式中：e表示应变张量；k和l取值为1，2，3，表示空间内x，y，z方向。
地铁车站混凝土材料采用Abaqus材料库中的混凝土塑形损伤（concrete damaged plasticity，缩写为CDP）模型进行模拟，钢筋采用理想弹塑性材料模拟；混凝土强度等级为C30，钢筋强度等级为HRB335。CDP模型的计算参数列于表2，混凝土材料屈服强度及损伤因子随非弹性应变的变化关系如图6所示；假定桥梁已经按照规范进行了完善的抗震设计，因此可以忽略桥梁桩基础、桥柱和梁板的强度破坏，模拟时桥梁采用弹性模型并施加阻尼比为5%的质量相关阻尼模拟材料耗能特性，桥梁材料密度为2 500 kg/m3，弹性模量为30 GPa，泊松比为0.2；桥梁支座和梁板间设置接触模拟支座，法向采用“硬接触”，切向采用摩擦接触，桥梁支座和梁板间的摩擦系数取为0.4。
表 2 CDP模型材料参数Table 2. Parameter of CDP model弹性模量/GPa 泊松比 剪胀角/° 偏心率 双轴与单轴抗压强度比值 屈服面形态的参数K 黏性系数 损伤发展 30 0.2 30 0.1 1.16 0.666 7 0.000 5 见图6 ![]()
图 6 混凝土拉伸（a）和压缩（b）损伤发展Figure 6. Development of tensile damage （a） and compression damage （b） for concrete1.3 模拟工况
选择1995年阪神地震的Kobe波的南北分量和竖直分量（图7）作为输入地震波地震；分别选取有限元模型中x向（地铁车站横断面方向）和z向（竖直方向）作为水平地震动和竖直地震动输入方向，将水平方向分量和竖直方向分量调幅为0.4g和0.3g，以此模拟在强震作用下地铁车站的地震响应。
![]()
图 7 阪神地震中神户海洋气象台记录Figure 7. Ground motion record at the Kobe marine observatory during the Hanshin earthquake针对单一车站模型和车站−桥梁模型，分别进行车站为弹性模型和弹塑性模型情况下的动力时程反应分析，并依据上述两组模型模拟结果对比研究桥梁和地铁车站间的相互作用机理及耦联破坏机制。弹塑性工况P-1模拟了桥梁和地铁车站耦联破坏的震害现象，并与单一车站（P-2）的破坏模式进行对比给出桥梁和地铁车站间的耦联破坏机制。弹性工况E-1和E-2用于分析在未发生破坏的条件下，邻近车站的桥梁基础与车站的相互作用，进一步解释桥梁−车站耦联破坏机制。数值模拟工况如表3所示。
表 3 数值模拟工况Table 3. simulated Working condition工况 计算模型 地铁车站
本构模型输入地震动 工况 计算模型 地铁车站
本构模型输入地震动 水平/g 竖向/g 水平/g 竖向/g P-1 车站−桥梁 CDP 0.4 0.3 E-1 车站−桥梁 黏弹性 0.4 0.3 P-2 单一车站 CDP 0.4 0.3 E-2 单一车站 黏弹性 0.4 0.3 2. 有限元模拟结果分析
2.1 地铁车站−桥梁耦联灾变模拟结果
车站−桥梁模型（P-1工况）4个关键时刻的水平位移与车站混凝土压损伤如图8所示，可以看出：当t＝0—7 s 时，地铁车站中柱混凝土在双向地震动作用下损伤逐渐发育，但地铁车站和桥梁均未发生破坏；当t＝7—8 s 时，地铁车站中柱混凝土损伤进一步发育，下层左侧柱子首先发生屈曲破坏；当t＝8—11 s 时，地铁车站下层左侧中柱失效后，顶板和楼板失去中柱约束发生破坏，顶板和楼板失效后，侧墙无法承受两侧土压力的作用而向车站内侧坍塌；当t＝11—15 s时，地铁车站破坏后，邻近桥梁基础随土层向车站方向产生较大变形，进而产生了“落梁”破坏。
![]()
图 8 车站−桥梁模型4个关键时刻的水平位移及车站混凝土压损伤因子云图Figure 8. Nephogram of horizontal displacement and compression damage of station concrete for station-bridge model at four key moments2.2 地铁车站−桥梁耦联破坏机理
地震荷载作用下地下结构受到周围土体的约束而发生剪切变形，地下结构的层间位移可以反应地下结构受到的地震荷载作用大小。图9给出弹性工况（E-1和E-2）和弹塑性工况（P-1和P-2）在结构破坏前（3—7 s）地铁车站上下两层在左侧、中部和右侧的层间相对位移。单一车站模型和车站−桥梁模型计算得到的车站层间位移基本一致。由此可知，在地铁车站未破坏时，桥梁对车站的地震响应影响很小。
![]()
图 9 地铁车站层间相对位移（a） 弹性工况；（b） 弹塑性工况（3—7 s）Figure 9. Relative displacement of subway station（a） Elastic condition；（b） Elastic-plastic condition （3−7 s）车站−桥梁模型（P-1工况）和单一车站模型（P-2工况）地铁车站破坏后引起邻近区域土层永久变形（变形放大比例为1），如图9所示。图中显示了车站楼板（D＝−11.3 m）位置水平位移等值线，单一车站模型左右两侧水平位移等值线与水平方向夹角分别为38.26°和42.94°，车站−桥梁模型左右两侧水平位移等值线与水平方向夹角分别为42.90°和43.00°。
两模型在地表（D＝0 m）、车站顶板（D＝−5.1 m）和车站楼板（D＝−11.3 m）位置的土层永久变形分布曲线，如图11所示。车站−桥梁模型得到的车站邻近区域土层永久位移大于单一车站模型，桥梁会加剧车站破坏引起的地基变形以及其影响范围。结合图10和图11的结果可知，桥梁不会改变地铁车站的破坏形式；但是由于桥梁的基础对土层存在一定的加固效果使得车站两侧土体整体向车站移动，进而加剧了车站−桥梁系统的震害。
![]()
图 10 车站邻近区域土层永久水平位移云图（a） 车站-桥梁模型 ；（b） 单一车站模型Figure 10. Nephogram of soil permanent horizontal displacement near subway station（a） Station-bridge model；（b） Single station model![]()
图 11 地表、车站顶板和车站楼板位置处的土层永久位移Figure 11. Permanent horizontal displacement of soil at the depths of soil surface，station roof and station floor3. 讨论与结论
本文通过对软土场地的单一两层三跨车站以及其与邻近桥梁耦联相互作用模型进行非线性动力时程反应分析，模拟了一种可能的地铁车站-桥梁模型耦联破坏模式。并通过对比分析桥梁对弹性地铁车站模型（E-1工况和E-2工况）和弹塑性地铁车站模型（P-1工况和P-2工况）地震响应的影响，对地铁车站与桥梁的耦联破坏机制进行了分析，得出如下结论：
1） 依托Abaqus强大的大变形模拟计算能力，并采用本文给出的土体非线性本构模型，可以实现对地下结构与地面结构复杂耦联破坏过程的仿真模拟，地铁车站发生破坏后会引起邻近区域土层发生较大的水平和竖向变形，地铁车站的破坏会引起邻近区域地基失效并使桥梁发生“落梁”破坏。
2） 分别根据弹性模型和弹塑性模型（未破坏阶段）模拟结果，对比单一地铁车站和地铁车站-桥梁模型得到车站响应，可知桥梁对车站的地震响应影响较小，车站发生破坏主要取决于地震作用以及其本身的强度。
3） 地铁车站-桥梁耦联破坏主要表现为地铁车站破坏后桥梁基础失效形式的“次生”灾害， 在车站破坏后，桥梁基础会带动临近区域土层整体向车站方向移动，增大周围土层永久变形的范围，加剧车站坍塌震害。
需要说明的是，对地下结构地震破坏过程的模拟是一个高度非线性的问题，此类问题的非线性主要表现为几何非线性和材料非线性两个方面，选取合适的数值模型描述上述两种非线性是仿真模拟破坏过程的关键。针对几何非线性问题，本文采用的Abaqus能够实现对单元网格大变形的模拟；针对材料非线性问题，本文采用基于边界面理论建立的土体弹塑性本构模型，在该模型中引入了模量随围压变化的经验公式，可以更合理的描述地下结构两侧土体围压改变引起的土体非线性响应。
综上所述，本研究利用了一种可以模拟地下结构和地面结构耦联灾变过程的方法，并通过对地铁车站-桥梁相互作用模型地震灾变过程的仿真模拟结果，得出了地铁车站和邻近桥梁的耦联灾变机理，研究成果对于城市地下结构选址及抗震设计具有一定的参考意义。
-
![]()
图 1 2011年至2019年国家地磁台网中心磁通门磁力仪数量
Figure 1. The number of fluxgate magnetometers in National Geomagnetic Networks Center of China from 2011 to 2019
![]()
图 2 2019年中国西部磁通门磁力仪（秒采样）分布
Figure 2. Distribution of fluxgate magnetometers （second sampling） in western Chinese mainland in 2019
![]()
图 3 喀什台2017年YZH1与YZH2的对比
Figure 3. Comparison of YZH1 with YZH2 results of Kashi station in 2017
![]()
图 4 沿经线和纬线分布的部分台站及其YZH和YZH1 2017年结果
Figure 4. Location of some stations and their YZH and YZH1 results in 2017
![]()
图 5 喀什台2017年YZH1序列与Kp （a）和Dst （b）指数曲线的对比
Figure 5. Comparison of YZH1 series with Kp （a） and Dst （b） curves in Kashgar station in 2017
![]()
图 6 YZH1
高值异常在12个月中的频次分布 Figure 6. The frequency distribution of YZH1 high value anomalies in 12 months
![]()
图 7 YZH1高值异常频次与地震频次图
Figure 7. The frequency of YZH1 high value anomalies and earthquake frequency
![]()
![]()
图 9 2015年至2016年磁场高值日YZH2的空间分布图
Figure 9. Spatial distribution of YZH2 on the magnetic field high value day from 2015 to 2016
![]()
图 10 2017年至2019年磁场
高值日YZH2的空间分布图 Figure 10. Spatial distribution of YZH2 on the magnetic field high value day from 2017 to 2019
![]()
图 11 喀什台、静海台和高邮台的空间分布（a）以及2017年极化YZH1时序曲线对比图（b）
图中红色竖线标示中国大陆西部YZH1出现高值异常的时间
Figure 11. The locations of the stations Kashi，Jinghai and Gaoyou （a） and YZH1 time series curves in 2017 （b）
The red vertical line indicates the time when the YZH1 high value anomaly appeared in western Chinese mainland
![]()
图 12 2015—2019年6组高值异常叠加前、后的YZH2与后续6个月内地震的空间关系
Figure 12. Spatial relationship of six groups of high YZH2 anomalies before and after superposition as well as the distribution of earthquakes in the following six months from 2015 to 2019
![]()
图 12 2015—2019年6组高值异常叠加前、后的YZH2与后续6个月内地震的空间关系
Figure 12. Spatial relationship of six groups of high YZH2 anomalies before and after superposition as well as the distribution of earthquakes in the following six months from 2015 to 2019
![]()
图 13 叠加后的YZH2高值区面积与后续6个月内地震震级的关系
Figure 13. The relationship between YZH2 high-value areas of the superposition of several high-value anomalies and the magnitude of strong earthquakes in six months
表 1 2015年至2020年1月中国大陆西部及周边M≥6.0强震信息
Table 1 Information of strong earthquakes （M≥6.0） during the period from 2015 to January of 2020 in western Chinese mainland and its surrounding areas
发震时刻
年-月-日 时:分:秒M 北纬/° 东经/° 震源深度/km 参考位置 2020-01-19 21:27:55 6.4 39.83 77.21 16 新疆伽师县 2019-11-21 07:50:45 6.0 19.50 101.35 10 老挝 2019-06-17 22:55:43 6.0 28.34 104.9 16 四川长宁县 2019-04-24 04:15:48 6.3 28.40 94.61 10 西藏墨脱县 2017-11-18 06:34:19 6.9 29.75 95.02 10 西藏米林县 2017-08-09 07:27:52 6.6 44.27 82.89 11 新疆精河县 2017-08-08 21:19:46 7.0 33.2 103.82 20 四川九寨沟县 2016-12-08 13:15:03 6.2 43.83 86.35 6 新疆呼图壁县 2016-11-25 22:24:30 6.7 39.27 74.04 10 新疆阿克陶县 2016-10-17 15:14:49 6.2 32.81 94.93 9 青海杂多县 2016-06-26 19:17:11 6.7 39.43 73.4 10 吉尔吉斯斯坦 2016-01-21 01:13:13 6.4 37.68 101.62 10 青海门源县 2015-12-07 15:50:03 7.4 38.2 72.9 30 塔吉克斯坦 2015-07-03 09:07:46 6.5 37.6 78.2 10 新疆皮山县 2015-04-25 14:11:26 8.1 28.2 84.7 20 尼泊尔 
下载: 导出CSV
表 2 2015年至2019年中国大陆西部
$Y_{ZH2}$ 高值异常信息Table 2 Information of high YZH2 value anomalies in western Chinese mainland from 2015 to 2019
序号 异常出现时间
年-月-日异常持续时间/d 异常台站最多日（异常日）
年-月-日异常台站数 总台数 YZH2 1 2015-01-16 6 2015-01-18 28 36 0.73 2 2015-02-09 9 2015-02-13 14 36 0.39 3 2015-07-17 5 2015-07-18 8 36 0.29 4 2015-09-26 5 2015-09-29 13 27 0.47 5 2015-11-18 11 2015-11-23 21 27 0.92 6 2016-01-27 8 2016-01-30 17 41 0.75 7 2016-08-14 9 2016-08-17 15 33 0.92 8 2016-09-11 9 2016-09-15 14 33 0.88 9 2016-12-01 6 2016-12-02 14 37 0.75 10 2017-01-13 5 2017-01-15 8 38 0.41 11 2017-02-13 5 2017-02-14 10 37 0.49 12 2017-03-11 9 2017-03-18 14 34 0.66 13 2017-08-24 6 2017-08-28 22 40 0.65 14 2017-10-29 6 2017-10-30 14 40 0.70 15 2018-04-16 3 2018-04-17 11 43 0.62 16 2018-06-10 6 2018-06-12 17 43 0.44 17 2019-09-22 3 2019-09-24 13 54 1.15 18 2019-11-12 5 2019-11-12 26 56 0.85
下载: 导出CSV
表 3 单次高值异常及短期内多次高值异常叠加后发生强震的高值区面积
Table 3 The high value areas of single high-value anomalies and the superposition of several high-value anomalies in a short period which followed with strong earthquakes
序号 出现高值日期
年-月-日震后出现高值区的
地震（提前的天数）半年内高值区强震 高值区面积①
/（104 km2）高值区总面积②
/（104 km2）1 2015-01-18 无 2015-07-03皮山M6.5 251.6 313 2 2015-02-13 2015-04-25尼泊尔M8.1 3 2015-07-18 2015-07-03皮山M6.5
（15）2015-12-07塔吉克斯坦M7.4 13 20.6 4 2015-09-29 无 2016-01-21门源M6.4 20 144 5 2015-11-23 2015-12-07塔吉克斯坦M7.4 123.7 6 2016-01-30 2016-01-21门源M6.4 （9） 2016-06-26吉尔吉斯斯坦M6.7 17.5 84.7 7 2016-08-17 无 2016-12-08呼图壁M6.2 24.1 82.2 8 2016-09-15 2016-10-17杂多M6.2 58 9 2016-12-02 无 2016-12-08呼图壁M6.2 18.2 80.8 10 2017-01-13 2016-12-08呼图壁M6.2 （37） 无 47.2 47.2 11 2017-02-14 无 2017-08-08九寨沟M7.0 41.6 84.4 12 2017-03-18 2017-08-09精河M6.6 30.5 13 2017-08-28 2017-08-08九寨沟M7.0 （19）
2017-08-09精河M6.6 （20）无 39.5 14 2017-10-31 无 16.1 15 2018-04-17 无 无 32.4 16 2018-06-12 17 2019-09-24 无 2019-11-21老挝M6.0 19.2 72.3 18 2019-11-12 注：① 指地震震中所在高值区的面积；② 指一次高值现象出现的多个高值区面积之和。
下载: 导出CSV
-
杜爱民,周志坚,徐文耀,杨少峰. 2004. 新疆和田ML7.1地震前ULF电磁辐射的激发机理[J]. 地球物理学报,47(5):832–837. doi: 10.3321/j.issn:0001-5733.2004.05.014 Du A M,Zhou Z J,Xu W Y,Yang S F. 2004. Generation mechanisms of ULF electromagnetic emissions before the ML=7.1 earthquake at Hotan of Xinjiang[J]. Chinese Journal of Geophysics,47(5):832–837 (in Chinese). doi: 10.3321/j.issn:0001-5733.2004.05.014（inChinese）
冯志生,李琪,卢军,李鸿宇,居海华,孙海军,杨福喜,张翼. 2010. 基于磁通门秒值数据的地震ULF磁场可靠信息提取研究[J]. 华南地震,30(2):1–7. doi: 10.3969/j.issn.1001-8662.2010.02.001 Feng Z S,Li Q,Lu J,Li H Y,Ju H H,Sun H J,Yang F X,Zhang Y. 2010. The seismic ULF geomagnetic reliable information exaction based on fluxgate magnetometer data of second value[J]. South China Journal of Seismology,30(2):1–7 (in Chinese).
何畅,冯志生. 2017. 极化方法在成都地磁台的应用研究[J]. 地震学报,39(4):558–564. doi: 10.11939/jass.2017.04.010 He C,Feng Z S. 2017. Application of polarization method to geomagnetic data from the station Chengdu[J]. Acta Seismologica Sinica,39(4):558–564 (in Chinese).
李琪,杨星,蔡绍平. 2015. 极化方法应用于地磁台阵的震例分析[J]. 震灾防御技术,10(2):412–417. doi: 10.11899/zzfy20150222 Li Q,Yang X,Cai S P. 2015. Case study of applying polarization method to geomagnetic array data[J]. Technology for Earthquake Disaster Prevention,10(2):412–417 (in Chinese).
马宗晋. 1980. 华北地壳的多(应力集中)点场与地震[J]. 地震地质,2(1):39–47. Ma Z J. 1980. Multipoints (concentrated) stress field and earthquakes in North China crust[J]. Seismology and Geology,2(1):39–47 (in Chinese).
梅世蓉. 1996. 地震前兆场物理模式与前兆时空分布机制研究(二):强震孕育时应力、应变场的演化与地震活动、地震前兆的关系[J]. 地震学报,18(1):1–10. Mei S R. 1996. Study on physical model of seismic precursors and time-space distribution mechanism of precursors(2)[J]. Acta Seismologica Sinica,18(1):1–10 (in Chinese).
潘晖,张建国,杨冬梅,李美. 2014. 地震地磁低点位移成因的一种可能机制[J]. 大地测量与地球动力学,34(4):83–87. doi: 10.14075/j.jgg.2014.04.029 Pan H,Zhang J G,Yang D M,Li M. 2014. A possible mechanism of the contributing factor of geomagnetic lowest shift[J]. Journal of Geodesy and Geodynamics,34(4):83–87 (in Chinese).
汤吉,詹艳,王立凤,徐建郎,赵国泽,陈小斌,董泽义,肖骑彬,王继军,蔡军涛,徐光晶. 2008. 5月12日汶川8.0级地震强余震观测的电磁同震效应[J]. 地震地质,30(3):739–745. doi: 10.3969/j.issn.0253-4967.2008.03.012 Tang J,Zhan Y,Wang L F,Xu J L,Zhao G Z,Chen X B,Dong Z Y,Xiao Q B,Wang J J,Cai J T,Xu G J. 2008. Coseismic signal associated with aftershock of the MS8.0 Wenchuan earthquake[J]. Seismology and Geology,30(3):739–745 (in Chinese). doi: 10.3969/j.issn.0253-4967.2008.03.012（inChinese）
汤吉,詹艳,王立凤,董泽义,赵国泽,徐建郎. 2010. 汶川地震强余震的电磁同震效应[J]. 地球物理学报,53(3):526–534. doi: 10.3969/j.issn.0001-5733.2010.03.006 Tang J,Zhan Y,Wang L F,Dong Z Y,Zhao G Z,Xu J L. 2010. Electromagnetic coseismic effect associated with aftershock of Wenchuan MS8.0 earthquake[J]. Chinese Journal of Geophysics,53(3):526–534 (in Chinese). doi: 10.3969/j.issn.0001-5733.2010.03.006（inChinese）
姚休义,滕云田. 2017. 2010年元谋MS5.2地震前地磁极化异常研究[J]. 地震研究,40(3):444–448. doi: 10.3969/j.issn.1000-0666.2017.03.021 Yao X Y,Teng Y T. 2017. Variation of geomagnetic polarization value associated with the 2010 Yuanmou MS5.2 earthquake[J]. Journal of Seismological Research,40(3):444–448 (in Chinese).
姚休义,冯志生. 2018. 地震磁扰动分析方法研究进展[J]. 地球物理学进展,33(2):511–520. doi: 10.6038/pg2018BB0063 Yao X Y,Feng Z S. 2018. Review on the recent development of analysis methods on magnetic disturbance associated with earthquakes[J]. Progress in Geophysics,33(2):511–520 (in Chinese).
袁桂平,张学民,吴迎燕,赵旭东. 2015. 汶川8.0级地震前地磁低点位移与内外源Sq等效电流体系关系的研究[J]. 地震,35(3):102–112. doi: 10.3969/j.issn.1000-3274.2015.03.011 Yuan G P,Zhang X M,Wu Y Y,Zhao X D. 2015. Minimum point shift of the geomagnetic vertical component in diurnal variation and the internal-external equivalent current system Sq before the 2008 Wenchuan MS8.0 earthquake[J]. Earthquake,35(3):102–112 (in Chinese).
Chang X T,Zou B,Guo J Y,Zhu G B,Li W,Li W D. 2017. One sliding PCA method to detect ionospheric anomalies before strong earthquakes:Cases study of Qinghai,Honshu,Hotan and Nepal earthquakes[J]. Adv Space Res,59(8):2058–2070. doi: 10.1016/j.asr.2017.02.007
Cianchini G,De Santis A,Barraclough D R,Wu L X,Qin K. 2012. Magnetic transfer function entropy and the 2009 MW=6.3 L’Aquila earthquake (Central Italy)[J]. Nonlinear Proc Geophys,19(4):401–409. doi: 10.5194/npg-19-401-2012
Currie J L,Waters C L. 2014. On the use of geomagnetic indices and ULF waves for earthquake precursor signatures[J]. J Geophys Res,119(2):992–1003. doi: 10.1002/2013JA019530
Gladychev V,Baransky L,Schekotov A,Fedorov E,Pokhotelov O,Andreevsky S,Rozhnoi A,Khabazin Y,Belyaev G,Gorbatikov A,Gordeev E,Chebrov V,Sinitsin V,Lutikov A,Yunga S,Kosarev G,Surkov V,Molchanov O,Hayakawa M,Uyeda S,Nagao T,Hattori K,Noda Y. 2001. Study of electromagnetic emissions associated with seismic activity in Kamchatka region[J]. Nat Hazards Earth Syst Sci,1(3):127–136. doi: 10.5194/nhess-1-127-2001
Han P,Hattori K,Xu G J,Ashida R,Chen C H,Febriani F,Yamaguchi H. 2015. Further investigations of geomagnetic diurnal variations associated with the 2011 off the Pacific coast of Tohoku earthquake (MW9.0)[J]. J Asian Earth Sci,114:321–326. doi: 10.1016/j.jseaes.2015.02.022
Hattori K. 2004. ULF geomagnetic changes associated with large earthquakes[J]. Terr Atmos Ocean Sci,15(3):329–360. doi: 10.3319/TAO.2004.15.3.329(EP)
Hattori K,Han P,Yoshino C,Febriani F,Yamaguchi H,Chen C H. 2013. Investigation of ULF seismo-magnetic phenomena in Kanto,Japan during 2000−2010:Case studies and statistical studies[J]. Surv Geophys,34(3):293–316. doi: 10.1007/s10712-012-9215-x
Hayakawa M,Kawate R,Molchanov O A,Yumoto K. 1996. Results of ultra-low-frequency magnetic field measurements during the Guam earthquake of 8 August 1993[J]. Geophys Res Lett,23(3):241–244. doi: 10.1029/95GL02863
Hayakawa M,Itoh T,Hattori K,Yumoto K. 2000. ULF electromagnetic precursors for an earthquake at Biak,Indonesia on February 17,1996[J]. Geophys Res Lett,27(10):1531–1534. doi: 10.1029/1999GL005432
Huang Q H,Ikeya M. 1998. Seismic electromagnetic signals (SEMS) explained by a simulation experiment using electromagnetic waves[J]. Phys Earth Planet Inter,109(3/4):107–114. doi: 10.1016/S0031-9201（98）00135-6
Ismaguilov V S,Kopytenko Y A,Hattori K,Voronov P M,Molchanov O A,Hayakawa M. 2001. ULF magnetic emissions connected with under sea bottom earthquakes[J]. Nat Hazards Earth Syst Sci,1(1/2):23–31. doi: 10.5194/nhess-1-23-2001
Karakelian D,Klemperer S L,Fraser-Smith A C,Thompson G A. 2002. Ultra-low frequency electromagnetic measurements associated with the 1998 MW5.1 San Juan Bautista,California earthquake and implications for mechanisms of electromagnetic earthquake precursors[J]. Tectonophysics,359(1/2):65–79.
Li Q,Zhu P Y,Mamatemin A,Xu X G. 2011. Detection of ULF electromagnetic emissions as a precursor to two earthquakes in China[J]. Earthquake Science,24(6):601–607. doi: 10.1007/s11589-011-0822-2
Masci F. 2011. On the seismogenic increase of the ratio of the ULF geomagnetic field components[J]. Phys Earth Planet Inter,187(1/2):19–32. doi: 10.1016/j.pepi.2011.05.001
Matsushima M,Honkura Y,Oshiman N,Bariş Ş,Tunçer M K,Tank S B,Çelik C,Takahashi F,Nakanishi M,Yoshimura R,Pektaş R,Komut T,Tolak E,Ito A,Iio Y,Işikara A M. 2002. Seismoelectromagnetic effect associated with the Izmit earthquake and its aftershocks[J]. Bull Seismol Soc Am,92(1):350–360. doi: 10.1785/0120000807
Molchanov O A,Kopytenko Y A,Voronov P M,Kopytenko E A,Matiashvili T G,Fraser-Smith A C,Bernardi A. 1992. Results of ULF magnetic field measurements near the epicenters of the Spitak (MS=6.9) and Loma Prieta (MS=7.1) earthquakes:Comparative analysis[J]. Geophys Res Lett,19(14):1495–1498. doi: 10.1029/92GL01152
Molchanov O A,Schekotov A,Fedorov E,Belyaev G,Gordeev E. 2003. Presismic ULF electromagnetic effect from observation at Kamchatka[J]. Nat Hazards Earth Syst Sci,3(3/4):203–209. doi: 10.5194/nhess-3-203-2003
Nagao T,Enomoto Y,Fujinawa Y,Hata M,Hayakawa M,Huang Q,Izutsu J,Kushida Y,Maeda K,Oike K,Uyeda S,Yoshino T. 2002. Electromagnetic anomalies associated with 1995 Kobe earthquake[J]. J Geodynam,33(4/5):401–411.
Prattes G,Schwingenschuh K,Eichelberger H U,Magnes W,Boudjada M,Stachel M,Vellante M,Wesztergom V,Nenovski P. 2008. Multi-point ground-based ULF magnetic field observations in Europe during seismic active periods in 2004 and 2005[J]. Nat Hazard Earth Syst Sci,8(3):501–507. doi: 10.5194/nhess-8-501-2008
Thomas J N,Love J J,Johnston M J S,Yumoto K. 2009. On the reported magnetic precursor of the 1993 Guam earthquake[J]. Geophys Res Lett,36(16):L16301. doi: 10.1029/2009GL039020
Uyeda S,Hayakawa M,Nagao T,Molchanov O,Hattori K,Orihara Y,Gotoh K,Akinaga Y,Tanaka H. 2002. Electric and magnetic phenomena observed before the volcano seismic activity in 2000 in the Izu island region,Japan[J]. Proc Natl Acad Sci USA,99(11):7352–7355. doi: 10.1073/pnas.072208499
Wang Z D,Chen B,Yuan J H,Yang F X,Jia L,Wang C. 2018. Localized geomagnetic field anomalies in an underground gas storage[J]. Phys Earth Planet Inter,283:92–97. doi: 10.1016/j.pepi.2018.08.009
Wessel P,Smith W H F. 1995. New version of the generic mapping tools[J]. Eos Trans AGU,76(33):329.
计量
- 文章访问数: 821
- HTML全文浏览量: 474
- PDF下载量: 62
