Geophysical analysis on the tectonic difference between northern and southern segments of Xianshuihe fault zone
-
摘要: 对鲜水河断裂带重磁异常进行向上延拓，通过计算观测面高度异常与延拓后异常之间的相关系数得出最佳向上延拓高度，该延拓高度所对应的延拓结果即为研究区构造背景产生的异常值。向上延拓结果显示布格重力异常值沿鲜水河断裂带自北西向南东逐渐增大，反映出下地壳底边界沿该方向呈升高趋势，可能由青藏高原地壳软弱物质 “东向逃逸” 所致。化极后ΔT磁异常延拓结果表明鲜水河断裂带南东段的康定—石棉以东为强磁性刚性基底。以鲜水河断裂带为分界的不同地块之间基底岩石及地层物性的不同是断裂带南北磁性差异的主导原因。对鲜水河断裂带两侧各约50 km范围内地震的震源深度进行统计，经投影至剖面及线性拟合求出鲜水河断裂带的三维几何形状。结果表明鲜水河断裂带总体倾向南西，倾角近乎直立，范围约为57°—88°。Abstract: Xianshuihe fault zone is one of the strongly active seismic belts in the eastern margin of Tibetan Plateau, with the features of large scale, high activity and seismic frequency. Four M≥7.0 earthquakes have occurred along Xianshuihe fault zone. It is indicated that the differences of geometry features, activity characteristics and intensity, surface rupture characteristics of earthquakes exist between northern and southern segments of Xianshuihe fault zone. The upward continuation result of gravity-magnetic anomaly in Xianshuihe fault zone is present. The optimum upward continuation height is determined by calculating correlation coefficient of anomaly at observed surface height and upward continuation height. The upward continuation result at this height is corresponding to the anomaly caused by structural setting. It shows that the bottom of lower crust is on the rise from NW to SE according to the consecutive increase of Bouguer gravity anomaly at such orientation along the fault, which probably results from eastward motion of crustal soft material in Tibetan Plateau area. Magnetic anomaly reduced to the pole reveals that there is strongly magnetic rigid base to east of Kangding-Shimian in the southern segment of Xianshuihe fault zone. The difference in physical property of basement rocks and stratum in different blocks bounded by Xianshuihe fault zone is the major reason for magnetic diversity between northern and southern segments of the fault zone. With the statistical results of focal depths of seismic events in the range of 50 kilometers on both sides of the fault zone, its three-dimensional geometry is determined by projection and linear fitting. It is indicated that the dip direction of Xianshuihe fault zone is SW and its dip angle is almost vertical in the range of 57° to 88°.
-
引言
我国是世界上岩溶最为发育的国家之一，岩溶面积达363万平方千米，主要分布在黔、桂、滇东、川南、湘西等省区，其中分布最广的是贵州省，约占全省面积的73%，其次是广西壮族自治区，约占全区面积的34% （覃子建等，1996）。近年来随着我国经济的飞速发展，工农业、居民生活用水迅速增加，大量地下水被抽取，导致地下水位不断下降，塌陷事件也逐年增多，岩溶塌陷已成为我国六大地质灾害之一（赵博超等，2015）。然而，天然地震与岩溶塌陷的波形具有较高的相似性，分析人员依靠经验根据波形特征进行二者的识别费时较多，缺乏时效性。因此亟需研究一种能够及时准确地识别天然地震与岩溶塌陷事件的方法，建立准确的事件目录，为快速开展塌陷灾害防治工作提供重要的科技支撑，为地方经济发展、国土资源规划提供理论基础。
近年来国内外研究人员针对地震事件类型的识别开展了广泛而深入的研究，总结了许多时域和频域方面的识别判据。时域的判据主要包括初动方向、震源深度、体波震级与面波震级比、P波初动振幅与P波和S波最大振幅比、P波与S波最大振幅比、P波和S波最大振幅与尾波持续时间比、体波和面波震相特征、P波视速度差异等（吴宝峰，闫峰，2009；Pasyanos，2010；Taylor，2011；刘莎等，2012；王婷婷，边银菊，2015）；而频域的判据主要有能量频谱分布、拐角频率、倒谱曲线、瞬态谱等（张萍等，2005；Yılmaz et al，2013 ；靳玉贞等，2015；Kortström et al，2016 ；叶庆东等，2017）。这些方法各有特色，均取得了很好的效果，对于建立准确的地震事件目录具有重要意义；但是这些方法多是单独使用时域特征或频域特征来进行识别，而且地质条件、仪器参数、波形传播路径等差异使得其普适性不足。因此，许多研究人员使用在时域和频域均具有很高分辨率的小波包变换方法来进行特征判据提取。和雪松等（2006）通过对地震和矿震波形进行小波包变换，再计算二者的奇异值，从而达到识别地震的目的；王婷婷等（2014）使用小波包对地震和爆破波形进行分解，根据P波和S波的发育频段不同，使用不同分解系数中的P波和S波能量进行对比，获得了地震和爆破的识别判据；Beccar-Varela等（2016）通过小波变换对美国亚利桑那州M3.0—3.5地震和爆破事件进行分解，利用分解得出的时频能量特征进行地震事件的识别。尽管上述利用小波包分解提取判据的方法能够有效地进行地震识别，但是这些方法主要是针对地震与爆破或矿震，至于地震与岩溶塌陷事件识别的研究则较少。此外，上述方法尽管提取了地震波形在时频域的特征，但均未涉及地震波形分形维数的特征。
为此，本文拟根据天然地震和岩溶塌陷在不同小波包分解频带内存在的差异，结合分形维数计算方法，对比各小波包分解系数的维数差异，以期提出一种能够体现出事件波形时频域特征和分形维数特征并且能够快速有效地识别天然地震与岩溶塌陷的新方法。
1. 资料
本文从广西地震台网中心2008—2015年记录到的塌陷和地震事件中，选取记录台站较多、波形较为清晰的广西境内的39个事件，包括17次塌陷事件和22次浅源地震事件，震源深度为5—10 km，震级范围为ML1.0—2.9，所有事件的类型均已由广西地震台网中心确认。从这39个事件中，提取出400条记录清晰的事件波形记录，其中天然地震事件和岩溶塌陷事件的波形各200条，最大、最小震中距分别为237 km和25 km。所选地震、塌陷事件和记录台站的分布如图1所示。
本文所选用事件的波形记录的采样率均为100 sps，波形记录的频率范围为0—50 Hz；所有事件只使用垂直向波形进行计算，事件波形的截取宽度均为120 s （即12 000个点），S波的初动位置尽量位于截取窗口的中部。
![]()
图 1 事件和地震台站位置示意图Figure 1. Location of earthquakes （green circles），karst collapses （red circles） and seismic stations （triangles） in this study2. 识别算法设计及基本原理
2.1 识别算法设计
一般来说，一套完整的自动识别算法主要由两个模块组成：其一为特征提取模块，用于从目标数据中提取出识别所需的特征参数；其二为识别判定模块，利用已提取出的特征参数对目标类型进行判定。本文在设计两个模块时考虑到：① 本文所使用的特征提取算法应能体现出天然地震和岩溶塌陷的震源机制在时频域的特点，同时还要能描述出两种事件在不同频带下波形的复杂性；② 由于测震仪器、地质构造、传播路径等差异，每一组事件所提取的特征数据可能存在一定差异，因此识别算法需要具有一定的鲁棒性和容错性，同时从现实应用的角度出发，算法还需具有较快的运算速度。
鉴于此，本文使用小波包分形方法来提取天然地震和岩溶塌陷的特征，再通过径向基函数(radial basis function，简写为RBF)神经网络进行事件分类的识别判定。小波包分形方法利用了小波包在时频域下高分辨的特点，并通过分形理论对不同频带下波形信号的复杂度进行描述，从而达到兼顾地震信号的细节分析与整体描述合理性的目的（曹茂森等，2005；谢全民等，2011；朱权洁等，2012）。而RBF神经网络除了神经网络普遍具有的自适应性强、能够实现各种非线性映射、有较高的鲁棒性和容错性、运算速度快等特点外，还具备结构简单、收敛快等优势。所以本文使用小波包分形作为特征提取模块，RBF神经网络作为识别判定模块。首先，对事件波形进行归一化处理，然后利用小波包分形提取出事件波形特征；其次，利用已知事件组成训练集完成对RBF神经网络的训练，并将训练好的神经网络用于对未知事件的类型进行识别。图2为本文算法的流程框图。
2.2 小波包变换
小波变换是一种重要的时频域波形分析方法，在数字信号处理、图像处理、声纹识别等方面均有广泛应用。而小波包变换是对小波变换的进一步扩展，弥补了小波变换对高频细节特征提取的不足。小波包变换方法是通过多组低通、高通正交滤波器将信号的各频段不断细分，以此来提高信号在整个频带内的分辨能力，进而获得信号中所包含的更多细节。也正因为此，小波包变换也被称作信号时频域的显微镜。
2.3 分形维数
分形维数是分形理论中用于描述曲线或者图形不规则程度的重要特征量，但是，在实际使用中往往难以直接计算得到不规则分形体的维数，而只能计算其近似值。因此，本文从快速计算的角度出发，选用计盒维数来计算维数。
计盒维数DB（F）定义为（张济忠，1995）：对于需要计算维数的目标图形F，构造足够多的边长均为δ的正方形盒子将其完全覆盖（图3），所用到的盒子数量为Nδ（F），那么F的计盒维数DB（F）则为当δ趋于0时Nδ（F）增加的对数速率，可以表示为
${D_{\rm{B}}}(F) {\text{＝}} \mathop {\lim }\limits_{\delta \to 0} \frac{{\lg {N_\delta }(F)}}{{ - \lg \delta }}.$
(1) 3. 特征提取
由于地震与岩溶塌陷的震源性质和传播路径有所不同，台站记录到两种事件的波形在时频谱中呈现不同的特征。因此，通过时频域分析来进行地震和塌陷事件的特征提取具有可行性。图4为归一化处理后的地震和塌陷事件的时频能量分布图，所示的两事件是由贞丰地震台记录到的地震事件（2010-01-25，ML2.9，震中距Δ=128 km）和塌陷事件（2015-09-09，ML1.8，Δ=124 km），两个事件的震中均位于河池市天峨县。可以看出，地震的优势频率处于0—20 Hz范围内，而塌陷的优势频率处于0—10 Hz范围内。
![]()
图 4 地震（a）、岩溶塌陷（b）事件的波形（上）及其时频谱（下）Figure 4. Waveforms （upper） and spectrograms （lower） of the earthquake （a） and the karst collapse （b）3.1 事件波形的小波包分解
使用小波包进行变换前首先需要选取小波包的基函数。常用的基函数有Haar，Daubechies，Mexican Hat，Morlet，Meyer，Symlet，Coiflet，双正交小波等，不同的基函数性质不同。考虑到基函数需具备快速计算的能力，以便增强算法的时效性，同时需要具备精确重构的性质，以便能够完整地反映原始波形的频率特性，为此，本文选取具有紧支撑正交性的Daubechies作为基函数。考虑到需要避免因支撑区间太大而产生边界问题或者因区间太小而使得能量无法集中，选择区间长度为7的4阶Daubechies小波，简写为db4。
对本文所选的400条事件波形记录进行db4小波包的4层分解，得到16个分解系数波形，对应的频带宽度分别为0—3.125 Hz，3.125—6.250 Hz，…，46.875—50 Hz。图5和图6分别给出了河池市天峨县的地震和塌陷事件的小波包分解结果。可以看出，虽然地震与塌陷事件的优势频率不同，但是在各频段内均有震动波形的记录。因此，为了获取事件在不同频段内的波形特征，提取特征时需对各频段的数据进行计算。
![]()
图 5 2010年1月25日ML2.9地震事件的第四层小波包分解Figure 5. The fourth-layer wavelet packet decomposition of the ML2.9 earthquake on 25 January 2010![]()
图 6 2015年9月9日ML1.8岩溶塌陷事件的第四层小波包分解Figure 6. The fourth-layer wavelet packet decomposition of the ML1.8 karst collapse on 9 September 20153.2 小波包系数的分形维数计算
本文采用计盒维数法计算分形维数（张济忠，1995），分别计算各个事件的小波包分解系数，最终获得由16个维数值组成的一维向量，此即为本文所要提取的地震和岩溶塌陷事件的特征。下面将以这16个维数值为输入向量，进行RBF人工神经网络训练和事件类型识别。图7给出了地震和岩溶塌陷事件小波包分解系数的分形维数计算结果，可以看出，岩溶坍塌事件在不同频段内波形的维数高于地震事件波形的维数。
![]()
图 7 地震与岩溶坍塌事件的小波包分形维数计算结果对比Figure 7. Comparison of wavelet packet fractal dimension of earthquake with that of collapse4. RBF神经网络
RBF神经网络是一种经典的前馈3层神经网络，根据人脑皮层神经网络结构提出，具备以任意精度逼近任意一个连续函数的能力，由输入层、隐含层、输出层组成，其中输入层将输入数据与网络连接起来，隐含层对输入层传递过来的数据进行非线性变换，输出层则负责对隐含层进行加权后输出（Haykin，2001）。图8给出了单输出的RBF神经网络结构图，由图可知整个网络的输出yn可以表示为
${y_n} \text{＝} \sum\limits_{i\text{＝}1}^m {{W_{ni}}{g_i}(x)} \text{，}\quad i \text{＝} 1, \; 2, \; 3, \; \cdots , \; m\text{，}$
(2) 式中，n为样本输入数据的数量，m为隐含层节点数，gi（x）为神经网络的径向基函数。
本文所使用的RBF神经网络结构为：① 输入层：神经网络的输入为小波包分形得出的16个特征值，所以输入层节点数为16，节点的权值均为1；② 隐含层：径向基函数选择应用广泛的高斯函数，其表达式为
${g_i}(x) \text{＝} \exp \left[ { - \frac{{||x - {C_i}|{|^2}}}{{{\sigma ^2}}}} \right]\text{，}\quad i \text{＝} 1, \; 2, \; 3, \; \cdots , \; m\text{，}$
(3) This page contains the following errors:
error on line 1 at column 1: Start tag expected, '<' not foundBelow is a rendering of the page up to the first error.
RBF神经网络是否优秀，关键在于隐含层的训练，即基函数中心的选取、节点数量的选择以及节点权值的训练。本文使用matlab神经网络工具箱来进行隐含层的训练，目标均方差设为0.01，高斯函数扩展步长为2.0。训练的具体步骤为：① 建立一个0神经元的网络；② 寻找最大方差样本；③ 增加一个神经元，并使用最大方差的输入样本作为基函数中心；④ 重新设计隐含层权值，使用能够实现最小输出误差的值作为隐含层权值；⑤ 检查输出的均方差是否小于所设定的目标值，若满足该条件则网络训练完毕，否则返回第二步，从剩余样本中继续寻找最大方差样本。
5. 识别能力检验
神经网络的训练和测试由训练集和测试集完成。本文按照发震时间的先后顺序，将前10次地震事件的100条波形数据和前9个岩溶塌陷事件的100条波形数据组成训练集，用于神经网络训练；测试集则是由余下12次地震事件的100条波形数据和8次岩溶塌陷事件的100条波形数据组成，用于验证分类能力。之后分别计算测试集和训练集的地震和岩溶塌陷事件的小波包分形维数，计算结果列于表1。
表 1 地震和岩溶塌陷事件的小波包分形维数均值和方差Table 1. Mean and variance of fractal dimension of earthquake and collapse events小波包
分解系数训练集（地震） 训练集（塌陷） 测试集（地震） 测试集（塌陷） 维数均值 维数方差 维数均值 维数方差 维数均值 维数方差 维数均值 维数方差 1 1.379 6 0.010 0 1.463 0 0.009 2 1.389 0 0.010 4 1.440 8 0.008 5 2 1.341 8 0.008 2 1.412 2 0.010 0 1.337 6 0.007 7 1.438 0 0.006 5 3 1.354 4 0.013 3 1.435 1 0.014 0 1.358 0 0.010 1 1.466 3 0.010 2 4 1.344 9 0.009 7 1.400 4 0.013 2 1.333 9 0.008 1 1.444 3 0.009 1 5 1.416 3 0.017 9 1.483 0 0.018 5 1.426 5 0.017 2 1.505 5 0.011 0 6 1.388 8 0.016 0 1.438 4 0.016 0 1.395 4 0.012 8 1.478 4 0.010 2 7 1.368 3 0.014 0 1.450 2 0.019 9 1.364 6 0.012 0 1.471 0 0.008 8 8 1.370 0 0.014 5 1.451 4 0.016 6 1.370 1 0.012 6 1.462 8 0.010 3 9 1.411 5 0.011 8 1.539 9 0.015 9 1.414 9 0.013 4 1.557 0 0.009 4 10 1.438 5 0.017 3 1.508 4 0.010 2 1.443 9 0.015 7 1.524 9 0.011 6 11 1.423 4 0.017 0 1.493 3 0.010 0 1.433 3 0.015 5 1.518 4 0.007 7 12 1.429 8 0.017 4 1.486 3 0.010 5 1.440 5 0.016 2 1.513 1 0.009 6 13 1.422 9 0.019 2 1.530 1 0.009 8 1.438 5 0.017 2 1.539 2 0.007 7 14 1.422 4 0.018 8 1.520 0 0.009 5 1.433 2 0.017 8 1.536 9 0.006 0 15 1.421 1 0.018 4 1.497 6 0.010 0 1.425 9 0.016 6 1.516 1 0.008 6 16 1.420 3 0.017 6 1.501 0 0.016 2 1.429 7 0.017 4 1.524 8 0.007 1 从表1可以看出，岩溶塌陷事件各小波包分解系数计算出的分形维数值均明显大于地震事件对应分解系数的维数值。使用训练集中200条波形小波包分形特征向量进行RBF神经网络的训练，训练完成后，分别使用训练集和测试集的波形进行分类检验。训练集和测试集的识别结果如表2所示，可以看出：对于训练集，本文的识别方法可以达到100%的识别率；而对于测试集，地震和岩溶塌陷的识别率可分别达到87%和92%，能够有效地识别出天然地震与岩溶塌陷事件。由表2还可以看出，本文错误识别的共有9个地震事件13条波形记录和5个岩溶塌陷事件的8条波形，其相关信息列于表3。通过分析错误识别的事件波形，本文认为导致错误识别的原因在于：地脉动等背景噪声幅值较大，导致事件波形信号不明显；事件波形存在零漂现象；事件波形中叠加了少量干扰，出现了多个波形信号；事件波形的毛刺多，信号不够平滑；事件波形中有断记，部分波形为直线。
表 2 小波包分形维数判据的识别结果Table 2. The recognition result of the criteria based on wavelet packet fractal事件类型 测试样本数 正确识别样本数 错误识别样本数 识别率 地震（训练集） 100 100 0 100% 塌陷（训练集） 100 100 0 100% 地震（测试集） 100 87 13 87% 塌陷（测试集） 100 92 8 92% 表 3 错误识别事件的相关信息Table 3. The directory of false recognition序号 发震时间
年-月-日事件类型 ML 地点 东经/° 北纬/° 台站 震中距/km 1 2010-08-26 地震 2.2 钦州 22.24 109.47 LNS 26 2 2010-09-29 地震 2.7 河池 25.24 106.99 TE 36 3 2010-09-29 地震 2.7 河池 25.24 106.99 YIZ 188 4 2011-02-28 地震 2.9 百色 24.42 108.11 YAX 101 5 2011-02-28 地震 2.9 百色 24.42 108.11 PGX 122 6 2011-02-28 地震 2.9 百色 24.42 108.11 GAL 136 7 2011-10-03 地震 2.9 百色 24.61 105.40 YTT 143 8 2012-03-28 地震 2.5 来宾 23.64 108.90 PGX 136 9 2013-03-28 地震 2.2 崇左 21.88 107.12 PGX 128 10 2013-06-17 地震 2.3 崇左 23.42 107.50 SHX 115 11 2014-02-06 地震 2.5 百色 23.38 107.51 SHX 118 12 2014-02-06 地震 2.5 百色 23.38 107.51 QIP 122 13 2015-01-09 地震 2.3 钦州 22.25 108.04 LNS 112 14 2010-06-01 塌陷 1.9 防城港 23.71 109.23 LIB 34 15 2015-01-14 塌陷 1.0 百色 23.14 106.39 JXX 25 16 2015-09-08 塌陷 2.4 来宾 23.88 108.89 XCT 35 17 2015-09-08 塌陷 2.1 来宾 23.85 108.89 DHX 82 18 2015-09-08 塌陷 2.1 来宾 23.85 108.89 TAP 95 19 2015-09-08 塌陷 2.1 来宾 23.85 108.89 JL 200 20 2015-09-08 塌陷 2.1 来宾 23.88 108.89 TAP 95 21 2015-09-09 塌陷 1.6 河池 25.24 107.00 GAL 90 天然地震事件和岩溶塌陷事件的波形差异的本质是两种事件产生的震源机制不同。天然地震从本质上来讲地壳介质在构造应力作用下的破裂过程，通常震源深度较深，破裂过程较长，同时地震波经地下介质多次反射、折射、散射、叠加，导致地震仪所记录到的地震波频率范围更宽，尾波持续时间更久；而岩溶塌陷则发生在地表或者地下几十米处，是由于塌陷体从周边岩体上脱落并撞击底部岩体而产生的瞬时压缩源，能量释放时间短，波的传播路径较简单，未经过多层介质干扰，频带范围窄，尾波持续时间短（靳玉贞等，2015）。
为了能够得到量化的地震波形小波包分解的特征，本文使用分形维数来计算特征参数。地震波形分形维数的大小能够在一定程度上反映出岩石所受加载应力的差异。当岩石受到不同大小应力作用时，所产生的瞬态弹性波（称为声发射技术）的分形维数变化规律也验证了这一点。尹贤刚等（2009）运用分形理论方法研究了不同加载应力比阶段岩石声发射的分形维数，结果表明加载应力达到岩石峰值应力的40%左右时，分形维数即可达到最大值，之后逐步减小直至发生破裂；张黎明等（2015）对大理岩卸围压破坏过程的研究结果表明，应力达到峰值前，声发射分形维数会大幅度降低。由此可见，由于岩溶塌陷事件震动所产生的加载力并未达到岩石破裂峰值，而天然地震事件的加载力已经达到破裂峰值，所以岩溶塌陷事件波形的分形维数要大于地震事件的分形维数，这可以作为识别两种事件的判据。
6. 讨论与结论
本文利用基于小波包的分形维数和RBF神经网络识别技术，对2008—2015年广西境内的17次岩溶塌陷事件和22次浅源地震事件进行了识别研究和识别效果的检验，得到以下几点认识：
1） 小波包分形维数识别方法对于训练集的识别能力能够达到100%，对于测试集的平均识别能力达89.5%，说明该方法能够有效地识别天然地震与岩溶塌陷事件。而且，一个地震事件通常会被多个台站记录到，综合多个台站的识别结果，可提高事件识别的准确率。
2） 时频谱图和小波包分解系数对比显示，天然地震波形的优势频段为0—20 Hz，岩溶塌陷的优势频段为0—10 Hz。尽管两种事件的优势频率不同，但是二者所产生的波形在各个频段上均有记录，因此进行时频分析时不能只计算优势频率，应将全频段考虑在内。
3） 地震事件的各小波包分形维数均小于塌陷事件的小波包分形维数，这可能与岩石受到的应力大小有关，岩石受力越接近临界状态，所产生震动的分形维数越小。
在实际工作中，可将本文获得的量化指标作为快速自动化识别地震与岩溶塌陷事件的判定依据之一。鉴于地震破裂过程、传播路径等因素的复杂性，地震事件性质的判定还需结合近台记录P波初动符号、P波与S波振幅比、现场地震地质调查、地震宏观现象等资料综合判定。
不同地区的地质背景、塌陷方式、传播路径等方面差异较大，波形记录的特征也有所不同。因此，本研究获得的指标仅具有一定参考价值，不一定适用于其它地区。本文研究样本为已知类型事件，下一步将收集更多的事件样本，进行更广泛的识别训练和研究，不断地完善和验证判定指标的普适性。
-
![]()
图 1 鲜水河断裂带及邻区地质构造简图
（a） 构造背景（阴影部分为本文研究区）（据徐锡伟等，2005）；（b） 区域地质
Figure 1. Generalized geological structure map of Xianshuihe fault zone and its adjacent areas
（a） Tectonic setting map （The dash area is the studied area of this article）（after Xu et al，2005）；（b） Regional geological map
![]()
图 2 相关系数、布格重力异常与化极后ΔT磁异常向上延拓结果
（a） 布格重力异常观测值与向上延拓值相关系数及其斜率随延拓高度的变化曲线；（b） 化极后ΔT磁异常观测值与向上延拓值相关系数及其斜率随延拓高度的变化曲线；（c） 布格重力异常观测值；（d） 化极后ΔT磁异常；（e） 向上延拓40 km后的布格重力异常；（f） 向上延拓20 km后化极后的ΔT磁异常
Figure 2. Upward continuation results of correlation coefficients，Bouguer gravity anomaly and polarized aeromagnetic anomaly
（a） Variation curve of correlation coefficients and its gradient calculated by observed Bouguer gravity anomaly and upward continuation anomaly；（b） Variation curve of correlation coefficients and its gradient calculated by observed polarized aeromagnetic anomaly and upward continuation anomaly；（c） Observed Bouguer gravity anomaly map；（d） Polarized aeromagnetic anomaly map；（e） Upward continuation result of Bouguer gravity anomaly；（f） Upward continuation result of polarized aeromagnetic anomaly
![]()
图 3 地震震源深度典型剖面反演及断裂三维几何形状拟合等值面图
（a） BB′剖面震源深度分布及断层线性拟合；（b） EE′剖面震源深度分布及断层线性拟合；（c） 鲜水河断裂带两侧地震深度分布图；（d） 鲜水河断裂带三维几何形状拟合等值面图
Figure 3. The inversion of focal depths of seismic events in typical profiles and the map of three-dimensional geometry fitting isosurface of Xianshuihe fault zone
（a） Distribution of focal depths of seismic events and linear fitting of the fault in profile BB′；（b） Distribution of focal depths of seismic events and linear fitting of the fault in profile EE′；（c） Map of focal depths distribution on both sides of Xianshuihe fault zone；（d） Map of three-dimensional geometry fitting isosurface of Xianshuihe fault zone
![]()
图 4 重磁异常延拓与鲜水河断裂带三维几何形状切片图
红色实心圆代表参与断裂带三维几何形状拟合的MS≥1.0地震事件。（a） 布格重力异常不同高度上延后水平总导数切片图（蓝色为低值，红色为高值，下同）；（b） 化极后ΔT磁异常不同高度上延后水平总导数切片图
Figure 4. The upward continuation result of gravity-magnetic anomaly and the slices of three-dimensional geometry about Xianshuihe fault zone
The red solid circles mean the MS≥1.0 earthquakes that contribute to the fitting. （a） The slices of total horizontal derivative results after upward continuation of Bouguer gravity anomaly at different depths （blue color means low value and red color means high value，the same below）；（b） The slices of total horizontal derivatives results after upward continuation of magnetic anomaly reduced to the pole at different depths
-
白登海,滕吉文,马晓冰,孔祥儒. 2011. 大地电磁观测揭示青藏高原东部存在两条地壳物质流[J]. 中国基础科学,13(1):7–10. doi: 10.3969/j.issn.1009-2412.2011.01.002 Bai D H,Teng J W,Ma X B,Kong X R. 2011. Crustal flow beneath eastern Tibetan Plateau revealed by magnetotelluric measurements[J]. China Basic Science,13(1):7–10 (in Chinese).
陈杰,陈宇坤,丁国瑜,田勤俭,王赞军,单新建,任金卫,赵瑞斌,王志才. 2003. 2001年昆仑山口西8.1级地震地表破裂带[J]. 第四纪研究,23(6):629–639. doi: 10.3321/j.issn:1001-7410.2003.06.006 Chen J,Chen Y K,Ding G Y,Tian Q J,Wang Z J,Shan X J,Ren J W,Zhao R B,Wang Z C. 2003. Surface rupture zones of the 2001 earthquake MS8.1 west of Kunlun pass,northern Qinghai-Xizang plateau[J]. Quaternary Sciences,23(6):629–639 (in Chinese).
陈立春,王虎,冉勇康,孙鑫喆,苏桂武,王继,谭锡斌,李智敏,张晓清. 2010. 玉树MS7.1级地震地表破裂与历史大地震[J]. 科学通报,55(13):1200–1205. Chen L C,Wang H,Ran Y K,Sun X Z,Su G W,Wang J,Tan X B,Li Z M,Zhang X Q. 2010. The MS7.1 Yushu earthquake surface rupture and large historical earthquakes on the Garzê-Yushu fault[J]. Chinese Science Bulletin,55(31):3504–3509. doi: 10.1007/s11434-010-4079-2
陈运泰. 2008. 汶川特大地震的震级和断层长度[J]. 科技导报,26(10):26–27. doi: 10.3321/j.issn:1000-7857.2008.10.007 Chen Y T. 2008. On the magnitude and the fault length of the great Wenchuan earthquake[J]. Science and Technology Review,26(10):26–27 (in Chinese).
邓起东,程绍平,马冀,杜鹏. 2014. 青藏高原地震活动特征及当前地震活动形势[J]. 地球物理学报,57(7):2025–2042. doi: 10.6038/cjg20140701 Deng Q D,Cheng S P,Ma J,Du P. 2014. Seismic activities and earthquake potential in the Tibetan Plateau[J]. Chinese Journal of Geophysics,57(7):2025–2042 (in Chinese).
高玲举,张健,董淼. 2015. 川西高原重磁异常特征与构造背景分析[J]. 地球物理学报,58(8):2996–3008. doi: 10.6038/cjg20150831 Gao L J,Zhang J,Dong M. 2015. The study of gravity-magnetic anomaly and tectonic background in Sichuan west region[J]. Chinese Journal of Geophysics,58(8):2996–3008 (in Chinese).
国家地震科学数据共享中心. 2016. 中国地震台网统一地震目录[EB/OL]. [2018−03−26]. http://data.earthquake.cn. China Earthquake Data Center. 2016. Unified Seismic Catalogue of China Seismic Network[EB/OL]. [2018−03−26]. http://data.earthquake.cn （in Chinese）. 何梅兴,方慧,王绪本,卢景奇,袁永真,白大为,杜炳锐,仇根根,高宝屯. 2017. 鲜水河断裂带南段深部电性结构特征研究[J]. 地球物理学报,60(6):2414–2424. doi: 10.6038/cjg20170630 He M X,Fang H,Wang X B,Lu J Q,Yuan Y Z,Bai D W,Du B R,Qiu G G,Gao B T. 2017. Deep conductivity characteris-tics of the southern Xianshuihe fault zone[J]. Chinese Journal of Geophysics,60(6):2414–2424 (in Chinese).
黄福明,杨智娴. 1987. 鲜水河断裂带的应力积累与释放[J]. 地震学报,9(2):18–32. Huang F M,Yang Z X. 1987. Stress accumulation and release in the Xianshuihe fault zone[J]. Acta Seismologica Sinica,9(2):18–32 (in Chinese).
黄金莉,赵大鹏,郑斯华. 2001. 川滇活动构造区地震层析成像[J]. 地球物理学报,44(增刊1):127–135. Huang J L,Zhao D P,Zheng S H. 2001. Seismic tomography of the Sichuan-Yunan active tectonic region[J]. Chinese Journal of Geophysics,44(S1):127–135 (in Chinese).
阚荣举,张四昌,晏凤桐,俞林胜. 1977. 我国西南地区现代构造应力场与现代构造活动特征的探讨[J]. 地球物理学报,20(2):96–109. Kan R J,Zhang S C,Yan F T,Yu L S. 1977. Present tectonic stress field and its relation to the characteristics of recent tectonic activity in southwestern China[J]. Acta Geophysics Sinica,20(2):96–109 (in Chinese).
李海兵,付小方,Van Der Woerd J,司家亮,王宗秀,侯立玮,邱祝礼,李宁,吴富峣,许志琴,Tapponnier P. 2008. 汶川地震(MS8.0)地表破裂及其同震右旋斜向逆冲作用[J]. 地质学报,82(12):1623–1643. doi: 10.3321/j.issn:0001-5717.2008.12.002 Li H B,Fu X F,van der Woerd J,Si J L,Wang Z X,Hou L W,Qiu Z L,Li N,Wu F Y,Xu Z Q,Tapponnier P. 2008. Co-seisimic surface rupture and dextral-slip oblique thrusting of the MS8.0 Wenchuan earthquake[J]. Acta Geologica Sinica,82(12):1623–1643 (in Chinese).
钱洪. 1988. 鲜水河断裂带上潜在震源区的地质学判定[J]. 四川地震,(2):22–30. Qian H. 1988. Identification of potential seismic area along the Xianshuihe zone[J]. Earthquake Research in Sichuan,(2):22–30 (in Chinese).
任金卫,王敏. 2005. GPS观测的2001年昆仑山口西MS8.1级地震地壳变形[J]. 第四纪研究,25(1):34–44. doi: 10.3321/j.issn:1001-7410.2005.01.006 Ren J W,Wang M. 2005. GPS measured crustal deformation of the MS8.1 Kunlun earthquake on November 14th 2001 in Qinghai-Xizang plateau[J]. Quaternary Sciences,25(1):34–44 (in Chinese).
单斌,李佳航,韩立波,房立华,杨嵩,金笔凯,郑勇,熊熊. 2012. 2010年MS7.1级玉树地震同震库仑应力变化以及对2011年MS5.2级囊谦地震的影响[J]. 地球物理学报,55(9):3028–3042. doi: 10.6038/j.issn.0001-5733.2012.09.021 Shan B,Li J H,Han L B,Fang L H,Yang S,Jin B K,Zheng Y,Xiong X. 2010. Coseismic coulomb stress change caused by 2010 MS=7.1 Yushu earthquake and its influence to 2011 MS=5.2 Nangqên earthquake[J]. Chinese Journal of Geophysics,55(9):3028–3042 (in Chinese).
四川省地质矿产局. 1982. 区域地质志[M]. 北京: 地质出版社: 306−309. Sichuan Bureau of Geology and Mineral Resources. 1982. Regional Geology[M]. Beijing: Geological Publishing House: 306−309 （in Chinese）.
宋鸿彪,罗志立. 1995. 四川盆地基底及深部地质结构研究的进展[J]. 地学前缘,2(3/4):231–237. Song H B,Luo Z L. 1995. The study of the basement and deep geological structures of Sichuan basin,China[J]. Earth Science Frontiers,2(3/4):231–237 (in Chinese).
孙洁,晋光文,白登海,王立凤. 2003. 青藏高原东缘地壳、上地幔电性结构探测及其构造意义[J]. 中国科学:D辑,33(增刊1):173–180. Sun J,Jin G W,Bai D H,Wang L F. 2003. Sounding of electrical structure of the crust and upper mantle along the eastern border of Qinghai-Tibet Plateau and its tectonic significance[J]. Science in China:Series D,46(S2):243–253.
唐文清,陈智粱,刘宇平,张清志,赵济湘,Burchfiel B C,King R W. 2005. 青藏高原东缘鲜水河断裂与龙门山断裂交会区现今的构造活动[J]. 地质通报,24(12):1169–1172. doi: 10.3969/j.issn.1671-2552.2005.12.013 Tang W Q,Chen Z L,Liu Y P,Zhang Q Z,Zhao J X,Burchfiel B C,King R W. 2005. Present-day tectonics activity in the intersection area of the Xianshuihe fault and Longmenshan fault on the eastern margin of the Qinghai-Tibet Plateau[J]. Geological Bulletin of China,24(12):1169–1172 (in Chinese).
王椿镛,王溪莉,苏伟,常利军,韩渭宾,吕智勇,戴仕贵. 2006. 青藏高原东缘下地壳流动的地震学证据[J]. 四川地震,(4):1–4. doi: 10.3969/j.issn.1001-8115.2006.04.001 Wang C Y,Wang X L,Su W,Chang L J,Han W B,Lü Z Y,Dai S G. 2006. Seismological evidence of the crust flowing under the eastern boundary part of Qinghai-Tibet Plateau[J]. Earthquake Research in Sichuan,(4):1–4 (in Chinese).
王辉,刘杰,申旭辉,刘勉,李青松,石耀霖,张国民. 2010. 断层分布及几何形态对川西及邻区应变分配的影响[J]. 中国科学: 地球科学,40(4):458–472. Wang H,Liu J,Shen X H,Liu M,Li Q S,Shi Y L,Zhang G M. 2010. Influence of fault geometry and fault interaction on strain partitioning within western Sichuan and its adjacent region[J]. Science China Earth Sciences,53(7):1056–1070. doi: 10.1007/s11430-010-3062-8
王卫民,郝金来,姚振兴. 2013. 2013年4月20日四川芦山地震震源破裂过程反演初步结果[J]. 地球物理学报,56(4):1412–1417. doi: 10.6038/cjg20130436 Wang W M,Hao J L,Yao Z X. 2013. Preliminary result for rupture process of Apr. 20,2013,Lushan earthquake,Sichuan,China[J]. Chinese Journal of Geophysics,56(4):1412–1417 (in Chinese).
闻学泽,Allen C R,罗灼礼,钱洪,周华伟,黄伟师. 1989. 鲜水河全新世断裂带的分段性,几何特征及其地震构造意义[J]. 地震学报,11(4):362–372. Wen X Z,Allen C R,Luo Z L,Qian H,Zhou H W,Huang W S. 1989. Segmentation,geometric features,and their seismotectonic implications for the Holocene Xianshuihe fault zone[J]. Acta Seismologica Sinica,11(4):362–372 (in Chinese).
闻学泽. 2000. 四川西部鲜水河-安宁河-则木河断裂带的地震破裂分段特征[J]. 地震地质,22(3):239–249. doi: 10.3969/j.issn.0253-4967.2000.03.005 Wen X Z. 2000. Character of rupture segmentation of the Xianshuihe-Anninghe-Zemuhe fault zone,western Sichuan[J]. Seismology and Geology,22(3):239–249 (in Chinese).
熊熊,滕吉文. 2002. 青藏高原东缘地壳运动与深部过程的研究[J]. 地球物理学报,45(4):507–515. doi: 10.3321/j.issn:0001-5733.2002.04.008 Xiong X,Teng J W. 2002. Study on crustal movement and deep process in eastern Qinghai-Xizang plateau[J]. Chinese Journal of Geophysics,45(4):507–515 (in Chinese).
徐锡伟,陈文彬,于贵华,马文涛,戴华光,张志坚,陈永明,何文贵,王赞军,党光明. 2002. 2001年11月14日昆仑山库赛湖地震(MS8.1)地表破裂带的基本特征[J]. 地震地质,24(1):1–13. doi: 10.3969/j.issn.0253-4967.2002.01.001 Xu X W,Chen W B,Yu G H,Ma W T,Dai H G,Zhang Z J,Chen Y M,He W G,Wang Z J,Dang G M. 2002. Characteris-tic features of the Hoh Sai Hu (Kunlunshan) earthquake (MS8.1),northern Tibetan Plateau,China[J]. Seismology and Geology,24(1):1–13 (in Chinese).
徐锡伟,闻学泽,郑荣章,马文涛,宋方敏,于贵华. 2003. 川滇地区活动块体最新构造变动样式及其动力来源[J]. 中国科学:D辑,33(增刊1):151–162. Xu X W,Wen X Z,Zheng R Z,Ma W T,Song F M,Yu G H. 2003. Pattern of latest tectonic motion and its dynamics for active blocks in Sichuan-Yunnan region,China[J]. Science in China:Series D,46(S2):210–226.
徐锡伟,张培震,闻学泽,秦尊丽,陈桂华,朱艾斓. 2005. 川西及其邻近地区活动构造基本特征与强震复发模型[J]. 地震地质,27(3):446–461. doi: 10.3969/j.issn.0253-4967.2005.03.010 Xu X W,Zhang P Z,Wen X Z,Qin Z L,Chen G H,Zhu A L. 2005. Features of active tectonics and recurrence behaviors of strong earthquakes in the western Sichuan Province and its adjacent regions[J]. Seismology and Geology,27(3):446–461 (in Chinese).
徐锡伟,陈桂华,王启欣,陈立春,任治坤,许冲,魏占玉,鲁人齐,谭锡斌,董绍鹏,石峰. 2017. 九寨沟地震发震断层属性及青藏高原东南缘现今应变状态讨论[J]. 地球物理学报,60(10):4018–4026. doi: 10.6038/cjg20171028 Xu X W,Chen G H,Wang Q X,Chen L C,Ren Z K,Xu C,Wei Z Y,Lu R Q,Tan X B,Dong S P,Shi F. 2017. Discussion on seismogenic structure of Jiuzhaigou earthquake and its implication for current strain state in the southeastern Qinghai-Tibet Plateau[J]. Chinese Journal of Geophysics,60(10):4018–4026 (in Chinese).
曾华霖. 2005. 重力场与重力勘探[M]. 北京: 地质出版社: 29. Zeng H L. 2005. Gravity Field and Gravity Exploration[M]. Beijing: Geological Publishing House: 29 （in Chinese）.
曾华霖,许德树. 2002. 最佳向上延拓高度的估计[J]. 地学前缘,9(2):499–504. doi: 10.3321/j.issn:1005-2321.2002.02.029 Zeng H L,Xu D S. 2002. Estimation of optimum upward continuation height[J]. Earth Science Frontiers,9(2):499–504 (in Chinese).
曾祥方,罗艳,韩立波,石耀霖. 2013. 2013年4月20日四川芦山MS7.0地震: 一个高角度逆冲地震[J]. 地球物理学报,56(4):1418–1424. doi: 10.6038/cjg20130437 Zeng X F,Luo Y,Han L B,Shi Y L. 2013. The Lushan MS7.0 earthquake on 20 April 2013: A high-angle thrust event[J]. Chinese Journal of Geophysics,56(4):1418–1424 (in Chinese).
张桂芳,屈春燕,单新建,刘云华,宋小刚. 2011. 2010年青海玉树MS7.1级地震地表破裂带和形变特征分析[J]. 地球物理学报,54(1):121–127. doi: 10.3969/j.issn.0001-5733.2011.01.013 Zhang G F,Qu C Y,Shan X J,Liu Y H,Song X G. 2011. The surface rupture and coseismic deformation characteristics of the MS7.1 earthquake at Qinghai Yushu in 2010[J]. Chinese Journal of Geophysics,54(1):121–127 (in Chinese).
张培震,徐锡伟,闻学泽,冉勇康. 2008. 2008年汶川8.0级地震发震断裂的滑动速率、复发周期和构造成因[J]. 地球物理学报,51(4):1066–1073. doi: 10.3321/j.issn:0001-5733.2008.04.015 Zhang P Z,Xu X W,Wen X Z,Ran Y K. 2008. Slip rates and recurrence intervals of the Longmenshan active fault zone,and tectonic implications for the mechanism of the May 12 Wenchuan earthquake,2008,Sichuan,China[J]. Chinese Journal of Geophysics,51(4):1066–1073 (in Chinese).
张旭,冯万鹏,许力生,李春来. 2017. 2017年九寨沟MS7.0级地震震源过程反演与烈度估计[J]. 地球物理学报,60(10):4105–4116. doi: 10.6038/cjg20171035 Zhang X,Feng W P,Xu L S,Li C L. 2017. The source-process inversion and the intensity estimation of the 2017 MS7.0 Jiu-zhaigou earthquake[J]. Chinese Journal of Geophysics,60(10):4105–4116 (in Chinese).
张勇,许力生,陈运泰. 2013. 芦山4.20地震破裂过程及其致灾特征初步分析[J]. 地球物理学报,56(4):1408–1411. doi: 10.6038/cjg20130435 Zhang Y,Xu L S,Chen Y T. 2013. Rupture process of the Lushan 4.20 earthquake and preliminary analysis on the disaster-causing mechanism[J]. Chinese Journal of Geophysics,56(4):1408–1411 (in Chinese).
周荣军,何玉林,黄祖智,黎小刚,杨涛. 2001. 鲜水河断裂带乾宁—康定段的滑动速率与强震复发间隔[J]. 地震学报,23(3):250–261. doi: 10.3321/j.issn:0253-3782.2001.03.004 Zhou R J,He Y L,Huang Z Z,Li X G,Yang T. 2001. The slip rate and strong earthquake recurrence interval on the Qianning-Kangding segment of the Xianshuihe fault zone[J]. Acta Seismologica Sinica,23(3):250–261 (in Chinese).
Abdelrahman E M,Bayoumi A I,Abdelhady Y E,Gobashy M M,El-Araby H M. 1989. Gravity interpretation using correlation factors between successive least-squares residual anomalies[J]. Geophysics,54(12):1614–1621. doi: 10.1190/1.1442629
Allen C R,Luo Z L,Qian H,Wen X Z,Zhou H W,Huang W S. 1991. Field study of a highly active fault zone: The Xianshuihe fault of southwestern China[J]. Geolog Soc Am Bull,103(9):1178–1199. doi: 10.1130/0016-7606(1991)103<1178:FSOAHA>2.3.CO;2
Bai D H,Unsworth M J,Meju M A,Ma X B,Teng J W,Kong X R,Sun Y,Sun J,Wang L F,Jiang C S,Zhao C P,Xiao P F,Liu M. 2010. Crustal deformation of the eastern Tibetan Plateau revealed by magnetotelluric imaging[J]. Nat Geosci,3(5):358–362. doi: 10.1038/ngeo830
Beaumont C,Jamieson R A,Nguyen M H,Lee B. 2001. Himalayan tectonics explained by extrusion of a low-viscosity crustal channel coupled to focused surface denudation[J]. Nature,414(6865):738–742. doi: 10.1038/414738a
BGI. 2012. Gravity databases[EB/OL]. [2017−09−06]. http://bgi.omp.obs-mip.fr/.
ISC. 2016. On-line bulletin[EB/OL]. [2018−03−26]. http://www.isc.ac.uk.
Royden L H,Burchfiel B C,King R W,Wang E,Chen Z L,Shen F,Liu Y P. 1997. Surface deformation and lower crustal flow in eastern Tibet[J]. Science,276(5313):788–790. doi: 10.1126/science.276.5313.788
-
期刊类型引用(9)
1. 孟上九，秦艺峰，王淼，孙志远，穆海龙. 矿震机理及其次生灾害监测技术研究进展. 世界地震工程. 2024(03): 103-118 . 
百度学术
2. 王忠忠，庄卓涵，胡飞跃，黄文龙. 广州北部丘陵区岩溶塌陷形成条件与易发性评价. 中国地质灾害与防治学报. 2024(04): 163-172 . 百度学术
3. 任涛，田国亮，宁志杰，周爱红，李宽，陈石. 基于地理探测器和随机森林的岩溶塌陷易发性评价. 灾害学. 2023(03): 227-234 . 百度学术
4. 徐玉聪，朱建，董建辉，周鲁. 三峡库区仙女山断裂周缘地震类型特征识别. 水利水电快报. 2022(02): 9-16 . 百度学术
5. 吴涛，庞聪，江勇，丁炜，廖成旺. 基于随机子空间和Ada Boost集成学习的地震事件性质辨识研究. 地球物理学进展. 2022(03): 981-988 . 百度学术
6. 尹欣欣，蔡润，陈文凯，彭立顺，李少华. 甘肃塌陷地震震相识别研究. 地震工程学报. 2021(02): 259-263 . 百度学术
7. 周少辉，蒋海昆，曲均浩，李健，郭宗斌，郑旭. 爆破、塌陷识别研究进展综述. 中国地震. 2021(02): 508-522 . 百度学术
8. 赵胤翔，赵金昌，马忠忠，常乐，李春明. 锚杆无损检测反射信号分形分析及其应用研究. 中国矿业. 2021(11): 101-108 . 百度学术
9. 梁飞，阎春恒，向巍，黄惠宁，原永东，李莎. 广西地区不同类型地震震源参数特征. 华北地震科学. 2020(S2): 33-38 . 百度学术
其他类型引用(3)
下载:
计量
- 文章访问数: 1819
- HTML全文浏览量: 1059
- PDF下载量: 107
- 被引次数: 12
