CN117993330A - 一种曲面河道非正交阻水范式计算方法 - Google Patents
一种曲面河道非正交阻水范式计算方法 Download PDFInfo
- Publication number
- CN117993330A CN117993330A CN202410400221.1A CN202410400221A CN117993330A CN 117993330 A CN117993330 A CN 117993330A CN 202410400221 A CN202410400221 A CN 202410400221A CN 117993330 A CN117993330 A CN 117993330A
- Authority
- CN
- China
- Prior art keywords
- local
- water
- engineering
- section
- orthogonal
- 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.)
- Pending
Links
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 87
- 230000000903 blocking effect Effects 0.000 title claims abstract description 49
- 238000004364 calculation method Methods 0.000 title claims abstract description 18
- 238000000034 method Methods 0.000 claims description 20
- 238000011144 upstream manufacturing Methods 0.000 claims description 3
- 238000010276 construction Methods 0.000 description 3
- 101100252016 Arabidopsis thaliana RMA2 gene Proteins 0.000 description 1
- WYTGDNHDOZPMIW-RCBQFDQVSA-N alstonine Natural products C1=CC2=C3C=CC=CC3=NC2=C2N1C[C@H]1[C@H](C)OC=C(C(=O)OC)[C@H]1C2 WYTGDNHDOZPMIW-RCBQFDQVSA-N 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000007599 discharging Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
Landscapes
- Measuring Volume Flow (AREA)
Abstract
本发明提供一种曲面河道非正交阻水范式计算方法，通过模拟河段水流流场，以河道内工程线路总长与线路起终点位置的流速偏角变化的比值作为划分依据，将非正交涉水工程沿线流向变化在一定范围内的单元格作为一个局部行洪断面，由此划分出多个局部行洪断面；计算局部行洪断面的平均流向，得到工程局部线路与局部行洪断面水流方向的夹角；根据各局部行洪断面内的涉水建筑物类型及结构尺寸，计算局部行洪断面内的工程阻水面积，最后得到工程总阻水面积及阻水比。
Description
技术领域
本发明涉及河道行洪技术领域，尤其涉及一种曲面河道非正交阻水范式计算方法。
背景技术
涉水工程的建设将占用河道原有的行洪断面，影响河道的行洪能力，对河道防洪安全造成影响。阻水比是衡量涉水工程建设对河道行洪影响的一个重要指标，阻水比是指设计洪水位条件下河道内工程建（构）筑物占用的过水断面面积（垂直于水流方向上的投影面积）与所在断面总面积的比值。对蜿蜒型、游荡型河道而言，设计洪水条件下，其过水断面为曲面，涉水工程线路基本无法做到正交于河道水流方向，基于此，阻水分析计算存在两个难点，一是河道行洪断面的确定，二是涉水建（构）筑物在河道行洪断面上的投影存在超出河道断面的问题。因此，如何计算非正交涉水工程对曲面河道局部行洪断面的阻水面积及阻水比，是本领域技术人员亟待解决的问题。
发明内容
本发明的目的在于针对上述现有技术的不足，提供了一种曲面河道非正交阻水范式计算方法，以正确评估工程建设对河道行洪能力的影响。
一种曲面河道非正交阻水范式计算方法，包括以下步骤：
S1、通过二维水动力模型模拟工程河段防洪标准条件下的水流流场，建立工程河段流场平面直角坐标系；
S2、将非正交涉水工程线路沿线流场均匀划分为个单元格，确定每个单元格的平均水深、平均流速；
S3、确定非正交涉水工程线路的局部行洪断面，计算局部行洪断面的行洪面积；
S4、计算各个局部行洪断面的平均流速，确定非正交涉水工程局部线路与局部行洪断面水流方向的夹角；
S5、确定各个局部行洪断面内涉水建筑物的类型及结构尺寸；
S6、计算各个局部行洪断面内的建筑物阻水面积、工程线路总阻水面积和工程阻水比。
进一步，所述S2中，确定每个单元格的平均水深和平均流速/>的组成集合为：
；
。
进一步，所述S3中，以河道内非正交涉水工程线路总长与线路起终点位置的流程偏角变化的比值作为确定局部行洪断面的依据，计算公式为：
；
其中，为非正交涉水工程线路的局部行洪断面宽度；/>为河道内非正交涉水工程线路总长；/>为线路终点位置的流程偏角；/>为线路起点位置的流程偏角；
当时，以/>为单位，将非正交涉水工程线路划分为多个局部行洪断面；
当时，以50m为单位，将非正交涉水工程线路划分为多个局部行洪断面。
进一步，提取所述工程沿线流场各个局部行洪断面的单元格数量，具体为：
；
其中，为局部行洪断面的单元格数量集；/>为局部行洪断面个数；/>为局部行洪断面的单元格数量；/>；
计算所述局部行洪断面的行洪面积，具体为：
；
其中，为第/>个局部行洪断面的行洪面积；/>为单元格的平均水深；/>为单元格的迎水面宽度；/>为/>个局部行洪断面的单元格数量；/>为第/>个局部行洪断面的单元格数量。
进一步，所述S4中，计算各个局部行洪断面的平均流速为：
；
其中，为第/>个局部行洪断面的平均流速；/>为第/>个局部行洪断面的单元格数量；/>为第/>个局部行洪断面的第1个网格编号；/>为第/>个局部行洪断面的最后一个网格编号。
进一步，所述S4中，确定非正交涉水工程局部线路与局部行洪断面水流方向的夹角具体为：
假设第个局部工程线路的方向向量为/>;由于第/>个局部工程线路能够看作直线，则有：
；
其中，分别为第/>个局部工程线路终点的横纵坐标；/>分别为第个局部工程线路起点的横纵坐标，第/>个局部工程线路为第/>个局部行洪断面对应的工程；
则：；
；
其中，为非正交涉水工程局部线路与局部行洪断面水流方向的夹角；分别为第/>个局部行洪断面平均流速和第/>个局部工程线路的模。
进一步，所述S5具体为：假设第个局部行洪断面中有/>个方形阻水建筑物和/>个圆形阻水建筑物；方形阻水建筑物尺寸为/>；圆形阻水建筑物直径为/>；
提取各局部行洪断面内方形阻水建筑物和圆形阻水建筑物所在位置的平均水深：
；
；
其中，分别为各个局部行洪断面内的方形、圆形建筑物所在位置的水深。
进一步，所述S6中，计算各个局部行洪断面内的建筑物阻水面积和工程线路总阻水面积，具体为：
；
；
根据局部工程线路与局部行洪断面的水流方向夹角，以及阻水建筑物的尺寸，分别计算局部行洪断面内的建筑物阻水面积：
；
令：；
；
；
则，工程线路总阻水面积和工程阻水比为：
；
；
其中，、/>为第/>个局部行洪断面内方形建筑物阻水面积；/>为第/>个局部行洪断面内圆形建筑物阻水面积；/>为非正交涉水工程阻水比；/>为非正交涉水工程总的阻水面积；/>为第/>个局部行洪断面与第/>个局部工程线路的夹角。
附图说明
图1为一种曲面河道非正交阻水范式计算方法的流程图；
图2为实施例中模拟的流场图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，下面结合附图，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本发明，并不用于限定本发明。
请参阅图1，以东荆河某工程为例，具体说明一种曲面河道非正交阻水范式计算方法，包括以下步骤：
S1、通过二维水动力模型模拟工程河段设计防洪标准条件下的水流流场，建立工程河段流场平面直角坐标系；
二维水动力学模型，其包括但不限于Delft 3D、MIKE 21、RMA2；所述的建立工程河段流速场平面直角坐标系，其包括但不限于2000国家大地坐标系、1954年北京坐标系、1980年西安坐标系。
具体的，模拟工程河段水水流流场可通过不同水动力模型模拟，本次实施例采用Mike21模型模拟流场，流场示例见图2，其中，流速流场平面坐标系采用2000国家大地坐标系。
S2、将非正交涉水工程线路沿线流场均匀划分为个单元格，确定每个单元格的平均水深、平均流速；
S3、确定非正交涉水工程线路的局部行洪断面，计算局部行洪断面的行洪面积；
S4、计算各个局部行洪断面的平均流速，确定非正交涉水工程局部线路与局部行洪断面水流方向的夹角；
S5、确定各个局部行洪断面内涉水建筑物的类型及结构尺寸；
S6、计算各个局部行洪断面内的建筑物阻水面积、工程线路总阻水面积和工程阻水比。
所述S2中，确定每个单元格的平均水深和平均流速/>的组成集合为：
；
；
其中，。
所述S3中，以河道内非正交涉水工程线路总长与线路起终点位置的流程偏角变化的比值作为确定非正交涉水工程线路的局部行洪断面的依据，计算公式为：
；
其中，为非正交涉水工程线路的局部行洪断面宽度；/>为河道内工程线路总长；为线路终点位置的流程偏角；/>为线路起点位置的流程偏角；
其中，； />，计算得到/>。
当时，以/>为单位，将非正交涉水工程线路划分为多个局部行洪断面；
当时，以50m为单位，将非正交涉水工程线路划分为多个局部行洪断面；
以50m为单位将工程沿线划分为多个局部行洪断面，其中。
提取所述工程沿线流场各个局部行洪断面的单元格数量，具体为：
；
其中，为局部行洪断面的单元格数量集；/>为局部行洪断面个数；/>为局部行洪断面的单元格数量；/>；
计算所述局部行洪断面的行洪面积，具体为：
；
其中，为第/>个局部行洪断面的行洪面积；/>为单元格的平均水深；/>为单元格的迎水面宽度；/>为/>个局部行洪断面的单元格数量；/>为第/>个局部行洪断面的单元格数量。
以第75个局部行洪断面为例说明，其中、/>、/>、。
所述S4中，计算各个局部行洪断面的平均流速为：
；
其中，为第/>个局部行洪断面的平均流速；/>为第/>个局部行洪断面的单元格数量；/>为第/>个局部行洪断面的第1个网格编号；/>为第/>个局部行洪断面的最后一个网格编号 。
以第75个局部行洪断面为例，其中、/>、/>、。
所述S4中，确定非正交涉水工程局部线路与局部行洪断面水流方向的夹角具体为：
假设第个局部工程线路的方向向量为/>；由于第/>个局部工程线路能够看作直线，则有：
；
其中，分别为第/>个局部工程线路终点的横纵坐标；/>分别为第个局部工程线路起点的横纵坐标，第/>个局部工程线路为第/>个局部行洪断面对应的工程；
则： ；
；
其中，为非正交涉水工程局部线路与局部行洪断面水流方向的夹角； 分别为第/>个局部行洪断面平均流速和第/>个局部工程线路的模。
以第75个局部行洪断面为例，其中，，/> 。
所述S5具体为：假设第个局部行洪断面中有/>个方形阻水建筑物和/>个圆形阻水建筑物；方形阻水建筑物尺寸为/>；圆形阻水建筑物直径为/>；提取各局部行洪断面内方形阻水建筑物和圆形阻水建筑物所在位置的平均水深/>：
；
；
其中，分别为各个局部行洪断面内的方形、圆形建筑物所在位置的水深。
其中，本实施例为长方形和圆形桥墩的组合，长方形桥墩总个数为10个，尺寸均为，圆形桥墩共131个，尺寸均为/>；以75个局部行洪断面为例，局部行洪断面内分布有2个圆形桥墩，其水深为/>。
所述S6中，计算各个局部行洪断面内的建筑物阻水面积和工程线路总阻水面积，具体为：
；
；
根据局部工程线路与局部行洪断面的水流方向夹角，以及阻水建筑物的尺寸，分别计算局部行洪断面内的建筑物阻水面积：
；
令：；
；
；
则，工程线路总阻水面积和工程阻水比为：
；
；
其中，、/>为第/>个局部行洪断面内方形建筑物阻水面积；/>为第/>个局部行洪断面内圆形建筑物阻水面积；/>为非正交涉水工程阻水比；/>为非正交涉水工程总的阻水面积；/>为第/>个局部行洪断面与第/>个局部工程线路的夹角。
本实施例中，、/>。
以上所述实施例仅表达了本发明的实施方式，其描述较为具体和详细，但并不能因此而理解为对本发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变形和改进，这些都属于本发明的保护范围。因此，本发明专利的保护范围应以所附权利要求。
Claims (8)
1.一种曲面河道非正交阻水范式计算方法，其特征在于，包括以下步骤：
S1、通过二维水动力模型模拟工程河段防洪标准条件下的水流流场，建立工程河段流场平面直角坐标系；
S2、将非正交涉水工程线路沿线流场均匀划分为个单元格，确定每个单元格的平均水深、平均流速；
S3、确定非正交涉水工程线路的局部行洪断面，计算局部行洪断面的行洪面积；
S4、计算各个局部行洪断面的平均流速，确定非正交涉水工程局部线路与局部行洪断面水流方向的夹角；
S5、确定各个局部行洪断面内涉水建筑物的类型及结构尺寸；
S6、计算各个局部行洪断面内的建筑物阻水面积、工程线路总阻水面积和工程阻水比。
2.根据权利要求1所述的一种曲面河道非正交阻水范式计算方法，其特征在于：所述S2中，确定每个单元格的平均水深和平均流速/>的组成集合为：
；
。
3.根据权利要求2所述的一种曲面河道非正交阻水范式计算方法，其特征在于：所述S3中，以河道内非正交涉水工程线路总长与线路起终点位置的流程偏角变化的比值作为确定局部行洪断面的依据，计算公式为：
；
其中，为非正交涉水工程线路的局部行洪断面宽度；/>为河道内非正交涉水工程线路总长；/>为线路终点位置的流程偏角；/>为线路起点位置的流程偏角；
当时，以/>为单位，将非正交涉水工程线路划分为多个局部行洪断面；
当时，以50m为单位，将非正交涉水工程线路划分为多个局部行洪断面。
4.根据权利要求3所述的一种曲面河道非正交阻水范式计算方法，其特征在于：提取所述工程沿线流场各个局部行洪断面的单元格数量，具体为：
；
其中，为局部行洪断面的单元格数量集；/>为局部行洪断面个数；/>为局部行洪断面的单元格数量；/>；
计算所述局部行洪断面的行洪面积，具体为：
；
其中，为第/>个局部行洪断面的行洪面积；/>为单元格的平均水深；/>为单元格的迎水面宽度；/>为/>个局部行洪断面的单元格数量；/>为第/>个局部行洪断面的单元格数量。
5.根据权利要求4所述的一种曲面河道非正交阻水范式计算方法，其特征在于：所述S4中，计算各个局部行洪断面的平均流速为：
；
其中，为第/>个局部行洪断面的平均流速；/>为第/>个局部行洪断面的单元格数量；/>为第/>个局部行洪断面的第1个网格编号；/>为第/>个局部行洪断面的最后一个网格编号。
6.根据权利要求5所述的一种曲面河道非正交阻水范式计算方法，其特征在于：所述S4中，确定非正交涉水工程局部线路与局部行洪断面水流方向的夹角具体为：
假设第个局部工程线路的方向向量为/>;由于第/>个局部工程线路能够看作直线，则有：
；
其中，分别为第/>个局部工程线路终点的横纵坐标；/>分别为第/>个局部工程线路起点的横纵坐标，第/>个局部工程线路为第/>个局部行洪断面对应的工程；
则：；
；
其中，为非正交涉水工程局部线路与局部行洪断面水流方向的夹角； />分别为第/>个局部行洪断面平均流速和第/>个局部工程线路的模。
7.根据权利要求6所述的一种曲面河道非正交阻水范式计算方法，其特征在于，所述S5具体为：假设第个局部行洪断面中有/>个方形阻水建筑物和/>个圆形阻水建筑物；方形阻水建筑物尺寸为/>；圆形阻水建筑物直径为/>；
提取各局部行洪断面内方形阻水建筑物和圆形阻水建筑物所在位置的平均水深：
；
；
其中，分别为各个局部行洪断面内的方形、圆形建筑物所在位置的水深。
8.根据权利要求7所述的一种曲面河道非正交阻水范式计算方法，其特征在于：所述S6中，计算各个局部行洪断面内的建筑物阻水面积和工程线路总阻水面积，具体为：
；
；
根据局部工程线路与局部行洪断面的水流方向夹角，以及阻水建筑物的尺寸，分别计算局部行洪断面内的建筑物阻水面积：
；
令：；
；
；
则，工程线路总阻水面积和工程阻水比为：
；
；
其中，、/>为第/>个局部行洪断面内方形建筑物阻水面积；/>为第/>个局部行洪断面内圆形建筑物阻水面积；/>为非正交涉水工程阻水比；/>为非正交涉水工程总的阻水面积；/>为第/>个局部行洪断面与第/>个局部工程线路的夹角。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410400221.1A CN117993330A (zh) | 2024-04-03 | 2024-04-03 | 一种曲面河道非正交阻水范式计算方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410400221.1A CN117993330A (zh) | 2024-04-03 | 2024-04-03 | 一种曲面河道非正交阻水范式计算方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN117993330A true CN117993330A (zh) | 2024-05-07 |
Family
ID=90900993
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202410400221.1A Pending CN117993330A (zh) | 2024-04-03 | 2024-04-03 | 一种曲面河道非正交阻水范式计算方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117993330A (zh) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111651885A (zh) * | 2020-06-03 | 2020-09-11 | 南昌工程学院 | 一种智慧型海绵城市洪涝预报方法 |
| CN112257352A (zh) * | 2020-10-21 | 2021-01-22 | 黄河水利委员会黄河水利科学研究院 | 一维水动力模型和二维水动力模型的耦合方法及系统 |
| CN112507419A (zh) * | 2020-11-19 | 2021-03-16 | 长江勘测规划设计研究有限责任公司 | 无资料、多碍洪构筑物的山区性河流洪水过程模拟方法 |
| CN114511995A (zh) * | 2022-01-05 | 2022-05-17 | 四川岷江港航电开发有限责任公司 | 一种基于cesium模型的洪水分级预警方法 |
| US20220413181A1 (en) * | 2021-06-29 | 2022-12-29 | One Concern, Inc. | Coupled pluvial, fluvial, and urban flood tool |
-
2024
- 2024-04-03 CN CN202410400221.1A patent/CN117993330A/zh active Pending
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111651885A (zh) * | 2020-06-03 | 2020-09-11 | 南昌工程学院 | 一种智慧型海绵城市洪涝预报方法 |
| CN112257352A (zh) * | 2020-10-21 | 2021-01-22 | 黄河水利委员会黄河水利科学研究院 | 一维水动力模型和二维水动力模型的耦合方法及系统 |
| CN112507419A (zh) * | 2020-11-19 | 2021-03-16 | 长江勘测规划设计研究有限责任公司 | 无资料、多碍洪构筑物的山区性河流洪水过程模拟方法 |
| US20220413181A1 (en) * | 2021-06-29 | 2022-12-29 | One Concern, Inc. | Coupled pluvial, fluvial, and urban flood tool |
| CN114511995A (zh) * | 2022-01-05 | 2022-05-17 | 四川岷江港航电开发有限责任公司 | 一种基于cesium模型的洪水分级预警方法 |
Non-Patent Citations (2)
| Title |
|---|
| JIANPING BING 等: "Hydrological changes along the Jing River and its causal analysis", 《2020 6TH INTERNATIONAL CONFERENCE ON HYDRAULIC AND CIVIL ENGINEERING》, 31 December 2020 (2020-12-31), pages 1 - 10 * |
| 陈勇 等: "城市河道生态治理及修复措施浅析", 《水利技术监督》, no. 3, 31 March 2024 (2024-03-31), pages 260 - 262 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108021780B (zh) | 一种基于无规则非结构网格模型的山洪动态仿真方法 | |
| CN108256193A (zh) | 一种河道溃决洪水一、二维数学模型耦合分析方法 | |
| CN112257352A (zh) | 一维水动力模型和二维水动力模型的耦合方法及系统 | |
| CN104361157A (zh) | 一种建筑物间风环境的评价方法 | |
| CN112036043B (zh) | 长距离重力流大管径供水管道破裂模型计算的方法 | |
| CN110532648A (zh) | 一种堰塞湖风险评估平台 | |
| CN107256312B (zh) | 一种基于潮流环境下海湾纳潮变化量计算方法 | |
| CN110147646B (zh) | 一种数值模拟框架下线性挡水构筑物的过流处理方法 | |
| Wang et al. | Three-dimensional flood routing of a dam break based on a high-precision digital model of a dense urban area | |
| CN113297753B (zh) | 山区河流透过性坝沿程水深预测方法 | |
| CN113095552B (zh) | 一种海洋平台泄漏可燃气云团体积预测方法及系统 | |
| CN117993330A (zh) | 一种曲面河道非正交阻水范式计算方法 | |
| CN113191054A (zh) | 一种基于显卡加速耦合管网的高精度城市雨洪模拟方法 | |
| CN116822180A (zh) | 一种模拟河渠输水过程的积分时滞模型参数辨识方法 | |
| CN111079330A (zh) | 一种考虑流固耦合效应的大跨度膜结构抗风设计方法 | |
| CN110412312A (zh) | 涌潮流速垂向分布的一种计算方法 | |
| CN103258085A (zh) | 一种污水深海排放扩散器上升管数量确定的方法 | |
| Fatimah et al. | Numerical simulation of groyne placement in minimising Krueng Aceh river bank erosion | |
| CN109724570B (zh) | 地下跌水的跌水量、跌水宽度、坎上水层厚度的计算方法 | |
| CN111125805A (zh) | 一种基于数字高程模型和河流流向的无资料地区河道断面信息识别方法 | |
| He et al. | A CFD Analysis of Wall-Effect Building Groups in a Curved Layout: A Study in Sha Tin, Hong Kong | |
| CN115408832B (zh) | 河岸生态护坡抗冲流速复核方法 | |
| CN114492238A (zh) | 一种河湖一维与平面二维水动力模型深度耦合方法 | |
| CN117150823A (zh) | 一种滨海城市洪水灾害的高分辨率快速模拟方法 | |
| CN114372427A (zh) | 一种水动力模型的建立方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination |