system: [SYSTEM]

user: 3. 项目对比分析 3.1. 项目实施 3.1.1. 外业生产 为验证引进的Sonic 2024型多波束测深系统在水深测量方面的功能、精度与效益等，单位组织了相关试验，并结合河道项目，选择了长江中下游荆江的郝穴河段、沙市河段等位置进行了相关试验与分析。 本试验项目采用1954北京坐标系，横轴墨卡托投影3˚分带。水下平面利用GNSS RTK进行定位，并采用七参数进行WGS-84坐标到BJ-54坐标的转换，每个试验河段的起算基准均一致。 试验外业生产主要包括：多波束系统设备安装、校正、施测及数据处理等步骤，单波束测深和RTK测量为传统精确测量方式。 3.1.2. 数据处理 单波束测深和RTK测量为常规测量，数据处理包含数据传输、水深处理、水位推算等。多波束系统数据处理包括数据预处理、数据成果输出和成图几个部分。 多波束系统数据预处理主要包括定位数据处理，声速剖面数据处理，潮位数据处理，姿态数据处理，深度数据处理和数据编辑、去噪、合并、清项；图型处理是对预处理后得到的水深数据进行网格化，生成数字地形模型(DTM)，最终形成河道地形图，数据处理流程如图2所示，可见其处理流程复杂，某一方面出错，则会影响整个数据的精度。 图2. CARIS处理多波束数据流程 经过数据处理，最终输出的部分河道3D影像如图3所示。 图3. 河道3D影像生成示意图 3.2. 数据统计 3.2.1. 数据统计原则 1) 为了有针对性的反映多波束系统在航道等重要区域的扫测结果，将比对区域内的所有数据都参与统计分析。 2) 多波束与河道固定护坡工程的对比，汛前采用传统RTK进行散点测量与汛期高水期的多波束扫射数据进行对比分析。 3) 点位对比的原则：多波束的水深数据按照单波束水道地形测量的计划线进行筛选，采用起点距相同的点计算高差的方式进行对比，或者多波束按照不同网格间距输出散点，通过程序寻找与单波束测点或RTK地形散点按照位置最接近点作对比。 4) 体积对比法原则：采用同一范围内的单波束与多波束测点分别成图后，再统计同面积范围内的体积法以及不同等深线包含体积的对比方式。 5) 多波束的水深数据筛选方式为不同点距间隔的方式抽稀输出。 3.2.2. 统计方法 为了更真实的反应多波束测量的精度，本文采用了不同类型设备、多种组合方式进行对比分析，水下高程点以单波束或RTK高程点为模板，统计方式为： 1) 水深互差统计，将多波束数据按照0.5 m间隔输出，然后将单波束按照亮点间距最小原则与其进行比较。 2) 体积法统计，主要有两种方式，选择多波束按照1 m、2 m、3 m、4 m、5 m不同点距与常规单波束进行同面积下的体积对比；同面积下31~29 m不同等深线的方式，进行比较两种仪器测量数据所包含的体积。 3) 固定的河道水泥护坡统计，汛期采用传统RTK散点测量方式，高水期采用多波束测量，再分斜坡与平坦两种区域进行统计分析。 3.3. 精度对比分析 3.3.1. 郝穴河段对比分析 本河段试验项目采用的比例尺为1:10，000中比例尺，试验长度约4.5 km。因本次试验为中比例尺测量，故单波束测量时未进行延时改正，且船速相对较快。 经过数据处理后，共统计了单波束HY1601同步试验数据23个断面456个典型的地形特征点。水深点高程互差统计以单波束测深数据为模版，统计数据见表1所示，通过对比分析最大互差为−3.03 m。 Table 1 Statistics of elevation difference of water depth points in 1:10，000 River topographic survey of Haoxue reach (m 互差分布 ≤±0.1 ±0.1≤ ∽<±0.2 ±0.2≤ ∽ <±0.5 0.52≤ ∽ <±1.0 ≥±1.0 统计个数 174 151 96 25 10 百分比 38.1% 33.1% 21.1% 5.5% 2.2% 表1. 郝穴河段1:10，000河道地形测量水深点高程互差统计单位(m) 取其中一个横断面，将单波束和多波束测量数据处理成水下河道测量横剖面图进行对比分析，见图4所示(红色为单波束测点，蓝色为多波束测点)。 通过表1和图4对比及分析可知： 1) 其误差分布符合正态分布规律。 2) 两者互差92.3%的数据小于0.5 m，且通过断面对比分析两者吻合较好，且在平坦地区，两种仪器水深数据互差较小，互差较大的基本出现在陡坡或水深突然剧变的位置。 3.3.2. 沙市河段对比分析 本河段试验项目采用的比例尺为1:1000大比例尺，试验过程中严格控制船行速度、严格进行延时改正，所有测量数据的起算与基准均一致。试验分三个区域，包括： 图4. 单波束和多波束同一断面测量对比图 第一区域：河道水下地形测量区域，进行单波束与多波束同步测深精度对比； 第二区域：固定护坡工程的斜坡处，高水期采用多波束扫射，低水期采用传统RTK散点测量模式； 第三区域：荆州长江三八滩平坦区域的固定维护工程处，测量方式为高水期采用多波束扫射，低水期采用传统RTK散点测量模式。 沙市河段三个区域试验数据对比分析如下。 1) 测深精度对比 第一区域数据经过处理后，共统计了单波束HY1601同步试验数据36个断面2344个典型的地形特征点。水深互差统计以单波束测深数据为模版，其统计结果见表2所示，通过对比分析最大互差为1.42 m。 Table 2 Statistical analysis of the difference between multi-beam and single-beam underwater elevation points in Shashi river reach (unit: m 互差分布 ≤±0.1 ±0.1≤ ∽ <±0.2 ±0.2≤ ∽ <±0.5 0.5 ≤ ∽ <±1.0 ≥±1.0 统计个数 705 613 910 107 9 百分比 30% 26.2% 38.8% 4.6% 0.4% 表2. 沙市河段与单波束水下高程点互差统计(单位：m) 通过分析对比，其误差的极值主要分布在岸坡区域。表2数据统计显示95%的互差在0.5 m以内，且极值分布优于郝穴河段，这主要因不同比例尺测量时，外业的设置、测量的船速、数据采集的时间同步性等有关。 2) RTK地形点对比 第二、三区域固定护坡工程RTK测量点与多波束扫射统计分析见表3、表4所示。分斜坡和平坦两种区域进行了统计分析。其中斜坡段最大互差为0.392 m，平坦区域最大互差为0.39 m。 Table 3 Statistical analysis of elevation difference of fixed slope protection project in Shashi river reach (slope) (unit: m 互差分布 ≤±0.1 ±0.1≤ ∽ <±0.2 ±0.2≤ ∽ <±0.5 0.5 ≤ ∽ <±1.0 ≥±1.0 统计个数 172 113 48 2 0 百分比 51.3% 33.7% 14.3% 0.7% 0% 表3. 沙市河段固定护坡工程高程点互差统计(斜坡) (单位：m) Table 4 Statistical analysis of elevation difference of fixed slope protection project in Shashi river reach (flat terrain) (unit: m 互差分布 ≤±0.1 ±0.1≤ ∽ <±0.2 ±0.2≤ ∽ <±0.5 0.5 ≤ ∽ <±1.0 ≥±1.0 统计个数 80 206 25 0 0 百分比 25.8% 66.2% 8% 0% 0% 表4. 沙市河段固定护坡工程高程点互差统计(平坦) (单位：m) 通过表3和表4的对比分析可知： 1、其误差分布均符合正态分布规律。 2、通过表3及表4的固定护坡工程对比，85%以上的互差在0.2 m以内，且平坦区域精度高于斜坡区域，这也表明采用测船进行水下地形测量时，其精度影响因素较多，而传统RTK测量模式则相对较少。 3) 体积法对比 根据沙市河段的比测数据，利用单波束测深仪数据绘制了5 m点距的地形图，多波束数据按1~5 m不同点距绘制了的地形图，分别建立DTM模型，计算相同面积下不同多波束点距或者多种等深线所围体积，并计算相对误差(多波束测深系统在不同点距，与单波束测深系统计算面积的百分比)，具体如表5、表6所示。 Table 5 Statistical table of volume difference between different points in the same area of 1:1000 River topographic survey in Shashi reac 方式 多波束测深计算体积 单波束测深仪 差值统计 点距 体积 测点间距5m体积 1 1.0 m 1，505，438.73 1，507，665.68 0.15% 2 2.0 m 1，512，611.21 0.33% 3 3.0 m 1，518，974.31 0.75% 4 4.0 m 1，523，209.32 1.03% 5 5.0 m 1，527，054.79 1.29% 表5. 沙市河段1:1000河道地形测量同面积不同点距的体积互差统计 通过表5对比及分析可知：若将单波束测量值作为真值，在进行多波束数据后处理时，点距取值越小则越接近真值。 Table 6 Statistical table of volume difference of different Isobaths in the same area of 1:1000 River topographic survey in Shashi reac 不同等深线 多波束测深计算体积 单波束测深仪计算体积 体积差 差值统计 33 m 1，092，619.87 1，093，578.32 −958.45 0.022% 32 m 969，966.34 970，615.24 −648.9 0.011% 31 m 848，915.08 849，300.29 −385.21 0.011% 30 m 730，759.64 730，843.03 −83.39 0.003% 29 m 615，346.49 615，143.73 202.76 0.008% 表6. 沙市河段1:1000河道地形测量同面积不同等深线的体积互差统计 通过表6的对比及分析可知：通过同面积不同等深线的体积互差分析，以及前面护坡的平坦与斜坡工程的对比分析可知，多波束在平坦区域测量与传统测量数据吻合良好，陡坡区域略差。