更全的杂志信息网
期刊发表 | 网站导航

多站载波相位差分GPS技术在风沙地貌台风响应模式研究中的应用与实现

更新时间：2009-03-28

0 前言

风沙运动是一个不断变化的过程，对其地形变化的观测是其研究的基本手段之一[1]。传统的海岸风沙观测方法通常有两种。第一种是：采用在观测区域内设定若干个标志（如埋石、插杆等）作为观测点进行定点观测，但由于这些点设定在沙质的地面上易发生位移、人为地破坏或被风沙刮掉，往往会导致观测失败，同时采取该设点方式的观测点密度不可能很大，也影响到观测结果的精度；第二种是：使用水准仪测量高程的变化，但该方法无法确定点的坐标。由于海岸风沙地形的变化量比较小，进行测量时的精度必须要达到厘米级[2]，这是传统观测技术难以达到的。而单站差分RTKGPS测量出来的坐标和高程精度只能达到分米级，本文根据风沙运动的特点和单站差分RTKGPS的特点，提出多站差分RTKGPS的测量方法，并通过对海岸风沙的实际观测检验其应用效果。

1 单站载波相位差分GPS技术

差分GPS通常分为伪距差分GPS和载波相位差分GPS[3]。单站载波相位差分GPS是至少用两台GPS接收机，一台作为基准站，另一台作为移动台，同步观测相同的GPS卫星，确定两台接收机天线之间的相对位置差（差分值），从而实现移动台的高精度定位[4]。

单站载波相位差分GPS的基本原理（图1）：将基准站采集的载波相位发给移动台，进行求解计算移动台的坐标、高程，也可以将基准站的载波相位修正值（差分值）发给移动台，改正移动台接收到的载波相位，再解算移动台的坐标、高程[5]。由于受地球表面曲率的影响，基准站和移动台之间是一种非线性关系，而单站载波相位差分GPS计算出来的差分值是一种线性关系，用线性关系的数值改正与基准站为非线性关系的移动台的坐标，其精度只有分米级，无法满足风沙地貌运动厘米级的精度要求。

  

图1 差分GPS测量原理图Fig.1 The principle of difference GPS

2 多站载波相位GPS技术在海岸风沙运动中的观测原理

由于海岸风沙运动导致的地形高程变化量相对很少，一般最少的变化为厘米级[6]，同时GPS测量出来的坐标和高程是在一定的范围内是变化的，要满足厘米级的精度要求，就必须对在多个已知点上测量出来的WGS-84坐标进行平滑处理，并对其高程进行曲面拟合，以实现已知点对应的椭球体曲面和WGS-84坐标对应的椭球体曲面之间的拟合，这是采用多站RTKGPS技术进行海岸风沙运动中地形变化测量的关键。

2.1 坐标平滑处理方法

GPS的坐标系统为WGS-84坐标系统，与我国的坐标系统（北京54或西安-80）不同[7]。GPS采用的椭球长半轴6378137±2 m，扁率1:298.257；我国的坐标系统（以北京54为例）采用的椭球长半轴6 378 245 m，扁率1:298.30[8]。两种坐标系统的原点不同，扁率和长半轴不同，在实际的测量过程中，必须要找出它们之间的转换关系，进行坐标平滑处理，才能实现高精度的测量[9]。具体采用七参数平滑处理的方法：

七参数平滑处理的方法是在至少已知三个点（多站点）的坐标前提下进行的。如图2所示：O1、X1、 Y1、 Z1 和 O2、 X2、 Y2、 Z2 为转换前后的坐标系统，两个坐标系相互间的旋转角为：ΩX、ΩY、ΩZ。

  

图2 三维直角坐标系间的转换Fig.2 The conversion between three-dimensional orthogonal coordinate

由于坐标轴之间的旋转，旋转角为ΩX、ΩY、ΩZ所引起的坐标变换如下：

 

由于GPS的卫星是运动型的卫星，即使在同一个位置利用不同的时间段进行观测得出来的坐标和高程也是不完全相同的，但在长时间的风沙运动观测中必须每次都要对相同的点进行观测，以求出风沙运动所引起的地形高程的变化。要满足这些条件，每次观测前必须要进行坐标的平滑处理和高程曲面拟合。

 

