GPS监测网位移和速率模型.docVIP

GPS监测网位移和速率模型 陶本藻(武汉测绘科技大学　430079) 　　【摘　要】　本文简述了已有的位移和速率模型，然后结合GPS监测网的特点，探讨了建立位移和速率模型问题，并给出了相应的检验位移和速率显著性的统计量及若干速率的拟合模型。http;//测绘信息网　　由于GPS测量的高精度、高速度，不受测边的长度和通视条件的限制等优点，国内外已广泛使用GPS技术来监测地壳运动和工程建筑物的变形。常规方法布设的形变监测网，考虑测量手段的精度和通视条件，监测网范围局限于小区域和短边长，而且平面位移与高程变化分开监测。GPS监测网是测定三维形变，而且所布的网点可视需要而定，网点的稳定性大致可先验预料。根据国内外文献报导［4］，采用一定方法，GPS测定形变，局部网的水平精度可达到±1～±2 mm，高程精度达到±5 mm左右。　　有关监测网的位移和速率模型，近一二十年来国内外已有许多研究，本文结合GPS监测网的特点，运用已有的成果，讨论其位移和速率模型的建立，以及位移速率显著性的检验方法等。　　　　一、位移速率模型　　设两期观测模型为： E(l1)=A1x1　　　　D(l1)＝σ20Q1＝σ20p-11　　　E(l2)=A2x2　　　　D(l2)＝σ20Q2＝σ20p-12　　(1) 分期平差模型为：http;//测绘信息网 即在基准权为px情况下作出xi的估计i。其平差结果为： 　　在两期时间间隔内的位移定义为： d=x2-x1 其估值为。位移的协因数为： 　　由式(4)知，位移d的基准采用了与确定i的基准相同，亦即位移量大小是在所选的基准权px下定义的。在自由网平差理论中已经知道，不论是固定基准，重心基准还是拟稳基准均可由px的特定形式表达［1］。　　两期观测位移模型也可写成下式： 　　(5) 相应的速率模型为：http;//测绘信息网 平差时x1与d(λ)一般是在同一基准下求得，如果选择不同基准，就可选取各自不同的px1和pd(pλ)，其相应的附加基准条件分别为： 此时求出的,其基准属于式(8)所确定。　　位移或速率显著性检验，可采用线性假设法［2］，原假设为： 　　 检验统计量为：http;//测绘信息网 式中，fd(fλ)为d(λ)中函数独立()的个数，单位权方差估计为： 　　(11) fi为i期多余观测数。http;//测绘信息网　　如果检验其中单点位移或速率，即H0∶dj=0或λj=0，则检验统计量为： 　　对于速率模型式(6)，很容易扩大至多期观测，计算单位时间形变量并进行其显著性检验。　　　　二、GPS网位移速率模型　　在三维直角坐标系，基线向量的误差方程为： 　　(13) 式中，xi为点i的三维坐标，lij为基线向量观测值。　　设GPS网有k期观测，则误差方程为： 为了检验k期观测两两位移总体显著性，需建立如下假设： H0∶ Hx=0 即线性假设条件为：http;//测绘信息网 按线性假设法，构成残差平方和： ΩH=Ω+R　　(15) 式中， 因为 式中， 由检验公式知，式(14)可单独分期平差。　　若要对其中任两期位移进行显著性检验，就可利用检验统计量(式9)。　　单点位移检验，考虑GPS点的三维坐标是相关的，一般最好以点的三维坐标的位移显著性一次性检验。http;//测绘信息网　　如果GPS网是变形监测的参考网，必须保证网点在监测期内是稳定的，上述检验如果原假设H0成立，表示参考网稳定。　　顺便指出，GPS网平差必须考虑网点坐标及其位移的基准。　　对于速率模型，基线向量的误差方程为： 　　(21) 式中，Δk=tk-t0，x0i为t0时刻i点的三维坐标向量。式(21)所描述的是，假定在t0至tk时刻内点的运动速率均匀，即单位时间位移相同。　　如果要考虑每期GPS监测网存在各自的尺度和方位的差异，为消除这种系统性影响，可在每期GPS网中分别设立一个尺度参数和3个方位参数，令其为：u=1+δu,αx,αy,αz，则误差方程为： 式中，Δx,Δy,Δz是由近似坐标算得的坐标差。　　全网的误差方程为：http;//测绘信息网 vk=Akxk+Bkuk-lk　　(23) 式中，xk为第k期坐标向量，uk为第k期尺度、方位参数向量，即 uk=［δuk　αx　αy　αz］T 文献［3］对分期平差基准问题作了详细讨论，提出固定位置、重心尺度、重心方位基准，固定位置和尺度、重心方位的基准，拟稳基准、重心基准等4种不同基准，给出了具体的附加基准条件。文中还对参数向量u的检验作了讨论，给出了检验公式。　　设有两期观测，即式(23)中k=1,2，选择基准进行平差，得和,求得位移量：http;//测绘信息网 　　检验H0∶ d=0，如果H0成立，可视为随机误差向量，由此可由计算单位权方差：