安新—宽城测深剖面的时域有限差分数值模拟

地震 ›› 2008, Vol. 28 ›› Issue (2): 54-64.

安新—宽城测深剖面的时域有限差分数值模拟

彭菲¹, 陈棋福¹, 刘澜波^1,2, 陈颙¹

1.中国地震局地震预测研究所, 北京 100036;
2.康涅狄格大学, 美国

收稿日期:2007-10-23 修回日期:2007-12-15 出版日期:2008-04-30 发布日期:2021-10-29
作者简介:彭菲(1981-), 女, 河北昌黎人, 硕士, 主要从事地震波场数值模拟等研究。
基金资助:
国际科技合作项目(2006DFA21650)资助

The finite difference time domain simulation on Anxin- Kuancheng deep seismic sounding profile

PENG Fei¹, CHEN Qi-fu¹, LIU Lan-bo^1,2, CHEN Yong¹

1. Institute of Earthquake Science, CEA, Beijing 100036, China;
2. Department of Civil and Environmental Engineering, University of Connecticut, Storrs, CT 06269, USA

Received:2007-10-23 Revised:2007-12-15 Online:2008-04-30 Published:2021-10-29

摘要/Abstract

摘要： 利用时域有限差分方法进行地震波模拟能够得到丰富的波场信息, 对研究复杂地质构造条件下的地震波波场十分有效。该文在对时域有限差分方法进行概述的基础上, 以2002年4月施测的安新—宽城地震测深剖面资料为基础, 采用时域有限差分方法对该剖面进行二维全波数值模拟, 同时以理想匹配层作为模拟的吸收边界条件。将模拟得到的合成地震图, 通过与野外实验得到的记录截面和射线追踪合成记录相比较和讨论, 验证了时域有限差分方法在模拟地震波传播过程的有效性。

关键词: 地震波场, 二维模拟, 时域有限差分, 吸收边界

Abstract: Using the method of finite difference time domain to simulate the seismic wave can obtain much more information of the wave field. It is very effective to investigate the seismic wave field on complex geological background. In this paper, the finite difference time domain (FDTD) was briefly introduced. For the Anxin-Kuancheng deep seismic sounding profile conducted on April 2002, the two-dimensional full wave numerical simulation was performed with FDTD method, as well as the perfectly matched layer as the absorbing boundary condition. Discussing and comparing the FDTD simulation synthetic seismograms with the recorded sections and the synthetic results given by ray tracing method, the effectiveness of the finite difference time domain method for the simulation of the seismic wave propagation has been proved.

Key words: Seismic wave propagation, 2-D simulation, Finite difference time domain, Absorbing boundary condition

中图分类号:

P315.7

彭菲, 陈棋福, 刘澜波, 陈颙. 安新—宽城测深剖面的时域有限差分数值模拟[J]. 地震, 2008, 28(2): 54-64.

PENG Fei, CHEN Qi-fu, LIU Lan-bo, CHEN Yong. The finite difference time domain simulation on Anxin- Kuancheng deep seismic sounding profile[J]. EARTHQUAKE, 2008, 28(2): 54-64.

参考文献

[1] 张永刚. 地震波场数值模拟方法[J]. 石油物探, 2000, 42(2): 143-148.
[2] Yee K S. Numerical solution of initial boundary value problems involving Maxwell equations in isotropic media[J]. IEEE Trans Antennas Propagat, 1966, 14(3): 302-307.
[3] Kelly K R, Ward R W, Alford R M. Synthetic seismograms: a finite-difference approach[J]. Geophysics, 1946, 41(1): 2-27.
[4] Berenger J P. A perfectly matcher layer for the absorption of electromagnetic waves[J]. J Comput Phys, 1994, 114(2): 185-200.
[5] Abramovici F, Alterman Z. Computations pertaining to the problem of propagation of a seismic pulse in a layered solid[J]. Methods of Computational Physics, 1965, Vol. IV, Academic Press: 349-379.
[6] Alterman Z, Karal F C. Propagation of elastic waves in layered media by finite-difference methods[J]. Bull Seism Soc Am, 1986, 58: 367-398.
[7] Madariaga R. Dynamics of an expanding circular fault[J]. Bull Seism Soc Am, 1976, 66(3): 639-666.
[8] Kaser M, Igel H. Numerical simulation of 2D wave propagation on unstructured grids using explicit differential operators[J]. Geophysical Prospecting, 2001, 49: 607-619.
[9] Lanbo Liu, Steve A. Arcone, Numerical simulation of the wave-guide effect of the near-surface thin layer on radar wave propagation[J]. JEEG, 2003, 8(2): 53-61.
[10] 田小波, 吴庆举, 曾融生. 弹性波场数值模拟的隐式差分多重网格算法[J]. 地球物理学报, 2004, 47(1): 81-87.
[11] 殷文, 印兴耀, 等. 高精度频率域弹性波方程有限差分方法及波场模拟[J]. 地球物理学报, 2006, 49(2): 562-568.
[12] 王雪秋. 利用有限积分变换—有限差分方法模拟地震波场[D]. 吉林大学硕士学位论文, 2002.
[13] 张文生. 科学计算中的偏微分方程有限差分法[M]. 北京: 高等教育出版社, 2006. 328-331.
[14] 郑宏兴, 黄文武, 谢洪波, 等. 弹性波时域仿真的两个问题[J]. 天津大学学报, 2005, 38(4): 338-342.
[15] 王夫运, 张先康, 陈棋福, 等. 张渤构造带中东段二维P波速度结构安新—香河—宽城剖面[J]. 地震学报, 2004, 26(增刊): 32-40.