依据上述原理，采用G10A载波相差分GPS测量系统，对如何应用多站载波相位差分GPS技术进行海岸风沙地貌受台风响应中地貌变化进行了野外实际测量，并采用徕卡0.3 mm精密水准仪进行验证。观测地点为福建平潭岛（图4）中的坛南湾白鸽屿后海岸沙丘。

在欠发达省份，农业是主要用水部门，水的利用效率和效益较低。但由于正处在加快发展阶段，水问题也日益复杂严峻。

 

因平潭岛海岸沙丘的典型性和受台风影响的代表性，选择平潭岛白鸽屿后海岸前丘（图4）进行海岸沙丘形态对台风响应的观测与研究。

 

式中：ΔX、ΔY、ΔZ是坐标原点平移引起的坐标变换；称为坐标系1变换到坐标系2的旋转矩阵；ΩX、ΩY、ΩZ是坐标轴旋转的角度；K是尺度因子（一个无量的值），它表明二坐标的长度标准稍有不同。

（2）PmGSTd1在雌、雄虫的各部位均有表达，且均表现为足部的表达量明显高于其他部位；雄虫各部位的表达量（残体除外）均显著高于雌虫，雄虫翅膀、触角、足部、残体中的表达分别是雌虫相应部位的4.49倍、2.11倍、1.68倍、1.33倍。

计算出 ΔX、 ΔY、 ΔZ， ΩX、 ΩY、 ΩZ和 K 等参数后，即可实现不同坐标之间的平滑处理。

XIA Q在他们的研究中简要地解决了医疗数据共享系统中的访问控制管理问题，主要设计了一个基于区块链的数据共享方案[1]，允许数据用户/所有者在身份验证和加密密钥验证后，从共享存储库访问电子病历。SIFAH E B等人也提出了基于区块链的共享医疗数据方案，重点在于提供数据访问控制、出处和审计的同时[2]，在云服务提供商之间共享医疗数据。SHAE Z提出了一个用于临床试验和精密医学的区块链平台架构，并讨论了各种设计方面问题，并对技术要求和挑战提供了一些见解[3]。

只得慰藉于苏珊娜与浮西努相同的身世，神似的气质，假借战争的苦难造成人物情感的偏执走向，并且借弗希里的口：“人在最无助的时候，爱情最容易入侵。在最危险的时候也是，这像嗜酒一样。”

2.2 高程拟合法

纵观人类历史的发展进程和历代王朝的兴亡更替，对于任何一个阶级来说，如果没有代表本阶级利益的领袖，是很难取得统治地位的。而以列宁为首的代表无产阶级利益的领袖在革命的实践活动中则非常重视领袖的权威，并且是符合当时的世情、国情和党情的，因而取得了非凡成就。

  

图3 GPS高程和正常高高程Fig.3 The altitude of GPS and the altitude of level

当旋转角ΩX、ΩY、ΩZ很小时，上式可简化成：

当GPS点布设成一定区域面时，可以应用数学曲面拟合法求待定点的正常高。根据测区中已知点的平面坐标X、Y（或大地坐标B、L）和ζ值（图3），用数值拟合法，拟合出测区似大地水准面，再内插出待求点的ζ，从而求出待求点的正常高。设点的ζ与平面坐标X、Y有以下关系：

 

选择在平潭岛海岸的两个典型海岸沙丘分布区设立野外定位观测点，即位于平潭岛东南部的白鸽屿后 （海岸前丘）典型点进行测量。

 

写成矩阵形式有：

二是水市场不发育，水权交易平台尚未构建。以往的水权交易全部是通过政府作为中介完成的，这不是规范的水市场，也不可能使市场在水资源配置中发挥决定性作用和更好发挥政府作用。水市场电子信息系统有待研发建设。水权交易需要的交易账户注册、登记、信息披露、金额结算、风险控制等功能的软件和硬件系统，以及与水行政主管部门监管实时对接系统等都有待研发建设。内蒙、河南、广东等地正在积极研究组建水权交易中心（所）。

 

