海洋论坛▏贾永刚等：基于数字孪生的海底悬浮物时空变化可视化技术实现与应用 - 海底地质过程模拟工具

点击上方“溪流之海洋人生”即可订阅哦

一、概述

近年来海上油气开采平台、海底管线、深海养殖等深海工程大量兴建，然而海洋地质灾害的发生对海洋工程会构成严重威胁。深海原位长期监测能够揭示灾害触发机制及演变过程，进而有效预报海底地质灾害并减少相关损失，但是由于深海环境复杂多变，监测面临高测量精度、耐压耐腐蚀、抗电磁干扰、长期电能供应等一系列挑战，导致缺乏长期、连续的海底现场观测数据，结合现有的分析方法，数据难以直观表现海洋工程地质环境变化，更无法精准识别海洋地质灾害，严重制约深海工程建设发展。

" X& s4 S5 D" Q' }; f, U

数字孪生技术是以数字化方式创建物理实体的虚拟模型，通过深海原位长期监测数据与虚拟模型进行交互，为海洋工程地质环境的可视化分析提供技术支撑。2003年美国密歇根大学Grieves教授提出数字孪生概念，2010年美国国家航空航天局（NASA）引入数字孪生技术对航天系统的运行状态进行诊断与预测，确保航天器的运行安全。近年来，数字孪生与大数据、人工智能等新兴技术融合快速发展，广泛应用于工业、互联网、智慧城市等新型场景，国内外团队也引入该技术开展海洋工程领域的数字孪生平台研发。Coraddu等研发数字孪生船，结合监测数据估算出海洋污垢造成的油耗损失；Mayani等构建海上油井的数字孪生系统进行作业优化，实现异常状态的早期检测；Wang等研发海上风力涡轮机支撑结构的数字孪生系统，实现支撑结构运行状态的实时监测及故障诊断；张磊等提出海上石油平台工艺数字孪生系统方案，包含远程监控、虚拟巡检等功能；王文明等构建海洋修井数字孪生系统，实现修井管柱分析、辅助决策等功能。总体而言，数字孪生技术在海上风电、油井等海洋工程上进行了初步应用，但是此类数字孪生平台在构建海洋环境方面多利用虚拟数据，而非实测数据，造成数字孪生技术在海洋环境数据可视化分析上存在诸多问题。

" N/ ?0 |; `# ]6 z$ ^4 T$ ], ?

本研究构建的数字孪生系统结合复杂深海工程地质原位长期观测设备（SEEGeo）的深海底实测数据，将SEEGeo搭载的RBR多参数水质仪监测的浊度数据、海水抽滤测量的悬浮物浓度和近底层悬浮物粒径等数据相结合，使用Unity3D技术构建虚拟环境粒子系统，实现了南海北部陆坡神狐海域近底层悬浮物浓度升高的时空变化分析，提供了一种新的数据可视化分析手段，为海洋工程地质环境监测、灾害预测提供技术支撑。

二、数字孪生系统框架

' U4 b, Z8 _6 v0 Y. d+ F; O

图1海洋工程地质环境数字孪生系统框架

陶飞等在2019年提出数字孪生五维模型理论，内容包括物理实体、虚拟实体、孪生数据、数据连接、系统服务，对卫星、飞机、船舶、智慧城市等多个领域的数字孪生应用提供理论参考。本研究在此理论的基础上，构建了海洋工程地质环境数字孪生系统，包括物理层、传输层、虚拟层、应用层，如图1所示。

⒈物理层

" p. w/ z# @/ n7 L

物理层是指作业现场的环境，由各种功能子系统组成。本研究的物理层基于SEEGeo构建，由海底观测平台、海面中继浮标、室内监控平台3部分组成，工作模式如图2所示。工作原理为海底观测平台监测海底边界层温度、盐度、深度、压力、水体流速流向、海水浊度等环境要素，数据通过水声通信传输至海面中继浮标，再通过卫星通信将数据输送到室内监控平台。

+ k5 E& m2 D- c' B/ M9 V. u

图2 SEEGeo工作模式

1 i) M+ Z5 P: ~

⒉传输层

0 y' `- f' m0 z3 q, I' L

