一种适用于径流潮汐河口段的河相关系及其推导方法技术

本发明专利技术公开了一种适用于径流潮汐河口段的河相关系及其推导方法，其特征在于，包括如下步骤：S1、结合O'Brien方法和窦国仁方法，初步推导河相关系；S2、根据河流的实测地形资料，建立水流数学模型，并采用有限体积法对其求解；S3、选取仅受径流作用的河段，拟合该径流河段的河相关系；S4、选取具有代表性的上游径流流量的工程河段；S5、将该河段分为一定数量的断面，若干等分各个断面的河宽，通过步骤S3拟合的径流河段的河相关系计算；S6、利用步骤S2中的水流数学模型，分别计算步骤S4中的河段的平均过水面积A与QT落；本发明专利技术的河相关系式在定性和定量描述径流潮汐河口段河流动力条件与河槽形态的关系时具有很高的精度。

【技术实现步骤摘要】

本专利技术涉及一种河相关系，具体涉及一种适用于径流潮汐河口段的河相关系。
技术介绍
径流潮汐河口段的河床形态是流域来水来沙、海洋潮汐动力等相互作用的结果，其间存在着一定关系，通常称为河相关系。河相关系的研究旨在寻求河床地貌形态与动力因子的定量关系，是认识河道特性的主要技术手段之一，一直为人们所关注。有关河相关系，国内外曾先后开展过大量研究，特别是有关潮量或流量与过水面积之间的关系，成果颇丰。1930年，O'Brien在调查美国太平洋砂质侵蚀海岸潮汐通道时，提出潮量和断面面积(相对于平均海平面)之间的表达式；1953年，Leopold和Maddock提出在准平衡状态的天然河流中，给定的河宽、平均水深、平均流速与流量之间存在着简单类似的函数关系，也就是河相关系；1964年，窦国仁将河床最小活动性原理应用到潮汐河流，提出平均水深、平均河宽和平均断面面积与落潮平均流量、落潮平均含沙量河相关系式；1988年，高抒对中国东海11个海湾通道资料进行统计后亦得到断面面积与潮量的关系；2010年，韩曾萃等在海湾通道断面面积与潮量关系的基础上，进一步将径流、落潮历时、落潮含沙量的河口河相关系应用于海湾，并验证了浙江省的杭州湾、象山港、三门湾、乐清湾等区域25个断面的河相关系；2014年，赵子辉利用最小活动性原理结合新的水流挟沙力公式，重构窦国仁河相关系公式，建立了一个新的适用于潮汐河口的河相关系公式；2015年，张玮和曹昊建立了灌河河口落潮流量与过水断面面积之间的关系，并对不同河道形态的理想宽深比及其在航道整治中的应用进行了讨论。上述研究中，主要针对径流或者潮流单独作用的情况，即使有些研究也涉及径流潮汐河口段，没有深入探讨径流和潮流的河相关系。
技术实现思路
为解决现有技术的不足，本专利技术的目的在于提供一种适用于径流潮汐河口段的河相关系。为了实现上述目标，本专利技术采用如下的技术方案：一种适用于径流潮汐河口段的河相关系的推导方法，其特征在于，包括如下步骤：S1、结合O'Brien方法和窦国仁方法，初步推导河相关系；S2、根据河流的实测地形资料，建立河流的一维潮流和二维潮流的耦合平面水流数学模型，并采用有限体积法对其求解；S3、选取上游仅受径流作用的河段，根据上述的水流数学模型，模拟该河段的一年水流量与过水断面面积的关系，拟合径流河段的河相关系；S4、选取计算工况：由于工程河段在一个潮汐周期中，上游径流量几乎不变，可以将其设为定值，因此可选取若干具有代表性的上游径流流量的该工程河段；S5、计算选取断面：根据步骤S2河流的实测地形资料，将该河段分为一定数量的断面，将各个断面的河宽分为若干等分，通过步骤S3拟合的径流河段的河相关系计算；S6、验证河相关系式：利用步骤S2中的水流数学模型，分别计算步骤S4中的代表性工况下该河段的平均过水面积A与QT落；上述步骤S1中的O'Brien方法，潮量和断面面积之间的表达式为：式中：P为平均大潮潮差下的潮量；A1为平均海平面下河口断面面积；k为由河口断面特性决定的常数；S为含沙量；所述窦国仁方法，如果潮汐河口在一个潮汐周期T内的冲淤数量能够相互抵偿，在此时期内上游来水量以及下游海口涨入的水量必然等于同一时期内落潮期间河床断面所能排泄的水量。令P′表示T时间内的上游来水量，同一时期内涨潮流期间下游涨入的水量为潮量P；为了落潮期间能够排泄上述水量，河床必须具有如下水力几何形态：P+P′＝BHvT落(19)式中：B为平均潮水位时的水面宽，H为相应的平均水深，v为平均落潮流速，T落为单个潮汐周期T内的总落潮历时；所述推导过程为：用Q表示落潮流(包括上游径流)的平均流量，用Q0表示落潮期间上游平均径流量，则有：P+Q0T＝QT落(20)，将(1)式代入(3)式，并经转换可得：其中：A1＝A-A0(22)，式中：A为平均过水断面面积，A0为径流产生的平均过水断面面积；A1为潮流产生的平均过水断面面积；对于径流潮流河段上游的冲积性河流，仅受径流作用，有如下关系式：式中：B′为平均河宽；H′为相应平均水深；k1、k2、m、n为常数；由于A0＝B′H′(25)，将(6)式、(7)式代入(8)式可得：式中：k3、e为率定参数；将(9)式代入(5)式得：将(10)式代入(4)式得：或者上述步骤S2中的一维潮流数学模型中，采用Saint-Venant方程组作为控制方程，其连续方程和运动方程可以表示为：式中：x和t分别为距离和时间的坐标，A为平均过水断面面积，Q为流量，h为水位，q为旁侧入流的流量，C为谢才系数，R为水力半径，α为动量校正系数，g为重力加速度；二维潮流数学模型中，连续方程为：运动方程为：式中：x,y为笛卡尔坐标，t为时间变量，η为水位，d代表总水深，且有h＝d+η，d为静水深，为x,y方向深度平均速度，τbx、τby为x,y方向底部应力，ρ0为水的密度，f为科氏力系数，且f＝2Ωsinψ，g为当地重力加速度，Txx、Tyy、Txy为水平粘滞应力项；在进行数值计算时，对于开边界，上、下游分别采用流量边界和水位边界条件；对于闭边界，根据不可入原理取法向流速为0；采用干湿判别技术进行动边界处理。上述步骤S3中径流河段的河相关系的拟合的方法：平均过水断面面积的计算方法是将断面的河宽20等分，由步骤S2的水流数学模型计算此河宽的21个坐标点在平均潮位下的水深，运用以下公式：式中：A0为径流产生的平均过水断面面积；h1，h21为该断面两端坐标点在平均潮位下的水深；hi为该断面上第i个坐标点在平均潮位下的水深；B为平均落潮水位时的河宽。上述步骤S6利用步骤S2中的数学模型，分别计算步骤S4中几个代表性工况下该工程河段的平均过水面积A与QT落的方法：其中：式中：t为某断面落潮期间里任意一个整点时刻；t0为该断面落潮起始时刻；t1为该断面落潮结束时刻；Qt为该断面t时刻的流量；将(QT落-Q0T)与(A-A0)1.1作为散点，点绘于平面(x,y)坐标系，进行线性回归分析，率定参数kS，并验证该河相关系式。本专利技术的有益之处在于：通过基于O'Brien提出的海湾通道断面面积与潮量关系式，考虑潮汐河口段上游径流量，参考窦国仁有关河口动力的处理方法，推导出一个适用于径流潮汐河口段河相关系式，该河相关系式可以分别考虑径流与潮流的贡献。利用数学模型模拟计算径流潮汐河口段的各物理量，验证推导的该河相关系式，该河相关系式在定性和定量描述径流潮汐河口段河流动力条件与河槽形态的关系时具有很高的精度，具有很强的实用性。附图说明图1是本专利技术中仅受径流作用下长江下游河段断面过水断面面积与流量关系图。图2是本专利技术的计算断面的选取的示意图。图3是本专利技术的澄通河段河相关系图(1#～15#)。具体实施方式以下结合附图和具体实施例对本专利技术作具体的介绍。以本专利技术提出的一种适用于径流潮汐河口段的河相关系应用于长江下游典型径流潮汐河口段澄通河段为例来详细介绍本专利技术，实施过程分为六个步骤。步骤S1：初步推导该河相关系参照O'Brien方法，潮量和断面面积之间的表达式为：式中：P为平均大潮潮差下的潮量；A1为平均海平面下河口断面面积；k为常数(河口断面特性决定)；S为含沙量。参照窦国仁方法，如果潮汐河口在一个潮汐周期T内的冲淤数量能够相互抵偿，在此时期内上游来本文档来自技高网...

【技术保护点】
一种适用于径流潮汐河口段的河相关系的推导方法，其特征在于，包括如下步骤：S1、结合O'Brien方法和窦国仁方法，初步推导河相关系；S2、根据河流的实测地形资料，建立河流的一维潮流和二维潮流的耦合平面水流数学模型，并采用有限体积法对其求解；S3、选取仅受径流作用的河段，根据上述的水流数学模型，模拟该河段的一年水流量与过水断面面积的关系，拟合径流河段的河相关系；S4、选取计算工况：由于工程河段在一个潮汐周期中，上游径流量几乎不变，可以将其设为定值，因此可选取若干具有代表性的上游径流流量的该工程河段；S5、计算选取断面：根据步骤S2河流的实测地形资料，将该河段分为一定数量的断面，将各个断面的河宽分为若干等分，通过步骤S3拟合的径流河段的河相关系计算；S6、验证河相关系式：利用步骤S2中的水流数学模型，分别计算步骤S4中的代表性工况下该河段的平均过水面积A与QT落；

【技术特征摘要】
1.一种适用于径流潮汐河口段的河相关系的推导方法，其特征在于，包括如下步骤：S1、结合O'Brien方法和窦国仁方法，初步推导河相关系；S2、根据河流的实测地形资料，建立河流的一维潮流和二维潮流的耦合平面水流数学模型，并采用有限体积法对其求解；S3、选取仅受径流作用的河段，根据上述的水流数学模型，模拟该河段的一年水流量与过水断面面积的关系，拟合径流河段的河相关系；S4、选取计算工况：由于工程河段在一个潮汐周期中，上游径流量几乎不变，可以将其设为定值，因此可选取若干具有代表性的上游径流流量的该工程河段；S5、计算选取断面：根据步骤S2河流的实测地形资料，将该河段分为一定数量的断面，将各个断面的河宽分为若干等分，通过步骤S3拟合的径流河段的河相关系计算；S6、验证河相关系式：利用步骤S2中的水流数学模型，分别计算步骤S4中的代表性工况下该河段的平均过水面积A与QT落；2.根据权利要求1所述的一种适用于径流潮汐河口段的河相关系的推导方法，其特征在于，所述步骤S1中的O'Brien方法，潮量和断面面积之间的表达式为：P=A11.1k1.1S1.1---(1),]]>式中：P为平均大潮潮差下的潮量；A1为平均海平面下河口断面面积；k为由河口断面特性决定的常数；S为含沙量；所述窦国仁方法，如果潮汐河口在一个潮汐周期T内的冲淤数量能够相互抵偿，在此时期内上游来水量以及下游海口涨入的水量必然等于同一时期内落潮期间河床断面所能排泄的水量。令P′表示T时间内的上游来水量，同一时期内涨潮流期间下游涨入的水量为潮量P；为了落潮期间能够排泄上述水量，河床必须具有如下水力几何形态：P+P′＝BHvT落(2)式中：B为平均落潮水位时的河宽，H为相应的平均水深，v为平均落潮流速，T落为单个潮汐周期T内的总落潮历时；所述推导过程为：用Q表示落潮流(包括上游径流)的平均流量，用Q0表示落潮期间上游平均径流量，则有：P+Q0T＝QT落(3)，将(1)式代入(3)式，并经转换可得：其中：A1＝A-A0(5)，式中：A为平均过水断面面积，A0为径流产生的平均过水断面面积；A1为潮流产生的平均过水断面面积；对于径流潮流河段上游的冲积性河流，仅受径流作用，有如下关系式：B′=k1Q0m---(6),]]>H′=k2Q0n---(7),]]>式中：B′为平均河宽；H′为相应平均水深；k1、k2、m、n为常数；由于A0＝B′H′(8)，将(6)式、(7)式代入(8)式可得：A0=k3Q0e---(9),]]>式中：k3、e为率定参数；将(9)式代入(5)式得：A1=A-k3Q0e---(10),]]>将(10)式代入(4)式得：或者3.根据权利要求1所述的一种适用于径流潮汐河口段的河相关系的推导方法，其特征在于，所述步骤S2中的一维潮流数学模型中，采用Saint-Venant方程组作为控制方程，其连续方程和运动方程可以表示为：∂A∂t+∂Q∂x=q∂Q∂t+&pa...

【专利技术属性】
技术研发人员：张玮，徐鹏，徐晶鑫，周凯，覃昌佩，吴彦颖，徐延琦，王超，
申请(专利权)人：河海大学，
类型：发明
国别省市：江苏;32