一种保持典型洪水形态的设计洪水过程解析推求方法及系统的制作方法
【技术领域】
[0001] 本发明涉及水利工程领域，特别是一种保持典型洪水形态的设计洪水过程解析推 求方法。
【背景技术】
[0002] 设计洪水过程线是防洪规划和防洪工程设计的基本依据。现行拟定设计洪水过程 线方法的程序是先进行洪水峰、量频率计算，分析设计流域的洪峰流量及时段洪量的分布 函数（或频率曲线），分别求得符合设计标准的设计洪峰流量值及设计时段的洪量值；再根 据洪水特性选择典型洪水，并考虑水工设计要求，选取其中一项或几项对防洪后果影响最 大的特征，以它们为控制，对典型洪水过程线进行放大，从而得到满足指定设计标准的设计 洪水过程。常用的放大方法有同倍比放大法和同频率放大法。
[0003] 然而，同频率放大法常用于峰量关系不够好、洪峰形状差别大的河流。在运用此法 进行典型洪水过程线放大时，缺陷之一是各时段按同倍比放大，改变了典型洪水过程线的 形态，有时甚至可能将矮胖型的洪水过程放大成尖瘦型。缺陷之二是由于在不同历时衔接 的地方放大倍比不一致，放大后在交界处产生不连续现象，此时需要人工修正不连续点，而 为了要满足时段洪量的设计要求，常常需要反复多次修改试算，工作量和任意性大。
【发明内容】
[0004] 发明目的：提供一种保持典型洪水形态的设计洪水过程解析推求方法，以解决现 有技术存在的上述问题。
[0005] 技术方案：一种保持典型洪水形态的设计洪水过程解析推求方法，包括以下步 骤：
[0006] S1、根据典型洪水和设计标准的各参数计算不同时段的抬升值；
[0007] S2、采用不连续修正策略，根据不同时段的抬升值计算交界处的抬升值；
[0008] S3、根据各时序的抬升值放大典型洪水过程线。
[0009] 所述不连续修正策略为：洪峰点以设计值为准，且在设计时段的不连续点处，取不 连续点两边抬升值的均值作为不连续点处的抬升值；将由于不连续修正而导致的水量均匀 分摊到不受影响的时段中。
[0010] 基于相似性原理，各设计时段平移的初始抬升值为：
【主权项】
1. 一种保持典型洪水形态的设计洪水过程解析推求方法，其特征在于，包括w下步 骤： 51、 根据典型洪水和设计标准的各参数计算不同时段的抬升值； 52、 采用不连续修正策略，根据不同时段的抬升值计算交界处的抬升值； 53、 根据各时序的抬升值放大典型洪水过程线。
2. 如权利要求1所述的保持典型洪水形态的设计洪水过程解析推求方法，其特征在 于，所述不连续修正策略为：洪峰点W设计值为准，且在设计时段的不连续点处，取不连续 点两边抬升值的均值作为不连续点处的抬升值；将由于不连续修正而导致的水量均匀分摊 到不受影响的时段中。
3. 如权利要求2所述的保持典型洪水形态的设计洪水过程解析推求方法，其特征在 于，所述步骤S1具体为： 基于相似性原理，各设计时段平移的初始抬升值为：
式中：*。为洪峰时刻所对应的时序;胖1。，胖3。，胖^分别为典型洪水过程对应的最大1日、 3日、7日洪量；9。。，￥1。，￥3。，￥7。分别为设计标准对应的洪峰流量和最大1日、3日、7日洪量； tk，ti。分别为最大1日洪量的开始和结束时序；t3,，t3。分别为最大3日洪量的开始和结束时 序；分别为最大7日洪量的开始和结束时序；Qd(0为典型洪水过程洪峰点的数值； AQ(tm)为洪峰点的抬升值；A幻日片）为最大1日洪量所对应的时段抬升值；么0如。为最大3 日洪量除去最大1日所对应的时段抬升值；A货(/)为最大7日洪量除去最大3日所对应的 时段抬升值；at为时段长。 采用下式衡量放大后的洪水过程线与典型洪水过程线形态差异的程度：
式中：f为形态差异性指标；Qp(i-l)，Qp(i)分别为设计洪水过程线第i-1时段的初、末 的流量吼(i-l)，Qi(i)分别为典型洪水过程线第i-1时段的初、末的流量；num为洪水过程 离散点的个数。 其约束条件为放大后的洪水过程线的洪峰及各种历时的洪量分别等于设计洪峰和设 计洪量： max(Qp(i), i G [I'num]) = Q哑；
4. 如权利要求3所述的保持典型洪水形态的设计洪水过程解析推求方法，其特征在 于，所述方法还包括： 在各区间重新计算抬升值， 步骤S1. 1、在[tk+l，tm-l] U [tm+l，tk-l]区间计算新的抬升值AQi(t)，依据不连续 修正策略，tw tie 处的抬升值为 0,(',,)-g',化、）二()"咕,)一g心|,.)= A 幻，A Q, (t) 为最大3日洪量除去最大1日所对应的时段的新的抬升值， Wip-Wid二{ A Q 1 (t) X (tie-tis_4) + [Qp (tis) -Qd (tis) + A Qi (t) ] /2+ 阳p (t J -Qd (t J + A Qi (t) ] /化 A Q (tm) + A Qi (t)} A T ； 步骤S1.2、在[t3,+l，tk-l] U [tk+l，t3e-l]区间计算新的抬升值AQ2(t)，依据不连续 修正策略，tw tse处的抬升值为(_),，a;.) - a 怡、）=(_),,,也,)-(2,也,）=聯，A Q3 (t) 为最大7日洪量除去最大3日所对应的时段的新的抬升值； (胖3。-胖3<1) - (Wlp-Wld) = { A (t) X (tis-tss+tse-tie-A) + [Qp (tis) ~Qd (tis) + A Q2 (t) ] /2+ [Qp (tig) -Qj (tig) + A Q2 (t) ] /2+ [Qp 咕日)-Qd 咕日)+ A Q2 (t) ] /化[Qp (tse) -Qd 咕6) + A Q2 (t) ] /w A T ; 步骤S1.3、在U [tgg+Lt,。]区间计算新的抬升值AQ3(t)， (胖冲-胖川)-(胖3口-胖3<1)二{ A Q3 (t) X (tss-th+tye-tse-S) + 阳p 咕曰)-Qd 咕曰)+ A Q3 (t) ] /2+ 阳p 咕6) -Qd 咕e) + A Q3 (t) ] /2} A T。 步骤S2、由不同时段的抬升值计算交界处的抬升值： tm处 Qp(tm) = Qmp，tss，tis，tie, tse处的计算公式如下：
十;其中 X = A Qi (t)，y = A (t)，z = A Qs (1：)。
5. 如权利要求4所述的保持典型洪水形态的设计洪水过程解析推求方法，其特征在 于，所述步骤S3具体为： 根据各时序的抬升值放大典型洪水过程线：
A Qi (t) = X t G (ti曰，tm) U (tm，t J ; A (t) = y t e (tss, tis) U (tie, tse); AQsCt) = Z t e 山曰，ts曰）U (tse'tj。
6. -种保持典型洪水形态的设计洪水过程解析推求系统，其特征在于，包括： 第一模块，用于根据典型洪水和设计标准的各参数计算不同时段的抬升值； 第二模块，用于采用不连续修正策略，根据不同时段的抬升值计算交界处的抬升值； 第H模块，用于根据各时序的抬升值放大典型洪水过程线。
【专利摘要】本发明公开了一种保持典型洪水形态的设计洪水过程解析推求方法及系统，其中所述方法包括以下步骤：根据典型洪水和设计标准的各参数计算不同时段的抬升值；采用不连续修正策略，根据不同时段的抬升值计算交界处的抬升值；根据各时序的抬升值放大典型洪水过程线。本发明提出的方法与传统同频率放大法相比设计洪水与典型洪水形态变形更小，多时段不连续性处理更规范，与现有的基于优化算法的放大方法相比，解析法计算更为简单，效率更高。
【IPC分类】G06F17-50
【公开号】CN104598676
【申请号】CN201510008167
【发明人】钟平安, 张阳, 张宇, 陈娟, 吴业楠
【申请人】河海大学
【公开日】2015年5月6日
【申请日】2015年1月7日