广义应变花方法及其在川滇地区地应变率场计算中的应用

doi:10.7523/j.issn.2095-6134.2005.3.006

摘要/Abstract

摘要：

提出将工程测量中的平面应变花分析方法推广到地球球面上的广义应变花方法;并通过对已知位移和应变场的特例下的数值实验,比较常用的最小二乘配置法和广义应变花方法计算应变的精度.结果显示广义应变花法应变的计算精度比最小二乘配置法高.最后以中国川滇地区为例,利用GPS观测资料计算了地应变率场,并对计算方法和结果进行了分析讨论.

关键词: 广义应变花方法, 最小二乘配置法, GPS观测, 地应变场, 川滇地区

Abstract:

GPS (global positioning system) velocity can be used to calculate the horiwntal strain rate，but it should be done with great care because digital differentiation is involved，which brings about a prominent problem in minimizing the error. In this paper,a general strain rosette method(GSRM) is proposed，in which the plane strain rosette method in engineering is widened to include spherical earth surface. Then the GSRM is compared with the least squares collocation method (LSCM) in specific case by a digital experiment in which displacements and the strain field are known. The result shows that the accuracy of strain rate computed by GSRM is higher than that of L SRM.Finally，with Chinese Sichuan-Yunnan area as an example,the strain rate is calculated by GPS data and some discussions of computation methods and results are made.

Key words: GSRM, LSCM, GPS measurements, geo-strain rate field, Sichuan-Yunnan area

中图分类号:

P541

朱守彪, 石耀霖. 广义应变花方法及其在川滇地区地应变率场计算中的应用[J]. 中国科学院大学学报, 2005, 22(3): 292-302.

ZHU Shou-Biao, SHI Yao-Lin. General Strain Rosette Method and its Application to the Computation of Geo-Strain Rate in Chinese Sichuan-Yunnan Area[J]. , 2005, 22(3): 292-302.

参考文献

[1] Savage JC,Gan WJ，Svarc JL. Strain accumulation and rotation in the Eastern California Shear Zone. J. Geophys. Res.，2001，106(B10) :21995 ～22007

[2] Murray JR， Segall P， Cervelli P， et al. Inversion of GPS data for spatially variable slip2rate on the San Andreas Fault near Parkfield， CA. Geophys. Res. Lett.， 2001， 28(2) :359～362

[3] DeMets C, Gordon RG, Argus DF， et al. Current plate motions. Geophys. J. Int.，1990， 101(3) :425～478

[4] Freymueller J，Bilham R，Bürgmann R， et al. Global Positioning System measurements of Indian plate motion and convergence across the lesser Himalaya. Geophys. Res. Lett.，1996，23(22) :3107～3110

[5] Larson KM, Freymueller J，Philipsen S. Global plate velocities from the Global Positioning System. J. Geophys. Res.， 1997， 102(B5) :9961～ 9982

[6] Larson KM,Bürgmann R，Bilham R， et al. Kinematics of the India2Eurasia collision zone from GPS measurements. J. Geophys. Res.，1999，104 (B1) :1077～1094

[7] Huang LR. QUAD method used for identifying relatively stable stations. Journal of Geodesy and Geodynamics. 2002，22(2) :10～15(in Chinese with English abstract)

[8] Shen ZK, David DJ， Bob XG. Crustal deformation across and beyond the Los Angeles basin from geodetic measurements. J. Geophys. Res.， 1996，101(B12) :27957～27980

[9] Kato T,El2Fiky GS, Oware EN. Crustal strains in the Japanese islands as deduced from dense GPS array. Geophys. Res. Lett.，1998，25(18) : 3445～3448

[10] Gu GH, Sun HR， Sun HJ， et al. Deduction of strain from GPS crustal deformation data in geodetic coordinate sysem. Crustal Deformation and Earthquake，1998，18(3) :26～31 (in Chinese with English abstract)

[11] Huang LR， WangM. Relative movement and strainof tectonic block. Crustal Deformation and Earthquake， 1999，19(2) :17～26 (in Chinese with English abstract)

[12] El2Fiky GS, Kato T,Oware EN. Crustal deformation and interpolate coupling in the Shikoku district， Japan， as seen from continuous GPS observation. Tectonophysics， 1999， 314 : 387～399

[13] Zhang X, Jiang ZS, Chen B， et al. Characteristics of recent horizontal movement and strain2stress field in the crust of north China. Chinese J. Geophys.，2000，43(5) :657～665(in Chinese with English abstract)

[14] Shen ZK, Zhao CK, Li YX, et al. Contemporary crustal deformation in east Asia constrained by global positioning system measurements. J. Geophys. Res.， 2000，105(B3) :5721～5734

[15] Wang Q， Zhang PZ, Freymueller J， et al. Present2day crustal deformation in China constrained by GPS measurements. Science， 2001， 294: 574 ～577

[16] Wahr J，Dam TV， Larson K, et al. Geodetic Measurements in Greenland and Their Implications. J. Geophys. Res.，2001，106:16567～16581

[17] HuJC, Yu SB. Active deformationof Taiwanfrom GPSmeasurements and numerical simulations. J. Geophys. Res.，2001，106(B2) :2265～2280

