CN108763829B - 基于水动力相似的区域降雨模拟系统降雨过程设计方法 - Google Patents

基于水动力相似的区域降雨模拟系统降雨过程设计方法 Download PDF

Info

Publication number
CN108763829B
CN108763829B CN201810652182.9A CN201810652182A CN108763829B CN 108763829 B CN108763829 B CN 108763829B CN 201810652182 A CN201810652182 A CN 201810652182A CN 108763829 B CN108763829 B CN 108763829B
Authority
CN
China
Prior art keywords
rainfall
watershed
model
sub
outdoor
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.)
Active
Application number
CN201810652182.9A
Other languages
English (en)
Other versions
CN108763829A (zh
Inventor
龚家国
杨苗
王友胜
冶运涛
翟正丽
彭玕
王英
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Institute of Water Resources and Hydropower Research
Original Assignee
China Institute of Water Resources and Hydropower Research
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China Institute of Water Resources and Hydropower Research filed Critical China Institute of Water Resources and Hydropower Research
Priority to CN201810652182.9A priority Critical patent/CN108763829B/zh
Publication of CN108763829A publication Critical patent/CN108763829A/zh
Application granted granted Critical
Publication of CN108763829B publication Critical patent/CN108763829B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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

基于水动力相似的区域降雨模拟系统降雨过程设计方法
技术领域
本发明涉及水利工程技术领域，具体涉及一种基于水动力相似的区域降雨模拟系统降雨过程设计方法。
背景技术
区域水文过程试验是当前水文试验研究的主要发展方向之一，开展室内区域水文过程试验首先需要根据试验研究目标科学设定室内区域降雨人工模拟系统的人工降雨过程。当前的人工降雨试验以点尺度的过程试验为主，其降雨过程的设定主要利用点尺度的降雨观测资料，选定典型降雨雨强和降雨历时作为人工降雨过程的设计依据。这些实验过程中一般强调降雨区域内的均匀度指标，而区域水文过程的研究对象强调指定区域内降雨、地形、径流等过程的相似性、一致性。由于两者之间的关注点不一致，当前科学设定区域人工降雨过程的方法还不多见，致使现有室内试验难以模拟野外研究区的人工降雨过程。
发明内容
针对现有技术中的上述不足，本发明提供的基于水动力相似的区域降雨模拟系统降雨过程设计方法通过构建的时段单位线能够模拟野外研究区的人工降雨。
为了达到上述发明目的，本发明采用的技术方案为：
提供一种基于地形比例尺的模型流域时段单位线，其包括：
获取研究流域的矢量信息，并将其按设定比例缩小得到的模型流域导入室内降雨区中；
按照区域降雨人工模拟系统的空间剖分结构将室内降雨区划分成若干区块，并将模型流域位于的区块标记为流域子块；
采用研究流域中雨量站测量的日降雨量及雨量站与流域子块的权重，计算流域子块的日降雨量；
采用流域子块的日降雨量、雨量雨力关系模型、降雨历时及每个流域子块中模型流域的占比，计算模型流域的室外小时降雨量；
计算模型流域在室外的净雨量R：
Figure BDA0001705122510000021
其中，I0为室外初损值；
Figure BDA0001705122510000022
为室外平均下渗率；tc为产流历时；P后为室外降雨后期不产流的雨量；D为模型流域的室外小时降雨量；
根据模型流域的室外净雨量和室内出口断面流量，计算模型流域在室外的瞬时单位线u(0，t)：
Figure BDA0001705122510000023
Figure BDA0001705122510000024
其中，n′为反应模型流域调蓄能力的参数；Γ(n′)为n′的伽马函数；k′为线性水库的调蓄系数；e为自然对数；t为时间变量；
Figure BDA0001705122510000025
分别为室外出口断面流量Q的一阶和二阶原点矩；
Figure BDA0001705122510000026
分别为室外净雨量R的一阶和二阶原点矩；
根据模型流域在室外的瞬时单位线，构建模型流域在室外净雨为z mm的时段单位线：
Figure BDA0001705122510000027
Figure BDA0001705122510000031
其中，q(Δt，t)为室外时段单位线；Δt为净雨时段；F为模型流域面积；
基于水动力过程相似原理，计算模型流域i时段末的室内出口断面流量Qi′：
Q′i＝K1K2 3/2Qi
其中，K1为模型流域的水平比例尺；K2为模型流域的垂直比例尺；Qi为模型流域i时段末的室外出口断面流量，i＝1，2，···，l为流量过程线时段数；
根据室外净雨的时段单位线，计算模型流域第j时段末的室内净雨量：
Figure BDA0001705122510000032
其中，
Figure BDA0001705122510000033
K5为q(Δt，t)与模型流域的室内时段单位线q′的换算比例尺，K5＜1；q′i-j+1为i-j+1时段的室内时段单位线，i-j+1＝1，2，···n为单位线时段数，j＝1，2，···，m为净雨时段数；
根据模型流域的模型流域第j时段末的室内净雨量R′j，计算模型流域的室内时段降雨量D′j：
Figure BDA0001705122510000034
其中，I′0为室内初损值；
Figure BDA0001705122510000035
为室内平均下渗率；tc为产流历时；P′后为室内降雨后期不产流的雨量；
根据计算的所有室内时段降雨量，生成室内模型流域正上方人工降雨装置的降雨过程控制文件。进一步地，所述流域子块的日降雨量的计算公式为：
Figure BDA0001705122510000036
其中，Hj为第j个流域子块的日降雨量；wij为第j个流域子块到第i个雨量站的权重；pi为第i个雨量站的日降雨量；c为雨量站的总数量。
进一步地，所述wij的计算公式为：
Figure BDA0001705122510000041
其中，rij为第j个流域子块到第i个雨量站的距离；b为权重指数；b＝0为平方法，b＝1为线性反比法，b＝2为RDS法。
进一步地，所述采用流域子块的日降雨量、雨量雨力关系模型、降雨历时及每个流域子块中模型流域的占比，计算模型流域的小时降雨量进一步包括：
采用流域子块的日降雨量、雨量雨力关系模型、降雨历时，计算流域子块在降雨历时内的小时降雨量为：
Figure BDA0001705122510000042
其中，Aj为第j个流域子块在降雨历时内的小时降雨量；Hj为第j个流域子块的日降雨量；Sj为第j个流域子块的雨量雨力关系模型；a和b均为雨量雨力关系参数；ε为残差；T为降雨历时；n为暴雨衰减系数；
根据流域子块的小时降雨量和每个流域子块中模型流域的占比，计算模型流域的室外小时降雨量为：
Figure BDA0001705122510000043
其中，D为模型流域的室外小时降雨量；q为每个流域子块中模型流域的占比；B为流域子块的面积；Bj为第j个流域子块中模型流域的面积。
进一步地，所述计算流域子块在降雨历时内的小时降雨量进一步包括：
根据模型流域的日降雨量，构建雨量雨力关系模型：
Sj＝aHj+b+ε；
根据计算的雨力Sj、流域子块的日降雨量Hj和率定的衰减系数n，计算一天内的降雨历时T：
Figure BDA0001705122510000051
根据雨力Sj和降雨历时T，流域子块在降雨历时内的小时降雨量为：
Figure BDA0001705122510000052
进一步地，所述
Figure BDA0001705122510000053
和
Figure BDA0001705122510000054
的计算公式分别为：
Figure BDA0001705122510000055
其中，Rl为第l时段的模型流域的室外净雨量；Ql为第l时段的模型流域的室内出口断面流量；
Figure BDA0001705122510000056
Δt为净雨时段。
进一步地，所述研究流域的矢量信息的生成方法包括：
获取研究流域的地形数据和雨量站的位置，并将地形数据和雨量站的位置导入ArcGis中生成研究流域的矢量信息。
本发明的有益效果为：提供了一种基于水动力过程的水文实验模型降雨过程的确定方法，此方法基于区域洪水过程所构建的水文模型，提供了一种科学设定室内区域水文试验中的人工模拟降雨过程的方法。该方法能较好的拟合野外研究区域，计算出比较符合野外实际情况的降雨过程，具有较好的实用性、相似性和一致性，根据水动力过程相似和研究区域的实际地形来确定水文模型的降雨过程，有较高的可靠度和准确性
附图说明
图1为基于水动力相似的区域降雨模拟系统降雨过程设计方法的流程图。
图2为模型流域导入室内降雨区的示意图。
图3为S曲线的示意图。
具体实施方式
下面对本发明的具体实施方式进行描述，以便于本技术领域的技术人员理解本发明，但应该清楚，本发明不限于具体实施方式的范围，对本技术领域的普通技术人员来讲，只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内，这些变化是显而易见的，一切利用本发明构思的发明创造均在保护之列。
参考图1，图1示出了基于水动力相似的区域降雨模拟系统降雨过程设计方法的流程图；如图1所示，该方法100包括步骤101至步骤111。
在步骤101中，获取研究流域的矢量信息，并将其按设定比例缩小得到的模型流域导入室内降雨区中，其中，区域降雨模拟系统为：有效降雨面积26m×40m，共分为110个独立控制单元(其中共分为10个大区，每个大区又分为11个小区)；每个降雨单元雨强变化范围10-200mm/h；单元内均匀度大于0.8；雨强变化频率可以做到5分钟调整一次；通过输入降雨控制文件，可以实现人工降雨。
如图2所示，40×26的区块为室内降雨区，其中位于室内降雨区中的弧形区域为模型流域。在进行缩小时，必须保证研究流域的矢量信息在x、y轴上的长宽比例不变的情况下对研究流域的矢量信息进行缩放得到模型流域，得到水平比例尺为K1，垂直比例尺为K2。
其中设定比例的选择主要跟室内降雨区的面积相关，在缩放过程中，为了更好地在室内模拟野外流域降雨过程，优选缩小后的模型流域尽量填充整个室内降雨区。
实施时，本方案优选研究流域的矢量信息的生成方法包括：获取研究流域的地形数据和雨量站的位置，并将地形数据和雨量站的位置导入ArcGis中生成研究流域的矢量信息。
在步骤102中，按照区域降雨人工模拟系统的空间剖分结构将室内降雨区划分成若干区块，并将模型流域位于的区块标记为流域子块；；实施时，本方案优选将室内降雨区共划分成110个区块，找到模型流域所在区域，并将有模型流域位于的区块标记为流域子块，并将所有的流域子块依次编上序号。
在步骤103中，采用研究流域中雨量站测量的日降雨量及雨量站与流域子块的权重，计算流域子块的日降雨量；具体地，流域子块的日降雨量的计算公式为：
Figure BDA0001705122510000071
其中，Hj为第j个流域子块的日降雨量；wij为第j个流域子块到第i个雨量站的权重；pi为第i个雨量站的日降雨量；c为雨量站的总数量。
所述wij的计算公式为：
Figure BDA0001705122510000072
其中，rij为第j个流域子块到第i个雨量站的距离；b为权重指数；b＝0为平方法，b＝1为线性反比法，b＝2为RDS法。
在步骤104中，采用流域子块的日降雨量、雨量雨力关系模型、降雨历时及每个流域子块中模型流域的占比，计算模型流域的小时降雨量；
在本发明的一个实施例中，所述采用流域子块的日降雨量、雨量雨力关系模型、降雨历时及每个流域子块中模型流域的占比，计算模型流域的小时降雨量进一步包括步骤201和步骤202：
在步骤201中，采用流域子块的日降雨量、雨量雨力关系模型、降雨历时，计算流域子块在降雨历时内的小时降雨量为：
Figure BDA0001705122510000073
其中，Aj为第j个流域子块在降雨历时内的小时降雨量；Hj为第j个流域子块的日降雨量；Sj为第j个流域子块的雨量雨力关系模型；a和b均为雨量雨力关系参数；ε为残差；T为降雨历时；n为暴雨衰减系数。
在步骤201中计算流域子块在降雨历时内的小时降雨量的方法包括：
首先通过对模型流域的日降雨量进行分析，构建雨量雨力关系模型：
Sj＝aHj+b+ε
由于雨量雨力关系模型中的流域子块的日降雨量Hj可以在103步骤中求解出来，则本方案通过构建的雨量雨力关系模型就可以直接计算出第j个流域子块的雨力Sj。
之后，根据计算的雨力Sj、流域子块的日降雨量Hj和率定的衰减系数n，计算一天内的降雨历时T：
Figure BDA0001705122510000081
之后，根据雨力Sj和降雨历时T，流域子块在降雨历时内的小时降雨量为：
Figure BDA0001705122510000082
在步骤202中，根据流域子块的小时降雨量和每个流域子块中模型流域的占比，计算模型流域的室内小时降雨量为：
Figure BDA0001705122510000083
其中，D为模型流域的室内小时降雨量；q为每个流域子块中模型流域的占比；B为流域子块的面积；Bj为第j个流域子块中模型流域的面积。
在步骤105中，计算模型流域在室外的净雨量R：
Figure BDA0001705122510000084
其中，I0为室外初损值；
Figure BDA0001705122510000085
为室外平均下渗率；tc为产流历时；P后为室外降雨后期不产流的雨量；D为模型流域的室外小时降雨量。
在步骤106中，根据模型流域的室外净雨量和室内出口断面流量，计算模型流域在室外的瞬时单位线u(0，t)：
Figure BDA0001705122510000091
Figure BDA0001705122510000092
其中，n′为反应模型流域调蓄能力的参数；Γ(n′)为n′的伽马函数；k′为线性水库的调蓄系数；e为自然对数；t为时间变量；
Figure BDA0001705122510000093
分别为室外出口断面流量Q的一阶和二阶原点矩；
Figure BDA0001705122510000094
分别为室外净雨量R的一阶和二阶原点矩。
实施时，本方案优选所述
Figure BDA0001705122510000095
和
Figure BDA0001705122510000096
的计算公式分别为：
Figure BDA0001705122510000097
其中，Rl为第l时段的模型流域的室外净雨量；Ql为第l时段的模型流域的室内出口断面流量；Δt为净雨时段。
在步骤107中，根据模型流域在室外的瞬时单位线，构建模型流域在室外净雨为zmm的时段单位线：
Figure BDA0001705122510000099
Figure BDA00017051225100000910
其中，q(Δt，t)为室外时段单位线；Δt为净雨时段；F为模型流域面积。
在本发明的一个实施例中，根据模型流域在室外的瞬时单位线，构建模型流域在室外净雨为z mm的时段单位线进一步包括：
首先对得到的模型流域在室外的瞬时单位线进行求导，得到S曲线：
Figure BDA0001705122510000101
由于n′、k′已知，以不同的t代入上式积分，就可以得到如图3所示的S曲线。以t＝0为起点的S(t)曲线向后平移一个Δt时段，即可得到S(t-Δt)曲线，两条曲线的纵坐标差：
U(Δt，t)＝S(t)-S(t-Δt)
上式U(Δt，t)即为时段为Δt的无因次时段单位线，将无因次时段单位线换算成时段为Δt，净雨为zmm的时段单位线为：
Figure BDA0001705122510000102
其中，q(Δt，t)为室外时段单位线；Δt为净雨时段；F为模型流域面积。
在步骤108中，基于水动力过程相似原理，计算模型流域i时段末的室内出口断面流量Q′i：
Qi′＝K1K2 3/2Qi
其中，K1为模型流域的水平比例尺；K2为模型流域的垂直比例尺；Qi为模型流域i时段末的室外出口断面流量，i＝1，2，···，l为流量过程线时段数；
在步骤109中，根据室外净雨的时段单位线，计算模型流域第j时段末的室内净雨量：
Figure BDA0001705122510000103
其中，
Figure BDA0001705122510000104
K5为q(Δt，t)与模型流域的室内时段单位线q′的换算比例尺，K5＜1；q′i-j+1为i-j+1时段的室内时段单位线，i-j+1＝1，2，···n为单位线时段数，j＝1，2，···，m为净雨时段数；
当i＝1的时候，j＝1；i＝2的时候，是j＝1和j＝2时两个结果的和；i＝3时，就是j＝1、j＝2、j＝3这三个的和。
在步骤110中，根据模型流域的模型流域第j时段末的室内净雨量R′j，计算模型流域的室内时段降雨量D′j：
Figure BDA0001705122510000111
其中，I′0为室内初损值；
Figure BDA0001705122510000112
为室内平均下渗率；tc为产流历时；P′后为室内降雨后期不产流的雨量；
在步骤111中，根据计算的所有室内时段降雨量，生成室内模型流域正上方人工降雨装置的降雨过程控制文件。
其中的控制文件为每个流域子块的人工降雨量，其具体获取方法为：按照分区情况将降雨(所有室内时段降雨量)分配到每一个区块上，利用面积权重来分配得到每个区域的降雨，并采用每一个降雨插值分区的室内降雨量实施人工降雨。
下面采用Nash效率、相关系数和相对误差对本方案构建的时段单位线的效果进行说明：
将假定的室内降雨量P′按照分区情况将降雨分配到每一个区块上，可以利用面积权重来分配得到每个区域的降雨量P′n。采用假定的每一个降雨插值分区的内室降雨量P′实施人工降雨，并对降雨过程中的雨强和指定观测点的流量过程进行监测，记录完整的数据。
采用以下计算公式计算模型流域的理论流量：
Figure BDA0001705122510000113
其中，Q′理i为i时段末的室内理论流量，i＝1，2，···，l为流量过程线时段数；q′i-j+1为i-j+1时段的室内时段单位线，i-j+1＝1，2，···n为单位线时段数；R′j为模型流域第j时段末的室内净雨量，j＝1，2，···，m为净雨时段数；
当i＝1的时候，j＝1；i＝2的时候，是j＝1和j＝2时两个结果的和；i＝3时，就是j＝1、j＝2、j＝3这三个的和。
下面采用Nash效率、相关系数和相对误差对Q′理i与实验得到实测流量过程线Q′n进行校验：
Nash效率的计算公式为：
Figure BDA0001705122510000121
其中，Qi为采用本方案的方法计算得到的降雨量，Qi等于上面的Q′理i，qi为室内实测降雨量，
Figure BDA0001705122510000122
为室内实测降雨量的平均值。
相关系数的计算公式为：
Figure BDA0001705122510000123
rxy为相关系数；n为系列的样本数；X、Y分别代表实测系列和模拟系列的数值。
相对误差的计算公式为：
Figure BDA0001705122510000124
其中，Dv为模相对误差(％)；F0为室内实测降雨量的均值；R为采用本方案的方法计算得到的降雨量均值。
将在模拟人工降雨过程中记录的雨强和指定观测点的流量带入Nash效率、相关系数和相对误差的计算公式中，可以得到Nash效率接近于1，相关系数rxy位于0.8-1.0之间；相对误差绝对值接近于零。
通过上述对对Q′理与实验得到实测流量过程线Q′n的校验可以看出，本方案构建的时段单位线能够很好的模拟模拟水沙过程及野外研究流域的人工降雨。

Claims (6)

1.一种基于水动力相似的区域降雨模拟系统降雨过程设计方法，其特征在于，包括：
获取研究流域的矢量信息，并将其按设定比例缩小得到的模型流域导入室内降雨区中；
按照区域降雨人工模拟系统的空间剖分结构将室内降雨区划分成若干区块，并将模型流域位于的区块标记为流域子块；
采用研究流域中雨量站测量的日降雨量及雨量站与流域子块的权重，计算流域子块的日降雨量；
采用流域子块的日降雨量、雨量雨力关系模型、降雨历时及每个流域子块中模型流域的占比，计算模型流域的室外小时降雨量；
计算模型流域在室外的净雨量R：
Figure FDA0002328456020000011
其中，I0为室外初损值；
Figure FDA0002328456020000012
为室外平均下渗率；tc为产流历时；P后为室外降雨后期不产流的雨量；D为模型流域的室外小时降雨量；
根据模型流域的室外净雨量和室外出口断面流量，计算模型流域在室外的瞬时单位线u(0，t)：
Figure FDA0002328456020000013
Figure FDA0002328456020000014
其中，n′为反应模型流域调蓄能力的参数；Γ(n′)为n′的伽马函数；k′为线性水库的调蓄系数；e为自然对数；t为时间变量；MQ (1)、MQ (2)分别为室外出口断面流量Q的一阶和二阶原点矩；MR (1)、MR (2)分别为室外净雨量R的一阶和二阶原点矩；
根据模型流域在室外的瞬时单位线，构建模型流域在室外净雨为zmm的时段单位线：
Figure FDA0002328456020000021
Figure FDA0002328456020000022
其中，q(Δt，t)为室外时段单位线；Δt为净雨时段；F为模型流域面积；
基于水动力过程相似原理，计算模型流域i时段末的室内出口断面流量Q′i：
Q′i＝K1K2 3/2Qi
其中，K1为模型流域的水平比例尺；K2为模型流域的垂直比例尺；Qi为模型流域i时段末的室外出口断面流量，i＝1，2，…，l为流量过程线时段数；
根据室外净雨的时段单位线，计算模型流域第j时段末的室内净雨量：
Figure FDA0002328456020000023
其中，
Figure FDA0002328456020000024
K5为q(Δt，t)与模型流域的室内时段单位线q′的换算比例尺，K5＜1；q′i-j+1为i-j+1时段的室内时段单位线，i-j+1＝1，2，…n为单位线时段数，j＝1，2，…，m为净雨时段数；
根据模型流域的模型流域第j时段末的室内净雨量R′j，计算模型流域的室内时段降雨量D′j：
Figure FDA0002328456020000025
其中，I′0为室内初损值；
Figure FDA0002328456020000026
为室内平均下渗率；tc为产流历时；P′后为室内降雨后期不产流的雨量；
根据计算的所有室内时段降雨量，生成室内模型流域正上方人工降雨装置的降雨过程控制文件；
所述研究流域的矢量信息的生成方法包括：
获取研究流域的地形数据和雨量站的位置，并将地形数据和雨量站的位置导入ArcGis中生成研究流域的矢量信息。
2.根据权利要求1所述的基于水动力相似的区域降雨模拟系统降雨过程设计方法，其特征在于，所述流域子块的日降雨量的计算公式为：
Figure FDA0002328456020000031
其中，Hj为第j个流域子块的日降雨量；wij为第j个流域子块到第i个雨量站的权重；pi为第i个雨量站的日降雨量；c为雨量站的总数量。
3.根据权利要求2所述的基于水动力相似的区域降雨模拟系统降雨过程设计方法，其特征在于，所述wij的计算公式为：
Figure FDA0002328456020000032
其中，rij为第j个流域子块到第i个雨量站的距离；b为权重指数；b＝0为平方法，b＝1为线性反比法，b＝2为RDS法。
4.根据权利要求1-3任一所述的基于水动力相似的区域降雨模拟系统降雨过程设计方法，其特征在于，所述采用流域子块的日降雨量、雨量雨力关系模型、降雨历时及每个流域子块中模型流域的占比，计算模型流域的小时降雨量进一步包括：
采用流域子块的日降雨量、雨量雨力关系模型、降雨历时，计算流域子块在降雨历时内的小时降雨量为：
Figure FDA0002328456020000033
其中，Aj为第j个流域子块在降雨历时内的小时降雨量；Hj为第j个流域子块的日降雨量；Sj为第j个流域子块的雨量雨力关系模型；a和b均为雨量雨力关系参数；ε为残差；T为降雨历时；n为暴雨衰减系数；
根据流域子块的小时降雨量和每个流域子块中模型流域的占比，计算模型流域的室外小时降雨量为：
Figure FDA0002328456020000041
其中，D为模型流域的室外小时降雨量；q为每个流域子块中模型流域的占比；B为流域子块的面积；Bj为第j个流域子块中模型流域的面积。
5.根据权利要求4所述的基于水动力相似的区域降雨模拟系统降雨过程设计方法，其特征在于，所述计算流域子块在降雨历时内的小时降雨量进一步包括：
根据模型流域的日降雨量，构建雨量雨力关系模型：
Sj＝aHj+b+ε；
根据计算的雨力Sj、流域子块的日降雨量Hj和率定的衰减系数n，计算一天内的降雨历时T：
Figure FDA0002328456020000042
根据雨力Sj和降雨历时T，流域子块在降雨历时内的小时降雨量为：
Figure FDA0002328456020000043
6.根据权利要求1所述的基于水动力相似的区域降雨模拟系统降雨过程设计方法，其特征在于，所述MQ (1)、MQ (2)、MR (1)和MR (2)的计算公式分别为：
Figure FDA0002328456020000044
Figure FDA0002328456020000045
其中，Rl为第l时段的模型流域的室外净雨量；Ql为第l时段的模型流域的室内出口断面流量；
Figure FDA0002328456020000046
Δt为净雨时段。
CN201810652182.9A 2018-06-22 2018-06-22 基于水动力相似的区域降雨模拟系统降雨过程设计方法 Active CN108763829B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810652182.9A CN108763829B (zh) 2018-06-22 2018-06-22 基于水动力相似的区域降雨模拟系统降雨过程设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810652182.9A CN108763829B (zh) 2018-06-22 2018-06-22 基于水动力相似的区域降雨模拟系统降雨过程设计方法

Publications (2)

Publication Number Publication Date
CN108763829A CN108763829A (zh) 2018-11-06
CN108763829B true CN108763829B (zh) 2020-04-14

Family

ID=63976643

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810652182.9A Active CN108763829B (zh) 2018-06-22 2018-06-22 基于水动力相似的区域降雨模拟系统降雨过程设计方法

Country Status (1)

Country Link
CN (1) CN108763829B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110659823B (zh) * 2019-09-21 2023-02-03 四川大学工程设计研究院有限公司 相似流域分析方法、模型、系统和计算机存储介质
CN111008259B (zh) * 2019-12-06 2020-08-11 中国水利水电科学研究院 一种流域降雨相似性搜索方法
CN115688622B (zh) * 2022-10-26 2023-06-13 中国长江三峡集团有限公司 一种库区区间来沙量的计算方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102930357A (zh) * 2012-11-20 2013-02-13 中铁第四勘察设计院集团有限公司 岩溶隧道暗河暴雨涌水洪峰峰值及峰值时间的预测方法
CN103886151A (zh) * 2014-03-27 2014-06-25 北京工业大学 不同暴雨情景下城市下凹式立交桥区简化水力模型构建方法
CN104281780A (zh) * 2014-10-11 2015-01-14 水利部交通运输部国家能源局南京水利科学研究院 线性水库滞留汇流及嵌套流域(多子流域)汇流方法
CN107730151A (zh) * 2017-11-21 2018-02-23 中国水利水电科学研究院 一种基于概念性水文模型的流域设计洪水推求方法
CN108018823A (zh) * 2017-12-15 2018-05-11 河海大学 基于瞬时单位线的流域地下水平均响应时间估计方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9772609B2 (en) * 2013-05-24 2017-09-26 Optirtc, Inc. Optimized hydromodification management with active stormwater controls

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102930357A (zh) * 2012-11-20 2013-02-13 中铁第四勘察设计院集团有限公司 岩溶隧道暗河暴雨涌水洪峰峰值及峰值时间的预测方法
CN103886151A (zh) * 2014-03-27 2014-06-25 北京工业大学 不同暴雨情景下城市下凹式立交桥区简化水力模型构建方法
CN104281780A (zh) * 2014-10-11 2015-01-14 水利部交通运输部国家能源局南京水利科学研究院 线性水库滞留汇流及嵌套流域(多子流域)汇流方法
CN107730151A (zh) * 2017-11-21 2018-02-23 中国水利水电科学研究院 一种基于概念性水文模型的流域设计洪水推求方法
CN108018823A (zh) * 2017-12-15 2018-05-11 河海大学 基于瞬时单位线的流域地下水平均响应时间估计方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
《基于流量_面积比值的小流域设计洪水计算方法对比研究》;叶永东;《中国水能及电气化》;20180531;第57-60页 *
《柳州市中小河流洪水预警预报存在的问题与解决方案探讨》;姜文;《广西水利水电》;20171030;第57-60页 *
《黄河呼和浩特段水动力与降雨径流耦合模型的构建》;冯利忠等;《中国农业大学学报》;20160715;第113-120页 *

Also Published As

Publication number Publication date
CN108763829A (zh) 2018-11-06

Similar Documents

Publication Publication Date Title
CN102289570B (zh) 基于降雨-径流-洪水演进计算的洪水预报方法
Gebre Application of the HEC-HMS model for runoff simulation of Upper Blue Nile River Basin
CN108022047A (zh) 一种海绵城市水文计算方法
CN108763829B (zh) 基于水动力相似的区域降雨模拟系统降雨过程设计方法
CN114997541B (zh) 一种基于数字孪生技术的城市内涝预测方法及预警平台
CN108984823B (zh) 一种合流制溢流调蓄池规模的确定方法
Li et al. Modeling the Quality and Quantity of Runoff in a Highly Urbanized Catchment Using Storm Water Management Model.
CN104281780A (zh) 线性水库滞留汇流及嵌套流域(多子流域)汇流方法
CN103886135A (zh) 基于二维非恒定流数值模型的电力工程选址方法
CN108829990B (zh) 一种区域降雨人工模拟系统的降雨过程设计方法
Yener et al. Modeling studies with HEC-HMS and runoff scenarios in Yuvacik Basin, Turkiye
Koutsovili et al. Numerical assessment of climate change impact on the hydrological regime of a small Mediterranean river, Lesvos Island, Greece
Bahrami et al. Estimation of flood hydrographs in the ungauged mountainous watershed with Gray synthetic unit hydrograph model
Gharib et al. Evaluation of ModClark model for simulating rainfall-runoff in Tangrah watershed, Iran.
CN114186424A (zh) 一种气候变化对水源涵养量影响的分析方法
Gyori et al. Unit hydrograph generation for ungauged subwatersheds. Case study: the Monoroştia river, Arad County, Romania
CN108824349B (zh) 基于水动力过程相似的模型流域单位线的构建方法
Adib et al. Comparison between GcIUH-Clark, GIUH-Nash, Clark-IUH, and Nash-IUH models
Chiew et al. Estimating groundwater recharge using a surface watershed modelling approach
CN106980764B (zh) 一种通用的直接利用遥感蒸发的水文模拟方法
Fathy et al. Runoff hydrograph modeling for arid regions: case study—Wadi Sudr-Sinai
CN114943490A (zh) 一种渠道衬砌对灌区生态耗水影响的评价方法及装置
Churiulin et al. Research of opportunities of combined use of the runoff formation ECOMAG model and mesoscale atmosphere circulation COSMO-Ru model (on the example of floods on the Sukhona River at the Velikiy Ustyug)
Koivusalo et al. Modelling runoff generation in a forested catchment in southern Finland
CN108951531B (zh) 基于地形比例尺的模型流域时段单位线的构建方法

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