层状介质中有限震源引起的地面运动的计算对点源情形的拓展

邸海滨　许力生

摘要: 从分层均匀介质中地震波传播的基本理论出发，参考已有的利用广义反射透射系数矩阵方法计算分层均匀介质中点源引起的地面运动的方法和程序，考虑了任意几何特征的有限震源及有限震源的震源机制随时间和空间发生变化的可能性，并考虑了并行计算的发展方向, 本文对点源情况下的计算流程进行了改进，重新设计编写了计算程序GRTMatSyn，使其不但适应于点源的情形，也适应于任意几何形状的、有限的、震源机制随时间和空间变化的有限震源的情形．为了验证该程序的可靠性和计算效率，设计了必要的有限震源模型和观测点分布，计算了地面运动，并与已有的被广泛认可的计算有限平面断层引起的地面运动的程序CompSyn的计算结果进行了对比分析．结果表明， GRTMatSyn的计算结果可靠、计算效率较高．

关键词:

Abstract: From fundamental theory of seismic wave propagation in stratified media, with reference to the techniques and applications already developed for calculating point-source-caused ground motion using generalized reflection and transmission matrices method, and with consideration of arbitrary geometry and tempo-spatially variable focal mechanisms of finite sources, as well as development trend of parallel computation, we improved computing-flow in the point-source case and developed a new program, GRTMatSyn, to be suitable for both point sources and finite sources with arbitrary geometries and tempo-spatially variable mechanisms. Numerical tests were carried out with necessary models of finite sources and observation points for reliability and computing efficiency of the application. The results were compared with those from CompSyn, a widely recognized program for calculating ground motion caused by a finite plane fault. The comparison indicated that the application developed in this study is reliable and rather more efficient.

邸海滨. 2011. 层状介质中有限震源引起的地面运动计算［D］. 北京：中国地震局地球物理研究所： 51——90.

李旭. 1993. 用地震波形资料反演1990年青海共和地震的震源过程［D］. 北京: 国家地震局地球物理研究所： 1——45.

李旭，陈运泰. 1996. 合成地震图的广义反射透射系数矩阵方法［J］. 地震地磁观测与研究， 17（2）： 1——20.

刘超. 2011. 断层厚度的地震效应和非对称矩张量［D］. 北京：中国地震局地球物理研究所： 1——110.

姚振兴, 郑天愉. 1985. 对二滩水电站坝区场地地面运动的估计［J］. 地球物理学报， 28(2): 170——180.

张勇. 2008. 震源破裂过程反演方法研究［D］. 北京：北京大学： 1——16.

张勇, 冯万鹏, 许力生, 周成虎, 陈运泰. 2008. 2008年汶川大地震的时空破裂过程［J］. 中国科学: D辑， 38(10): 1186——1194.

张勇, 许力生, 陈运泰. 2010. 2010年青海玉树地震震源破裂过程［J］. 中国科学： D辑， 40(7): 819——821.

Bouchon M. 1981. A simple method to calculate Green′s functions for elastic layered media［J］. Bull Seism Soc Amer，71（4）: 957——971.

Chapman C H. 1974. Generalized ray theory for an inhomogeneous medium［J］. Geophys J R astr Soc， 36(3): 673——704.

Chen X F. 1999. Seismogram synthesis in multi——layered half——space part Ⅰ: Theoretical formulations［J］. Earthquake Research in China, 13(2): 149——174.

Gilbert F, Helmberger D V. 1972. Generalized ray theory for a layered sphere［J］. Geophys J Inter，27(1): 57——80.

Helmberger D V. 1968. The crust——mantle transition in the Bering Sea［J］. Bull Seism Soc Amer， 58(1): 179——214.

Kennett B L N, Kerry N J. 1979. Seismic waves in a stratified half space［J］. Geophys J R astr Soc，58: 179——214.

Kennett B L N. 1980. Seismic waves in a stratified space Ⅱ: Theoretical seismograms［J］. Geophys J R astr Soc，61(1): 1——10.

Kennett B L N. 1983. Seismic Wave Propagation in a Stratified Media［M］. Cambridge: Cambridge University Press: 1——342.

Olson A H, Orcutt J A, Frazier G A. 1984. The discrete wave number/finite element method for synthetic seismograms［J］. Geophys J R astr Soc，77(2): 421——460.

Spudich P, Archuleta R. 1987. Techniques for earthquake ground motion calculation with applications to source parameterization of finite faults［G］∥Bolt B A ed. Seismic Strong Motion Synthetics. Orlando, Florida: Academic Press: 205——265.

Spudich P, Xu L S. 2003. Software for calculating earthquake ground motion from finite faults in vertically varying media［G］∥Lee W H K, Kanamori H, Jennings P C, Kisslinger C eds. International Handbook of Earthquake and Engineering Seismology: 1633——1634.

Yao Z X, Harkrider D G. 1983. A generalized reflection transmission coefficient matrix and discrete wave number method for synthetic seismograms［J］. Bull Seism Soc Amer，73(6A): 1685——1699.