有限元方法在地震孕育过程与前兆机理研究中的应用进展

摘要/Abstract

摘要： 有限元方法的基本思想是将连续的求解区域离散为一组有限个、且按一定方式相互联结在一起的单元的组合体。有限元方法是目前应用最为广泛的数值模拟方法, 由于能够将复杂介质的力学本构关系及边界条件问题转化成为常规问题的计算程序, 所以越来越多的地震学者将其应用于地震成因与前兆机理研究中。该文着重介绍了有限元数值分析方法及其模拟在地震孕育过程与前兆机理研究中的应用进展, 并提出开展小尺度规模如岩石中矿物颗粒之间的力学行为的有限元模拟研究等几点建议。

关键词: 有限元方法, 孕育过程, 前兆机理, 数值模拟方法

Abstract: The finite element method, the basic idea of which is the decomposition of a continuous solution field into simpler components in a given way, is a numerical modeling method used widely. It is used in studying earthquake formation and mechanism of earthquake precursors by more and more researchers because of its advantage that it can treat the problems of mechanical constitutive properties of complicated media and boundary conditions in a simple way. This paper introduced the development of the finite element numerical modeling in study of the seismogenic process and mechanism of earthquake precursors. Several suggestions were presented for better use of the method in earthquake researches, for example, numerical modeling mechanical behaviors among mineral grains in rocks.

Key words: Finite element method, Seismic process, Mechanism of earthquake precursors, Numerical modeling method

中图分类号:

P315.7

陈玉香, 杜建国, 刘红. 有限元方法在地震孕育过程与前兆机理研究中的应用进展[J]. 地震, 2007, 27(4): 99-109.

CHEN Yu-xiang, DU Jian-guo, LIU Hong. Application of finite element method to study the seismic process and mechanism of earthquake precursors[J]. EARTHQUAKE, 2007, 27(4): 99-109.

参考文献

