Deep structure of eastern margin of Bayan Har block and its adjacent areas by using teleseismic P-wave tomography
-
摘要: 为了揭示巴颜喀拉地块东缘及邻区的壳幔速度结构差异，获取2017年九寨沟MS7.0地震的深部构造背景，本文收集了2009年5月至2016年8月期间四川及邻区数字测震台网的203个地震台站所记录到的远震P波走时数据，应用有限频体波走时层析成像方法，反演得到了巴颜喀拉地块东缘及邻区50—600 km深度范围内的三维壳幔P波速度结构。反演结果表明：巴颜喀拉地块东缘及邻区的壳幔速度结构具有明显的横向不均匀性和分区特征，松潘—甘孜地槽褶皱系、西秦岭和祁连山褶皱系的整体速度异常较低，研究区东部具有克拉通性质的四川盆地西北缘和鄂尔多斯地块南缘则呈明显的高速异常。上地幔P波速度结构特征差异表明松潘—甘孜地块的抬升可能与地幔上涌有关，巴颜喀拉地块东缘九寨沟震区及周边50—250 km深度范围内的上地幔存在低速异常，在400—600 km地幔过渡带深度范围内表现为明显的高速异常特征。巴颜喀拉地块向东南方向运移受到东部高速、高强度的扬子克拉通地块对青藏高原物质东向挤出的强烈阻挡，而九寨沟震区处于松潘—甘孜地块重要的北东边界断裂交会处附近，应力容易在此集中，这些因素均可能是东昆仑断裂塔藏段与岷江断裂北段交会处附近发生九寨沟MS7.0地震的深部动力学背景。
-
关键词:
- 有限频走时层析成像 /
- 巴颜喀拉地块 /
- 九寨沟MS7.0地震 /
- P波速度结构 /
- 松潘—甘孜地槽褶皱系
Abstract: This paper collected the teleseismic P-wave travel data recorded by 203 broadband stations in digital seismic network of Sichuan and its neighboring areas from May 2009 to August 2016. And by using the finite-frequency tomography this paper carried out the inversion for P-wave velocity structure of the crust and upper mantle in the depth of 50−600 km in the eastern margin of Bayan Har block and its surrounding areas. The results show that the velocity structure of the Bayan Har block in the eastern part of the Qinghai-Xizang （Tibetan） Plateau and the surrounding crustal-mantle is characterized by obviously lateral inhomogeneity and zoning. Specifically, the Songpan-Garze trough fold system, the western Qinling and Qilianshan fold systems in the eastern margin of the Tibetan Plateau displayed low velocity, whereas the northwestern margin of the Sichuan basin and the southern edge of the Ordos block in the eastern part of the studied area showed obvious high velocity anomalies. P wave velocity structure characteristics of upper mantle suggested that the uplift of Songpan-Garze block is related to upper mantle upwelling. There are low velocity anomalies in the depth range of 50−250 km in Jiu-zhaigou earthquake area and the surrounding upper mantle in the east margin of Bayan Har block. The three-dimensional P-wave velocity structure also reveals that the source region of Jiuzhaigou MS7.0 earthquake shows a high-velocity anomaly at the mantle transition zone with depth range of 400−600 km. Therefore, it was deduced that this anomaly may be caused by upwelling of hot asthenosphere material. With the crustal thickening and strong uplift of the Tibetan Plateau, the Sichuan-Qinghai block in the west of Longmenshan fault zone slipped and was pushed to the SE direction, which was influenced by the blocking of high velocity and high strength craton blocks to the eastern Qinghai-Tibet Plateau extrusion, its hard upper crust should form a huge stress accumulation, and Jiuzhaigou is located in the vicinity of the north-east boundary fault intersection at the Songpan-Garze block, where the stress is likely to be concentrated. These factors may result in the occurrence of the Jiuzhaigou MS7.0 earthquake near the East Kunlun fault zone and the northern section of Minjiang fault. -
引言
与地震相关的热辐射研究始于20世纪80年代（Gorny et al，1988；Tronin，1996），自此国内外不少研究人员进行了与此相关的研究（强祖基等，1990；刘德富等，1999；Tramutoli et al，2001；Saraf，Choudhury，2003；Choudhury et al，2006）。地表温度场中存在着与活动构造带密切相关的热信息，地震前后温度场的变化与同震形变观测结果相吻合（马瑾等，2010；陈顺云等，2014），与地震活动相关的热异常也已得到论证（邓志辉等，2003；Ouzounov，Freund，2004）。Tronin等（2004）研究了2003年12月26日伊朗巴姆MW6.6地震得出，震前4天出现陆面温度（land surface temperature，简写为LST）增温，幅度为7°C—10°C；Saraf等（2012）分析2011年1月19日达尔本丁MS7.2地震得出震前两天LST增温异常幅度达8°C—10°C，增温面积约为6×104 km2。近年来国内研究人员在此方面的研究取得了一定进展（魏乐军等，2008；康春丽等，2011；解滔等，2015），大量震例的研究结果表明，强震前确实存在热辐射异常。张元生等（2010）根据地表热辐射影响因素所具有的不同周期，利用小波变换法去除影响因素，从而有效地提取地震热异常，并对2008年汶川、于田、乌恰地震进行了热辐射变化分析，分析结果显示这3次地震在发震前均呈现明显的热异常特征，进而总结出发震区域一般位于异常区边缘及区内的活动断层上，发震时间在功率谱达到峰值之后的35天内。其它的相关研究还有很多（张璇等，2013；李青梅等，2015），这些研究成果在发震地点和发震时间上意义明确，也可为短期地震预测提供经验。
在大量震例研究的基础上，目前已形成了多种热异常提取方法，这些方法主要是通过对比地震前后卫星热红外亮温图像的变化来识别异常。随着该领域研究的不断深入，又发展出一些新的研究方法。例如：Tramutoli等（2001）提出的“稳健卫星技术”（robust satellite techniques，简写为RST）反映了时空域上的地震热辐射异常信息；陈梅花等（2007）将“断层带内外温差法”用于研究1997年玛尼MS7.5地震，尽管受到降雪等因素的干扰，但还是有效地提取到了地震前的增温信息，震后相对亮温差恢复至正常状态；郭晓等（2010）通过研究强震长波辐射异常变化时提出“功率谱相对变化法”，并使用该方法对中国大陆6次MS6.5地震进行分析，结果表明震前均出现了不同程度的功率谱幅值增强现象；温少妍（2011）利用“历年同期亮温偏移指数法”研究2008年汶川MS8.0地震和2010年玉树MS7.1地震，得出了二者亮温异常的时空演化规律；解滔等（2015）使用“连续小波变换法”对于田地震的亮温异常进行了研究。地震热红外异常成因机理方面的探索也取得了一定进展，但目前尚未形成统一认识，主要有以下几种观点：① 地球放气-电磁增温说。震前应力积累变化使地下气体（如CO2，CH4）沿微裂隙溢出，辐射增强，另外地壳活动强烈时电磁场会出现异常，和太阳辐射一并促进温室气体的热辐射，进而在震源区附近出现大面积、高幅度热异常（强祖基等， 1992，1998）；② 应变能转换说。吴立新等（2001）的撞击瞬态监测试验显示撞击瞬间出现明显的升温现象，且震源区岩层蠕动使断裂带两侧摩擦生热也会引起地面升温；刘培洵等（2004）研究得出岩石破裂前在未来断层处出现红外热像，岩石弹性形变压力与温度表现为正相关，即压力越大温度增加越快；③ 地球温室效应-大气耦合说。地球放气温室效应为地震热辐射异常的基本成因，间接成因与大气中水汽含量有关，地球放气形成相对高温区，加速水汽蒸发，在高空水蒸气遇冷凝结，这一机理可以解释多数大地震发生后震中附近下大雨或下大雪的自然现象（张元生等，2010）。
本文基于中国静止气象卫星FY-2C/E/G亮温数据产品，使用小波变换和功率谱估计法对2017年九寨沟MS7.0地震前的资料进行扫描，发现沿龙门山断裂带存在一条热红外增温带。九寨沟地震的发震断层为东昆仑断裂的一个分支，即介于岷江断裂与塔藏断裂之间的树正断裂。从构造上来说，汶川地震与九寨沟地震是处于不同构造上的两次地震，但对比汶川地震的热异常区域，为何两次地震的热红外增温区域表现出了一定的相似性？本文基于该问题，拟分析讨论九寨沟MS7.0地震前的热红外亮温变化过程，类比分析盆地周缘MS≥6.5地震的热辐射异常形态，结合地球放气理论尝试性地解释盆地周缘地震的热辐射异常现象。
1. 数据与方法
本文从中国静止气象卫星FY-2C/E/G亮温小时数据产品中选取午夜（凌晨1时至5时）5个时次的数据，以避免太阳直接辐射造成的影响，并对5个时次数据求均值后采用补窗法进行简单的去云处理，形成亮温日值数据。
数据处理方法为小波变换和功率谱估计法。小波变换是在傅里叶变换的基础上发展演变而来的，傅里叶变换适用于稳态信号的频域分析，而自然界中的信号多为非稳态信号，信号中不同频率成分的出现时间无法判断，短时傅里叶变换便应运而生，该变换在一定程度上克服了上述不足，但一旦加窗，又无法适应所有的频率成分，对任一时间点上的频率分量依然无从知晓。小波变换则弥补了上述不足，其最大的特点是具有自适应性和数学显微镜性，在几乎不损失原始信号的情况下对原始信号进行重构，能够聚焦到信号的任何细节，是当前信号时频分析的有效工具。有限时间序列的小波变换公式为
$${W_\psi }f(a, b) {\text{＝}} \frac{1}{{\sqrt a }}\int_{ {\text{－}} \infty }^\infty {f(t)} \psi \left(\frac{{t {\text{－}} b}}{a}\right){\rm{d}}t{\text{，}}$$ (1) 式中：
$f(t)$ 为原始信号；${W_\psi }f(a, b)$ 为小波变换后的各频段分量；a为尺度因子，控制小波函数的伸缩，对应于变量的频率；b为时间平移因子，控制小波函数的平移；$\displaystyle\frac{1}{{\sqrt a }}\psi \left(\frac{{t - b}}{a}\right)$ 为小波母函数，文中以Daubecheies小波系中的db8小波基函数对原始亮温数据进行小波变换处理。以九寨沟地震为例，图1给出了震中（33.2°N，103.8°E）处2014—2017年的原始亮温（brightness temperature，简写为BT）值（图1a）及小波变换后的信息（数据更新至2017年9月10日，图1b为2阶小波尺度部分）。从图1可以看出，7阶小波尺度部分（图1d）的年变特征明显，因此可将地球年变温度场等长周期成分通过去除7阶小波尺度部分来消除其影响，云雨、寒热气流引起的温度变化时间较短，一般为几小时至几天，这部分影响因素可通过去除2阶小波细节部分（图1c）来消除，处理后的数据为亮温波形数据，单位为°C （图1e）。
![]()
图 1 热红外原始亮温（BT）数据及其经小波变换后的信息（a） 原始亮温数据；（b） 2阶小波尺度部分；（c） 2阶小波细节部分；（d） 7阶小波尺度部分；（e） 2阶小波尺度部分与7阶小波尺度部分之差Figure 1. The results of wavelet transformation of thermal infrared brightness temperature （BT） data（a） Original record of brightness temperature；（b） Results with 2nd-order scale analysis；（c） Details of 2nd-order wavelet；（d） Results with 7th-order scale analysis；（e） Difference between 2nd-order scale analysis and 7th-order scale analysis功率谱估计分为经典谱估计和现代谱估计，而经典谱估计主要有周期图和自相关两种方法，Welch算法是在周期图法的基础上改进而来的。本文中的功率谱估计参考了Welch算法，先将长度为N的数据分为L段，每段M个点，再对每段数据加窗求傅里叶变换，最后计算各段功率谱的平均值，其计算式为
$$ \overline {{P_x}(w)} {\text{＝}} \frac{1}{L}\sum\limits_{i {\text{＝}}1}^L {\frac{1}{{MU}}} {\left| {\sum\limits_{n {\text{＝}} 0}^{M {\text{－}} 1} {{x_i}(n)w(n)} {{\rm{e}}^{ {\text{－}} jwn}}} \right|^2}{\text{，}} \quad\quad i {\text{＝}} 1{\text{，}} 2{\text{，}} \cdot \cdot \cdot {\text{，}} L{\text{，}} $$ (2) 式中，
$U {\text{＝}} \displaystyle\frac{1}{M}\sum\limits_{n {\text{＝}} 0}^{M {\text{－}} 1} {{w^2}} (n)$ 为归一化因子，$w(n)$ 为窗函数，$\overline {{P_x}(w)} $ 为功率谱平均值的最后计算结果。计算上述小波变换及去除影响因素之后数据的功率谱，以n＝64天为窗长，m＝1天为滑动窗长作傅里叶变换，时间约定为窗内数据的最新时间，对每个像元的时程数据每滑动一次可得一组功率谱。为了对比地震前后功率谱的异同，对其进行幅值相对处理，最后得到的结果为各像元（0.05°×0.05°）的相对功率谱幅值。功率谱计算及相对处理过程见郭晓等（2010），在此基础上本文计算了相对功率谱背景值和标准差，详细计算过程见张丽峰等（2016）。
2. 结果分析
2.1 九寨沟地震热红外时空变化
2017年8月8日四川阿坝州九寨沟县发生了MS7.0地震，震源深度为20 km，震中为（33.2°N，103.8°E）。在九寨沟地震前的日常资料扫描分析中，发现沿龙门山断裂带存在一条热红外增温带（图2）。
![]()
图 2 2017年7月23日中国中西部的热红外异常空间分布图F1：塔臧断裂；F2：树正断裂；F3：岷江断裂；F4：虎牙断裂；F5：雪山梁子断裂Figure 2. Spatial distribution of thermal infrared anomaly in central-western China on July 23，2017F1：Tazang fault；F2：Shuzheng fault；F3：Minjiang fault；F4：Huya fault；F5：Xueshanliangzi fault对上述区域进行时空扫描，观察到在7月初沿龙门山断裂带即有微弱的热辐射异常存在（图3）；随着时间的推移，热辐射的幅度和面积都逐渐增大，至7月23日辐射面积达到最大，大于6倍平均值的面积大概为6×104 km2。本次地震与张元生等（2010）对汶川地震前的分析结果在空间上存在相似的辐射异常区，汶川地震在功率谱达到峰值后11天发震。鉴于在汶川地震中积累的经验，持续跟进分析资料（更新至7月30日），观察到辐射增温区域相比7月中旬有所收缩，但仍然清晰可见。随后8月8日九寨沟县发生了MS7.0地震，该地震处于增温区域边缘。震后持续补充资料，结果显示：8月2日辐射增温区依旧沿龙门山断裂带展布，辐射面积有所增大，之后收缩，直至8月6日辐射面积缩至最小，8月7日辐射面积又有回返增大趋势，8月8日发震，之后辐射面积持续增大，直到8月14日转为减小；辐射增温区后续向南收缩，面积持续减小，直到9月10日热辐射异常几乎消失。
从构造空间看，本次热辐射异常最初出现在龙门山断裂带附近，于7月20日局部扩展至树正断裂附近并逐渐增强。异常整体显示为增强趋势，其展布与龙门山断裂带走向一致，主要分布于四川盆地西北边缘（图3）。7月23日以后异常辐射面积收缩至龙门山断裂带南段，到8月6日辐射面积减至近几日的最小值，7日略有增大，8日发生地震，之后异常持续向四川盆地内部及龙门山断裂带北段扩展。8月14日以后异常呈减弱趋势，收缩于龙门山断裂带南段与鲜水河断裂的交会区域，9月中旬基本消失。本次九寨沟地震前后的热辐射异常的空间演化与汶川地震极其类似，均经历了从龙门山断裂带北段扩展到四川盆地最后收缩于龙门山断裂带南段这一过程。异常随时间的演变可能与应力的积累变化过程有关，应力持续变化导致地热资源丰富的区域逸出气体增多（郭晓等，2014），这可能是上述相似性出现的原因。
![]()
图 3 九寨沟地震前后相对功率谱的时空变化过程Figure 3. Spatio-temporal evolution of relative power spectrum for Jiuzhaigou earthquake2.2 九寨沟地震时序曲线特征
提取异常显著区域0.5°×0.5°范围内（图3中2017−07−20小图的黑色方形）的相对功率谱均值，在此基础上计算其功率谱背景值和标准差，生成时序曲线。图4中的功率谱背景值指每年相同天数的相对功率谱均值，标准差在此基础上计算而得。
![]()
图 4 九寨沟地震区的相对功率谱时序曲线Figure 4. Time-series curves of relative power spectrum of anomaly area for Jiuzhaigou earthquake图4显示出相对功率谱在地震前出现了3次峰值（蓝色曲线），最近一次出现在震前6天（8月2日），相对功率谱幅值为平均值的11.5倍，截至目前幅值大于8倍的持续时间为40天，震后两天幅值则变为转折上升。九寨沟地震前后的相对功率谱幅值变化起伏较大，在往年无地震的时候是否也存在这种功率谱幅值大幅度变化的现象？为了对比研究这个问题，文中计算了相对功率谱的两组背景值及标准差，其中一组背景值及标准差的计算结果包含九寨沟地震前后的数据(黄色和紫色曲线)，另一组由2007—2016年10年数据计算得到(绿色和红色曲线)。结果显示：① 7月之前，两组背景值和标准差曲线各自几乎完全重合，说明2017年1—7月之间的数据对背景值的贡献不大；② 两组曲线在5—6月有所上升，这表示往年资料在此期间其相对功率谱幅值也有所上升，但是上升幅度与7—8月数据（2008—2017年10年数据）相比较低，这说明往年资料的上升幅度与九寨沟地震期间出现的这种大幅度变化不同，这也使得两组曲线在7—8月各自出现分离（紫色与红色曲线分离，黄色与绿色曲线分离）；③ 这种偏离直到9月初又趋于一致，由九寨沟地震引起的相对功率谱幅值的显著变化消失。
3. 对盆地效应的探讨
3.1 盆地与其周缘地震热辐射关系
九寨沟地震的热辐射异常边界经过龙门山断裂带，与龙门山断裂带走向一致，主要分布在四川盆地西北边缘，异常在南部和西北部有小的分支，西北部分支恰好覆盖了岷江、虎牙及雪山梁子等断裂，九寨沟地震即发生在西北部异常边缘（图5e）。汶川地震的热辐射以龙门山断裂带为边界，主要分布于四川盆地北部及以外，汶川地震发生在异常最西部的顶角上（图5b）。这两次地震的热辐射异常在空间上均与龙门山断裂带走向一致，形态相似。汶川地震发生在龙门山推覆构造带中段（徐锡伟等，2013），而九寨沟地震的发震断层为东昆仑断裂的一个分支，即介于岷江断裂与塔藏断裂之间的树正断裂，从构造上来说汶川地震与九寨沟地震是两次独立的地震，但两者的热辐射异常区域表现出了一定的相似性。
![]()
图 5 6次盆地周缘地震的热辐射异常空间分布（a） 2008年3月21日于田MS7.3地震；（b） 2008年5月12日汶川MS8.0地震；（c） 2012年6月30日新源MS6.6地震； （d） 2015年7月3日皮山MS6.5地震；（e） 2017年8月8日九寨沟MS7.0地震； （f） 2017年8月9日精河MS6.6地震Figure 5. Spatial distributions of thermal infrared anomalies for six earthquakes around the basins（a） Yutian MS7.3 earthquake on March 21，2008；（b） Wenchuan MS8.0 earthquake on May 12，2008；（c） Xinyuan MS6.6 earthquake on June 30，2012；（d） Pishan MS6.5 earthquake on July 3，2015； （e） Jiuzhaigou MS7.0 earthquake on August 8，2017；（f） Jinghe MS6.6 earthquake on August 9， 2017基于这种现象，本文分析了2008年以来新疆地区盆地周边地震所对应的热辐射状态，结果显示2008年于田MS7.3地震、2012年新源MS6.6地震和2015年皮山MS6.5地震的热辐射异常形态均与盆地相关，如图5a，c，d所示，热辐射异常均贯穿了塔里木盆地，但从发震构造上来说这是三次独立地震，可见新疆塔里木盆地在大震前也存在类似现象。九寨沟地震后一天就发生了新疆精河MS6.6地震，该地震处于准噶尔盆地边缘，图5f为精河地震前的热红外辐射空间展布。表1统计了上述地震与相关盆地之间的一些参数，从各震例所对应的热辐射异常面积的分配来看，盆地一侧的面积占比均超过50%，其中于田、皮山及精河地震的占比更是高达100%，似乎异常面积越大，处于盆地内部的异常占比也越高。
表 1 盆地周缘地震的热辐射与对应盆地之间的相关性参数Table 1. The correlation parameters between the thermal radiation anomalies for six earthquakes and the corresponding basins发震日期 地点 MS 特征
周期/d发震构造 响应盆地 盆地性质 异常最大
面积/（104 km2）盆地一侧
占比相对盆地
位置2008−03−21 于田 7.3 64 郭扎错断裂 塔里木盆地 含油气田 55 100% 盆地内及边缘 2008−05−12 汶川 8.0 13 龙门山断裂带 四川盆地 含油气田 12 66% 盆地北部及以外 2012−06−30 新源 6.6 11 伊犁盆地北缘断裂 塔里木盆地 含油气田 60 83% 盆地中东部及以外 2015−07−03 皮山 6.5 32 泽普断裂 塔里木盆地 含油气田 44 100% 盆地内及边缘 2017−08−08 九寨沟 7.0 11 树正断裂 四川盆地 含油气田 8 50% 盆地西北及以外 2017−08−09 精河 6.6 21 库松木契克山前断裂 准噶尔盆地 含油气田 25 100% 盆地及以外 这些与异常相关的盆地均蕴藏着丰富的天热气，含油气盆地的大量天然气对震前应力的变化较为敏感，当应力积累到一定程度，盆地周缘的活动构造带及一些微裂隙均成为天然气上涌的通道，溢出地表的甲烷、二氧化碳等温室气体的辐射增温效果明显。大震前，盆地出现大片区域的辐射增温异常可能与天然气外泄有关，这或许可以解释为何不同发震构造上的地震却呈现出相似的辐射增温异常分布。这是否是与地震有关的“盆地效应”？为了说明地震盆地效应的可能性，本文分析了2008年以来四川和新疆盆地内的热辐射异常变化情况，结果显示研究区域内共出现8次类似盆地形态的异常，对应于盆地周缘MS6.5以上地震，8个震例中的6次地震存在此类异常，占比为75%。以往地震热辐射研究均从地震出发寻找异常，反过来由异常来对应地震的研究是对异常的一种判断，上述研究结果表明地震盆地效应是有一定可信度的。
此外，强祖基等（1994）的研究结果表明，油气聚集带与地震前的卫星热红外增温异常有较好的对应关系；黄福林等（1998）认为地震前油气盆地区低空大气中的甲烷浓度高于平均值两倍多，二氧化碳高于平均值2.6倍；卢振权等（2005）经过长期的观察研究认为中国近海海域临震前的卫星热红外增温异常与已知油气盆地、热水盆地等的分布区一致；郭卫英等（2006）关于2003年伽师—巴楚MS6.8地震的研究显示，震前的红外增温区就出现在塔里木盆地西北部边缘，并指出该增温现象可能与孕震机制有关。这些研究均表明地震前的卫星热辐射增温异常与含地热资源丰富的盆地有一定关联。
3.2 对地震预测的启示
图6为图5中震中附近黑色方块所示区域一年时间尺度相对功率谱幅值的变化曲线，由表2所列地震的相对功率谱时序曲线总结出一些特征：相对功率谱峰值均在7倍以上；峰值大都出现在震前，最长为震前8天，最近为震前4天。值得指出的是，地震孕育及临震阶段的应力变化是一个复杂的过程，可以通过GPS等手段来观测地表形变，但是对于地下微弱细小的变化和破裂仍难以捕捉。另外，究竟应力积累到什么程度才会发震，目前尚难以判断。盆地在震前所表现出的热辐射变化信息或许是对地下复杂精细变化的宏观反映，因此分析盆地周缘地震的热辐射异常时序变化，对地震短临预测的作用更显著一些。至于发震地点与其对应热辐射异常的位置关系，邓志辉等（2003）的研究认为强震震中常位于红外增温异常区域的边缘或附近，张元生等（2010）的研究也得出相似的结论，从本文中给出的6个震例（图5）来看，震中也均处在热辐射异常区域的边缘。
![]()
图 6 6次盆地周缘地震的相对功率谱时序曲线Figure 6. Time-series curves of relative power spectrum for six earthquakes around basin表 2 地震相对功率谱时序曲线特征Table 2. The features of time-series curves of relative power spectrum for six earthquakes发震日期 地点 MS 相对功率谱峰值 峰值出现时间 相对功率大于2的持续时间/d 2008−03−21 新疆于田 7.3 9.8 3月13日（震前8天） 105 2008−05−12 四川汶川 8.0 7.4 5月8日（震前4天） 81 2012−06−30 新疆新源 6.6 16.0 6月25日（震前5天） 75 2015−07−03 新疆皮山 6.5 13.8 8月25日（震后53天） 33 2017−08−08 四川九寨沟 7.0 11.5 8月2日（震前6天） 73 2017−08−09 新疆精河 6.6 13.5 8月1日（震前8天） 52 4. 讨论与结论
对2017年8月8日四川九寨沟地震热辐射异常的研究结果显示：其空间展布与龙门山断裂带走向一致，主要分布在四川盆地边缘，与汶川地震前热辐射异常的展布特征相似，热辐射异常面积最大可达6×104 km2，地震发生在增温区域的边缘；其时序演化特征显示，相对功率谱幅值于7月初迅速升高后，长时间处于高幅值状态，震前最后一次峰值为均值的11.5倍，尽管有所转折，但幅值均在9倍以上，出现3次波峰值后发震，与最后一次波峰之间相隔6天；7月之前，两组相对功率谱背景值和标准差曲线各自几乎完全重合，7月之后，由于九寨沟地震的发生，两组背景值和标准差数据明显偏离，即九寨沟地震的相对功率谱幅值与过去十年的数据相比较为突出。
文中通过6个震例分析了盆地周缘地震的热辐射异常与盆地的相关性，从各震例所对应的热辐射异常面积来看，盆地一侧的面积占比均超过50%，其中于田、皮山及精河地震的占比更是高达100%，似乎异常面积越大，处于盆地内部的异常占比也越高。为了说明地震盆地效应的可能性，文中分析了自2008年以来四川和新疆盆地内的热辐射异常变化情况，结果显示研究区域内共出现8次类似盆地形态的异常，对应于盆地周缘MS≥6.5的地震，共8个震例中的6次地震存在此类异常，占比为75%。这一结果表明地震盆地效应具有一定可信度。
地震盆地效应可以解释盆地周缘具有不同发震构造的地震却表现出类似的盆地形态异常。地震盆地效应是从盆地周缘地震的热辐射异常现象出发，为地震热红外成因机理积极寻求答案，这明确地支持了地球放气理论。尽管如此，目前对于热红外成因机理的认识尚处于探索阶段，后续仍需进行建模和定量化计算研究，以进一步推进此方面的工作。
感谢中国地震局兰州地震研究所张元生研究员提供了热红外数据分析计算软件。
-
![]()
图 1 巴颜喀拉地块及邻区的构造背景及地震分布
Figure 1. Tectonic settings and distribution of earthquakes around the Bayan Har block and its surroundings
![]()
图 3 检测板测试结果
（a） 不同深度h上的水平向分辨率；（b） 沿不同经纬度的垂向分辨率
Figure 3. Test results of checkboard test
（a） Resolution in the horizontal direction at different depths h；（b） Resolution in the vertical direction along different longitudes and latitudes
![]()
图 3 检测板测试结果
（a） 不同深度h上的水平向分辨率；（b） 沿不同经纬度的垂向分辨率
Figure 3. Test results of checkboard test
（a） Resolution in the horizontal direction at different depths h；（b） Resolution in the vertical direction along different longitudes and latitudes
![]()
图 4 剩余走时残差与速度模型间的权衡曲线
图中数字为反演过程中采用的阻尼值
Figure 4. The trade-off curve between traveltime residual variances and velocity anomaly
The figures denote the damping values used in inversion
![]()
图 5 反演前（a）、后（b）走时残差对比图
Figure 5. Distribution of the travel-time residuals before （a） and after （b） inversion
![]()
图 6 不同深度h上的P波速度扰动结果
Figure 6. P wave velocity disturbance results at different depth h
（a） h＝50 km；（b） h＝150 km；（c） h＝200 km；（d） h＝250 km；（e） h＝300 km；（f） h＝400 km；（g） h＝500 km；（h） h＝600 km
-
邓起东,高翔,陈桂华,杨虎. 2010. 青藏高原昆仑-汶川地震系列与巴颜喀喇断块的最新活动[J]. 地学前缘,17(5):163–178. Deng Q D,Gao X,Chen G H,Yang H. 2010. Recent tectonic activity of Bayankala fault-block and the Kunlun-Wenchuan earthquake series of the Tibetan Plateau[J]. Earth Science Frontiers,17(5):163–178 (in Chinese).
丁志峰,何正勤,孙为国,孙宏川. 1999. 青藏高原东部及其边缘地区的地壳上地幔三维速度结构[J]. 地球物理学报,42(2):197–205. doi: 10.3321/j.issn:0001-5733.1999.02.007 Ding Z F,He Z Q,Sun W G,Sun H C. 1999. 3-D crust and upper mantle velocity structure in eastern Tibetan Plateau and its surrounding areas[J]. Chinese Journal of Geophysics,42(2):197–205 (in Chinese).
杜方,龙锋,阮祥,易桂喜,宫悦,赵敏,张致伟,乔慧珍,汪智,吴江. 2013. 四川芦山7.0级地震及其与汶川8.0级地震的关系[J]. 地球物理学报,56(5):1772–1783. doi: 10.6038/cjg20130535 Du F,Long F,Ruan X,Yi G X,Gong Y,Zhao M,Zhang Z W,Qiao H Z,Wang Z,Wu J. 2013. The M7.0 Lushan earthquake and the relationship with the M8.0 Wenchuan earthquake in Sichuan,China[J]. Chinese Journal of Geophysics,56(5):1772–1783 (in Chinese). doi: 10.6038/cjg20130535
杜方,闻学泽,冯建刚,梁明剑,龙锋,吴江. 2018. 六盘山断裂带的地震构造特征与强震危险背景[J]. 地球物理学报,61(2):545–559. doi: 10.6038/cjg2018L0181 Du F,Wen X Z,Feng J G,Liang M J,Long F,Wu J. 2018. Seismo-tectonics and seismic potential of the Liupanshan fault zone (LPSFZ),China[J]. Chinese Journal of Geophysics,61(2):545–559 (in Chinese). doi: 10.6038/cjg2018L0181
郭飚,刘启元,陈九辉,刘立申,李顺成,李昱,王峻,齐少华. 2009. 川西龙门山及邻区地壳上地幔远震P波层析成像[J]. 地球物理学报,52(2):346–355. Guo B,Liu Q Y,Chen J H,Liu L S,Li S C,Li Y,Wang J,Qi S H. 2009. Teleseismic P-wave tomography of the crust and upper mantle in Longmenshan area,west Sichuan[J]. Chinese Journal of Geophysics,52(2):346–355 (in Chinese).
郭慧丽,丁志峰,徐小明. 2017. 南北地震带北段的远震P波层析成像研究[J]. 地球物理学报,60(1):86–97. Guo H L,Ding Z F,Xu X M. 2017. Upper mantle structure beneath the northern South-Nouth Seismic Zone from teleseismic traveltime data[J]. Chinese Journal of Geophysics,60(1):86–97 (in Chinese). doi: 10.6038/cjg20170108
嘉世旭,张先康. 2008. 青藏高原东北缘深地震测深震相研究与地壳细结构[J]. 地球物理学报,51(5):1431–1443. doi: 10.3321/j.issn:0001-5733.2008.05.016 Jia S X,Zhang X K. 2008. Study on the crust phases of deep seismic sounding experiments and fine crust structures in the northeast margin of Tibetan Plateau[J]. Chinese Journal of Geophysics,51(5):1431–1443 (in Chinese).
嘉世旭,林吉焱,郭文斌,赵楠,邱勇. 2017. 巴颜喀拉块体地壳结构多样性探测[J]. 地球物理学报,60(6):2226–2238. doi: 10.6038/cjg20170616 Jia S X,Lin J Y,Guo W B,Zhao N,Qiu Y. 2017. Investigation on diversity of crustal structures beneath the Bayan Har block[J]. Chinese Journal of Geophysics,60(6):2226–2238 (in Chinese). doi: 10.6038/cjg20170616
李松林,张先康,张成科,赵金仁,成双喜. 2002. 玛沁—兰州—靖边地震测深剖面地壳速度结构的初步研究[J]. 地球物理学报,45(2):210–217. doi: 10.3321/j.issn:0001-5733.2002.02.007 Li S L,Zhang X K,Zhang C K,Zhao J R,Cheng S X. 2002. A preliminary study on the crustal velocity structure of Maqin-Lanzhou-Jingbian by means of deep seismic sounding profile[J]. Chinese Journal of Geophysics,45(2):210–217 (in Chinese).
李永华,田小波,吴庆举,曾融生,张瑞青. 2006a. 青藏高原INDEPTH-Ⅲ剖面地壳厚度与泊松比:地质与地球物理含义[J]. 地球物理学报,49(4):1037–1044. Li Y H,Tian X B,Wu Q J,Zeng R S,Zhang R Q. 2006a. The Poisson ratio and crustal structure of the central Qinghai-Xizang inferred from INDEPTH-Ⅲ teleseismic waveforms:Geological and geophysical implications[J]. Chinese Journal of Geophysics,49(4):1037–1044 (in Chinese).
李永华,吴庆举,安张辉,田小波,曾融生,张瑞青,李红光. 2006b. 青藏高原东北缘地壳S波速度结构与泊松比及其意义[J]. 地球物理学报,49(5):1359–1368. Li Y H,Wu Q J,An Z H,Tian X B,Zeng R S,Zhang R Q,Li H G. 2006b. The Poisson ratio and crustal structure across the NE Tibetan Plateau determined from receiver functions[J]. Chinese Journal of Geophysics,49(5):1359–1368 (in Chinese).
李大虎,丁志峰,吴萍萍,郑晨,叶庆东,梁明剑. 2015. 鲜水河断裂带南东段的深部孕震环境与2014年康定MS6.3地震[J]. 地球物理学报,58(6):1941–1953. doi: 10.6038/cjg20150610 Li D H,Ding Z F,Wu P P,Zheng C,Ye Q D,Liang M J. 2015. The deep seismogenic environment of the southeastern section of the Xianshuihe fault zone and the 2014 Kangding MS6.3 earthquake[J]. Chinese Journal of Geophysics,58(6):1941–1953 (in Chinese).
汤吉,詹艳,赵国泽,邓前辉,王继军,陈小斌,赵俊猛,宣飞. 2005. 青藏高原东北缘玛沁—兰州—靖边剖面地壳上地幔电性结构研究[J]. 地球物理学报,48(5):1205–1216. doi: 10.3321/j.issn:0001-5733.2005.05.032 Tang J,Zhan Y,Zhao G Z,Deng Q H,Wang J J,Chen X B,Zhao J M,Xuan F. 2005. Electrical conductivity structure of the crust and upper mantle in the northeastern margin of the Qinghai-Tibet Plateau along the profile Maqên-Lanzhou-Jingbian[J]. Chinese Journal of Geophysics,48(5):1205–1216 (in Chinese).
闻学泽,杜方,张培震,龙锋. 2011. 巴颜喀拉块体北和东边界大地震序列的关联性与2008年汶川地震[J]. 地球物理学报,54(3):706–716. doi: 10.3969/j.issn.0001-5733.2011.03.010 Wen X Z,Du F,Zhang P Z,Long F. 2011. Correlation of major earthquake sequences on the northern and eastern boundaries of the Bayan Har block,and its relation to the 2008 Wenchuan earthquake[J]. Chinese Journal of Geophysics,54(3):706–716 (in Chinese).
吴福元,黄宝春,叶凯,方爱民. 2008. 青藏高原造山带的垮塌与高原隆升[J]. 岩石学报,24(1):1–30. Wu F Y,Huang B C,Ye K,Fang A M. 2008. Collapsed Himalayan-Tibetan orogen and the rising Tibetan Plateau[J]. Acta Petrologica Sinica,24(1):1–30 (in Chinese).
徐果明,姚华建,朱良保,沈玉松. 2007. 中国西部及其邻域地壳上地幔横波速度结构[J]. 地球物理学报,50(1):193–208. doi: 10.3321/j.issn:0001-5733.2007.01.026 Xu G M,Yao H J,Zhu L B,Shen Y S. 2007. Shear wave velocity structure of the crust and upper mantle in western China and its adjacent area[J]. Chinese Journal of Geophysics,50(1):193–208 (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. Characteristic features of the Hoh Sai Hu (Kunlunshan) earthquake (MS8.1),northern Tibetan Plateau,China[J]. Seismology and Geology,24(1):1–13 (in Chinese).
徐锡伟,闻学泽,陈桂华,于贵华. 2008. 巴颜喀拉地块东部龙日坝断裂带的发现及其大地构造意义[J]. 中国科学:D辑,38(5):529–542. Xu X W,Wen X Z,Chen G H,Yu G H. 2008. Discovery of the Longriba fault zone in eastern Bayan Har block,China and its tectonic implication[J]. Science in China:Series D,51(9):1209–1223. doi: 10.1007/s11430-008-0097-1
杨峰,黄金莉,杨挺. 2010. 应用远震有限频率层析成像反演首都圈上地幔速度结构[J]. 地球物理学报,53(8):1806–1816. doi: 10.3969/j.issn.0001-5733.2010.08.006 Yang F,Huang J L,Yang T. 2010. Upper mantle structure beneath the Chinese capital region from teleseismic finite-frequency tomography[J]. Chinese Journal of Geophysics,53(8):1806–1816 (in Chinese).
杨歧焱,吴庆举,盛艳蕊,高家乙,宋键,邸龙. 2018. 张渤地震带及邻区近震体波成像及孕震环境分析[J]. 地球物理学报,61(8):3251–3262. doi: 10.6038/cjg2018L0628 Yang Q Y,Wu Q J,Sheng Y R,Gao J Y,Song J,Di L. 2018. Regional seismic body wave tomography and deep seismogenic environment beneath Zhangbo seismic belt and its adjacent area[J]. Chinese Journal of Geophysics,61(8):3251–3262 (in Chinese). doi: 10.6038/cjg2018L0628
詹艳. 2008. 青藏高原东北缘地区深部电性结构及构造涵义[D]. 北京: 中国地震局地质研究所: 59−70. Zhan Y. 2008. Deep Electric Structures Beneath the Northeastern Margin of the Tibetan Plateau and Its Tectonic Implications[D]. Beijing: Institute of Geology, China Earthquake Administration: 59−70 （in Chinese）.
詹艳,赵国泽,王立凤,王继军,陈小斌,赵凌强,肖骑彬. 2014. 西秦岭与南北地震构造带交汇区深部电性结构特征[J]. 地球物理学报,57(8):2594–2607. doi: 10.6038/cjg20140819 Zhan Y,Zhao G Z,Wang L F,Wang J J,Chen X B,Zhao L Q,Xiao Q B. 2014. Deep electric structure beneath the intersection area of West Qinling orogenic zone with North-South Seismic tectonic zone in China[J]. Chinese Journal of Geophysics,57(8):2594–2607 (in Chinese).
张风雪,吴庆举,李永华. 2013. 中国东北地区远震P波走时层析成像研究[J]. 地球物理学报,56(8):2690–2700. doi: 10.6038/cjg20130818 Zhang F X,Wu Q J,Li Y H. 2013. The traveltime tomography study by teleseismic P wave data in the Northeast China area[J]. Chinese Journal of Geophysics,56(8):2690–2700 (in Chinese).
张风雪,吴庆举,丁志峰. 2018. 青藏高原东部P波速度结构及其对高原隆升的启示[J]. 科学通报,63(19):1949–1961. Zhang F X,Wu Q J,Ding Z F. 2018. A P-wave velocity study beneath the eastern area of Tibetan Plateau and its implication for plateau growth[J]. Chinese Science Bulletin,63(19):1949–1961 (in Chinese). doi: 10.1360/N972018-00337
张风雪,吴庆举. 2019. 松辽盆地北缘的上地幔速度结构及该区火山成因探讨[J]. 地球物理学报,62(8):2918–2929. doi: 10.6038/cjg2019M0665 Zhang F X,Wu Q J. 2019. Velocity structure in upper mantle and its implications for the volcanism nearby the north edge of Songliao Basin[J]. Chinese Journal of Geophysics,62(8):2918–2929 (in Chinese). doi: 10.6038/cjg2019M0665
张培震,邓起东,张国民,马瑾,甘卫军,闵伟,毛凤英,王琪. 2003. 中国大陆的强震活动与活动地块[J]. 中国科学:D辑,33(增刊1):12–20. Zhang P Z,Deng Q D,Zhang G M,Ma J,Gan W J,Min W,Mao F Y,Wang Q. 2003. Active tectonic blocks and strong earthquakes in the continent of China[J]. Science in China:Series D,46(S2):13–24.
Dahlen F A,Hung S H,Nolet G. 2000. Fréchet kernels for finite-frequency traveltimes: Ⅰ . Theory[J]. Geophys J Int,141(1):157–174. doi: 10.1046/j.1365-246X.2000.00070.x
Densmore A L,Ellis M A,Li Y,Zhou R J,Hancock G S,Richardson N. 2007. Active tectonics of the Beichuan and Pengguan faults at the eastern margin of the Tibetan Plateau[J]. Tectonics,26(4):TC4005.
Huang Z C,Wang P,Xu M J,Wang L S,Ding Z F,Wu Y,Xu M J,Mi N,Yu D Y,Li H. 2015. Mantle structure and dynamics beneath SE Tibet revealed by new seismic images[J]. Earth Planet Sci Lett,411:100–111. doi: 10.1016/j.jpgl.2014.11.040
Hung S H,Dahlen F A,Nolet G. 2000. Fréchet kernels for finite-frequency traveltimes: Ⅱ . Examples[J]. Geophys J Int,141(1):175–203. doi: 10.1046/j.1365-246X.2000.00072.x
Hung S H,Shen Y,Chiao L Y. 2004. Imaging seismic velocity structure beneath the Iceland hot spot:A finite frequency approach[J]. J Geophys Res,109(B8):B08305.
Hung S H,Chen W P,Chiao I Y. 2011. A data-adaptive,multiscale approach of finite-frequency,traveltime tomography with special reference to P and S wave data from central Tibet[J]. J Geophys Res,116(B6):B06307.
Kirby E,Whipple K X,Burchfiel B C,Tang W Q,Berger G,Sun Z M,Chen Z L. 2000. Neotectonics of the Min Shan,China:Implications for mechanisms driving Quaternary deformation along the eastern margin of the Tibetan Plateau[J]. Geol Soc Am Bull,112(3):375–393. doi: 10.1130/0016-7606(2000)112<375:NOTMSC>2.0.CO;2
Laske G,Masters G,Ma Z T,Pasyanos M. 2013. Update on CRUST1.0:A 1-degree global model of Earth’s crust[J]. Geophys Res Abstr,15:EGU2013-2658.
Lei J S,Zhao D P. 2016. Teleseismic P-wave tomography and mantle dynamics beneath eastern Tibet[J]. Geochem Geophys Geosyst,17(5):1861–1884. doi: 10.1002/2016GC006262
Li D H,Liao H,Ding Z F,Zhan Y,Wu P P,Xu X M,Zheng C. 2018. Joint inversion of the 3D P wave velocity structure of the crust and upper mantle under the southeastern margin of the Tibetan Plateau using regional earthquake and teleseismic data[J]. Acta Geologica Sinica,92(1):16–33. doi: 10.1111/1755-6724.13492
Liu Z,Tian X B,Gao R,Wang G C,Wu Z B,Zhou B B,Tan P,Nie S T,Yu G P,Zhu G H,Xu X. 2017. New images of the crustal structure beneath eastern Tibet from a high-density seismic array[J]. Earth Planet Sci Lett,480:33–41. doi: 10.1016/j.jpgl.2017.09.048
Meyer B,Tapponnier P,Bourjot L,Métivier F,Gaudemer Y,Peltzer G,Guo S M,Chen Z T. 1998. Crustal thickening in Gansu-Qinghai,lithospheric mantle subduction,and oblique,strike-slip controlled growth of the Tibet Plateau[J]. Geophys J Int,135(1):1–47. doi: 10.1046/j.1365-246X.1998.00567.x
Paige C C,Saunders M A. 1982a. Algorithm 583:LSQR:Sparse linear equations and least squares problems[J]. ACM Trans Math Softw,8(2):195–209. doi: 10.1145/355993.356000
Paige C C,Saunders M A. 1982b. LSQR:An algorithm for sparse linear equations and sparse least squares[J]. ACM Trans Math Softw,8(1):43–71. doi: 10.1145/355984.355989
Rawlinson N,Kennett B L N. 2004. Rapid estimation of relative and absolute delay times across a network by adaptive stacking[J]. Geophys J Int,157(1):332–340. doi: 10.1111/j.1365-246X.2004.02188.x
Shen X Z,Mei X P,Zhang Y S. 2011. The crustal and upper-mantle structures beneath the northeastern margin of Tibet[J]. Bull Seismol Soc Am,101(6):2782–2795. doi: 10.1785/0120100112
Tian Y,Zhao D P,Sun R M,Teng J W. 2009. Seismic imaging of the crust and upper mantle beneath the North China Craton[J]. Phys Earth Planet Inter,172(3/4):169–182.
van der Woerd J,Tapponnier P,Ryerson F J,Meriaux A S,Meyer B,Gaudemer Y,Finkel R C,Caffee M W,Zhao G G,Xu Z Q. 2002. Uniform postglacial slip-rate along the central 600 km of the Kunlun fault (Tibet),from 26Al,10Be,and 14C dating of riser offsets,and climatic origin of the regional morphology[J]. Geophys J Int,148(3):356–388. doi: 10.1046/j.1365-246x.2002.01556.x
Wang C Y,Wu J,Lou H,Lü Z Y. 2003. P-wave crustal velocity structure in the western Sichuan and eastern Tibetan region[J]. Scice in Chin:Series D,46(Suppl):254–265.
Yao H J,Beghein C,van der Hilst R D. 2008. Surface wave array tomography in SE Tibet from ambient seismic noise and two-station analysis: Ⅱ . Crustal and upper-mantle structure[J]. Geophys J Int,173(1):205–219. doi: 10.1111/j.1365-246X.2007.03696.x
-
期刊类型引用(4)
1. 赵一帆，王毅. 基于三重震相波形模拟的东北亚地区上地幔S波速度结构研究. 物探化探计算技术. 2024(02): 137-145 . 
百度学术
2. 苏慧，魏荣强，周元泽，崔清辉，李国辉. 东北亚边缘地区地幔过渡带内滞留太平洋板片上界面的三重震相研究. 地球物理学报. 2023(06): 2431-2444 . 百度学术
3. 崔辉辉，周元泽，石耀霖，王晓冉，李国辉. 华北克拉通东部滞留板块下方低速异常的地震三重震相探测. 地球物理学报. 2016(04): 1309-1320 . 百度学术
4. 李国辉，眭怡，周元泽. 基于P波三重震相的下扬子克拉通地幔转换带顶部低速层初探. 地球物理学报. 2014(07): 2362-2371 . 百度学术
其他类型引用(11)
下载:
计量
- 文章访问数: 2129
- HTML全文浏览量: 888
- PDF下载量: 161
- 被引次数: 15
