一种确定入海河口生态植被水流阻力的方法与流程

1.本发明属于水利工程领域，具体涉及一种确定入海河口生态植被水流阻力及阻力系数的方法。


背景技术：

2.在水体流经有生态植被分布区域时，水体中的微团相互碰撞、混掺，此时水体流动的主要阻力是由植被产生的阻力。有生态植被分布的入海河口水域，由植被产生的水流阻力往往是主要阻力，这种由生态植被带来的水流阻力，会减缓水流流速，滞留水体中的泥沙和形成促淤等，对入海河口的沙体沙洲发育和生态植被的生长发展有明显的影响。因此，生态植被水流阻力的确定对入海河口的湿地规划、管理及影响研究具有重要意义。
3.目前现有水流阻力计算方程需要给定阻力系数，一维水流条件下已有多种确定阻力系数的方法，而二维非恒定流条件下的阻力系数取值都是基于经验选取，尚没有基于理论和观测试验来确定二维非恒定流条件下阻力系数的方法。
4.申请号为：201410193863.5的发明申请，公开了“一种含漂浮植物河道的模拟系统、水流阻力计算方法及装置”，其中，含漂浮植物河道的水流阻力计算方法包括：建立用于模拟河道的循环可变坡水槽系统，并在所述循环可变坡水槽系统中建立用于模拟漂浮植物的漂浮植物固定系统；调节所述循环可变坡水槽系统及漂浮植物固定系统，使其物理量参数满足设计要求，所述物理量参数包括水流量、水位及没入水中的植物长度；测量所述循环可变坡水槽系统中的水槽底坡及流速统计变量；根据所述测量的水槽底坡及流速统计变量，通过预设的计算模型，计算得到漂浮植物水流阻力。
5.申请号为：201510524660.4的发明申请，公开了“一种计算不同状态挺水植被水流阻力的方法”，主要包括四个步骤：计算单株挺水植被的阻力系数，计算单株挺水植被的阻力，计算单断面挺水植被群的水流阻力，计算多断面挺水植被群落的水流阻力。不仅能够直接计算挺水植被保持直立状态时的水流阻力，而且能够计算挺水植被处于摇摆状态或倒伏状态下的水流阻力。
6.申请号为：201610900304.2的发明申请，公开了“一种确定城市生态河道水流阻力的模拟系统及水流阻力计算方法”，所述的模拟系统包括城市生态河道模拟系统和河道水流阻力测量系统，所述城市生态河道模拟系统包括模拟城市生态河道和模拟植被带，所述模拟植被带置于所述模拟城市生态河道中。通过建立城市生态河道模拟系统和河道水流阻力测量系统，结合水力学阻力理论知识得出不同水力条件下的水流阻力，并利用力矩平衡原理计算水流阻力，建立的模拟系统简单实用，填补了城市生态河道水流阻力研究的空白，为城市生态河道规划设计提供技术支持。


技术实现要素：

