地震地质 ›› 2006, Vol. 28 ›› Issue (4): 612-622.

• 研究论文 • 上一篇    下一篇

近场强地震动数值模拟的简化计算方法

张冬丽1,2, 陶夏新3, 周正华4   

  1. 1. 中国地震局地质研究所, 北京, 100029;
    2. 中国地震局地震预测研究所兰州科技创新基地, 兰州, 730000;
    3. 哈尔滨工业大学, 哈尔滨, 150083;
    4. 中国地震局工程力学研究所, 哈尔滨, 150083
  • 收稿日期:2006-05-24 修回日期:2006-08-20 出版日期:2006-12-21 发布日期:2009-08-27
  • 作者简介:张冬丽,女,1974年生,副研究员,2000年在中国地震局兰州地震研究所获理学硕士学位,构造地质学专业,2005年在中国地震局工程力学研究所获博士学位,岩土工程专业,现在中国地震局地质研究所做博士后研究工作,主要从事工程地震和近场强地震动数值模拟预测等研究,电话:010-62009416,E-mail:zdl516@163.com.
  • 基金资助:

    国家发展与改革委员会发改投资“城市活断层试验探测”项目(20041138);国家自然科学基金(40604005);地震科学联合基金(A07027);中国地震局“十五”重大建设项目(J105B16)共同资助

SIMPLIFIED METHOD FOR NUMERICAL SIMULATION OF NEAR-FIELD STRONG GROUND MOTION

ZHANG Dong-li1,2, TAO Xia-xin3, ZHOU Zheng-hua4   

  1. 1. Institute of Geology, China Earthquake Administration, Beijing 100029, China;
    2. Lanzhou Base of Institute of Earthquake Prediction, China Earthquake Administration, Lanzhou 730000, China;
    3. Harbin Institute of Technology, Harbin 150083, China;
    4. Institute of Engineering Mechanics, China Earthquake Administration, Harbin 150083, China
  • Received:2006-05-24 Revised:2006-08-20 Online:2006-12-21 Published:2009-08-27

摘要:

近场强地震动除受场地条件的影响外,还受到震源破裂面上子源的空间分布特点、子源破裂先后顺序的强烈控制,基于数值格林函数法的近场强地震动数值模拟方法可以综合考虑震源、传播途径及局部场地条件的影响,对计算过程进行合理简化,分2步完成地震动模拟:第1步,在介质均匀区采用矩张量的解析解计算所有子源在盖层底面的位移,形成下一步有限元计算的输入场;第2步,在盖层介质不均匀区,结合局部人工透射边界技术,采用时、空解耦的波动显式有限元方法计算地表强地震动.在有限断层模型中,采用具有9个力偶的等效地震矩张量表达断层产状、滑动方向等的影响,采用Brune模型定义各子源的滑动时间函数,描述滑动的时、空不均匀分布特征,从而细化震源模型.通过对Northridge地震中4个基岩台站地表地震动的模拟结果和强震记录,验证了此简化计算方法的可行性.

关键词: 近场强地震动, 数值模拟, 数值格林函数法, Northridge地震

Abstract:

Besides the site condition,the distribution of sub-sources and expansion of sub-source rupture on finite fault effect the near-field ground motions strongly,too.In order to simulate strong ground motion effectively,we develop a simplified numerical simulation method by numerical Green's function's computation.The whole process can be divided into two steps.At first,in infinite uniform elastic medium,take the finite fault as an exterior source under asymmetry layers,the displacement time histories of overburden layers' bottom caused by all sub-sources can be synthesized by analytic solution.The results can be used as the input for the finite element numerical simulation.Secondly,in asymmetry medium,with local artificial transmitting boundary condition,the explicit decoupled finite element simulation approach is applied to model the near-field strong ground motions.Comparing with former source studies,two ideas are added into the source model.One is that,instead of simple couples source,the earthquake moment tensor is introduced to describe fault attitude and the rupture mechanism,which could express the relationships among force,slip direction and slip-time function of every sub-sources in detail.The dislocation could by equivalent to three-order earthquake moment tensor with 9 force couples and every couple expression includes source parameters,such as fault strike,dip angle and rake angle.The other idea is about the slip-time function.On a finite fault,the inhomogeneous characteristics can be expressed by the nonuniform dislocation distributions and the nonuniform rupture expanding in temporal and spatial space.In this paper,the Brune model is used to characterize the slip-time function.For a certain sub-source,three parameters can define the function,average dislocation,rising time and the rupture time delay.It is feasible of the simplified method of strong ground motion simulation by an example test.According to the source parameters and velocity structures of the 1994 Northridge MW 6.7 earthquake,long period ground motions of four rock stations are simulated and compared with their records.The simulation acceleration histories and velocity histories are coincident well with the records in durations,amplitudes and envelops.The large acceleration pulse and velocity pulse in the simulation are clear.The simulation results also express engineering characters of hanging wall effect and rupture directivity effect distinctly.All above are well in accord with the ground motion records and earthquake investigation which prove that the approach is reasonable and applicable for seismic hazard evaluation and prediction.

Key words: near-field strong ground motion, numerical simulation, numerical Green’s function method, Northridge earthquake

中图分类号: