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

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)基于二维潮波数值模拟结果计算海湾纳潮量；2d)按照成果与报告要求编写海湾纳潮量数值计算报告。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纳潮量计算37.1平均纳潮量计算对潮波计算结果采用最小二乘法进行调和分析，得到潮汐调和常数，半日潮海湾平均潮差可依据潮汐调和常数进行计算。模型第i个网格的平均潮差ΔH计算见公式（1）：式中：ΔH——第i个网格的平均潮差，单位为米（m）；H2——第i个网格半日分潮M2的振幅，单位为米（m）；H2——第i个网格半日分潮S2的振幅，单位为米（m）；H1——第i个网格日分潮K1的振幅，单位为米（m）；振幅，单位为米（m）。中：平均纳潮量，单位为立方米（m3ΔH——第i个网格的平均潮差，单位为米（m）；si——第i个网格的面积，单位为平方米（m2n——计算域网格数。7.2潮周期纳潮量计算在湾口断面上布设M个控制点，相邻两个控制点间距为Δli，涨潮或落潮期设为N个时间步，时间步长为Δt。从潮波计算结果中提取断面上M个控制点N个时刻的水深、流速和流向，则流经断面第i个控制点j时刻的单宽流量q计算见公式（3 点j时刻的单宽流量q计算见公式（3式中：q——第i个控制点j时刻的单宽流量，单位为平方米每秒（m2/s）；ℎ——第i个控制点j时刻的水深，单位为米（m）；v——第i个控制点j时刻的流速，单位为米每秒（m/s）；a——第i个控制点j时刻的流向，单位为弧度（rad）；Y——N向与断面起点至终点方向顺时针的夹角，单位为弧度（rad），面向海湾，断面左端点为起点，右端点为终点，见图1。4图1Y角定义示意图j时刻湾口断面流量ΔQj计算见公式（4）：式中：ΔQj——j时刻湾口断面流量，单位为立方米每秒（m3/s）；M——湾口断面控制点数；q——第i个控制点j时刻的单宽流量，单位为平方米每秒（m2/s）；q+1——第i+1个控制点j时刻的单宽流量，单位为平方米每秒（m2/s）；Δli——相邻两个控制点间距，单位为米（m）。潮周期纳潮量Q计算见公式（5）：式中：Q——潮周期纳潮量，单位为立方米（m3N——涨潮或落潮期时间步数；ΔQj——j时刻湾口断面流量，单位为立方米每秒（m3/s）；ΔQj+1——j+1时刻湾口断面流量，单位为立方米每秒（m3/s）；Δt——时间步长，单位为秒（s）。8成果与报告8.1成果要求计算成果应包括潮周期纳潮量和平均纳潮量，潮周期纳潮量至少包括大、小潮期纳潮量，计算时长不少于60个连续潮周期，计算结果保留3位小数，计量单位采用立方千米（km3）。海湾纳潮量数值计算结果汇总表参照附录A。8.2报告编写成果报告主要内容包括前言、海湾概况、纳潮量计算、结论与建议等内容。海湾纳潮量数值计算报告格式和内容参照附录B。5（资料性）海湾纳潮量数值计算结果汇总表表A.1给出了海湾纳潮量数值计算结果的内容与格式。表A.1海湾纳潮量数值计算结果汇总表6（资料性）海湾纳潮量数值计算报告格式与内容B.1文本格式B.1.1文本规格海湾纳潮量数值计算报告文本外形尺寸为A4（210mm×297mm）。B.1.2封面格式海湾纳潮量数值计算报告封面格式如下：第一行书写：×××海湾纳潮量（居中）；第二行书写：数值计算报告（居中）；落款书写：编制单位全称（居中）；第四行书写：××××年××月（居中）；以上内容字体字号应适宜，各行间距应适中，保持封面美观。B.1.3封里1内容封里1应写明报告编制单位的全称、通讯地址、邮政编码、联系电话、传真、E-mail等。B.1.4封里2内容应分行写明：报告编制单位全称（加盖公章）、报告技术负责人（姓名、职务、职称）、项目负责人（姓名、职务、职称）、主要编制人员（姓名、职务、职称）、审核人（姓名、职务、职称）等，并签名。B.2报告章节内容报告应包括以下全部或部分内容。根据各海湾特点和纳潮量计算具体内容，可对下列章节及内容适当增减。B.2.1前言B.2.1.1目的意义B.2.1.2计算内容、范围B.2.1.3总体技术路线B.2.2海湾概况B.2.2.1地理位置及区位条件B.2.2.2自然环境条件B.2.2.3开发利用现状B.2.3