海岸海洋潮流模拟可视化与虚拟现实建模

《海岸海洋潮流模拟可视化与虚拟现实建模》由会员分享，可在线阅读，更多相关《海岸海洋潮流模拟可视化与虚拟现实建模（8页珍藏版）》请在装配图网上搜索。

1、海岸海洋潮流模拟可视化与虚拟现实建模实验室论文编号：SCIEL21101102基金项目：国家自然科学基金资助项目(49701013)；“九五”国家科技攻关项目(96-922-03-01) 作者简介：马劲松(1969-),男(回族),南京人,南京大学副教授,博士,2000年赴加拿大滑铁卢大学作为访问学者,著有地 理信息系统与空间可视化主要从事GIS空间数据结构理论、技术及其在海岸海洋方面的应用研究。马劲松，朱大奎南京大学大地海洋科学系，南京大学海岸与海岛开发教育部重点实验室，南京，210093电话：86-25-3593770(0), 86-25-3593756(H), Email： majs)V

2、isualization Modeling and Virtual Reality for CoastalOcean Tidal Current Numerical SimulationMA Jingsong,ZHU Dakui(Dept. of Geo & Ocean Sciences, Nanjing University,The Key Laboratory of Coast and Island Development of Ministry ofEducation P.R.China Nanjing, 210093,Tel: 86-25-3593770(O), 3593756(H),

3、 Email majs )Abstract: The technology of scientific visualization and virtual reality (VR) for tidal current numerical simulation is one of the most important modern methods utilized in coastal ocean science. This paper firstly illustrates a data model used for tidal current simulation from the view

4、point of Geographical Information Systems (GIS), presents an object-oriented modeling methodology suitable for hydrodynamic numerical simulation. Secondly, discusses the field visualization technique of tidal current velocity and describes the utilization of VR in simulation. The authors also design

5、ed and implemented a prototypic software system-VR0cean, used in visualization of coastal ocean tidal simulation as well as oceanographic studies, and verified the results in the research program of tidal current numerical simulation on the sand ridges in South Yellow Sea of China.Keywords: Tidal Cu

6、rrent Numerical Simulation; Virtual Reality; Visualization Model; Coastal 0cean摘 要 : 潮流数值模拟的科学可视化与虚拟现实是海岸海洋科学研究重要的现代技术手段之 一。本文首先从地理信息系统的角度对潮流数值模拟系统的数据模型进行了论述,提出了适 合水文数值模拟的对象模型方法；其次探讨了空间流场的可视化技术,阐述了虚拟现实在潮 流数值模拟中的作用；最后设计并实现了一个海岸海洋潮流模拟虚拟现实的原型软件系统VROcean,并在南黄海辐射沙洲等的潮流数值模拟试验中进行了实际对比检验。关键词：潮流数值模拟；虚拟现实；可视化

7、模型；海岸海洋1概述潮流是海岸海洋主要的水动力因素之 一,它影响着诸如泥沙、盐分、各类污染物 及能量的输运过程12。潮流的数值模拟是 海岸海洋科学研究中一个重要且经济快速 的技术方法,它通过离散化描述潮流运动的 控制方程,用近似的求解方法模拟潮流的运 动规律,达到对海岸海洋潮流场特征及其变 化的科学认识,并指导对海岸海洋的管理、 开发和利用5。然而,数值模拟方法在具有经济快速特 点的同时,也存在形象性与直观性差的问 题。潮流数值模拟的计算结果往往是大量枯 燥的数据,以及根据这些数据所绘制的静态 流场图和水位过程曲线,这就阻碍了进一步进行潮流运动规律性认识的研究工作。针对 这一问题，科学可视化和

8、虚拟现实技术的兴 起，无疑是改变这一状况的契机。科学可视化和虚拟现实为研究者提供 了直观地处理研究结果的技术方法，越来越 被公认为是科学研究过程的重要组成部分。 它在计算机空间(Cyberspace)中为研究者 开辟了一个具有沉浸感的虚拟环境，实现了 三维空间和时态数据的可视化，并使研究者 既能够在虚拟环境中交互地操控研究对象， 更可以在仿真模拟等科学计算过程中实时 地得到正在处理的动态过程的反馈。如果将 科学可视化与虚拟现实技术运用于海岸海 洋潮流数值模拟中，将极大地提高海洋研究 者的研究成效。建立海岸海洋潮流数值模拟的可视化 与虚拟现实系统，关键技术在于数据模型的 建立。结合海洋学、水文学

