用（DDA+FEM）方法数值模拟1975年海城、1999年岫岩地震发生的过程l

用非连续变形分析方法（DDA+FEM）， 数值模拟了华北地区各地块相互制约构造环境中发生的1975年海城地震及1999年岫岩地震的过程. 通过对岫岩前震及近期余震分布图象的研究， 以及前人对海城地震的研究， 提出了海城、 岫岩地震发震构造块体模型. 数值模拟结果给出了两次地震释放的主应力变化、 最大剪应力变化等值线图， 地震前后位移变化矢量图， 及发震断层滑移随时间的变化， 分别与相应地震的震源机制、 宏观等震线、 发震断层的走滑性质等的结果基本一致.      

Numerical simulation of dynamic process of the Tangshan earthquake

In this paper, the new LDDA (Lagrangian Discontinuous Deformation Analysis) method is used in modeling the dynamic process of the MS=7.8 Tangshan earthquake on July 28, 1976 and obtain directly the dynamic and quasi-static dislocations, shear stress drops, fracture velocities of the Tangshan earthquake fault. The simulation shows that the slip history at each point of the fault is different. The displacement vectors at the concave side of the fault is greater than that of the convex side of the fault. The "over shoot" of the fault slip is greatest at the middle part of the fault and attenuates to its ends. The rupture velocities of the fault from the epicenter towards south-west and towards north-east are 3.08 m/s and 1.18 m/s, respectively, the average one is 2.13 m/s. The maximum dynamic and quasi-static dislocations are 7.1 m and 6.2 m respectively, the average quasi-static one on the fault is 4.5 m. The maximum dynamic and quasi-static shear stress drops are 8.1 MPa and 5.4 MPa, respectively, the average quasi-static shear stress drop is 3.9 MPa.We found that the rupture velocities and shear stress are related to the initial stress states of the fault.     

1997年玛尼地震对青藏川滇地区构造块体系统稳定性影响的三维DDA+FEM方法数值模拟

本文用三维非连续变形与有限元相结合（DDA+FEM）的方法,在青藏川滇地区三维构造块体相互制约的大背景中,通过用GPS资料做位移速率边界约束和震源机制约束,计算得到研究区的初始位移场和应力场与该地区GPS测量结果和震源机制分布结果基本一致.在此基础上进一步数值模拟1997年玛尼7.9级大震的发生过程,研究大震引起研究区各块体边界断层应力状态变化的特征.（1）发震断层两侧发生左旋走滑错动,最大水平位错大约7 m;（2）深部位错面上位错分布与用地震波资料震源反演的结果类似;（3）最大差应力变化等值线图与由星载D\|INSAR技术获取的地表形变场图像相似;（4）地表垂直位移表明地震断层面略向北逆冲.计算模拟得到了玛尼地震发生引起青藏川滇地区构造块体系统各边界断层上库仑破裂应力变化的分布,表明玛尼大震的发生除了使其发震断层的两端库仑破裂应力增大,应力进一步集中外,位于上地壳层上东昆仑断裂中段的2001年昆仑山8.1大震（H=11 km）发震断层段的库仑破裂应力增加约2 MPa,位于中地壳层上喀拉昆仑断裂带中的2008年改则6.9级地震（H=30 km）发震断层段的库仑破裂应力也增加约0.7 MPa,可见这两个已接近破裂强度地段的失稳对发生大震起了一定促进作用.研究结果也表明：作者发展的三维DDA+FEM方法能有效地用于大震活动与各构造块体相互作用关系的研究.     

Numerical simulations of deformation and movement of blocks within North China in response to 1976 Tangshan earthquake