其中：

 

对于每个已知点，都可以列出以上方程，在趋向于零的条件下，解出各ai，再按ζ=XB+ε式求出待求点的ζ，从而求出Hr。

3 多站差分RTKGPS技术在风沙运动观测中的应用

可以得到用矩阵形式表示的变换方程：

GPS高程系统采用的是大地高程系统，它是与参考椭球体作为参照系的高程系统，与我国以拟大地水准面作为基准面的高程系统不同。如图3所示：H84是GPS的高程，属于大地高程；Hr是以拟大地水准面作为基准面正常高高程。

3.1 观测区域

把坐标轴的平移量ΔX、ΔY、ΔZ和尺度因子K一起考虑时得出公式（6）：

  

图4 观测区域、观测点位置图及沙丘景观Fig.4 Observation area,observation location map and sand dunes landscape

3.2 技术方法

其中， f（X，Y）为ζ趋势值， ε为误差。

图5中黑色圆点为载波相位差分GPS第一次观测的点位，基准点J0、Y1、Y2为已知点，其坐标值、高程已知。在基准点设立基准站，移动台分别到Y1点和Y2点进行观测，并测量出该两点的WGS-84坐标，通过地面的已知的J0、Y1、Y2点和对应的GPS测量出来的坐标、高程，求解出GPS的WGS-84坐标与已知坐标J0、Y1、Y2的转换7参数及高程异常值参数（表1、表2）。

第二次及后面多次进行测量时，把在J0、Y1、Y2求出的七参数及高程异常值输入到移动台的控制器内，应用GPS的导航功能，每次观测都沿着这些点进行观测。经过不同时段的观测，对各次测量相同点的高程进行相互比较，即可分析出各时期的高程变化，从而实现了观测不同环境、不同时间下的海岸风沙地貌台风响应的地貌变化。

3.3 观测结果

根据已知点坐标（J0、Y1、Y2）和 GPS 在这些已知点上测定的WGS-84坐标和高程，经计算后得出的拟合参数如表1。

从2014年7月14日开始至2017年7月30日，共观测了20次，每次测量相同315个点，表3为第20次（2017年7月30日）观测315个点中的1001号点至1080号点共80个点的坐标数据（Y、X）和高程数据（Hr），利用精密水准仪测量对应点的高程h，利用相同的点的高程比较，即可计算出其误差值。

2004年8月份我国加入了非物质文化遗产保护公约，1998年10月我省设立了第一座生态博物馆——六枝梭嘎生态博物馆，该馆的设立标志着我省非物质文化保护领域进入了新的发展阶段。历经20余年的快速发展，当前我省非物质文化博物馆主要分为三类：（1）区域性的非物质文化博物馆，即以各地市命名的非物质文化博物馆，该类博物馆大部分是国有，由地方财政负责运营养护等其他非要，对经济效益的要求不高；（2）以专题为主要展览内容的博物馆，该类博物馆以专题性文化为主，性质有国营也有民办；（3）传统民俗博物馆，该类博物馆将民族习俗保留下来，将艺人和他们的手艺传承保护下来，供游客参观和鉴赏。

  

图5 白鸽屿后（海岸前丘）观测位置图Fig.5 Observation location map of BaiGeYuHou

 

表1 福建平潭岛已知点坐标Table 1 The known coordinate of Pingtan island,Fujian province

  

坐标 J 0点 Y 1点 Y 2点X 2 8 1 2 5 3 9.0 4 5 2 8 1 3 6 5 7.3 2 8 2 8 1 4 5 2 3.5 0 3 Y 4 7 4 9 6 0.5 1 4 4 7 4 6 3 5 7.1 6 4 7 4 8 4 6 9.1 0 6 H 2 0.8 0 7 1 4.6 3 0 1 3.3 3 7

 

表2 福建平潭岛拟合参数Table 2 The fitting parameters of Pingtan island,Fujian province

  

