基于多波束技术的内陆水下地形三维测绘方法.pdf

理论研究水利技术监督 年第 期:/基于多波束技术的内陆水下地形三维测绘方法温 旋 孔德博 杜丽雯(江苏省工程勘测研究院有限责任公司 江苏 扬州)摘要:为解决现有水利工程地形三维测绘方法误差较大、适用性较差的问题 设计基于多波束技术的水利工程地形三维测绘方法 经过回波信号聚集和滤波匹配处理 得到测绘目标声呐图像 经过点云数据提取和散乱点云数据平滑处理 使用 三维软件构建测绘目标模型 实现水利工程地形三维测绘 实验结果表明:该方法测绘水下地形点位精度较高 适用性较强 呈现水下地形特征性能较好 可有效模拟水下三维地形 有效缩短工期为 推动工程的顺利进行关键词:多波束技术 水利工程 地形三维测绘 激励信号 三维建模中图分类号:文献标识码:文章编号:()收稿日期:作者简介:温 旋(年)男 工程师:我国领土面积较大 多山多水且东西向横跨地形复杂 使河流落差较高 为我国建设水电站提供了优良的条件 同时我国海岸线长达 万 海底油气储存量巨大 非常利于海洋勘探但在进行海洋油气勘探和构建水电站时 受河道和海底复杂环境影响 水利工程建设难度较大目前研究水利工程地形测绘方法很多 如张建芳设计的船载激光地形测绘方法 该方法使用激光扫描仪获取海岸地形范围 以海岸边缘信息为基础 使用摄像机拍摄海底图像 经过坐标系转换和点云数据配准后 获取海底等高线实现海底地形测绘 但该方法使用摄像机拍摄的海底图像存在大量噪声干扰 且该方法并未对图像噪声进行预处理导致其获取的海底地形等高线存在误差 精度不足 仲照东则设计了三维地形测量系统 该系统通过建立 地理信息数据库 使用无人机搭载激光雷达方式获取水底地形信息 经过数据融合后实现水底地形测绘信息 但该系统仅适用于 操作系统 不具备兼容性 导致其适用范围较小为了解决上述文献成果的不足 提高测绘误差和适用性 本文提出了基于多波束技术的水利工程地形三维测绘方法 多波束测量技术是通过透镜、反射等途径形成的多个元波束 利用该波束包围的区域获取测绘目标的技术 该技术被应用于地形测绘、水底基质分析及建筑工程定点等 具备较高的精度 三维技术则是通过对测绘目标进行多方位数据提取 确定其三维坐标的技术 该技术通过构建目标模型 为用户呈现可视化的测绘目标 为水利工程建设提供立体的数据基础本文将多波束技术和三维技术相结合 设计基于多波束技术的水利工程地形三维测绘方法 实验结果表 明:本 文 方 法 不 仅 将 测 绘 误 差 降 低 到 还使测绘目标要素三维符号化 可以在任何操作系统上完成测绘目标三维建模过程 具备兼容性 同时测量的水下地形距离不受环境温度、水密度及其他参量影响 适用性较强 为水利工程建设提供可视化的数据基础 水利工程地形三维测绘 多波束声呐成像模型构建 激励信号函数多波束声呐成像为视声呐成像技术的一种 其通过多波束声呐仪发射多波束声波 声波在水体内扩散时 遇到障碍物会形成散射波 散射波自其来向返回至多波束声呐仪接收器内 形成若干目标声呐回波 通过向测绘目标区域发射数个多波束声波 依据其若干目标回波信号形成测绘目标图像 在向测绘目标区域发射声波时 需设置其激励信号函数 以提升多波束成像声呐仪的成像质量 年第 期水利技术监督理论研究多波束声呐成像仪发射的声波包括单频矩形脉冲信号和线性调频脉冲信号 设置单频矩形脉冲信号和线性调频脉冲信号的数学函数表达式分别为:()其他()()()其他()式中 自然常数 单频矩形脉冲信号中心频率 脉冲宽度 参考点输出波形 脉冲周期 /为脉冲信号变化率 声波调频宽度利用公式()和()约束单频矩形脉冲信号和线性调频脉冲信号的频率与宽度 其回波信号更为均匀与稳定 生成的测绘目标图像更为精准 水底数学模型在上节生成的测绘目标图像基础上测绘水底目标时 将水底目标看作刚性球体 声呐基阵参考阵元由()表示 则测绘目标的数学表达式如下:()()()()式中 探测距离()多波束声呐基阵内参考点()测绘目标三维坐标原点设、分别为回波相位差和波束角 则探测目标坐标值()的表达式为:()多波束声呐基阵内参考点与测绘目标上任意点直线表达式为:()结合公式()和()得到联合测绘目标和直线方程的参数 表达式为:()()()()()()()()()()利用公式()获取声波辐射测绘目标区域 公式()内 值为 可获取到测绘目标的阴影覆盖区域 回波信号模型水底测绘目标阻挡声波形成反向声呐回波后由于声波过长或者测绘目标较大 导致声呐图像分辨率不够清晰 为此依据声呐图像分辨率将其划分为若干较小单元 依据较小单元内的声呐散射点描述声波受测绘目标阻挡后的回波强度(如图 所示)图 多波束声呐回波方程示意图由图 可知 在声呐传播时 为辐射声源级 为水面下方的声级 由于传输介质不同产生能量损失 为水面上方的损失声级 为反射点 结合声波、测绘目标和多波束声呐发射仪参数等 对当前多波束声呐仪接收声呐回波状态进行预判 回波信号损失后返回到基阵时的表达式如下:()式中 测绘目标强度利用声呐方程获取到测绘目标的声压数 令 为反射系数 为信号功率 测绘目标的反射回波声压表达式如下:()()()()式中 测绘目标的声压 扩展参数 介质吸收系数 校正系数 信号辐射区域跨度()指向性函数 阵元间隔 线阵列 波长 多波束声呐发射仪的预判角度预判角度正弦值经过上述过程 获取到测绘目标的反射回波声压 生成射频信号 作为回波信号 其中存在大量干扰噪声 需对回波信号进行聚集和滤波匹配处理 回波信号处理对多波束成像声呐模型获取到的测绘目标回波信号进行聚集和滤波匹配处理 使测绘目标声呐图理论研究水利技术监督 年第 期像特征更为清晰令 为基阵线阵列阵元 该阵元接收到的测绘目标区域分辨单元为()则该分辨单元的子单元的回波信号输出表达式如下:()()()()()()()式中()子单元的回波信号 分辨单元内子单元数量()测绘目标区域高分辨单元 的子阵元回波函数()阵元接收到的回波信号()测绘目标区域分辨单元的滤波匹配函数基阵所有单元接收到的回波信号聚集为多波束测深信号 其中 进行离散傅里叶变换后 通过相移方法获取第 个波束 第 个波束的含噪回波信号表达式为:()()()在回波信号相位噪声内引入模数转换器 分析离散梯度 进行聚集处理 输出积累后的回波信号为:()()()受回波信号内相位噪声影响 对聚集处理后的回波信号进行滤波匹配处理将回波信号输入至匹配滤波器内 利用该匹配滤波器输出的回波信号即为经过压缩处理的回波信号 令 为存在干扰噪声的复杂信号 其线阵列数值 此时该回波信号长度可由 /表示 回波信号的功率由 表示 则信号功率 和长度 与输出信号功率 和长度 关系如下:()通过上述公式可知 输出信号功率 约等于 复杂的回波信号经过匹配滤波处理后 信号长度得到缩减 而信号功率则得到提升 由此实现多波束声呐呈现的距离分辨率提高 得到增强后的测绘结果 测绘目标云数据提取经过 小节和 小节获取到测绘目标图像后 以测绘目标声呐图像内的像素点作为点云数据 并将该图像内的测绘目标轮廓像素点换算成三维空间坐标由于测绘目标的声呐图像为水体成像 在提取测绘目标时 仅能提取其二维坐标 此时利用新坐标参数 将测绘目标的二维转换为三维 该转换过程称为点云数据拼接 其步骤如下:在点云空间三维坐标系内 横轴与垂直轴均对应任意单一测绘目标声呐图像 垂直坐标为垂直方向 纵坐标为法向量方向 则纵坐标方向为多波束声呐发射方向 此时多波束声呐发射方向表达公式如下:()()式中 多波束声呐发射方向的纵坐标轴 声呐移动速度 多波束声呐获取测绘目标图像间隔时间利用公式()可将测绘目标的二维信息转换为三维信息 完成测绘目标点云数据提取 散乱点云数据平滑处理测绘目标点云数据提取完成后 测绘目标边缘的点云数据点存在散乱情况 影响测绘结果精度 需对测绘目标边缘的散乱点云数据进行平滑处理 其步骤如下 特征区域判定令 为散乱点云数据内的测绘点()为该测绘点的单位法矢()为测绘点的邻域中点集合 该集合内第 个元素 与测绘点 的法矢夹角由 表示 该法矢夹角的变化率表达公式如下:()/()()式中 元素的法矢夹角通过设定法矢夹角阈值 对比 与阈值大小当 小于该阈值时 则判断点位于特征区域内 测绘点法矢量修正计算测绘点法矢夹角变化率 后 其中 计算公式如下:()()()()式中()测绘点顶点修正函数()测绘点法矢夹角计算函数 调整测绘点位置利用公式()结果对测绘点位置进行调整其表达公式如下:()/()式中 未调整前测绘点位置 三角片 年第 期水利技术监督理论研究测绘点位置调整函数 三 角 片 面 积()测绘点投影 经过上述步骤可实现散乱点云数据平滑处理 水底地形三维建模将经过平衡处理后的点云数据输入至 软件内 建立测绘目标三维模型并将其保存为 格式文件 使用 三维软件即可调取该 格式文件创建点云测绘目标三维模型 经过连接点、线使测绘目标要素三维符号化 可以在任何操作系统上完成测绘目标三维建模过程 具备兼容性 利用该模型实现水利工程地形三维测绘 实验分析以某省适合构建水电站河流水域为实验对象使用 仿真软件模拟水域环境 使用本文方法对其水下地形展开测绘 测绘目标距离测试以多波束回波信号的延时和距离数值为基础对多波束回波信号的距离和延时关系进行拟合 测试本文方法的测绘水下地形距离能力 结果如图 所示图 水底地形距离测绘结果拟合曲线分析图 可知 初始回波信号数据点的距离随着其延时时间的增加而增加 初始数据分布为散点状态 经过本文方法处理后 拟合曲线呈现平滑状态 该结果表明本文方法测量的水下地形距离精度高的同时 不受环境温度、水密度以及其他参量影响 对水底地形距离测绘结果较为精准 且环境适用性较强 多波束点位误差测试使用多波束测绘水下地形时 需对水下地形定点 在 仿真软件内设置 个水下地形定点位置 编号为 使用本文方法对其点位进行测绘 结果见表 表 多波束点位误差单位:定点编号已知点位坐标测绘点位坐标互差横坐标纵坐标横坐标纵坐标横坐标纵坐标 分析表 可知 使用本文方法对水下点位测绘时 测绘坐标与已知点位坐标相差不大 其中最大偏差出现在第 个定点位置的横轴方向 但偏差仅为 测绘误差百分比约为 该数值表明本文方法测绘水下地形误差较低 可有效作为水利工程建设技术含量较高的水下地形数据 散乱点云平滑处理散乱点云平滑处理是提高测绘数据质量的一种有效手段 以一组该水域水下地形声呐图像为实验对象 测试经过本文方法对其散乱点云平滑处理后 声呐图像呈现效果 结果如图 所示图 声呐图像平滑处理前(左)与处理后(右)的结果对比由图 可知 未经过平滑处理的声呐图像整体较为模糊 水下地形边缘线不够清晰 而经过本文方法对该声呐图像进行平滑处理后 该声呐图像的清晰度得到一定程度提升 且水下地形的边缘轮廓点均已被标出 经过平滑处理后的水下地形声呐图像更能呈现其自身特征 上述结果表明:本文方法处理水下声呐图像平滑效果较好 也从侧面说明其水下地形测绘能力好 水下地形三维模拟效果使用本文方法模拟该水域某区域水下地形 模理论研究水利技术监督 年第 期拟结果如图 所示 可有效模拟水下地形 且依据颜色标准可分析出该水域水下地形起伏状况 为水利工程建设提供良好的参考条件图 水下地形三维模拟为更清晰呈现本文方法测绘的水下地形效果依据图 的水下地形三维模型 绘制该水域水下地形剖面图(如图 所示)图 水下地形剖面图分析图 可知 该水域水下地形在 轴方向表现为深谷 峰 谷的形态 深度数值最大出现在 轴方向约为 处 谷深度约为 自 轴方向延伸 在 轴方向得 位置表现为峰状 而后又呈现谷状 通过对比图 内等高颜色分布情况 水下地形剖面图与三维模拟图表现的水下地形情况较为相符 该结果表明:本文方法不仅可以有效模拟水下地形 还可呈现水下地形起伏情况 具备较好的应用性从工期角度测试本文方法应用效果 以水利工程一个桥墩施工工期为例 测试预计工期和使用本文方法对水下地形测绘后的施工工期变化情况(如图 所示)图 施工工期分析图 可知 该水利工程在钢围堰建设阶段、抽水阶段以及注浆阶段总工期为 而应用本文方法后 该水利工程桥墩建设工期仅为 有效缩短工期天数为 该结果表明:本文方法应用后 缩短水利工程建设工期效果较好 结论本文将多波束技术和三维建模技术相结合 设计基于多波束技术的水利工程地形三维测绘方法该方法环境适用性较强 能够为水利工程建设提供重要的数据 且稳定性良好 可以有效地平滑掉噪声 提供清晰的水下声呐图像 提高水下地形测绘水平 准确且完整地呈现水下地形起伏情况 测绘水下地形 能够有效减少工期天数 大大提高了水利工程建设的效率 但由于多波束声呐工作原理的局限性 测量范围受到一定的限制 在水深较大的情况下 可能无法获取到足够的数据 在后续研究中 可以尝试新型多波束声呐技术 如宽带多波束声呐等 以拓展测量范围和提高测量精度参考文献 李铁 周丰年 赵建虎.多波束系统安装偏差整体校准的地形特征匹配方法.测绘学报 ():.赵文超 刘满杰 王国岗.智慧水利中地质与测绘技术信息化提升与 应 用 .水 利 规 划 与 设 计 ():.张建芳.基于船载激光扫描技术的海岸地形测绘方法.舰船科学技术 ():.仲照东.基于 技术和激光机载技术的三维地形测量系统.激光杂志 ():.王皓冉 汪双 陈永灿 等.无人机测量系统在消力池水下地形测量中的应用.水力发电学报 ():.曹斌 朱述龙 邱振戈 等.影像双介质摄影测量的浅海地形测绘试验.遥感学报 ():.王晨 林杰.基于多波束测试系统对桥梁水下基础、河床检测方 年第 期水利技术监督理论研究法研究.公路 ():.张超 吕连港 姜莹 等.基于多波束测深仪的西太平洋声散射层测量.海洋科学 ():.刘尚蔚 周钰航 魏群.基于 文件的水利工程地形模型构建方法研究.人民黄河 ():.杜九博 杨士瑞 吴先敏 等.基于三维精细模型的水利工程高边坡稳定分析.人民黄河 ():.贾正元 崔斌 张贵宾.基于数字高程模型数据绘制立体地貌图方法.科学技术与工程 ():.梁冲.基于机载 与无人船多波束的水陆一体化测量技术研究.水利技术监督():.牛作鹏 李国杰 郭涛.技术在航道水下地形空间监测中的应用.水运工程 ():.何亮 李端有 刘源 等.基于 和多波束技术的库容计算及应用.人民长江 ():.张静远 徐振烊 王新鹏.基于 算法的水下地形辅助导航误差研究.海军工程大学学报 ():.(上接第 页)率 二次絮凝还可以改善水质的 值和电导率 有利于水中污染物的去除 因此新型硅系絮凝剂二次絮凝对、等指标的去除率比一次絮凝效果好 可能是因为二次絮凝增强了新型硅系絮凝剂的桥连作用、中和作用和卷扫作用 使得水中更多的污染物被吸附或沉淀出来图 絮凝次数对絮凝效果的影响 不同絮凝剂对絮凝沉淀速度的影响不同絮凝剂形成的絮体的沉降速度 见表 新型硅系絮凝剂的沉降速度最大 三氯化铁次之 最小 且新型硅系絮凝剂的沉降速度远大于规范推荐值 新型硅系絮凝剂在处理臭黑水时效果非常好表 不同絮凝剂形成的絮体的沉降速度项目新型硅系絮凝剂三氯化铁沉降时间/沉降速度/(/)规范推荐值/(/)结论文章以某黑臭河道为研究对象 分析不同类型和投加量的絮凝剂对黑臭水体处理效果 得出以下结论()为保证、的去除效果最好投入量应控制在最低/最高/采用新型硅系絮凝剂处理后 臭黑水中的 和 均达到一级 标()尽管在高 条件下硅系絮凝剂具有更好的去除磷效果 但是过高的 值可能会导致硅系絮凝剂分解或水解 从而降低其去除效率 因此最优 范围为 ()在 时、和 的去除效果比较好 污泥回流比控制在 时 的去除效果较好 二次絮凝可以增加絮体的密度和强度提高沉降速度和固液分离效率参考文献 吴卓阳.黑臭水体成因及相关修复技术的探究.辽宁化工():.张进武.水处理絮凝剂研究现状与前景.山西化工():.刘明源 钱佳 沙赟颖.胞外多糖生物絮凝剂与化学絮凝剂复配处理印染废水的效果.化学与粘合 ():.姚航 胡真虎.浅析国内黑臭水体修复方法进展.绵阳师范学院学报 ():.林聪.生态景观堰在河道生态治理工程中的应用探讨.水利规划与设计():.袁存明.生物修复技术治理河道底泥污染应用研究.水利技术监督():.