CN108875222A - 基于水动力过程相似的水文模型流域比例尺确定方法 - Google Patents
基于水动力过程相似的水文模型流域比例尺确定方法 Download PDFInfo
- Publication number
- CN108875222A CN108875222A CN201810652924.8A CN201810652924A CN108875222A CN 108875222 A CN108875222 A CN 108875222A CN 201810652924 A CN201810652924 A CN 201810652924A CN 108875222 A CN108875222 A CN 108875222A
- Authority
- CN
- China
- Prior art keywords
- basin
- model
- research
- watershed
- scale
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种基于水动力过程相似的水文模型流域比例尺确定方法,其包括获取研究流域的矢量信息,并将室内自定义坐标的矢量文件与研究流域的边界矢量文件进行匹配,得到模型流域的水平比例尺;基于水动力相似原则,采用研究流域的长和宽及水平比例尺分别计算研究流域和模型流域的水文参数;获取阻力相似模型,并采用研究流域和模型流域的水文参数及阻力相似模型得到垂直比例尺。
Description
技术领域
本发明涉及水利工程技术领域,具体涉及一种基于水动力过程相似的水文模型流域比例尺确定方法。
背景技术
区域水文过程试验是当前水文试验研究的主要发展方向之一,开展室内区域水文过程试验首先需要根据试验研究目标科学确定室内区域水文过程试验的物理模型与研究流域之间的转换比例尺。由于水循环伴生过程涉及到水文、泥沙、环境、生态等多个方面,目前以泥沙和污染物运移等为主要研究目标的、以水动力过程为基础的模型试验方法中,确定室内物理模型比尺方法多以单一过程或点尺度为主。由于两者之间的关注点不一致,当前科学设定区域人工降雨过程的方法还不多见,这样可能会致使计算的与降雨直接相关参数的准确性,最终导致模拟的人工降雨难以达到预期的效果。
发明内容
针对现有技术中的上述不足,本发明提供的基于水动力过程相似的水文模型流域比例尺确定方法能够通过阻力相似模型和水文参数准确地计算出模型流域的垂直比例尺。
为了达到上述发明目的,本发明采用的技术方案为:
提供一种基于水动力过程相似的水文模型流域比例尺确定方法,其包括:
获取研究流域的矢量信息,并将室内自定义坐标的矢量文件与研究流域的边界矢量文件进行匹配,得到模型流域的水平比例尺;
基于水动力相似原则,采用研究流域的长和宽及水平比例尺分别计算研究流域和模型流域的水文参数;
获取阻力相似模型,并采用研究流域和模型流域的水文参数及阻力相似模型得到垂直比例尺。
进一步地,所述水文参数包括出口断面面积、湿周、水力半径、自分水岭到出口断面主河道的平均比降和断面出口流速。
进一步地,所述研究流域的水文参数的计算公式为:
A=LH;WP=2H+L;
其中,L为研究流域的长;H为研究流域的宽;WP为室外研究流域的湿周;R为室外研究流域的水力半径;J为研究流域自分水岭到出口断面主河道的平均比降;n为研究流域的糙率;
h1、h2、h3······hn为研究流域河道纵断面各转折点的河底高程与设计断面处河底高程的差值;L1、L2、L3······Ln为研究流域各转折点之间的分段河长,L′为主河长;v为室外研究流域的断面出口流速。
进一步地,所述模型流域的水文参数的计算公式为:
A′=K1K2A=K1K2LH;W′P=2K2H+K1L;
其中,A′为室内模型流域的出口断面面积;WP′为室内模型流域的湿周;R′为室内模型流域的水力半径;J′为模型流域自分水岭到出口断面主河道的平均比降;n′为模型流域的糙率;K1为水平比例尺;K2为垂直比例尺;
h′1、h′2、h′3······h′n为模型流域河道纵断面各转折点的河底高程与设计断面处河底高程的差值,且h′n=K2hn;L′1、L′2、L′3······L′n为模型流域各转折点之间的分段河长,L′n=K1Ln;v′为室内模型流域的断面出口流速。
进一步地,所述阻力相似模型的计算公式为:
其中,αu为流速比尺,αL为水平几何比尺,αL=K1;αn为糙率比尺,αR是水力半径比尺,
进一步地,根据模型流域和研究流域的水文参数的计算公式及阻力相似计算公式,得到垂直比例尺的计算公式:
进一步地,所述研究流域的矢量信息包括研究流域形状矢量信息和雨量站位置的矢量信息;
所述研究流域的矢量信息提取方法包括:获取研究流域的地形数据和雨量站位置;将所述地形数据与雨量站位置导入ArcGis中得到研究流域形状和雨量站位置的矢量信息。
本发明的有益效果为:本方案基于水动力过程相似和阻力相似,通过计算得到研究流域和模型流域的水文参数,可以知道两者水文参数之间关系,再结合构建的阻力相似模型就能够快速地计算出垂直比例尺。
由于构建的模型流域的垂直比例尺计算公式直接与水平比例尺及研究流域的高和宽相关,一旦得到水平比例尺后,通过该函数能够快速地得到不同地理特性的研究流域在室内研究区域进行人工模拟时的垂直比例尺。
另外,通过垂直比例尺和水平比例尺对研究流域缩放得到的模型流域计算的降雨控制参数更能够反应研究流域降雨的真实情况,从而保证了模拟的人工降雨的预期的效果。
附图说明
图1为基于水动力过程相似的水文模型流域比例尺确定方法的流程图。
图2为研究流域的示意图。
图3为模型流域导入室内降雨区的示意图。
具体实施方式
下面对本发明的具体实施方式进行描述,以便于本技术领域的技术人员理解本发明,但应该清楚,本发明不限于具体实施方式的范围,对本技术领域的普通技术人员来讲,只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内,这些变化是显而易见的,一切利用本发明构思的发明创造均在保护之列。
参考图1,图1示出了基于水动力过程相似的水文模型流域比例尺确定方法的流程图;如图1所示,该方法包括步骤101至步骤103。
在步骤101中,获取研究流域的矢量信息,并将室内自定义坐标的矢量文件与研究流域的边界矢量文件进行匹配,得到模型流域的水平比例尺;
在实施时,本方案优选所述研究流域的矢量信息包括研究流域形状矢量信息和雨量站位置的矢量信息;
所述研究流域的矢量信息提取方法包括:获取研究流域的地形数据和雨量站位置;将所述地形数据与雨量站位置导入ArcGis中得到研究流域形状和雨量站位置的矢量信息。
图2示出的研究流域根据水平比例尺和垂直比例尺导入室内降雨区后可以得到如图3所示的示意图。其中,室内降雨区(区域降雨人工模拟系统)为:有效降雨面积26m×40m,共分为110个独立控制单元(其中共分为10个大区,每个大区又分为11个小区);每个降雨单元雨强变化范围10-200mm/h;单元内均匀度大于0.8;雨强变化频率可以做到5分钟调整一次;通过输入降雨控制文件,可以实现人工降雨。
在步骤102中,基于水动力相似原则,采用研究流域的长和宽及水平比例尺分别计算研究流域和模型流域的水文参数。
其中的水文参数包括出口断面面积、湿周、水力半径、自分水岭到出口断面主河道的平均比降、断面出口流速和研究流域单位时间通过出口断面的流量,本方案主要考虑前面五个水文参数。
在本发明的一个实施例中,所述研究流域的水文参数的计算公式为:
A=LH;WP=2H+L;
其中,L为研究流域的长;H为研究流域的宽;WP为室外研究流域的湿周;R为室外研究流域的水力半径;J为研究流域自分水岭到出口断面主河道的平均比降;n为研究流域的糙率;
h1、h2、h3······hn为研究流域河道纵断面各转折点的河底高程与设计断面处河底高程的差值;L1、L2、L3······Ln为研究流域各转折点之间的分段河长,L′为主河长;v为室外研究流域的断面出口流速。
所述模型流域的水文参数的计算公式为:
A′=K1K2A=K1K2LH;WP′=2K2H+K1L;
其中,A′为室内模型流域的出口断面面积;WP′为室内模型流域的湿周;R′为室内模型流域的水力半径;J′为模型流域自分水岭到出口断面主河道的平均比降;n′为模型流域的糙率;K1为水平比例尺;K2为垂直比例尺;
h′1、h′2、h′3······h′n为模型流域河道纵断面各转折点的河底高程与设计断面处河底高程的差值,且h′n=K2hn;L′1、L′2、L′3······L′n为模型流域各转折点之间的分段河长,L′n=K1Ln;v′为室内模型流域的断面出口流速。
在步骤103中,构建阻力相似模型,并采用研究流域和模型流域的水文参数及阻力相似模型得到垂直比例尺。其中,所述阻力相似模型为:
其中,αu为流速比尺,αL为水平几何比尺,αL=K1;αn为糙率比尺,αR是水力半径比尺,
最后,根据模型流域和研究流域的水文参数的计算公式及阻力相似计算公式,得到垂直比例尺的计算公式:
综上所述,本方案通过研究流域的矢量信息可以快速的模型流域在X、Y轴的水平比例K1,之后通过研究流域的宽和长,就能快速地求取模型流域的垂直比例尺K2,这样采用通过垂直比例尺K2和水平比例K1对研究流域进行缩放后的模型流域求得的室内降雨强度进行人工降雨时,更能真实地反应研究区域的实际降雨情况。
Claims (7)
1.基于水动力过程相似的水文模型流域比例尺确定方法,其特征在于,包括:
获取研究流域的矢量信息,并将室内自定义坐标的矢量文件与研究流域的边界矢量文件进行匹配,得到模型流域的水平比例尺;
基于水动力相似原则,采用研究流域的长和宽及水平比例尺分别计算研究流域和模型流域的水文参数;
构建阻力相似模型,并采用研究流域和模型流域的水文参数及阻力相似模型得到垂直比例尺。
2.根据权利要求1所述的基于水动力过程相似的水文模型流域比例尺确定方法,其特征在于,所述水文参数包括出口断面面积、湿周、水力半径、自分水岭到出口断面主河道的平均比降和断面出口流速。
3.根据权利要求2所述的基于水动力过程相似的水文模型流域比例尺确定方法,其特征在于,所述研究流域的水文参数的计算公式为:
其中,L为研究流域的长;H为研究流域的宽;WP为室外研究流域的湿周;R为室外研究流域的水力半径;J为研究流域自分水岭到出口断面主河道的平均比降;n为研究流域的糙率;
h1、h2、h3……hn为研究流域河道纵断面各转折点的河底高程与设计断面处河底高程的差值;L1、L2、L3……Ln为研究流域各转折点之间的分段河长,L′为主河长;v为室外研究流域的断面出口流速。
4.根据权利要求3所述的基于水动力过程相似的水文模型流域比例尺确定方法,其特征在于,所述模型流域的水文参数的计算公式为:
其中,A′为室内模型流域的出口断面面积;WP′为室内模型流域的湿周;R′为室内模型流域的水力半径;J′为模型流域自分水岭到出口断面主河道的平均比降;n′为模型流域的糙率;K1为水平比例尺;K2为垂直比例尺;
h′1、h′2、h′3……h′n为模型流域河道纵断面各转折点的河底高程与设计断面处河底高程的差值,且h′n=K2hn;L′1、L′2、L′3……L′n为模型流域各转折点之间的分段河长,且L′n=K1Ln;v′为室内模型流域的断面出口流速。
5.根据权利要求4所述的基于水动力过程相似的水文模型流域比例尺确定方法,其特征在于,所述阻力相似模型为:
其中,αu为流速比尺,αL为水平几何比尺,αL=K1;αn为糙率比尺,αR是水力半径比尺,
6.根据权利要求5所述的基于水动力过程相似的水文模型流域比例尺确定方法,其特征在于,根据模型流域和研究流域的水文参数的计算公式及阻力相似计算公式,得到垂直比例尺的计算公式:
7.根据权利要求1-6任一所述的基于水动力过程相似的水文模型流域比例尺确定方法,其特征在于,所述研究流域的矢量信息包括研究流域形状矢量信息和雨量站位置的矢量信息;
所述研究流域的矢量信息提取方法包括:获取研究流域的地形数据和雨量站位置;将所述地形数据与雨量站位置导入ArcGis中得到研究流域形状和雨量站位置的矢量信息。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810652924.8A CN108875222B (zh) | 2018-06-22 | 2018-06-22 | 基于水动力过程相似的水文模型流域比例尺确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810652924.8A CN108875222B (zh) | 2018-06-22 | 2018-06-22 | 基于水动力过程相似的水文模型流域比例尺确定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108875222A true CN108875222A (zh) | 2018-11-23 |
CN108875222B CN108875222B (zh) | 2022-11-29 |
Family
ID=64340530
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810652924.8A Active CN108875222B (zh) | 2018-06-22 | 2018-06-22 | 基于水动力过程相似的水文模型流域比例尺确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108875222B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117034527A (zh) * | 2023-09-22 | 2023-11-10 | 水利部交通运输部国家能源局南京水利科学研究院 | 一种全域河网物理模型局部变比尺阻力相似设计方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110306022A1 (en) * | 2010-06-11 | 2011-12-15 | Adam Wcislo | Manufacturing of three-dimensional topographical scale models |
CN103473809A (zh) * | 2013-09-27 | 2013-12-25 | 武汉大学 | 基于3d打印技术的流域水文模拟方法 |
US20140324541A1 (en) * | 2013-04-30 | 2014-10-30 | International Business Machines Corporation | Using real-time online analytics to automatically generate an appropriate measurement scale |
CN104978451A (zh) * | 2015-08-12 | 2015-10-14 | 河海大学 | 一种基于ArcGIS的复杂平原河网地区水文模型处理方法 |
-
2018
- 2018-06-22 CN CN201810652924.8A patent/CN108875222B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110306022A1 (en) * | 2010-06-11 | 2011-12-15 | Adam Wcislo | Manufacturing of three-dimensional topographical scale models |
US20140324541A1 (en) * | 2013-04-30 | 2014-10-30 | International Business Machines Corporation | Using real-time online analytics to automatically generate an appropriate measurement scale |
CN103473809A (zh) * | 2013-09-27 | 2013-12-25 | 武汉大学 | 基于3d打印技术的流域水文模拟方法 |
CN104978451A (zh) * | 2015-08-12 | 2015-10-14 | 河海大学 | 一种基于ArcGIS的复杂平原河网地区水文模型处理方法 |
Non-Patent Citations (1)
Title |
---|
张妙仙等: "小流域产流汇流模型试验的应用研究", 《山西水利科技》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117034527A (zh) * | 2023-09-22 | 2023-11-10 | 水利部交通运输部国家能源局南京水利科学研究院 | 一种全域河网物理模型局部变比尺阻力相似设计方法 |
CN117034527B (zh) * | 2023-09-22 | 2024-03-19 | 水利部交通运输部国家能源局南京水利科学研究院 | 一种全域河网物理模型局部变比尺阻力相似设计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108875222B (zh) | 2022-11-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108829986A (zh) | 基于能量过程相似的水文实验模型比例尺确定方法 | |
CN106884405B (zh) | 一种无资料地区溃堤型山洪灾害分析评价方法 | |
CN109543275B (zh) | 一种城区地表径流二维数值模拟方法 | |
CN111369059B (zh) | 基于雨洪模拟耦合模型的城市内涝快速预测方法及系统 | |
CN108022047A (zh) | 一种海绵城市水文计算方法 | |
CN106529176B (zh) | 一种双核双驱洪水预报方法 | |
CN109815305A (zh) | 一种无资料地区场次洪水径流过程反演的方法 | |
CN107563139A (zh) | 一种点源突发性水污染源事故溯源贡献程度计算方法 | |
CN103886135B (zh) | 基于二维非恒定流数值模型的电力工程选址方法 | |
CN102289570A (zh) | 基于降雨-径流-洪水演进计算的洪水预报方法 | |
CN106909731B (zh) | 一种膨胀性土壤产汇流过程模拟方法 | |
Isidoro et al. | Influence of wind-driven rain on the rainfall-runoff process for urban areas: Scale model of high-rise buildings | |
Kazama et al. | Evaluation of groundwater resources in wide inundation areas of the Mekong River basin | |
CN109033589B (zh) | 一种基于swmm模型考虑lid净化作用的污染物去除率计算方法 | |
CN100580692C (zh) | 一种基于航空影像的水体与居民地变化检测方法 | |
Li et al. | Slope spectrum variation in a simulated loess watershed | |
CN108897940A (zh) | 基于矩形网格的分布式水文模型和二维水动力模型的单向耦合方法 | |
CN108875222A (zh) | 基于水动力过程相似的水文模型流域比例尺确定方法 | |
CN108763829B (zh) | 基于水动力相似的区域降雨模拟系统降雨过程设计方法 | |
CN108920571A (zh) | 基于分布式水文模型的室内水文实验模型比例尺确定方法 | |
CN108897715A (zh) | 一种基于洪水过程的水文实验模型比例尺确定方法 | |
CN109308375A (zh) | 一种基于地貌参数的流域最优流速的测算方法 | |
CN108824349A (zh) | 基于水动力过程相似的模型流域单位线的构建方法 | |
Hassan et al. | AL-Dibdiba formation basin hydrological aspects extraction using GIS techniques and quantitative morphometric analysis | |
Gautam | Flow routing with Semi-distributed hydrological model HEC-HMS in case of Narayani River Basin. |
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 |