CN107256312A - 一种基于潮流环境下海湾纳潮变化量计算方法 - Google Patents
一种基于潮流环境下海湾纳潮变化量计算方法 Download PDFInfo
- Publication number
- CN107256312A CN107256312A CN201710443352.8A CN201710443352A CN107256312A CN 107256312 A CN107256312 A CN 107256312A CN 201710443352 A CN201710443352 A CN 201710443352A CN 107256312 A CN107256312 A CN 107256312A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mfrac
- mover
- bay
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01F—MEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
- G01F22/00—Methods or apparatus for measuring volume of fluids or fluent solid material, not otherwise provided for
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Fluid Mechanics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Other Liquid Machine Or Engine Such As Wave Power Use (AREA)
Abstract
本发明公开了一种基于潮流环境下海湾纳潮变化量计算方法，涉及海湾围垦工程。本发明其包括海湾围垦工程实际占用的海水体积的变化量；再以海湾围垦工程实施后造成的海湾水体潮波变形所导致的海水体积变化量来修正海湾纳潮变化量的误差。本发明给出基于数值模拟计算方法的新方法，为数值模拟计算方法与传统计算方法的改进方法。它真实有效的揭示了围垦工程实施造成的纳潮变化量。在实践中并证明了相关理论，为整个海湾进行围垦工程的设计与研究奠定了基础。
Description
技术领域
本发明涉及海湾围垦工程，具体涉及一种基于潮流环境下海湾纳潮变化量计算方法。
背景技术
近年来，随着我国改革开放的不断推进，经济建设的热潮一浪高过一浪，工程建设行业也得到了前所未有的大发展。沿海地区，尤其是沿海港口城市，土地资源日超紧缺，海涂围垦工程蓬勃发展。然而在河口地区，尤其是海湾内进行围垦会造成海湾纳潮量的改变。
纳潮量是海湾可以接纳的潮水体积，是一个潮周期中的最大水交换量，为海湾高潮水量与低潮水量之差。其大小直接影响着海湾水深的维持以及海湾与外海的水体交换和污染物的迁移扩散，制约着海湾的自净能力和环境容量，对维护海湾良好的生态环境至关重要。因此，研究围垦工程实施后，海湾内纳潮量的变化至关重要。
现有纳潮量的计算方法主要有两种：
(1)传统计算方法。即纳潮量计算可通过海湾高潮时的最大纳水体积与港湾低潮时的最小纳水体积之差求出，公式为：
W＝0.5×(S1+S2)(h1-h2)
式中：W为纳潮量；h1、h2分别为高、低潮位；S1、S2分别为高、低潮潮位对应的水域面积。
(2)ADCP测量方法。该方法通过在半封闭海湾的湾口设计一个封闭断面直接走航测量得到进出该断面的水量(纳潮量)。
(3)数值模拟方法。通过建立并验证工程海域潮流数学模型，计算出封闭湾口断面的潮位及流速流向，得到通过湾口断面的流量变化序列，再将这些流量值积分可得到海湾的纳潮量。
对于传统计算方法，该方法将高低潮位之间水体默认为坡度均匀一致的斜面，存在一定的局限性(如图1所示)。同时，在计算围垦工程实施后海湾纳潮变化量时，无法考虑围垦工程造成的潮波变形情况，且当围垦工程顶高程低于高潮位或围垦工程位置水深高于低潮位时则无法计算围垦工程对纳潮量变化的影响。
对于ADCP测量方法。该方法需要进行现场实测工作，除费用高之外，因海湾繁忙的船行以及较差的海况经常难于实施，即便实施测得的结果也难以做到同步而失真，需要做一定的数据修正工作。同时，ADCP测量方法只能测量现状条件下海湾纳潮量，并不能预测围垦工程实施后造成的海湾纳潮量变化。
对于数值模拟计算方法，其计算误差直接取决于模型研究人员的水平。
发明内容
本发明所要解决的技术问题是提供一种基于潮流环境下海湾纳潮变化量的计算方法，旨在给出基于数值模拟计算方法的新方法，为数值模拟计算方法与传统计算方法的改进方法。它真实有效的揭示了围垦工程实施造成的纳潮变化量。
为解决上述技术问题，本发明所采取的技术方案是：一种基于潮流环境下海湾纳潮变化量计算方法，其特征在于：其包括海湾围垦工程实际占用的海水体积的变化量；再以海湾围垦工程实施后造成的海湾水体潮波变形所导致的海水体积变化量来修正海湾纳潮变化量的误差。
进一步的技术方案在于，对海湾围垦工程实施后造成的海湾水体潮波变形所导致的水体体积变化量，是基于建立海湾潮流运动数学模型进行数值模拟，计算出封闭湾口断面以内水体高低潮位变化，进而得到水体体积变化量。
进一步的技术方案在于，所述海湾潮流运动数学模型为：
控制方程如下：
其中：h＝η+d，η和d分别表示水面高度和静水深，x和y分别表示横轴和纵轴坐标，t为时间，g为重力加速度，和分别为沿x和y方向的深度平均流速，f为柯氏力系数，ρ为流体密度，ρ0为参考密度，S为点源流量，us与vs为点源流速，Tij为应力项，包括粘性应力、紊流应力和对流等，根据水深平均的流速梯度计算；
底部应力由下式计算：
其中：cf是拖曳力系数，是水深平均的流速；拖曳力系数可以根据Chezy系数C或Manning系数M计算，
Manning系数可以根据底部糙率计算；
风应力计算公式为：
其中：ρa是空气密度，cd是空气拖曳力系数，是海面上10m高处的风速。
进一步的技术方案在于，所述海湾纳潮变化量为：
ΔW＝S工(h工1-h工2)+S河1×Δh河1+S河2×Δh河2
式中：ΔW为海湾纳潮变化量；S工为围填海工程面积；h工1为围填海工程附近高潮位(当高潮位低于围垦工程堤顶高程)或围垦工程堤顶高程(当高潮位高于围垦工程堤顶高程)；h工2为围填海工程附近低潮位(当低潮位高于围垦工程附近平均水深)或围垦工程附近水深(当低潮位低于围垦工程附近平均水深)；S河1、S河2分别为高、低潮潮位对应的河口水域面积(除去围垦工程面积)；Δh河1、Δh河2分别数值模拟计算出的高、低潮潮位变化值(海湾平均)；Δh河1增加为正值，减小为负值；Δh河2增加为负值，减小为正值。
采用上述技术方案所产生的有益效果在于：本发明给出基于数值模拟计算方法的新方法，为数值模拟计算方法与传统计算方法的改进方法。它真实有效的揭示了围垦工程实施造成的纳潮变化量。在实践中并证明了相关理论，为整个海湾进行围垦工程的设计与研究奠定了基础。
附图说明
下面结合附图和具体实施方式对本发明作进一步详细的说明。
图1是传统计算方法断面示意图；
图2是围垦工程前示意图(图中灰色为海域，白色为陆域或岛屿)；
图3是围垦工程后示意图(图中灰色为海域，白色为陆域或岛屿)；
图4是围垦工程占用海湾纳潮水体体积示意图(图中灰色水体)；
图5是实测潮位及流速流向测站示意图；
图6是网格示意图；
图7是潮位验证图；
图8是流速流向验证图；
图9是表1的特征点布置图。
具体实施方式
下面结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明的一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
在下面的描述中阐述了很多具体细节以便于充分理解本发明，但是本发明还可以采用其他不同于在此描述的其它方式来实施，本领域技术人员可以在不违背本发明内涵的情况下做类似推广，因此本发明不受下面公开的具体实施例的限制。
本发明提供了一种基于潮流环境下海湾纳潮变化量计算方法，其特征在于：其包括海湾围垦工程实际占用的海水体积的变化量；再以海湾围垦工程实施后造成的海湾水体潮波变形所导致的海水体积变化量来修正海湾纳潮变化量的误差。
优选的，对海湾围垦工程实施后造成的海湾水体潮波变形所导致的水体体积变化量，是基于建立海湾潮流运动数学模型进行数值模拟，计算出封闭湾口断面以内水体高低潮位变化，进而得到水体体积变化量。
优选的，所述海湾潮流运动数学模型为：
控制方程如下：
其中：h＝η+d，η和d分别表示水面高度和静水深，x和y分别表示横轴和纵轴坐标，t为时间，g为重力加速度，和分别为沿x和y方向的深度平均流速，f为柯氏力系数，ρ为流体密度，ρ0为参考密度，S为点源流量，us与vs为点源流速，Tij为应力项，包括粘性应力、紊流应力和对流等，根据水深平均的流速梯度计算；
底部应力由下式计算：
其中：cf是拖曳力系数，是水深平均的流速；拖曳力系数可以根据Chezy系数C或Manning系数M计算，
Manning系数可以根据底部糙率计算；
风应力计算公式为：
其中：ρa是空气密度，cd是空气拖曳力系数，是海面上10m高处的风速。
优选的，所述海湾纳潮变化量为：
ΔW＝S工(h工1-h工2)+S河1×Δh河1+S河2×Δh河2
式中：ΔW为海湾纳潮变化量；S工为围填海工程面积；h工1为围填海工程附近高潮位(当高潮位低于围垦工程堤顶高程)或围垦工程堤顶高程(当高潮位高于围垦工程堤顶高程)；h工2为围填海工程附近低潮位(当低潮位高于围垦工程附近平均水深)或围垦工程附近水深(当低潮位低于围垦工程附近平均水深)；S河1、S河2分别为高、低潮潮位对应的河口水域面积(除去围垦工程面积)；Δh河1、Δh河2分别数值模拟计算出的高、低潮潮位变化值(海湾平均)；Δh河1增加为正值，减小为负值；Δh河2增加为负值，减小为正值。
实施例：
一、纳潮变化量分解
纳潮量为海湾高潮时最大纳水体积与海湾低潮时最小纳水体积之差。以某海湾内围垦工程为例(图2和图3)，整个海湾面积可分为两部分，即：
S＝S工+S河
其纳潮量也可分为两部分：
W＝W工+W河
则纳潮变化量也可分为两部分：
ΔW＝ΔW工+ΔW河
二、围垦工程占用海湾水体体积
根据上述公式可见，整个海湾进行围垦工程后，海湾纳潮变化量包括两部分，其中一部分为围垦工程实际占用的海水体积，一部分为围垦工程实施后造成的海湾水体潮波变形所导致的水体体积变化量。其中，ΔW工可采用实测资料进行计算，如图4所示。
ΔW工＝S工(h工1-h工2)
其中：S工为围填海工程面积(已知)；h工1为围填海工程附近高潮位(当高潮位低于围垦工程堤顶高程时)，或围垦工程堤顶高程(当高潮位高于围垦工程堤顶高程时)；h工2为围填海工程附近低潮位(当低潮位高于围垦工程附近平均水深时)或围垦工程附近水深(当低潮位低于围垦工程附近平均水深时)。
其中，围垦工程附近水深可根据实测水深得到。高、低潮位可根据工程附近潮位站实测潮位资料得到，对应到本工程即为工程附近T4潮位站实测潮位资料(如图5所示)。
本工程围填海面积约为2.95km2，实测大潮两次涨潮过程潮差分别为5.20m和5.66m(工程附近T4站)。则围垦工程所占的潮棱体约0.32×108m3。
三、潮波变形造成纳潮变化量
ΔW河为围垦工程实施后海湾水域潮波变形导致的纳潮变化量，可采用数值模拟方法进行计算。可采用自主开发程序或比较成熟的MIKE21、FVCOM等软件计算。
3.1潮流运动数学模型
控制方程如下：
其中：h＝η+d，η和d分别表示水面高度和静水深，x和y分别表示横轴和纵轴坐标，t为时间，g为重力加速度，和分别为沿x和y方向的深度平均流速，f为柯氏力系数，ρ为流体密度，ρ0为参考密度，S为点源流量，us与vs为点源流速，Tij为应力项，包括粘性应力、紊流应力和对流等，根据水深平均的流速梯度计算。
底部应力由下式计算：
其中：cf是拖曳力系数，是水深平均的流速。拖曳力系数可以根据Chezy系数C或Manning系数M计算，
Manning系数可以根据底部糙率计算。
风应力计算公式为：
其中：ρa是空气密度，cd是空气拖曳力系数，是海面上10m高处的风速。
3.2计算域选取及网格剖分
以上述某海区为例进行计算：计算区域西边界取至上游水闸，东西方向长约97km，北方向长约62km。模型模拟区域总面积达6000km2。
如图6所示，采用三角形网格进行剖分，距离工程区较远的海区网格较稀疏，对工程海域所在的海湾进行网格加密。计算域采用的三角形网格边长最大为2000m，最小网格边长为10m，总共有13669个网格节点，25028个三角形单元。
3.3边界条件及概化方法
外海开边界采用潮位过程控制，由中国海域潮汐模型提供。上游主要径流根据水文站实测流量给出。模型中可考虑漫滩、潜堤等动边界的处理。
3.4模型验证
为验证模型的合理性，采用实测水文全潮资料对模型进行率定和验证(测站参见图5)。部分潮位及流速流向验证曲线参见图7和图8，精度满足规程要求。
3.5潮位变化统计
在工程所在海湾内均匀地选取尽量多的特征点(图9)。表1给出了数值模拟计算出来的围填海工程实施前后特征点位置高低潮位变化情况
表1特征点位置高低潮位变化值(m)
3.6潮波变形造成纳潮变化量
潮波变形造成的纳潮变化量为高潮水体变化量与低潮水体变化量之和：
ΔW河＝S河1×Δh河1+S河2×Δh河2
S河1、S河2分别为高、低潮潮位对应的河口水域面积(除去围垦工程)，其可根据海湾实测水深及工程附近实测高低潮位统计得到。
Δh河1、Δh河2分别为数值模拟计算出的高、低潮潮位变化值(海湾平均)。
其中：Δh河1增加为正值，减小为负值；Δh河2增加为负值，减小为正值。
本文实例计算可得，由于潮波变形造成的纳潮量变化约为0.043×108m3。
四、海湾纳潮变化量计算方法总结
针对围垦工程实施后海湾纳潮变化量，作者提出了一种新的计算方法，该方法基于潮流数值模拟，公式如下,为上述两部分之和：
ΔW＝ΔW工+ΔW河
ΔW工＝S工(h工1-h工2)
ΔW河＝S河1×Δh河1+S河2×Δh河2
ΔW＝S工(h工1-h工2)+S河1×Δh河1+S河2×Δh河2
式中：ΔW为由于围填海工程造成的海湾纳潮变化量；S工为围填海工程面积；h工1为围填海工程附近高潮位(当高潮位低于围垦工程堤顶高程)或围垦工程堤顶高程(当高潮位高于围垦工程堤顶高程)；h工2为围填海工程附近低潮位(当低潮位高于围垦工程附近平均水深)或围垦工程附近水深(当低潮位低于围垦工程附近平均水深)；S河1、S河2分别为高、低潮潮位对应的河口水域面积(除去围垦工程)。Δh河1、Δh河2分别数值模拟计算出的高、低潮潮位变化值(海湾平均)。Δh河1增加为正值，减小为负值；Δh河2增加为负值，减小为正值。
综上所述，本为实例中围填海工程实施后纳潮变化量为0.363×108m3。
Claims (4)
1.一种基于潮流环境下海湾纳潮变化量计算方法，其特征在于：其包括海湾围垦工程实际占用的海水体积的变化量；再以海湾围垦工程实施后造成的海湾水体潮波变形所导致的海水体积变化量来修正海湾纳潮变化量的误差。
2.根据权利要求1所述的一种基于潮流环境下海湾纳潮变化量计算方法，其特征在于：对海湾围垦工程实施后造成的海湾水体潮波变形所导致的水体体积变化量，是基于建立海湾潮流运动数学模型进行数值模拟，计算出封闭湾口断面以内水体高低潮位变化，进而得到水体体积变化量。
3.根据权利要求2所述的一种基于潮流环境下海湾纳潮变化量计算方法，其特征在于：所述海湾潮流运动数学模型为：
控制方程如下：
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>h</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>h</mi>
<mover>
<mi>u</mi>
<mo>&OverBar;</mo>
</mover>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>h</mi>
<mover>
<mi>v</mi>
<mo>&OverBar;</mo>
</mover>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>y</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mi>h</mi>
<mi>S</mi>
</mrow>
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>h</mi>
<mover>
<mi>u</mi>
<mo>&OverBar;</mo>
</mover>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>h</mi>
<msup>
<mover>
<mi>u</mi>
<mo>&OverBar;</mo>
</mover>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>h</mi>
<mover>
<mrow>
<mi>v</mi>
<mi>u</mi>
</mrow>
<mo>&OverBar;</mo>
</mover>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>y</mi>
</mrow>
</mfrac>
<mo>=</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>f</mi>
<mover>
<mi>v</mi>
<mo>&OverBar;</mo>
</mover>
<mi>h</mi>
<mo>-</mo>
<mi>g</mi>
<mi>h</mi>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&eta;</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mfrac>
<mrow>
<msup>
<mi>gh</mi>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<mn>2</mn>
<msub>
<mi>&rho;</mi>
<mn>0</mn>
</msub>
</mrow>
</mfrac>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&rho;</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<msub>
<mi>&tau;</mi>
<mrow>
<mi>s</mi>
<mi>x</mi>
</mrow>
</msub>
<msub>
<mi>&rho;</mi>
<mn>0</mn>
</msub>
</mfrac>
<mo>-</mo>
<mfrac>
<msub>
<mi>&tau;</mi>
<mrow>
<mi>b</mi>
<mi>x</mi>
</mrow>
</msub>
<msub>
<mi>&rho;</mi>
<mn>0</mn>
</msub>
</mfrac>
<mo>+</mo>
<mfrac>
<mo>&part;</mo>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>hT</mi>
<mrow>
<mi>x</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mfrac>
<mo>&part;</mo>
<mrow>
<mo>&part;</mo>
<mi>y</mi>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>hT</mi>
<mrow>
<mi>x</mi>
<mi>y</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>hu</mi>
<mi>s</mi>
</msub>
<mi>S</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>h</mi>
<mover>
<mi>v</mi>
<mo>&OverBar;</mo>
</mover>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>h</mi>
<mover>
<mi>v</mi>
<mo>&OverBar;</mo>
</mover>
<mover>
<mi>u</mi>
<mo>&OverBar;</mo>
</mover>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>h</mi>
<msup>
<mover>
<mi>v</mi>
<mo>&OverBar;</mo>
</mover>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>y</mi>
</mrow>
</mfrac>
<mo>=</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mi>f</mi>
<mover>
<mi>u</mi>
<mo>&OverBar;</mo>
</mover>
<mi>h</mi>
<mo>-</mo>
<mi>g</mi>
<mi>h</mi>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&eta;</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>y</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mfrac>
<mrow>
<msup>
<mi>gh</mi>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<mn>2</mn>
<msub>
<mi>&rho;</mi>
<mn>0</mn>
</msub>
</mrow>
</mfrac>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&rho;</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>y</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<msub>
<mi>&tau;</mi>
<mrow>
<mi>s</mi>
<mi>y</mi>
</mrow>
</msub>
<msub>
<mi>&rho;</mi>
<mn>0</mn>
</msub>
</mfrac>
<mo>-</mo>
<mfrac>
<msub>
<mi>&tau;</mi>
<mrow>
<mi>b</mi>
<mi>y</mi>
</mrow>
</msub>
<msub>
<mi>&rho;</mi>
<mn>0</mn>
</msub>
</mfrac>
<mo>+</mo>
<mfrac>
<mo>&part;</mo>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>hT</mi>
<mrow>
<mi>x</mi>
<mi>y</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mfrac>
<mo>&part;</mo>
<mrow>
<mo>&part;</mo>
<mi>y</mi>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>hT</mi>
<mrow>
<mi>y</mi>
<mi>y</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>hv</mi>
<mi>s</mi>
</msub>
<mi>S</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
其中：h＝η+d，η和d分别表示水面高度和静水深，x和y分别表示横轴和纵轴坐标，t为时间，g为重力加速度，和分别为沿x和y方向的深度平均流速，f为柯氏力系数，ρ为流体密度，ρ0为参考密度，S为点源流量，us与vs为点源流速，Tij为应力项，包括粘性应力、紊流应力和对流等，根据水深平均的流速梯度计算；
底部应力由下式计算：
<mrow>
<mfrac>
<mover>
<msub>
<mi>&tau;</mi>
<mi>b</mi>
</msub>
<mo>&RightArrow;</mo>
</mover>
<msub>
<mi>&rho;</mi>
<mn>0</mn>
</msub>
</mfrac>
<mo>=</mo>
<msub>
<mi>c</mi>
<mi>f</mi>
</msub>
<mover>
<msub>
<mi>u</mi>
<mi>b</mi>
</msub>
<mo>&RightArrow;</mo>
</mover>
<mrow>
<mo>|</mo>
<mover>
<msub>
<mi>u</mi>
<mi>b</mi>
</msub>
<mo>&RightArrow;</mo>
</mover>
<mo>|</mo>
</mrow>
</mrow>
其中：cf是拖曳力系数，是水深平均的流速；拖曳力系数可以根据Chezy系数C或Manning系数M计算，
<mrow>
<msub>
<mi>c</mi>
<mi>f</mi>
</msub>
<mo>=</mo>
<mfrac>
<mi>g</mi>
<msup>
<mi>C</mi>
<mn>2</mn>
</msup>
</mfrac>
</mrow>
<mrow>
<msub>
<mi>c</mi>
<mi>f</mi>
</msub>
<mo>=</mo>
<mfrac>
<mi>g</mi>
<msup>
<mrow>
<mo>(</mo>
<msup>
<mi>Mh</mi>
<mrow>
<mn>1</mn>
<mo>/</mo>
<mn>6</mn>
</mrow>
</msup>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mfrac>
</mrow>
Manning系数可以根据底部糙率计算；
风应力计算公式为：
<mrow>
<msub>
<mi>&tau;</mi>
<mi>s</mi>
</msub>
<mo>=</mo>
<msub>
<mi>&rho;</mi>
<mi>a</mi>
</msub>
<msub>
<mi>c</mi>
<mi>d</mi>
</msub>
<mrow>
<mo>|</mo>
<msub>
<mi>u</mi>
<mi>w</mi>
</msub>
<mo>|</mo>
</mrow>
<mover>
<msub>
<mi>u</mi>
<mi>w</mi>
</msub>
<mo>&RightArrow;</mo>
</mover>
</mrow>
1
其中：ρa是空气密度，cd是空气拖曳力系数，是海面上10m高处的风速。
4.根据权利要求2一种基于潮流环境下海湾纳潮变化量计算方法，其特征在于：所述海湾纳潮变化量为：
ΔW＝S工(h工1-h工2)+S河1×Δh河1+S河2×Δh河2
式中：ΔW为海湾纳潮变化量；S工为围填海工程面积；h工1为围填海工程附近高潮位(当高潮位低于围垦工程堤顶高程)或围垦工程堤顶高程(当高潮位高于围垦工程堤顶高程)；h工2为围填海工程附近低潮位(当低潮位高于围垦工程附近平均水深)或围垦工程附近水深(当低潮位低于围垦工程附近平均水深)；S河1、S河2分别为高、低潮潮位对应的河口水域面积(除去围垦工程面积)；Δh河1、Δh河2分别数值模拟计算出的高、低潮潮位变化值(海湾平均)；Δh河1增加为正值，减小为负值；Δh河2增加为负值，减小为正值。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710443352.8A CN107256312B (zh) | 2017-06-13 | 2017-06-13 | 一种基于潮流环境下海湾纳潮变化量计算方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710443352.8A CN107256312B (zh) | 2017-06-13 | 2017-06-13 | 一种基于潮流环境下海湾纳潮变化量计算方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107256312A true CN107256312A (zh) | 2017-10-17 |
| CN107256312B CN107256312B (zh) | 2020-09-25 |
Family
ID=60024043
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710443352.8A Expired - Fee Related CN107256312B (zh) | 2017-06-13 | 2017-06-13 | 一种基于潮流环境下海湾纳潮变化量计算方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107256312B (zh) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111695250A (zh) * | 2020-06-04 | 2020-09-22 | 哈尔滨工程大学 | 一种内潮波特征提取方法 |
| CN108595831B (zh) * | 2018-04-22 | 2021-05-07 | 天津大学 | 一种计算实时潮致混合的动态模拟方法 |
| CN115100537A (zh) * | 2022-06-27 | 2022-09-23 | 中国科学院地理科学与资源研究所 | 一种基于遥感影像的潮汐能资源评估方法 |
| CN116796517A (zh) * | 2023-05-30 | 2023-09-22 | 山东省海洋资源与环境研究院(山东省海洋环境监测中心、山东省水产品质量检验中心) | 一种湾口断面纳潮量计算方法、系统、设备及介质 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103774605A (zh) * | 2014-01-10 | 2014-05-07 | 河海大学 | 一种提高环抱式港池水体交换能力的设计方法 |
| CN106168991A (zh) * | 2016-06-24 | 2016-11-30 | 珠江水利委员会珠江水利科学研究院 | 一种基于水动力数值模拟的感潮河网潮位预报方法 |
| CN106709262A (zh) * | 2017-01-11 | 2017-05-24 | 河海大学 | 基于动态概化的潮流能水轮机阵列水动力影响数值模拟方法 |
-
2017
- 2017-06-13 CN CN201710443352.8A patent/CN107256312B/zh not_active Expired - Fee Related
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103774605A (zh) * | 2014-01-10 | 2014-05-07 | 河海大学 | 一种提高环抱式港池水体交换能力的设计方法 |
| CN106168991A (zh) * | 2016-06-24 | 2016-11-30 | 珠江水利委员会珠江水利科学研究院 | 一种基于水动力数值模拟的感潮河网潮位预报方法 |
| CN106709262A (zh) * | 2017-01-11 | 2017-05-24 | 河海大学 | 基于动态概化的潮流能水轮机阵列水动力影响数值模拟方法 |
Non-Patent Citations (4)
| Title |
|---|
| YU-HAI WANG: "Combined effects of channel dredging,land reclamation and long-range jetties upon the long-term evolution of channel-shoal system in Qinzhou bay,SW China", 《OCEAN ENGINEERING》 * |
| 孔俊 等: "风暴潮模拟中潮位对风拖曳力系数的影响研究", 《海洋预报》 * |
| 张坤: "钦州湾茅尾海纳潮量和水交换能力的数值模拟研究", 《中国优秀硕士学位论文全文数据库 基础科学辑》 * |
| 杨世伦 等: "半封闭海湾潮间带部分围垦后纳潮量计算的商榷——以胶州湾为例", 《海洋科学》 * |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108595831B (zh) * | 2018-04-22 | 2021-05-07 | 天津大学 | 一种计算实时潮致混合的动态模拟方法 |
| CN111695250A (zh) * | 2020-06-04 | 2020-09-22 | 哈尔滨工程大学 | 一种内潮波特征提取方法 |
| CN111695250B (zh) * | 2020-06-04 | 2022-12-13 | 哈尔滨工程大学 | 一种内潮波特征提取方法 |
| CN115100537A (zh) * | 2022-06-27 | 2022-09-23 | 中国科学院地理科学与资源研究所 | 一种基于遥感影像的潮汐能资源评估方法 |
| CN116796517A (zh) * | 2023-05-30 | 2023-09-22 | 山东省海洋资源与环境研究院(山东省海洋环境监测中心、山东省水产品质量检验中心) | 一种湾口断面纳潮量计算方法、系统、设备及介质 |
| CN116796517B (zh) * | 2023-05-30 | 2024-02-02 | 山东省海洋资源与环境研究院(山东省海洋环境监测中心、山东省水产品质量检验中心) | 一种湾口断面纳潮量计算方法、系统、设备及介质 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107256312B (zh) | 2020-09-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111950211B (zh) | 基于海流设计流速的海床基础局部冲刷深度确定方法及系统 | |
| CN107256312A (zh) | 一种基于潮流环境下海湾纳潮变化量计算方法 | |
| CN107895059A (zh) | 一种淤泥质海岸高浊度海区岛礁促淤工程模拟方法 | |
| Li et al. | Responses of water environment to tidal flat reduction in Xiangshan Bay: Part I hydrodynamics | |
| CN103397615B (zh) | 一种环抱式港区建设期防淤减淤的设计方法 | |
| Kuang et al. | A two-dimensional morphological model based on next generation circulation solver II: Application to Caofeidian, Bohai Bay, China | |
| CN104091065A (zh) | 一种求解浅水问题模拟间断水流数值的方法 | |
| CN115034468B (zh) | 一种淤泥质海岸丁坝或防波堤拆除后海床冲淤变化预测方法 | |
| CN108625337B (zh) | 一种确定潮流界以下沙质河床河段整治水位的方法 | |
| CN114329950B (zh) | 基于动态概化的斜坡式潜堤波浪水动力影响数值模拟方法 | |
| CN115840975A (zh) | 一种风暴潮增水漫堤预警方法、系统、装置和存储介质 | |
| Daly et al. | Modelling accretion at nha trang beach, vietnam | |
| Hou et al. | The morphodynamic responses to deposition-promoting projects in island and reef coasts of the Zhoushan Archipelago, China | |
| Deng et al. | Wind influence on the orientation of estuarine landforms: An example from Lake Illawarra in southeastern Australia | |
| CN111611641B (zh) | 一种考虑余流的港口建设方法 | |
| Uzaki et al. | Numerical and field study of sediment budgets on an intertidal flat at the mouth of the Shirakawa River | |
| Bao et al. | Tidal characteristics in the Wenzhou offshore waters and changes resulting from the Wenzhou Shoal Reclamation Project | |
| Zhou et al. | Cross‐Shelf Penetrating Fronts of Buoyant Coastal Currents Around the Headland | |
| Yang et al. | Potential causes of salinity trend in seawater overlaying a mudflat inside a semi-enclosed bay | |
| CN112648960A (zh) | 一种计算湖岸远岸坝与近岸坝间距的方法 | |
| Xie et al. | Promoting siltation effects and impacts of Hengsha East Shoal on the Yangtze River Estuary | |
| Hu et al. | Sediment transport in the nearshore area of Phoenix Island | |
| Ding et al. | Three-dimensional numerical simulation of tidal flows in the Yangtze River Estuary | |
| Zhang et al. | A numerical study on storm surge and sediment Resuspending in Modaomen Estuary during Typhoon Hagupit | |
| Jiang et al. | Waterway design criteria of land reclamation based on an environmental fluid dynamics computer model |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200925 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |