CN106909797B - 一种基于马斯京根法反推区间入流的方法 - Google Patents

一种基于马斯京根法反推区间入流的方法 Download PDF

Info

Publication number
CN106909797B
CN106909797B CN201710134508.4A CN201710134508A CN106909797B CN 106909797 B CN106909797 B CN 106909797B CN 201710134508 A CN201710134508 A CN 201710134508A CN 106909797 B CN106909797 B CN 106909797B
Authority
CN
China
Prior art keywords
downstream
flow
sequence
value
website
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
CN201710134508.4A
Other languages
English (en)
Other versions
CN106909797A (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 CN201710134508.4A priority Critical patent/CN106909797B/zh
Publication of CN106909797A publication Critical patent/CN106909797A/zh
Application granted granted Critical
Publication of CN106909797B publication Critical patent/CN106909797B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
    • G16Z99/00—Subject matter not provided for in other main groups of this subclass

Landscapes

  • Measuring Volume Flow (AREA)

Abstract

本发明公开了一种基于马斯京根法反推区间入流的方法，涉及水文资料整编领域。所述方法：获得上游站点流程序列和下游站点流量序列，获取上、下游站点实测流量过程，同时，使用马斯京根法计算得到的下游站点估计流量过程，然后判断得到上、下游站点实测流量极值序列和下游站点估计流量极值序列；设计目标函数，当目标函数f达到最小值fmin时，所对应的qsim即为下游站点出流估计流量过程，将同一时刻下游站点实测流量过程与下游站点出流估计流量过程相减，得到该目标河段估计区间入流过程。本发明实现了在掌握资料较少的情况下通过较简单的操作仍旧可以得到较高精度的区间入流过程及马斯京根系数。

Description

一种基于马斯京根法反推区间入流的方法
技术领域
本发明涉及水文资料整编领域，尤其涉及一种基于马斯京根法反推区间入流的方法。
背景技术
当前为研究水库调度、水资源管理等问题，河道演进及区间入流资料的收集/推求是必不可少的步骤之一。现广泛使用的河道演进方法可分为水动力学方法和水文学方法。水动力学方法计算精度较高，物理意义明确，但所需资料较多，计算复杂，操作性较差，且求解时间较长。而，作为众多研究者采用的水文学河道演进方法----马斯京根法，计算过程简便，所需资料较少，且精度能满足一般科研、工程需求。但，马斯京根法建立在河段入流水量、出流水量平衡的假定上提出，多数情况下不考虑河道上、下游断面之间的区间入流影响，对有支流汇入的情况，往往采取“先演后合”或“先合后演”方法，不过“先演后合”方法或“先合后演”方法使用的前提是需要收集到支流的流量过程资料。
然而，在一些小流域，水文站点建设并不完善，一些支流往往缺少观测资料；或者在山区流域，山洪暴发时区间入流量等于或超出河道径流流量。此时，直接应用马斯京根方法进行河道演算得到的结果的准确性低，需要反推区间入流。
发明内容
本发明的目的在于提供一种基于马斯京根法反推区间入流的方法，从而解决现有技术中存在的前述问题。
为了实现上述目的，本发明1.一种基于马斯京根法反推区间入流的方法，其特征在于，所述方法包括：
S1，资料收集与整编
收集所计算目标河段上游站点、下游站点同时段流量资料，并分别将上游站点资料、下游站点资料整编为相同长度、相同计算时间步长的流量序列，分别得到上游站点流量序列和下游站点流量序列；
S2，洪峰自动识别
在上游站点流量序列和下游站点流量序列的基础上，获取上、下游站点实测流量过程，同时，使用马斯京根法计算得到的下游站点估计流量过程；根据预先设置的识别条件分别识别出两个流量过程中流量出现的各个极大值或各个极小值；所述识别条件包括：
判断第t时刻的流量qt是否为极大值的识别条件为公式(1)：
qt-1<qt且qt+1<qt (1)；
判断第t时刻的流量qt是否为极小的识别条件为公式(2)：
qt-1>qt且qt+1>qt (2)；
其中，qt-1、qt、qt+1分别为同一个流量过程中第t-1、t、t+1时刻的流量值；
将识别得到的任意一个流量过程的各个极大值按流量大小排序；或将识别得到的任意一个流量过程的各个极小值按流量大小排序，得到上、下游站点实测流量极值序列和下游站点估计流量极值序列；
S3，设计目标函数
设定马斯京根法参数优化的目标函数f为公式(3)：
f＝α×f1(Tobs,Tsim)+β×f2(qobs,qsim) (3)；
其中，f1为考虑洪峰峰现时间的分式，f2为考虑洪水水量拟合度的分式，α为f1的系数，β为f2的系数，Tobs为下游站点实测流量过程极大值或极小值出现的计算量序列，Tsim为下游站点估计流量过程极大值或极小值的出现的计算量序列，qobs为下游站点实测流量极值序列中的实际流量，qsim为下游站点估计流量极值过程中的估计流量；
其中，f1通过公式(4)计算：
公式(4)中，Tobs,j为下游站点实测流量极值序列中第j标记点出现的计算量，Tsim,j为下游站点估计流量极值序列中第j标记点出现的计算量，所述标记点为极大值或极小值，所述计算量为出现时间或出现序数；M为人为选定的下游站点实测流量极值序列中流量为极大值或极小值时洪峰出现的个数，M≥1；
f2通过公式(5)计算：
公式(5)中，qobs,i为下游站点实测流量极值序列中i时刻的实际流量，qsim,i为下游站点估计流量极值过程中i时刻的估计流量，表示下游站点实际流量过程中的实际流量均值，N为下游站点实测流量极值序列或下游站点估计流量极值过程中的时刻总数；
S4，求目标函数f的最小值fmin；
当目标函数f达到最小值fmin时，所对应的qsim即为下游站点出流估计流量过程，得到该qsim的马斯京根法参数即为目标河段马斯京根法参数；
S5，区间入流反推
将同一时刻的下游站点实测流量过程与下游站点出流估计流量过程相减，得到该目标河段估计区间入流过程。
优选地，步骤S3中，系数α和系数β的限定条件为：0<β<α<1，α+β＝1。
优选地，步骤S1中，所述流量资料为上游站点和下游站点在某个时间段内不同时刻的流量值，在所述流量资料中包括流量值和所述流量值对应的时间。
优选地，所述流量序列按照时间的先后顺序排列。
优选地，上游站点流量序列作为所述方法的输入变量，下游站点流量序列作为步骤S3中所述目标函数f的输入变量之一。
优选地，利用遗传算法或粒子群算法求目标函数f的最小值fmin。
在公式(4)中需要将下游站点估计流量序列和下游站点实测流量序列中的排序位置相同的极大值点或极小值点进行计算，例如：将两个流量序列中排序为1的极大值点相减，将两个流量序列中排序为2的极大值相减等等。
本发明的有益效果是：
本发明所述基于马斯京根法反推区间入流的方法，在计算马斯京根参数时，将计算流量过程与实测流量过程的峰现时间差加入目标函数当中，并取较大权重值。本发明所述方法，在掌握资料较少的情况下通过较简单的操作仍旧可以得到较高精度的区间入流过程及马斯京根系数。
附图说明
图1为2015年9月18日8时至2015年9月25日23时西江龙滩水电站出库流量过程和岩滩水电站入库流量过程；
图2为经过遗传算法优化的马斯京根法演进得到的岩滩水电站入库估计流量过程与实际流量过程；
图3为计算得到的西江龙滩水电站、岩滩水电站区间入流过程。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施方式仅仅用以解释本发明，并不用于限定本发明。
实施例
本实施例以广西西江流域龙滩水电站至岩滩水电站间的水量演进为实例，以表现本发明达到的效果。
参照图1，龙滩、岩滩水电站均位于广西壮族自治区西北部，同属红水河梯级电站群。龙滩水电站位于广西天峨县城上游15公里处，设计蓄水位375米，装机容量630万千瓦，属多年调节水库。岩滩水电站位于大化瑶族自治县境内，正常蓄水位223米，装机容量121万千瓦，属季调节水库。龙滩至岩滩区间干流长166公里，区间总年际8080平方公里，占坝址控制面积的7.6％。实施例以起止时间为2015年9月18日8时至2015年9月25日23时，时间步长为1小时的龙滩水电站出库资料和岩滩水电站入库资料为基础，计算龙滩、岩滩水电站间红水河干流的马斯京根演进参数，并计算其区间入流过程。本实施例基于马斯京根法反推区间入流的方法的步骤如下：
步骤一：收集龙滩、岩滩水电站2015年9月18日8时至2015年9月25日23时的逐小时入库流量过程资料和出库流量过程资料，见图2。
步骤二：由于该时段两水库流量过程极大值较复杂，因此选择识别极小值点。
步骤三：设计目标函数，初选α、β的值分别为0.7和0.3。
步骤四：使用经典遗传算法，计算龙滩至岩滩河段马斯京根法参数及岩滩水库入库估计流量过程。
步骤五：将下游岩滩水库入库流量实测过程和步骤四得到的下游岩滩的估计流量过程相减，得到龙滩水电站至岩滩水电站之间的区间入流过程。得到的结果图见图3和图4。
通过采用本发明公开的上述技术方案，得到了如下有益的效果：
本发明所述基于马斯京根法反推区间入流的方法，在计算马斯京根参数时，将计算流量过程与实测流量过程的峰现时间差加入目标函数当中，并取较大权重值。本发明所述方法，在掌握资料较少的情况下通过较简单的操作仍旧可以得到较高精度的区间入流过程及马斯京根系数。
以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视本发明的保护范围。

Claims (6)

1.一种基于马斯京根法反推区间入流的方法，其特征在于，所述方法包括：
S1，资料收集与整编
收集所计算目标河段上游站点、下游站点同时段流量资料，并分别将上游站点资料、下游站点资料整编为相同长度、相同计算时间步长的流量序列，分别得到上游站点流量序列和下游站点流量序列；
S2，洪峰自动识别
在上游站点流量序列和下游站点流量序列的基础上，获取上、下游站点实测流量过程，同时，使用马斯京根法计算得到的下游站点估计流量过程；根据预先设置的识别条件分别识别出两个流量过程中流量出现的各个极大值或各个极小值；所述识别条件包括：
判断第t时刻的流量qt是否为极大值的识别条件为公式(1)：
qt-1<qt且qt+1<qt (1)；
判断第t时刻的流量qt是否为极小的识别条件为公式(2)：
qt-1>qt且qt+1>qt (2)；
其中，qt-1、qt、qt+1分别为同一个流量过程中第t-1、t、t+1时刻的流量值；
将识别得到的任意一个流量过程的各个极大值按流量大小排序；或将识别得到的任意一个流量过程的各个极小值按流量大小排序，得到上、下游站点实测流量极值序列和下游站点估计流量极值序列；
S3，设计目标函数
设定马斯京根法参数优化的目标函数f为公式(3)：
f＝α×f1(Tobs,Tsim)+β×f2(qobs,qsim) (3)；
其中，f1为考虑洪峰峰现时间的分式，f2为考虑洪水水量拟合度的分式，α为f1的系数，β为f2的系数，Tobs为下游站点实测流量过程极大值或极小值出现的计算量序列，Tsim为下游站点估计流量过程极大值或极小值出现的计算量序列，qobs为下游站点实测流量极值序列中的实际流量，qsim为下游站点估计流量极值过程中的估计流量；
其中，f1通过公式(4)计算：
公式(4)中，Tobs,j为下游站点实测流量极值序列中第j标记点出现的计算量，Tsim,j为下游站点估计流量极值序列中第j标记点出现的计算量，所述标记点为极大值或极小值，所述计算量为出现时间或出现序数；M为人为选定的下游站点实测流量极值序列中流量为极大值或极小值时洪峰出现的个数，M≥1；
f2通过公式(5)计算：
公式(5)中，qobs,i为下游站点实测流量极值序列中i时刻的实际流量，qsim,i为下游站点估计流量极值过程中i时刻的估计流量，表示下游站点实际流量过程中的实际流量均值，N为下游站点实测流量极值序列或下游站点估计流量极值过程中的时刻总数；
S4，求目标函数f的最小值fmin；
当目标函数f达到最小值fmin时，所对应的qsim即为下游站点出流估计流量过程，得到该qsim的马斯京根法参数即为目标河段马斯京根法参数；
S5，区间入流反推
将同一时刻的下游站点实测流量过程与下游站点出流估计流量过程相减，得到该目标河段估计区间入流过程。
2.根据权利要求1所述方法，其特征在于，步骤S3中，系数α和系数β的限定条件为：0<β<α<1，α+β＝1。
3.根据权利要求1所述方法，其特征在于，步骤S1中，所述流量资料为上游站点和下游站点在某个时间段内不同时刻的流量值，在所述流量资料中包括流量值和所述流量值对应的时间。
4.根据权利要求1所述方法，其特征在于，所述流量序列按照时间的先后顺序排列。
5.根据权利要求1所述方法，其特征在于，上游站点流量序列作为所述方法的输入变量，下游站点流量序列作为步骤S3中所述目标函数f的输入变量之一。
6.根据权利要求1所述方法，其特征在于，利用遗传算法或粒子群算法求目标函数f的最小值fmin。
CN201710134508.4A 2017-03-08 2017-03-08 一种基于马斯京根法反推区间入流的方法 Active CN106909797B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710134508.4A CN106909797B (zh) 2017-03-08 2017-03-08 一种基于马斯京根法反推区间入流的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710134508.4A CN106909797B (zh) 2017-03-08 2017-03-08 一种基于马斯京根法反推区间入流的方法

Publications (2)

Publication Number Publication Date
CN106909797A CN106909797A (zh) 2017-06-30
CN106909797B true CN106909797B (zh) 2019-05-24

Family

ID=59186382

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710134508.4A Active CN106909797B (zh) 2017-03-08 2017-03-08 一种基于马斯京根法反推区间入流的方法

Country Status (1)

Country Link
CN (1) CN106909797B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107491643A (zh) * 2017-08-18 2017-12-19 中国水利水电科学研究院 一种确定设计洪峰流量的方法、装置及系统
CN107688702B (zh) * 2017-08-24 2020-11-17 西安理工大学 一种基于狼群算法的河道洪水流量演进规律模拟方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102663267A (zh) * 2012-05-15 2012-09-12 南京大学 一种半湿润区流域面源污染负荷的确定方法
CN104462863A (zh) * 2015-01-07 2015-03-25 武汉大学 一种推求河道区间入流的计算方法
CN106056267A (zh) * 2016-05-12 2016-10-26 中国水利水电科学研究院 一种并联水库群优化调度方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105631168B (zh) * 2016-03-25 2018-12-28 中国水利水电科学研究院 一种实时高效的流域洪水演进可视化仿真方法
CN106202790B (zh) * 2016-07-20 2018-05-01 中国水利水电科学研究院 一种分布式河北模型构建方法及其应用

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102663267A (zh) * 2012-05-15 2012-09-12 南京大学 一种半湿润区流域面源污染负荷的确定方法
CN104462863A (zh) * 2015-01-07 2015-03-25 武汉大学 一种推求河道区间入流的计算方法
CN106056267A (zh) * 2016-05-12 2016-10-26 中国水利水电科学研究院 一种并联水库群优化调度方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
Muskingum river routing with dynamic bank storage;Birkhead A L,James C S;《Journal of Hydrology》;20021231;第264卷(第1期);第113-132页
on the subject of flood propagation computation method(muskingum method）;Cunge J A;《Journal of hydraulic Research》;19690731;全文
基于改进二分法的马斯京根法入流过程反演分析;张悦，李国芳，吴洁，郭财秀;《水电能源科学》;20150331;第33卷(第3期);第38-41页
基于马斯京根法推求入库洪水计算方法的改进;唐文涛，等;《水电能源科技》;20121231;第30卷(第12期);第51-54页
马斯京根模型改进新思路;王家彪，雷晓辉，廖卫红，王浩;《南水北调与水利科技》;20160430;第14卷(第2期);第87-92页

Also Published As

Publication number Publication date
CN106909797A (zh) 2017-06-30

Similar Documents

Publication Publication Date Title
CN111401666B (zh) 一种利用预报误差开展上游水库群影响下径流预报的方法
Zhang et al. Effectiveness of rainwater harvesting in runoff volume reduction in a planned industrial park, China
CN108108838A (zh) 一种高水量利用率的季调节水库优化调度方法
Yuting et al. Impacts of temperature and precipitation on runoff in the Tarim River during the past 50 years
CN108018823B (zh) 基于瞬时单位线的流域地下水平均响应时间估计方法
CN110852526B (zh) 一种基于雨洪过程相似性判别的实时洪水预报方法
CN108256780A (zh) 一种基于迭代库区回水计算的梯级水库防洪调度计算方法
CN106909797B (zh) 一种基于马斯京根法反推区间入流的方法
CN107392509A (zh) 一种基于物元分析的河道生态径流过程评价方法
CN109785979B (zh) 一种场次洪水降雨径流过程界定方法
CN104462863A (zh) 一种推求河道区间入流的计算方法
CN105354416B (zh) 一种基于代表电站的流域降雨径流电量宏观预测方法
Chaube et al. Synthesis of flow series of tributaries in Upper Betwa basin
Li et al. Quantifying the effects of channel change on the discharge diversion of Jingjiang Three Outlets after the operation of the Three Gorges Dam
CN104615883A (zh) 基于降雨过程理论的城市内河洪峰流量的计算方法
CN106383997A (zh) 一种反推三峡水库区间入流过程的计算方法
CN104594282B (zh) 一种水库运行调度图绘制方法
Wang et al. Combined use of volume-area and volume-length scaling relationships in glacio-hydrological simulation
CN106383935B (zh) 一种水库坝址天然日均水沙过程还原方法
Nguyen et al. Rainfall-storage-utilization-discharge model for flood mitigation and water conservation
Rahimi et al. Rooftop rainwater harvesting and its efficiency in Kabul new city
CN111047213B (zh) 基于多要素联合判别的中长期水资源量临界预警指标方法
Boroto Limpopo River: steps towards sustainable and integrated water resources management
Marinović et al. Karst groundwater quantity assessment and sustainability: the approach appropriate for river basin management plans
Ansorge Comparing various methods for determining the water footprint of electricity generation at Orlík hydroelectric power station–case study

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