福建仙游震群的应力场及其成因分析

余海琳, 万永革, 崔华伟, 王晓山, 黄少华

余海琳,万永革,崔华伟,王晓山,黄少华. 2024. 福建仙游震群的应力场及其成因分析. 地震学报,46(4):600−613. DOI: 10.11939/jass.20220204
引用本文: 余海琳,万永革,崔华伟,王晓山,黄少华. 2024. 福建仙游震群的应力场及其成因分析. 地震学报,46(4):600−613. DOI: 10.11939/jass.20220204
Yu H L,Wan Y G,Cui H W,Wang X S,Huang S H. 2024. Analysis of stress field and its causes of the Xianyou earthquake swarm in Fujian. Acta Seismologica Sinica46(4):600−613. DOI: 10.11939/jass.20220204
Citation: Yu H L,Wan Y G,Cui H W,Wang X S,Huang S H. 2024. Analysis of stress field and its causes of the Xianyou earthquake swarm in Fujian. Acta Seismologica Sinica46(4):600−613. DOI: 10.11939/jass.20220204

福建仙游震群的应力场及其成因分析

基金项目: 国家自然科学基金(42174074, 41674055)、中央高校科研业务费专项(ZY20215117)、河北省地震科技星火计划(DZ20200827053)和河北省地震动力学重点实验室开放基金(FZ212105)共同资助
详细信息
    作者简介:

    余海琳,助理工程师,主要从事构造应力场方面的研究,e-mail:1534541141@qq.com

    通讯作者:

    万永革,博士,研究员,主要从事构造应力场、地震应力触发等方面的研究,e-mail:wanyg217217@vip.sina.com

  • 中图分类号: P315.2

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.

  • 时间同步技术是多种行业的基础支撑技术,它所提供的高精度同步时间是科学实验和工程技术等方面的基本物理量。随着高精度同步授时在各个行业的需求愈加广泛,科研、生产和生活中时间的同步性和稳定性也越来越重要,例如电力网的时间精确和统一就是电力系统自动化和安全运行的要素之一,时分同步码分多址(time division-synchronous code division multiple access,缩写为TD-SCDMA)等通讯网络系统一般需要微秒级的时间同步,准确统一的时间也是利用各类地震观测系统进行地震科学研究的前提。因此,自主研究简单实用高精度的同步授时系统并将其应用于各类科研、生产、生活中具有重要意义,特别是对于地震研究而言,在各类地震动观测台阵、强震动观测台阵及边坡滑坡震动监测系统的建设与长期运行中,获取更高质量的同步地震动监测数据,对我国地震监测预报和地球科学研究等都非常重要。

    地震动监测系统对系统内仪器的时间同步都具有一定的精度要求,其中小孔径地震台阵和密集台阵等对地震仪器时间同步的精度要求更高。小孔径地震台阵是在与所观测地震波波长相当的孔径范围内有规则排列的一组地震计,采用独特的地震数据处理方法对这些地震计的输出信号进行各种组合分析。将各子台的数据时间对齐后进行抑制地面噪声、压低干扰背景、提高信噪比等方面的处理,从而提升地震监测能力,实现提取更微弱地震信号的目的。故台阵内各地震计时间同步误差将极大地影响微弱信号的提取效果。Meng等(2014)提出了利用近断层密集台阵波束形成技术反演震源方位角的方法来实时估测震源的破裂长度,用美国的UPSAR强震动台阵记录对多次地震进行破裂长度的估测,取得了较好的结果。该方法需要小孔径地震台阵提供精确的走时,但我国现有的一些强震动观测台阵中各观测子台的观测仪器均为独立授时,各子台记录可能存在时钟同步误差问题,因此对密集台阵中台站之间的时间延迟及其精度的确定有一定影响(刘辰等,2018)。此外,水库大坝等地震动监测系统的地震动监测仪器安装于密闭廊道内,不具备使用传统的全球定位系统和北斗卫星导航系统(Global Positioning System/BeiDou Navigation Satellite System,缩写为GPS/BDS)授时设备的条件,因此水库大坝、水电站等结构设施的地震动监测系统也需要一种更实用的授时系统。鉴于上述情况,本文拟提出一种基于实时光纤通信的多通道同步授时系统,以期解决各类分布式地震数据采集系统的高精度同步授时问题。

    目前各类地震观测仪器及地震动观测台阵常用的时间同步技术主要有短波(长波)授时技术、卫星授时技术、网络授时技术和IRIG-B (Inter-Range Instrumentation Group-B)码授时技术等,其中短波(长波)授时技术是利用不同波长信号通信实现地震动观测台阵中多个观测设备时间同步的方法,该种方法授时简单但授时精度有限,短波授时精度一般为毫秒级,长波授时精度为微秒级,一般仅用于野外特殊场地的部分观测设备。多数地震动观测台阵采用GPS/BDS卫星授时技术的地震观测设备,其内配置GPS/BDS接收模块,当GPS/BDS卫星天线接收到足够多的空间卫星信号后(王向军等,2016郭长发,2018),地震观测设备就能够获取到准确而无累积误差的高精度时间,进而实现同一区域的多台地震观测设备的时间同步。

    网络时间同步授时技术是基于网络传输的一种授时技术,如基于网络时间协议(network time protocol,缩写为NTP)、简单网络时间协议(simple network time protocol,缩写为SNTP)、网络测量和控制系统的精密时钟同步协议标准等标准或协议的授时技术(魏丰,孙文杰,2009彭栋,郭伟,2018)。网络时间协议是通过NTP协议和时间源进行时间校准,根据服务器与客户端往返报文来确定两者之间的时钟差值和报文在网路传输中的延迟向网络内的其他客户端提供精准时间服务。这种技术是可以跨越广域网、局域网的复杂时间协议,通常可以获得毫秒级的同步精度。网络授时技术也可以用硬件辅助解决网络协议栈和以太网中交换器、路由器产生的时延稳定问题(马红皎等,2014翟学明等,2016),用支持网络协议的硬件芯片在以太网口的物理层设计实时时钟、打时间戳,从而将时延降低到极小的数量级(小于1 µs)。

    虽有如上多种同步授时方法,但部分地震观测设备所处地理位置特殊,无法架设GPS/BDS接收天线,且其它授时方法的授时精度无法达到高精度同步地震观测的需求,因此本文将提出以基于光纤通信的多通道同步授时技术实现多通道分布式地震数据采集系统的高精度同步授时,并将其应用于小孔径地震台阵监测系统的授时。

    光信号在光纤中传播的速度约为2×108 m/s,长度为l (单位为m),光纤中单程光信号传播的时间为tl/(2×108) s。若光信号从光纤初始端发送并在光纤中匀速传播,该光信号到达光纤接收端即刻沿光纤原路返回至光纤初始端,那么在此传播过程中,光信号所需的传输时间为t=2l/(2×108) s。当所使用的光纤长度固定时,光信号在该段光纤中的传播时间即可确定,因此可以设计高精度时间测量电路来测量该传输时间。若一种设备的光信号发送端的发送时间可控、光信号接收端识别与返回电路的时间可控,则可以精确地测量光信号在不同长度光纤中的传输时间,并可以利用延时测量的时间值对授时操作进行时间校正,所以基于此原理设计一种多通道时间同步授时系统来实现分布式地震数据采集设备的高精度远距离同步授时是可行的。

    通过研究授时秒脉冲发送和接收时间测量方法、纳秒级定时器设计方法、授时秒脉冲及对应时间数据的编译与解析方法,设计一套高精度时间同步授时系统,具体步骤如下:① 通过高速可编程电路设计时间测量模块,实时测量授时秒脉冲的发送时间与返回时间,计算出该授时通道的授时秒脉冲全路径时间的延迟值(图1a);② 以内部高精度高频时钟信号为源,用硬件逻辑编程设计纳秒级精度的定时器,并在下一个授时秒脉冲到达前提前ta s 完成延迟校正的秒脉冲发送(图1b),从而达到补偿该授时通道的系统延时及光纤传输线路延时的目的;③ 图1atd为某授时通道的授时秒脉冲传输全路径的时间延迟值,图1b中第一组秒脉冲信号为同步授时端发送出去的、经过时间校正的授时秒脉冲信号,第二组秒脉冲信号为时间接收模块端接收到的授时秒脉冲信号,其中ta=0.5td。此秒脉冲信号理论上与该授时通道的初始发送秒脉冲时间同步。

    图  1  同步授时秒脉冲(a)及延时校正后的授时秒脉冲(b)
    Figure  1.  Synchronous timing second pulse (a) and timing second pulse after delay correction (b)

    本文设计的高精度同步授时系统基于GPS卫星授时技术,GPS接收模块在卫星条件及气候良好条件下,与GPS时间同步精度为纳秒级别,这种精度的时钟作为本系统的时间基准能够满足需求。高精度授时系统的结构如图2所示,包含同步授时中心、同步授时模块、光纤和时间接收模块。该系统进行同步授时脉冲传输的延迟时间测量及延迟时间校正、同步授时脉冲发送及时间信息编译、光电转换及光信号传输数据、时间接收模块识别授时脉冲及时间信息,自动完成多路时间接收模块端的同步授时。

    图  2  时间同步授时系统结构图
    Figure  2.  Structure chart of synchronous timing system

    在同步授时系统中,同步授时中心通过串口连接多个同步授时模块,每路同步授时模块由光纤连接对应的时间接收模块。多个同步授时模块从同步授时中心模块获得高度同步的授时秒脉冲信号,之后通过串口获得秒脉冲对应的时间信息,再根据如上时间信息同步授时模块完成与同步授时中心模块的时间同步,并进一步与时间接收模块配合完成传输延迟时间测量、同步授时及校正、通信数据发送等工作。

    同步授时中心模块是该同步授时系统的时间源,它管理高精度时间源、多路授时秒脉冲输出、多路串口通信端口,通过串口实现对多个同步授时模块的并行管理,即通过并行控制将时间信号以同步授时脉冲、时间信息发送至多路同步授时模块,从而完成整个系统的同步授时脉冲管理、时间数据发送等功能。同步授时中心模块的结构见图3

    图  3  同步授时中心模块结构图
    Figure  3.  Structure chart of synchronous time center mudule

    同步授时模块是同步授时系统中的一个关键模块,负责完成授时通道的授时脉冲传输时间测量及延时校正等功能,具体包括同步授时脉冲发送时间测量、经授时通道返回的同步脉冲时间测量、纳秒级定时器、同步授时脉冲及对应时间数据的编译、正常通信数据的收发管理等功能。

    同步授时模块内部设计分为逻辑控制区、数据发送区、数据接收区和光电转换区,其原理图见图4。逻辑控制区包括控制器、定时器和时间测量,主要负责逻辑控制、同步授时脉冲发送和返回时间的测量;发送区包括通信数据缓存区和发送控制器,负责将正常通信信息对外发送,其中控制器具有发送使能功能(见TxE控制线);接收区包括数据缓存区和接收控制器,负责时间接收模块端发送来的通信数据的管理,其中控制器有接收缓存的控制使能(见RxE控制线);光电转换区为高度集成的光电收发一体接收器,支持高速远距离数据传输。同步授时模块电路由高速ARM控制器、可编程器件、时钟管理模块、串口模块和光电转换模块组成。通过高速可编程电路设计时间测量模块,实时测量同步授时脉冲的发送时间与返回时间,计算出该授时通道的全路径时间延迟值(图1a)。具体步骤如下:① 以内部高精度高频时钟信号为源,用硬件逻辑编程设计纳秒级的定时器,在下一个授时秒脉冲到达前自动完成延迟校正的授时脉冲发送(图1b),从而达到补偿该授时通道的系统延时及光纤传输延时等;② 用硬件编程设计纳米级的定时器,对于两个授时秒脉冲之间的1 s,采用同步授时模块的定时器及可编程器件进行精确延时,生成一个初值为(1-ta) s的定时器;③ 启动定时器工作,当定时时间到达,由同步授时模块发出该通道同步授时脉冲;④ 该同步授时脉冲发出后,所对应串口发出该同步授时脉冲对应的时间值;⑤ 将同步授时脉冲及其时间数据嵌入正常通信数据,遵循通信协议命令经光电转换模块和光纤发送出去,且以时间同步脉冲电信号形式发送;⑥ 在同步授时脉冲发送前通过内部缓存器暂存正常通信数据,同步授时脉冲发送结束、停止数据缓存、恢复数据发送。

    图  4  同步授时模块原理图
    Figure  4.  Schematic diagram of synchronous timing module

    当时间同步脉冲信号到来,时间信息依次写入控制器、定时器和时间测量等三个模块。控制器通过TxE使能暂停其它通信数据发送,定时器关闭开关K1和开关K2,时间同步脉冲经由K1线路直接发送。控制器通过TxE使能失效而恢复发送,定时器打开开关K1,发送时间数据和通信数据,带有同步脉冲标识的同步授时数据帧编译完毕。

    当返回时间同步脉冲信号到来,经由开关K2直接送入时间测量模块,控制器通过RxE使能暂停数据接收。返回时间同步脉冲信号结束,定时器控制开关K2打开,控制器通过RxE使能恢复数据接收。启动时间同步脉冲信号与返回时间同步脉冲信号的比较过程,即完成传输延迟时间测量,该时间延迟值存入定时器。

    时间接收模块是应用在远端待授时采集设备端的硬件模块,它自动识别同步授时脉冲信号和带有同步脉冲标志的时间数据帧,并提取同步脉冲信号、时间信息和数据信息,自动向光纤发送返回的同步脉冲信号,并向采集设备发送授时脉冲、时间数据及管理通信数据等。同步时间接收模块必须与同步授时模块成对使用,完成一路数据采集设备的时间授时。

    同步时间接收模块电路采用高速可编程器件及相关配件设计完成,主要包括控制器、接收区(缓存和接收)、发送区(发送和缓存)和光电转换区,同步时间接收模块的内部元件有ARM控制器、逻辑控制器、串口控制器、脉冲转发电路及高速光电收发一体接收器。其中控制器可以通过使能信号RxE和TxE控制通信数据接收和发送的启停,也控制开关KaKb的操作,接收同步授时脉冲信号以及返回同步授时脉冲信号。同步时间接收模块的设计结构见图5

    图  5  时间接收模块原理图
    Figure  5.  Schematic diagram of time receiving module

    同步授时脉冲信号经光电转换后到达控制器并完成授时脉冲识别后,控制器关闭脉冲控制开关Ka和发送数据控制开关Kb,同步授时脉冲经由开关Ka线路被置入待授时采集器的秒脉冲输入端,同步经过接收串口将对应授时脉冲的时间数据送入时间数据传输串口。同时经控制开关Kb线路,直接通过光纤将该授时脉冲沿原传输路径返回,即该授时脉冲经电光转换、光纤、同步授时模块形成延时测量脉冲(返程的时间同步授时脉冲信号)。

    同步授时脉冲信号到达控制器时,控制使能信号RxE和TxE暂停数据发送与接收。

    同步授时脉冲信号发送结束后,控制器关闭,打开开关KaKb,控制使能信号RxE和TxE恢复数据发送与接收,一组接收时间同步脉冲工作时序结束。

    本文的时间同步授时系统虽由多个分立电路模块组成,但进行同步授时操作需要多组模块协同操作共同完成,各组模块分别完成自身通道的同步授时操作。

    同步授时模块工作分为测量时钟校准(clock correction)、线路延时校准(delay correction)、同步信息发送和微帧数据发送几个步骤,其中:测量时钟校准是指利用时间同步信息调整和校准模块内的时钟误差;线路延时校准是指同步授时模块进行测量线路的延时,然后将其存储为下发同步授时秒脉冲信号时的补偿线路延时值;同步信息发送是指收到同步授时秒脉冲信号时,对其执行测量时钟校准算法和线路延时校准算法后再发送授时秒脉冲,以保证远端待授时设备收到的同步授时秒脉冲信号与同步授时中心的授时秒脉冲信号一致,其中同步授时秒脉冲信号嵌入授时秒脉冲信号,并且同步授时秒脉冲信号以硬件信号的形式进行收发;微帧数据发送则是指根据发送或接收信号的类型来确定串口缓存内部的通信数据发送与否。

    时间接收模块工作包括配合同步授时模块校准线路延时(delay correction of synchronous time module,缩写为DCS)、同步信息接收和微帧数据接收。配合同步授时模块校准线路延时(DCS)是指当同步授时模块申请测量线路延时的时候,时间接收模块无延时返回接收到的授时秒脉冲信号;同步信息接收是指时间接收模块识别出授时秒脉冲信号后,以硬件信号的形式无延时地把授时秒脉冲信号输出到授时秒脉冲输出端口,提取时间脉冲对应时间信息数据并将其发送到串口1;当上述工作完成后,时间接收模块恢复微帧数据接收,并把收(发)的正常数据由串口2对应收发处理。同步授时系统工作流程图见图6,在此流程下同步授时系统中各个模块中的授时脉冲时序见图7

    图  6  同步授时工作流程图
    Figure  6.  Flow chart of synchronous timing program
    图  7  校正前(a)、后(b)的各模块时间脉冲时序图
    Figure  7.  Uncorrected (a) and corrected (b) timing pulse sequence diagrams of each module

    此处以一路授时通道的同步授时工作时序为例介绍各模块中授时秒脉冲信号的工作时序。未进行同步授时操作时各个模块的时间脉冲工作时序见图7a,图中t1为同步授时模块发送脉冲延迟时间,该延迟时间内同步授时模块将完成上一微帧数据通信,然后停止数据通信;t2为同步授时模块授时秒脉冲信号在光纤中正向传输的延迟时间值;t3为授时秒脉冲从时间接收模块返回后在光纤中反向传输的时间延迟值,一般情况下t2t3。当完成一次同步授时操作后,各个模块的时间脉冲将按照图7b所示的时序进行工作,从该时序图可以看出,时间接收模块端的秒脉冲已经与同步授时模块输出的秒脉冲信号达到同步。

    基于以上技术研究,本文采用授时中心模块、同步授时模块、时间接收模块及通信光纤搭建了一套同步授时系统,并且以GPS授时模式作为授时中心模块的时间源,授时通道的通信光纤均采用2 km长的单模通信光纤。启动该同步授时系统连续运行,采用高精度时间测量设备连续记录授时通道的授时脉冲及对应时间信息。测试试验所使用通用计数器的频率范围为10 Hz—1 300 MHz,其频率稳定度为3×10−10/s,电力专用卫星接收机的工作频率为5 MHz和10 MHz,其稳定度为1×10−11/s。同步授时系统工作时,定时器的时钟频率为100 MHz,试验时实测得到光纤传输时间的延迟值td约为20 μs,时间延迟校正值ta约为10 μs。

    当授时中心模块的GPS信号有效后,该时间同步授时系统的同步授时工作启动。此时记录授时中心模块输出的授时脉冲和各路通道时间接收模块的时间同步脉冲的连续测试时长为1 h,在同次试验中相同条件下,同时采集三路通道的同步时间接收模块的时间同步脉冲,分别统计授时中心模块输出的授时脉冲与各路通道时间接收模块的时间脉冲之差,以此时间差值作为各路授时通道的授时时间误差。图8给出了三路授时通道的授时时间误差,其中通道一、二和三的授时时间误差均方根值分别为126.22 ns,129.19 ns和129.85 ns,三路通道的授时误差均方根值均低于200 ns。

    图  8  授时中心模块的授时脉冲与第一(a)、第二(b)和第三(c)通道时间接收模块的授时脉冲的时间差值数据曲线
    Figure  8.  The time difference data curves between the timing pulse of the timing central module and the timing pulse of the first (a),second (b) and third (c) channel time receiving module

    本文采用时间同步授时技术研制了一套高精度时间同步授时系统,并搭建了一套同步授时实验系统进行了连续授时观测试验。对三路通道的同步授时脉冲的连续观测记录的分析计算结果表明,本文所研制的时间同步授时系统能够连续稳定工作,且在通讯光纤为2 km长时,其1 h连续观测记录的同步授时误差均小于200 ns。该同步授时误差主要来源于两种硬件模块内部ARM控制器指令执行时间误差,具体为同步授时模块内部定时器计数溢出时向ARM控制器发出中断请求,当ARM控制器正在执行指令或中断服务程序时,该ARM控制器响应定时器溢出中断会产生延迟,对同步授时脉冲的发送造成影响。定时器计数溢出中断与ARM控制器响应该中断的时间误差具有非固定性特点,这种非固定性给误差分析、误差补偿也带来了困难。为控制发送同步授时脉冲时间的不确定性,后续还要选用FPGA模块替代ARM控制器设计同步授时模块及时间接收模块的硬件电路,以降低模块操作时间不确定性所带来的同步误差,从而进一步降低同步授时系统的时间偏差。

  • 图  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%的置信水平下应力场的最大剪应力节面;黄色箭头表示最大剪应力的滑动方向。S3S2S1及标注周边的闭合曲线分别表示主压应力轴、中间应力轴和主张应力轴的置信范围。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 S3S2S1 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

  • 期刊类型引用(6)

    1. 张志宏,侯作亮,杨士超,李梦莹,黄明威. 基于AMT的地电阻率观测场地优化选择研究——以辽宁新城子深井地电阻场地为例. 大地测量与地球动力学. 2025(04): 436-440 . 百度学术
    2. 王丽,李伟,沈红会,孙君嵩,李鸿宇. 地电阻率受高压输电线路架线施工干扰的判定. 华北地震科学. 2024(02): 102-108 . 百度学术
    3. 曹原. 高导电率地质区间地铁杂散电流控制技术研究. 机电信息. 2024(11): 85-89 . 百度学术
    4. 张宇,王兰炜,胡哲,张世中,张兴国,娄晓宇. 加权平均算法在地电阻率日均值计算中的应用. 地震. 2023(02): 25-37 . 百度学术
    5. 李雪浩,何思源,李国超,刘华姣,廖绍欢,赵乃千. 成都地震监测中心站地电阻率交、直流观测系统对比分析. 地震地磁观测与研究. 2023(05): 118-124 . 百度学术
    6. 刘建波,雷生学,张明东,刘文兵,刘金城,马义山. 用自然电位估算地铁干扰大小——以天津塘沽地震台为例. 地震地磁观测与研究. 2023(05): 107-117 . 百度学术

    其他类型引用(1)

图(8)  /  表(3)
计量
  • 文章访问数:  186
  • HTML全文浏览量:  53
  • PDF下载量:  72
  • 被引次数: 7
出版历程
  • 收稿日期:  2022-11-02
  • 修回日期:  2023-03-15
  • 网络出版日期:  2023-09-27
  • 刊出日期:  2024-07-14

目录

/

返回文章
返回