-
摘要: 基于1999年以来的GPS观测数据，重点分析了川滇地区近期的变形特征，计算给出了该区域主要断裂带的滑动速率，并结合GPS剖面初步识别了断裂带所处的孕震阶段。结果表明，红河断裂带始终处于弱形变区域，值得进一步关注；汶川地震后，由于研究区域应力场的调整和应变的重新分配，使得安宁河断裂带的活动变化较小，逐步趋缓，走滑分量存在明显趋缓现象，并且安宁河断裂带的闭锁深度在一定程度上有所增加；小江断裂带的走滑分量存在明显趋缓现象，显示存在强闭锁现象，小江断裂带平行于断层方向GPS速度剖面显示断裂带两侧变形平缓且变形范围较大；昭通、莲峰断裂的挤压和剪切变形积累均存在弱化现象，可能预示该区域挤压应变积累背景较高。Abstract: Based on the GPS observation data since 1999, recent deformation characteristics of Sichuan-yunnan region is analyzed, slip rates of main faults are calculated, meanwhile, GPS profiles are combined to analyze the periods of the earthquake generating. The results show that, East edge of Sichuan-Yunnan block:Xianshuihe fault, Anninghe fault, Zemuhe fault, Xiaojiang fault have intense activity, and the movement of the faults are given priority to left-lateral strike-slip. Honghe fault is still in the weak deformation area and the middle section of the Honghe fault zone is the maximum shear strain rate region, while the shear deformation in the northwest and southeast segment of the Honghe fault is higher, which deserves further attention. After Wenchuan earthquake, due to the adjustment of regional stress field and strain redistribution, the activity of the Anning river fault zone is relatively small, easing gradually, and the sliding component has the obvious tendency of slowing down, and the blocking depth of the Anning river fault zone has increased to a certain extent. With the south of the material, the speed of the south segment of the Xiaojiang fault is slowed down and the locking depth increases. GPS velocity profile in parallel to the fault direction of the Xiaojiang fault line shows that the deformation of the fault zone is smooth and the deformation scale is large, showing that there are strong blocking phenomenon. After the Ludian earthquake, the strain release of Zhaotong and Lianfeng fault zone was not significant, compression strain and shear deformation accumulation of Zhaotong-lianfeng is weakening, may indicate that the compressive strain accumulating background is higher.
-
Keywords:
- Sichuan-Yunnan region /
- fault /
- velocity field /
- strain rate /
- slip rate /
- GPS velocity profile
-
数据融合最初仅仅局限于军事领域，是针对异质传感器产生的海量数据进行一致性估计的方法[1-3]。20世纪90年代后，数据融合的概念被引入中国[4]。文献[5]将数据融合的结构分为3层，即数据层、特征层和决策层，在不同层进行数据融合会得到不同的结果，其差异体现在融合精度和处理效率等方面，此结构也作为基本的数据融合结构应用于其他领域。
在大地测量领域，文献[6]从函数模型误差和随机模型误差两方面详细讨论了不同的融合结果，具有较大的理论意义，相关的数据融合应用也较为广泛，体现在GPS与INS组合导航[7]、GPS与InSAR[8]、GPS不同系统间的融合[9]、遥感图像及地图数据融合[9, 10]以及3S集成[11]等方面。在工业测量中，相关的数据融合理论方法并不常见，文献[12]介绍了在盾构自动导向中综合利用全站仪和倾斜仪数据联合平差解算盾构姿态的方法，体现了采用多源异质传感器进行数据融合的优势。
在mm级精度的工业测量中，是采用全站仪进行特定点的三维坐标的采集，全站仪是测距仪(获取斜距)和经纬仪(获取角度)的集成，其获取的坐标实际是测距数据和测角数据的融合结果；而在船舶、飞机制造等工业测量领域，由于一些大型构件无法一站式获取全部待测点位三维信息，需采用移站测量(moving station method, MSM)的方法[13]，通过在不同测站上重复测量公共点，将不同测站的坐标系统一起来。但是，随着移站次数的增加，测量误差会不断积累，影响整体测量的精度。
本文以工业测量整体平差(overall adjustment, OA)为例，讨论了两种不同数据融合方法，即基于原始观测值(水平角、竖直角、斜距)的平差模型和基于坐标观测值(直角坐标系中的三维坐标)的平差模型，从解算精度、多余观测值的敏感性以及随机模型偏差对结果的影响等3方面，详细对比了两者的区别并和传统的移站测量结果进行比较，得到在工业测量中融合策略选取的相关结论。
1 工业测量整体平差函数模型
1.1 基于原始观测值的模型
在传统的工业测量中，采用全站仪获取数据是使用最多的方法之一。在每站测量前经过整平后，全站仪的原始观测值为测站点到待测点的水平方向、竖直角和斜距，这3个观测值互相独立，并在全站仪出厂时会给出角度和距离测量的标称精度。通常所有的待测点无法在一个测站全部观测，故采用移站测量的方法将不同测站的坐标系通过重复测量公共点进行统一，一般都以起始测站的坐标系作为全局坐标系。
在基于原始观测值的整体平差函数模型中，设第j测站在全局坐标系中的坐标为$ \left( {{{\hat{X}}}_{\text{stat}-j}}, \ \ {{{\hat{Y}}}_{\text{stat}-j}}, \ \ {{{\hat{Z}}}_{\text{stat}-j}} \right) $，测量了n个全局坐标系中的待测点$ \left( {{{\hat{X}}}_{i}}, \ \ {{{\hat{Y}}}_{i}}, \ \ {{{\hat{Z}}}_{i}} \right) $，其中i=1, 2, …, n，则有：
$$ \left\{ \begin{array}{l} {L_{\rm{s}}} + {V_{\rm{s}}} = \sqrt {{{\left( {{{\hat X}_i} - {{\hat X}_{{\rm{stat}} - j}}} \right)}^2} + {{\left( {{{\hat Y}_i} - {{\hat Y}_{{\rm{stat}} - j}}} \right)}^2} + {{\left( {{{\hat Z}_i} - {{\hat Z}_{{\rm{stat}} - j}}} \right)}^2}} \\ {L_{\rm{H}}} + {V_{\rm{H}}} = \arctan \frac{{{{\hat Y}_i} - {{\hat Y}_{{\rm{stat}} - j}}}}{{{{\hat X}_i} - {{\hat X}_{{\rm{stat}} - j}}}} - {{\hat A}_{{\rm{azi}} - j}}\\ {L_{\rm{V}}} + {V_{\rm{V}}} = \arctan \frac{{{{\hat Z}_i} - {{\hat Z}_{{\rm{stat}} - j}}}}{{\sqrt {{{\left( {{{\hat X}_i} - {{\hat X}_{{\rm{stat}} - j}}} \right)}^2} + {{\left( {{{\hat Y}_i} - {{\hat Y}_{{\rm{stat}} - j}}} \right)}^2}} }} \end{array} \right. $$ (1) 式中，Ls、LH、LV分别为斜距、水平角、竖直角观测值；对应Vs、VH、VV为观测值的改正数；$ {{\hat{A}}_{\text{azi}-j}} $为第j测站在全局坐标系中水平角初始零方向参数。将式(1)线性化后有：
$$ \left\{ \begin{array}{l} {V_{\rm{s}}} = {p_1}\delta {{\hat X}_{{\rm{stat}} - j}} + {p_2}\delta {{\hat Y}_{{\rm{stat}} - j}} + {p_3}\delta {{\hat Z}_{{\rm{stat}} - j}} - {p_1}\delta {{\hat X}_i} + {p_2}\delta {{\hat Y}_i} + {p_3}\delta {{\hat Z}_i} - {l_s}\\ {V_H} = - {p_4}\delta {{\hat X}_{{\rm{stat}} - j}} + {p_5}\delta {{\hat Y}_{{\rm{stat}} - j}} + {p_4}\delta {{\hat X}_i} - {p_5}\delta {{\hat Y}_i} - \delta {{\hat A}_{{\rm{az}} - j}} - {l_H}\\ {V_V} = {p_6}\delta {{\hat X}_{{\rm{stat}} - j}} + {p_7}\delta {{\hat Y}_{{\rm{stat}} - j}} + {p_8}\delta {{\hat Z}_{{\rm{stat}} - j}} - {p_6}\delta {{\hat X}_i} + {p_7}\delta {{\hat Y}_i} + {p_8}\delta {{\hat Z}_i} - {l_V} \end{array} \right. $$ (2) $$ \left\{ \begin{array}{l} {p_1} = \frac{{{X_{{\rm{stat}} - {j_0}}} + {X_{{i_0}}}}}{{{S_{{i_0}}}}}\\ {p_2} = \frac{{{Y_{{\rm{stat}} - {j_0}}} + {Y_{{i_0}}}}}{{{S_{{i_0}}}}}\\ {p_3} = \frac{{{Z_{{\rm{stat}} - {j_0}}} + {Z_{{i_0}}}}}{{{S_{{i_0}}}}}\\ {p_4} = \frac{{{Y_{{\rm{stat}} - {j_0}}} + {Y_{{i_0}}}}}{{D_{{i_0}}^2}}\\ {p_5} = \frac{{{X_{{\rm{stat}} - {j_0}}} + {X_{{i_0}}}}}{{D_{{i_0}}^2}}\\ {p_6} = \frac{{\left( {{Z_{{\rm{stat}} - {j_0}}} - {Z_{{i_0}}}} \right) \cdot \left( {{X_{{\rm{stat}} - {j_0}}} - {X_{{i_0}}}} \right)}}{{{D_{{i_0}}} \cdot \left( {D_{{i_0}}^2 + {{\left( {{Z_{{\rm{stat}} - {j_0}}} - {Z_{{i_0}}}} \right)}^2}} \right)}}\\ {p_7} = = \frac{{\left( {{Z_{{\rm{stat}} - {j_0}}} - {Z_{{i_0}}}} \right) \cdot \left( {{Y_{{\rm{stat}} - {j_0}}} - {Y_{{i_0}}}} \right)}}{{{D_{{i_0}}} \cdot \left( {D_{{i_0}}^2 + {{\left( {{Z_{{\rm{stat}} - {j_0}}} - {Z_{{i_0}}}} \right)}^2}} \right)}}\\ {p_8} = \frac{{ - D_{{i_0}}^2}}{{{D_{{i_0}}} \cdot \left( {D_{{i_0}}^2 + {{\left( {{Z_{{\rm{stat}} - {j_0}}} - {Z_{{i_0}}}} \right)}^2}} \right)}}\\ {l_{\rm{s}}} = {L_{\rm{s}}} - {S_{{i_0}}};\\ {l_H} = {L_H} + {A_0} - {A_1}\\ {l_V} = {L_V} - \arctan \frac{{{Z_i} - {Z_{{\rm{stat}} - {j_0}}}}}{{{D_{{i_0}}}}} \end{array} \right. $$ (3) 式中，Xstat-j0、Ystat-j0、Zstat-j0为在全局坐标系中第j测站坐标的近似值；Xi0、Yi0、Zi0为在全局坐标系中第i目标点的坐标近似值；A0为初始零方向近似值；A1为通过测站和目标点近似值解算的初始零方向，两者的范围都为0°~360°；Si0和Di0分别为通过测站点和目标点近似值求得的点间距及其水平投影。由于式(1)为非线性函数，求解通过最小二乘法迭代，这些近似值在不断迭代的过程中经求解的改正数调整，直到求解的改正数小于某一给定的阈值δ。
在某次测量任务中，假设共测了s个测站，目标点一共有t个，以第一站的坐标系作为全局坐标系，则需要解得参数共有(s－1)×3+t×3+(s－1)个。需要注意的是，在t个目标点中某些点是可以在多个测站中重复观测的，也就是在移站测量中的公共点。参数解为：
$$ \begin{array}{*{20}{c}} {\mathit{\boldsymbol{\hat X}} = }\\ {\left( {\underbrace {\begin{array}{*{20}{c}} {\delta {{\hat X}_{{\rm{stat}} - 2}}}&{\delta {{\hat Y}_{{\rm{stat}} - 2}}}&{\delta {{\hat Z}_{{\rm{stat}} - 2}}}& \cdots &{\delta {{\hat X}_{{\rm{stat}} - s}}}&{\delta {{\hat Y}_{{\rm{stat}} - s}}}&{\delta {{\hat Z}_{{\rm{stat}} - s}}} \end{array}}_{s - 1个测站坐标}} \right.}\\ {\underbrace {\begin{array}{*{20}{c}} {\delta {{\hat X}_1}}&{\delta {{\hat Y}_1}}&{\delta {{\hat Z}_1}}& \cdots &{\delta {{\hat X}_t}}&{\delta {{\hat Y}_t}}&{\delta {{\hat Z}_t}} \end{array}}_{t\;个目标点坐标}}\\ {{{\left. {\underbrace {\begin{array}{*{20}{c}} {\delta {{\hat a}_{{\rm{azi}} - 2}}}& \cdots &{\delta {{\hat a}_{{\rm{azi}} - s}}} \end{array}}_{s - 1个测站初始零方向参数}} \right)}^{\rm{T}}}} \end{array} $$ (4) 由于第一个测站的测站坐标和初始零方向都为0，因此该测站列立的误差方程仅有目标点的坐标参数。
观测方程的随机模型由仪器标定精度给定，测距精度σs=(a+b×Ls), a为固定误差，b为比例误差, Ls为测站点到目标点的斜距；由于一般测量目标点的角度值都是半测回测一次的结果，而标称精度一般为一测回方向的精度，所以实际的测角精度σangle为给定测角精度的$ \sqrt{2} $倍，在空间上带来的误差影响为σangle×Ls。此外，距离观测与角度观测之间是独立的，因此基于原始观测值的整体平差权阵为对角阵。对于一次观测目标点的3个原始观测值权阵为：
$$ {\mathit{\boldsymbol{P}}_1} = {\rm{diag}}\left[ {\begin{array}{*{20}{c}} {\frac{{\sigma _0^2}}{{\sigma _s^2}}}&{\frac{{\sigma _0^2}}{{{{\left( {{\sigma _{{\rm{angle}}}} \times {L_s}} \right)}^2}}}}&{\frac{{\sigma _0^2}}{{{{\left( {{\sigma _{{\rm{angle}}}} \times {L_s}} \right)}^2}}}} \end{array}} \right] $$ (5) 式中，σ0为给定的验前单位权中误差，一般设定为1 mm。验后单位权中误差由式(6)求得：
$$ \hat \sigma _0^2 = \frac{{{\mathit{\boldsymbol{V}}^{\rm{T}}}{\mathit{\boldsymbol{P}}_1}\mathit{\boldsymbol{V}}}}{{n - \left( {s - 1} \right) \times 4 - t \times 3}} $$ (6) 式中，n为总的观测值个数; V为式(2)中的残差向量，则求解参数$ \mathit{\boldsymbol{\hat{X}}} $的精度：
$$ {\mathit{\boldsymbol{D}}_{\hat X\hat X}} = \hat \sigma _0^2{\mathit{\boldsymbol{Q}}_{\hat X\hat X}} = \hat \sigma _0^2{\left( {{\mathit{\boldsymbol{B}}^{\rm{T}}}{\mathit{\boldsymbol{P}}_1}\mathit{\boldsymbol{B}}} \right)^{ - 1}} $$ (7) 式中，矩阵B为式(2)中的系数阵。
1.2 基于坐标观测值的模型
基于坐标观测值的整体平差中，在第j测站，观测i号点的函数模型为：
$$ \left( {\begin{array}{*{20}{c}} {{L_X}}\\ {{L_Y}}\\ {{L_Z}} \end{array}} \right) + \left( {\begin{array}{*{20}{c}} {{V_X}}\\ {{V_Y}}\\ {{V_Z}} \end{array}} \right) = \left( {\begin{array}{*{20}{c}} {\Delta {{\hat X}_j}}\\ {\Delta {{\hat Y}_j}}\\ {\Delta {{\hat Z}_j}} \end{array}} \right) + \left( {\begin{array}{*{20}{c}} {\cos {{\hat \gamma }_j}}&{ - \sin {{\hat \gamma }_j}}&0\\ {\sin {{\hat \gamma }_j}}&{\cos {{\hat \gamma }_j}}&0\\ 0&0&1 \end{array}} \right)\left( {\begin{array}{*{20}{c}} {{{\hat X}_i}}\\ {{{\hat Y}_i}}\\ {{{\hat Z}_i}} \end{array}} \right) $$ (8) 式中，$ \Delta {{\hat{X}}_{j}} $、$ \Delta {{\hat{Y}}_{j}} $、$ \Delta {{\hat{Z}}_{j}} $、$ {\hat{\gamma }}_{j} $为第j测站相对于全局坐标系的平移参数和绕Z轴的旋转角；$ {{\hat{X}}_{i}} $、$ {{\hat{Y}}_{i}} $、$ {{\hat{Z}}_{i}} $为i号待测点在全局坐标系中的三维坐标值。
式(8)中，求解站与站之间转换参数采用的观测值为坐标，而坐标值是由笛卡尔坐标系中原始观测的斜距、水平角和竖直角推算求得的：
$$ \left\{ \begin{array}{l} {L_X} = {L_{\rm{s}}} \cdot \cos {L_H} \cdot \sin {L_V}\\ {L_Y} = {L_{\rm{s}}} \cdot \sin {L_H} \cdot \sin {L_V}\\ {L_Z} = {L_{\rm{s}}} \cdot \cos {L_V} \end{array} \right. $$ (9) 可以看出，虽然原始观测值之间是独立的，但是坐标观测值是相关的，由误差传播可知，坐标观测值的方差阵为：
$$ \begin{array}{*{20}{c}} {{\mathit{\boldsymbol{D}}_{XYZ}} = \mathit{\boldsymbol{A}} \cdot {\rm{diag}}}\\ {\left[ {\begin{array}{*{20}{c}} {\sigma _s^2}&{{{\left( {{\sigma _{{\rm{angle}}}} \times {L_s}} \right)}^2}}&{{{\left( {{\sigma _{{\rm{angle}}}} \times {L_s}} \right)}^2}} \end{array}} \right] \cdot {\mathit{\boldsymbol{A}}^{\rm{T}}}} \end{array} $$ (10) 式中，
$$ \mathit{\boldsymbol{A}} = \left[ {\begin{array}{*{20}{c}} {\cos {L_H}\sin {L_V}}&{\sin {L_H}\sin {L_V}}&{\cos {L_H}\cos {L_V}}\\ {\sin {L_H}\sin {L_V}}&{\cos {L_H}\sin {L_V}}&{\sin {L_H}\cos {L_V}}\\ {\cos {L_V}}&0&{\sin {L_V}} \end{array}} \right] $$ (11) 仍设定验前单位权中误差为σ0，则对应权阵为：
$$ {\mathit{\boldsymbol{P}}_2} = {\left( {\frac{{{\mathit{\boldsymbol{D}}_{XYZ}}}}{{\sigma _0^2}}} \right)^{ - 1}} $$ (12) 与式(1)一样，式(8)也是非线性的函数，需进行线性化，可以参考大角度七参数三维基准转换模型[14]。需要注意的是，由于以第一测站的坐标系作为全局坐标系，其误差方程更为简单：
$$ \left( {\begin{array}{*{20}{c}} {{L_X}}\\ {{L_Y}}\\ {{L_Z}} \end{array}} \right) + \left( {\begin{array}{*{20}{c}} {{V_X}}\\ {{V_Y}}\\ {{V_Z}} \end{array}} \right) = \left( {\begin{array}{*{20}{c}} {{{\hat X}_i}}\\ {{{\hat Y}_i}}\\ {{{\hat Z}_i}} \end{array}} \right) $$ (13) 传统移站测量方法中，站与站之间通过观测公共点将坐标系统一并传递下去。实际公共点在两套坐标系的观测值中都含有误差，而实际上一般仅采用最小二乘方法处理，忽略了系数阵的误差，因此随着移站的增多，误差累积越来越大。
对比§1.1和§1.2可以看出，两者采用不同的函数模型，但最终都能通过整体平差求解出待测点在第一测站坐标系下的坐标，而且没有传统移站测量中出现的误差累积的问题；此外，两者的系数阵中都不存在含观测误差的观测值，因此通过传统最小二乘算法和准确的随机模型即可获得准确的参数估计，而避免采用整体最小二乘等较复杂的方法来提高解算精度。
2 算例与分析
为了对比研究两种数据融合模型的优劣，笔者利用索佳Net05高精度全站仪，围绕某旗杆墩架设了4个测站，并在墩上的4个面观测了24个贴有反射片的待测点，在每个测站上都可以观测到两个墩面共12个目标点，测站1观测1~12号点，测站2观测7~18号点，测站3观测13~24号点，测站4观测19~24号点和1~6号点，如图 1所示。
对于索佳Net05，当反射物为反射片时，测距标称精度为σs = 0.5+1× ppm mm，测角精度为0.5″，则半测回某方向的精度为σangle=0.5×$ \sqrt{2} $≈0.7″，采用§1.1中的融合方法及式(7)的精度评定公式，以第一测站坐标系作为全局坐标系，得到该坐标系下除测站1的其他3个测站的测站坐标及初始零方向参数和24个待测点的坐标，解算结果见表 1、表 2。同理可以得到按照§1.2的融合方法的解算结果，见表 3、表 4。
表 1 基于原始观测值的整体平差结果Table 1. Result of Overall Adjustment Based on Original Data类型 编号 X/m Y/m Z/m $ \hat{\sigma } $X/mm $ \hat{\sigma } $Y/mm $ \hat{\sigma } $Z/mm 待测点 1 -1.847 0 -11.186 5 0.764 1 0.2 0.2 0.1 2 -1.483 9 -10.830 7 0.779 2 0.2 0.1 0.1 3 -1.064 8 -10.431 0 0.796 8 0.1 0.2 0.1 4 -1.872 2 -11.153 5 0.230 3 0.2 0.2 0.1 5 -1.521 4 -10.806 9 0.196 2 0.2 0.1 0.1 6 -1.098 7 -10.396 3 0.173 7 0.1 0.2 0.1 7 -0.786 9 -10.499 8 0.823 1 0.1 0.2 0.1 8 -0.454 7 -10.842 2 0.828 7 0.2 0.2 0.2 9 -0.108 5 -11.190 8 0.835 5 0.2 0.2 0.2 10 -0.768 9 -10.449 1 0.260 0 0.1 0.2 0.1 11 -0.417 1 -10.811 3 0.272 4 0.2 0.2 0.1 12 -0.082 2 -11.157 0 0.318 2 0.2 0.2 0.2 13 -0.093 4 -11.504 2 0.759 5 0.2 0.2 0.2 14 -0.439 7 -11.837 9 0.779 4 0.3 0.2 0.2 15 -0.826 0 -12.211 0 0.767 9 0.3 0.2 0.2 16 -0.054 3 -11.531 9 0.179 0 0.2 0.2 0.2 17 -0.442 7 -11.904 3 0.191 9 0.3 0.2 0.2 18 -0.798 3 -12.253 4 0.168 8 0.3 0.2 0.2 19 -1.116 3 -12.211 5 0.733 2 0.3 0.2 0.1 20 -1.468 5 -11.854 8 0.745 2 0.3 0.2 0.1 21 -1.845 4 -11.474 2 0.735 1 0.2 0.2 0.1 22 -1.150 9 -12.234 0 0.235 4 0.3 0.2 0.1 23 -1.493 8 -11.888 1 0.222 7 0.3 0.2 0.1 24 -1.877 1 -11.502 8 0.192 6 0.2 0.2 0.1 2 11.670 7 -12.423 4 -0.018 8 0.3 2.5 0.1 测站 3 -1.140 5 -23.007 8 -0.013 3 2.1 0.1 0.1 4 -11.095 8 -11.317 5 0.039 3 0.1 1.8 0.1 表 2 测站其他参数平差结果Table 2. Result of Other Parament编号 初始零方向
/rad解算精度/(″) 验后单位权中误差/mm 2 -2.922 588 41.6 0.5 3 0.741 673 37.4 4 5.704 257 37.9 表 3 基于坐标观测值的整体平差结果Table 3. Result of Overall Adjustment Based on Coordinate Data类型 编号 X/m Y/m Z/m $ \hat{\sigma } $X/mm $ \hat{\sigma } $Y/mm $ \hat{\sigma } $Z/mm 待测点 1 -1.847 1 -11.186 5 0.764 1 0.1 0.4 0.0 2 -1.483 8 -10.831 0 0.779 3 0.1 0.3 0.1 3 -1.064 9 -10.431 4 0.796 8 0.1 0.3 0.1 4 -1.872 2 -11.153 6 0.230 3 0.1 0.4 0.0 5 -1.521 4 -10.807 0 0.196 2 0.1 0.3 0.0 6 -1.098 7 -10.396 7 0.173 7 0.1 0.3 0.0 7 -0.786 8 -10.499 2 0.823 0 0.1 0.3 0.1 8 -0.455 0 -10.841 7 0.828 7 0.1 0.3 0.1 9 -0.108 4 -11.190 4 0.835 5 0.1 0.4 0.1 10 -0.768 7 -10.448 7 0.260 0 0.1 0.3 0.1 11 -0.417 4 -10.811 0 0.272 4 0.1 0.3 0.1 12 -0.082 0 -11.156 6 0.318 2 0.1 0.4 0.1 13 -0.093 3 -11.504 0 0.759 5 0.5 0.4 0.1 14 -0.439 8 -11.837 5 0.779 5 0.5 0.3 0.1 15 -0.825 9 -12.210 5 0.767 9 0.6 0.3 0.1 16 -0.054 5 -11.531 7 0.179 1 0.5 0.4 0.1 17 -0.442 8 -11.903 8 0.191 9 0.5 0.3 0.1 18 -0.798 1 -12.252 8 0.168 8 0.6 0.3 0.1 19 -1.116 1 -12.211 9 0.733 2 0.6 0.3 0.1 20 -1.468 5 -11.855 1 0.745 2 0.5 0.3 0.1 21 -1.845 6 -11.474 2 0.735 1 0.5 0.4 0.1 22 -1.150 8 -12.234 4 0.235 4 0.6 0.3 0.1 23 -1.494 0 -11.888 4 0.222 7 0.5 0.3 0.1 24 -1.877 4 -11.502 9 0.192 6 0.5 0.4 0.1 2 8.697 1 -14.663 0 0.018 8 5.2 1.5 0.0 测站 3 16.385 9 16.191 1 0.013 3 3.7 4.0 0.0 4 3.090 5 15.541 9 -0.039 3 4.4 2.2 0.0 表 4 测站其他参数平差结果Table 4. Result of Other Parament编号 绕Z轴旋转角/rad 解算精度/(″) 验后单位权中误差/mm 2 2.922 953 100.6 3 -0.742 009 98.3 1.7 4 0.579 411 91.0 从表 1、表 2可以看出，待测点的点位解算精度可达0.3 mm左右，且精度分布均匀。从表 3、表 4可以看出，点位解算精度基本在0.5 mm左右，且测站的定位参数解算精度也不如表 1、表 2的结果高。虽然从数学角度看，§1.2的融合方法中，观测值是由§1.1的原始观测值等价推导得到的，随机模型也是从§1.1的误差模型通过误差传播得来的，但是融合结果却不是完全等价的。
按照传统移站测量的思路，采用§1.1的函数模型进行分站平差，即第一测站观测的1~12点都为全局坐标系的坐标；第二测站观测7~12点来解算测站参数，根据解算得到的测站参数和13~18点的原始观测值解算出13~18点的全局坐标系的坐标；同理，第三测站也可求解出19~24点的全局坐标系坐标。此方法称为基于原始观测值的移站测量法。同理采用§1.2的函数模型进行传统移站测量的方法称为基于坐标观测值的移站测量法。
通过不同的方法得到24个待测点的全局坐标，有如下4种方案：①基于原始观测值的整体平差法；②基于坐标观测值的整体平差法；③基于原始观测值的移站测量法；④基于坐标观测值的移站测量法。分别对比表 1、表 2和表 3、表 4，可以看出方案①的参数解算精度较高。因此，以方案①解算的待测点坐标为基准，可以得到方案①与方案②、方案①与方案③以及方案②与方案④的坐标点位偏差，如图 2所示。
可以看出，方案①与方案②的各点点位偏差较均匀，范围在0.1~0.6 mm之间；而在传统的MSM中，无论是基于原始观测值或是坐标观测值，其解算的点位偏差随着搬站的增加而递增。
图 2中MSM采用6个公共点进行搬站，在经过2次搬站后测得的19~24号点与整体平差的点位偏差可达3~4 mm左右，而实际测量中一般不会采用这么多公共点进行移站测量。如若在24个待测点中选取12个点，如选取1、3、5、…、19、21、23号待测点，每两个测站之间只有3个公共点，即减少了多余观测数，可以得到两种OA结果与对应传统MSM的点位偏差，如图 3所示。
![]()
图 3 3个公共点时整体平差与移站测量的点位偏差Figure 3. Coordinate Differences Between OA and MSM with 3 Common Points从图 2和图 3可以看出，基于原始观测值的OA和MSM受多余观测值数量的影响较小，而基于坐标观测值的OA和MSM受多余观测数的影响非常明显。当搬站公共点仅有3个的时候，经两次搬站后测量的目标点与整体平差点位偏差可达将近3 cm，因此基于坐标观测值的MSM的点位偏差更加难以控制。
为了进一步比较当随机模型与实测数据存在偏差，即存在随机模型误差时，方案①和方案②的解算稳定性，在所有斜距原始观测值分别添加中误差为1 mm和5 mm的白噪声(标称精度仅为0.5 mm)，而随机模型仍然采用仪器的标称精度计算，可以得到两种方案在相同额外白噪声影响下解算的坐标偏差，如图 4所示。
![]()
图 4 添加不同白噪声对两种解算点位偏差影响Figure 4. Influence on the Coordinate Deviation when Adding Different White Noise方案①在外加1 mm斜距误差后解算的24个待测点位与不加误差时偏差在0.2~1.5 mm，在外加5 mm斜距误差后解算的24个待测点位与不加误差时偏差在1.1~5.9 mm；同样的条件，方案②的偏差分别在0~0.8 mm和1.1~2.9 mm。从以上偏差的范围和稳定性都可以看出，当随机模型含有误差时，方案②参数解算稳定性更高。
3 结语
1) 当随机模型的误差不明显，即仪器实际测量精度和标称精度相差不大时，采用基于原始观测值的整体平差融合算法来解算待定点坐标的方法精度最高。在传统的移站测量解算中，采用基于原始观测值的方法也比基于坐标观测值的结果要稳定。
2) 当多余观测数较少，即公共点个数较少时，若采用传统移站测量法解算点位坐标，不宜采用基于坐标观测值的融合方法，采用基于原始观测值的融合方法会更可靠。
3) 当随机模型含有明显误差时，基于坐标观测值的整体平差法的融合结果要比基于原始观测值的整体平差融合结果稳定可靠。
致谢: 感谢中国地震局第一监测中心武艳强博士在计算过程中给予的无私帮助和指导。 -
![]()
图 5 安宁河、小江断裂带剖面反正切位错模型分析
Figure 5. The Arctangent Dislocation Model Analyze of Anninghe Fault and Xiaojiang Fault Profiles
表 1 主要断裂带滑动结果/(mm·a-1)
Table 1 Slip Rate of Principle Faults/(mm·a-1)
断裂带 1999~2007年(左旋为正，拉张为正) 2009~2013年(左旋为正，拉张为正) 走滑速率 拉张/挤压速率 走滑速率 拉张/挤压速率 安宁河断裂 4.49±0.32 -0.14±0.37 2.33±0.77 -0.41±0.85 则木河断裂 3.73±1.65 0.74±1.44 2.33±0.58 2.33±0.61 大凉山断裂 2.91±0.54 -2.26±0.55 3.83±0.87 -1.97±0.85 丽江小金河断裂南段 5.96±0.8 -1.11±0.81 3.62±0.86 -1.98±0.86 丽江小金河断裂北段 1.83±0.67 -5.90±0.63 -0.34±0.96 -3.45±0.95 程海断裂 2.41±0.72 4.63±0.67 0.87±2.1 4.2±2.28 红河断裂 -0.38±0.9 1.92±0.85 -0.44±1.1 2.72±1.05 小江断裂北段 3.25±0.1 2.05±0.1 4.14±0.83 2.51±0.69 小江断裂中段 7.16±0.76 0.72±0.67 4.99±0.87 2.84±0.75 小江断裂南段 6.72±0.94 1.68±0.86 4.5±0.66 2.79±0.74 昭通鲁甸断裂 -6.22±0.1 -5.89±0.1 -2.22±0.99 -4.57±0.99 莲峰断裂 -3.85±0.1 -2.58±0.1 -2.9±0.66 -4.58±0.66 
下载: 导出CSV
-
[1] Bayer R, Chery J, Tatar M, et al. Active Deformation in Zagros-Makran Transition Zoneinferred from GPS Measurements[J]. Geophys J Int, 2006, 165:373-381 doi: 10.1111/gji.2006.165.issue-1
[2] Wu Y Q, Jiang Z S, Wang M, et al. Preliminary Results of the Co-seismic Displacement and Pre-seismic Strain Accumulation of the Lushan Ms7.0 Earthquake Reflected by the GPS Surveying[J]. Chin Sci Bull, 2013, 58, doi: 10.1007/s11434-013-5998-5
[3] 乔学军, 王琪, 杜瑞林.川滇地区活动地块现今地壳形变特征[J].地球物理学报, 2004, 47(5):805-811 http://kns.cnki.net/KCMS/detail/detail.aspx?filename=dqwx200405010&dbname=CJFD&dbcode=CJFQ Qiao Xuejun, Wang Qi, Du Ruilin. Characteristics of Current Crustal Deformation of Active Blocks in the Sichuan-Yunnan Region[J].Chinese J.Geophys(in Chinese), 2004, 47(5):805-811 http://kns.cnki.net/KCMS/detail/detail.aspx?filename=dqwx200405010&dbname=CJFD&dbcode=CJFQ
[4] 丁开华, 许才军, 邹蓉, 等.利用GPS分析川滇地区活动地块运动与应变模型[J].武汉大学学报·信息科学版, 2013, 38(7):822-827 http://ch.whu.edu.cn/CN/abstract/abstract2694.shtml Ding Kaihua, Xu Caijun, Zou Rong, et al. Crustal Movement and Strain Model of Active Blocks Analyzed by GPS in Sichuan-Yunnan Region[J]. Geomatics and Information Science of Wuhan University, 2013, 38(7):822-827 http://ch.whu.edu.cn/CN/abstract/abstract2694.shtml
[5] 魏文薪, 江在森, 武艳强, 等.利用GPS数据研究川滇块体东边界主要断裂带运动特性[J].武汉大学学报·信息科学版, 2012, 37(9):1041-1044 http://ch.whu.edu.cn/CN/abstract/abstract315.shtml Wei Wenxin, Jiang Zaisen, Wu Yanqiang, et al. Motion Characteristics of Major Faults in East Boundary of Sichuan-Yunnan Block Obtained with GPS Data[J]. Geomatics and Information Science of Wuhan University, 2012, 37(9):1041-1044 http://ch.whu.edu.cn/CN/abstract/abstract315.shtml
[6] 赵静, 刘杰, 牛安福, 等.大凉山次级块体周边断层的闭锁特征[J].地震地质, 2014, 36(4):1136-1144 http://www.docin.com/p-1202435378.html Zhao Jing, Liu Jie, Niu Anfu, et al. Study on Fault Locking Characteristic Around the Daliangshan Sub-Block[J]. Seismology and Geology, 2014, 36(4):1136-1144 http://www.docin.com/p-1202435378.html
[7] 方颖, 江在森, 邵志刚, 等.利用GPS连续资料分析川滇地区的地壳变形特征[J].地震研究, 2014, 37(2):204-209 doi: 10.3969/j.issn.1000-0666.2014.02.006 Fang Ying, Jiang Zaisen, Shao Zhigang, et al. Analysis on Crustal Deformation Characteristics in Sichuan-Yunnan Region with GPS Continuous Data[J].Journal of Seismological Research, 2014, 37(2):204-209 doi: 10.3969/j.issn.1000-0666.2014.02.006
[8] 张培震, 邓起东, 张竹琪.中国大陆的活动断裂、地震灾害及其动力过程[J].中国科学:地球科学, 2013, 43(10):1607-1620 http://www.wenkuxiazai.com/doc/81d94fc776c66137ee0619f7.html Zhang Peizhen, Deng Qidong, Zhang Zhuqi, et al. Active Faults, Earthquake Hazards and Associated Geodynamic Processes in Continental China[J].Scientia Sinica Terrae, 2013, 43(10):1607-1620 http://www.wenkuxiazai.com/doc/81d94fc776c66137ee0619f7.html
[9] 杨博、张风霜、韩月萍.球面水平应变场无偏差计算与滤波[J].地震研究, 2011, 34(1):59-66 http://kns.cnki.net/KCMS/detail/detail.aspx?filename=dzyj201101011&dbname=CJFD&dbcode=CJFQ Yang Bo, Zhang Fengshuang, Han Yueping. Non-Deviation Calculation and Filtering of Spherical Horizontal Stain Field[J]. Journal of Seismological, 2011, 34(1):59-66 http://kns.cnki.net/KCMS/detail/detail.aspx?filename=dzyj201101011&dbname=CJFD&dbcode=CJFQ
[10] 王阎昭, 王恩宁, 沈正康, 等.基于GPS资料约束反演川滇地区主要断裂现今活动速率[J].中国科学D辑:地球科学, 2008, 38(5):582-597 http://www.oalib.com/paper/4150550 Wang Yanzhao, Wang Enning, Shen Zhengkang, et al. Inversion of Current Movement Velocity of Sichuan-Yunnan Main Rupture Constrained by GPS Data[J].Science in China(Series D:Earth Science), 2008, 38(5):582-597 http://www.oalib.com/paper/4150550
[11] 蒋锋云, 朱良玉, 王双绪, 等.川滇地区地壳块体运动特征研究[J].地震研究, 2013, 36(3):263-268 http://www.cnki.com.cn/Article/CJFDTotal-DZYJ201303001.htm Jiang Fengyun, Zhu Liangyu, Wang Shuangxu, et al. Research on the Characteristics of Crustal Block Movement in Sichuan-Yunnan Area[J].Journal of Seismlolgical Research, 2013, 36(3):263-268 http://www.cnki.com.cn/Article/CJFDTotal-DZYJ201303001.htm
[12] Wu Y Q, Jiang Z S, Yang G H, et al. Comparison of GPS Strain Rate Computing Methods and Their Reliability[J]. Geophys. J. Int., 2011, 185:703-717 doi: 10.1111/gji.2011.185.issue-2
[13] 武艳强, 江在森, 杨国华, 等.利用最小二乘配置在球面上整体解算GPS应变场的方法及应用[J].地球物理学报, 2009, 52(7):1707-1714 http://kns.cnki.net/KCMS/detail/detail.aspx?filename=dqwx200907006&dbname=CJFD&dbcode=CJFQ Wu Yanqiang, Jiang Zaisen, Yang Guohua, et a1.The Application and Method of GPS Strain Calculation in Whole Mode Using Least Square Collocation in Sphere Surface[J].ChineseJ.Geophys, 2009, 52(7):1707-1714 http://kns.cnki.net/KCMS/detail/detail.aspx?filename=dqwx200907006&dbname=CJFD&dbcode=CJFQ
[14] 武艳强. 高性能计算与GPS地壳形变信息提取[D]. 北京: 中国地震局地震预测研究所, 2008 http://cdmd.cnki.com.cn/Article/CDMD-85405-2008129224.htm Wu Yanqiang. High Performance Computing and Distilling Crustal Deformation Information from GPS Data[D]. Beijing:Institute of Earthquake Science, China Earthquake Administration, 2008 http://cdmd.cnki.com.cn/Article/CDMD-85405-2008129224.htm
[15] 武艳强, 江在森, 杨国华, 等.南北地震带北段近期地壳变形特征研究[J].武汉大学学报·信息科学版, 2012, 37(9):1045-1048 http://ch.whu.edu.cn/CN/abstract/abstract316.shtml Wu Yanqiang, Jiang Zaisen, Yang Guohua, et al.Deformation Characteristics of North Section of North-South Seismic Zone in Recent Period[J]. Geomatics and Information Science of Wuhan University, 2012, 37(9):1045-1048 http://ch.whu.edu.cn/CN/abstract/abstract316.shtml
[16] 周荣军, 何玉林, 杨涛, 等.鲜水河-安宁河断裂带磨西-冕宁段的滑动速率与强震位错[J].中国地震, 2001, 17(3):253-262 http://mall.cnki.net/magazine/Article/ZGZD200103003.htm Zhou Rongjun, He Yulin, Yang Tao, et al. Slip Rate and Strong Earthquake Rupture on the Moxi-Mianning Segment Along the Xianshuihe-Anninghe Fault Zone[J].Earthquake Research in China, 2001, 17(3):253-262 http://mall.cnki.net/magazine/Article/ZGZD200103003.htm
[17] 徐锡伟, 闻学泽, 郑荣章, 等.川滇地区活动块体最新构造变动样式及其动力来源[J].中国科学D辑:地球科学, 2003, 33(增刊):210-216 doi: 10.3321/j.issn:1006-9267.2003.z1.017 Xu Xiwei, Wen Xueze, Zheng Rongzhang, et al. Pattern of Latest Tectonic Motion and its Dynamics for Active Blocks in Sichuan-Yunnan Region, China[J]. Science in China(Series D), 2003, 46(Suppl.):210-216 doi: 10.3321/j.issn:1006-9267.2003.z1.017
[18] Shen Z K, Lv J, Wang M, et al. Contemporary Crustal Deformation Around the Southeast Borderland of the Tibetan Plateau[J]. J Geophys Res, 2005, 110:B11409, doi: 10.1029/2004JB003421
[19] 闻学泽, 杜方, 易桂喜, 等.川滇交界东段昭通、莲峰断裂带的地震危险背景[J].地球物理学报, 2013, 56(10):3361-3372 doi: 10.6038/cjg20131012 Wen Xueze, Du Fang, Yi Guixi, et al. Earthquake Potential of Zhaotong and Lianfeng Fault Zones of the Eastern Sichuan-Yunnan Border Region[J]. Chinese J.Geophys, 2013, 56(10):3361-3372 doi: 10.6038/cjg20131012
-
期刊类型引用(4)
1. 朱杰，郑加柱，陈红华，杨静，胡平昌，陆敏燕. 结合POI数据的南京市商业中心识别与集聚特征研究. 现代测绘. 2022(06): 34-39 . 
百度学术
2. 金澄，安晓亚，陈占龙，马啸川. 矢量居民地多边形多级图划分聚类方法. 武汉大学学报(信息科学版). 2021(01): 19-29 . 百度学术
3. 张铭龙，何贞铭. 基于因子分析法的城市商业中心抽取研究. 地理空间信息. 2021(08): 58-60+64+5 . 百度学术
4. 李卫东，张铭龙，段金龙. 基于POI数据的南京市空间格局定量研究. 世界地理研究. 2020(02): 317-326 . 百度学术
其他类型引用(3)
计量
- 文章访问数:
- HTML全文浏览量:
- PDF下载量:
- 被引次数: 7
