The numerical simulation of thesource development of the Tangshan earthquake

Ｔｈｅｎｕｍｅｒｉｃａｌｓｉｍｕｌａｔｉｏｎｏｆｔｈｅｓｏｕｒｃｅｄｅｖｅｌ┐ｏｐｍｅｎｔｏｆｔｈｅＴａｎｇｓｈａｎｅａｒｔｈｑｕａｋｅＺＨＩ－ＲＥＮＮＩＵ（牛志仁）ＤＡＮＧ－ＭＩＮＣＨＥＮ（陈党民）ＢＩＮＧ－ＱＩＡＮＬＩ（李炳乾）Ｓｅｉｓｍｏｌｏ...     

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

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

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.     

数值模拟方法在应力场演化及地震科学中的研究进展

本文对数值模拟方法(主要针对有限单元法)在应力场演化及地震科学中的应用和发展做了回顾与总结,主要包括构造应力场模拟、地震活动性模拟、地震危险性模拟等部分,并对数值模拟方法的发展趋势、方法的优势及其所存在的主要问题做了探讨.     

Mechanical features of Shanxi Datong earthquake swarm sequence and their numerical simulation by using finite element method

The seismicity characteristics of the earthquake swarm occurred in Datong area in Shanxi Province during 1989–1993 period are studied in detail. The focal mechanisms and rupture processes of the principal earthquakes in it are also estimated in this paper. A mechanical model of fracture system of asymmetrical conjugate array is suggested from the view point of fracture propagation in order to simulate the mechanical features in the seismic source respect of the Datong swarm sequence. The computation results of displacement field, stress field and strain energy distribution from this model by using finite element method are presented. In addition, some recorded data of geothermal precursors before these three climaxes of the Datong earthquake swarm sequence are discussed mainly in this paper. Contribution No. 94A0093, Institute of Geophysics, SSB, China.     

识别震源性质的一种新方法

提出了利用地震事件的震源差异在倒谱域所呈现的不同特征识别震源性质的新方法。识别结果表明新方法是有效的 ,展示出良好的应用前景     

Mechanical model of precursory source proces for strong earthquakes in western Yun－nanarea——conceptandnumericalsimula┐tion

Ｍｅｃｈａｎｉｃａｌｍｏｄｅｌｏｆｐｒｅｃｕｒｓｏｒｙｓｏｕｒｃｅｐｒｏ┐ｃｅｓｆｏｒｓｔｒｏｎｇｅａｒｔｈｑｕａｋｅｓｉｎｗｅｓｔｅｒｎＹｕｎ┐ｎａｎａｒｅａ——ｃｏｎｃｅｐｔａｎｄｎｕｍｅｒｉｃａｌｓｉｍｕｌａ┐ｔｉｏｎＹＯＵ－ＪＩＮＳＵ（苏有...     

Strong motion simulation by the composite source modeling： A case study of 1679 M8.0 Sanhe-Pinggu earthquake

In this study, a composite source model has been used to calculate the realistic strong ground motions in Beijing area, caused by 1679 MS8.0 earthquake in Sanhe-Pinggu. The results could provide us the useful physical parame-ters for the future seismic hazard analysis in this area. Considering the regional geological/geophysical background, we simulated the scenario earthquake with an associated ground motions in the area ranging from 39.3°N to 41.1°N in latitude and from 115.35°E to 117.55°E in longitude. Some of the key factors which could influence the characteristics of strong ground motion have been discussed, and the resultant peak ground acceleration (PGA) distribution and the peak ground velocity (PGV) distribution around Beijing area also have been made as well. A comparison of the simulated result with the results derived from the attenuation relation has been made, and a suf-ficient discussion about the advantages and disadvantages of composite source model also has been given in this study. The numerical results, such as the PGA, PGV, peak ground displacement (PGD), and the three-component time-histories developed for Beijing area, have a potential application in earthquake engineering field and building code design, especially for the evaluation of critical constructions, government decision making and the seismic hazard assessment by financial/insurance companies.     

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.     

Comparative analysis of the characteristics of global seismic wave propagation excited by the Mw9.0 Tohoku earthquake using numerical simulation

