海湾纳潮量数值计算技术规程.docxVIP

1

海湾纳潮量数值计算技术规程

1范围

本文件确立了海湾纳潮量数值计算的基本程序，规定了资料收集、二维潮波数值模型、纳潮量计算、成果与报告等内容。

本文件适用于海域内的海湾纳潮量计算。

2规范性引用文件

下列文件中的内容通过文中的规范性引用而构成本文件必不可少的条款。其中，注日期的引用文件，仅该日期对应的版本适用于本文件；不注日期的引用文件，其最新版本（包括所有的修改单）适用于本文件。

GB/T12763.2海洋调查规范第2部分：海洋水文观测

HY/T147.7海洋监测技术规程第7部分：卫星遥感技术方法

JTS/T231水运工程模拟试验技术规范

3术语和定义

下列术语和定义适用于本文件。

3.1

湾口baymouth

海湾两个岬角之间的出口。

3.2

纳潮量tidalprism

由于潮汐涨落引起的海湾水体体积的变化量。

3.3

潮差tidalrange

在一个潮汐周期内，相邻的高潮与低潮之潮位差。

[来源：GB/T39619—2020，6.2.5]

3.4

单宽流量flowdischargeperunitwidth

单位时间通过单位宽度断面的水体体积。

4基本程序

纳潮量计算宜按下列程序进行：

a)收集模型建立与验证资料；

b)进行海湾二维潮波数值模拟与验证；

c)基于二维潮波数值模拟结果计算海湾纳潮量；

2

d)按照成果与报告要求编写海湾纳潮量数值计算报告。

5资料收集

5.1岸线与水深

岸线资料来源于遥感影像的解译，遥感影像的空间分辨率、成像时间和图像质量等要求按照HY/T147.7的规定执行。水深资料来源于数字化海图或现场实测水深资料，水深测量年份应与岸线资料年份相近。

5.2海流与水位

海流与水位资料来源于海上现场观测，海流、水位观测原则上同步进行。海流与水位观测方法和资料处理按照GB/T12763.2的规定执行。

6二维潮波数值模型

6.1控制方程与数值算法

模型的控制方程按照JTS/T231的规定执行。方程的离散求解方法可采用有限体积法、有限差分法或有限元法等。

6.2计算域的确定及网格剖分

模型计算域应足够大，并使开边界附近的水文要素不受海湾岸线、水深变化等的影响。海湾内部网格应加密，在疏密网格镶嵌界线上的加密度不超过3倍，并能够概化海湾岸线边界，网格结点水深应能够反映水下地形特征。

6.3边界条件

模型的边界条件应满足下列要求：

a)固边界条件：垂直于陆地边界的法向速度为零；

b)动边界条件：潮间带宜采用“干、湿”滑动边界；

c)开边界条件：开边界水位采用潮位或潮汐调和常数形式给出，在采用潮汐调和常数形式给出时，至少包括K1、O1、M2、S2等分潮。

6.4基本参数确定

曼宁糙率系数由实测潮位、流速和流向资料确定，应考虑水深、海床组成变化等因素引起的糙率平面分布不同。水平涡流粘度宜采用Smagorinsky公式计算。

6.5验证计算及精度控制

模型应不少于2个站位的潮位、3个站位潮流观测资料验证。模型应选取与岸线数据年份对应的观测资料进行验证。潮位验证时长不少于1个月，潮流验证时长不少于25h。模型验证精度控制按照JTS/T231的规定执行。

7纳潮量计算

3

7.1平均纳潮量计算

对潮波计算结果采用最小二乘法进行调和分析，得到潮汐调和常数，半日潮海湾平均潮差可依据潮汐调和常数进行计算。模型第i个网格的平均潮差ΔH计算见公式（1）：

?H=2.02H2+0.58+0.08·································(1)

式中：

ΔH——第i个网格的平均潮差，单位为米（m）；

H2——第i个网格半日分潮M2的振幅，单位为米（m）；

H2——第i个网格半日分潮S2的振幅，单位为米（m）；

H1——第i个网格日分潮K1的振幅，单位为米（m）；

振幅，单位为米（m）。

=

中：平均纳潮量，单位为立方米（m3）；

ΔH——第i个网格的平均潮差，单位为米（m）；

si——第i个网格的面积，单位为平方米（m2）；

n——计算域网格数。

7.2潮周期纳潮量计算

在湾口断面上布设M个控制点，相邻两个控制点间距为Δli，涨潮或落潮期设为N个时间步，时间步长为Δt。从潮波计算结果中提取断面上M个控制点N个时刻的水深、流速和流向，则流经断面第i个控制

点j时刻的单宽流量q计算见公式（3）：q=??v?