球坐标系中图形单元应变与旋转张量及其误差解算 - 中国地震信息网.pdf

维普资讯 第 3O卷 第 1期 地 震 学 报 Vo1．30，NO．1 2008年 1月 (67-75) ACTA SEISM OLOGICA SINICA Jan．，2008 孟 国杰 ，任金卫 ，伍吉仓 ，申旭辉．2008．球坐标系 中图形单元应变与旋转张量及其误差解算．地震学报 ，30(1)：67—75． MengGuojie，RenJinwei，WuJicang，ShenXuhui．2008．Computationofstrainandrotationtensorandtheiruncertainty forsmallarraysinsphericalcoordinatesystem．ActaSeismologicaSinica，30(1)：67-75． 球坐标系中图形单元应变与旋转 张量及其误差解算 孟国杰”’’ 任金卫 ¨ 伍吉仓 申旭辉 1)中国北京 100036中国地震局地震预测研究所 2)中国上海 200092同济大学测量与国土信息工程系 摘要 利用泰勒级数展开和弹性力学的几何方程 ，推导了球坐标系中由GPS位移数据解算 图形单元应变和旋转张量 的解析公式．通过线性化处理 ，并利用误差传播定律 ，详细推导 了 图形单元应变和旋转张量 的误差公式．利用川滇地 区最新 的GPS测站位移速率数据 ，采用 图 形单元法解算 了该区的面应变率及最大剪切应变率分布 ，并对其进行 了初步的分析．同时阐 述 了图形单元应变数学模型的局限性 ，讨论 了图形单元法计算应变的意义以及 图形单元 的选 择 问题 ，进一步分析 了应变计算 的定权方法 ，还讨论 了GPS观测 网图形单元 的地心半径差异 与应变 的关系 问题． 关键词 球坐标系 图形单元 应变分析 误差分析 文章编号 ：0253—3782(2OO8)01—0067—09 中图分类号 ：P315．727 文献标识码 ：A Computation ofstrainandrotationtensorandtheir uncertaintyforsmallarraysinspherical coordinatesystem MengGuojie’’ RenJinwei WuJicang ShenXuhui’ 1)InstituteofEarthquakeScience，ChinaEarthquakeAdministration，Beijing100036，China 2)DepartmentofSurveyingandGeo-informatics，TongjiUniversity，Shanghai200092，China Abstract：BasedonTaylorseriesexpansionandstraincomponentexpressionsin elasticmechanics，wederivetheformulaofstrain androtation tensorforsmal1 arraysinspherica1coordinatesystem．By1inearization processoftheformula， wealsoderivetheexpressionsofstraincomponentsandEulervectoruncertain— ty，respectively forsubnetsusing thelaw ofuncertaintypropagation．Taking GPSvelocityfieldin Sichuan-Yunnan region asan example，wecompute