传输层实现监测数据与虚拟层的交互，构建物理层、虚拟层和应用层的联系。研究采用MySQL作为监测数据存储平台，在MySQL中设计节点数据表的逻辑结构，将数据在垂向上分成海洋环境数据、海底地形数据和海底地层数据3部分，通过Unity3DC#脚本将MySQL中的数据批量传输给虚拟层，在UI界面中实现数据的三维可视化分析，分析结果传输至应用层实现人机交互、辅助决策等功能。

9 F) x) \( v2 ]& ^

⒊虚拟层

2 v9 w4 X4 w0 b- R( M3 R

虚拟层是由几何模型、物理模型、行为模型构成的多维度模型，反映设备的实际工作状态。研究使用Unity3D技术构建设备几何模型，根据现实设备重力、浮力等力场属性构建设备物理模型，基于海洋工程学与流体力学的方法构建设备行为模型。通过虚拟现实技术（VR）增强实测数据在虚拟层的多维度表达，建立具有交互功能的海洋工程地质环境数字孪生系统。

3 B3 M, ^8 b1 b

⒋应用层

! A. z0 p- a$ l$ g0 \3 {, i

应用层具备智能监测、数据分析、人机交互和辅助决策等功能。智能监测功能监控设备的运行状况；数据分析功能实现监测数据的可视化分析；人机交互功能实现设备作业模拟和海底环境漫游；辅助决策功能帮助用户快速掌握海底环境异常信息，辅助用户做出决策，优化设备作业。

+ R9 a8 ^, n2 A* L0 A

三、虚拟层场景构建

虚拟层场景构建是数字孪生系统的关键，其精细化程度决定系统对现实环境的映射质量。本研究构建的虚拟层融合设备三维模型、虚拟仿真模型和虚拟现实交互模型，设备三维模型依据真实设备尺寸和材料构建，虚拟仿真模型依据设备实际工作状态构建，虚拟现实交互模型依据用户对现实环境的主观感受构建，最大程度实现虚拟场景对现实环境的映射，实现过程如图3所示。

图3虚拟层场景实现过程

) n: b% @$ z# ^, U7 k+ Y

⒈设备三维模型构建

本研究使用SolidWorks构建设备三维模型，包括设备平面几何、零件与装配体等信息，通过采集设备的平面几何信息构建各个零件模型，组装零件模型实现设备模型的装配及仿真演示。通过转换装配体模型格式、重新计算法线、删除冗余顶点，降低模型复杂度进而优化模型，最后导入Unity3D建立虚拟仿真模型。

. S+ v5 i% w* @5 [2 T' Y9 |

⒉虚拟仿真模型构建

本研究使用Unity3D构建虚拟仿真模型，包括设备模型和海洋环境仿真模型两部分。设备模型导入Unity3D后，根据设备的材质信息，在软件中设置不同材质球的着色器（Standard）参数赋予模型材质属性，在场景中添加定向光（DirectionalLight）和区域光（AreaLight）增强虚拟模型的显示效果，添加刚体（RigidBody）、碰撞体（Collider）等组件模拟设备在海底环境中受重力、浮力、碰撞等作用，如图4所示。

图4设备虚拟仿真模型

本研究使用实测的多波束数据建立海底三维地形模型，在Unity3D中使用Heightmaps方法创建地形模型，使用平滑高度（SmoothHeight）工具对模型进行平滑处理，进而实现海底三维地形模型的构建，如图5所示。

图5海底地形模型

6 _" z& [2 p2 _! e& p6 h# r

通过粒子系统（ParticleSystem）模拟浊度计现场记录的近底层浊度升高事件，设置粒子发射状态、运动周期、初始速度、重力等参数，模拟粒子在虚拟空间的运动行为，效果如图6a所示。使用天空、海洋等动态环境模块表现虚拟海洋环境细节，增强环境真实感，如图6b所示。

图6海洋环境仿真模型

$ R5 K% ^1 k: J) H, b, _) J1 m

⒊虚拟现实模型构建

为增强监测数据可视化分析效果，本研究采用虚拟现实技术进行人机交互。通过虚拟现实模型将监测数据可视化分析，加强用户对海洋环境信息的理解。虚拟现实技术结合计算机和现场监测数据生成具有交互功能的数字化环境，实现框架如图7所示。研究使用HTCVive设备，包括头戴式VR眼镜和操控手柄，将HTCVive连接至终端SteamVR，同时将SteamVRPlugin和VRTK插件导入Unity3D虚拟现实空间，通过编写VisualStudioC#脚本控制对象。用户使用操控手柄控制VR视角，触控进行空间移动、抓取特征、选中对象等UI交互，实现人机交互功能。