Ω X Ω Y Ω Z Δ X Δ Y Δ Z K ζ 0 0 0°0 7 9 5 0 0.8 0 4 0 0 0°0 9 9 3 6 0.2 2 7 0 0 0°0 6 9 0 3 0.3 5 1 4 7.3 1 5 6 7 2 7 1.0 2 4 1 6 2 3 7.2 7 5 3 8 6 0.9 9 1 2 5 4 7 5 1 4 4.5 2 7

 

表3 2017年7月30日载波相位差分GPS测量的高程以精密水准仪测量高程比较Table 3 The comparison of carrier phase difference GPS measurement elevation and the leveling elevation on July 30th,2017

  

载波相位差分G P S测量结果观测点序号 Y X H r 1 0 0 1 4 7 4 7 3 2.4 4 2 2 8 1 2 3 6 0.9 7 1 3.4 5 9 1 0 0 2 4 7 4 7 3 0.9 7 6 2 8 1 2 3 6 4.3 6 1 3.4 4 1 1 0 0 3 4 7 4 7 2 9.9 7 6 2 8 1 2 3 6 7.7 8 1 3.4 3 6 1 0 0 4 4 7 4 7 2 9.3 1 8 2 8 1 2 3 6 9.7 1 1 3.4 3 4 1 0 0 5 4 7 4 7 2 8.3 5 1 2 8 1 2 3 7 3.2 0 1 3.4 2 3 1 0 0 6 4 7 4 7 2 7.4 7 1 2 8 1 2 3 7 6.9 6 1 3.4 0 4 1 0 0 7 4 7 4 7 2 6.6 0 3 2 8 1 2 3 7 9.7 1 1 3.3 8 2 1 0 0 8 4 7 4 7 3 9.7 3 8 2 8 1 2 3 8 4.2 2 1 3.5 9 1 1 0 0 9 4 7 4 7 4 0.1 2 4 2 8 1 2 3 8 0.2 9 1 3.5 9 5 1 0 1 0 4 7 4 7 4 0.2 1 1 2 8 1 2 3 7 6.8 4 1 3.5 9 3 1 0 1 1 4 7 4 7 4 0.4 0 1 2 8 1 2 3 7 4.3 9 1 3.6 0 1 1 0 1 2 4 7 4 7 4 0.8 0 0 2 8 1 2 3 7 0.1 0 1 3.5 8 7精密水准仪测量/m h误差/m h-H r 1 3.4 3 1 1 3.4 1 4 1 3.4 2 2 1 3.4 2 1 3.3 8 5 1 3.3 9 3 1 3.3 8 1 1 3.5 9 2 1 3.5 9 4 1 3.5 8 5 1 3.5 8 4 1 3.6-0.0 2 8-0.0 2 7-0.0 1 4-0.0 1 4-0.0 3 8-0.0 1 1-0.0 0 1 0.0 0 1-0.0 0 1-0.0 0 8-0.0 1 7 0.0 1 3 1 0 1 3 1 0 1 4 1 0 1 5 4 7 4 7 4 1.2 7 8 4 7 4 7 4 2.6 3 1 4 7 4 7 5 9.5 4 8 2 8 1 2 3 6 6.3 4 2 8 1 2 3 6 1.1 4 2 8 1 2 3 6 3.0 1 1 3.6 0 3 1 3.6 1 2 1 3.8 6 5 1 3.5 9 1 1 3.6 1 1 3.9 7 8-0.0 1 2-0.0 0 2 0.1 1 3

 

（接表 3）

  