[18] Chen J， Tao BZ. Geodetic Deformation Survey. Beijing: Seismological Press， 1980 (in Chinese)

[19] Yang GH. A description and realization of horizontal deformation field by the multi2kernel function. Crustal Deformation and Earthquake，1992，12 (1) :45～51 (in Chinese with English abstract)

[20] Zhang X, Jiang ZS,Zhang SX. Whole calculation of the crustal visual strain field with the least square collocation. Crustal Deformation and Earthquake，1998，18(2) :57～62 (in Chinese with English abstract)

[21] Zhang X, Jiang ZS. Study on some questions of dynamic pictures of crustal deformation and strain field obtained by the least square collocation. Crustal Deformation and Earthquake， 1999，19(3) :32～39 (in Chinese with English abstract)

[22] Moritz H. Interpolation and predication of point gravity anomalies. Annals Academiace Scientiarum Fennicae，1963，69(ⅢA) :1～32

[23] Moritz H. High Physics Geodetic Survey.NingJS, et al. (Trans.) Beijing: Survey Publishing House，1984 (in Chinese)

[24] Fujii Y, Xia S. Estimation of distributionof the ratesof vertical crustal movement rates in the Tokai district with the aidof least squares prediction. J. Phys. Earth， 1993，41:239～256

[25] Yang GH, Xie JM, Han YP. Current horizontal deformation and kinematic mechanismof major tectonic units and boundary zones in north China. Chinese J. Geophys.，2001，44(5) :645～653 (in Chinese with English abstract)

[26] MA ZJ，et al. Study on the Recent Deformation and Dynamics of the Lithosphere of Qinghai2Xizang Plateau. Beijing: Seismological Press， 2001

[27] Li YX, Yang GH, Li Z,et al. Movement and strain conditions of active blocks in the Chinese mainland. Science in China， 2003(supp.)，46:82 ～117

[28] Jiang ZS, Zhang X, Chen WS,et al. Studyon scale dependence of strain value obtained from crustal deformation data. Acta Seismologica Sinica， 2000，22(4) :352～359

[29] Chinese lithospheric kinetics charts editorial board of CSB (MA XY, chief editor). Lithosphetic Kinetics Charts of Chinese Mainland and Its Adjoining Sea Area. Beijing: Geological Publishing House， 1986 (in Chinese)

[30] Zhu SB， Shi YL. Genetic algorithm2finite element inversion of drag forces produced by the lower crust flow to the upper crust in Sichuan2Yunnan area. Chinese J. Geophys.，2004，47(2) :258～267

[31] Huang LR， Ma Q. Statistical test for tectonic location of GPS stations. Crustal Deformation and Earthquakes，1999， 19(4) :27～31 (in Chinese with English abstract)

[32] Shi YL， Zhu SB. Method for division of present active crustal blocks by GPS survey data. J. Geodesy and Geodynamics，2004， 24(2) :1～5(in Chinese with English abstract)

[7] 黄立人. 用于相对稳定点组判别的 QUAD 法. 大地测量与地球动力学，2002，22(2) :10～15

[10] 顾国华，孙汉荣，孙惠娟，等. 利用GPS 地形变资料在大地坐标系中计算应变. 地壳形变与地震，1998，18(3) : 26～31

[11] 黄立人，王敏. 构造块体的相对运动和应变. 地壳形变与地震，1999，19(2) :17～26

[13] 江在森，张希，陈兵，等. 华北地区近期地壳水平运动与应力应变场特征. 地球物理学报，2000，43(5) :657～665

[18] 陈健，陶本藻. 大地形变测量学. 北京:地震出版社，1980

[19] 杨国华. 水平形变场的函数描述及其实现方法. 地壳形变与地震，1992，12(1) :45～51

[20] 张希，江在森，张四新. 借助最小二乘配置整体解算地壳视应变场. 地壳形变与地震，1998，18(2) :57～62

[21] 张希，江在森. 用最小二乘配置获得地形变场动态图象的几个问题研究. 地壳形变与地震，1999， 19(3) :32～39

[23] [奥]赫尔墨特?莫里兹著. 宁津生、管泽霖等译. 高等物理大地测量学. 北京:测绘出版社，1984

[25] 杨国华,谢觉民,韩月萍. 华北主要构造单元及边界带现今水平形变与运动机制. 地球物理学报，2001，44(5) :645～653

[26] 马宗晋等主编. 青藏高原岩石圈现今变动与动力学研究，北京:地震出版社，2001，8

[28] 江在森，张希，陈文胜，等. 地形变资料求解应变值的尺度相对性问题研究. 地震学报，2000，22(4) :352～359

[29] 国家地震局中国岩石圈动力学图集编委会(马杏垣主编). 中国及邻近海域岩石圈动力学图. 北京:地质出版社，1986

[31] 黄立人，马青. GPS测站所处构造位置的统计检验. 地壳形变与地震，1999，19(4) :27～31

[32] 石耀霖，朱守彪. 利用 GPS观测划分现今地壳活动块体的方法及其应用. 大地测量与地球动力学， 2004， 24(2) : 1～5