莱州湾数值模拟的快速建立研究.pdfVIP

莱州湾数值模拟的快速建立研究 李宝山，王红川 河海大学交通学院，南京（210098 ） E-mail 摘 要：本文主要针对莱州湾工程研究为主，当繁杂的数据资料充斥在工程建模当中，如果 希望充分利用数据资料而又不浪费大量的时间和精力在数据的处理上，那么，一套目的明确， 集成度高并且方便的方法便是不可多得的。介绍了工程前、后处理的高效方法和批处理程序 的应用，以及后处理中图形统一化的应用。其中涉及Visual Basic 和 Script 程序在工程的应 用。 关键词：数值模拟；批处理；Visual Basic ；Surfer，莱州湾 随着经济的发展，海岸及近海工程的日益扩大，实施要求的不断完善，工程领域对波浪 数学模型的建立提出了更高的标准。数值模拟在工程中从无到有，它的存在价值就是非常明 确的——运用数学方法去解决实际问题，即用数学的语言、方法去近似地模拟实际问题的过 程。而快速的建模就是提高工作效率的最直接的方法。 以下均以莱州湾的地形图为研究对象，分析、研究此方法的实用性。此海图比例为 1 ： 2 150000 的CAD 文件的海图，海图覆盖范围为 120000 ⅹ100000(m ) 。 1. 海图快速数字化 初步地数字化过程，就是将海图中的水深点、等深线、岸线等有用的参数以北京坐标的 形式确切地提取出来。首先，我们要利用大型扫描仪将纸质海图转化为图片格式，再将其调 整为 1：1000 的比例以便数据的读取。然后，以点和线的形式进行模拟。其中，由于方法不 同，工作者处理工作的效率会有天壤之别。 1）一般方法是：将海图的图片用CAD 或者 MAPinfo 等常用的软件进行读取，再对海 图进行微调；先经过角度的旋转，比例的缩放，坐标的平移，然后把每个水深点一一放到海 图对应得位置并输入水深值。平均每分钟可录入20-30 个。数字化完成后，还需用 Surfer 将 人工录入的数据按照随机插值并用等深线图来校验其正确性，对突变的地形数据进行核对和 修改。 2 ）快速方法是：配合 Surfer 的图形数字化功能进行快速提取。其过程为：首先，根据 一个参照点的坐标，计算出其海图四个轮廓点的坐标，并把这四个点做成一个数据文件并调 入 Surfer 中，我们便可以看见页面上有四个点了；然后，导入这张海图的图片，并将其拉伸 成四个点的大小，通过这个方法把海图放大到工程所需比例；最后，再用“数字化”命令可以 调出相应点击位置的坐标。我们按照等深线的方向把相同的水深点放到一组中，再整体赋值。 这样就实现了同一水深数据点的批处理。平均输出速度为一般方法的 2-3 倍。这样不仅输入 的更快，同时，不需要数据校验。因为批处理的点位置都是根据等深线的走势来分组处理的， 所以不会存在等深线间的突变的情况。 2. 插值方法 由于数值求解的精度主要取决于差分法的插值精度、网格尺度和迭代步长的选取，因此， 在制作地形的时候，合适的插值网格对数值求解有效果好，稳定性强等特点[3] 。在实际工程 中波浪是个复杂的随机过程，如果波浪数模计算值略偏大，那么需要充分考虑波浪的衰减， - 1 - 如：地形变化引起的波浪折射；受海底摩擦影响的波浪衰减；水深变浅而波浪进入极浅水域 的破碎；小风区的影响。由此可见，忽略了任何的一方面的考虑都会在最终的数值结果中反 映出来，而全面的对地形数据进行正确高效地处理可以有效地保证工程的质量。 在数据点分布不均，范围大的情况下，利用插值提高网格精度是工程数模中矩形网格 常采用的方便方法之一[2] 。目前，大量应用于工程的网格形式有：矩形网格和三角网格。对 于三角网格而言，在数据密集，不同形状的区域均可很好的模拟工程区域，特别是狭长的河 流水域工程。但对于一个长达数十公里的工程范围进行计算对于三角网格而言则需要更