计算断层自发破裂过程的离散网格划分研究

英文篇名：Study on the meshing scheme in the computation of spontaneous fault rupture process
中文刊名：地震学报
英文刊名：Acta Seismologica Sinica
作者：刘启明 ; 陈晓非
英文作者：Liu Qiming1) Chen Xiaofei2) ; 1) School of Earth and Space Science ; Peking University ; Beijing 100871 ; China 2) School of Earth and Space Science ; University of Science and Technology of China ; Hefei 230026 ; China
中文关键词：震源动力学 ; 有限震源破裂 ; 边界积分方程 ; 网格划分
英文关键词：seismic source dynamics ; finite fault source rupture ; boundary integral equation ; meshing scheme
出版日期：2008-10-15
机构：北京大学地球与空间科学学院;中国科学技术大学地球和空间科学学院;
年：2008
期：05
出版单位：地震学报

摘要

离散化网格的空间步长选取在各种数值算法中都是一个很受关注的问题,在全空间均匀介质模型和简单离散化方案下,利用边界积分方程方法研究在自发破裂求解过程中动力学参数组合Dc和Te(Dc为临界滑动弱化位移,Te为有效的断层破裂面上的初始应力)对计算网格划分的影响,初步得到了Dc和Te参数空间中的有效计算网格的选取规律,对合理、有效地运用边界积分方程方法计算地震震源的破裂过程具有重要的指导意义.
The choice of spatial grid size has been being a crucial issue in all kinds of numerical algorithms. As using BIEM(Boundary Integral Equation Method) to calculate the rupture process of a planar fault embedded in an isotropic and homogeneous full space with simple discretization scheme, this paper focuses on what grid size should be applied to control the error as well as to keep the computing efficiency for different parameter combinations of (Dc ,Te ), where Dc is the critical slip-weakening distance and Te is the initial stress on the fault plane. We have preliminarily found the way of properly choosing the spatial grid size, and this is of great significance in the computation of seismic source rupture process with BIEM.

引文

Aagaard B,Heaton T,Hall J.2001.Dynamic earthquake rupturesinthe presence of lithostatic normal stresses:I mplica-tions for friction models and heat production[J].Bull Seism Soc Amer,91:1765--1796.
    Aki K,Richards P.2002.Quantitative Seismology:Theory and Methods[M].Second edition.Sausalito,California:University Science Books:38--60.
    Andrews D.1976.Rupture propagation with finite stress in antiplane strain[J].J Geophys Res,81:3575--3582.
    Aochi H,Fukuyama E,Matsu’ura M.2000.Spontaneous rupture propagation on a non-planar fault in3-Delastic medi-um[J].Pure Appl Geophys,157:2003--2027.
    Fukuyama E,Madariaga R.1995.Integral equation method for plane crack with arbitrary shape in3D elastic medium[J].Bull Seism Soc Amer,85:614--628.
    Ida Y.1972.Cohesive force across the tip of a longitudinal-shear crack and Griffith’s specific surface energy[J].J Geo-phys Res,77:3796--3805.
    Madariaga R.1976.Dynamic of an expanding circular fault[J].Bull Seism Soc Amer,66:639--667.
    Madariaga R,Olsen K,Archuleta R.1998.Modeling dynamic rupture in a3D earthquake fault model[J].Bull SeismSoc Amer,88:1182--1197.
    Oglesby D,Archuleta R,Nielson S.2000.The three-di mensional dynamics of dipping faults[J].Bull Seism Soc Amer,90:616--628.
    Zhang H M,Chen X F.2006.Dynamic rupture on a planar fault in three-di mensional half space:Ⅰ.Theory[J].Geo-phys J Int,164(3):633--652.