9、、地理信息系统 和计算机图形学等理论和技术，该系统可由 图 1 所示的三个子系统模型组成：图 1 海岸海洋潮流数值模拟可视化与虚拟现实系统的组成Fig.1 The components of coastal ocean tidal current numerical simulation visualization and VR system连续方程：针对上述系统模型的建立，目前的国际 和国内研究者已分别在潮流数值模拟数学 模型、面向对象空间建模、以及科学计算可 视化与虚拟地理环境等方面都作了许多卓 有成效的研究和探索。例如 Bartlett6,7、 McCalMB和SherinMio用于海洋

10、GIS的线 性参照系统；Gold13,14等的用Voronoi 图和其对偶 Delaunay 三角剖分建立的三维 海岸海洋数字地形模型；王颖、朱大奎等的 南黄海辐射沙洲潮流通道边坡稳定性 4GIS 与遥感综合集成研究，利用南黄海海底沙脊 2潮流数值模拟系统的对象模型 潮流的数学模型：直角坐标系下三维 潮流数学模型基本方程是 Navier-Stokes 方 程组，水流间的紊动应力采用 Boussinesq 假 区的地形、潮流、沉积物的大量数据，以此 基础建立潮流泥沙模型；马劲松、王颖等用 于海岸动态研究的基于关系数据库的空间 时态数据模型3；林珲等的东中国海潮波系 统与海岸演变模拟；LiW181