Using methods of discontinuous deformation analysis and finite element (DDA+FEM), this paper simulates dynamic processes of the Tangshan earthquake of 1976, which occurred in the northern North China where its internal blocks apparently interacted. Studies focus upon both the movement and deformation of the blocks, in particular, the Ordos block, and variations of stress states on the boundary faults. The Tangshan earthquake was composed of three events: slipping motions of NNE-striking major fault, NE-striking fault near the northeastern end of the NNE-striking fault, and NW-striking fault on the southeastern side of the NNE-striking fault. Compared with previous studies, our model yields a result that is more agreeable with the configuration of aftershock distributions. A number of data are presented, such as the principle stress field during the earthquake, contours of the maximum shear stress, the strike-slip deformation between blocks near the earthquake focus, time-dependent variations of slips of earthquake-triggered faulting, the maximum slip distance, and stress drops. These results are in accord with the earthquake source mechanism, basic parameters from earthquake wave study, macro-isoseismic line, observed horizontal displacement vectors, etc. The Tangshan earthquake exerted different influences on the adjacent blocks and boundary faults between them, thus resulting in differential movement and deformation. The Ordos block seems to have experienced the small-scale counterclockwise rotation and deformation, but its northeast part, bounded on the east by the Taihangshan and on the north by the Yanshan and Yinshan belts, underwent relatively stronger deformation. The Tangshan earthquake also changed the stress state of boundary faults of the North China, leading to an increase in shear stress and a decrease in normal stress in the NW-trending Zhangjiakou-Penglai fault through Tangshan City and the northern border faults of the Ordos block, and therefore raises the potential risk of earthquake occurrence. This result is supported by the facts that a series of Ms≥ 6 earthquakes took place at the northern margin of the Ordos block after the Tangshan earthquake.     

液化场地地震反应的离散元数值模拟

采用一个描述饱和砂土中固相颗粒和液相流体耦合的离散元数值模型,固相颗粒采用离散元的颗粒流理论(PFC)模拟,液相流体通过求解平均Navier-Stokes方程的计算流体动力学(CFD)技术计算.采用上述模犁数值模拟了在地震荷载作用液化场地的加速度、剪应力、剪应变时程以及超孔隙水压力和土层表面位移的发展过程.数值模拟结果...     

Database application platform for earthquake numerical simulation

Introduction In recent years, all kinds of observation networks of seismology have been established, which have been continuously producing numerous digital information. In addition, there are many study results about 3D velocity structure model and tectonic model of crust （Huang and Zhao, 2006; Huang et al, 2003; Li and Mooney, 1998）, which are valuable for studying the inner structure of the earth and earthquake preparation process. It is badly needed to combine the observed data, experimental study and theoretical analyses results by the way of numerical simulation and develop a database and a corresponding application platform to be used by numerical simulation, and is also a significant way to promote earthquake prediction.     

拐角频率在营海岫地区地震序列判定中的应用

利用拐角频率方法,研究营、海、岫地区1999年以来发生的显著震群性质,结果表明,前震序列与其他震群拐角频率曲线变化形态不同,拐角频率的大小能较好反映震源区应力场的大小。     

地震次生毒气泄漏与扩散数值模拟的参数分析

本文在破坏性地震中可能发生的次生毒气泄漏与扩散灾害的数值模拟与动态仿真研究的基础上，分析了影响毒气扩散的主要参数。在设定的几种工况下，给出相应的图表结果，直观地描述了各参数对扩散过程的影响。进一步验证了所采用毒气扩散简化模型的准确性，对大震应急和防御与减轻城市和大型企业的地震次生毒气泄漏与扩散灾害都有重要的意义。     

Numerical simulation of precursors of Gonghe M 7.0 earthquake in Qinghai Province

Based on a simplified media model of Gonghe area, the precursor characteristics of Gonghe M7.0 earthquake in 1990 are simulated in this paper by using the constitutive relationship of binary medium (solid and water). The results show that the simulated state distribution and extension variation of media are identical with the spatial and temporal distribution of reliable anomalies before Gonghe earthquake. The study also suggests that the development of the Gonghe earthquake has experienced a series of processes such as elastic deformation of large scale, early nonelastic dilatation, strain softening and elastic recovery in the neighbour region of Xining, nonelastic dilatation of high dense block near seismic source and earthquake occurrence. Therefore, it can be concluded that this earthquake development is neither a simple process for fissures developing and linking up with one another, nor a process of nonelastic volume expansion and water flowing into medium around seismic source, contrary, it is a complicated medium state changing process, and the style of such changing depends on geological structure environment of seismic source and its neighbour region. It is considered that different earthquakes are accompanied by different geological conditions, the spatial and temporal behavior of their precursors are certainly different.     