载波相位差分G P S测量结果观测点序号 Y X H r 1 0 1 6 4 7 4 7 5 9.6 1 1 2 8 1 2 3 6 3.5 2 1 3.8 5 6 1 0 1 7 4 7 4 7 5 9.1 6 4 2 8 1 2 3 6 7.3 4 1 3.8 4 0 1 0 1 8 4 7 4 7 5 9.0 6 2 2 8 1 2 3 7 0.2 7 1 3.7 9 9 1 0 1 9 4 7 4 7 5 8.8 5 1 2 8 1 2 3 7 3.8 7 1 3.8 4 9 1 0 2 0 4 7 4 7 5 7.9 3 8 2 8 1 2 3 7 7.6 7 1 3.8 1 6 1 0 2 1 4 7 4 7 5 7.0 1 1 2 8 1 2 3 8 1.9 1 1 3.8 1 6 1 0 2 2 4 7 4 7 5 5.8 0 1 2 8 1 2 3 8 6.6 9 1 3.8 0 5 1 0 2 3 4 7 4 7 6 6.8 3 7 2 8 1 2 3 9 2.1 8 1 4.0 3 2 1 0 2 4 4 7 4 7 6 7.6 8 0 2 8 1 2 3 8 8.2 5 1 4.0 6 2 1 0 2 5 4 7 4 7 6 7.7 8 6 2 8 1 2 3 8 4.0 2 1 4.0 5 0 1 0 2 6 4 7 4 7 6 8.0 5 9 2 8 1 2 3 8 0.9 4 1 4.0 6 2 1 0 2 7 4 7 4 7 6 8.4 2 7 2 8 1 2 3 7 7.0 1 1 4.0 5 8精密水准仪测量/m h误差/m h-H r 1 3.9 6 9 1 3.8 6 2 1 3.8 4 5 1 3.8 5 3 1 3.8 5 2 1 3.8 4 5 1 3.8 0 9 1 4.0 7 1 1 4.0 7 4 1 4.0 6 1 1 4.0 7 2 1 4.0 8 5 0.1 1 3 0.0 2 2 0.0 4 6 0.0 0 4 0.0 3 6 0.0 2 9 0.0 0 4 0.0 3 9 0.0 1 2 0.0 1 1 0.0 1 0 0.0 2 7 1 0 2 8 1 0 2 9 1 0 3 0 1 0 3 1 1 0 3 2 1 0 3 4 1 0 3 5 1 0 3 6 1 0 3 7 1 0 3 8 1 0 3 9 1 0 4 0 1 0 4 1 1 0 4 2 1 0 4 3 1 0 4 4 1 0 4 5 1 0 4 6 1 0 4 7 1 0 4 8 1 0 4 9 1 0 5 0 1 0 5 1 1 0 5 2 1 0 5 3 1 0 5 4 1 0 5 5 1 0 5 6 1 0 5 7 1 0 5 8 4 7 4 7 6 8.8 5 1 4 7 4 7 6 9.3 8 6 4 7 4 7 7 0.6 9 1 4 7 4 7 8 2.3 8 9 4 7 4 7 8 1.8 9 4 4 7 4 7 8 1.1 7 8 4 7 4 7 8 0.0 7 3 4 7 4 7 7 9.6 8 5 4 7 4 7 7 9.0 0 3 4 7 4 7 7 8.2 3 1 4 7 4 7 7 7.7 1 3 4 7 4 7 8 6.3 0 3 4 7 4 7 8 6.8 9 9 4 7 4 7 8 7.3 8 5 4 7 4 7 8 7.8 5 4 4 7 4 7 8 8.5 8 3 4 7 4 7 8 9.1 9 8 4 7 4 7 9 0.0 0 8 4 7 4 7 9 9.2 2 7 4 7 4 7 9 8.4 9 1 4 7 4 7 9 7.7 5 4 7 4 7 9 7.5 8 2 4 7 4 7 9 7.5 8 2 4 7 4 7 9 7.9 8 5 4 7 4 7 9 8.4 7 5 4 7 4 7 9 7.6 5 1 4 7 4 7 9 7.2 7 4 4 7 4 7 9 6.2 8 7 4 7 4 7 9 5.2 9 5 4 7 4 7 9 9.6 7 7 2 8 1 2 3 7 3.2 2 2 8 1 2 3 6 9.4 4 2 8 1 2 3 6 3.9 8 2 8 1 2 3 6 7.4 8 2 8 1 2 3 7 1.2 9 2 8 1 2 3 7 5.2 4 2 8 1 2 3 7 9.5 5 2 8 1 2 3 8 3.4 3 2 8 1 2 3 8 7.1 4 2 8 1 2 3 9 2.2 9 2 8 1 2 3 9 5.7 9 2 8 1 2 4 0 0.6 4 2 8 1 2 3 9 6.4 8 2 8 1 2 3 9 1.4 8 2 8 1 2 3 8 6.9 4 2 8 1 2 3 8 3.5 1 2 8 1 2 3 7 9.6 7 2 8 1 2 3 7 5.2 2 2 8 1 2 3 7 6.8 8 2 8 1 2 3 7 9.3 4 2 8 1 2 3 8 2.1 3 2 8 1 2 3 8 5.1 0 2 8 1 2 3 8 5.1 0 2 8 1 2 3 8 7.5 1 2 8 1 2 3 8 9.7 4 2 8 1 2 3 9 3.0 9 2 8 1 2 3 9 6.2 6 2 8 1 2 3 9 9.6 2 2 8 1 2 4 0 3.4 4 2 8 1 2 4 0 6.4 0 1 4.0 6 5 1 4.0 8 1 1 4.0 9 3 1 4.5 5 7 1 4.4 9 1 1 4.4 5 8 1 4.4 4 0 1 4.4 3 2 1 4.4 0 5 1 4.3 8 7 1 4.3 8 2 1 4.8 5 3 1 4.8 6 6 1 4.8 9 6 1 4.9 3 0 1 5.0 4 4 1 5.0 4 8 1 5.1 4 7 1 6.7 6 1 1 6.7 3 9 1 6.6 0 6 1 6.6 2 0 1 5.5 4 5 1 5.5 4 8 1 5.6 8 9 1 5.4 3 9 1 5.4 5 3 1 5.4 0 7 1 5.3 1 5 1 5.6 7 9 1 4.1 0 1 1 4.1 1 2 1 4.1 1 5 1 4.5 8 1 1 4.5 7 1 1 4.5 0 3 1 4.5 2 0 1 4.5 2 5 1 4.4 1 5 1 4.3 7 9 1 4.3 6 6 1 4.8 4 5 1 4.9 0 1 1 4.8 5 1 1 4.9 2 2 1 5.0 2 3 1 5.0 7 1 1 5.0 7 9 1 6.7 5 4 1 6.7 7 6 1 6.6 1 3 1 6.6 5 8 1 5.5 6 6 1 5.5 7 7 1 5.7 2 2 1 5.4 5 4 1 5.4 2 9 1 5.4 1 5 1 5.3 3 9 1 5.6 5 2 0.0 3 6 0.0 3 1 0.0 2 2 0.0 2 4 0.0 8 0 0.0 4 5 0.0 8 0 0.0 5 3 0.0 1 0-0.0 0 8-0.0 1 6-0.0 0 8 0.0 3 5-0.0 4 5-0.0 0 8-0.0 2 1 0.0 2 3-0.0 6 8-0.0 0 7 0.0 3 7 0.0 0 7 0.0 3 8 0.0 2 1 0.0 2 9 0.0 3 3 0.0 1 5-0.0 2 4 0.0 0 8 0.0 2 4-0.0 2 7 1 0 5 9 1 0 6 0 1 0 6 1 4 7 4 8 0 2.0 6 4 4 7 4 8 0 1.9 9 9 4 7 4 8 0 2.4 4 9 2 8 1 2 4 0 6.2 0 2 8 1 2 4 0 4.3 2 2 8 1 2 4 0 2.2 4 1 5.9 0 5 1 6.5 3 5 1 6.7 5 5 1 5.9 2 1 1 6.5 5 8 1 6.7 8 1 0.0 1 6 0.0 2 3 0.0 2 6

 