11、9基于超图 理论的三维海岸空间数据模型；以及 Su2o 的海洋时态序列三维可视化系统研究与开 发等。本文结合上述研究成果，提出了基于UML 统一建模语言的潮流数值模拟可视化 与虚拟现实系统的面向对象系统模型。定，则基本方程为：dudvdw小+ + = 0dx dy dz动量守恒方程：duf 1 dPd (+ u V u f +dtH vdy 丿 dz v VM dz 丿)d(dv d(dv A IA丿dyIv H dy 丿dzv VMdz丿(dw)d(dw )AAV H dy 丿dz v VMdz丿p dxdx v H dx 丿dv1 dPd ( . dv+ v + f + I A dtH u

12、 p dydx v H dxdw r 1 dPd ( . dw+ u V w + f + g I A dtH w p dzdx v H dx 丿密度守恒方程：竺+W -dtt Qx(Qg )Q(QgQQg )A tAtAt H Qx 丿QyI H Qy 丿Qz(vgQz丿流体静力学方程：QP其中：p=p0+qp：水密度。卩=(u v, w)：速度矢量。f: Coriolis系数，f=23sin申g :地 转角速度，件纬度。P：静压力。g：重力 加速度。Ah：水平涡旋扩散系数，Avm:垂HVM直涡旋粘性系数，AV :垂直涡旋扩散系数。 潮流数值模拟模型：潮流数值模拟分 为二维和三维的模型。对于二

13、维模型，按差 分网格形状可以分为三角形、正方形、矩形、 多边形、曲线坐标网格以及各种形状网格的 组合等；按计算方法可分为显式法、半隐半 显式法和隐式法；对三维的模型，计算方法 有分层二维法、谱方法、流速分解法和坐标 变换法等。潮流数值模拟数学模型简化的对 象类层次可以用以下图2来表达。图2潮流数值模拟数学模型的简化UML类图Fig. 2 UML class diagram of tidal current numerical simulation hydrodynamic model3潮流数值模拟空间数据对象模型四维时空数据模型：潮流数值模拟的 空间与时态数据构成四维的多分辨率时空 数据结构，

14、可由三维空间模型与一维时间序 列组成。二维平面的网格结构可以采用栅格 表示矩形或正方形网格，不规则三角网 TIN 表示平面的三角剖分。多个二维网格组合成 三维序列用以表达海水的垂直分层结构，再 加上时间序列，就可以形成四维的潮流时空 数据模型。空间参照系作为时空的一部分与 空间数据模型相关联。实例化的边界条件、 水下地形、水位场和流速场等需要与可视化 模型相应部分相关联，以生成静态三维和动 态四维的潮流场虚拟环境。潮流数值模拟空 间数据对象模型的简化类图如图3 所示。图3简化的潮流数值模拟空间时态数据对象UML类图Fig. 3 Simplified object class diagram o

15、f spatiotemporal data model空间多分辨率存储的实现可以借助分 维的方法。海底地形的分维数反映了地形的 粗糙程度，运用空域到频域的离散傅立叶变 换求得分维数，再利用分维数控制地形的内 插和采样，就可以获得多分辨率的地形数 据。图 4 显示了同一海底地形不同空间分辨 率的模型。空间矢量的对象数据模型：模型中水 下地形与水位场都是数量场，而流速场是矢 量场。同样，潮流虚拟现实对象模型中，虚 拟环境的透视投影、表面渲染及模拟飞行的 实现也需要空间矢量对象的支持，所以在模 型中加入空间矢量的对象模型，如图5所示。 空间矢量及方向的操作函数如下表1 定义。高分辨低分辨图 4 多分

16、辨率多维时空数据模型Fig. 4 Multi-resolution and multi-dimensional spatiotemporal data modelVector+Position: Point +Magnitude() +Addition(Vector) +Scale() +DotProduct(Vector) +CrossProduct(Vector) 0, c 0, s.t.d c d + c d1 2 1 2 112 2Among丄厶丄厶丄丄厶厶 ii*=【1=dAmongd , d , d , if3c 0, c 0, c 0, s.t.d c d + c d + c d

17、123123112233方向vPoint, Front, Right, Above平移Translate (O,v) = Translate (OP，V),frOnt,right,ab两旋转f*Rotate(O,RotationMatrix) (OP, front ,right ,above Jnnn f(front , right ,above ) - (front, right,above) - RotationMatrix 空间参照系：潮流模拟计算网格的空 间定位需要用空间参照系来确定。空间参照 系是地球的数学模型，定义了地球参数化三 维表面的空间坐标系统，及各坐标系统间的 数学关系。在

18、较大空间尺度的海岸海洋进行 潮流数值模拟，地球形状和曲率是影响计算 精度的重要因素。同时，空间参照系也是将 潮流模拟模型、地球空间数学模型和可视化 虚拟现实模型结合在一起的空间定位基准。海岸海洋潮流数值模拟采用的空间参 照系中，为了建立统一的海陆空间定位基 准，需要解决以下几个方面的数学建模问 题。a.确定旋转椭球体的参数，包括长短半 径、扁率，以及大地测量的定位参数等等。 确立区域坐标系统与全球坐标系统 WGS84 等的转换数学模型；b.确定海岸带高程起算 面、水下深度基准面和平均海面等的空间数 学关系；c.解决国家标准地形图的高斯投影、 海图的等角圆柱（墨卡托）投影及国际常用 的UTM （

19、通用横轴墨卡托）投影之间的变 换关系。海岸海洋潮流模拟空间参照系可由 图 6 所示。图 6 海岸海洋潮流数值模拟空间参照系 UML 类图Fig.6 Object class diagram of spatial reference system4潮流数值模拟可视化与虚拟现实对象模型图 7 简化的潮流数值模拟可视化与虚拟现实对象模型 UML 类图Fig.7 Simplified object class diagram of tidal current numerical simulation visualization and VR model潮流数值模拟的可视化与虚拟现实系 统是“虚拟地理环

20、境”和“数字地球”的组 成部分，通过计算机系统的多媒体特性，能 向研究者提供对海岸海洋潮流对象直观空 间认知的手段，使研究者产生虚拟现实系统 的三个特征，即主体的沉浸感(Immersion), 运用对虚拟客体的交互操作( Interaction) 以及认识反映研究客体规律性的自治特性(Autonomy)。对于潮流数值模拟的可视化 和虚拟现实主要体现在具有真实感的动态 视觉模型和交互操作潮流参数而获得的虚 拟潮流客体自治的变化特性，其对象模型由 图7 所示。由于需要显示水下的地形起伏，海水必 须呈现蓝色半透明的视觉效果，所以，水位 场的显示需要用到色彩融合的技术。实现海 面和水下三维飞行模拟，需

21、要设定飞行器实 时的空间坐标及其空中姿态。这可以通过上 述的空间方向矢量 Orientation 来表达，其 中， Position 确定飞行器的空间位置，并假 设 Front 方向矢量确定飞行器的飞行方向并 控制飞行器的滚动角度；Right方向矢量控 制飞行器对地的俯仰角度； Above 方向矢量 控制飞行器的偏航角度；根据人机交互操作 发生的动作，实时刷新飞行器的空间位置和 空中姿态，就可实现实时的三维模拟飞行。综合上述可视化与虚拟现实技术，我们 设计并实现了一个海岸海洋潮流数值模拟 虚拟现实的原型软件系统VROcean,其中实 现了潮流水位场和流速场的三维实时动态 模拟，交互控制的三维空

22、间飞行，遥感正射 影像的三维地形映射，多分辨率DEM数据 的分形压缩等功能。图 8和图 9是用该软件 生成的黄海辐射沙洲潮流数值模拟中水下 数字地形模型和潮流场的虚拟场景。进一步 的工作将在于设计并实现交互控制潮流参 数获得虚拟潮流的自治特性的研究。图 8 黄海辐射沙洲水下数字地形模型的虚拟场景Fig.8 Virtual scene of sand ridges under the water of South Yellow Sea of China图 9 二维潮流场矢量的虚拟场景（涨潮流）Fig.9 Virtual scene of two-dimensional tidal current

23、 velocity field5 结论本文所述的研究将可视化和虚拟现实 的建模技术手段运用到海岸海洋潮流数值 模拟研究之中，其动态图像演示了潮流场的 运动演化过程，及海底地形与潮流场的相互 演化关系。这与海底地貌与沉积学研究得出 的结果可互为印证4，表明该项技术具有实 际应用价值。海岸海洋潮流数值模拟的可视化与虚 拟现实是海岸海洋潮流研究的重要技术方 法，借助于地理信息系统、计算机图形图像 学、以及科学可视化等理论和技术实现的虚 拟现实系统，可以很好地满足海岸海洋科学 的深入研究。在此基础之上，实现虚拟地理 环境、数字地球等网络化、全球化的系统将 是今后的发展方向。参考文献：1 WANG Yi

24、ng, REN Mei-e, ZHU Da-kui. Sediment Supply to the Continental Shelf by the Major Rivers of China J. Journal of Geological Society, 1986, 143: 935-944.2 WANG Ying, REN Mei-e, SYVITSKI J. Sediment Transport and Terrigenous Fluxes J. The Sea, 1998, 10: 252-292.3 MA Jing-song, WANG Ying. A Spatiotempora

25、l Data Model on Relational Database for Coastal Dynamic Research J. Marine Geodesy, 1999, 22(2): 105-114.4 WANG Ying, ZHU Da-kui, et al. The Sediment Characteristics and Evolution of Radial Sand Ridges in South Yellow Sea of China J. Science in China, 1998, D : 385-393. (in Chinese)5 LIN Hui, et al.

26、 Study on Simulation of Tidal System in East China Sea Based on GIS J. ACTA GEOGRAPHICA SINICA, 1997, 52 : 160-169. (in Chinese)6 BARTLETT D. Space, Time Chaos and Coastal GIS C. 16th International Cartographic Conference and the 42nd Deutscher Kartographentag, Koln, Germany, 1993. 539-551.7 BARTLET

27、T D, DEVOY R, MCCALL S OCNNNOR I. A Dynamically Segmented Linear Data Model of the Coast J. Marine Geodesy, 1997, 20: 137-151.8 MCCALL S. The Application of Dynamic Segmentation in the Development of a Coastal Geographic Information System C. CoastGIS 95 International Symposium on GIS and Computer M

28、apping for Coastal Zone Management, Cork Ireland, 1995, 305-312.9 SHERIN A G, EDWARDSON K A. Using GIS and Dynamic Segmentation to Build a Digital Coastal Information Database C. Coastal Zoon Canada, Halifax, Nova Scotia, 1994.10 SHERIN A G. Linear Reference Data Models and Dynamic Segmentation: App

29、lication to Coastal and Marine Data C. Marine and Coastal Geographical Information Systems, Taylor & Francis, 2000, 95-115.11 MANDELBROT B. How Long is the Coast of Britain? Statistical self-similarity and fractional dimension J. Science, 1967, 156: 636-638.12 MANDELBROT B. The Fractal Geometry of N

30、ature M. San Francisco: Freeman, 1982.13 GOLD C M. Problems with Handling Spatial Data: The Voronoi Approach J. CISM Journal, 1991, 45: 65-80.14 GOLD C M, CONDAL A R. A Spatial Data Structure Integrating GIS and Simulation in a Marine Environment J. Marine Geodesy, 1995, 18: 213-228.15 GOLD C M. An

31、Algorithmic Approach to Marine GIS C. Marine and Coastal Geographical Information Systems, Taylor & Francis, 2000, 37-52.16 WRIGHT D J. Spatial Reasoning for Marine Geology and Geophysics C. Marine and Coastal Geographical Information Systems, Taylor & Francis, 2000, 37-52.17 LI R, QIAN L, BLAIS J A

32、 R. A Hypergraph-based Conceptual Model for Bathymetric and Related Data Management J. Marine Geodesy, 1995, 18: 172-182.18 LI R. Shape from Shading: A Method of Integration of Sonar Images and Bathymetric Data for Ocean Mapping J. Marine Geodesy, 1993, 15: 115-127.19 LI R. Data Models for Marine and Coastal Geographical Information Systems C. Marine and Coastal Geographical Information Systems , Taylor & Francis, 2000, 25-36.20 SU Y. A User-friendly Marine GIS for Multi-dimensional Visualization C. Marine and Coastal Geographical Information Systems , Taylor & Francis, 2000, 227-236.