震源动力学破裂过程数值模拟研究

首先阐述了震源动力学过程研究的重要意义, 在此基础上, 研究了复杂的断层几何形态及介质模型对动力学破裂过程的影响, 并对常用的有限元方法、 离散元方法、 有限差分方法和边界积分方程方法等进行了相应介绍. 讨论了这些数值模拟方法各自的优缺点, 建议在方法的选择上应视具体问题及计算的精度而定. 最后对动力学数值模拟的关键部分, 滑动摩擦准则进行了论述. 常用的滑动摩擦准则有滑动弱化准则、 速率弱化准则和速率-状态依赖摩擦准则. 在单纯考虑某个地震的动力学破裂传播过程时, 滑动弱化准则较为常用, 其中滑动弱化距离的选取至关重要. 但若考虑整个地震循环, 速率-状态依赖摩擦准则更为合适.      

地震次生毒气泄漏与扩散的数值模拟与动态仿真

文章对破坏性地震中可能发生的次生毒气泄漏与扩散灾害进行了数值模拟与动态仿真研究，设定了几种有实际意义的数值分析模型，分别进行了结构地震破坏，有毒有害气体泄漏，在一定泄漏和气象条件下的向周围空间扩散的动态仿真研究，给出了典型毒气扩散过程中可致人员死亡，严重危险和有感区域的范围的动态变化结果，作者认为地震次生毒气泄漏与扩散研究对大地震应急、城市和大型企业防灾都有重要的意义和广阔的应用前景。     

Current state of integrated earthquake simulation for earthquake hazard and disaster

This paper presents the current state of integrated simulation for earthquake hazard and disaster. This simulation takes advantage of the macro–micro analysis method; this method estimates an earthquake’s strong motion with high spatial and temporal resolution, using the bounding medium theory to obtain optimistic and pessimistic estimates of expected strong motion distribution and the singular perturbation expansion that results in an efficient multi-scale analysis. Integrated earthquake simulation calculates seismic responses for all structures in a target area, inputting simulated strong ground motion to a structure analysis method that is plugged into the system by means of a wrapper; a suitable method, linear or nonlinear, is chosen depending on the type of the structure. The results of all simulations are visualized so that residences and government officials can share a common recognition of earthquake hazard and disaster. Two examples of this integrated earthquake simulations are presented; one is made by plugging nonlinear structure analysis methods into the system, and the other is made for an actual city, the computer model of which is constructed with the help of available geographical information systems.     

冲绳海槽地震海啸的数值试验研究

在分析东海水深与地震情况基础上,我们认为冲绳海槽地震具有引发海啸的条件,并对该区海啸传播过程采用单侧破裂方式的有限移动源模型和高阶Boussinesq方程做了一维有限差分法数值模拟,给出了海啸传播过程波形和速度变化图.我们从模拟结果得到初步结论:冲绳海槽特大地震海啸传播到上海近海在5小时后,近海波高在几十厘米范围内.     

Numerical simulation on the influences of Wenchuan earthquake on the surrounding faults