（接表 3）

  

载波相位差分G P S测量结果观测点序号 Y X H r 1 0 6 2 4 7 4 8 0 2.2 3 9 2 8 1 2 3 9 9.9 5 1 6.4 2 6 1 0 6 3 4 7 4 8 0 2.2 5 0 2 8 1 2 3 9 8.0 8 1 6.4 5 6 1 0 6 4 4 7 4 8 0 2.6 2 5 2 8 1 2 3 9 5.3 7 1 6.4 3 5 1 0 6 5 4 7 4 8 0 2.3 1 2 2 8 1 2 3 9 3.9 1 1 6.1 7 1 0 6 6 4 7 4 8 0 0.8 9 0 2 8 1 2 3 9 2.2 1 1 6.3 4 2 1 0 6 7 4 7 4 8 0 0.1 0 2 2 8 1 2 3 8 9.6 9 1 6.7 5 6 1 0 6 8 4 7 4 8 0 1.0 0 1 2 8 1 2 3 8 8.6 8 1 6.5 2 7 1 0 6 9 4 7 4 8 0 0.8 9 0 2 8 1 2 3 8 7.1 6 1 6.7 3 8 1 0 7 0 4 7 4 8 0 0.6 4 6 2 8 1 2 3 8 5.5 0 1 7.4 6 8 1 0 7 1 4 7 4 8 0 0.5 2 3 2 8 1 2 3 8 3.6 3 1 6.8 1 5 1 0 7 2 4 7 4 8 0 1.1 1 0 2 8 1 2 3 8 1.2 8 1 6.7 3 1 1 0 7 3 4 7 4 8 0 1.6 0 7 2 8 1 2 3 7 9.8 4 1 6.6 2 6精密水准仪测量/m h误差/m h-H r 1 6.4 0 3 1 6.4 6 1 1 6.4 0 1 1 6.1 7 6 1 6.3 2 7 1 6.7 1 7 1 6.5 7 2 1 6.7 8 5 1 7.4 1 9 1 6.8 0 1 1 6.7 1 2 1 6.6 4 1-0.0 2 3 0.0 0 5-0.0 3 4 0.0 0 6-0.0 1 5-0.0 3 9 0.0 4 5 0.0 4 7-0.0 4 9-0.0 1 4-0.0 1 9 0.0 1 5 1 0 7 4 1 0 7 5 4 7 4 8 0 2.2 6 7 4 7 4 8 0 6.9 1 7 2 8 1 2 3 7 8.1 7 2 8 1 2 3 7 8.2 8 1 6.6 1 9 1 6.8 8 8 1 6.6 9 9 1 6.8 6 5 0.0 8-0.0 2 3 1 0 7 6 1 0 7 7 1 0 7 8 1 0 7 9 1 0 8 0 4 7 4 8 0 7.2 5 0 4 7 4 8 0 5.3 5 1 4 7 4 8 0 8.4 2 8 4 7 4 8 0 9.7 3 1 4 7 4 8 1 0.3 9 0 2 8 1 2 3 8 0.2 2 2 8 1 2 3 8 2.6 6 2 8 1 2 3 8 2.2 2 2 8 1 2 3 8 5.4 6 2 8 1 2 3 8 7.1 8 1 7.1 2 9 1 6.8 6 8 1 7.0 7 2 1 7.3 1 5 1 7.2 9 6 1 7.1 2 5 1 6.8 7 8 1 7.0 6 6 1 7.3 1 9 1 7.2 8 8-0.0 0 4 0.0 1-0.0 0 6 0.0 0 4-0.0 0 8