Giant earthquakes generate rich signals that can be used to explore the characteristics of the hierarchical structure of the Earth’s interior associated with the eigenfrequencies of the Earth.We employ the spectral element method,incorporated with large-scale parallel computing technology,to investigate the characteristics of global seismic wave propagation excited by the2011 Mw9.0 Tohoku earthquake.The transversely isotropic PREM model is employed as a prototype of our numerical global Earth model.Topographic data and the effect of the oceans are taken into consideration.Wave propagation processes are simulated by solving three-dimensional elastic wave governing equations with the seismic moment tensor obtained from the Global Centroid Moment Tensor Catalog.Three-dimensional visualization of our computing results displays the nature of the global seismic wave propagation.Comparative analysis of our calculations with observations obtained from the Incorporated Research Institutions for Seismology demonstrates the reliability and feasibility of our numerical results.We compare synthetic seismograms with incorporated and unincorporated ocean models.First results show that the oceans have obvious effects on the characteristics of seismic wave propagation.The peak displacement and peak velocity of P waves become relatively small under the effect of the ocean.However,the effect of the ocean on S-waves is complex.The displacement and velocity of S waves decrease rapidly over time using an unincorporated ocean model.Therefore,the effects of the ocean should be incorporated when undertaking quantitative earthquake hazard assessments on coastal areas.In addition,we undertake comparative analysis on the characteristics of the Earth’s oscillation excited by the 2004 Sumatra-Andaman,2008 Wenchuan,and 2011Tohoku earthquakes that incorporate the effect of the Earth’s gravitational potential.A comparison of the amplitude spectra of the numerical records indicates that energy released by the three big earthquakes is different.Our comparative analysis realizes that the computing results can accurately reproduce some eigenfrequencies of the Earth,such as toroidal modes 0T2 to 0T13and spheroidal modes 0S7 to 0S31.These results demonstrate that numerical simulations can be successfully used to investigate the Earth’s oscillations.We propose that numerical simulations can be used as one of the major tools to further reveal how the Earth’s lateral heterogeneities affect the Earth’s oscillations.     

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.     

热源激发重力波特征以及波流作用的数值模拟研究

本文在二维等温可压大气中引入了一个随时间和空间变化的热源扰动,分别以静止风和中纬1月份月平均向东的纬向风急流为背景,对不同背景下热源激发的重力波的传播详细过程及其特性进行了数值模拟研究.热源激发出来的重力波在初始阶段有很宽的频谱范围,随后由于重力波的传播效应,水平波长和垂直波长分布范围随时间都有所减小.顺风传播的重力波的小尺度和低频部分会容易被急流吸收,从而加强了对流层急流,而逆风传播的重力波更容易上传,会导致中间层区域向西的背景风增强.这体现了低层大气急流对中间层大气风场结构的影响.热源的尺度直接决定激发波的尺度;激发波的垂直尺度和时间尺度对热源的变化比其水平尺度更敏感.     

震源参数对强地面震动模拟结果的影响

基于运动学震源模型,进行了不同震源参数情形下强地面震动数值模拟.结果表明,不同的破裂过程会产生差别甚大的强地面运动分布,一次确定性震源参数的模拟结果不能作为活动断层地震危害性评价的指标,只有通过大量三维地震动场模拟计算,给出地面震动评估的统计结果,才是比较合理的发展方向.由于一次三维地震动场计算耗时很大,因此解决问题的关键是如何考虑合理的震源参数.     

Spatial sensitivity of seismic hazard results to different background seismic activity and temporal earthquake occurrence models

Spatial sensitivity of seismic hazard results to different models with respect to background seismic activity and earthquake occurrence in time is investigated. For the contribution of background seismic activity to seismic hazard, background area source with uniform seismicity and spatially smoothed seismicity models are taken into consideration. For the contribution of faults, through characteristic earthquakes, both the memoryless Poisson and the time dependent renewal models are utilized. A case study, involving the assessment of seismic hazard for the Bursa province in Turkey, is conducted in order to examine quantitatively the influence of these models on seismic hazard results. The spatial variation of the difference in Peak Ground Acceleration (PGA) values obtained from these different models is presented in the form of difference maps for return periods of 475 and 2475 years. Best estimate seismic hazard maps for PGA and Spectral Accelerations (SA) at 0.2 and 1.0 s are obtained by using the logic tree method.     

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

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

A phased approach to enable hybrid simulation of complex structures