On 12 May 2008, the devastating Wenchuan earthquake struck the Longmenshan fault zone, which comprised the eastern margin of the Tibetan Plateau, and this fault zone was predominantly a convergent boundary with a right-lateral strike-slip component. After such a large-magnitude earthquake, it was crucial to analyze the influences of the earthquake on the surrounding faults and the potential seismic activity. In this paper, a complex viscoelastic model of western Sichuan and eastern Tibet regions was constructed including the topography. Based on the findings of co-seismic static slip distribution, we calculated the stress change caused by the Wenchuan earthquake with the post-seismic relaxation into consideration. Our preliminary results indicated that: (1) The tectonic stressing rate was relatively high in Kunlun mountain pass-Jiangcuo, Ganzi-Yushu, Xianshuihe and Zemuhe faults; while in the east Kunlun and Longriba was medium; also the value was less in the Minjiang, Longmenshan, Anninghe and Huya faults. As to the Longmenshan fault, the value was 0.28×10^-3 MPa/a to 0.35×10^-3 MPa/a, which is coincident with the previous long recurrence interval of Wenchuan earthquake; (2) The Wenchuan earthquake not only caused the Coulomb stress decrease in the source region, but also the stress increase in the two terminals, especially the northeastern segment, which is comparatively consistent with the aftershock distribution. Meanwhile, the high concentration areas of the static slip distribution were corresponding to the Coulomb stress reductions; (3) The Coulomb stress change caused by Wenchuan earthquake showed significant increase on five major faults, which were northwestern segment of Xianshuihe fault, eastern Kunlun fault, Longriba fault, Minjiang fault and Huya fault respectively; also the Coulomb stress on the fault plane of the Yushu earthquake was faintly increased; (4) We defined the recurrence interval as the time needed to accumulate the magnitude of the stress drop, and the recurrence interval of Wenchuan earthquake was estimated about 1 714 a to 2 143 a correspondingly.     

Role of fluids in surface deformation caused by the 1999 Chi‐Chi earthquake in Taiwan

The 1999 Chi‐Chi earthquake significantly altered the landscape of central Taiwan. Surface deformation produced by the earthquake along the trace of the Chelungpu thrust can be classified into two styles: (1) uplift without significant surface rupture, and (2) uplift accompanied by surface rupture. Here we examine areas that exhibited the first style of deformation (e.g. Wufeng). Seismic stress at the time of the main shock may have been relieved by high pore‐fluid pressure in a 300‐m‐thick sand and gravel aquifer. Along the thrust fault, frictional heating of these sediments resulted in thermal expansion and an increase in pore‐fluid pressure. High pore‐fluid pressure damped seismic‐wave energy and enhanced intergranular slips of unconsolidated sandy and gravel sediments, which were possibly assisted by sulphuric acid corrosion, leading to a high sulphate content in the groundwater (c. 70 mg L^?1). These changes permitted surface folding and terrace‐style uplifting to occur without significant rupture. In contrast, other areas in which the second style of deformation is dominant (e.g. Fengyuen‐Shihkang) have thin (0–10 m) sand and gravel deposits and lower concentrations of sulphate (c. 30 mg L^?1) in groundwater. In these areas, sediments were heated but not sufficiently to produce significant thermal expansion and increase in pore‐fluid pressure; accumulation of stress in these locations led to rupture at the ground surface, with the formation of steep fault scarps. The areas exhibiting the first deformation style are characterized by the presence of high pore‐fluid pressure, frictional heat conduction, and possibly chemical corrosion related to sulphuric acid attack and formation of sulphate, in contrast to those involving significant uplift and surface rupture. The areal distribution of these two surface deformation styles suggests that the aforementioned fluid‐related subsurface processes may have altered the characteristics of sediments and caused diverse responses to the quake. Copyright © 2002 John Wiley & Sons, Ltd.     

Deformation patterns associated with the M5.9 Athens (Greece) earthquake of 7 September 1999