图7虚拟现实模型实现框架

8 A6 |* F: M5 G$ l& V4 b0 p

四、系统功能

! X4 {! _9 y1 N$ P7 e* A

本研究设计的数字孪生系统，使用Unity3D技术构建场景漫游、分析工具、告警记录等UI界面，结合MATLAB、ArcGIS数据分析技术实现智能监测、数据分析、人机交互、辅助决策等功能。

6 }7 C- L$ [3 z8 d

⒈智能监测功能

智能监测功能实现对设备运行状态的实时监控。应用倾斜姿态测量传感器数据判别设备的坐底姿态，应用北斗卫星通信数据实时定位设备，综合判断海底设备的作业安全状态。智能监测界面如图8a、图8b所示，用户通过UnityGUI图形控件切换地形、设备视角，远程监控设备作业过程，查看设备基本信息、当前状态、作业任务、执行情况等，查询设备部件实时运转情况及历史运行数据，实现设备的智能化监测。

% q. J1 [" v6 A# p* M

图8智能监测功能

⒉数据分析功能

系统数据分析功能在MATLAB、ArcGIS等数据分析方法的基础上，建立海洋环境、海底地形、地层等多空间维度数据处理模式。对于海洋环境数据，系统使用快速傅里叶变换频谱分析方法对压力数据进行滤波处理，采用平滑/高斯滤波器对浊度、温度数据进行降噪处理；对于海底地形数据，使用等值线、坡度、坡向分析等工具对海底DEM数字高程模型进行等高线加密、高程统计、坡度及坡向分析；对于海底地层数据，采用路径剖面工具快速提取DEM断面高程数据并生成剖面图。分析结果使用Unity3DEcharts分析工具进行可视化展示，如图9a、图9b所示。

. ]. m3 z# S/ v8 n& C, @! I+ m

图9数据分析功能

' M- G& Z: k1 M0 ?% O

⒊人机交互功能

+ ^ g8 I4 N8 U& \+ `: @

系统具有虚拟现实人机交互功能。通过在Unity3DCamera上编写VisualStudioC#脚本，用户能够在场景中进行旋转、平移、缩放等操作，点击相应控件切换场景，实现海底场景漫游功能。系统使用MySQL数据库实时记录监测数据及设备运行信息，用户通过访问数据库进行数据可视化分析、判断设备运行状态。人机交互功能实现用户在虚拟场景中分析数据，从用户角度出发，增强监测数据的可视化分析效果。

⒋辅助决策功能

辅助决策功能根据可视化分析结果辅助用户进行决策，系统将异常监测数据传输至虚拟层，虚拟层根据可视化分析结果做出辅助决策，将海底环境异常情况进行弹窗警告，使用户及时对监测设备下达应急操作指令并告知相关工程单位，提出应对措施以减少工程损失。

% S8 E3 _8 }; f$ ~5 |5 @3 i3 n

五、系统应用

/ O1 s- c+ B/ A; p1 N1 d8 f {

⒈研究区概况

" w$ q3 u( t5 M; \

本研究开展上述数字孪生系统在南海神狐峡谷区的应用设计与实践，研究区位置如图10所示。神狐海域位于南海北部陆坡的中部，构造上隶属于珠江口盆地最大的沉积凹陷——白云凹陷，该区域海底地貌复杂，发育17条SSE走向的陆坡限制型海底峡谷，平均坡度约为2°。研究区水深约为1300m，沉积物类型以细粒沉积为主，主要为粉砂、含泥质粉砂、泥质粉砂，该区域陆源碎屑来源丰富，沉积速率较高，且由于洋流、等深流和沉积物漂积等对碎屑物质的长途搬运，使得研究区内普遍发育滑塌沉积、浊流沉积。

3 {4 F/ h7 M$ ~* r: A) W

图10研究区位置

⒉研究方法

2020年9月，SEEGeo在南海神狐海域（水深1342m）由“海洋地质六号”船进行布放并开展200余天深海原位长期观测，搭载的RBR多参数水质仪监测到浊度、温度、盐度、压力、溶解氧浓度、电阻率等近底层环境要素变化。其中，水质仪在2021年3月7−21日期间监测到近底层沉积物悬浮导致浊度升高事件，但根据传统分析方法，数据只能以二维图表形式表现出浊度值随时间的变化，研究者难以将浊度值与近底层悬浮物的空间变化特征相联系，本研究使用数字孪生技术将浊度值变化与虚拟仿真模型进行交互分析，提供一种新的数据可视化分析手段，使沉积物悬浮事件在数字孪生虚拟平台进行分析展示。

7 ^* S5 h/ {/ _, Y' [

研究使用Unity3D技术构建虚拟环境粒子系统，对沉积物悬浮事件进行可视化分析。水质仪浊度数据量化悬浮物浓度的单位为NTU，使用最小二乘法建立浊度与抽滤测定悬浮物浓度的经验关系，将浊度数据反演为悬浮颗粒物浓度（单位：mg/L）；根据悬浮物的平均粒径（30μm），计算单个悬浮物的体积为14130μm3；根据悬浮物的平均密度（1.48g/cm3），计算在1m3空间范围内悬浮物的体积为3943mm3（以1mg/L为例）；根据单个悬浮物的体积，估算1m3空间范围内约有悬浮物数量279021059个；根据传感器的有效探测范围（125cm3）与海底摄像机的拍摄视野范围按一定比例（1∶4096000）进行放大，设置512m3的粒子体积空间；由于悬浮颗粒很少以单个颗粒的形式输运，而是以胶体、絮凝等复合形式输运，因此参照ROV拍摄近底层悬浮物直径（25.6mm）将单个粒子体积按比例（1∶621378370）放大，相关设置参数如表1所示。

. e) w! F1 I8 i" n5 t

表1近底层悬浮物颗粒虚拟仿真粒子信息表

在此基础上建立粒子在虚拟空间的运动行为，包括粒子发射状态、运动周期、初始速度、重力等。为了更好地描述近底层悬浮物的空间运动特性，根据Woodward研究可知悬浮物很少以单个颗粒的形式存在，而是以胶体、絮凝或聚集的形式存在，在系统中设置粒子碰撞器，使粒子间碰撞形成微团用于描述悬浮物颗粒的聚集过程。Gillard等的研究表明，沉积是海底悬浮物最重要的运移去向，在系统中设置碰撞检测，当检测到粒子与海底面碰撞时，使粒子消失，用于描述近底层悬浮物的沉降过程。编写VisualStudioC#脚本建立实测浊度值与虚拟环境悬浮粒子的对应关系，实现近底层浊度值的三维可视化分析。

+ S9 E' i }/ E' }( \; ?

⒊结果分析

1 k" E/ U* p% u

⑴浊度数据分析

在整个观测期内（图11a，2020年9月15日至2021年6月27日），RBR多参数水质仪监测到为期15d（2021年3月7–21日）的近底层浊度升高事件，如图11b所示。在近底层浊度升高事件发生之前，浊度值维持在0.56NTU左右，在2021年3月7日后的第1至9日，浊度整体呈锯齿状周期上升，上升周期在20h左右，每个周期浊度值增加近10NTU，然后快速下降到10NTU或12NTU，下降时间为1～2h，在此期间浊度值的峰值为49.69NTU，约是最初水平的88倍。观测区域位于峡谷底部，峡谷头部陆源碎屑不断向底部进行垂向输运，可能是造成此次近底层浊度升高事件的原因。后第10至15日浊度开始下降，最终回归之前的均值水平。

* T* I' X) k* z' ?; e9 S+ l0 H( A+ I

图11浊度随时间变化关系

⑵海底悬浮物浓度的反演

- ^6 w' t2 `- T6 A- k% z( S% W+ o

为实现数字孪生技术在此次近底层浊度升高事件的可视化应用，对目标区域悬浮物浓度及粒径开展研究，于2022年9月，使用“东方红3”号船CTD采水器（SBE911PLUS）对研究区5m、50m、100m、200m、300m、500m、1000m深度的海水进行抽滤分析，结果如图12所示，使用最小二乘法建立浊度与抽滤测定悬浮物浓度的经验关系，将浊度（Tu）反演为悬浮颗粒物浓度（SSC，mg/L），得到表达式为SSC=38.405Tu−32.570，相关性达到0.815。

图12浊度与悬浮物颗粒浓度之间的线性关系

. W+ X- x8 R9 D

⑶悬浮颗粒物粒径分析

7 P- |% H- z2 K( D# H X; M3 V

在CTD采水器上使用LISST-deep设备获取研究区近底层悬浮物平均粒径，图13显示LISST-deep所测悬浮物体积浓度在32组均匀粒径（1.36～230μm）下的垂直分布情况，并给出悬浮物的中值粒径（Md）曲线。发现研究区悬浮物粒径范围在16.458～57.479μm之间，平均粒径为30μm，这与Lokhov等在巴伦支海使用LISST-deep测量近底层悬浮物粒径的结果相似。在800～900m深度处的悬浮物体积浓度结果显著，存在一个明显的底部高浊度层，此时悬浮物体积浓度达到峰值，粒径为36μm，可能是该环境富存来自陆源输送的悬浮物、有机聚合物和生物碎屑。

图13悬浮颗粒物粒径的垂直分布及中值曲线

⑷海底悬浮物时空变化规律

* F- r h$ Y# {3 \: e4 p

为进一步探讨海底悬浮物的时空变化规律，本文对2021年3月15-16日的近底层悬浮物浓度升高事件使用数字孪生技术进行可视化分析。结果发现在虚拟环境粒子系统中，悬浮物浓度、粒子数量与聚集程度具有较大的时空差异。如图14所示，在0～12h期间，悬浮物浓度变化呈上升趋势，图14d中粒子数量相比图14a增加1.63倍，过程中粒子间相互碰撞发生聚集，在视野中粒子密集程度增高，其中前4h悬浮物浓度以每小时137.9mg/L的速率迅速上升，图14b视野中的粒子出现较为明显的增加，数量是原来的2.2倍，4h后上升速率减缓，以每小时39.5mg/L左右的速率波动上升，图14b至图14d中粒子总数维持较高水平，视野空间中粒子数目增加，形成多个粒子微团，悬浮物浓度达到峰值，Droppo认为悬浮物微团是矿物颗粒、有机物和生物有机体（如硅藻、细菌和藻类）的混合物，与悬浮物颗粒在物理、化学和构成上具有明显差异，在本研究观测到微团的形成可能是有机质与海底悬浮物强烈相互作用的结果。通常，悬浮物微团的密度较低，不易沉降，在图14d视野中粒子之间出现较明显的重叠，对设备的覆盖面积达到最大，峰值后悬浮物浓度以每小时240.2mg/L的速率波动下降，经碰撞器检测，大部分粒子与海底面碰撞后消失，如图14e所示，粒子数量与之前相比减少了64%，视野中粒子密集程度减小，空间内粒子相互重叠的现象减弱，微团生成数目减少，恢复到最初水平。

图14悬浮物浓度时间变化曲线及对应数字孪生虚拟粒子系统可视化分析图

⑸虚拟粒子系统可行性分析

& D! r4 _9 O; Q! f' y1 a6 |

为验证上述可视化分析方法的可行性，本研究将ROV海底摄像影片和数字孪生虚拟粒子系统可视化分析图进行对比分析。图15a至图15c为ROV于2021年6月22日挂缆回收设备时拍摄到的近底层悬浮物画面，使用二值化方法处理图像，通过直方图寻找二值化阈值，增强悬浮物在图像中的显示，对比分析ROV摄像影片和虚拟粒子系统可视化分析图中颗粒在图像面积上的占比，如图15a至图15c所示，虚拟粒子系统分析误差在0.16%～2.80%范围内，与ROV拍摄的情况基本一致，证明该分析方法具有较高的可行性。本研究通过建立数字孪生虚拟粒子系统，与深海原位长期监测的浊度数据相结合，能够展示近底层悬浮物的时空变化特征。

图15ROV海底摄像影片和数字孪生虚拟粒子系统可视化分析对比图

六、结论

⑴本文研发了海洋工程地质环境数字孪生系统，阐述了数字孪生系统的组成，基于深海长期原位观测数据，使用Unity3D技术构建了设备三维模型、虚拟仿真模型和虚拟现实模型等多模型融合的虚拟海洋环境，实现了对现实海洋环境的同步映射，同时结合MATLAB、ArcGIS数据分析技术开发了智能监测、数据分析、人机交互、辅助决策等功能。

6 y3 `/ K) L! y$ F) `! I, {. [

⑵本研究基于数字孪生技术对南海北部陆坡神狐海域峡谷区的沉积物悬浮事件展开研究，基于RBR多参数水质仪监测的浊度数据、海水抽滤测量的悬浮物浓度及LISST-deep获取的近底层悬浮物粒径等数据，使用Unity3D技术构建虚拟环境粒子系统，将SEEGeo监测到的近底层悬浮物浓度升高事件在虚拟环境粒子系统上进行可视化分析，提供了一种新的数据可视化分析手段。研究发现在整个观测期内悬浮物浓度、粒子数量与聚集程度具有较大的时空差异，其中，悬浮物浓度变化呈上升趋势时，视野中粒子间相互碰撞发生聚集，密集程度增高；悬浮物浓度维持较高水平时，视野中粒子数目显著增加，并发现粒子间形成多个微团，出现较明显的重叠现象，加剧了粒子对设备的覆盖程度。

y" C2 Y4 s& Q! ~9 ^( L

⑶经进一步验证，虚拟粒子系统分析误差在0.16%～2.80%范围内，具有较高的可行性。本研究通过建立数字孪生虚拟粒子系统，与深海原位长期监测的浊度数据相结合，能够展示近底层悬浮物的时空变化特征，为海洋工程监测数据可视化分析提供了新的方法，是数字孪生技术在海洋工程领域的一个应用实例。数字孪生技术在揭示地质环境对海底工程的影响方面优势突出，将温度、压力、流速等监测数据应用于系统进行综合分析，对探索海底地质灾害监测预警方面能够起到重要作用。

0 e) z# t4 M' @2 f; p6 x

【致谢】本研究的数据及样品采集得到国家自然科学基金重大仪器专项（41427803）和国家自然科学基金委员会共享航次计划项目（42149905）的资助。该航次（NORC2022-05）由“东方红3”号科考船实施，感谢徐景平、胡聪、孙中强、冯学志、孙均楷对本论文给予的帮助和支持，在此一并致谢。

; }* F+ W9 e' }% y

1

+ c! L- |/ H- Y: r0 I8 z

END

1

' P% l3 e2 e7 ^

【作者简介】文/乔玥单红仙王宏威朱超祁孙志文贾永刚，来自中国海洋大学山东省海洋环境地质工程重点实验室、中国海洋大学海洋环境与生态教育部重点实验室、中国海洋大学海底科学与探测技术教育部重点实验室。第一作者乔玥，1998年出生，女，江苏省徐州市人，研究方向为基于数字孪生技术的海底地质灾害监测预警；通讯作者贾永刚，1965年出生，男，吉林省伊通满族自治县人，教授，博士，研究方向为海洋工程地质与环境。本文受基金项目赞助，国家自然科学基金重点项目（41831280基金）、国家自然科学基金面上项目（41877223基金）、国家自然科学基金青年项目（42207173）。文章来自《海洋学报》（2023年第8期），参考文章略，版权归出版单位与作者所有，用于学习与交流，转载也请备注由“溪流之海洋人生”微信公众平台编辑与整理。

相关阅读推荐

海洋论坛▏贾永刚等：深海滑坡研究进展科技动态▏中国海洋大学贾永刚团队在内孤立波悬浮深海沉积物研究方面取得新进展政务信息▏第十五届光华工程科技奖评选启动，贾永刚、荆少东等47名专家在列科技动态▏国家重点研发计划“海底沉积物力学特性的原位测试装置”顺利通过海试验收海洋论坛▏深海可视化移动地质取样系统的应用与展望海洋论坛▏海底地形数据组织与可视化选取技术 1 s: |) Q! _; y; B4 L7 q

公众号

溪流之海洋人生

9 S/ Y; Z3 s) C; V2 s; j

微信号▏xiliu92899

用专业精神创造价值

( n; Q# O9 n: B3 E! h1 [0 y

用人文关怀引发共鸣

您的关注就是我们前行的动力 

投稿邮箱▏191291624@qq.com

9 I+ @8 K# a" t" X; |- Z ' v b4 n0 y; T; O/ ?& T& T