分析表3得出：精密水准仪测量的高程和载波相位差分GPS测量出来高程之差值，最大为5 cm，最小为1 mm。

 

表4 19次不同时间观测的精密水准仪与载波相位GPS测量高程值比较Table 4 Comparison of elevation values measured by carrier phase GPS and precision level at 19 different times

  

观测时间 h-H r/c m最大值 最小值2 0 1 4-0 7-1 4 3.3 1.0 2 0 1 4-0 9-0 7 1.9 0.9 2 0 1 4-1 2-1 4 2.7 1.0 2 0 1 5-0 1-2 5 2.8 0.8 2 0 1 5-0 3-0 3 2.1 1.1 2 0 1 5-0 5-0 1 3.7 0.8 2 0 1 5-0 7-1 8 1.9 0.6 2 0 1 5-0 8-2 7 2.4 0.8 2 0 1 5-1 0-1 4 2.7 1.4 2 0 1 5-1 2-1 0 3.1 0.9 2 0 1 6-0 2-2 8 2.4 0.8 2 0 1 6-0 5-0 7 2.0 0.6 2 0 1 6-0 7-1 3 2.1 0.5 2 0 1 6-0 9-2 0 2.7 0.8 2 0 1 6-1 0-0 3 2.8 0.8 2 0 1 6-1 2-0 4 2.2 0.7 2 0 1 7-0 2-1 0 2.1 0.7 2 0 1 7-0 5-0 2 2.8 0.7 2 0 1 7-0 7-2 9 1.7 1.1

