Analysis of stress field and its causes of the Xianyou earthquake swarm in Fujian
-
摘要:
利用双差定位法对2010年8月1日至2014年12月31日福建仙游地区中国地震台网中心记录到的地震群进行重定位，并采用P波初动的方法求解出该震群的159个震源机制解，显示该震群的震源机制解结果与背景构造应力场表现明显不一致。为了进一步探究差异性存在的原因，将得到的震群震源机制解结果分区域和分深度求解应力场，结果显示：金钟水库的东南部附近应力方向更远离背景构造应力方向，并且应力形因子小于0.5，中间轴呈现为拉张应力；随着深度的增加，应力场的方向也更接近背景构造应力场的方向。此外，根据前人研究推断该区域应力方向发生偏转是由于太平洋板块在深部的低角度俯冲挤压导致浅部张开，应力呈现为东西向拉张所致。在上述两原因的共同作用下，使得经震源机制解反演出的应力场结果与背景应力场差别较大。
Abstract:This study employed the double-difference location method to relocate earthquakes recorded by the seismic network center in the Xianyou region of Fujian from August 1, 2010, to December 31, 2014. It also used the initial motion of P waves to determine the focal mechanism solutions for 159 events within this seismic cluster. The results reveal a significant inconsistency between the focal mechanism solutions of the cluster and the background tectonic stress field. To further investigate the reasons behind this discrepancy, the obtained focal mechanism solutions were analyzed to calculate the stress field by region and depth. The findings indicate that in the southeastern part, near the Jinzhong Reservoir, the direction of stress deviates more from the background tectonic stress direction, with a stress shape factor less than 0.5, and the intermediate axis showing extensional stress. As depth increases, the direction of the stress field also becomes closer to that of the background tectonic stress field. Additionally, based on previous research, it is inferred that the deviation in stress direction in this area is due to the low-angle subduction and compression of the Pacific Plate at depth, leading to extension at shallow levels, with the stress manifesting as east-west extension. Under the combined influence of these two factors, the stress field results deduced from the focal mechanism solutions differ significantly from the background stress field.
-
引言
与地震研究相关的井孔深部微小温度变化的高精度观测（分辨率优于10−3℃）始于20世纪70年代（Nakamura，Wakita，1984；Shimamura et al，1984 ）。我国自1984年开始在云南省和首都圈地区布设区域性台网，获得了1988年澜沧—耿马MS7.6等典型的中强震震例（付子忠，1988），自此之后，井水温度动态观测得到了大力发展。目前，地震系统内已拥有300多口井泉的水温观测点，而水温则已成为地下流体监测台网中继水位之后的第二大测项（冯恩国等，2012）。
国内外关于井水温度的地震前兆、同震及震后的动态变化等研究已有不少报道：例如Mogi等（1989）报道了日本伊豆半岛东北温泉区一口自流泉的温度变化，发现温度变化与地震的发生存在密切关系，水温观测记录到了一些中强震的前兆和同震变化；Cicerone等（2009）对包括水温在内的各类前兆进行了统计整理，认为地下水的流动能快速有效地传递热异常信息，是水温前兆异常变化的主要原因；Orihara等（2014）报道了2011年日本MW9.0地震3个月前，位于震中NW 155 km的岩手县一口2 000 m深的井，其井孔水位和水温同时出现大幅度下降的现象。
王瑜青等（1994）分析了云南地区9口井自1986年至1993年的水温观测资料，总结出水温正常动态可分为6种类型，前兆变化动态可分为5种类型，认为水温变化具有显著的中强震前兆映震能力；刘耀炜等（2008）报道了2007年宁洱MS6.4地震前云南地区的群体性水温异常变化；赵刚等（2009）分析整理了2008年四川汶川MS8.0地震前后全国277个观测台站的水温观测资料，发现有134个观测台站的水温观测数据存在同震变化和震后调整，有37个观测台站的水温观测数据存在“疑似”前兆变化；张彬等（2014）分析研究了2007年以前云南地区水温异常持续时间、异常空间位置、异常幅度与后续地震的关系，结果表明，水温异常主要为短临异常，但强震前也存在水温中期异常。一般情况下，地震震级越大，异常范围越大，发震地点通常出现在水温异常集中的区域。
在大量的水温前兆异常报道和研究中，大部分是在地震发生后对资料进行回溯研究的过程中发现水温异常变化，仅有少量是在地震发生前即根据包括井水温在内的异常变化提出了短临预测意见（车用太等，1998；付虹，赵小艳，2013；付虹等，2015），说明对中强震发生前水温微异常变化的分析仍然面临较大的困难。水温变化形态复杂多样，如何识别具有多样性和个性特征的异常变化是前兆分析中面临的关键问题之一。
丽江地区位于川滇块体中部，云南省的西北部，是红河、丽江—剑川、龙蟠—乔后、鹤庆—洱源、程海等多条断裂的复合交汇区（国家地震局地质研究所，云南省地震局，1990）（图1），中强震发生频繁。本文拟以丽江党校井为例，分析水温基本动态和同震变化特征，并对前兆异常变化进行讨论，得到的结果和认识可以对未来水温动态分析、异常核实跟踪以及井孔中温度探头的放置提供一定的依据和参考。
1. 观测背景
党校井位于丽江盆地东边缘，与黑龙潭泉群接近（图2），地下水类型为碳酸盐类裂隙溶洞水。补给区位于盆地北端的九子海洼地及其周边山系，该区域海拔2 800—3 588 m，高出丽江盆地400—1 188 m （图3），翟军伟等通过投放钼酸铵示踪剂证实了九子海与黑龙潭和下游泉点存在直接的水力联系，是这些泉水的补给区（康晓波等，2013；曾成等，2013）。九子海溶蚀区侵蚀作用强烈，溶蚀洼地和落水洞沿九子海环状断裂带分布。降雨产生的地表水通过洼地和落水洞灌入式补给后，沿岩溶管道、节理裂隙由北向南径流，在黑龙潭泉群出口处（海拔2 420 m）受地形切割和第四系相对隔水层阻隔，形成地下水富集带，浅层径流呈股状涌出地表，深层径流侧向补给丽江盆地底部第四系含水层后再向盆地南部作深远程径流，从九子海到黑龙潭泉群运移距离约为20 km。
![]()
图 1 云南丽江周边活动断裂构造和部分地震震源机制解活动断裂引自邓起东等 （2003）；震源机制解引自Ekström et al （2012）Figure 1. Active faults and focal mechanism solutions of some major earthquakes in Lijiang area and its vicinity，Yunnan ProvinceThe active faults are referred to Deng et al （2003）; the mechanism solutions are referred to Ekström et al （2012）![]()
图 2 党校井区域水文地质和井点位置图（修改自康晓波等，2013）Figure 2. The hydrogeological map of the Dangxiao well region （modified from Kang et al，2013 ）![]()
图 3 党校井地下水系统补径排过程示意图 （改自康晓波等，2013）Figure 3. The sketch map of recharge-flow-discharge process of the Dangxiao well （modified from Kang et al，2013 ）党校井井孔结构及观测示意图如图4a所示，井深为347.3 m，其中167.5—310 m安装滤水管，含水层为三叠系中统灰岩岩溶裂隙承压水。党校井管高出观测室地面6 m，在距地面0.4 m处安装一直径7 cm的泄流口，当水位高出泄流口时则自动泄流，称之为自流状态或动水位状态；低于泄流口停止自流，称之为非自流或静水位状态。井水位观测值为水面至泄流口的距离。
党校井自1992年3月1日起使用浮子式模拟水位仪进行观测。2007年6月改造为数字水位仪观测，分辨率为1 mm，同时增加了一套高精度石英温度仪，温度探头放置在距井口300 m深度处，采样率为分钟值。党校井分别于2006年11月24日和2015年7月17日进行了水温梯度测试（图4b）。
![]()
图 4 党校井井孔柱状图与观测示意图 （a） 和井温梯度图 （b）Figure 4. Histogram of strata and schematic diagram of observation （a） as well as the geothermal gradients map （b） of the Dangxiao well2. 基本动态特征分析
图5显示，2007—2015年党校井水温时值动态曲线以2011年末为界，可划分为两个阶段（图5a）：① 2007—2011年出现多次上升或下降的台阶式变化；② 2012—2015年期间表现为缓慢的上升—下降变化。与此相对，水位变化也可划分为两个阶段（图5b）：① 2007—2011年变化规律清晰，一般情况下每年4—6月水位下降直至低于泄流口不自流，8月之后水位逐渐升高，直至高于泄流口则开始自流至下一年度断流，年变化幅度不大，一般约为4 m；② 2012年开始水位低于分界线不再自流，保持静水位状态，水位动态表现为大幅度升降变化。
![]()
图 5 党校井水温（a）和水位（b）时值图Figure 5. Hourly values of water level （a） and temperature （b） for the Dangxiao well对水温水位动态曲线进行局部放大，可以看出，在2007—2011年期间，水温阶梯状的上升和下降变化与水位自流、非自流状态的转换具有明显的相关性（图6）。
当井水位由自流状态转换至非自流状态时，水温动态变化表现出以下特征（图6a）：① 即刻出现明显上升，例如2008年5月4—7日，一旦水位的潮汐波动低于泄流口，水温立即出现上升，使得水位由自流向非自流时断时续地过渡的情况下，水温出现大幅震荡；② 从自流到稳定的非自流状态一般要经历3—7日，在此前后温度上升幅度约为0.004—0.007℃；③ 自流状态下水温波动幅度大，转变为非自流状态后，水温动态更为平稳。
当井水位由静水位转换至自流状态时，水温动态变化表现出以下特征（图6b）：① 出现阶梯式下降；② 除2008年8月水位上升和下降均缓慢以外，其余3次当水位状态呈自流后，温度自下降至稳定状态只需7—16个小时，水位自流前后温度下降幅度约为0.008—0.011℃。
![]()
图 6 党校井由自流转换为断流 （a）和断流转换为自流 （b）状态时的水位水温时值图Figure 6. The hourly values of water level and temperature for the Dangxiao well at the state from artesian to non-artesian （a） and from non-artesian to artesian （b）水温与水位动态之间良好的实时相关性显示，温度探头位于井孔中水体与含水层连通的主要部位。同时，根据水温与水位的反向动态变化对比分析也可推断出，来自补给区的水源温度较低，这可以通过不同时期的两次井水温梯度测试结果 （图4b）进行验证。
第一次井水温梯度测试期间，井水位处于自流状态，井水温度自上而下变化幅度小 （图4b），50 m处温度值最低，190 m处温度值最高，50—190 m区段内水温缓慢增加，190—300 m区段内温度随深度增加略微降低，越接近温度探头放置处水温越低。
第二次测试期间，井水位处于非自流状态，井水温度自上而下变化幅度远大于自流状态下的观测值（图4b），40 m处温度值最低，226 m处温度值最高，40—196 m区段内温度随深度增加较快，196—226 m区段内增温缓慢，226—250 m区段内表现出微小的温度降低。
水温基本动态特征主要受水位泄流状态变化的影响。自流期间，含水层中较低温度的水在水头压力作用下不断向井孔内部流动上升，并通过泄流口溢出排泄，水温梯度曲线相对平缓，当含水层水头压力逐渐降低，较低温度的水向井孔中流动，通过泄流口排泄的水量也逐渐减少，直至转变为接近断流时，探头放置处的水温缓慢上升 （图5a），断流后水温迅速上升（图5a，6a）。断流状态下，当含水层水头压力逐渐增加，水位不断上升至超过泄流口时，较低温度的水加速向井孔中流动溢出，探头所在区域的水温很快出现下降（图6b）。
3. 同震变化特征分析
水温水位同震变化能够直接地、有效地揭示地下流体对应力-应变的响应，国内外近年来针对此问题已有较多的研究（石耀霖等，2007；杨竹转等，2008；Wang，Manga，2010；Cox et al，2015 ；Shi et al，2015 ）。本文使用分钟值观测资料来分析水温同震变化特征。
表1列出了2007—2015年引起党校井显著水温同震变化的8次强震及其相关参数，相应的震中分布图如图7所示。表中地震能量密度e（r）与震中距r和震级M之间的关系式为（Wang，Manga，2010）
$e\left( r \right) = {10^{1.45M - 3.03\lg r - 4.24}}.$
(1) 表 1 引起党校井显著水温同震变化的地震及其相关参数Table 1. The earthquakes caused obvious water temperature coseismic variations in the Dangxiao well and the corresponding parameters序号 发震时刻
东经/° 纬度/° MW 发震地点 震中距/km 地震能量
密度/（J·m–3）泄流
状态水位振荡
幅值/m水温同震
变化类型水温变化
幅值/10−4℃年-月-日 时：分 1 2007−09−12 19:11 101.40 4.46 （S） 8.5 印尼苏门答腊 3 487 0.002 2 自流 0.208 上升—下降 47 2 2007−09−13 07:49 100.73 2.62 （S） 7.9 印尼苏门答腊 3 281 0.000 4 自流 0.166 上升—下降 35 3 2008−05−12 14:28 103.37 31.06 （N） 7.9 四川汶川 556 0.079 0 静水位 0.119 下降—上升 −30 4 2011−03−11 13:46 142.50 38.30 （N） 9.1 日本 4 119 0.010 1 自流 0.711 上升—下降 48 5 2012−04−11 16:38 93.01 2.24 （S） 8.6 印尼苏门答腊 2 846 0.005 8 静水位 0.228 下降—上升 −41 6 2012−04−11 18:43 92.43 0.77 （S） 8.2 印尼苏门答腊 3 021 0.001 3 静水位 0.156 下降—上升 −30 7 2015−04−25 14:11 84.79 28.28 （N） 7.9 尼泊尔 1 529 0.003 7 静水位 0.163 下降—上升 −74 8 2015−12−07 15:50 72.78 38.21 （N） 7.2 塔吉克斯坦 2 850 0.000 1 静水位 0.016 下降—上升 −22 注：水位震荡幅值指最高值与最低值的差值；水温同震变化幅值指地震前的温度观测值与震后的水温上升或者下降最大值的差值；地震相关参数来自于ISC （International Seismological Centre，2017）。 ![]()
党校井水位水温同震变化分钟值曲线图（图8）显示，8次强震引起的水温同震变化形态可以分为两种类型：1，2，4号这3次地震（表1）为上升—下降—恢复型，其余5次地震为下降—上升—恢复型。两种类型的变化经过几十分钟至几小时后均能恢复到变化前的观测值附近。8次强震引起的水位同震变化类型均为振荡型，振荡很快达到峰值，峰值过后又快速衰减，振荡持续时间通常为几分钟至几十分钟，之后水位振荡逐渐减小直至恢复到变化前的观测值附近。中强震引起的水位振荡一般认为是在长周期瑞雷波的作用下，含水层膨胀和压缩引起孔隙压力的变化，导致含水层中的水进出井孔而造成井水面波动（Wang，Manga，2010）。
在井水位同震变化均为振荡的情况下，水温的同震变化表现出两种不同类型的变化形态。从水位观测值上看，1，2，4号这3次地震发生时井水位的观测值均大于0，显示水位为自流状态；其余5次地震发生时井水位的观测值均小于0，显示水位为静水位状态，表明水温同震变化明显受到泄流状态的影响，与震中的方位和震源机制解没有明确的相关性（图7）。
比较表1中水位振荡幅值与地震能量密度，可以观察到：自流状态下，两者之间基本呈正相关关系；而在静水位状态下，除汶川地震外，对于其余的4次地震而言，两者之间也保持了正相关关系。实际上，汶川地震是8次地震中地震能量密度最大的地震，但由于在地震前4—5天刚刚断流，水位处于泄流口以下0.05 m （图8），地震时振荡产生的水震波在向上运动时超过泄流口的部分通过泄流口流失，受此影响，振荡幅度偏小。即无论是在自流或非自流状态下，水位振荡幅值与地震能量密度均大致呈正相关关系。
比较表1中水温变化幅值与地震能量密度，可以看出，两者之间无明显的相关性，水温变化幅值的影响因素更为复杂。
自流期间，当地震波引起含水层水头压力波动时，加速了探头所在下部区域较低温度的水体与井孔上部较高温水体的对流混合，较高温的水混入到低温的水体中，使得探头所在处的温度值短暂相对升高，水位震荡平息后水温迅速恢复。非自流期间，井孔内部水体与含水层中水体交换微弱，当地震波引起含水层水位震荡时，含水层中较低温度的水体进出井孔中，混水作用引起了探头附近温度值的迅速下降，震荡平息后很快恢复。因此放置在含水层与井孔水体连通主要部位的温度探头能够迅速地反映以上温度的同震变化特征。
水温同震变化特征反映出的温度探头放置位置和低温水源补给与基本动态分析中获得的认识相一致。
![]()
图 8 党校井水温水位同震变化分钟值曲线图数字与表1中的8次地震相对应，虚线表示发震时间Figure 8. Coseismic variations of minute values of water temperature and water level for the Dangxiao wellDashed lines denote the earthquakes occurrence time；the numbers are consistent with the serial numbers in Table 14. 前兆变化分析
4.1 水温前兆异常机理
关于水温前兆异常机理，一般认为，由于岩石热传导率低，震源体附近的热量或者断层表面的摩擦热传导到地面需要很长时间，引起温度变化的异常量也非常小，所以要在岩石中观测到温度前兆异常是极为困难的；但是深部产生的热通过地下水的流动可以在较短时间传输到地表，特别是在导水导气的断裂带中，由于存在热对流，很有可能在有限的时间内影响到地壳表层而且变化量相对较大（Cicerone et al，2009 ）。
与中强震孕育有关的应力变化，一方面可引起地下水流动状态或者流向的变化，进而导致水温的变化（车用太等，2008）；另一方面，又可导致岩体孔隙结构的变化，产生新裂隙，例如震前岩体膨胀产生的新裂隙和原有裂隙的加宽使得流体活动增强（刘耀炜等，2015），使受到限制的地下水循环到更深的、更热的地方。深部岩体的热量通过地下水和气体的流动传导上来引起近地表水体的升温变化，若裂隙的改变堵塞了高温水与地表水的连通，则会引起降温的异常变化（Cicerone et al，2009 ）。
4.2 降雨对党校井水位的影响
考虑到党校井与黑龙潭泉群具有相同的区域水文地质背景，我们将1992年以来党校井的水位日均值（图9a）与黑龙潭泉群的泄流情况（图9b）进行对比，可以看出党校井与黑龙潭泉群的断流时段有很好的对应性。多年的观测值显示，党校井水位下降至泄流口以下0.6—0.7 m时黑龙潭泉群口会出现断流，两者良好的相关性也进一步证实它们具有相同的补给来源。
图9c为1992年以来丽江的年降雨量，在1992，1994，1997，2003，2005，2006和2009年降雨量均小于920 mm的情况下，黑龙潭泉群在次年均发生了断流（康晓波等，2013），党校井水位除2004年外也均出现了断流的情况。2011—2015年在年降雨量连续小于920 mm的情况下，党校井水位自2011年11月20日起断流，至今仍为静水位，黑龙潭泉群从2012年1月21日起断流。降雨量最少的年份为2012年，只有634 mm，党校井水位于次年，即2013年7月12日下降至17.738 m的历史最低值。上述分析显示，党校井水位泄流状态的变化主要受大气降雨的补给影响，降雨对党校井水位的影响在时间上有滞后性，表现为趋势性的影响。
![]()
图 9 党校井水位 （a）、黑龙潭泉群泄流状态 （b）及年降雨量统计直方图 （c）Figure 9. Water level of the Dangxiao well （a），the discharging state of Heilongtan spring groups （b） and the annual rainfall histogram （c）4.3 党校井水温前兆变化分析
据文献资料显示，黑龙潭泉群多年平均泄流量为1.23 m3/s （康晓波等，2013），党校井的泄流量为1.181×10–3 m3/s，前者约为后者的1 000倍。泄流期间，水头压力变化的大部分信息通过流量的变化流失掉了。
2007—2011年期间，党校井水位每年出现两次自流与非自流的转换，相应的水温每年出现两次上升或下降的台阶式变化（图5）。在此期间，受地下水体剧烈活动的影响，水温动态难以反映与地震孕育相关的变化信息。
2012—2015年期间，党校井一直处于断流状态，黑龙潭泉群除了在2014年9月下旬至2015年1月当井水位的观测值大于−0.6 m时出现泄流外，其余时间也处于断流状态。这一时期内，受年度降雨的趋势性影响，党校井水位表现出夏低冬高的年变形态（图10），水位自年初缓慢下降，至7—8月达到最低值，之后逐渐上升，至11—12月达到年度最高值。水温则在趋势上升的背景上叠加夏高冬低的大致年变形态（图10），一般在8—9月达到最高值。与水位年度最低值时间相比，水温年度最高值一般情况下要滞后10—20天，显示含水层中水头压力上升，低温水源增多时探头附近水温缓慢降低。党校井水温2014年9月加速下降，而2015年1月快速上升，分析认为可能与黑龙潭泉群在这段时间的泄流和断流有关。总体来说，水温年动态变化受到水位动态规律的明显影响。对照周边区域发生的中强震（图1），从目前的分析结果来看，水温前兆变化特征不明显，其中是否含有前兆变化信息有待于资料的进一步积累和更多的分析研究。
![]()
图 10 党校井2012—2015年水温水位时值曲线图Figure 10. Hourly values of water level and temperature for the Dangxiao well from 2012 to 20155. 讨论与结论
本研究以丽江党校井为例，结合同井水位动态、区域水文地质、井孔结构和井温梯度，分析水温动态的基本特征，并根据分钟值观测资料记录到的８次中强震的水温水位同震变化分析了水温的动态特征，最后对水温的前兆变化进行了讨论，得到如下结论和认识：
1） 水温水位对比和井温梯度结构分析显示，党校井温度探头放置在井孔水体与含水层连通的主要部位，补给水源的水温较低，水温基本动态特征主要受水位泄流状态和水位变化的影响；
2） 水位同震变化为振荡，水温同震变化形态在不同的泄流状态下存在明显差异，自流期间为上升—下降—恢复型，非自流期间为下降—上升—恢复型，同震变化分析反映出的温度探头放置位置和低温水源补给的认识与基本动态分析获得的认识相一致；
3） 前兆变化方面，党校井水位动态受降雨的趋势性影响，2007—2011年水位每年有两次泄流和非泄流的状态转换，地下水体活动剧烈，水温动态难以反映与地震孕育相关的变化信息，2012—2015年非自流状态期间，水温动态前兆变化特征不明显。
研究分析得到的结果和认识可以为未来水温动态分析和异常落实提供一定的依据和参考。
井水温度对地下水的运动非常敏感，因此水温动态分析首先需要结合同井水位的动态变化来研究，对比长时间水位水温的观测值，分析不同时期两者之间的变化特征差异，必要时，需要对水位变化进行深入研究，通过区域地质构造、水文地质资料、井孔结构、区域降雨资料等分析井孔所在含水层的补给和降雨情况；其次，结合井温梯度结构分析了补给水源的温度特性和井孔与含水层连通的主要补给部位；最后，水温水位的同震变化是地下流体对应力应变的直接反映，使得我们可以从另一个角度分析水温水位动态变化的相互关系。
党校井水温动态分析对温度探头的放置也具有一定的启示意义。关于温度探头放置位置的选择在 《地震地下流体观测方法：井水和泉水温度观测》 （DB/T 49—2012）中已有说明（中国地震局，2012），但每口井的情况仍需要具体对待。当水位补给来源主要为浅部降雨时，水温探头的放置应避开井孔与含水层连通的主要部位。
实际观测中，井孔通常在某一深度以下均与含水层相通，要寻找避开井孔水体与含水层连通的主要段落，本文研究显示，可以通过短期水位水温对比、温度梯度结构特性、水温水位同震变化记录等方面帮助寻找。
由于井孔位置、结构、水文地质情况的差异，每一口地震观测井的水温均有其独特的变化特征，本文的分析可能仍然带有很大的局限性，未来需要通过对更多的井孔加以研究来弥补。
云南省地震局付虹研究员、中国地震局地壳应力研究所刘耀炜研究员在稿件撰写过程中给予了建议、帮助和支持，审稿专家提出了有益的修改建议，作者在此一并表示衷心的感谢。
-
![]()
图 1 研究区周边断裂及台站分布
Figure 1. Geological structure and station distribution in the study area and surroundings
![]()
图 2 精定位前（未填充颜色）、后（填充颜色）的地震分布
Figure 2. Earthquake distribution before （uncolored） and after precisely locating （colored）
![]()
图 3 地震序列1段（区域2）的断层面拟合
（a） 地震在水平面投影；（b） 地震在断层面的投影；（c） 地震在垂直于断层的横断面上的投影；（d） 地震与断层面的距离
Figure 3. Fitting of fault surface in the first section （region 2） of earthquake sequence
（a） The earthquake in the horizontal projection；（b） The earthquake in the faultplane projection；（c） The projection on a cross section perpendicular to the fault；（d） The distance between the earthquake and the fault plane
![]()
图 4 地震序列2段（区域1)的断层面拟合
（a）地震在水平面投影；（b）地震在断层面的投影；（c）地震在垂直于断层的横断面上的投影；（d）地震与断层面的距离
Figure 4. Fitting of fault surface in the second section （region 1） of earthquake sequence
（a） The earthquake in the horizontal projection of the middle panel；（b） The earthquake in the fault plane projection；（c） The projection on a cross section perpendicular to the fault；（d） The distance between the earthquake and the fault plane
![]()
图 5 震源机制解的P轴（a）和T轴（b）统计分布及总体震源机制的等面积投影（c）
蓝色表示压缩，红色表示拉张
Figure 5. P-axis （a） and T-axis （b） statistical distribution of the focal mechanism solution and the equal area projection of the overall focal mechanism （c）
Blue means compression，red means expansion
![]()
图 6 本研究中区域1和区域2的应力场反演结果
图（a）为区域整体应力示意图；图（b）、（d）为应力张量反演的等面积投影图。蓝色箭头和红色箭头分别表示“可能断层面”的实际测量的滑动方向和理论滑动方向；绿色表示在95%的置信水平下应力场的最大剪应力节面；黄色箭头表示最大剪应力的滑动方向。S3，S2和S1及标注周边的闭合曲线分别表示主压应力轴、中间应力轴和主张应力轴的置信范围。U和D分别为上和下两个方向。图（c）、（e）为两个子区应力张量反演结果的三维表示，红色表示主压应力的大小和方向；蓝色表示主张应力的大小和方向。下同
Figure 6. Stress field inversion result in region 1 and region 2
Figs.（b） and （d） are equal-area projections of the stress tensor inversion. The blue and red arrows indicate the actual measured sliding direction and the theoretical sliding direction of the “possible fault plane”，respectively；the green indicates the maximum shear stress nodal plane of the stress field at 95% confidence level； the yellow arrow indicates the sliding direction of the maximum shear stress. The S3， S2， S1 and the closed curves marked around them respectively represent the confidence ranges of the principal compressive stress axis, the intermediate stress axis and the principal tensile stress axis.U and D are the top and bottom directions respectively. Figs.（c） and （e） are three-dimensional representations of the results of the inversion of the stress tensor in the two sub-areas，with the magnitude and direction of the principal compressive stress in red；the magnitude and direction of the asserted stress in blue. The same below
![]()
图 7 应力场方向随深度变化示意图
图中阴影表示板块交界面，黑色箭头表示构造应力场方向
Figure 7. Schematic diagram of stress field changing with depth
The shadow zone in the figure indicates the plate interface，and the black arrows indicate the compressional axis direction of tectonic stress
![]()
图 8 2013年2月10日M1.4地震 （a）和2014年3月14日M3.8 地震（b）的P波初动分布和震源机制的等面积投影
Figure 8. Equal area projection of the P-wave initial motion sign and estimated source mechanism for the February 10，2013，M1.4 earthquake （a） and March 14，2014，M3.8 earthquake （b）
表 1 福建地区地壳速度模型（张路，2008）
Table 1 Crustal velocity model in Fujian area
地壳厚度/km vp/（km·s−1） vS/（km·s−1） 密度/（kg·m−3） 0—4.0 5.6 3.23 2.5 4.0—12.0 6.1 3.52 2.7 12.0—18.0 6.0 3.46 2.6 18.0—25.0 6.4 3.70 2.8 25.0—30.0 6.9 3.98 2.9 ≥30.0 8.0 4.57 3.3 
下载: 导出CSV
表 2 应力张量反演数值结果
Table 2 Results of the inverted stress tensor
分区 震源
机制解
个数主压应力轴S1 中间应力轴S2 主张应力轴S3 R 走向/° 倾伏角/° 走向/° 倾伏角/° 走向/° 倾伏角/° 区域1 138 182.38°—188.46 10.16°—24.28° 329.00°— 335.00° 61.00°—78.00° 89.58°—95.58° −0.07°—16.93° 0.61 区域2 20 337.02°—344.43° 1.31°—7.04° 107.50°— 115.00° 79.50°—88.00° 246.07°—253.57° 7.03°—15.53° 0.41 深度8—9 km 31 178.62°—197.62° 13.57°— 39.02° 328.00°—347.00° 46.00°—75.00° 78.18°—97.18° 9.08°—38.08° 0.49 深度10—11 km 80 174.89°—181.25° 19.29°—26.97° 339.00°—343.00° 62.00°—70.00° 83.04°—87.04° −1.68°—6.32° 0.45 深度12—15 km 30 334.41°—346.41° 8.75°—15.41° 228.00°—240.00° 39.00°—60.00° 72.22°—84.22° 18.53°—39.53° 0.81 深度16—21 km 17 344.03°—347.28° 1.31°—4.80° 99.50°—103.00° 79.00°—87.00° 253.22°—256.72° 6.08°—14.08° 0.53 前人研究 50 −0.44°—3.56° 12.42°—16.24° 162.00°—166.00° 73.00°—77.00° 268.45°—272.45° 2.34°—6.34° 0.32
下载: 导出CSV
1 P波初动求得的震源机制解
序号 发震时间
年−月−日 时:分:秒东经
/°北纬
/°深度
/km震级
M节面Ⅰ 节面Ⅱ P 轴 T 轴 B 轴 矛盾比 P波初动个数 走向/° 倾角/° 滑动角/° 走向/° 倾角/° 滑动角/° 方位角/° 倾伏角/° 方位角/° 倾伏角/° 方位角/° 倾伏角/° 1 2 014−07−19 00:47:06 118.742 25.633 10 1.4 330 81 −177 240 87 −9 195 9 285 4 40 80 0 10 2 2 014−04−06 10:12:33 118.743 25.633 10 1.5 174 53 −115 32 44 −60 24 69 282 5 190 2 0 0.09 11 3 2 014−01−12 03:55:36 118.747 25.633 9 1.7 136 64 164 233 76 27 3 7 97 29 260 60 0.07 14 4 2 014−03−14 20:25:11 118.743 25.635 10 1.7 134 67 −160 36 71 −24 354 30 86 3 180 60 0 14 5 2 014−03−14 19:53:36 118.744 25.635 11 3.8 140 90 170 230 80 0 185 7 95 7 320 80 0 36 6 2 014−07−19 01:00:22 118.739 25.636 9 1.7 315 81 −120 211 31 −17 194 45 69 30 320 30 0.06 18 7 2 014−07−19 00:45:27 118.742 25.636 10 1.6 322 44 120 104 53 65 212 5 314 69 120 2 0 0.07 14 8 2 014−03−14 19:52:20 118.742 25.636 9 1.3 340 82 174 70 84 8 2 05 1 295 10 110 80 0 10 9 2 014−01−17 23:35:43 118.741 25.636 9 1.5 23 31 −109 226 61 −78 163 72 307 15 40 10 0.09 11 10 2 014−01−07 06:49:43 118.742 25.636 10 1.7 146 64 164 243 76 27 13 7 107 29 270 60 0.20 10 11 2 014−01−07 06:47:17 118.744 25.636 10 1.4 174 64 −124 51 41 −41 38 57 288 13 190 30 0.09 11 12 2 014−07−23 23:52:37 118.741 25.637 9 1.8 334 22 −117 184 70 −79 110 63 265 25 0 10 0 12 13 2 014−05−12 18:28:38 118.740 25.637 10 1.2 320 80 178 50 88 10 185 6 275 8 60 80 0.09 11 14 2 014−03−14 20:05:29 118.743 25.637 10 1.7 310 80 180 220 90 170 175 7 265 7 40 80 0 11 15 2 014−03−12 14:50:09 118.741 25.637 9 1.8 330 81 −177 240 87 −9 195 9 285 4 40 80 0 11 16 2 014−02−27 14:32:29 118.740 25.637 10 2.5 310 80 180 220 90 170 175 7 265 7 40 80 0 2 0 17 2 014−08−01 22:43:04 118.740 25.638 9 2.3 308 75 167 42 77 16 175 2 265 2 0 80 70 0 12 18 2 014−07−25 06:44:18 118.741 25.638 10 1.8 334 61 −102 177 31 −71 217 72 73 15 340 10 0.06 16 19 2 014−05−24 00:18:51 118.738 25.638 10 1.2 160 10 −90 340 80 −90 250 55 70 35 340 0 0 10 2 0 2 014−05−12 18:24:19 118.739 25.638 10 3.0 325 81 −120 221 31 −17 2 04 45 79 30 330 30 0.07 30 21 2 014−03−12 05:37:06 118.744 25.638 9 1.7 132 44 120 274 53 65 22 5 124 69 290 2 0 0.09 11 22 2 014−03−12 04:24:19 118.743 25.638 9 1.6 160 88 170 250 80 2 2 05 6 115 8 330 80 0 12 23 2 014−01−07 06:43:42 118.740 25.638 10 2.9 156 64 164 253 76 27 23 7 117 29 280 60 0.07 27 24 2 014−01−02 13:56:38 118.740 25.638 10 2.5 150 44 158 257 75 48 17 19 126 44 270 40 0 18 25 2 014−07−19 01:50:15 118.739 25.639 9 3.3 150 2 0 180 60 90 110 349 42 131 42 240 2 0 0.03 31 26 2 014−07−03 17:14:26 118.738 25.639 11 2.9 320 87 171 50 81 3 5 4 275 9 120 80 0 26 27 2 014−07−03 04:56:06 118.738 25.639 13 2.1 316 67 160 54 71 24 184 3 276 30 90 60 0.11 19 28 2 014−06−29 03:41:30 118.739 25.639 10 1.9 130 90 170 220 80 0 175 7 85 7 310 80 0 13 29 2 014−06−29 02:17:36 118.739 25.639 10 2.0 152 61 −138 37 54 −37 8 50 273 4 180 40 0.04 23 30 2 014−05−24 00:14:12 118.737 25.639 11 2.8 140 90 110 230 2 0 0 211 42 69 42 320 2 0 0.03 29 31 2 014−07−19 00:32:21 118.738 25.640 10 2.3 312 48 −149 2 00 67 −46 157 48 260 11 0 40 0.07 27 32 2 014−05−24 00:34:06 118.738 25.640 10 1.3 309 62 −113 171 36 −54 178 65 55 14 320 2 0 0.09 11 33 2 014−02−10 23:29:21 118.739 25.640 10 1.8 138 61 174 231 85 30 1 17 99 24 240 60 0.10 2 0 34 2 014−02−02 00:03:47 118.737 25.640 10 2.0 316 64 164 53 76 27 183 7 277 29 80 60 0.04 28 35 2 014−01−27 13:50:18 118.737 25.640 12 3.1 130 90 170 220 80 0 175 7 85 7 310 80 0 30 36 2 014−01−27 01:02:39 118.737 25.640 10 1.9 320 82 −174 230 84 −8 185 10 275 1 10 80 0 16 37 2 014−01−02 23:34:08 118.738 25.640 10 2.3 321 83 −161 229 71 −7 186 18 94 8 340 70 0.04 27 38 2 014−10−10 04:39:32 118.738 25.641 8 1.7 292 61 −138 177 54 −37 148 50 53 4 320 40 0 15 39 2 014−08−13 13:29:12 118.735 25.641 10 1.9 130 90 100 220 10 0 210 44 50 44 310 10 0.07 15 40 2 014−07−27 20:52:40 118.736 25.641 9 2.3 130 87 −171 40 81 −3 355 9 265 4 150 80 0.04 23 41 2 014−07−05 22:15:24 118.738 25.641 12 1.9 140 73 −122 24 36 −31 13 52 253 21 150 30 0.10 2 0 42 2 014−05−25 05:04:11 118.736 25.641 10 1.8 150 90 100 240 10 0 230 44 70 44 330 10 0.11 18 43 2 014−05−24 00:36:56 118.736 25.641 10 1.9 130 90 100 220 10 0 210 44 50 44 310 10 0 15 44 2 014−05−02 20:13:57 118.735 25.641 10 1.9 130 90 100 220 10 0 210 44 50 44 310 10 0 17 45 2 014−01−27 13:53:04 118.738 25.641 10 2.5 140 90 100 230 10 0 220 44 60 44 320 10 0 18 46 2 014−01−21 00:53:13 118.740 25.641 10 1.9 160 90 170 250 80 0 2 05 7 115 7 340 80 0.08 13 47 2 014−01−17 23:07:55 118.738 25.641 9 1.6 137 62 169 233 80 28 2 12 98 27 250 60 0.14 14 48 2 014−12−02 02:13:14 118.735 25.642 8 2.0 140 90 0 230 90 180 5 0 95 0 0 90 0.25 16 49 2 014−07−27 21:02:03 118.735 25.642 9 1.5 150 87 171 240 81 3 195 4 105 9 310 80 0.09 11 50 2 014−07−27 20:28:59 118.734 25.642 9 2.6 324 82 −130 225 41 −12 197 39 84 26 330 40 0.07 29 51 2 014−07−07 23:44:34 118.734 25.642 11 1.5 339 70 176 71 87 20 2 03 11 297 16 80 70 0.20 10 52 2 014−05−29 05:35:08 118.737 25.642 10 1.9 151 83 −161 59 71 −7 16 18 284 8 170 70 0.14 14 53 2 014−02−01 15:43:52 118.735 25.642 10 1.9 130 90 100 220 10 0 210 44 50 44 310 10 0 11 54 2 014−01−31 13:31:55 118.736 25.642 10 1.5 150 90 110 240 2 0 0 221 42 79 42 330 2 0 0 11 55 2 014−01−17 22:32:14 118.738 25.642 9 1.3 314 64 −124 191 41 −41 178 57 68 13 330 30 0 11 56 2 014−10−12 06:20:45 118.737 25.643 9 2.1 120 90 170 210 80 0 165 7 75 7 300 80 0 21 57 2 014−10−10 04:37:57 118.738 25.643 8 1.0 111 31 163 215 81 60 329 30 94 45 220 30 0.10 10 58 2 014−10−08 12:52:44 118.736 25.643 10 1.8 328 41 −105 168 51 −77 133 79 249 5 340 10 0 17 59 2 014−10−08 04:32:36 118.736 25.643 9 1.6 70 10 −90 250 80 −90 160 55 340 35 250 0 0 10 60 2 014−07−27 21:44:02 118.736 25.643 8 1.6 130 90 0 220 90 180 355 0 85 0 0 90 0 12 61 2 014−05−27 10:13:04 118.735 25.643 10 1.7 130 90 100 220 10 0 210 44 50 44 310 10 0 13 62 2 014−02−05 04:06:32 118.733 25.643 10 1.7 70 2 0 −90 250 70 −90 160 65 340 25 250 0 0 12 63 2 014−12−13 19:58:19 118.736 25.644 9 1.4 89 85 150 182 61 6 139 17 41 24 260 60 0.09 11 64 2 014−12−05 12:35:40 118.737 25.644 8 2.5 344 61 −102 187 31 −71 227 72 83 15 350 10 0.29 14 65 2 014−11−10 06:54:39 118.735 25.644 9 1.5 134 22 −117 344 70 −79 270 63 65 25 160 10 0 10 66 2 014−11−10 06:05:07 118.734 25.644 9 1.3 299 62 −113 161 36 −54 168 65 45 14 310 2 0 0 11 67 2 014−10−15 01:22:03 118.737 25.644 9 2.0 300 80 −178 210 88 −10 165 8 255 6 20 80 0 19 68 2 014−07−12 06:45:45 118.733 25.644 12 2.1 330 30 180 240 90 120 177 38 303 38 60 30 0 2 0 69 2 014−07−08 22:08:55 118.733 25.644 11 1.5 137 81 −110 23 22 −26 24 50 243 33 140 2 0 0.10 10 70 2 014−11−26 07:00:14 118.732 25.645 10 1.3 287 22 −154 173 81 −70 106 50 247 33 350 2 0 0 11 71 2 014−11−18 22:30:16 118.732 25.646 10 1.7 130 90 130 220 40 0 187 33 73 33 310 40 0 17 72 2 014−09−22 14:44:57 118.733 25.647 10 1.9 320 81 175 50 85 9 185 3 275 10 80 80 0 13 73 2 013−09−04 06:53:04 118.761 25.616 11 2.4 330 80 180 240 90 170 195 7 285 7 60 80 0 19 74 2 013−09−05 04:21:47 118.758 25.617 11 1.8 153 71 −111 23 28 −43 34 58 259 23 160 2 0 0.11 18 75 2 013−11−06 21:31:24 118.756 25.619 11 2.8 137 76 153 234 64 16 187 7 93 29 290 60 0 26 76 2 013−11−04 20:34:40 118.755 25.619 11 3.0 158 73 170 252 80 18 24 5 116 19 280 70 0.12 26 77 2 013−09−06 01:13:03 118.757 25.619 10 1.6 316 71 −144 213 56 −23 179 38 82 10 340 50 0.08 12 78 2 013−09−14 03:13:50 118.756 25.620 10 1.4 318 80 163 52 73 10 6 5 274 19 110 70 0 10 79 2 013−08−30 05:22:13 118.756 25.620 11 1.9 308 80 −100 175 14 −45 2 06 54 47 34 310 10 0.10 2 0 80 2 013−11−04 20:34:28 118.755 25.621 11 1.4 340 60 −90 160 30 −90 250 75 70 15 160 0 0.19 16 81 2 013−10−31 20:50:19 118.756 25.621 11 2.2 159 48 −132 33 56 −53 0 60 98 4 190 30 0 2 0 82 2 013−09−27 11:45:15 118.756 25.621 11 2.1 318 80 163 52 73 10 6 5 274 19 110 70 0 12 83 2 013−09−14 02:59:50 118.756 25.622 10 3.4 133 75 132 240 44 22 193 19 84 44 300 40 0.03 38 84 2 013−08−25 22:10:29 118.756 25.622 11 1.8 149 83 161 241 71 7 196 8 104 18 310 70 0 13 85 2 013−11−19 03:03:58 118.755 25.623 10 3.2 133 75 132 240 44 22 193 19 84 44 300 40 0 26 86 2 013−09−13 14:22:14 118.752 25.623 11 2.1 150 88 170 240 80 2 195 6 105 8 320 80 0.12 17 87 2 013−10−03 16:49:27 118.757 25.624 10 2.3 320 73 −122 2 04 36 −31 193 52 73 21 330 30 0.05 19 88 2 013−09−13 15:12:20 118.749 25.624 11 2.0 329 70 176 61 87 20 193 11 287 16 70 70 0.09 11 89 2 013−09−08 22:56:38 118.750 25.627 9 2.5 311 87 −160 219 70 −4 177 16 83 11 320 70 0.15 26 90 2 013−08−24 00:49:57 118.752 25.627 11 3.1 130 90 110 220 2 0 0 2 01 42 59 42 310 2 0 0.06 31 91 2 013−08−09 13:37:12 118.750 25.628 11 3.1 310 80 178 40 88 10 175 6 265 8 50 80 0.10 10 92 2 013−08−09 13:53:24 118.749 25.629 12 2.9 132 84 −140 37 51 −8 2 32 258 22 140 50 0.04 25 93 2 013−08−07 05:30:32 118.746 25.629 11 1.8 331 85 −150 238 61 −6 199 24 101 17 340 60 0.09 11 94 2 013−08−03 03:38:44 118.748 25.629 12 2.0 327 56 −127 2 01 48 −48 180 60 83 4 350 30 0.06 16 95 2 013−08−19 17:36:20 118.748 25.630 10 3.8 140 90 170 230 80 0 185 7 95 7 320 80 0 34 96 2 013−08−09 13:37:07 118.749 25.630 12 2.5 320 90 170 50 80 0 5 7 275 7 140 80 0.05 21 97 2 013−08−02 21:10:17 118.747 25.630 10 1.8 139 83 161 231 71 7 186 8 94 18 300 70 0.09 11 98 2 013−01−20 23:33:58 118.744 25.630 10 1.7 330 81 −177 240 87 −9 195 9 285 4 40 80 0 11 99 2 013−10−18 14:05:23 118.744 25.631 10 3.3 140 90 170 230 80 0 185 7 95 7 320 80 0.03 35 100 2 013−08−09 13:38:40 118.748 25.631 12 3.5 137 75 −132 30 44 −22 6 44 257 19 150 40 0.04 23 101 2 013−08−03 05:23:36 118.749 25.631 9 2.8 324 77 −142 225 53 −16 191 36 90 16 340 50 0.04 24 102 2 013−02−11 01:05:50 118.740 25.633 12 2.4 139 87 160 231 70 4 187 11 93 16 310 70 0.03 29 103 2 013−02−10 03:25:46 118.742 25.633 14 2.5 320 84 −172 230 82 −6 185 10 95 1 0 80 0.03 32 104 2 013−12−09 07:21:06 118.746 25.634 10 2.0 320 80 180 230 90 170 185 7 275 7 50 80 0 17 105 2 013−08−19 18:05:28 118.744 25.634 10 2.7 140 90 0 230 90 180 5 0 95 0 0 90 0.12 32 106 2 013−08−03 02:43:55 118.750 25.634 9 2.9 320 85 171 50 81 5 5 3 275 10 110 80 0.06 17 107 2 013−12−09 23:25:13 118.746 25.635 10 1.8 316 82 130 55 41 12 16 26 263 39 130 40 0 17 108 2 013−11−03 00:13:46 118.749 25.635 9 2.3 127 76 153 224 64 16 177 7 83 29 280 60 0 24 109 2 013−10−17 03:29:21 118.745 25.635 9 1.7 320 80 −178 230 88 −10 185 8 275 6 40 80 0 10 110 2 013−08−03 02:43:56 118.746 25.635 10 4.2 320 80 180 230 90 170 185 7 275 7 50 80 0 29 111 2 013−03−23 02:39:47 118.742 25.635 10 1.6 320 80 178 50 88 10 185 6 275 8 60 80 0.08 13 112 2 013−02−10 05:59:11 118.741 25.635 14 1.9 320 80 −178 230 88 −10 185 8 275 6 40 80 0 21 113 2 013−02−10 04:56:09 118.740 25.635 12 1.4 140 90 170 230 80 0 185 7 95 7 320 80 0.05 21 114 2 013−02−10 03:29:35 118.744 25.635 10 1.0 306 64 164 43 76 27 173 7 267 29 70 60 0.18 11 115 2 013−01−20 18:59:31 118.739 25.636 13 1.9 330 80 −178 240 88 −10 195 8 285 6 50 80 0 11 116 2 013−12−17 04:56:45 118.740 25.637 10 2.1 320 84 −172 230 82 −6 185 10 95 1 0 80 0 16 117 2 013−02−23 23:23:14 118.742 25.637 13 3.0 320 80 −178 230 88 −10 185 8 275 6 40 80 0 30 118 2 013−01−19 22:13:09 118.739 25.637 10 1.7 310 80 180 220 90 170 175 7 265 7 40 80 0 12 119 2 013−01−12 22:57:21 118.740 25.637 10 2.4 310 80 180 220 90 170 175 7 265 7 40 80 0 22 120 2 013−12−30 08:13:46 118.740 25.638 10 2.5 320 80 −178 230 88 −10 185 8 275 6 40 80 0 19 121 2 013−02−24 20:39:22 118.738 25.638 11 1.9 334 71 −156 236 67 −20 196 30 104 3 10 60 0 16 122 2 013−02−24 07:25:52 118.736 25.638 11 1.6 320 80 −178 230 88 −10 185 8 275 6 40 80 0 12 123 2 013−01−27 10:00:27 118.736 25.638 12 1.8 130 90 110 220 2 0 0 2 01 42 59 42 310 2 0 0.09 11 124 2 013−01−11 18:16:55 118.742 25.638 11 2.0 159 87 160 251 70 4 2 07 11 113 16 330 70 0.12 16 125 2 013−12−16 19:01:36 118.739 25.639 10 3.1 140 90 0 230 90 180 5 0 95 0 0 90 0.09 32 126 2 013−12−20 23:40:48 118.739 25.640 10 1.8 320 80 180 230 90 170 185 7 275 7 50 80 0 11 127 2 013−09−04 06:23:26 118.750 25.640 10 4.8 141 83 −161 49 71 −7 6 18 274 8 160 70 0.02 80 128 2 013−05−28 12:30:01 118.738 25.641 12 1.8 310 80 180 220 90 170 175 7 265 7 40 80 0 12 129 2 013−05−28 13:06:38 118.740 25.642 13 1.7 320 80 −178 230 88 −10 185 8 275 6 40 80 0 10 130 2 012−11−30 14:59:13 118.749 25.627 18 3.3 330 88 170 60 80 2 15 6 285 8 140 80 0.10 21 131 2 012−11−25 08:44:31 118.747 25.627 17 2.6 320 85 171 50 81 5 5 3 275 10 110 80 0 13 132 2 012−11−25 17:42:39 118.749 25.629 17 3.4 144 77 −142 45 53 −16 11 36 270 16 160 50 0.09 22 133 2 012−11−25 08:04:03 118.748 25.629 15 2.2 320 90 170 50 80 0 5 7 275 7 140 80 0 12 134 2 012−11−25 07:48:49 118.748 25.631 17 3.8 310 80 180 220 90 170 175 7 265 7 40 80 0 22 135 2 012−11−15 23:44:32 118.747 25.631 19 2.4 329 48 −132 2 03 56 −53 170 60 268 4 0 30 0 16 136 2 012−11−25 00:10:52 118.750 25.632 16 1.9 329 36 −126 191 62 −67 142 65 265 14 0 2 0 0 11 137 2 012−11−12 12:26:56 118.747 25.632 15 2.0 137 71 111 267 28 43 211 23 76 58 310 2 0 0 11 138 2 012−11−25 07:55:45 118.744 25.634 13 2.3 142 80 100 275 14 45 223 34 64 54 320 10 0.08 13 139 2 012−01−21 08:11:14 118.726 25.649 12 2.0 247 54 143 2 61 42 123 4 218 50 30 40 0.10 10 140 2 012−01−19 18:45:47 118.730 25.651 17 2.7 290 85 171 20 81 5 335 3 245 10 80 80 0 16 141 2 012−05−24 18:02:17 118.724 25.653 16 2.8 302 61 152 47 66 33 174 3 267 40 80 50 0.08 13 142 2 012−05−11 22:39:22 118.721 25.653 18 3.2 140 82 174 230 84 8 5 1 95 10 270 80 0.04 23 143 2 012−04−15 12:02:02 118.717 25.654 14 3.4 140 80 −178 50 88 −10 5 8 95 6 220 80 0 21 144 2 012−01−19 05:51:18 118.727 25.654 16 2.2 295 53 164 34 77 38 160 16 261 36 50 50 0.07 14 145 2 012−01−17 22:51:27 118.723 25.654 14 2.4 310 88 170 40 80 2 355 6 265 8 120 80 0.08 13 146 2 012−06−10 16:17:54 118.717 25.657 17 2.6 149 41 −139 26 64 −56 342 57 92 13 190 30 0 13 147 2 012−04−02 05:56:39 118.712 25.659 13 2.5 320 85 171 50 81 5 5 3 275 10 110 80 0 13 148 2 012−03−22 05:10:09 118.717 25.659 12 2.5 120 80 −180 30 90 −10 345 7 75 7 210 80 0 14 149 2 012−07−22 02:27:33 118.709 25.660 14 2.0 280 10 180 190 90 100 110 44 270 44 10 10 0.09 11 150 2 012−04−15 11:57:01 118.708 25.660 14 4.1 147 81 −110 33 22 −26 34 50 253 33 150 2 0 0.06 32 151 2 012−04−04 02:15:41 118.716 25.660 13 2.6 140 80 −178 50 88 −10 5 8 95 6 220 80 0 11 152 2 012−02−02 23:06:17 118.715 25.660 15 2.1 300 90 90 165 0 135 30 45 210 45 300 0 0.15 13 153 2 012−04−22 10:22:46 118.719 25.661 18 2.5 290 30 180 2 00 90 120 137 38 263 38 20 30 0 14 154 2 012−04−01 18:39:46 118.712 25.661 16 2.5 130 87 −171 40 81 −3 355 9 265 4 150 80 0.09 11 155 2 012−03−30 03:34:43 118.713 25.661 18 2.4 302 84 −140 2 07 51 −8 172 32 68 22 310 50 0.07 14 156 2 012−04−18 13:11:33 118.717 25.662 19 2.3 310 81 175 40 85 9 175 3 265 10 70 80 0 13 157 2 012−04−02 04:44:40 118.713 25.662 18 2.0 129 62 −113 351 36 −54 358 65 235 14 140 2 0 0 11 158 2 011−12−30 02:38:38 118.722 25.657 21 2.9 130 80 −178 40 88 −10 355 8 85 6 210 80 0 15 159 2 010−08−06 11:46:26 118.721 25.623 13 1.7 121 62 113 259 36 54 195 14 72 65 290 2 0 0 10
下载: 导出CSV
-
陈翰林,赵翠萍,修济刚,陈章立. 2009. 龙滩库区水库地震震源机制及应力场特征[J]. 地震地质,31(4):686–698. doi: 10.3969/j.issn.0253-4967.2009.04.012 Chen H L,Zhao C P,Xiu J G,Chen Z L. 2009. Study on the characteristics of focal mechanisms of reservoir induced earthquakes and stress field in the Longtan reservoir area[J]. Seismology and Geology,31(4):686–698 (in Chinese).
郭增建. 1958. 由地震波初动求断层面方法的一些推广和改进[J]. 地球物理学报,7(1):20–30. Kuo T C. 1958. Notes on the fault-plane determinations by means of initial motions of earthquakes[J]. Acta Geophysica Sinica,7(1):20–30 (in Chinese).
李强,李军,袁丽文,邱毅,李锋,秦双龙,曹轶. 2015. 福建仙游震群序列ML≥4.0事件震源机制与序列活动特征[J]. 地震,35(4):147–156. doi: 10.3969/j.issn.1000-3274.2015.04.016 Li Q,Li J,Yuan L W,Qiu Y,Li F,Qin S L,Cao Y. 2015. Focal mechanisms and activity of ML≥4.0 events in the 2013 Xianyou earthquake swarm sequence in Fujian Province[J]. Earthquake,35(4):147–156 (in Chinese).
毛松林,蔡欣欣. 2014. 仙游震群型地震的特点及地震机理分析[J]. 华南地震,34(1):86–93. Mao S L,Cai X X. 2014. Analysis of the characteristics and mechanism of xianyou earthquake swarm[J]. South China Journal of Seismology,34(1):86–93(in Chinese).
马志江,钟羽云,韩用兵,叶建青,张震峰,林胜法,徐梦林,刘倩倩. 2016. 温州珊溪水库诱发地震构造条件[J]. 地球科学,41(8):1413–1423. Ma Z J,Zhong Y Y,Han Y B,Ye J Q,Zhang Z F,Lin S F,Xu M L,Liu Q Q. 2016. The tectonic conditions of Shanxi reservoir induced earthquake in Wenzhou[J]. Earth Science,41(8):1413–1423 (in Chinese).
秦双龙,李强,邱毅,廖丽霞,陈莹. 2015. 福建仙游震群重定位及序列特征分析[J]. 地震地磁观测与研究,36(3):67–71. doi: 10.3969/j.issn.1003-3246.2015.03.013 Qin S L,Li Q,Qiu Y,Liao L X,Chen Y. 2015. Relocation and sequence analysis for Xianyou earthquake swarm in Fujian Province[J]. Seismological and Geomagnetic Observation and Research,36(3):67–71 (in Chinese).
秦双龙,邱毅. 2018. 福建仙游震群序列发震构造分析研究[J]. 地震工程学报,40(6):1306–1311. doi: 10.3969/j.issn.1000-0844.2018.06.1306 Qin S L,Qiu Y. 2018. Analysis of the seismogenic structure of the Xianyou earthquake swarm sequence in Fujian Province[J]. China Earthquake Engineering Journal,40(6):1306–1311 (in Chinese).
邱毅,李军,康兰池,袁丽文. 2014. 福建仙游地震序列的震源机制解[J]. 中国地震,30(2):280–288. doi: 10.3969/j.issn.1001-4683.2014.02.016 Qiu Y,Li J,Kang L C,Yuan L W. 2014. Focal mechanism of the Xianyou earthquake sequence,Fujian[J]. Earthouake Research In China,30(2):280–288 (in Chinese).
苏珊,韩立波,郭祥云. 2020. 溪洛渡水库近场区蓄水前后震源机制及应力场研究[J]. 地震研究,43(2):402–411. doi: 10.3969/j.issn.1000-0666.2020.02.024 Su S,Han L B,Guo X Y. 2020. Study on focal mechanism and stress field before and after the impoundment in the vicinity of the Xiluodu reservoir[J]. Journal of Seismological Research,43(2):402–411 (in Chinese).
万永革,沈正康,兰从欣. 2006. 根据走滑大地震前后应力轴偏转和应力降求取偏应力量值的研究[J]. 地球物理学报,49(3):838–844. doi: 10.3321/j.issn:0001-5733.2006.03.028 Wan Y G,Shen Z K,Lan C X. 2006. Deviatoric stress level estimation according to principle axes rotation of stress field before and after large strike-slip type earthquake and stress drop[J]. Chinese Journal of Geophysics,49(3):838–844 (in Chinese).
万永革,沈正康,刁桂苓,王福昌,胡新亮,盛书中. 2008. 利用小震分布和区域应力场确定大震断层面参数方法及其在唐山地震序列中的应用[J]. 地球物理学报,51(3):793–804. doi: 10.3321/j.issn:0001-5733.2008.03.020 Wan Y G,Shen Z K,Diao G L,Wang F C,Hu X L,Sheng S Z. 2008. An algorithm of fault parameter determination using distribution of small earthquakes and parameters of regional stress field and its application to Tangshan earthquake sequence[J]. Chinese Journal of Geophysics,51(3):793–804 (in Chinese).
万永革,吴逸民,盛书中,沈正康,万迪. 2011. P波极性数据所揭示的台湾地区三维应力结构的初步结果[J]. 地球物理学报,54(11):2809–2818. doi: 10.3969/j.issn.0001-5733.2011.11.011 Wan Y G,Wu Y M,Sheng S Z,Shen Z K,Wan D. 2011. Preliminary result of Taiwan 3-D stress field from P wave polarity data[J]. Chinese Journal of Geophysics,54(11):2809–2818 (in Chinese).
万永革. 2019. 同一地震多个震源机制中心解的确定[J]. 地球物理学报,62(12):4718–4728. doi: 10.6038/cjg2019M0553 Wan Y G. 2019. Determination of center of several focal mechanisms of the same earthquake[J]. Chinese Journal of Geophy sics,62(12):4718–4728 (in Chinese).
万永革. 2024. 震源机制水平应变花面应变的地震震源机制分类方法及序列震源机制总体特征分析[J]. 地球科学,49(7):2675–2684. Wan Y G. 2024. Focal mechanism classification based on areal strain of the horizontal strain rosette of focal mechanism and characteristic analysis of overall focal mechanism of the earthquake sequence[J]. Journal of Earth Science,49(7):2675–2684 (in Chinese).
王鹏,孙冬军,任烨,于俊谊,温燕林,朱艾斓. 2017. 浙江珊溪水库2014年震群活动发震机理研究[J]. 地震学报,39(5):648–658. doi: 10.11939/jass.2017.05.002 Wang P,Sun D J,Ren Y,Yu J Y,Wen Y L,Zhu A L. 2017. Investigation on the seismogenic mechanism of the 2014 earthquake swarm in Shanxi reservoir,Zhejiang Province[J]. Acta Seismologica Sinica,39(5):648–658 (in Chinese).
许忠淮,汪素云,黄雨蕊,高阿甲. 1989. 由大量的地震资料推断的我国大陆构造应力场[J]. 地球物理学报,32(6):636–647. doi: 10.3321/j.issn:0001-5733.1989.06.004 Xu Z H,Wang S Y,Huang Y R,Gao A J. 1989. The tectonic stress field of Chinese continent deduced from a great number of earthquakes[J]. Acta Geophysica Sinica,32(6):636–647 (in Chinese).
阎春恒,周斌,陆丽娟,孙学军,文翔. 2015. 龙滩水库蓄水后库区中小地震震源机制[J]. 地球物理学报,58(11):4207–4222. doi: 10.6038/cjg20151127 Yan C H,Zhou B,Lu L J,Sun X J,Wen X. 2015. Focal mechanisms of moderate and small earthquakes occurred after reservoir recharge in the Longtan reservoir region[J]. Chinese Journal of Geophysics,58(11):4207–4222 (in Chinese).
杨福平,余刚群,侯林峰,周昕. 2019. 大地电磁测深在珊溪水库诱发地震研究中的应用[J]. 工程地球物理学报,16(2):203–210. doi: 10.3969/j.issn.1672-7940.2019.02.011 Yang F P,Yu G Q,Hou L F,Zhou X. 2019. The application of magnetotelluric sounding (MT) to research of reservoir-induced earthquake in Shanxi region[J]. Chinese Journal of Engineering Geophysics,16(2):203–210 (in Chinese).
余海琳,万永革,黄少华,崔华伟. 2021. 利用P波初动数据研究2021年云南漾濞MS6.4地震序列震源机制解及应力场[J]. 地震研究,44(3):338–347. Yu H L,Wan Y G,Huang S H,Cui H W. 2021. Study on focal mechanism solution and stress field of the 2021 Yangbi,Yunnan MS6.4 earthquake sequence using P-wave first motion data[J]. Journal of Seismological Research,44(3):338–347 (in Chinese).
袁丽文,李强,陈彩虹. 2016. 仙游地震序列中小地震震源机制解特征[J]. 防灾减灾学报,32(3):99–104. Yuan L W,Li Q,Chen C H. 2016. Focal mechanism solution characteristics of moderate and small earthquakes in Xianyou area[J]. Journal of Disaster Prevention and Reduction,32(3):99–104 (in Chinese).
张帆,朱新运,钟羽云. 2013. 珊溪水库地震小震震源机制解特征研究[J]. 地震工程学报,35(3):684–691. doi: 10.3969/j.issn.1000-0844.2013.03.0684 Zhang F,Zhu X Y,Zhong Y Y. 2013. Study on the characteristic of focal mechanism solutions of the Shanxi reservoir earthquake[J]. China Earthquake Engineering Journal,35(3):684–691 (in Chinese).
张路. 2008. 福建东南沿海盆地第四纪构造运动模式与动力学成因[D]. 北京:中国地震局地质研究所;141−151. Zhang L. 2008. A Kinematic Model and Dynamic Cause of Quaternary Tectonic Movement of Southeastern Coastal Basins in Fujian Province[D]. Beijing:Institute of Geology,China Earthquake Administration;141−151 (in Chinese).
Byerly P. 1938. The earthquake of July 6,1934:Amplitudes and first motion[J]. Bull Seismol Soc Am,28(1):1–13. doi: 10.1785/BSSA0280010001
Hodgson J H,Storey R S. 1953. Tables extending byerly’s fault-plane techniques to earthquakes of any focal depth[J]. Bull Seismol Soc Am,43(1):49–61. doi: 10.1785/BSSA0430010049
Huang J L,Zhao D P. 2006. High-resolution mantle tomography of China and surrounding regions[J]. J Geophys Res:Solid Earth,111(B9):B09305.
Waldhauser F,Ellsworth W L. 2000. A double-difference earthquake location algorithm:Method and application to the northern Hayward fault,California[J]. Bull Seismol Soc Am,90(6):1353–1368. doi: 10.1785/0120000006
Wan Y G. 2010. Contemporary tectonic stress field in China[J]. Earthq Sci,23(4):377–386. doi: 10.1007/s11589-010-0735-5
Wan Y G,Sheng S Z,Huang J C,Li X,Chen X. 2016. The grid search algorithm of tectonic stress tensor based on focal mechanism data and its application in the boundary zone of China,Vietnam and Laos[J]. Journal of Earth Science ,27(5):777–785. doi: 10.1007/s12583-015-0649-1
计量
- 文章访问数: 186
- HTML全文浏览量: 53
- PDF下载量: 72