7.本发明提供了一种确定入海河口生态植被水流阻力的方法，其技术方案具体如下：
8.一种确定入海河口生态植被水流阻力的方法，其特征在于包括如下步骤：
9.s1：基于阻力方程，并结合守恒型变形微分动量方程，反算出实时水流情况对应的阻力系数；
10.s2：据此阻力系数，完成相应的水流阻力确定。
11.根据本发明的一种确定入海河口生态植被水流阻力的方法，其特征在于：
12.步骤s1中的阻力方程为一维恒定流或一维非恒定流阻力方程；反算出一维恒定流或一维非恒定流阻力系数。
13.根据本发明的一种确定入海河口生态植被水流阻力的方法，其特征在于：
14.步骤s1中的阻力方程为二维恒定流或二维非恒定流阻力方程；反算出二维恒定流或二维非恒定流阻力系数。
15.根据本发明的一种确定入海河口生态植被水流阻力的方法，其特征在于：
16.步骤s1具体包括如下步骤：
17.s11：根据阻力方程形成阻力系数表达；
18.s12：离散化守恒型变形微分动量方程；
19.s13：通过联立阻力系数表达与离散化的守恒型变形微分动量方程，并结合测量确定的数值，形成时间序列上的非恒定流阻力系数确定。
20.根据本发明的一种确定入海河口生态植被水流阻力的方法，其特征在于：
21.所述的守恒型变形微分动量方程为：
[0022][0023]
其中，
[0024]
u：x方向的流速，单位：m/s；
[0025]
v：y方向的流速，单位：m/s；
[0026]
ρ：水密度，单位：kg/m3；
[0027]
g：重力加速度，单位：m/s2；
[0028]
t：计算时间，单位：s；
[0029]
η：水位，单位：m；
[0030]
h：总水深，单位：m；
[0031]
τ
x
：x方向的阻力，单位：n；
[0032]
τ
y
：y方向的阻力，单位：n。
[0033]
根据本发明的一种确定入海河口生态植被水流阻力的方法，其特征在于：
[0034]
所述的阻力方程为：
[0035]
[0036]
其中，
[0037][0038]
u：x轴向的平均流速分量，单位：m/s；
[0039]
v：y轴向的平均流速分量，单位：m/s；
[0040]
ρ：水密度，单位：kg/m3；
[0041]
g：重力加速度，单位：m/s2；
[0042]
n：阻力系数，为无量纲数；
[0043]
h：总水深，单位：m；
[0044]
τ
x
：x方向的阻力，单位：n；
[0045]
τ
y
：y方向的阻力，单位：n。
[0046]
根据本发明的一种确定入海河口生态植被水流阻力的方法，其特征在于：
[0047]
步骤s12的离散化通过有限差分法与有限体积法综合完成。
[0048]
根据本发明的一种确定入海河口生态植被水流阻力的方法，其特征在于：
[0049]
步骤s13中的联立阻力系数表达与离散化的守恒型变形微分动量方程，具体如下：
[0050][0051]
其中，
[0052]
m：计算时间分步编号；
[0053]
δx
j
：控制体单元第j条边的x方向增量，单位：m；
[0054]
δy
j
：控制体单元第j条边的y方向增量，单位：m；
[0055]
k：控制体边个数；
[0056]
f：控制体单元的中心；
[0057]
：时间编号m对应时间的边中心的总水深，单位：m；
[0058]
：时间编号m+1对应时间的边中心的总水深，单位：m；
[0059]
：时间编号m对应时间的x轴向的边中心的平均流速分量，单位：m/s；
[0060]
：时间编号m+1对应时间的x轴向的边中心的平均流速分量，单位：m/s；
[0061]
：时间编号m对应时间的y轴向的边中心的平均流速分量，单位：m/s；
[0062]
：时间编号m+1对应时间的y轴向的边中心的平均流速分量，单位：m/s；
[0063]
：时间编号m对应时间的单元中心的总水深，单位：m；
[0064]
：时间编号m对应时间的x轴向的单元中心的平均流速分量，单位：m/s；
[0065]
：时间编号m对应时间的y轴向的单元中心的平均流速分量，单位：m/s；
[0066]
：时间编号m对应时间的单元中心的水位，单位：m；
[0067]
da：控制单元面积，单位：m2。
[0068]
根据本发明的一种确定入海河口生态植被水流阻力的方法，其特征在于：
[0069]
步骤s13中所述的测量确定，具体为：
[0070]
先组建与布放观测设备，而后根据观测设备实时采集相应数据、计算完成相应数值的测量确定。
[0071]
根据本发明的一种确定入海河口生态植被水流阻力的方法，其特征在于：
[0072]
所述观测设备包括有：
[0073]
用于测量近底高分辨率流速场的多普勒声学流速剖面仪及脉冲耦合模式多普勒声学流速剖面仪；
[0074]
用于测量近底高精度流速过程与紊动过程的多普勒点流速仪；
[0075]
用于获取水体中悬浮物浊度的光学后向散射浊度仪；
[0076]
用于测量波浪和潮高的浪潮仪。
[0077]
本发明的一种确定入海河口生态植被水流阻力的方法，通过阻力方程与守恒型变形微分动量方程的联立，并结合实际测量，反算出时间序列上的二维非恒定流阻力系数，并据此阻力系数完成阻力值的确定，使得基于此的阻力值更能贴合与表征实际情况下的二维非恒定流水流阻力(该技术方案同样适用于二维恒定流以及一维恒定或一维非恒定流)；针对现况中的入海河口区域的水动力受径潮流影响，具有明显的非恒定性及二维特征，且入海河口的生态植被具有复杂的生物特征，导致入海河口生态植被水流阻力通常难以精确给定的问题，提供了一种便于获取二维非恒定流阻力系数的现场测量方案以及基于现场观测的二维非恒定流阻力系数的计算与确定方法，通过此阻力系数，提高入海河口生态植被水流阻力的计算准确度，为确定入海河口生态植被水流阻力提供了新的便捷、高精度和科学的方法，对入海河口的湿地规划、管理及影响研究具有重要意义。
附图说明
[0078]
图1为本发明的步骤示意图；
[0079]
图2为本发明中的阻力系数确定步骤示意图；
[0080]
图3为本发明实施例中的生态植被水流阻力计算示意图；
[0081]
图4为本发明实施例中的控制单元结构及观测示意图。
具体实施方式
[0082]
下面，根据说明书附图和具体实施方式对本发明的一种确定入海河口生态植被水流阻力的方法作进一步具体说明。
[0083]
如图1所示的一种确定入海河口生态植被水流阻力的方法，其特征在于包括如下步骤：
[0084]
s1：基于阻力方程，并结合守恒型变形微分动量方程，反算出实时水流情况对应的阻力系数；
[0085]
s2：据此阻力系数，完成相应的水流阻力确定。
[0086]
其中，
[0087]
步骤s1中的阻力方程为一维恒定流或一维非恒定流阻力方程；反算出一维恒定流或一维非恒定流阻力系数。
[0088]
其中，
[0089]
步骤s1中的阻力方程为二维恒定流或二维非恒定流阻力方程；反算出二维恒定流或二维非恒定流阻力系数。
[0090]
其中，
[0091]
步骤s1具体包括如下步骤，如图2所示：
[0092]
s11：根据阻力方程形成阻力系数表达；
[0093]
s12：离散化守恒型变形微分动量方程；
[0094]
s13：通过联立阻力系数表达与离散化的守恒型变形微分动量方程，并结合测量确定的数值，形成时间序列上的非恒定流阻力系数确定。
[0095]
其中，
[0096]
所述的守恒型变形微分动量方程为：
[0097][0098]
其中，
[0099]
u：x方向的流速，单位：m/s；
[0100]
v：y方向的流速，单位：m/s；
[0101]
ρ：水密度，单位：kg/m3；
[0102]
g：重力加速度，单位：m/s2；t：计算时间，单位：s；
[0103]
η：水位，单位：m；
[0104]
h：总水深，单位：m；
[0105]
τ
x
：x方向的阻力，单位：n；
[0106]
τ
y
：y方向的阻力，单位：n。
[0107]
其中，
[0108]
所述的阻力方程为：
[0109][0110]
其中，
[0111][0112]
u：x轴向的平均流速分量，单位：m/s；
[0113]
v：y轴向的平均流速分量，单位：m/s；
[0114]
ρ：水密度，单位：kg/m3；
[0115]
g：重力加速度，单位：m/s2；
[0116]
n：阻力系数，为无量纲数；
[0117]
h：总水深，单位：m；
[0118]
τ
x
：x方向的阻力，单位：n；
[0119]
τ
y
：y方向的阻力，单位：n。
[0120]
其中，
[0121]
步骤s12的离散化通过有限差分法与有限体积法综合完成。
[0122]
其中，
[0123]
步骤s13中的联立阻力系数表达与离散化的守恒型变形微分动量方程，具体如下：
[0124][0125]
其中，
[0126]
m：计算时间分步编号；
[0127]
δx
j
：控制体单元第j条边的x方向增量，单位：m；
[0128]
δy
j
：控制体单元第j条边的y方向增量，单位：m；
[0129]
k：控制体边个数；
[0130]
f：控制体单元的中心；
[0131]
：时间编号m对应时间的边中心的总水深，单位：m；
[0132]
：时间编号m+1对应时间的边中心的总水深，单位：m；
[0133]
：时间编号m对应时间的x轴向的边中心的平均流速分量，单位：m/s；
[0134]
：时间编号m+1对应时间的x轴向的边中心的平均流速分量，单位：m/s；
[0135]
：时间编号m对应时间的y轴向的边中心的平均流速分量，单位：m/s；
[0136]
：时间编号m+1对应时间的y轴向的边中心的平均流速分量，单位：m/s；
[0137]
：时间编号m对应时间的单元中心的总水深，单位：m；
[0138]
：时间编号m对应时间的x轴向的单元中心的平均流速分量，单位：m/s；
[0139]
：时间编号m对应时间的y轴向的单元中心的平均流速分量，单位：m/s；
[0140]
：时间编号m对应时间的单元中心的水位，单位：m；
[0141]
da：控制单元面积，单位：m2。
[0142]
其中，
[0143]
步骤s13中所述的测量确定，具体为：
[0144]
先组建与布放观测设备，而后根据观测设备实时采集相应数据、计算完成相应数值的测量确定。
[0145]
其中，
[0146]
所述观测设备包括有：
[0147]
用于测量近底高分辨率流速场的多普勒声学流速剖面仪及脉冲耦合模式多普勒声学流速剖面仪；
[0148]
用于测量近底高精度流速过程与紊动过程的多普勒点流速仪；
[0149]
用于获取水体中悬浮物浊度的光学后向散射浊度仪；
[0150]
用于测量波浪和潮高的浪潮仪。
[0151]
实施例
[0152]
目前已有的水流阻力及阻力系数的获取方法描述如下：
[0153]
(1)一维恒定流条件下的阻力计算方法
[0154]
现有的一维恒定流条件下的阻力τ计算方程如下：
[0155][0156]
其中ρ为水密度，g为重力加速度，n即为阻力系数，u为垂线平均流速，h为垂线水深，这里n一般为实验结果或按经验取值。
[0157]
(2)二维恒定和非恒定流条件下的阻力计算方法
[0158]
二维恒定流和非恒定流的阻力计算一般借用式(1)的形式，分解成x和y两个方向的阻力，其计算可描述如下：
[0159][0160]
这里u和v分别为x和y轴两个方向的平均流速分量，n值一般直接选取一维恒定流时的实验室结果或者直接按经验取值。
[0161]
(3)阻力系数n试验室获取方法
[0162]
①
一维恒定流条件下阻力系数获取方法
[0163]
式(1)和(2)中最为关键是确定阻力系数的取值。一维情况下的n可以通过实验室获取，其主要通过变坡水槽或者平底水槽，其计算式如下：
[0164]
①
当变坡水槽时：
[0165][0166]
其中，i为水面坡降，当水面坡降和底坡一致时上式可通过试验测量数据求出阻力及阻力系数。
[0167]
②
当用平底水槽时：
[0168][0169]
其中表示偏导数，v＝(v1+v2)/2，h1、h2、u1、u2分别为图3中1、2断面的水深和平均流速，i为底坡，δl为计算断面长度，τ为阻力。
[0170]
②
二维非恒定流条件下的阻力系数获取方法
[0171]
目前采用的方法是直接用一维恒定条件下的阻力系数进行。
[0172]
(4)现有的阻力系数获取方法的存在的问题
[0173]
目前关于植被水流阻力系数的确定方法研究基本都是基于一维恒定的水流条件，但入海河口的实际水流具有非恒定性及二维特征，导致常用的水流阻力系数估算得到的阻力与实际会有较大的误差。
[0174]
另外，生态植被本身具有枝叶茂密枯萎、植被密度以及种类的多样性等生态特征，采用经验确定阻力系数的方法误差比较大。
[0175]
3、本实施例的技术方案
[0176]
(1)提出了一种基于现场观测的二维非恒定流阻力系数的计算新公式
[0177]
在二维非恒定流条件下，略去高阶项和科氏力等相对小值后，常用的守恒型动量方程可写成如下格式：
[0178][0179]
其中，u、v分别为x和y方向的流速，t为计算时间，η为水位，h为总水深，τ
x
、τ
y
为x和y方向的阻力。
[0180]
改写式(4)，则含有生态植被的水流阻力计算式描述成如下形式：
[0181][0182]
采用有限差分和有限体积法在控制体上进行式(5)的离散，可获得控制体单元的水流阻力，如下计算式描述：
[0183]
[0184]
其中，m为计算时间分步编号，δx
j
和δy
j
为控制体单元第j条边的x和y方向增量，k为控制体边个数，f表示控制体单元的中心，da为单元面积；
[0185]
根据式(2)和(6)，可得二维非恒定流的阻力系数n的计算式，描述如下：
[0186][0187]
因此，确定二维非恒定流阻力系数的计算公式的最终形式可描述如下：
[0188][0189]
(8)
[0190]
式(8)中计算需要的主要参数，式(8)中计算需要的主要参数，及δx
j
、δy
j
等，都可以很容易的通过现场测量的方式获取，详细测量方法参见下述。
[0191]
(2)设计了一套便于获取二维非恒定流阻力系数的现场测量方案
[0192]
本实施例先选用观测测量装置进行测量。所述观测测量装置由支架及设在支架上的多普勒声学流速剖面仪、脉冲耦合模式多普勒声学流速剖面仪、多普勒点流速仪、光学后向散射浊度仪及浪潮仪组成；
[0193]
所述支架呈正棱台状，分三层设置，其中，所述多普勒声学流速剖面仪及脉冲耦合模式多普勒声学流速剖面仪用于测量近底高分辨率流速场；所述多普勒点流速仪用于测量近底高精度流速过程与紊动过程；所述光学后向散射浊度仪用于获取水体中悬浮物浊度；所述浪潮仪用于测量波浪和潮高。
[0194]
通过对上述观测装置中各仪器的调制、布放与采集数据及数据处理等步骤，完成本实施例需要的观测点流速、水深和水位的观测。
[0195]
具体方案设计如下：为了获取式(8)中计算需要的节点(实心黑圆)位置处的数据，布置观测点位置如图4所示(三角形单元的m＝3，即测量4个点)，并把观测点包围形成的三角形区域作为计算生态植被水流阻力系数的计算单元，通过距离测量可获取δx
j
、δy
j
。
[0196]
本发明的一种确定入海河口生态植被水流阻力的方法，通过阻力方程与守恒型变形微分动量方程的联立，并结合实际测量，反算出时间序列上的二维非恒定流阻力系数，并据此阻力系数完成阻力值的确定，使得基于此的阻力值更能贴合与表征实际情况下的二维非恒定流水流阻力(该技术方案同样适用于二维恒定流以及一维恒定或一维非恒定流)；针对现况中的入海河口区域的水动力受径潮流影响，具有明显的非恒定性及二维特征，且入海河口的生态植被具有复杂的生物特征，导致入海河口生态植被水流阻力通常难以精确给定的问题，提供了一种便于获取二维非恒定流阻力系数的现场测量方案以及基于现场观测的二维非恒定流阻力系数的计算与确定方法，通过此阻力系数，提高入海河口生态植被水流阻力的计算准确度，为确定入海河口生态植被水流阻力提供了新的便捷、高精度和科学
的方法，对入海河口的湿地规划、管理及影响研究具有重要意义。