2014年7月14日上午第一次测量至2017年7月29日19次观测，载波相位差分GPS各次观测高程数据以精密水准仪测量高程数据比较的最大值和最小值（表4）。

分析表4得出：2014年7月14日至2017年7月29日19次观测中，相同点高程误差最大为3.3 cm，最小为0.5cm，全部在cm范围内。

4 结语

由于海岸风沙地貌受台风响应的变化量比较少，使用传统的观测方法无法对其点位坐标和相应的高程进行精密的观测，经过多站载波相位差分GPS技术求解坐标拟合参数及高程异常值的方法,观测的精度达到了厘米甚至毫米级的精度，通过用精密水准测量方法对相同的观测点测定其高程，并以拟合后的载波相位差分GPS测量的相同点位的高程进行比较，可以得出，其点位误差均在厘米级范围内，表明多站载波相位GPS测量完全满足了海岸风沙运动观测的精度要求。

企业作为纳税人需要依法履行缴纳税务的义务，同样也依法享有诸多税法所规定的权利。例如当企业收入不佳、资金周转困难的情况下，企业可以利用缓期缴纳税款的权力，保持企业的稳定，当企业摆脱困难时再行缴纳税务，类似的还具有陈述权、申辩权、行政复议权、请求国家赔偿权等权力。通过专门的税务会计可以切实有效得为煤炭企业争取合法权益，保证企业稳定发展，为我国经济发展做出贡献。

多站RTKGPS测量能够全天候实时地测定出空间任一点的三维坐标，基准站与移动台之间无须通视，同时在测量区域内无须为待测点做标志，直接采取导航的方式引导到待测点进行测量，无论在精度还是在工作方式上都是传统测量方法无法比拟的，实践证明也是可行的。

参考文献：

[1]吴正.风沙地貌与治沙工程学[M].北京：科学出版社，2003：1-40

[2]韩致文，董治宝，王涛.塔克拉玛干沙漠风沙运动若干特征观测研究[J]．中国科学（D 辑）2003，33（3）：255－263.

[3]KEE C，PARKINSON B W， AXELRAD P．Wide Area Differential GPS[J]．Journal of the Institute of Navigation，1991，38（2）：123-135．

[4]王广运.载波相位差分GPS定位技术.测绘工程[J].1999，8（1）：12-14.

[5]PAIELLI B，MENALLY D，BBCH R E．Carrier phase diffferential gPS of approach and landing：algorithms and preliminary results[C].Salt Lake City：Proc．of the sixth International Technical Meeting，1993：831-840．

[6]刘贤万.实验风沙物理与风沙工程学[M].北京：科学出版社，1995：1-20.

[7]LANGLEY R．Time Clock and GPS [J].GPS World Magazine，1991，11（2）：38－42．

[8]国家质量技术监督局.海道测量规范GB12327-1998[S].北京：中国标准出版社，1999.

[9]黄德全，董玉祥，哈斯，等.多站差分RTKGPS技术在海岸风沙观测中的应用[J]．中山大学学报（自然科学版），2007，46（4）：121-124.

 
杨国辉
《华南地震》 2018年第01期
《华南地震》2018年第01期文献

服务严谨可靠 7×14小时在线支持 支持宝特邀商家 不满意退款

本站非杂志社官网，上千家国家级期刊、省级期刊、北大核心、南大核心、专业的职称论文发表网站。
职称论文发表、杂志论文发表、期刊征稿、期刊投稿，论文发表指导正规机构。是您首选最可靠，最快速的期刊论文发表网站。
免责声明：本网站部分资源、信息来源于网络，完全免费共享，仅供学习和研究使用，版权和著作权归原作者所有
如有不愿意被转载的情况，请通知我们删除已转载的信息 粤ICP备2023046998号