[1] 王勖成. 有限单元法[M]. 北京: 清华大学出版社, 2005.
[2] 张国民, 傅征祥, 桂燮泰, 等. 地震预报引论[M]. 北京: 科学出版社, 2001.
[3] 傅承义, 陈运泰, 祁贵仲, 等. 地球物理学基础[M]. 北京: 科学出版社, 1985.
[4] Richardson R M. Finite element modeling of stress in the Nazca plate: driving forces and plate boundary earthquakes[J]. Tectonophysics, 1978, 50: 223-248.
[5] 张东宁, 曾融生. 冀中坳陷滑脱构造动力的数值模拟[J]. 地震学报, 1995, 17(4): 414-421.
[6] Matsu'ura M, Sato T. Loading mechanism and scaling relation of large interplate earthquakes[J]. Tectonophysics, 1997, 277: 189-198.
[7] Hashimoto C, Matsu'ura M. 3-D simulation of earthquake generation cycles and evolution of fault constitutive properties[J]. Pure Appl Geophys, 2002, 159: 2175-2199.
[8] Hashimoto C, Matsu'ura M. 3-D physical modeling of stress accumulation and release processes at transcurrent plate boundaries[J]. Pure Appl Geophys, 2000, 157: 2125-2147.
[9] Aochi H, Matsu'ura M. Slip and time-dependent fault constitutive law and its significance in earthquake generation cycles[J]. Pure Appl Geophys, 2002, 159: 2029-2046.
[10] 李大鹏, 林邦慧, 谢云. 地震剪切破裂的三维断裂力学有限元分析—1966年邢台地震的数值模拟[J]. 中国地震, 1997, 13(2): 151-163.
[11] 郑文俊, 袁道阳, 张冬丽, 等. 1927年古浪8级地震的破裂习性及破裂机制的数值模拟[J]. 中国地震, 2004, 20(4): 353-363.
[12] 梅世蓉, 梁北援. 唐山地震孕震过程的数值模拟[J]. 中国地震, 1989, 5(3): 9-17.
[13] 冯德益, 刘喜兰, 舒立德, 等. 唐山7.8级地震孕育背景及长期前兆异常特征的三维数值模拟[J]. 地球物理学进展, 1996, 11(2): 105-113.
[14] 李步云. 昌马地震应力场及发震机制研究[J]. 西北地震学报, 1994, 16(3): 14-23.
[15] 宋惠珍, 兰印刚, 徐鹏飞, 等. 1303 年洪洞8级地震震源过程研究[J]. 地震地质, 1995, 17(1): 13-18.
[16] 周新伟. 北天山活动构造区的深浅构造特征与孕震模拟[D]. 北京: 中国地震局地质研究所, 2000.
[17] 李启成. 地震场源应变前兆的区别及其在地震预报中的应用[J]. 西北地震学报, 1995, 17(1): 1-9.
[18] Piatanesi A, Tinti S, Bortoluui E. Finite-element simulations of the 28 December 1908 Messina Straits(Southern Italy) tsunami[J]. Phys Chem Earth, 1998, 23(9): 933-938.
[19] 王妙月, 底青云, 张美根, 等. 地震孕育、发生、发展动态过程的三维有限元数值模拟[J]. 地球物理学报, 1999, 42(2): 218-227.
[20] 王继存, 续春荣. 唐山地震断层破J坏及其力学过程的数值模拟[J]. 地震地质, 1989, 11(4): 71-76.
[21] 陶振宇, 张黎明. 与唐山地震有关的深部岩石力学特性的数值模拟[J]. 西北地震学报, 1991, 13(2): 22-28.
[22] 孟国杰, 于莉, 高锡铭. 唐山地震的非均匀错动模型-水准资料的反演结果[J]. 地壳形变与地震, 1995, 15(2): 7-14.
[23] 王成宝, 张之立. 1966年邢台地震序列破裂过程的有限元计算[J]. 华北地震科学, 1990, 8(4): 31-40.
[24] 陈祖安, 白武明, 林邦慧, 等. 1966年以来华北地区一系列七级大震破裂过程的数值模拟[J]. 地球物理学报, 2003, 46(3): 373-381.
[25] 任春, 罗奇峰. 断层破碎带对非发震断层场地地震动的影响[J]. 地震研究, 2005, 28(2): 162-166.
[26] 白超英, 秦保燕. 深部剪切形变带对浅源地震的控制-立交模式有限元的模拟计算[J]. 西北地震学报, 1990, 12(1): 1-11.
[27] 李步云. 甘宁交界地区弧形断裂失稳理论的模拟研究[J]. 高原地震, 1999, 11(2): 22-30.
[28] Oglesby D D, Archuleta R J, Nielsen S B. The three-dimensional dynamics of dipping faults[J]. Bull Seism Soc Am, 2000, 90: 616-628.
[29] Oglesby D D, Archuleta R J. The three-dimensional dynamics of a nonplanar thrust fault[J]. Bull Seism Soc Am, 2003, 93(5): 2222-2235.
[30] Oglesby D D, Archuleta R J., Nielsen S B. Earthquakes on dipping faults: The effect of broken symmetry[J]. Science, 1998, 280: 1055-1059.
[31] Oglesby D D, Archuleta R J, Nielsen S B. The dynamics of dip-slip faulting: explorations in two dimensions[J]. J Geophys Res, 2000, 105(B6): 13643.
[32] 刘启方, 丁海平, 袁一凡, 等. 三维地震断层动力破裂的显式并行有限元解法[J]. 地震工程与工程振动, 2003, 23(4): 22-28.
[33] 孙荀英, 刘激扬, 王仁. 1976年唐山地震震时和震后变形的模拟[J]. 地球物理学报, 1994, 37(1): 45-55.
[34] 黄立人, 马青, 郭良迁, 等. 华北部分地区水平变形的力学机制-三维有限单元计算和GPS复测结果的分析[J]. 地震学报, 1999, 21(1): 50-56.
[35] 白武明, 柳邦慧, 陈祖安. 1976年唐山大震发生对华北地区各地块运动与变形影响的数值模拟研究[J]. 中国科学(D辑), 2003, 33(增刊): 99-107.
[36] 陈连旺, 陆远忠, 张杰, 等. 1966年邢台地震引起的华北地区应力场动态演化过程的三维粘弹性模拟[J]. 地震学报, 2001, 23(5): 480-491.
[37] 陈连旺, 陆远忠, 郭若眉, 等. 华北地区断层运动与三维构造应力场的演化[J]. 地震学报, 2004, 23(4): 349-361.
[38] 冯向东, 魏东平, 陈棋福. 基于观测应力场的大华北地区动力学机制探讨[J]. 地震学报, 2005, 27(1): 1-10.
[39] 李志才, 许才军, 赵少荣, 等. 基于地壳分层的唐山地震断层震后变形分析[J]. 地球物理学进展, 2005, 20(4): 961-968.
[40] Hu Y, Wang K, He J, et al. Three-dimensional viscoelastic finite element model for postseismic deformation of the great 1906 Chile earthquake[J]. J Geophys Res, 2004, 109, B12403.
[41] Alik Ismail-Zadeh, Birgit Mueller, Gerald Schubert. Three-dimensional numerical modeling of contemporary mantle flow and tectonic stress beneath the earthquake-prone Southeastern Carpathians based on integrated analysis of seismic, heat flow, and gravity Data[J]. Physics of the Earth and Planetary Interiors, 2005, 149: 81-98.
[42] Ismail-Zadeh A T, Keilis-Borok V I, Soloviev A A. Numerical modeling of earthquake flows in the southeastern Carpathians(Vrancea): effect of a sinking slab[J]. Phys Earth Planet Int, 1999, 111: 267-274.
[43] Ismail-Zadeh A T, Korotkii A I, Naimark B M. Numerical modelling of three-dimensional viscous flow under gravity and thermal effects[J]. Comput Math Phys, 2001, 41(9): 1331-1345.
[44] 王仁, 何国琦, 殷有泉, 等. 华北地区地震迁移规律的数学模拟[J]. 地震学报, 1980, 1(1): 32-42.
[45] 王仁, 孙荀英, 蔡永恩. 华北地区近年地震序列的数学模拟[J]. 中国科学(B辑), 1982, 18(8): 715-753.
[46] 焦明若, 张国民, 刘杰. 中国大陆及其邻区地震活动的数值模型研究[J]. 地震学报, 1999, 21(6): 583-590.
[47] 陈化然, 蒋淳, 王俊国, 等. 强震活动有限元模型研究及其初步应用-以西南地区为例[J]. 中国地震, 2001, 18(3): 249-258.
[48] 陈化然, 李轶群, 汪翠芝, 等. 基于断层相互作用的地震活动有限元模型[J]. 大地测量与地球动力学, 2002, 22(2): 86-90.
[49] 王凯英, 马瑾. 川滇地区断层相互作用的地震活动证据及有限元模拟[J]. 地震地质, 2004, 26(2): 259-272.
[50] 陈忠辉, 傅宇方, 唐春安. 组合震源模式的数值模拟研究[J]. 西北地震学报, 1997, 19(1): 15-21.
[51] 刘杰, 石耀霖, 张国民. 基于准三维有限元方法建立的地震活动模型[J]. 地球物理学报, 2001, 44(6): 814-824.
[52] 李平, 卢良玉, 卢造勋, 等. 辽宁及邻区地壳构造应力场及其与地震活动关系的三维有限元数值模拟研究[J]. 地震学报, 2001, 23(1): 24-35.
[53] 陆明勇, 郑文衡. 地壳块体动应力变化特征及其与地震关系的初步研究[J]. 地震学报, 2005, 27(3): 282-291.
[54] 郑勇, 傅容珊, 熊熊. 中国大陆及周边地区现代岩石圈演化动力学模拟[J]. 地球物理学报, 2006, 49(2): 415-427.
[55] 丁原章, 黄新辉. 华南滨海断裂带地震危险性的数学模拟[J]. 华南地震, 1995, 15(1): 12-20.
[56] 闻则刚, 王正尚, 杨马陵. 广东地区现代构造应力场的模拟计算[J]. 华南地震, 2001, 21(3): 29-36.
[57] 刘洁, 宋惠珍. 用数值模拟方法评估华北北部地震危险性[J]. 地震地质, 1999, 21(3): 221-227.
[58] 钟羽云, 姚立珣. 浙北及邻近地区应力场数值模拟与地震危险区预测[J]. 地震学刊, 2000, 20(3): 9-17.
[59] 李步云, 杨斌, 刘洪春, 等. 三维有限元数值模拟在宁夏大坝电厂地震安全性评价中的应用[J]. 西北地震学报, 1995, 17(3): 80-87.
[60] 文彦君, 郑文涛. 延怀盆地现代构造应力场模拟与官厅水库大坝地震安全性初探[J]. 地震地质, 2004, 26(1): 141-152.
[61] 王爱国, 杨斌, 马巍, 等. 数值模拟在黑山峡两坝址构造稳定性评价中的应用[J]. 地震研究, 2004, 27(4): 350-356.
[62] Bird P. Testing hypotheses on plate driving mechanisms with global lithosphere models including topography, thermal structure, and faults[J]. J Geophys Res, 1998, 103: 10115-10129.
[63] Bird P, Zhen Liu. Global finite-element model makes a small contribution to intraplate seismic hazard estimation[J]. Bull Seis Soc Am, 1999, 6: 1642-1647.
[64] 张永仙, 石耀霖, 张国民. 孕震过程中波速比及S波分裂时空变化的数值模拟[J]. 地震,1994, 14(6): 26-33.
[65] 张永仙, 石耀霖, 张国民. 孕震过程中孔隙压及地下水位变化的数模拟[J]. 地震, 1994, 14(1): 65-72.
[66] 石耀霖, 张永仙, 张国民. 孕震过程中岩石膨胀及断层蠕动的数值模拟研究[J]. 地震学报, 1994, 16(增刊): 10-17.
[67] 梁子斌, 钱家栋. 用二相介质有限元方法对共和7.0级地震前兆的数值模拟[J]. 西北地震学报, 1997, 19(3): 69-75.
[68] 董治平, 郭守年, 蒋梅, 等. 天水地震区深部地温场的有限元计算[J]. 西北地震学报, 1996, 18(3): 48-53.
[69] 张永仙. 热应力三维数值模拟研究及热物质的运移在地震孕育过程中的作用探讨[D]. 北京: 中国地震局地球物理研究所, 1999.
[70] Mark D Behn, Margaret S Boettcher, Greg Hirth. Thermal structure of oceanic transform faults[J]. Geology, 2007, 35(4): 307-310.
[71] 郑文衡, 王乘, 陈湘鹏. 地壳块体脉冲响应的时间反转效应研究[J]. 地震学报, 2001, 23(4): 370-379.
[72] 张杰卿, 何巧云, 陈化然. 地震应力触发数值模拟的初步研究[J]. 东北地震研究, 2004, 20(3): 47-53.
[73] 胡志奇, 郑文衡, 陆明勇. 岩体动态应力在地震机理中的作用研究[J]. 大地测量与地球动力学, 2005, 25(1): 108-112.
[74] 杜建国, 康春丽. 强地震前兆异常特征与深部流体作用探讨[J]. 地震, 2000, 20(3): 95-101.
[75] 徐常芳. 中国大陆壳内与上地幔高导层成因及唐山地震机理研究[J]. 地学前缘, 2003, 10(增刊): 101-111.
[76] 车用太, 金鱼子. 地震地下流体学[M]. 北京: 气象出版社, 2006.
[77] Gabrielov A M, Keilis-Borok V I, Pinsky V, et al. Fluids migration and dynamics of a blocks-and-faults system[J]. Tectonophysics, 2007, 429: 229-251.
[78] Loo Huan-yen, Gao W A, Yuan Y G. Three-dimensional modeling of the driving mechanism of the 1976 Tangshan earthquake in China[J]. J of Geodynamics, 1989, 11: 243-265.
[79] 刘武英, 李龙海. 唐山地震与海城地震之共性特征及有关问题的讨论[J]. 地震研究, 1995, 18(2): 161-167.
[80] 叶正仁, 王建. 中国大陆现今地壳运动的动力学机制[J]. 地球物理学报, 2004, 47(3): 456-461.
[81] 张国民, 李丽. 地壳介质的流变性与孕震模型[J]. 地震地质, 2003, 25(1): 1-10.
[82] Kari Moisio. Numerical lithospheric modelling: rheology, stress and deformation in the central Fennoscandian Shield[D]. Finland: University of Oulu, 2005.
[83] 刘武英, 李龙海, 吴建华, 等. 从地球化学角度讨论1976年唐山7.8级地震的成因[J]. 地质地球化学, 1996, 5: 66-69.
[84] 刘巍, 杜建国, 丁天才. 地球内部矿物相变及其地质意义[J]. 现代地质, 2001, 15(2): 211-215.
[85] 何昌荣, 周永胜, 桑祖南. 攀枝花辉长岩半脆性-塑性流变的实验研究[J]. 中国科学(D辑), 2002, 32(9): 717-726.
[86] Majbritt D Lund, Sandra Piazolo, Simon L Harley. Ultrahigh temperature deformation microstructures in felsic granulites of the Napier Complex, Antarctica[J]. Tectonophysics, 2006, 427: 133-151.