数值模拟在地球动力学中的研究进展
详细信息 本馆镜像全文    |  推荐本文 | | 获取馆网全文
摘要
本文主要针对有限单元法,对近10几年来数值模拟在地球动力学中的应用和发展作了回顾,分为构造应力场模拟;地幔热对流模拟;板块碰撞模拟;岩石圈流变学模拟;以及地震机制与预测模拟等部分.简单阐述了当前地球动力学数值模拟的发展趋势并论述了其面临的主要问题.
Just aimed at finite element method,this paper gives a general review on the development and application of the numerical modeling in the Geodynamics research since the nineties.It contains sections on: numerical simulation of stress field and deformation field in structural geology;simulation on thermal convection of the mantle;simulation on the plate collision and convergence,including Continent-Continent collision and Ocean-Continental collision;modeling of Lithospheric Rheology;simulation on the mechanical of earthquake,and so on.At the end,this paper discusses the development tendency and mainly facing questions of numerical modeling.
引文
[1]Love A E H.Some problems of Geodynamics,Cambridge Uni-versity press.1911.
    [2]王仁,地球动力学的历史和近期发展[J].地球物理学进展,1996,11(1),1~10.
    [3]Clough R W.The finite element method in plane stress analy-sis[A].Proceedings of 2nd ASCE Conference on ElectronicComputation[C].Pittsburgh:ASCE,1960.
    [4]王仁,有限单元等数值方法在我国地球科学中的应用和发展[J].地球物理学报,1994,37增刊1,128~138.
    [5]石耀霖,地幔对流研究的一些新进展[J].地球科学进展,2001,16(4):496~500.
    [6]Pysklywec R N,Beaumont C,Fullsack P,Modeling the be-havior of the continental mantle lithosphere during plate con-vergence[J].Geology,2000,28(7):655~658.
    [7]Beaumont C,Jamiesont R A,Nguyen M H,et al.Himalayantectonics explained by extrusion of a low~viscosity crustalchannel coupled to focused surface denudation[J].Nature,2001,414:738~742.
    [8]Houseman G,England P C.A lithospheric-thickening modelfor the Indo-Asian collision in tectonic evolution of Asia[A],in A.Yin&T.M.Harrison(eds)[C],Cambridge universitypress,New York,1996,pp3~17.
    [9]Vanderhaeghe O,Medvedev S,Fullsack P,et al,Evolution oforogenic wedges and continental plateaus:insights from ther-mal~mechamical models with subduction basal boundary con-ditions[J].Geophys.J.Int,2003,153(1):27~51.
    [10]Kasapoglu K E,Toksoz M N,Tectonic consequences of thecollision of the Arabian and Eurasian plates:finite elementmodels[J].Tectonophysics,1983,100:71~95.
    [11]Farhad S,Jafar A H,Thermo-mechanical modeling of sub-duction of continental lithosphere[J].Physics of the Earthand Planetary Interiors,2002,131:185~203.
    [12]Magali I.Billen,Michael G,A low viscosity wedge in sub-duction zones[J].Earth and Planetary Science Letters,2001,193,227~236.
    [13]Shi Y L,Wang C Y.Two-demensional modeling of the P-T-tpaths of regional metamorphism in simple overthrust terrains[J].Geology,1987,15:1048~1051.
    [14]Liu M,Yang Y Q,Extensional collapse of the Tibetan Plat-eau:Results of three-dimensional finite element modeling[J].J.Geophys.Res.,2003,108(8):1~15.
    [15]Pysklywec R N.Evolution of subducting mantle lithosphereat a continental plate boundary[J].Geophysical Reseach Let-ters,2001,28(23):4399~4402.
    [16]Pysklywec R N,Beaumont C,Fullsack P.Lithospheric de-formation during the early stages of continental collision:Nu-merical experiments and comparison with South Island,NewZealand[J].J.Geophys.Res.,2002,107(B7):ETG 1~19..
    [17]England P,Molnar P,Active deformation of Asia:from kine-matics to dynamics[J].Science,1997,278:647~650.
    [18]Zhong S J,Zuber M T.Crustal compensation during moun-tain-building[J].Geophysical Reseach Letters,2000,27(18):3009~3012.
    [19]Soofi M A,King S D.Oblique convergence between India andEurasia[J].J.Geophys.Res.,2002,107(b5):ETG 1~8.
    [20]Pysklywec R N,Shahnas M H.Time-dependent surface to-pography in a coupled crust-mantle convection model[J].Geo-phys.J.Int.,2003,154:268~278.
    [21]Zhong S J,Michael G,Dynamic interaction between tectonicplates,subducting slabs,and the mantle[J].Earth Interac-tions,1997,1(3):1~18.
    [22]C.P.Conrad,B.H.Bager,Mantle convection with strongsubduction zones[J].Geophys.J.Int.,2001,144:271~288.
    [23]Robl J,Stüwe K.The role of indenter rheology during conti-nental collision:Concepts and application to the EasternAlps.(GeoMod Abstr.vol.);Bollettino di Geofisica.2004,45(suppl.1):133~135.
    [24]Hunen J v,A.P.van den Berg,N.J.Vlaar,Variousmechanisms to induce present-day shallow flat subduction andimplications for the younger Earth:a numerical parameterstudy[J].Physics of the Earth and Planetary Interiors,2004,146:179~194.
    [25]Michael A.Eberle,Olivier Grasset,Christophe Sotin,Anumerical study of the interaction between the mantle wedge,subducting slab,and overriding plate[J].Physics of theEarth and Planetary Interiors,2002,4158:1~12.
    [26]Gerbault M,Willingshofer E.Lower crust indentation orhorizontal ductile flow during continental collision?[J].Tec-tonophysics,2004,387:169~187.
    [27]Willett S,Beaumont C,Fullsack P.Mechanical model for thetectonics of doubly vergent compressional orogens[J].GEOL-OGY,1993,21:371~374.
    [28]Sean D.Willett,Christopher Beaumont,Subduction of Asianlithospheric mantle beneath Tibet inferred from models ofcontinental collision[J].Nature,1994,369:642~645.
    [29]Liu Z,Peter B.Finite element modeling of neotectonics inNew Zealand[J].J.Geophys.Res,2002,107(B12):ETG 1~18.
    [30]汪素云,陈培善.中国及邻区现代构造应力场的数值模拟[J].地球物理学报,1980,23(1):35~45.
    [31]焦明若,张国民,车时,刘杰.中国大陆及其周边地区构造应力场的数值计算及其在地震活动性解释上的应用[J].地震学报,1999,21(2):123~132.
    [32]梅世蓉,车时.我国现代构造应力场方向与大小的数值模拟及其与地震活动的关系[A].见:梅世蓉地震科学研究论文选集(1993)[C],北京:地震出版社42:8~447.
    [33]肖兰喜,朱元清,陶九庆,杜宪宋.岩石圈流变强度与中国大陆构造运动关系的探讨[J].西北地震学报,2003,25(4):304~311.
    [34]安美建,石耀霖,李方全.用遗传有限单元反演法研究东亚部分地区现今构造应力场的力源和影响因素[J].地震学报,1998,20(3):225~231.
    [35]尹京苑,梅世蓉.地壳应力状态及其在三维模拟计算中的意义[J].地震,1998,18(3):226~232.
    [36]傅容珊,黄建华,常筱华,曹抗震.东亚大陆变形应力场格局演化的数值模拟[J].地壳变形与地震,2002,20(3):1~10.
    [37]陈开平,马谨.印度与欧亚大陆碰撞构造变形数值分析[J].地震地质,1995,17(3):277~284).
    [38]张东宁,许忠淮.青藏高原现代构造应力状态及构造运动的三维粘弹性数值模拟[J].中国地震,1994,10(2):136~143.
    [39]张东宁,许忠淮.中国大陆岩石层动力学数值模型的边界条件[J].地震学报,1999,21(2):133~139.
    [40]张健,石耀霖.中国西部地区重力位能与板内变形动力[J].中国科学院研究生院学报,2001,18(1):43~50.
    [41]张健,石耀霖.青藏高原隆升及伸展变形中的重力位能[J].地球物理学报,2002,45(2):226~232.
    [42]熊熊,刘孙君,许厚泽.青藏高原地壳东-西向拉张及力学机制[J].武汉大学学报,2003,28特刊:150~154.
    [43]张东宁,许忠淮.青藏高原南部上地壳正断层地震活动的一种可能解释[J].地震学报,1995,17(2):188~195.
    [44]张东宁,许忠淮.西藏南部地堑构造成因的数值模拟[J].中国地震,1997,13(4):349-357.
    [45]宋惠珍,曾海容,DenisHeliot,剡公瑞.逆冲断层应力场的数值模拟[J].地震地质,1999,21(3):275~282.
    [46]武红岭,陈柏林.郯庐断裂带构造演化的数值分析[J].地质力学学报,1996,2(2):66~74.
    [47]Cai Y,He T,Wang R,Numerical simulation of Dynamicprocess of the Tangshan earthquake by a new method-LDDA[J].Pure Apple Geophys,2000,157:2083~2104.
    [48]熊熊,许厚泽.地幔对流诱发的莫霍面起伏与下地壳塑性流动[J].测绘学报,2000,29增刊:21~25.
    [49]熊熊,许厚泽,傅容珊.大陆岩石圈的增厚及对流剥离对青藏高原隆升的影响[J].地壳形变与地震,1998,18(3):1~3.
    [50]傅容珊,黄建华,李力刚,常筱华.青藏高原隆升三阶段模型的数值模拟[J].地学前缘,2000,7(4):588~596.
    [51]熊熊,傅容珊,许厚泽,李平.增厚大陆岩石圈热边界层对流剥离的数值模拟[J].地球物理学报,1998,41增刊:33~40.
    [52]石耀霖,朱元清,沈显杰.青藏高原构造热演化的主要控制因素[J].地球物理学报,1992,35(6):710~72.
    [53]董文杰,汤懋苍.青藏高原隆升和夷平过程的数值模型研究[J].中国科学(D辑),1997,27(1):65~69.
    [54]沈显杰,Ladislaus Rybach.喜马拉雅-青藏地区的运动学-抬升模型[J].地震学报,1994,16(3):326~335.
    [55]傅容珊,徐耀民,黄建华,李力刚.青藏高原挤压隆升过程的数值模拟[J].地球物理学报,2000,43(3):346~355.
    [56]熊熊,滕吉文.青藏高原东缘地壳运动与深部过程的研究,地球物理学报[J].2002,45(4):507~515.
    [57]熊熊,许厚泽,滕吉文.青藏高原物质东流力学背景探讨[J].地壳形变与地震,2001,21(2):1~7.
    [58]金文,石耀霖.高喜马拉雅P-T-t轨迹热模拟[J].中国科学院研究生院学报,1994,11(2):187~195.
    [59]石耀霖,王其允.高喜马拉雅淡色花岗岩形成的热模拟[J].地球物理学报,1997,40(5):667~676.
    [60]孙荀英,王仁,王其允.海沟倒退对地幔对流的影响[J].地球物理学报,1994,37(6):738~748.
    [61]张健,石耀霖.东亚陆缘带构造扩张的深部热力学机制[J].大地构造与成矿学,2003,27(3):222~227.
    [62]臧绍先,宁杰远,陈玉文,周丽梅.两种地幔对流模式下俯冲带的热结构[J].地球物理学报,1994,37(4):448~455.
    [63]张健,石耀霖.海洋和大陆岩石圈下的混合地幔对流模型[J].中国科学院研究生院学报,1998,15(2):157~164.
    [64]宁杰远,臧绍先.对俯冲带深震成因的探讨[J].地震学报,1999,21(5):523~532.
    [65]宁杰远,臧绍先.对俯冲带波速结构的数值模拟[J].地球物理学报,2001,44(2):190~198.
    [66]Scott D.King,Subduction zones:observations and geody-namic models[J].Physics of the Earth and Planetary Interi-ors,2001,127:9~24.
    [67]金维浚,石耀霖.大别山超高压变质带P~T~t轨迹数值模拟[J].地学前缘,1998,5增刊:127~134.
    [68]范桃园,石耀霖.大别-苏鲁超高压变质带P~T~t轨迹的动力学模拟[J].地球物理学报,2001,44(5):627~635.
    [69]武红岭,董树文.大别山陆—陆点碰撞和构造超压的形成[J].地球科学,2001,26(5):457~463.
    [70]Tapponier P,Peltzer G,Le Dain A Y,Armijo R,P,C.Propagating extrusion tectonics in Asia:new insights fromsimple experiments with plasticine[J].Geology,1982,10:611~616..
    [71]England PC,Houseman GA.Finite strain calculations ofcontinental deformation;1,Method and general results forconvergent zones[J].J.Geophys.Res.,1986,91(3):3664~3676.
    [72]Houseman G A,McKenzie D P,Molnar P.Convective insta-bility of a thickened boundary layer and its relevance for thethermal evolution of continental onvergent belts[J].J.Geo-phys.Res.,1981,86,6115~6132,.
    [73]Beaumont C,Kamp P J J,Hamilton J,Fullsack P.The con-tinental collision zone,South Island,New Zealand:compari-son of geodynamical models and observations[J].J.Geo-phys.Res.,1996,101,3333~3359,.
    [74]England P,Houseman G.Extension during continental con-vergence,with application to the Tibetian Plateau[J].J.Geophys.Res.,1989,94,17561~17579..
    [75]张贵宾,高锐,肖序常.横过西昆仑和塔里木结合带的山隆盆降机制的动力学模拟[J].地球学报,2001,22(6):541~546.
    [76]石耀霖,王其允.斜俯冲板块边界变形分配的力学分析[J].地球物理学报,1994,37(5):606~612.
    [77]Joel L,King S D.The influence of thermodynamic formula-tion on simulations of subduction zone geometry and history[J].Geophysical Reseach Letters,1998,25(9):1463~1466.
    [78]Royden L H,Clark B B,Robert W.K,et al.Surface deform-ation and lower crustal flow in eastern Tibet[J].Science,1997,276:788~790.
    [79]石耀霖,朱守彪.中国大陆震源机制深度变化反映的地壳-地幔流变特征[J].地球物理学报,2003,46(3):359~365.
    [80]王仁,何国琦,殷有泉,等.华北地区地震迁移的数学模拟[J].地震学报,1980,2(1):32~42.
    [81]曹雪峰,尹京苑,杨立明.南北地震带应力场演变的数值模拟研究,西北地震学报,2003,25(3):221~225.
    [82]Shi Y l,et al,Nonlinear Dynamic Model of Earthquake Pre-diction[J].Geophysics Series,1994,18:81~90.
    [83]张国民,耿鲁明,石耀霖.中国大陆强震轮回活动的计算机模型研究[J].中国地震,1993,9(1):20~32.
    [84]耿鲁明,张国民,石耀霖.地震孕育发生的场源关系初步研究[J].中国地震,1994,9(4):310~319.
    [85]王妙月,底青云,张美根,等.地震孕育、发生、发展动态过程的三维有限元数值模拟[J].地球物理学报,1999,42(2):218~227.
    [86]张永仙,石耀霖,刘桂萍.热物质上涌与震前重力异常关系的初探[J].地震,2000,20(增刊):135~143.
    [87]何秀琴.地壳变形与地应力变化关系的研究[J].高原地震,1994,6(2):30~38.

版权所有:© 2023 中国地质图书馆 中国地质调查局地学文献中心