GPS地壳形变监测网应变解算方法与程序设计  

METHODS AND PROGRAM DESIGN OF STRAIN COMPUTATION FROM GPS DEFORMATION DATA

作　　者：孟国杰[1,2] 申旭辉[3] 伍吉仓[1] 沈云中[1] 张秋文[4] 

机构地区：[1]同济大学测量与国土信息工程系 [2]中国地震局地震预测研究所,北京100036 [3]中国地震局地震预测研究所 [4]华中科技大学水电与数字化工程学院

出　　处：《大地测量与地球动力学》2006年第2期70-74,共5页Journal of Geodesy and Geodynamics

基　　金：中国博士后科学基金(20040350495);地震科学联合基金(605032);科技部重大基础研究前期研究专项(2003CCB00600)

摘　　要：分析了在平面直角坐标系和大地坐标系内解算应变的方法及其特点,并以Matlab为开发平台,研制出一套在平面直角坐标系和大地坐标系中解算应变的程序GPSstrain。利用GPS监测网测站位移(速率),该程序能够解算剪切应变(γ1、γ2)、面应变、块体旋转、主应变及其方向、最大剪切应变及其方向,同时也可解算出这些应变分量的误差。另外该程序还具有高斯投影解算功能,根据各测站坐标能够自动形成所有的Delaunay三角形,并以图形方式显示。最后以实际算例给出了该解算程序的部分结果,并进行了初步的分析和讨论。This paper analyzes the ordinate system and their characteristics methods for calculating strains in plane coordinate system and geodetic corespectively. We have designed a software, GPS strain based on the Matlab platform for computing strain in the two coordinate systems, which can be used to compute shear strain components T1 ,T2, area strain, block rotation, principal strains and their directions, maximum shear strain and its direction, and the uncertainties of all the components can be computed as well. Besides, Gaussian projection is included in the software package, and thus Delaunay triangles can be formed automatically from all GPS geodetic coordinates, and displayed in graphics mode. Taking the GPS station velocities in northeastern China as an example, this paper analyzes and discuses preliminarily the parts of computational results.

关 键 词：GPS 地壳形变 应变场分析 应变解算方法 程序设计 

分 类 号：P227[天文地球—大地测量学与测量工程]