A Monte Carlo adapted finite element method for dislocation simulation of faults with uncertain geometry

Dislocation modelling of an earthquake fault is of great importance due to the fact that ground surface response may be predicted by the model. However, geological features of a fault cannot be measured exactly, and therefore these features and data involve uncertainties. This paper presents a Monte Carlo based random model of faults with finite element method incorporating split node technique to impose the effects of discontinuities. Length and orientation of the fault are selected as random parameters in the domain model, and hence geometrical uncertainties are encountered. Mean and standard deviation values, as well as probability density function of ground surface responses due to the dislocation are computed. Based on analytical and numerical calculation of dislocation, two approaches of Monte Carlo simulations are proposed. Various comparisons are examined to illustrate the capability of both methods for random simulation of faults.

[1]  S. Day,et al.  Comparison of Fault Representation Methods in Finite Difference Simulations of Dynamic Rupture , 2006 .

[2]  Vito Volterra,et al.  Sur l'équilibre des corps élastiques multiplement connexes , 1907 .

[3]  H. Melosh,et al.  Anelastic response of the earth to a dip slip earthquake , 1983 .

[4]  Y. Kagan.  Stochastic model of earthquake fault geometry , 1982 .

[5]  Sumio Sawada,et al.  Simulation of Fault Rupture Process by the Stochastic Finite Element Method , 1987 .

[6]  Luis A. Dalguer,et al.  Staggered-grid split-node method for spontaneous rupture simulation , 2007 .

[7]  Salvatore Barba,et al.  Normal-fault stress and displacement through finite-element analysis , 2005 .

[8]  A. V. Mikhailova,et al.  Finite element modelling of elastic intraplate stresses due to heterogeneities in crustal density and mechanical properties for the Jabalpur earthquake region, central India , 2008 .

[9]  Arthur Raefsky,et al.  A simple and efficient method for introducing faults into finite element computations , 1981 .

[10]  Muneo Hori,et al.  Introduction To Computational Earthquake Engineering , 2006 .

[11]  Static deformation of a uniform half-space due to a long dip-slip fault , 1992 .

[12]  D. J. Andrews,et al.  Test of two methods for faulting in finite-difference calculations , 1999 .

[13]  Bala Neeru,et al.  Deformation of a two-phase medium due to a long buried strike-slip fault , 2013 .

[14]  H. Melosh,et al.  Mechanics of graben formation in crustal rocks - A finite element analysis , 1989 .

[15]  G. A. Frazier,et al.  Finite element analysis of the residual displacements for an earthquake rupture: Source parameters for the San Fernando earthquake , 1973 .

[16]  Tsuyoshi Ichimura,et al.  ANALYSIS METHODS OF STOCHASTIC MODEL: APPLICATION TO STRONG MOTION AND FAULT PROBLEMS , 2003 .

[17]  E. H. van Brummelen,et al.  Discontinuities without discontinuity : the Weakly-enforced Slip Method , 2013, 1304.7810.

[18]  Y. Okada.  Internal deformation due to shear and tensile faults in a half-space , 1992, Bulletin of the Seismological Society of America.

[19]  R. Hanssen,et al.  Overview of a range of solution methods for elastic dislocation problems in geophysics , 2013 .

[20]  G. Brietzke,et al.  Comparison of deterministic and stochastic earthquake simulators for fault interactions in the Lower Rhine Embayment, Germany , 2013 .