《Ocean Modelling》：OCSMesh and an Automated Creek-to-Ocean Mesh Generation Workflow
编辑推荐：
尽管在自动化与可复现性方面已有大量努力，非结构网格（Unstructured Mesh/Grid）生成仍是海岸海洋数值模拟中的一大"灰箱（gray area）"。本文介绍OCSMesh——一个开源、并行化的Python框架，用于自动化非结构网格生成。OCSMe
尽管在自动化与可复现性方面已有大量努力，非结构网格（Unstructured Mesh/Grid）生成仍是海岸海洋数值模拟中的一大"灰箱（gray area）"。本文介绍OCSMesh——一个开源、并行化的Python框架，用于自动化非结构网格生成。OCSMesh利用几何对象与尺寸函数（sizing function）对象定义计算域与分辨率，主要面向SCHISM（Semi-implicit Cross-scale Hydroscience Integrated System Model，半隐式跨尺度水科学集成系统模式）建模系统的应用。与传统网格生成器强制要求单元正则性及水深（Bathymetric）平滑不同，OCSMesh优先考虑数字高程模型（Digital Elevation Model, DEM）保真度，可在不限制单元形状的前提下解析小尺度水力特征。研究人员通过两个算例验证OCSMesh的适用性：其一，重建一套基于SCHISM的三维业务运行预报系统所用的跨尺度网格，结果表明基于OCSMesh的模型可在极短网格开发时间内复现业务模式的模拟技巧（skill）；其二，展示高效的OCSMesh网格合并流程可用于可迁移（relocatable）预报场景——研究人员可在合理时间（＜2小时）内为大河流域生成高分辨率网格并将其无缝并入较粗的业务网格。这些应用表明OCSMesh可提供更为真实的洪水预报，对及时的跨尺度洪灾减灾至关重要。
本文解读对象为Cassalho F, Mani S, Moghimi S, Ye F, Zhang Y J发表于《Ocean Modelling》的论文"OCSMesh and an Automated Creek-to-Ocean Mesh Generation Workflow"。
一、研究背景与立项缘由
非结构网格（Unstructured Grid, UG）模型因可按区域重要性灵活局部加密或粗化以更好捕捉地形，日益广泛地应用于溪-海（creek-to-ocean / riverine-to-ocean）跨尺度海岸海洋模拟。然而网格生成是海洋环流模式开发中最耗时且影响深远的步骤，常依赖人工操作、图形用户界面（GUI）及主观判断，阻碍自动化与可复现性。现有自动化UG生成软件包虽在一定程度上缓解了此问题，但仍面临"三角悖论（trilemma）"——需在精度（真实反映地形水深特征）、计算成本（节点与单元数量）及网格质量（如三角形接近等边程度）间取舍。因数值稳定性与计算成本通常不可妥协，最终常以牺牲网格对DEM的保真度为代价。此外许多自动化包为保证目标海洋模式的数值稳定而内置等边三角形约束和平滑过渡，甚至在水深插值前对DEM做平滑处理，这不仅引入不必要细化推高计算开销，还因偏离真实DEM可能在模式校准时被错误补偿，影响模拟物理真实性。因此亟需一种优先保证DEM保真度而非严格几何正则性的自动化网格框架。
SCHISM（Semi-implicit Cross-scale Hydroscience Integrated System Model，半隐式跨尺度水科学集成系统模式；Zhang and Baptista, 2008; Zhang et al., 2016）作为社区开源模式，允许倾斜（skewed）水平网格单元及三角/四边形混合单元，可接受DEM原始起伏而无需水深平滑与人为网格操控，其水平网格宽容度配合可配置垂向分层方案使单一网格可支持1D/2D/3D多形态配置。依托SCHISM对单元变形的容忍能力，NOAA海岸测量发展实验室（Coast Survey Development Laboratory, OCS）开发并维护OCSMesh——面向SCHISM应用的自动化非结构网格生成Python工具包，强调DEM保真、跨尺度（小溪系统至开阔海域）连续网格构建及不同子域（河流、洪泛平原、海洋）独立生成与无缝合并，适用于业务化模式框架与可迁移海洋模拟平台。
二、主要关键技术方法概述
研究人员开发开源Python包OCSMesh（v2.0+解耦网格引擎依赖，支持Gmsh等），以DEM驱动为主（亦支持shapefile或已有网格），通过几何（geometry）对象和尺寸函数（sizing function / Hfun）对象定义计算域边界与空间分辨率分布；利用SCHISM对skewed三角形及四边形混合单元的兼容性放宽单元正则约束以保DEM原貌；实现独立子域网格分别生成后基于高效算法无缝合并（merge），形成连续跨尺度非结构网格。验证采用美国NOAA Surge and Tide Operational Forecast System——STOFS-3D-Atlantic业务预报系统对应的大西洋岸域为基准：例1用OCSMesh自动重建STOFS-3D-Atlantic所用跨尺度网格并对比模拟技巧；例2以Pearl River Basin（路易斯安那–密西西比州间）为对象，借助RiverMapper获取精细河道划分产品，用Hfun加密生成5万米级洪泛平原高分辨率网格（22万余节点）并在<2小时内合并入较粗业务网格，演示可迁移高分辨率嵌入能力。算例流程复现教程与Jupyter Notebook公开于OCSMesh的GitHub仓库。
三、研究结果
Description of Application Cases（应用案例说明）
复合洪水（compound flooding）由河流（fluvial）、降雨（pluvial）与风暴增水/潮汐（coastal surge/tide）共同作用，精确模拟需能无缝覆盖溪-海尺度的网格。研究人员据此设计两案例：基于STOFS-3D-Atlantic域评估OCSMesh重建业务网格能力，以及在大域中局部嵌入高分辨率流域网格的可迁移性。
Mesh Engines（网格引擎）
OCSMesh v2.0起解除对特定网格引擎绑定，用户可选适配的求解器（当前v2.1+支持Gmsh等），保持DEM主导流程并兼容skewed三角与四边形单元，体现框架灵活性。
Model Pre-Processing（模式前处理）
为使结果可复现，研究人员编写逐步Jupyter Notebook教程涵盖洪泛平原、河道及海洋网格的生成与合并全流程（符合文中设计准则与工作流程），教程公布于OCSMesh的GitHub供用户参照。
Relocatable Mesh Generation（可迁移网格生成）
针对同计算域内局部提升河网分辨率的需求，研究人员用RiverMapper产出精细河道划分，为Pearl River Basin建高分辨率河道网格（223,349节点），并利用Hfun加密洪泛平原至约50 m分辨率；演示OCSMesh可将此高分辨率子域网格在<2 h内合并入原有较粗业务网格，验证了可迁移高分辨率嵌套方案的可行性及效率。
Discussion and Conclusions（讨论与结论）
研究人员指出OCSMesh是数据驱动、并行、面向对象、面向SCHISM应用的自动化非结构网格生成包，优先DEM保真并允许skewed单元。自2020年发布以来已有活跃用户群用于真实业务网格构建（Mani et al., 2022; Sun et al., 2022; Martins et al., 2024）。两案例证明：①OCSMesh能以远少于手工方式的时间重建跨尺度STOFS-3D-Atlantic网格并使SCHISM模式维持与原业务网格相当的模拟技巧；②OCSMesh的高效子域独立生成与合并功能支持在大域中快速嵌入局部高分辨率河-洪泛平原网格（如Pearl River Basin），耗时<2 h，适合可迁移业务预报部署。结论为OCSMesh通过减少人为干预、保留DEM原始特征及支持跨尺度无缝合并，可提升复合洪水模拟真实性，为及时洪灾风险评估与减缓提供实用工具。OCSMesh开源发布于https://github.com/noaa-ocs-modeling/OCSMesh，示例数据与流程可于该仓库获取；SCHISM源码见https://github.com/schism-dev/schism。
四、研究意义
该研究表明，在SCHISM这类对skewed单元具包容性的非结构网格模式支撑下，放弃传统对单元正则性与DEM平滑的过度追求、转而优先保持DEM保真度，既可降低网格开发主观性与耗时，又避免因人工修饰地形带来的模拟偏差。OCSMesh的自动化、并行化及子域合并机制为跨尺度（creek-to-ocean）海岸海洋模拟、尤其是复合洪水与业务化可迁移预报系统的网格构建提供了可复现、高效率的解决方案，对推动海岸数值模拟工作流标准化与灾害预警能力提升具有重要实用价值。
