Ground motion prediction of MW≥7.5 on Xiadian fault
-
摘要: 针对夏垫断裂开展了MW≥7.5地震动预测研究。首先基于全破裂模式设定震源（使其尽可能涵盖夏垫断裂的未知信息）模拟得到夏垫断裂发生MW≥7.5地震时研究区域内的地面地震动场，进而依据分位数筛选出各场点的地震动空间分布，讨论了包含不确定震源下的加速度峰值和速度峰值的分布特征，结果显示当夏垫断裂发生MW7.9地震时，通州城区、北京中心城区均会发生强烈的运动。之后对比讨论了仿真震源下MW7.5地震所引起的地面运动场的空间变化，结果显示对于同等震级而言，两种震源的模拟结果可以相互印证。
-
关键词:
- 夏垫断裂 /
- NNSIM随机有限断层法 /
- 震源设定 /
- 地震动预测
Abstract: The prediction of ground motion fields for potential earthquakes is important for urban planning and regional seismic hazard assessment. In this paper, the prediction of MW≥7.5 ground motion is carried out for the Xiadian fault. Firstly, the earthquake source is set based on the full rupture patterns to cover the unknown information of the Xiadian fault as much as possible. Then the ground motion field in the study area is simulated due to the MW≥7.5 earthquake on the Xiadian fault. Furthermore, the spatial distribution of ground motion at each site is filtered based on the quantile method, and the distribution characteristics of peak ground acceleration and peak ground velocity for the inclusion of uncertain sources are discussed. The results show that an MW7.9 earthquake on the Xiadian fault will produce the strong ground motion in Tongzhou district and Beijing center zone. After that, the spatial distribution of the ground motion due to the MW7.5 earthquake from the simulated source is discussed. The results illustrate that the simulated ground motion from the two types of sources can corroborate each other for the same magnitude. Our study provides a method for the earthquake hazard predictions due to the potential sources with some unknowns. -
引言
对于地震学而言，高频全球导航卫星系统（global navigation satellite system，简写为GNSS）是能够记录到地震波传播信息的一种新手段，其在强震中可直接提供时变位移场观测信息，弥补了传统观测手段在强震观测时出现诸如震级饱和、漂移等现象的不足。但由于高频GNSS数据存在信噪比低、信号稳定性差等问题，在预警中需要与测/强震观测手段联合使用。地震波信号往往是一种非线性、非平稳信号，故研究其频率的局部特性显得至关重要，这就需要用时频分析方法将一维时域信号转换到二维的时频平面。而傅里叶变换和反变换虽然能从时间域和频率域对平稳信号进行分析，却不能同时保留时间和频率信息。故此，本文拟以2008年汶川地震时获取的高频GNSS同震信号为例，通过对其进行广义S变换，系统分析变换后同震信号的时间域和频率域，以期能够更加准确地判断P波和S波到时。同时，广义S变换的使用也为高频GNSS在地震学中的应用提供了一条从多层面提取地震波形到时的新途径。
2. S变换
美国地球物理学家Stockwell等（1996）综合了小波变换和短时傅里叶变换的特点，提出了一种分析非平稳信号的新方法—S变换。其时频随着频率的不同而变化，具有多分辨率、逆变换唯一等特点，同时获得的二维时频谱也是基于傅里叶变换的基本理论。该方法已经在地震勘探、电能信号分析、故障检测等多个领域得到了广泛的应用（齐春燕等，2010；周竹生，陈友良，2011；郑成龙，王宝善，2015）。
S变换的表示形式为
$$ {{S}}(\tau {\text{，}}f) {\text{＝}} \int _{ {\text{－}} \infty }^{ {\text{＋}} \infty }h(t)\left\{ {\frac{{|f|}}{{\sqrt {2\pi } }}\exp \left[\frac{{ {\text{－}} {f^2}{{(t {\text{－}} \tau )}^2}}}{2}\right]\exp ( {\text{－}} {\rm i}2\pi ft)} \right\}{\rm d}t{\text{，}} $$ (1) 式中：τ为窗函数的中心点，控制窗函数在时间轴上的位置；f为采样频率；h（t）为原始信号。S变换使用的高斯窗函数gf （t）和基本小波ωf （t）分别定义为
$$ {{g}_f}(t) {\text{＝}} \frac{{|f|}}{{\sqrt {2\pi } }}\exp \left({\text{－}} \frac{{{t^2}{f^2}}}{2}\right){\text{，}} $$ (2) $$ {\omega _f}(t) {\text{＝}} \frac{{|f|}}{{\sqrt {2\pi } }}\exp \left( {\text{－}} \frac{{{t^2}{f^2}}}{2} {\text{－}} {\rm i}2\pi ft\right) {\text{＝}} {g_f}(t)\exp ( {\text{－}} {\rm i}2\pi ft){\text{，}} $$ (3) 从上式可以看到，S变换在分析高频信号或低频信号时具有良好的特性：对于高频信号，设计窗口比较窄，可满足高时间分辨率的要求；对于低频信号，设计窗口比较宽，可满足高频率分辨率的要求。
相比于单分辨率的短时傅里叶变换，S变换属于多分辨率时频分析方法，可以更加有效地对复杂的地震波信号进行分析。同时由于信号x（t）可以由S变换S（τ，f ）通过逆变换进行重构，S（τ，f ）的S逆变换为
$$ {\rm{x}}(t) {\text{＝}} \int_{ {\text{－}} \infty }^{ {\text{＋}} \infty } {\left[ {\int_{ {\text{－}} \infty }^{ {\text{＋}} \infty } {S(\tau{\text{，}} f\;)\rm d\tau } } \right]} \exp ({\rm i}2\pi ft){\rm {d}}f{\text{，}} $$ (4) S变换所具有的无损逆变换的特性，使得信号在经过时频域处理后可逆变换为时间信号。
由于S变换的基本高斯窗函数固定，因此在处理复杂的地震信号时，不能根据需求对频率域和时间域的分辨率进行调节，缺乏必要的灵活性。随着后期的不断发展，多位研究人员对S变换中的高斯窗函数和基本小波进行改变，构建出了广义的S变换（Mansinha等，1997；McFadden et al，1999；高静怀等，2003；Pinnegar，Mansinha ，2003；陈学华等，2006；齐春燕等，2010；周竹生，陈友良，2011）。
本文使用的广义S变换是通过两个调节因子对高斯窗进行改造（陈学华等，2006），以便根据实际要求对高斯窗作出调整，从而更好地分析复杂频率的地震波形，其表达式为
$$ {{S}}(\tau {\text{，}} f) {\text{＝}} \int\nolimits_{ {\text{－}} \infty }^\infty {x(t)\frac{{|{\lambda _{\rm a}}||f{|^p}}}{{2\pi }}} {{\rm exp}({\frac{{ {\text{－}} {\lambda ^2_{\rm{a}}}{f^{2p}}{{(t {\text{－}} \tau )}^2}}}{2}}}){{\rm exp}({ {\text{－}} {\rm i}2\pi ft}}){\rm d}t{\text{，}} $$ (5) 式中λa和p均为调节因子。选定p值后，当λa＞1时，时窗宽度随频率呈反比变化的速度加快；当λa＜1时，时窗宽度随频率呈反比变化的速度减慢；选定λa之后，通过调节p值可以对频率域的分辨率进行调整。引入λa和p两个参数之后，可以根据实际信号的频率分布特点和时频分析的侧重点，改变高斯窗函数随频率变化的规律，对时频分辨率进行调整，但由于其运算过程与标准的S变换相似，因此并不需要额外增加运算量。当两个调节因子λa和p均设置为1时，则为标准S变换。
3. 到时分析中的S变换
本文收集了2008年汶川MS8.0地震中重庆台网和昆明台网共28个站点的高频GNSS信号（采样率为1 Hz），采用GAMIT/GLOBK软件包（Bock，2000）中的高频GNSS处理模块TRACK （Chen，1998）对附近站点的高频GNSS数据进行解算。由于所选取数据按小时文件进行解算，时间最长为1小时，故未考虑昼夜变化等周期性较长的噪音。按照GNSS站点的位置及分布情况，选取15个站点的数据使用广义S变换进行分析，并挑选其中6个站点的同震波形进行展示，如图1所示。
根据距离震中位置的远近选取了bana （距震中约346 km）、djxm （距震中381 km）、fdlh （距震中477 km）和ghxx （距震中661 km）等4个站点的同震波形，对其进行广义S变换后识别出P波、S波到时并进行展示（图1），在计算P波和S波的理论到时中使用了频率-波数法（Zhu，2002），计算到时的速度模型为一维速度模型（吴建平等，2009；殷海涛等，2013；邓文哲等，2014；王小娜等，2015；杨彦明等，2016）。
bana站点距离震中346 km，该站点记录高频GNSS信号的频率为1 Hz，根据奈奎斯特（Nyquist）采样定理，其可记录到的地震波信号的最大频率为0.5 Hz。bana站点东西向同震信号广义S变换的结果表明（图2），地震信号的能量主要集中于0—0.3 Hz，所研究各站点的P波和S波理论到时均在300 s以内。为了更加直观地显示所研究的到时区域，在后续分析中，按照S变换后信号频率为0—0.3 Hz，取地震发生后的300 s进行分析。结果显示，bana站点的P波理论到时为56.7 s，S波理论到时为96.1 s。如图1所示，经过广义S变换后的图像显示，在P波和S波理论到时（图2）的位置，均有着较明显的能量汇集。在本文分析的地震中，P波和S波周期均较长、频率均相近，故P波的识别周期为20—50 s左右的初至拾取、高于背景噪声的能量汇集，S波能量为P波能量的2.5倍以上（水平向）。
![]()
图 2 bana站点东西向同震信号的S变换结果Figure 2. S transform of coseismic signal in the EW directionin at the bana station在使用广义S变换进行频率分析时，通过对两个调节因子λa和p进行调整，以获得最优的分辨率。以bana站东西向为例进行分析，结果如图3所示。
![]()
图 3 通过微调调节因子λa和p所产生的同震信号广义S变换结果展示Figure 3. S transform results of the same coseismic signal with different adjustment factor λa and p（a） λa＝1，p＝1.1；（b） λa＝1，p＝1.2；（c） λa＝1.05，p＝1.05；（d） λa＝1.05，p＝1.05在λa不变的情况下，提高因子p值，可以增加时间域的分辨率，详见图2和图3a。但是p值过高，会导致频率域的分辨率增加，时间域的分辨率降低，从而降低了对初至波的识别能力（图3b）。当调节因子λa＝1.05，p＝1.05时，对于P波和S波的识别均较为清晰。图3d将图3c进行三维展示，更加直观地展示了震后随着时间的变化，不同周期的波形所汇聚的能量大小。
对于震中距为381 km的djxm站点，其P波理论到时为62.5 s，S波理论到时为105.9 s。将P波和S波理论到时分别标注于S变换结果上（选用的调节因子为λa＝1.05，p＝1.05），可以看出，在同震信号中，djxm站记录的到时要比理论到时延迟约4—5 s （图4a），可能与本文使用的速度结构有一定关系，而djxm站在垂向上很难识别S波到时（图4b）。在原始波形记录中，无论垂直向还是水平向均已很难判定P波和S波到时。这表明通过微调广义S变换的两个调节因子λa和p，能够得到较好的到时分辨率。鉴于目前GNSS预警中的触发方式为阈值触发（Allen，Ziv，2011；Melgar et al，2012），很难通过反演的方式进行精确的震源定位。而相比于人工或者阈值判断，通过对地震信号进行广义S变换，能够对到时作出更为精准的判断。
![]()
图 4 各高频GNSS站点同震信号S变换结果展示（a） djxm站点EW向同震信号S变换；（b） djxm站点UD向同震信号S变换；（c） fdlh站点EW向同震信号S变换；（d） fdlh站点UD向同震信号S变换；（e） ghxx站点EW向同震信号S变换；（f） ghxx站点UD向同震信号S变换Figure 4. S treansform of coseismic signals for several stations（a） S transform of EW coseismic signal in djxm；（b） S transform of UD coseismic signal in djxm；（c） S transform of the EW coseismic signal in fdlh ；（d） S transform of the UD coseismic signal in fdlh；（e） S transform of the EW coseismic signal in ghxx；（f） S transform of the UD coseismic signal in ghxx对于震中距为477 km的fdlh站，P波和S波理论到时分别为79 s和133 s。东西向同震信号的广义S变换（选用的调节因子为λa＝1.05，p＝1.05）波形识别基本与理论到时一致，相差2—3 s （图4c）。fdlh站垂直向上的分辨率也可以满足波形识别的需求（图4d）。其结果较好，可能与fdlh站的背景噪声较低有关。三维图像中，面波的识别更加清晰（图4c）。
对于震中距为661 km的ghxx站，在同震信号中很难看到P波到时和S波到时。经过计算，其P波理论到时为108.3 s，S波理论到时为183 s。在图4e中可以看到，通过广义S变换后对S波的识别较好，对P波的识别稍显模糊。而对于垂直向（图4f），通过广义S变换之后，P波已经无法识别，S波的识别时间较理论地震波形要晚将近20 s。这可能与GNSS背景噪声较高，且垂直向的噪音水平高于水平向2—3倍有关。而同震信号的分析中，P波和S波的到时基本均无法分辨。这说明通过广义S变换，即使是震中距大于600 km，也能有效地提取到地震波到时。
4. 讨论与结论
对于实时解算的高频GNSS数据，其背景噪声相对较大，使用广义S变换，可以更加准确地得到P波到时。即使远距离的站点记录，在时频域中也可以很好地指示出P波初至到时。这为高频GNSS数据在地震学中的应用和分析提供了一种新的途径。其优势在于：
1） 广义S变换集合了短时傅里叶变换和小波变换的优点，在地球物理领域及电能分析领域取得了很多成果，但依旧存在很多应用前景。本文通过对高频GNSS同震信号进行S变换，从以往的只针对单独的时间域或频率域分析高频GNSS信号，转变为从时间域和频率域同时进行分析。从频率域的角度看，P波和S波的能量均汇聚于0.025 Hz左右，面波的能量汇聚于0.04—0.25 Hz并能够明显地看出面波是分为多个频率区间；从时间域的角度看，P波和S波的理论到时均有明显的能量汇集，而面波到时则更加明显。通过对不同频率的信号进行分析，认为近震范围内P波和S波实际到时与理论到时的分辨率相差仅为1—2 s，中远震其分辨率相差3—5 s，该方法可以作为高频GNSS信号的分析及各类波形到时的提取手段。
2） 通过对高频GNSS信号使用S变换可见，有效信号与噪音信号的分离较为明显，说明S变换在降噪去噪方面具有独特的优势。鉴于信号可以进行无损S变换和S逆变换，因此可以使用该方法提高高频GNSS信号的信噪比，尤其是在压制噪音方面具有很大的优势，下一步研究中将根据各波形的到时按时间提取波形并用于专门性的分析。
感谢两位审稿专家对本文提出的意见和建议。
① Wang Z Y，Li Y L，Wang W Q，Zhang W Q，Zhang Z G. 2022. Revisiting paleoearthquakes with computational seismology：A case study of the 1679 Sanhe-Pinggu earthquake （submitted）. -
![]()
图 7 模拟MW7.5加速度峰值及其对应的不同周期T的反应谱分布
Figure 7. The distribution of simulated acceleration peaks and response spectra with different period T for MW7.5
![]()
图 1 MW7.8地震多破裂模式下50%，60%，70%和80%分位数的地震动峰值（PGA）分布
Figure 1. Distribution of 50%，60%，70% and 80% quantile PGAs at MW7.8 by multiple rupture patterns
![]()
图 2 MW7.9地震多破裂模式下50%，60%，70%和80% 分位数的PGA分布
Figure 2. Distribution of 50%，60%，70% and 80% quantile PGAs at MW7.9 by multiple rupture patterns
![]()
图 2 MW7.9地震多破裂模式下50%，60%，70%和80%分位数的PGA分布
Figure 2. Distribution of 50%，60%，70% and 80% quantile PGAs at MW7.9 by multiple rupture patterns
![]()
图 3 断层走向为30°时50%，60%，70%和80%分位数下的地震动PGA分布
Figure 3. PGA distribution of ground motion at 50%，60%，70% and 80% quartile due to the fault with strike 30°
![]()
图 4 断层走向为45°时50%，60%，70%和80%分位数下的地震动PGA分布
Figure 4. PGA distribution of ground motion at 50%，60% ，70% and 80% quartile due to the fault with strike 45°
![]()
图 5 断层走向为55°时50%，60%，70%和80%分位数下的地震动PGA分布
Figure 5. PGA distribution of ground motion at 50%，60%，70%和80% quartile due to the fault with strike 55°
![]()
图 6 三个走向的加速度合成后50%，60%，70%和80%分位数的下地震动PGA分布
Figure 6. Distribution of PGAs of the synthesis records at 50%，60%，70% and 80% quartiles for the faults with three strikes
-
巴振宁,赵靖轩,吴孟桃,梁建文. 2022. 基于CPU-GPU异构并行的复杂场地近断层地震动谱元法模拟[J]. 地震学报,44(1):182–193. doi: 10.11939/jass.20210076 Ba Z N,Zhao J X,Wu M T,Liang J W. 2022. Simulation of near-fault ground motions in complex sites based on CPU-GPU heterogeneous parallelism by spectral element method[J]. Acta Seismologica Sinica,44(1):182–193 (in Chinese).
邓起东,张培震,冉勇康,杨晓平,闵伟,陈立春. 2003. 中国活动构造与地震活动[J]. 地学前缘,10(增刊1):66–73. Deng Q D,Zhang P Z,Ran Y K,Yang X P,Min W,Chen L C. 2003. Active tectonics and earthquake activities in China[J]. Earth Science Frontiers,10(S1):66–73 (in Chinese).
付长华,高孟潭,陈鲲. 2012. 北京盆地结构对长周期地震动反应谱的影响[J]. 地震学报,34(3):374–382. doi: 10.3969/j.issn.0253-3782.2012.03.009 Fu C H,Gao M T,Chen K. 2012. A study on long-period response spectrum of ground motion affected by basin structure of Beijing[J]. Acta Seismologica Sinica,34(3):374–382 (in Chinese).
高景华,徐明才,荣立新,柴铭涛,王广科,王小江,刘冠军. 2007. 探测城市断裂活动性的浅层地震方法技术[J]. 物探与化探,31(增刊1):4–8. Gao J H,Xu M C,Rong L X,Chai M T,Wang G K,Wang X J,Liu G J. 2007. The shallow seismic method for detecting city fault activity[J]. Geophysical and Geochemical Exploration,31(S1):4–8 (in Chinese).
高景华,徐明才,荣立新,柴铭涛,王广科,王小江,刘冠军. 2008. 利用地震剖面研究夏垫断裂西南段的活动性[J]. 地震地质,30(2):497–504. doi: 10.3969/j.issn.0253-4967.2008.02.014 Gao J H,Xu M C,Rong L X,Chai M T,Wang G K,Wang X J,Liu G J. 2008. Activity of the southwest segment of Xiadian fault investigated by seismic reflection profiling[J]. Seismology and Geology,30(2):497–504 (in Chinese).
高孟潭,俞言祥,张晓梅,吴健,胡平,丁彦慧. 2002. 北京地区地震动的三维有限差分模拟[J]. 中国地震,18(4):356–364. doi: 10.3969/j.issn.1001-4683.2002.04.005 Gao M T,Yu Y X,Zhang X M,Wu J,Hu P,Ding Y H. 2002. Three-dimensional finite-difference simulations of ground motions in the Beijing area[J]. Earthquake Research in China,18(4):356–364 (in Chinese).
高清武,李霓. 1998. 夏垫断裂氢释放与地震活动的关系[J]. 地震学报,20(4):432–436. Gao Q W,Li N. 1998. Relationship between H2 release and seismicity on Xiadian fault[J]. Acta Seismologica Sinica,11(4):501–505. doi: 10.1007/s11589-998-0097-4
高振寰. 1979. 唐山地震引起北京地区烈度异常区的原因分析[J]. 地震地质,1(2):74–82. Gao Z H. 1979. On distribution of intensity anomalies in the Beijing region caused by the Tangshan earthquake[J]. Seismology and Geology,1(2):74–82 (in Chinese).
何付兵,白凌燕,王继明,刘予,蔡向民,孙永华,张磊,方同明,郭高轩. 2013. 夏垫断裂带深部构造特征与第四纪活动性讨论[J]. 地震地质,35(3):490–505. doi: 10.3969/j.issn.0253-4967.2013.03.004 He F B,Bai L Y,Wang J M,Liu Y,Cai X M,Sun Y H,Zhang L,Fang T M,Guo G X. 2013. Deep structure and Quaternary activities of the Xiadian fault zone[J]. Seismology and Geology,35(3):490–505 (in Chinese).
何宏林,闵伟,原口强. 2008. 1679年三河—平谷8级地震破裂带的大地切片实验研究[J]. 地震地质,30(1):289–297. doi: 10.3969/j.issn.0253-4967.2008.01.021 He H L,Min W,Haraguchi T. 2008. Testing geo-slicer on the rupture of the M8 Sanhe-Pinggu earthquake of 1679[J]. Seismology and Geology,30(1):289–297 (in Chinese).
黄秀铭,汪良谋,徐杰,方仲景,张裕明,向家翠,王辉. 1991. 北京地区新构造运动特征[J]. 地震地质,13(1):43–51. Huang X M,Wang L M,Xu J,Fang Z J,Zhang Y M,Xiang J C,Wang H. 1991. Characteristics of neotectonic movement in Beijing area[J]. Seismology and Geology,13(1):43–51 (in Chinese).
江娃利,侯治华,肖振敏,谢新生. 2000. 北京平原夏垫断裂齐心庄探槽古地震事件分析[J]. 地震地质,22(4):413–422. doi: 10.3969/j.issn.0253-4967.2000.04.010 Jiang W L,Hou Z H,Xiao Z M,Xie X S. 2000. Study on paleoearthquakes of Qixinzhuang trench at the Xiadian fault,Beijing plain[J]. Seismology and Geology,22(4):413–422 (in Chinese).
焦青, 邱泽华. 2006. 北京平原地区主要活动断裂带研究进展[G]//地壳构造与地壳应力文集（18）. 北京: 地震出版社: 72–84. Jiao Q, Qiu Z H. 2006. Research progress of major active faults in Beijing plain area[G]//Bulletin of the Institute of Crustal Dynamics （18）. Beijing: Seismological Press: 72–84 （in Chinese）.
李赫,董一兵,王熠熙,彭研枫,刘双庆,吴博洋. 2020. 廊固凹陷深部剪切破裂构造的地震学证据[J]. 地球物理学报,63(2):492–504. doi: 10.6038/cjg2020N0192 Li H,Dong Y B,Wang Y X,Peng Y F,Liu S Q,Wu B Y. 2020. Seismological evidence for a deep-seated shear zone in the Langgu depression[J]. Chinese Journal of Geophysics,63(2):492–504 (in Chinese).
刘保金,胡平,孟勇奇,酆少英,石金虎,姬计法. 2009. 北京地区地壳精细结构的深地震反射剖面探测研究[J]. 地球物理学报,52(9):2264–2272. doi: 10.3969/j.issn.0001-5733.2009.09.010 Liu B J,Hu P,Meng Y Q,Feng S Y,Shi J H,Ji J F. 2009. Research on fine crustal structure using deep seismic reflection profile in Beijing region[J]. Chinese Journal of Geophysics,52(9):2264–2272 (in Chinese).
刘博研,史保平,张健. 2007. 复合地震源模拟强地面运动:以1679年三河—平谷MS8.0地震为例[J]. 地震学报,29(3):302–313. doi: 10.3321/j.issn:0253-3782.2007.03.009 Liu B Y,Shi B P,Zhang J. 2007. Strong motion simulation by the composite source modeling:A case study of 1679 M8.0 Sanhe-Pinggu earthquake[J]. Acta Seismologica Sinica,29(3):302–313 (in Chinese).
毛昌伟, 丁锐, 龚正, 张世民. 2010. 1679年三河—平谷8级地震地表断层陡坎的GPS测量[G]//地壳构造与地壳应力文集（22）. 北京: 地震出版社: 11–18. Mao C W, Ding R, Gong Z, Zhang S M. 2010. GPS survey of the surface fault scarp of 1679 Sanhe-Pinggu M8 earthquake[G]//Bulletin of the Institute of Crustal Dynamics （22）. Beijing: Seismological Press: 11–18 （in Chinese）.
潘波,许建东,关口春子,何宏林. 2006. 北京地区近断层强地震动模拟[J]. 地震地质,28(4):623–634. doi: 10.3969/j.issn.0253-4967.2006.04.010 Pan B,Xu J D,Sekigguchi H,He H L. 2006. Simulation of the near-fault strong ground motion in Beijing region[J]. Seismology and Geology,28(4):623–634 (in Chinese).
潘波,许建东,刘启方. 2009. 1679年三河-平谷8级地震近断层强地震动的有限元模拟[J]. 地震地质,31(1):69–83. doi: 10.3969/j.issn.0253-4967.2009.01.007 Pan B,Xu J D,Liu Q F. 2009. Simulations of the near-fault strong ground motion of the 1679 Sanhe-Pinggu M8 earthquake[J]. Seismology and Geology,31(1):69–83 (in Chinese).
冉勇康,邓起东,杨晓平,张晚霞,李如成,向宏发. 1997. 1679年三河—平谷8级地震发震断层的古地震及其重复间隔[J]. 地震地质,19(3):2–10. Ran Y K,Deng Q D,Yang X P,Zhang W X,Li R C,Xiang H F. 1997. Paleoearthquakes and recurrence interval on the seismogenic fault of 1679 Sanhe Pinggu M8 earthquake,Hebei and Beijing[J]. Seismology and Geology,19(3):2–10 (in Chinese).
冉志杰,李皓,吕国军,温超,杨柳,杨歧炎,孟立鹏,彭远黔,骆艳欣. 2013. 夏垫断裂夏垫段浅部构造特征地震探测[J]. 地震工程学报,35(3):656–663. doi: 10.3969/j.issn.1000-0844.2013.03.0651 Ran Z J,Li H,Lü G J,Wen C,Yang L,Yang Q Y,Meng L P,Peng Y Q,Luo Y X. 2013. Seismic detecting for the shallow tectonic features of the Xiadian fault[J]. China Earthquake Engineering Journal,35(3):656–663 (in Chinese).
沈正康,万永革,甘卫军,李铁明,曾跃华. 2004. 华北地区700年来地壳应力场演化与地震的关系研究[J]. 中国地震,(3):2–19. Shen Z K,Wan Y G,Gan W J,Li T M,Zeng Y H. 2004. Crustal stress evolution of the last 700 years in North China and earthquake occurrence[J]. Earthquake Research in China,(3):2–19 (in Chinese).
吴晶,高原,陈运泰,黄金莉. 2007. 首都圈西北部地区地壳介质地震各向异性特征初步研究[J]. 地球物理学报,50(1):209–220. doi: 10.3321/j.issn:0001-5733.2007.01.027 Wu J,Gao Y,Chen Y T,Huang J L. 2007. Seismic anisotropy in the crust in northwestern capital area of China[J]. Chinese Journal of Geophysics,50(1):209–220 (in Chinese).
向宏发,方仲景,徐杰,李如成,贾三发,郝书俭,王景钵,张晚霞. 1988. 三河—平谷8级地震区的构造背景与大震重复性研究[J]. 地震地质,10(1):15–28. Xiang H F,Fang Z J,Xu J,Li R C,Jia S F,Hao S J,Wang J B,Zhang W X. 1988. Seismotectonic background and recurrence interval of great earthquakes in 1679 Sanhe-Pinggu M=8 earthquake area[J]. Seismology and Geology,10(1):15–28 (in Chinese).
肖修来,李莹甄,许同旭,沈军. 2017. 三河—平谷地区b值和地震复发周期的时间变化分析[J]. 防灾科技学院学报,19(4):23–29. doi: 10.3969/j.issn.1673-8047.2017.04.004 Xiao X L,Li Y Z,Xu T X,Shen J. 2017. Analysis on the time variation of b-value and earthquake recurrence period in Sanhe-Pinggu area[J]. Journal of Institute of Disaster Prevention,19(4):23–29 (in Chinese).
徐锡伟,张培震,闻学泽,秦尊丽,陈桂华,朱艾斓. 2005. 川西及其邻近地区活动构造基本特征与强震复发模型[J]. 地震地质,27(3):446–461. doi: 10.3969/j.issn.0253-4967.2005.03.010 Xu X W,Zhang P Z,Wen X Z,Qin Z L,Chen G H,Zhu A L. 2005. Features of active tectonics and recurrence behaviors of strong earthquakes in the western Sichuan Province and its adjacent regions[J]. Seismology and Geology,27(3):446–461 (in Chinese).
杨晓平,曹景虎,陈献程. 2012. 夏垫活动断裂两盘岩心氧化铁变化[J]. 地震地质,34(4):659–671. doi: 10.3969/j.issn.0253-4967.2012.04.010 Yang X P,Cao J H,Chen X C. 2012. The iron oxide changes in drilling cores from the two walls of Xiadian active fault[J]. Seismology and Geology,34(4):659–671 (in Chinese).
雍凡,蒋正中,罗水余,李颜贵,刘建生,刘子龙. 2014. 夏垫断裂北段浅部构造的高分辨率地震反射勘探研究[J]. 工程地球物理学报,11(6):832–836. doi: 10.3969/j.issn.1672-7940.2014.06.016 Yong F,Jiang Z Z,Luo S Y,Li Y G,Liu J S,Liu Z L. 2014. The seismic reflection study on high-resolution profile of shallow structure in north part of Xiadian fault[J]. Chinese Journal of Engineering Geophysics,11(6):832–836 (in Chinese).
余中元,潘华,沈军,李金臣,张萌,戴训也. 2020. 夏垫断裂荣家堡探槽揭示的断裂活动特征及未来地震危险性概率评价[J]. 地震地质,42(3):688–702. doi: 10.3969/j.issn.0253-4967.2020.03.010 Yu Z Y,Pan H,Shen J,Li J C,Zhang M,Dai X Y. 2020. The activity features of Xiadian fault zone revealed by Rongjiabao trench and its probabilistic seismic hazard evaluation[J]. Seismology and Geology,42(3):688–702 (in Chinese).
张培震,张会平,郑文俊,郑德文,王伟涛,张竹琪. 2014. 东亚大陆新生代构造演化[J]. 地震地质,36(3):574–585. doi: 10.3969/j.issn.0253-4967.2014.03.003 Zhang P Z,Zhang H P,Zheng W J,Zheng D W,Wang W T,Zhang Z Q. 2014. Cenozoic tectonic evolution of continental Eastern Asia[J]. Seismology and Geology,36(3):574–585 (in Chinese).
张晚霞,向宏发,李如成. 1995. 夏垫隐伏断裂土壤气氡分布特征的初步研究[J]. 西北地震学报,17(2):46–50. Zhang W X,Xiang H F,Li R C. 1995. Preliminary study on soil-radon distribution along the Xiadian buried fault[J]. Northwestern Seismological Journal,17(2):46–50 (in Chinese).
张先康,赵金仁,刘国华,宋文荣,刘保金,赵成斌,成双喜,刘建达,顾梦林,孙振国. 2002. 三河—平谷8.0级大震区震源细结构的深地震反射探测研究[J]. 中国地震,18(4):326–336. doi: 10.3969/j.issn.1001-4683.2002.04.002 Zhang X K,Zhao J R,Liu G H,Song W R,Liu B J,Zhao C B,Cheng S X,Liu J D,Gu M L,Sun Z G. 2002. Study on fine crustal structure of the Sanhe-Pinggu earthquake (M8.0) region by deep seismic reflection profiling[J]. Earthquake Research in China,18(4):326–336 (in Chinese).
赵金仁,张先康,张成科,张建狮,杨卓欣,刘宝峰,刘保金,赵成斌. 2004. 利用宽角反射/折射和深反射探测剖面揭示三河—平谷大震区深部结构特征[J]. 地球物理学报,47(4):646–653. doi: 10.3321/j.issn:0001-5733.2004.04.015 Zhao J R,Zhang X K,Zhang C K,Zhang J S,Yang Z X,Liu B F,Liu B J,Zhao C B. 2004. Deep structural features of the Sanhe-Pinggu great earthquake area imaged by wide-angle and deep seismic reflection profiling[J]. Chinese Journal of Geophysics,47(4):646–653 (in Chinese).
周红. 2018. 基于NNSIM随机有限断层法的7.0级九寨沟地震强地面运动场重建[J]. 地球物理学报,61(5):2111–2121. doi: 10.6038/cjg2018L0651 Zhou H. 2018. Reconstruction of strong ground motion of Jiuzhaigou M7.0 earthquake based on NNSIM stochastic finite fault method[J]. Chinese Journal of Geophysics,61(5):2111–2121 (in Chinese).
周红,李亚南,常莹. 2021. 云南漾濞6.4级地震强地面运动的模拟和空间分布特征分析[J]. 地球物理学报,64(12):4526–4537. doi: 10.6038/cjg2021P0421 Zhou H,Li Y N,Chang Y. 2021. Simulation and analysis of spatial distribution characteristics of strong ground motions by the 2021 Yangbi,Yunnan Province MS6.4 earthquake[J]. Chinese Journal of Geophysics,64(12):4526–4537 (in Chinese).
Andrews D J. 1980. A stochastic fault model:1. Static case[J]. J Geophys Res,85(B7):3867–3877. doi: 10.1029/JB085iB07p03867
Beresnev I A,Atkinson G M. 1997. Modeling finite-fault radiation from the ωn spectrum[J]. Bull Seismol Soc Am,87(1):67–84.
Harris R A. 2000. Earthquake stress triggers,stress shadows,and seismic hazard[J]. Curr Sci,79(9):1215–1225.
Karakostas V,Papadimitriou E,Jin X S,Liu Z H,Paradisopoulou P,He Z. 2013. Potential of future seismogenesis in Hebei Province (NE China) due to stress interactions between strong earthquakes[J]. J Asian Earth Sci,75:1–12. doi: 10.1016/j.jseaes.2013.06.015
Li S Z,Zhao G C,Santosh M,Liu X,Dai L M,Suo Y H,Tam P Y,Song M C,Wang P C. 2012. Paleoproterozoic structural evolution of the southern segment of the Jiao-Liao-Ji belt,North China Craton[J]. Precambrian Res,200/201/202/203:59–73. doi: 10.1016/j.precamres.2012.01.007
Liu M A, Yang Y Q, Shen Z K, Wang S M, Wang M, Wan Y G. 2007. Active tectonics and intracontinental earthquakes in China: The kinematics and geodynamics[G]//Continental Intraplate Earthquakes: Science, Hazard, and Policy Issues. Boulder: Geological Society of America: 13–24.
Robinson R,Zhou S Y. 2005. Stress interactions within the Tangshan,China,earthquake sequence of 1976[J]. Bull Seismol Soc Am,95(6):2501–2505. doi: 10.1785/0120050091
Somerville P,Irikura K,Graves R,Sawada S,Wald D,Abrahamson N,Iwasaki Y,Kagawa T,Smith N,Kowada A. 1999. Characterizing crustal earthquake slip models for the prediction of strong ground motion[J]. Seismol Res Lett,70(1):59–80. doi: 10.1785/gssrl.70.1.59
Stein R S,Barka A A,Dieterich J H. 1997. Progressive failure on the North Anatolian fault since 1939 by earthquake stress triggering[J]. Geophys J Int,128(3):594–604. doi: 10.1111/j.1365-246X.1997.tb05321.x
Sun R,Vaccari F,Marrara F,Panza G F. 1998. The main features of the local geological conditions can explain the macroseismic intensity caused in Xiji-Langfu (Beijing) by the MS7.7 Tangshan 1976 earthquake[J]. Pure Appl Geophys,152(3):507–521. doi: 10.1007/s000240050164
Wang X S,Feng X D,Xu X W,Diao G L,Wan Y G,Wang L B,Ma G Q. 2014. Fault plane parameters of Sanhe-Pinggu M8 earthquake in 1679 determined using present-day small earthquakes[J]. Earthquake Science,27(6):607–614. doi: 10.1007/s11589-014-0099-3
Yu Z Y, Pan H, Xi H, Zhang Y H, Chen H. 2019. Late Quaternary paleoseismicity of the Xiadian fault in the North China Plain with implications for earthquake potential[J]. J Asian Earth Sci, 184. https://doi.org/10.1016/j.jseaes.2019.103997.
Zhang X, Wang M F. 2017. Research on the application of change detection technology based on multisource remote sensing images in the dynamic change of land use[C]//2017 3rd International Conference on Computational Systems and Communications （ICCSC 2017）. Beijing: Int. Conf. Comput. Syst. Commun: 34–40.
Zhou H, Chang Y. 2019. Stochastic finite-fault method controlled by the fault rupture process and its application to the MS7.0 Lushan Earthquake[J]. Soil Dyn Earthq Eng, 126: 1–14. https://doi.org/10.1016/j.soildyn.2019.105782.
-
期刊类型引用(2)
1. 董龙飞，李艳艳，龙海立. 联合CEEMD和改进修正能量比的高频GNSS同震位移P波震相识别方法. 全球定位系统. 2025(01): 25-32+59 . 
百度学术
2. 张环宇，林吉航，陈庭. 利用高频GPS数据反演地震要素. 全球定位系统. 2020(06): 37-45 . 百度学术
其他类型引用(1)
下载:
计量
- 文章访问数: 523
- HTML全文浏览量: 263
- PDF下载量: 158
- 被引次数: 3
