用户名: 密码: 验证码:
球面数值流形方法及其在地壳运动中的应用研究
详细信息    本馆镜像全文|  推荐本文 |  |   获取CNKI官网全文
摘要
以甚长基线干涉测量(VLBI)、卫星激光测距(SLR)、全球定位系统(GPS)和差分合成孔径雷达干涉(D-InSAR)等新兴的空间对地观测技术为代表的现代大地测量为监测和研究地壳运动提供了新途径。大地测量参与地壳形变的主要技术GPS可以10~(-9)的精度监测点距几十至几百公里不同尺度范围的地壳地表运动和形变,用VLBI和SLR技术可以测定相距数千里的全球板块运动;现代精密重力仪可以微伽的精度测定重力场的变化。通过这些运动信息的数据处理,可以直观地描述地壳运动规律。然而地表的地壳运动是地球深部动力学过程的一种表象,研究地壳运动,需要结合地质、地球物理等资料,探讨构造运动的动力学机制,研究构造运动的力源。
     地球物理、地质等学科的观测结果提供的是几十年至上百万年的空间上不连续的资料;虽然GPS观测网络越来越多,但相对于广袤的的大陆来说测站数仍严重不足,同时由于GPS观测主要以地震监测和预报为目的,测站分布多位于一些主断裂或构造活跃地区,地理分布极不均匀,许多地域受观测条件限制而成为观测盲区。因此,要获得特定块体或区域的地壳位移场和应力场的连续图像,必须要借助数理方法,由已知测站的观测得来的速度或位移来推求未被GPS观测覆盖地区的形变位移,并进而作应力分析。
     在已有数值模拟工具当中,有限元法因其适用性广、解算过程规范而得到广泛的应用,已成为解决数值正反演问题的主要工具。非连续变形分析方法(DDA)中所有单元是被不连续缝所包围的隔离块体,对每个块体允许有位移和变形,块体间允许有滑移。非连续变形分析方法侧重于处理不连续的问题,它要求区域内所有块体都是被完全分割开的。
     本文提出了利用数值流形方法进行地壳形变的数值模拟。数值流形方法是利用现代数学“流形”的有限覆盖技术建立起来的。有限覆盖是由物理覆盖和数学覆盖组成,通过它可以处理连续问题和非连续问题,解决了有限元只能计算连续和非连续变形分析方法侧重于不连续问题的不足,这两种方法都是数值流形方法的特例。
     本文最主要的工作是将数值流形方法从平面扩展到了球面,以适应模拟大范围地壳运动和形变的要求。
With the rapid progress of Very Long Baseline Interferometry(VLBI), Satellite Laser Ranger, Global Positioning System (GPS) and Differential Interferometry Synthetic Aperture Radar(D-InSAR), modern geodesy provides wholly new ways to monitor and research crustal movement. The relative accuracy 10~(-9) of baseline measurement can be achieved between GPS stations several hundred kilometers apart. Techniques such as VLBI and SLR are used to monitor global plate's movement. The results from GPS, VLBI and SLR can directly describe the crustal movement very well. However, if we want to explore the geodynamic mechanism and the driving forces of the plates, related geological and geophysical materials should be combined with geodetic data.Results from geophysical and geological data are in terms of hundred years and millions years respectively. They can not provide continuous current crustal movement field. The coverage of GPS stations is not enough compared with the big area of plate and block. For the observation networks are mainly set up to monitor earthquake and tectonic activity, GPS stations are mostly deployed along faults. So the distribution of GPS stations is not as even as desired. In order to map the velocity field and stress field of crustal movement, numerical tools are used to simulate the movement and deformation of the district without geodetic data and geophysical observation. Among the simulation tools Finite Element Method, which is based on the idea of dividing a complicated domain into small and manageable pieces, has been used in many fields. It can deal with continuous deformation very well. Discontinuous Deformation Analysis is capable of dealing with discontinuity, it requires that all blocks be totally separated by faults.In this dissertation Numerical Manifold Method is proposed to simulate crustal deformation. Here, the most important work is developing spherical Numerical Manifold Method from plane Numerical Manifold Method. A set of formulae is deduced in order that large scale crustal movement could be properly simulated with this method.
    Main research work is as following:1) The advantage and disadvantage of formerly used numerical simulation techniques are well discussed. FEM and DDA are good at dealing with continuum and discontinuity respectively. They are two special cases of NMM, which incorporate continua and discontinuity in a single model.2) Spherical numerical manifold method is firstly developed from NMM in order that it can deal with large scale crustal movement. Mathematical cover, physical cover and manifold element for general finite covers are explained; the usual forms of displacement function on physical cover are presented. The global displacement function is the weighted averages of local independent cover functions on the common part of several covers.3) Based on the simplified relation formula of displacement and strain in the spherical coordinate system, all relative formulae for spherical general finite covers are deduced. The relation equation of the manifold elements on spherical surface is established according to the theory of least Potential Energy. The potential energy expressions are separately worked out, including stiffness stress, initial stress, point loading, body force loading, inertia matrix, velocity matrix and fixed point matrix. By differentiating the above expressions, the correspondent coefficient sub-matrix are deduced and added to the equilibrium equation.4) The kinematics equation of the manifold elements system is established. A local plane projection coordinate system is set up in order to simplify the calculation of spherical invasion distance and so on. Normal contact matrix, shear contact matrix and friction force matrix for general finite covers are deduced in details.
    5) Based on the simplified relation formulation of displacement and strain in the spherical coordinate system, all relative formulations for spherical triangle finite covers are deduced. The relation equation of the manifold elements on spherical surface is established according to the theory of least Potential Energy. The potential energy expressions are separately worked out, including the stiffness stress, initial stress, point loading, body force loading, inertia matrix, velocity matrix and fixed point loading. By differentiating the above expressions, the correspondent coefficient sub-matrix are deduced and added to the equilibrium equation.6) Normal Contact matrix, Shear Contact matrix and friction force matrix for triangle mesh covers are deduced. They help to connect the individual discontinuous boundaries into a system.7) To meet the needs of simulating crustal movement, Spherical Numerical Manifold Method software is developed with VC++ language. Time step based on large displacement analysis, open-close iterations and others are introduced as well.8) Spherical numerical manifold method and Spherical discontinuous deformation analysis are used to calculate models with one and three faults respectively. The results indicate that the SNMM can deal with Continua and discontinuity very well.9) GPS station velocities relative to Euroasia plate in Chuandian district are derived from twice observations of Crustal Movement Observation Network of China and NUVEL1A model. The main faults framework is established according to Chuandian's geological materials. The fields of velocity, maximum and minimum stress change and maximum shear stress change of whole Chuandian district are simulated with spherical numerical manifold
    method. Finally result analysis is fully made.