The static displacement field of the Athens 1999 earthquake has been numerically modeled by a BEM method and analysed from SAR interferometry images with compatible results: (a) for a fault that reaches the surface the subsidence field coincides with the hangingwall domain of the Fili neotectonic normal fault with maximum amplitude, d _max, 5.5–7 cm, which is consistent with the possibly co-seismic displacement of 6–10 cm observed in the field, the average fault dislocation of 5–8 cm found by the application of circular source models, and the displacement up to 6 cm predicted by empirical relations between magnitude and displacement; the field of uplift covers the footwall domain of the fault with d _max1.5 cm;d gradually decreases with distance from the fault at a gradient of 0.4 cm/km, (b) for a blind fault d _max is only 1.8 and 0.3 cm in the hangingwall and footwall, respectively, and the decay gradient becomes 0.15 cm/km, (c) the total deformation area is 15 km × 15 km and the Fili fault, with a preferred mean dip of 60°, constitutes the natural boundary between the subsidence and uplift areas. The macroseismic field pattern is similar with that of the static ground deformation. The majority of intensity values VI (MM and EMS-98 scales), are distributed within the hangingwall of the Fili fault, while the highest intensities (VIII and IX) concentrate very close to the Fili fault within its hangingwall domain. A gradual decrease of the intensities with the distance from the Fili fault is evident. Because of the similarity between the intensity distribution pattern and the static ground deformation pattern, we make the hypothesis that the latter predicts well enough the main characteristics of the former although the ground displacement is dominated by relatively low frequency as compared to the ground acceleration.     

强震孕育、发生及其复发循环过程的有限单元法模拟

强震孕育发生及其复发过程的定量研究对于预测预报地震有着重要的科学意义;合理连续地计算地震孕育、同震破裂过程及其复发循环特征将有助于我们更好地认识强震发生的时空分布规律.为此,本文基于Newmark隐式时间积分法,根据模拟孕震-同震循环过程的特殊要求,发展了一种新的有限单元计算方法.新方法具有以下特点：（1）在不改变时间积分方法的情况下,实现对时间步长进行自动平滑地缩放,进而可以连续的模拟准静态、动态不同力学状态下的孕震-同震循环过程;（2）模型的初始应力场非人为指定,计算时通过施加重力及缓慢的构造加载获得;（3）地震破裂的成核区域及成核方式等亦非人为给定,模拟中破裂成核是自然形成的,这样更加符合实际地震地质情况.通过大量的模型计算,其结果表明,新的计算方法可以连续稳定的模拟断层孕震-同震及其循环过程,计算结果不仅可以给出强震的复发间隔,同时还可以给出地震时断层破裂行为的详细过程.此外,本文还考察了摩擦系数对地震准周期性的影响,发现断层上的静、动摩擦系数差值直接影响着断层的强震复发周期,两者差值越小,复发周期越短;差值越大,复发间隔越长.

     

Studies and experiments on earthquake prediction during 1999～2002

This paper briefly reviewed the research progresses of earthquake prediction and/or forecasting in China during 1999～2002, especially focused on mid-short term prediction methods with approaches of seismicity, crustal de-formation, electromagnetism, ground water and the analysis by synthesis, and the application of the methods to the practice of earthquake prediction.     

云南武定M6.5地震动态库仑破裂应力变化的数值模拟及其与余震活动的关系

在离散波数法（DWN）基础上,计算了武定M6.5地震断层破裂在周围介质中产生的位移时程（位移理论地震图）和动态位移场;进行弹性动力学转换后,求得应变时程和动态应变场;最后得到了武定M6.5地震所产生的动态库仑破裂应力变化量和动态库仑破裂应力变化场,进而研究其与后续余震的关系. 结果表明,动态应力最大峰值和静态应力的正区均呈非对称性分布,两者的分布特征与余震的分布特征基本一致. 在动态应力峰值为正的确定区应力值超过了0.1 MPa触发阈值, 在静态应力值为正的确定区应力值超过了0.01 MPa触发阈值. 这说明动态应力和静态应力均有助于余震的发生.      

Numerical simulation of dynamic Coulomb stress changes induced by M6.5 earthquake in Wuding, Yunnan and its relationship with aftershocks

Based on the discrete wavenumber method, we calculate the fields of dynamic Coulomb rupture stress changes and static stress changes caused by M6.5 earthquake in Wuding, and study their relationship with the subsequent after- shocks. The results show that the spatial distribution patterns of the positive region of dynamic stress peak value and static stress peak value are similarly asymmetric, which are basically identical with distribution features of aftershock. The dynamic stress peak value and the stati...