Hybrid simulation has been shown to be a cost-effective approach for assessing the seismic performance of structures. In hybrid simulation,critical parts of a structure are physically tested,while the remaining portions of the system are concurrently simulated computationally,typically using a finite element model. This combination is realized through a numerical time-integration scheme,which allows for investigation of full system-level responses of a structure in a cost-effective manner. However,conducting hybrid simulation of complex structures within large-scale testing facilities presents significant challenges. For example,the chosen modeling scheme may create numerical inaccuracies or even result in unstable simulations; the displacement and force capacity of the experimental system can be exceeded; and a hybrid test may be terminated due to poor communication between modules(e.g.,loading controllers,data acquisition systems,simulation coordinator). These problems can cause the simulation to stop suddenly,and in some cases can even result in damage to the experimental specimens; the end result can be failure of the entire experiment. This study proposes a phased approach to hybrid simulation that can validate all of the hybrid simulation components and ensure the integrity largescale hybrid simulation. In this approach,a series of hybrid simulations employing numerical components and small-scale experimental components are examined to establish this preparedness for the large-scale experiment. This validation program is incorporated into an existing,mature hybrid simulation framework,which is currently utilized in the Multi-Axial Full-Scale Sub-Structuring Testing and Simulation(MUST-SIM) facility of the George E. Brown Network for Earthquake Engineering Simulation(NEES) equipment site at the University of Illinois at Urbana-Champaign. A hybrid simulation of a four-span curved bridge is presented as an example,in which three piers are experimentally controlled in a total of 18 degrees of freedom(DOFs). This simulation illustrates the effectiveness of the phased approach presented in this paper.     

青海湖湖陆风的数值研究

用一个含植被参数化的二维中尺度模式研究青海湖夏季的湖陆风,模拟结果再次表明该模式具有重现局地中尺度环境特征,如冷（热）岛环境和湿岛等。同时表明,模式对植被覆盖度是敏感的,即有稀疏植被覆盖（σ＝３０％）情况的湖陆风白天比裸地情况要小,夜间则大,湖岸水平温度（或湿度）梯度随植被覆盖度增加而减小。     

轨道不平顺随机过程的数值模拟

轨道不平顺是列车-轨道系统的动力响应的主要影响因素,不平顺样本是列车-轨道系统动力分析模型中不可缺少的参数.根据随机理论原理推导了轨道不平顺谱的空间、时间域的相互转化关系,得到了时域功率谱的表达式.分析了目前广泛应用几种轨道不平顺随机过程的模拟方法的原理及作法,通过对功率谱快速数值算法的分析给出了一种新三角级数模拟方法,该方法与逆傅氏变换法具有等效性,同时也论证了逆傅氏变换方法样本是高斯过程.最后通过实际算例的分析论证了新三角级数法与傅氏逆变换法在简便、合理等方面的优点.     

2016年台湾美浓MW6.4地震震源参数的InSAR和GPS反演

2016年2月6日台湾西南部高雄市美浓区发生了M_W6.4地震.本文结合ALOS2卫星升降轨、Sentinel-1A升轨SAR数据,采用两轨差分干涉技术获取了该区域的同震形变场,形变结果表明震中西北部以抬升为主,最大视线向形变量约为11.2 cm.基于均匀位错模型和多峰值粒子群（MPSO）算法,利用InSAR和GPS形变数据联合反演了美浓地震的断层几何参数,结果表明震源中心位于22.920°N,120.420°E,深度约12 km,发震断层长度约15 km,走向角307°,倾角16.5°,平均滑动角为51.5°,此次地震是以逆冲倾滑兼左旋走滑的破裂模式.利用格网迭代搜索法得到最优倾角为15.7°,GPS和InSAR最优权比为18：1,最优平滑因子为0.06.基于非均匀位错模型,利用非负最小二乘方法进行线性反演,结果显示最大倾滑和走滑量分别为51.7 cm和55.3 cm,对应矩震级为M_W6.38,略小于GCMT（M_W6.4）的结果.通过与已有文献的比较和对该区域断层构造的分析,发现美浓地震的发震断层为单一断层的解释更为合理,我们推测发震断层是位于左镇、后甲里等断层之间的一条东南-西北走向往东北倾斜的盲断层,并初步推测2010年M_W6.3甲仙地震也同该断层有关.

     

地铁车站结构大型振动台试验与数值模拟的比较研究

根据可液化土层上土-地铁车站结构动力相互作用大型振动台模型试验结果,以软件ABAQUS为平台,将地基土-地铁车站结构体系视为平面应变问题,采用记忆型嵌套面黏塑性动力本构模型模拟土体的动力特性,采用混凝土动塑性损伤模型模拟车站结构混凝土的动力特性,建立了土-地铁车站结构非线性动力相互作用的有限元分析模型,对各种试验工况下地基土-地铁车站结构体系的地震反应进行了数值模拟,并与试验结果进行了对比,结果表明:数值模拟与振动台模型试验结果基本一致,体现出了相似的规律性,相互印证了计算分析的力学建模和振动台试验结果的正确性。