引文
[1] 安美建、石耀霖、李方全.用遗传有限单元反演法研究东亚部分地区现今构造应力场的力源和影响因素[J].地震学报,1998,vol20,No.3,P225—231
    [2] 陈兵、江在森、王双绪等.利用非连续变形数值方法研究块体运动及其应力场初探,地壳形变与地震,2000,vol20,No.1,P38—42
    [3] 陈兵、江在森、张希、王双绪.块体系统非连续变形数值方法——DDA及其在地壳应力分析中的应用,高原地震,1999,vol11,No.1,P1—5
    [4] 陈兵、江在森、赵振才.中国西部断层形变趋势异常特征与地震活动关系研究,中国地震,2000,vol16,No.1,P77—85
    [5] 陈化然、蒋淳、李轶群等.强震活动有限元模型研究及其初步应用,中国地震,2002,vol18,No.3,P249—258
    [6] 陈开平、马瑾.印度与欧业大陆碰撞构造变形数值分析,地震地质,1995,vol 17,No.3,P277—284
    [7] 曹文贵、速宝玉.流形元覆盖系统自动形成方法之研究,岩土工程学报,2001,vol23,No.2,P187—190
    [8] 程宗颐、朱文耀.APRGP97和APRGP98的GPS联测资料确定亚太地区地壳运动[J].中国科学院上海天文台年刊,2001,No.22,P15—25
    [9] 程万正、刁桂苓、吕弋培、张永久、李桂芳、陈天长.川滇地块的震源力学机制、运动速率和活动方式,地震地质,2003,vol25,No.1,P72—87
    [10] 程万正、杨永林.动态地壳形变与川滇地区成组强震活动,地震研究,2003,vol26,No.6,P62—69
    [11] 蔡永昌、朱合华、夏才初.流形方法覆盖系统自动生成算法,同济大学学报,2004,vol32,No.5,P585—590
    [12] 蔡永昌、张湘伟、骆少明.数值流形方法在连续体数值分析中的应用,力学与实践,1999,vol21,No.6,P53—54
    [13] 陈智梁、张选阳、沈风等.中国西南地区地壳运动的GPS监测[J].验科学通报,1999,vol44,No.8,P851—854
    [14] 邓起东、张培震、杨晓平等.中国活动构造基本特征,中国科学(D辑),2002,vol32,No.12,P1020—1030
    [15] 党亚民.基于反演理论的大地测量形变分析与解释的理论和方法[博士论文],武汉测绘科技大学,1998
    [16] 党亚民、陈俊勇、晁定波.中国及其邻区地壳应力场数值特征的研究[J].测绘科学,2001,vol26,No.2,P11—14
    [17] 党亚民,陈俊勇、刘经南等.利用国家GPSA级网资料对中国大陆现今水平形变场的初步分析,测绘学报,1998,No.3
    [18] 党亚民、李繁麟.川滇GPS地壳形变监测网的观测与数据处理,地壳形变与地震,1996,vol16,No.4,P31—36
    [19] 独知行.基于力学模式的大地测量反演理论及应用[博士论文],中国科学院,2001
    [20] 傅容珊、黄建华、常筱华、曹抗震.东亚大陆形变应力场格局演化的数值模拟,地壳形变与地震,2000,vol20,No.3,P1—10
    [21] 傅容珊、黄建华、徐耀明等.印度与欧亚板块碰撞的数值模拟和现代中国大陆形变[J].地震学报,2000,vol22,No.1,P1—7
    [22] 顾国华、孙汉荣.利用GPS地形资料在大地坐标系中计算应变[J].地壳形变与地震,1998,vol18,No.3,P26—31
    [23] 顾国华,申旭辉、王敏等,中国大陆现今地壳水平运动基本特征[J].地震学报,2001,vol23,No.3,P362—369
    [24] 郭良迁.对中国大陆地壳水平变形的初步探索[J].中国地震,2000,vol16,No.2,P126—134
    [25] 郭良迁、应绍奋.从GPS复测结果初步研究华北北部的水平形变及构造活动[J].中国地震,1998,vol14,No.4,P11—19
    [26] 甘卫军、沈正康、张培震.青藏高原地壳水平差异运动的GPS观测研究,大地测量与地球动力学,2004,vol24,No.1,P29—35
    [27] 胡明城.现代大地测量学的理论及其应用,测绘出版社,2003
    [28] 胡明城、鲁福.现代大地测量学。北京:测绘出版社,1994
    [29] 黄立人、郭良迁.华北(北部)地区的地壳水平运动结征及其模型—GPS测量结果的初步分析[J].地壳形变与地震,1998,vol18,No.2,P20—27
    [30] 黄立人、符养、段五杏等.由GPS观测结果推断中国大陆活动构造边界,地球物理学报,2003,vol46,No.5,P609—615
    [31] 黄立人、王敏.中国大陆现今地壳水平运动,地震学报,2000,vol22,No.3,P257—262
    [32] 黄立人、张培震.中国大陆的GPS速度场、地块运动与应变场,工程地球物理学报,2004,vol1,No.3,P204—212
    [33] 韩渭宾、蒋国芳.川滇地区强震活动分布特征及其与地壳块体构造背景关系的研究,地震学报,2004,vol26,No.2,P211—222
    [34] 何玉梅、姚振兴.中国台湾南部及其周边岛屿现今地壳形变的位错模型.地球物理学报,2002,vol45,No.5,P638—645
    [35] 华中工学院.现代工程数学手册(Ⅰ、Ⅱ、Ⅲ、Ⅳ、Ⅴ),华中工学院出版社
    [36] 焦明若、张国民、刘杰等.中国大陆及其周边地区构造应力场的数值计算及其在地震活动性解释上的应用[J].地震学报,1999,vol21,No.2,P123—132
    [37] 敬少群、吴云、殷志山.中国大陆及其周边体的现时运动、变形及动力学解释[J].地壳形变与地震,2001,vol21,No.4,P8—16
    [38] 江在森、丁平、王双绪等.中国西部大地形变监测与地震预测,地震出版社,2001
    [39] 江在森、马宗晋、牛安福、张晓亮、王双绪、陈兵.GPS技术应用于中国地南运动研究的方法及初步结果,地学前缘(中国地质大学.北京),2003,vol10,No.1,P71—79
    [40] 江在森、王双绪、赵振才.南北地震带和青藏块体东部近期大地形变与地震特征,中国地震,1997,vol13,No.2,P139—149
    [41] 裴觉民.数值流形方法与非连续变形分析,岩石力学与工程学报,1997,vol16,No.3,P279—292
    [42] 李桂芳、程万正.川滇交界地区构造应力场及强震趋势的统计预测,四川地震,1994,No.4,P1—7
    [43] 李国和、王思敬、张强等.川滇交界地区地壳结构及现代地壳活动模式,地质力学学报,vol6,No.2,P82—90
    [44] 李铁明、邓志辉、吕弋培.川滇地区现今地壳形变及其与强震时空分布的相关性研究,中国地震,vol19,No.6,P132—147
    [45] 牛安福、李旭东、郑文衡、柯丽君.丽江地震前后川滇地区地壳形变异常强度场及其演化,地震,1998,vol18,No.3,P303—307
    [46] 牛安福、王琪、江在森等.伽师地震与区域断层活动研究,大地测量与地球动力学,2003,vol23,No.1,P17—21
    [47] 刘本培、吕弋培、李建中、廖华、彭长虹.川滇断块东界中段地区现代地壳形变和断裂现今活动,四川地震,1994,No.4,P53—63
    [48] 吕江宁、沈正康、王敏.川滇地区现代地壳运动速度场和活动块体模型研究,地震地质,2003,vol25,No.4,P543—554
    [49] 吕弋培、廖华、苏琴、王兰.川滇菱形块边界的现今地壳形变,中国地震,2002,vol18,No.1,P28—37
    [50] 刘经南、葛茂荣.92中国GPS会战(A级网)数据处理分析,武汉测绘科技大学学报,1996,vol20,No,1,P40—45
    [51] 刘经南、施闯、陈俊勇.92、96国家高精度GPSA级网整体平差结果分析与我国块体运动模型研究,武汉测绘科技大学学报,1998,vol23,No.4,P314—319
    [52] 刘经南、施闯、许才军等.利用局域复测GPS网研究中国大陆块体现今地壳运动速度场[J].武汉大学学报信息科学版,2001,vol26,No.3,P189—195
    [53] 刘经南、施闯、姚宜斌、陶本藻.多面函数拟合法及其在建立中国地壳平面运动速度场模型中的应用研究,武汉大学学报.信息科学报,2001,vol26,No.6,P500—508
    [54] 刘经南、许才军、宋成骅等.精密全球卫星定位系统多期复测研究青藏高原现今地壳运动与应变[J].科学通报,2002,vol45,No.24,P2658—2663
    [55] 刘经南、姚宜斌、施闯.中国地壳运动整体速度场模型的建立方法研究,武汉大学学报.信息科学版,2002,vol27,No.4,P331—336
    [56] 刘经南、姚宜斌、施闯.中国大陆现今垂直形变特征的初步探讨,大地测量与地球动力学,2002,vol22,No.3,P1—5
    [57] 宁津生,晁定波,金标仁.现代大地测量在研究相关地学问题中的作用,地球物理学进展,1997,No.1
    [58] 刘澜波.用DDA模拟华北地区的无震断层滑动和断块形变,世界地震译丛,2001,No.5,P56—63
    [59] 李延兴、胡新康、帅平等.中国大陆地壳水平运动统一速度场的建立与分析[J].地震学报,2001,vol23,No.5,P453—459
    [60] 骆少明、蔡永昌、张湘伟.数值流形方法中的网格重分技术及其应用,重庆大学学报(自然科学版),2001,vol24,No.4,P34—41
    [61] 骆少明、张湘伟、姜东茹.三维数值流形的面向对象设计,重庆大学学报,2003,vol26,No.12,P10—14
    [62] 骆少明、张湘伟、蔡永昌.数值流形方法的变分原理与应用,重庆出版社,2001,vol22,No.6,P587—592
    [63] 骆少明、张湘伟、蔡永昌.非线性数值流形方法的变分原理与应用,重庆出版社,2000,vol21,No.12,P1265—1270
    [64] 赖锡安、黄立人、徐菊生等.中国大陆现今地壳运动,地震出版社2004
    [65] 赖锡安、徐菊生、晁定波,利用GPS基线变化预测滇西地区的地震危险性[J].震情研究,1993,vol19,No.3,P35—40
    [66] 马杏垣(主编).中国岩石圈动力学地图集,北京:中国地图出版社,1989
    [67] 马宗晋、陈鑫连、叶叔华等.中国大陆区现今地壳运动的GPS研究,科学通报,2001,vol46,No.13,P1118—1120
    [68] 牛之俊、马宗晋、陈鑫连等.中国地壳运动观测网络,大地测量与地球运动力学,2002,vol22.No.3.P1—7
    [69] 裴锡瑜.晚第四纪安宁河活动断裂了分段的基本特征,四川地震,1998,vol4,P52—61
    [70] 彭自强、葛修润.数值流形方法在有限元三维二十结点单元上的实现,岩石力学与工程学报,2004,vol23,No.15,P2622—2627
    [71] 彭自强、王水林、葛修润.单位分解法、无网格法、数值流形方法这形函数的内在联系,岩石力学与工程学报,2002,vol21,No.2,P2429—2431
    [72] 秦小军、刘宏先、赵齐乐等.断层带及其周围块体系统的非连续运动变形模拟,地壳形变与地震,2001,vol21,No.2,P17—24
    [73] 乔学军、王琪、杜瑞林.川滇地区活动地块现今地壳形变特征,地球物理学报,2004,vol47,No.5,P805—811
    [74] 乔学军、王晓强、王琪等.新疆伽师地区现今构造运动和地形变GPS监测的初步研究[J].地壳形变地震,2000,vol20,No.3,P49—56
    [75] 秦小军、周硕愚、乔党军等.依据GPS数据研究伽师及帕米尔东北侧1994~1998年地壳水平运动及其与强震群的关系,地壳形变与地震,2001,vol21,№.4,P1—7
    [76] 任金卫.利用GPS观测结果反演中国大陆及邻区现代构造变形运动[J].地球物理学报,2002,vol45,(增刊),P123—135
    [77] 申重阳、王琪、吴云、游新兆、甘家思、杨少敏.川滇菱形块体主要边界运动模型的GPS数据反演分析,地球物理学报,2002,vol45,No.3,P352—361
    [78] 石根华.块体系统不连续变形分析[M].任放等译,北京:科学出版社,1993
    [79] 石根华、裴觉民.数值流形方法与非连续变形分析,清华大学出版,1997
    [80] 帅平、吴云、周硕愚.用GPS测量数据模拟中国大陆现今地壳水平速度场及应变场[J].地壳形变与地震,2001,vol21,No.4,P1—8
    [81] 苏有锦、晏凤桐.滇西地区强震前兆源过程的力学模型与数值摸拟,地震学报,1996,vol18,No.3,P309—311
    [82] 苏有锦、刘祖荫、蔡民军.去南地区强震分布的深部地球介质背景,地震学报,1999,vol21,No.3,P313—322
    [83] 时振梁.我国强震活动和板内构造,地质科学,vol4,1973
    [84] 王凯英、马瑾、顾国华、刁桂苓.1996年川滇地区发生的一次应力转向事件,地学前缘(中国地质大学.北京),2003,vol10,P233—239
    [85] 王敏、沈正康、甘卫军等.现今中国大陆地壳运动与活动块体模型,中国科学(D辑),2003,vol33,增刊,P21—32
    [86] 王琪.中国大陆现今地壳运动研究,地震学报,2003,vol25,No.5,P541—547
    [87] 王琪,丁国瑜、乔学军等.用GPS研究南天山(伽师)地区现今地壳变形[J].地震学报,2000,vol22,No.3,P263—270
    [88] 王琪、赖锡安、游新兆、王启梁.红河断裂的GPS监测与现代构造应力场[J].地壳形变与地震,1998,vol18,No.2,P49—56
    [89] 王琪、游新兆、王文颖等.跨喜马拉雅的GPS观测与地壳形变[J].压器地壳形变与地震,1998,vol18,No.3,P43—50
    [90] 王琪、张培震、马宗晋.中国大陆现今构造变形GPS观测数据与速度场[J].地学前缘,2002,vol9,No.2,P416—429
    [91] 王琪、张培震、牛之俊、J.T.Freymueller等.中国大陆现今地壳运动和构造变形,中国科学(D辑),2001,vol31,No.7,P529—536
    [92] 王仁.我国地球动力学的研究进展与展望,地球物理学报,1997,vol40,增刊
    [93] 王仁、丁中一、殷有泉.固体力学基础,地震出版社,1979
    [94] 王仁、梁海华.用叠加法反演东亚地区现代构造应力场,见:马宗晋主编,国际交流地质学术论文集2—为二十七届国际地质大会撰写,北京:地震出版社,1985,P29—35
    [95] 王书法.岩体不连续非线性变形分析的数值流形方法研究[博士论文],中国科学院,2000
    [96] 王绍晋、龙晓帆.云南三次强震前后环境剪应力场动态变化,地震学报,1999,vol21,No.2,P150—158
    [97] 王绍晋、秦嘉政、关蕾等.云南强震活动的多层次动力源分析,地震研究,2003,vol26,No.2,P142—150
    [98] 王绍晋、张俊伟、李中华等.云南及邻区环境剪应力场时间分布与强震活动,地震研究,2004,vol27,No.2,P165—172
    [99] 汪素云、陈培善.中国及邻区现代构造应力场的数值模拟,地球物理学报,1980,vol23,No.1,P35—45
    [100] 王勖成等.有限单元法基本原理和数值方法,清华大学出版社,1997
    [101] 王小亚、朱文耀等.GPS监测的中国及其周边现时地壳形变,地球物理学、报,2002,vol45,No.2,P198—209
    [102] 王勇、许厚泽.青藏高原印度板块向欧亚大陆俯冲速率的研究-GPS观测资料的反演结果.地球物理学报,2003,vol46,No.2,P185—190
    [103] 吴云、王辉、周硕愚、施顺英.球面上的块体系统非连续变形反分析法,大地测量与地球动力学,2002,vol22,No.3,P15—18
    [104] 王泽民.非连续变形分析与现代地壳运动[博士论文],武汉大学,2001
    [105] 王泽民.非连续变形分析与现代地壳运动研究,武汉大学学报.信息科学版,2001,vol26,No.2,P122—126
    [106] 王泽民、伍岳、刘经南、崔天鹏.利用GPS资料采用非连续变形分析模型确定中国大陆地壳运动速度场,武汉大学学报.信息科学报,2002,vol27,No.3,P241—245
    [107] 王泽民、刘经南.球面坐标系非连续变形分析(DDA)的数学模型,武汉大学学报(信息科学版),2001,vol26,No.3
    [108] 王芝银、李云鹏.数值流形方法及其研究进展,力学进展,2003,vol33,No.2,P261—266
    [109] 王芝银、李云鹏.数值流形方法中的几点改进,岩石力学与工程学报,1998,vol20,No.6,P33—34
    [110] 王书法、朱维申、李术才、陈胜宏.岩体弹塑性分析的数值流形方法,岩石力学与工程学报,2002,vol21,No.6,P900—904
    [111] 许才军.青藏高原地壳运动模型与构造应力场,测绘出版社,2001
    [112] 许才军.大地测量联合反演理论和方法研究进展,武汉大学学报.信息科学版,2001,vol26,No.6,P555—561
    [113] 许才军、董立祥、李志才.华北地区壳形变的GPS及地震矩张量反演分析[J].武汉测绘科技大学学报,2000,vol25,No.6,P471—475
    [114] 许才军、董立祥、施闯等.华北地区GPS地壳应变能密度变化率场及其构造运动分析[J].地球物理学报,2002,vol45,No.4,P497—506
    [115] 许才军、晁定波.用大地测量资料反演青藏高原构造应力场的初步尝试[J].测绘学报,1997,vol26,No.2,P95—100
    [116] 徐菊生、赖锡安、卓力格图等.中国大陆西部地区现今的块体运动GPS监测的初步研究[J].地壳运动与地震,2000,vol20,No.2,P15—23
    [117] 熊熊、滕吉文、郑勇.中国大陆地壳运动的GPS观测及相关动力学研究,地球物理学 进展,2004,vol19,No.1,P16—25
    [118] 徐锡伟、闻学泽、郑荣章、马文涛、宋方敏、于贵华.川滇地区活动块体最新构造变动样式及其动力来源,中国科学(D辑),2003,vol33,P151—162
    [119] 许忠淮、汪素云、黄雨蕊等.由大量的地震资料推断的我国大陆构造应力场,地球物理学报,1989,vol32,No.6,P636—647
    [120] 尹京苑、梅世蓉.地壳应力状态及其在三维模拟计算中的意义,地震,1998,vol18,No.3,P226—232
    [121] 尹祥础.固体力学基础,地震出版社,1985
    [122] 杨红心、王乘、张秋文.青藏高原隆升的非线性动态有限元仿真研究[J].大地测量与地球动力学,2002,vol22,No.2,P75—81
    [123] 杨凯.中国测绘学科发展蓝皮书,测绘出版社,2003
    [124] 叶叔华、赖锡安、刘经南.现代地壳运动与地球动力学研究,中国大陆主要活动带现今地壳运动及动力学研究,地震出版社,2001
    [125] 杨少敏、游新兆、杜瑞林等.用双三次样条函数和GPS资料反演现今中国大陆构造形变场,大地测量与地球动力学,2002,vol22,No.1,P68—75
    [126] 游新兆、杜瑞林、王琪等.中国大陆地壳现今运动的GPS测量结果与初步分析,地壳形变与地震,2001,vol21,No.3,P1—8
    [127] 游新兆、王启梁、吴星华.青藏高原1993年GPS观测成果的精度分析[J].地壳形变与地震,1994,vol14,No.3,P27—33
    [128] 姚宜斌、刘经南、施闯、陶本藻.ITRF97参考框架下的中国大陆区域地壳板块运动背景场的建立及应用,武汉大学学报.信息科学报,2002,P363—376
    [129] 殷有泉.有限单元方法及其在地学中的应用,地震出版社,1987,
    [130] 张崇立、马宗晋.我国地震的现今地球动力学研究的进展与方向,中国地震,1995,vol11,No.1,P38—48
    [131] 张大林、栾茂田、杨庆、杨新辉.数值流形方法的网格自动剖分技术及其数值方法,岩石力学与工程学报,2004,vol23,No.11,P1836—1840
    [132] 张东宁、许忠淮.中国大陆岩石层动力学数值模型的边界条件,地震学报,1999,vol21,No.2,P133—139
    [133] 张东宁、许忠淮.青藏高原现代构造应力状态及构造运动的三维弹粘性数值模拟[J].中国地震,1994,vol10,No.2,P136—143
    [134] 张东宁、许忠淮.中国大陆大岩石层动力学数值模型,地震学报,1999,vol21,No.2,P133—139
    [135] 张培震、王琪、马宗晋.中国大陆现今构造运动的GPS速度场与活动地块[J].地学前缘,2002,vol9,No.2,P431—441
    [136] 张培震、王琪、马宗晋.青藏高原现今构造变形特征与GPS速度场[J].地学前缘,2002,vol9,No.2,P442—450
    [137] 张强、朱文耀、程宗颐等.中国地壳各构造块体运动模型的初建[J].科学通报,2000,vol45,No.9,P960—967
    [138] 赵齐乐.构造块体及边界带系统现今年运动及变形的数值模拟研究[硕士论文],中国地震局,2001
    [139] 赵齐乐、秦小军、乔学军等.伽师地区强震群前后地壳运动变形的数值模拟,地壳形变与地震,2000,vol20,No.4,P29—35
    [140] 朱守彪.球面有限单元法分析与应用,防灾技术高等专科学校学报,2004,vol6,No.1,P1—7
    [141] 朱守彪、石耀霖.用遗传有限单元法反演川滇下地壳流动对上地壳的拖曳作用,地球物理学报,2004,vol47,No.2,P232—239
    [142] 赵森林、胡建国、李毓麟.川滇断块的现代地壳运动,测绘通报,1998,、P4—6
    [143] 赵少荣.动态大地测量反演及物理解释的理论与应用,博士论文,武汉测绘科技大学,1991
    [144] 赵少荣.基于力学模式的动态大地测量数据反演研究[J].测绘学报,1994,vol23,No.2,P90—97
    [145] 臧绍先、吴忠良、宁杰远.中国周边板块的相互作用及其对中国应力场的影响——(Ⅱ)印度板块的影响,地球物理学报,1992,vol35,No.4,P428—440
    [146] 朱文耀、程宗颐、王小亚等.中国大陆地壳运动的初步研究[J].院中国科学(D辑),1997,vol27,No.5,P385—389
    [147] 朱文耀、程宗颐、熊永清.中国大陆地壳运动的背景场[J].科学通报,1999,vol44,No.14,P1537—1540
    [148] 周硕愚.走向21世纪的地壳形变学——对大陆动力学与地震预测的新推动,地壳形变与地震,1999,vol19,No.1,P2—13
    [149] 周硕愚.板内地壳形变系统动力学与地震预报,地壳形变与地震,1990,vol10,No.4,P1—7
    [150] 周硕愚、吴云、秦小军等.基于多种GPS数据研究福建及其邻近海域199421997年地壳水平运动.地球物理学报,2000,vol43,No.4,P471—479
    [151] 周硕愚、吴云、王若柏等.地震孕育过程中的地壳形变场图像动力学参量的研究,地震学报,1994,vol16,№.3,P336—340
    [151] 周硕愚、张跃刚、丁国瑜等.依据GPS数据建立中国大陆板内块体现时运动模型的初步研究,地震学报,1998,vol20,№.4,P347—355
    [152] 赵少荣.基于力学模式的动态大地测量数据反演研究[J].测绘学报,1994,vol23,No.2,P90—97
    [153] 臧绍先、吴忠良、宁杰远.中国周边板块的相互作用及其对中国应力场的影响——(Ⅱ)印度板块的影响,地球物理学报,1992,vol35,No.4,P428—440
    [154] 朱文耀、程宗颐、王小亚等.中国大陆地壳运动的初步研究[J].院中国科学(D辑),1997,vol27,No.5,P385—389
    [155] 朱文耀、程宗颐、熊永清.中国大陆地壳运动的背景场[J].科学通报,1999,vol44,No.14,P1537—1540
    [156] 朱文耀、王小亚、程宗颐等.利用GPS技术检测中国大陆地壳运动的初步结果[J].中国科学(D辑),2000,vol30,No.4,P394—400
    [157] 张湘伟、蔡永昌、骆少明.数值流形方法的对象设计,计算力学学报,2001,vol.18,No.1,P8—14
    [158] 张晓亮、江在森、王双绪等.青藏块体东北缘及其周围地区现今时空运动变形特征,中国地震,2003,vol19,No.4,P377—385
    [159] 张希、江在森、王琪、王双绪、张晓亮.川滇地区地壳水平运动特征与强震关系研究,大地测量与地球动力学,2003,vol23,No.3,P35—41
    [160] Abdrakhmatov K Ye, Aldazhanow S A, Hager B H, Hamburger M W. Relatively recent construction of the Tien Shah inferred from GPS measurements of present day crustal deformation rates [J]. Nature, 1996, vol384, P450—453
    [161] Armijo R, Tapponnier P, Mercier J L, Han JT, Han TL. Quaternary extension in southern Tibet: Field observations and tectonic implications[J]. J Geophys Res., 1986, vol91, P13803—13872
    [162] Bird P. Thin plate and thin shell finite element programs for for2 ward dynamic modeling of plate deformation and faulting [J]. Computers & Geosciences, 1999, vol25, P383—394
    [163] Ball A W. Dynamics of plate interiors [M]. American Geophysical Union, Geological Society of America. 1980.
    [164] Bird P. New finite element technique for modeling deformation histories of continents of continents with stratified temperature dependent rheology [J]. J Geophys Res., 1989, vol94, P3967—39990
    [165] Bendick R, Bilham R, Freymueller J, Larson K, Yin G H. Geodetic evidence for a low slip rate in the Altyn Tagh fault sys2tem [J]. Nature, 2000, vol404, P69 —72
    [166] Bilham R, Larrson K, Freymueller J, Project Iedylhim members. GPS measurements of present day convergence across the Nepal Himalaya [J]. Nature, 1997, vol386, P61—64
    [167] Chen Z, Burchfiel B C, Liu Y, King R W, Royden, L H, Tang W, Wang E, Zhao J, Zhang X. Global Positioning System measurements from eastern Tibet and their implications for India / Eurasia intercontinental deformation [J]. J Geophys Res, 2000, voll05, No. 7 P16215—16227
    [168] England P, McKenzie D P. A thin viscous sheet model for continental deformation [J]. Geophys. J R Astron Soc., 1982, vol70, P295—321
    [169] England P, Molnar P. Active deformation of Asia: From kinematics to dynamics [J]. Science, 1997, vol278, P647—650
    [170] England P, Molnar P. The field of crustal velocity in Asia calcu lated from Quaternary rates of slip on faults [J]. Geophys. J. Int., 1997, voll30, P551 —582
    [171] Flesh L M, Haines AJ, Holt W E. Dynamics of the India Eurasia collision zone [J]. J Geophys Res., 2001, voll06, No. 8, P16435—16460
    [172] D. M. Doolin, N. Sitar. DDAML — discontinuous deformation analysis markup language, International Journal of Rock Mechanics and Mining Sciences, 2001, P1—8
    [173] HAN Wei-bin、 JIANG Guo-fang. Study on distribution characteristics of strong earthquakes in Sichuan-Yunnan area and their geological tectonic background, ACTA SEISMOLOGICA SINICA, voll7, No. 2, P230—243
    [174] Haines A J. Calculating velocity field across plate boundaries from observed shear rates from observed shear rates [J]. Geophys. J. R . Astro. Soc. London, 1982, vol68, P203—209
    [175] Haines A J, Holt W E. A procedure for obtaining the complete horizontal motions within zones of distributed deformation from the inversion of strain fates data [J]. J Geophys Res, 1993, vol98, No. 7, P12057—12082
    [176] Holt W E, Chamot Rooke N, Le Picheon X, Haines AJ. Velocity field in Asia inferred from Quaternary fault slip tares and Global Positioning System observations [J]. J Geophys Res. , 2000, voll05 (B8), P19,185—19,209
    [177] Hager B H, King, R W, Murray M H. Measurement of crustal deformation using the Global Positioning System [J]. Annu. Rev. Earth Planet. Sci., 1991, voll9, P351 —382
    [178] Holt W E, Haines AJ. The kinematics of northern South Island, New Zealand, determined from geological strain rates [J]. J. Geophys. Res., 1995, vol100, P17991—18010
    [179] Holt W E, Haines AJ. Velocity fields in deforming Asia from the inversion of earthquake2released strains [J]. Tectonics, 1993, voll2, P1—12
    [180] Jackson J A, McKenzie D P. The relationship between plate motions and seismic moment tensors, and the rates or active deformation on Mediterranean and Middle East [J]. Geophys. J R Astron Soc. , 1988, vol93, P45—73
    [181] Jeen-Shang Lin, Member ASCE and Cheng-Yu Ku. SIMULATION OF SLOPE FAILURE USING A MESHED BASED PARTITON OF UNITY METHOD, 15th ASCE Engineering Mechanics Conference, 2002, P1—6
    [182] Kostov V V. Seismic moment, energy of earthquakes, and the seismic flow of rock [J]. IZV Acad Sci. USSR Phys. Solid Earth, Engl. Transl, 1974, vol10, P23— 44
    [183] Kong X H, Bird P. SHELLS: A thin2shell program from modeling neotectonic of regional or global lithosphere with faults [J]. J Geophys Res. , 1995, vollOO, P22, 129—22, 131
    [184] Larson K M, Burgmann R, Bilham R, Freymueller J T. Kinematics of the India2Eurasia collision zone from GPS measurements [J]. J Geophys Res. , 1999, vol104, No. 1, P1077—1093
    [185] Larson KM, Freymueller J T, Phiipen S. Global plate velocities form the Global Positioning System [J]. J Geophys. Res., 1997, voll02, P9951—9981
    [186] LUO Shao-ming, ZHANG Xiang-wei, CAL Yong-china. THE VARIA TIONAL PRINCIPLE AND APPLICATION OF NUMERICAL MANIFOLD METHOD, Variational Principle of Numerical Manifold Method, 2001, vol22, No. 6, P658—663
    [187] M. MacLaughlin, N. Sitar, D. Doolin, T. Abbot. Investigation of slope-stability kinematics using discontinuous deformation analysis, International Journal of Rock Mechanics and Mining Sciences, 2001, P1—10
    [188] Molnar P, Lyon2Caen H. Fault plane solution of earthquakes and active tectonics of the Tibetan Plateau and its margins [J]. Geo2phys J Int., 1989, vol99, P123—153
    [189] N. Bicanic, C. Stirling, C. J. Pearce. Discontinuous Modelling of Structural Masonry, Fifth World Congress on Computational Mechanics, 2002, P1—18
    [190] Peltzer G, Saucier F. Present day kinematics of Asia derived from geological fault rates [J]. J Geophys Res., 1996, vollOl, No. 12, P27943—27956
    [191] Qin C Y, Papazachos C, Papadimitriou E. Velocity field for crustal deformation in China derived from seismic moment tensor summation of earthquakes [J]. Tectonophys. , 2002, vol359, P29—46
    [192] S. Li. Cheng. Y. -F. Wu. Numerical manifold method based on the method of weighted residuals, 2004
    [193] S hen Z K, Zhao C, Li Y, Jackson D, Fang P, Dong D, Yin A. Contemporary crustal deformation in east Asia constrained by global positioning system measurements [J]. J. Geophys. Res., 2000, voll05, P5721—5734
    [194] Shen Z K, Wang M, Li YX. Crustal deformation along the AltTagh fault system, western China from GPS [J]. Geophys. Res., 2001, voll06, No. 12, P30607—30621
    [195] Tapponnier P, Molnar P. Active faulting and tectonics of China [J]. J Geophys Res., 1977, vol82, P2905—2930
    [196] Wang Q, Zhang P Z, Freymueller J T, Biham R, Larson KM, Lai X A, You X Z. Niu ZZ, Wu J C, Li Y X, Liu J N, Yang Z Q, Chen Q Z. Present2day crustal deformation in China constrained by Global Positioning System measurements [J]. Science, 2001, vol294, P574—577
    [197] ZHANG Ming, YANG Heqing, LI Zhongkui. A Coupling Model of the Discontinuous Deformation Analysis Method and the Finite Element Method, tsinghua science and technology, 2005, vol10, No. 2, P221—226

© 2004-2018 中国地质图书馆版权所有 京ICP备05064691号 京公网安备11010802017129号

地址:北京市海淀区学院路29号 邮编:100083

电话:办公室:(+86 10)66554848;文献借阅、咨询服务、科技查新:66554700