DOI： 10.11817/j.issn.1672-7207.2020.09.004
基于滑动分析法的黄河下游造床流量变化研究
来志强1, 2，赵连军1，江恩慧1，吴国英1，张向1，程家林3
(1. 黄河水利委员会黄河水利科学研究院 水利部黄河泥沙重点实验室，河南 郑州，450003；
2. 中国科学院 地理科学与资源研究所，北京，100101；
3. 武汉大学 水资源与水电工程科学国家重点实验室，湖北 武汉，430072)
摘要：采用滑动分析法修正经典马卡维也夫法造床能力表达式，避免在分析流量频率时人为划分流量级带来的误差，基于黄河下游花园口水文站1950—2015年实测资料，分析黄河下游造床流量变化规律及新水沙关系。研究结果表明：改进后的造床流量计算方法中流量级长度参数在200~600 m3·s-1范围内、流量级滑动距离参数在1~600 m3·s-1范围内选取，均可准确计算出花园口第一和第二造床流量；1986—1999年，黄河下游第一和第二造床流量显著减少，洪水频次、洪峰流量和多年平均流量明显下降，下游河道淤积萎缩严重；1999年自小浪底水库运用及调水调沙以来，第二造床流量增加，下游河道主槽过流输沙能力提高。目前黄河下游造床能力已由第一造床流量主导转变为第二造床流量主控。
关键词：黄河下游；造床流量；滑动分析方法；水沙关系；马卡维也夫方法
中图分类号：TV143 文献标志码：A 开放科学(资源服务)标识码(OSID)
文章编号：1672-7207（2020）09-2396-09
Study on dominant discharge change in the lower Yellow River based on moving analysis method
LAI Zhiqiang1, 2, ZHAO Lianjun1, JIANG Enhui1, WU Guoying1, ZHANG Xiang1, CHENG Jialin3
(1. Key Laboratory of Yellow River Sediment Research, Yellow River Institute of Hydraulic Research,Zhengzhou 450003, China;
2. Institute of Geographic Sciences and Natural Resources Research, Chinese Academy of Sciences,Beijing 10010, China;
3. State Key Laboratory of Water Resources and Hydropower Engineering Science, Wuhan University,Wuhan 430072, China)
Abstract: The moving analysis method was employed to revise the bed-forming capacity equation based on the typical Makayev method. The inaccuracy from the man-made discharge stages can be avoided while analyzing the discharge appearance frequency. According to the observed data from the Huayuankou hydrological station from 1950 to 2015, the dominant discharge evolution in the lower Yellow River was analyzed. The new water and sediment relationship was explored. The results show that when the length of the discharge stage is chosen between 200 and 600 m3·s-1 and the moving distance of the discharge stage is chosen between 1 and 600 m3·s-1, the first and second dominant discharge of the observed data from the Huayuankou hydrological station can be obtained accurately. The first and second dominant discharge, flood-peak frequency, flood-peak discharge and average annual discharge in the lower Yellow River from 1986 to 1999 decrease significantly. The lower Yellow River channels also shrunken notably due to the sediment deposition. Since the operation of Xiaolangdi reservoir and water and sediment regulation in 1999, the second dominant discharge increases, and the discharge and sediment transport capability of the main channel in the lower Yellow River improves. The current dominant bed-forming capacity in the lower Yellow River changes from the first dominant discharge to the second one.
Key words: the lower Yellow River; dominant discharge; moving analysis method; water and sediment relationship; Makayev method
泥沙淤积是黄河下游洪水危害的本质原因。黄河下游河床宽浅、散乱，主流摆动不定，河势变化剧烈。泥沙不断淤积及复杂多变的河势使得防洪负担日益加剧。黄河下游河道冲淤演变特性很大程度上取决于来水来沙条件，但是天然河道的水沙过程随时间发生变化，因此，在进行河床演变分析和河道整治时，往往需要确定与其多年流量过程综合造床作用相当的特征流量，即造床流量。目前确定造床流量的理论并不成熟，常用的方法有平滩流量法[1]、输沙率法[2]、输沙能力法等。其中，输沙能力法在工程实践中应用最多，最经典的是马卡维也夫提出的计算方法[3]，通常也称为马氏造床流量方法。该方法通常会得到2个造床流量，其中，第一造床流量相当于多年平均最大洪水流量，第二造床流量稍大于多年平均流量。张红武等[4]采用马氏造床流量方法分析了黄河“92·8”高含沙洪水异常现象成因，并针对多沙河流造床特点，引入水流挟沙能力修正了该计算方法[5]。随后，韩其为等[6-9]定义塑造河床纵剖面的流量为第一造床流量，塑造横断面的流量为第二造床流量，并将其应用于黄河下游调水调沙对造床流量影响研究中。江恩惠等[10]针对黄河下游来水来沙人为调控特征与变化趋势，提出了洪、中、枯兼容的“三级流路”整治模式，即整治导线按中水(相当于马氏第一造床流量)流路设计。基于黄河下游各水文站1950—2003年实测资料，陈建国等[11]系统分析了50多年来黄河下游平滩流量和造床流量变化过程，提出增大河道的造床流量是治理河道萎缩的重要措施之一。陈绪坚等[12]通过计算黄河下游第一和第二造床流量，分析了黄河下游河槽的萎缩过程。胡春宏等[13]研究了1950—2003年黄河下游平滩流量、造床流量与洪水特征值的响应关系，阐明了黄河水沙过程调控与下游河道中水河槽塑造内在联系。孙东坡等[14]提出了针对黄河下游水沙变异特点的造床流量确定方法，指出维持黄河下游河道功能的河床塑造必须考虑适当的中常洪水，而不是长时间小流量。还有一些学者利用水文站长期观测数据探讨了造床流量与河床启冲流速[15]、平滩河宽[16]和洪峰流量重现期[17]之间的关系，以及通过物理试验研究了砾石河床[18]、砂质河床[19]演变规律与不同流量级冲刷之间的内在联系。造床流量对于分析天然河道水沙过程及河床演变特性十分重要。目前，大部分造床流量确定方法均以马氏造床流量方法为基础，流量级的人为界定往往带来误差。近10年来，随着黄河水沙调控工程体系的运用，黄河下游水沙关系发生了很大变化。因此，本文作者基于黄河下游花园口水文站1950—2015年实测数据，通过滑动分析法对马氏造床能力表达式进行修正，避免人为划分流量级导致的计算误差，研究近年来黄河下游造床流量变化规律，旨在探明人为调控作用下黄河下游水沙过程与河床的适应性。
1 黄河下游造床流量确定方法
1.1 马氏造床流量
经典马氏造床流量认为某个流量的造床作用既与该流量的输沙能力有关，也与该流量所持续时间有关。与流量Q输沙能力有关的造床作用与Q的m次方及比降J的乘积成正比，其中m为双对数坐标上实测断面日均输沙率Gs与相应日均流量Q关系曲线的斜率，与流量Q所持续时间有关的造床作用可用该流量出现的频率P来表示[1,3]，P的求解方式为将河段断面随时间变化的流量过程分为若干流量级，统计各流量级在各流量级内出现的频率。马氏造床能力的表达式为QmJP。诸多学者[20-27]在运用马氏造床流量计算方法时发现，人为划分流量级差会给造床流量的确定带来很大影响。
1.2 滑动分析方法
滑动分析方法通过顺序逐期增减新旧数据求取移动平均值或累加值，借以消除偶然变动和人为主观因素的影响，近年来该方法已被广泛地应用于黄河来沙来水量的计算分析[28-29]。本文基于滑动分析方法原理，对马氏造床流量进行修正，以消除人为划分流量级(参数P)的不良影响，具体计算步骤如下。
1) 定义一系列流量级为
(1)
式中：i为流量级的序号；n为流量级的总数；
和
分别为第i级流量的上、下界；
为每组流量级长度；
为流量级滑动距离；
，
，具体定义如图1所示，其中
和
分别为第i+1级流量的上、下界。
图1 流量级、
和
定义示意图
Fig. 1 Sketch map of definitions of discharge stages, and
2) 定义第1流量级
的初始值为所分析流量数据中最小值
，则该系列流量级
可表述为
；
(2)
对于流量级
，定义该流量级的表征流量
为
；统计所分析流量落入该流量级的流量
，对其造床能力进行累加以求取流量级的表征造床能力
，式中，
为落入该流量级内
的数目；当流量级上界滑动至所分析流量中的最大值
时，计算结束。
3) 绘制
图，积分求取其围成的面积
；对于某个水文站某段时间内水沙过程，S为常数；为了便于对比不同时间段流量的造床能力，对进行归一化处理，即Φi=φi/S；-Φi变化规律与
变化规律是一致的，两者横坐标相同，区别仅在于纵坐标的值不同(Φi纵坐标值是
纵坐标值的1/S倍)；绘制
Φi图，在图中查出Φi的最大值和次大峰值，分别为第一造床能力和第二造床能力。2个流量中，流量较大值为第一造床流量，较小值为第二造床流量。
2 黄河下游水沙关系研究
2.1 参数敏感分析
基于黄河下游花园口水文站1950—2015年实测数据(地理位置如图2所示)，采用滑动分析法修正马氏造床流量计算方法，分析黄河下游造床流量变化规律。首先在双对数坐标系中点绘花园口1950—2015年日均流量Q和相应日均输沙率Gs，如图3所示，对数据进行拟合，得到Q和Gs关系为Gs=5e-5Q2.1(e为自然底数；R2为函数拟合可决系数，R2=0.9，其越接近1，表明函数拟合程度越好)。
图2 黄河部分水文站和水库地理位置
Fig. 2 Geographic locations of some hydrological stations and reservoirs in the Yellow River
图3 花园口1950—2015年流量Q与输沙率Gs关系
Fig. 3 Relationship between discharge Q and sediment transport rate Gs based on data from the Huayuankou hydrological station from 1950 to 2015
马氏造床流量方法中流量级长度取值直接影响其组数，根据文献[20-27]，马氏方法选取流量级组数从十几组到近百组不等。根据实测资料，花园口水文站1986—2015年实测日均流量中最小值为7.58 m3·s-1，最大值为7 270 m3·s-1。当流量级取100 m3·s-1时，流量级组数为73个；当其取800 m3·s-1时，流量级组数为10个。因此，选取100和800 m3·s-1以分别代表流量级组数较多和较少时情况，进而对比分析2种流量级划分对造床流量的影响，计算结果如图4所示，其中比降J取花园口实测比降资料的平均值0.000 2。由图4(a)可知：当流量级为100 m3·s-1时，计算得到第一造床流量(粉色标记)和第二造床流量(红色标记)分别为1 260和860 m3·s-1。由图4(b)可知：当流量级为800 m3·s-1时，第一造床流量和第二造床流量分别为5 050和1 060 m3·s-1。可以看出，人为划分流量级对马氏造床流量计算结果影响较大；若流量级取得小，则会造成某一小流量级因出现频率过高而掩盖真实值；若流量级取得大，则计算精度不高。
图4 花园口1986—2015年不同流量级下马氏造床流量分析
Fig. 4 Analysis of dominant discharge under different man-made discharge stages using Makayev method based on data from the Huayuankou hydrological station from 1986 to 2015
分析滑动分析法中和这2个参数对计算结果的影响。以花园口水文站1986—2015年实测数据为研究对象，首先选定=1 m3·s-1(最小的正整数流量)，这样，流量级滑动时既可确保流量数据得到充分利用，又不致于计算量太大，然后分别计算当为200，300，400，500和600 m3·s-1时Qi-Φi关系曲线，如图5所示。从图5可以看出：随着增大，Φi整体分布略降低，但规律上几乎没有变化，第一和第二造床流量分布没有变化。
图5 花园口1986—2015年不同ΔQ1下Qi-Φi关系曲线
Fig. 5 Relationship curves between Qi and Φi with different ΔQ1 based on data from the Huayuankou hydrological station from 1986 to 2015
另一方面，选定=500 m3·s-1(由图5可知在200~600 m3·s-1范围内对造床流量的确定没有影响)，设定为1，50，100，300和600 m3·s-1，以分析对Qi-Φi关系曲线的影响，如图6所示。由图6可知：当增大时，由于数据点减少，Qi-Φi关系曲线变得越来越不光滑(特别是当大于时)，但整体分布规律几乎没有变化，第一和第二造床流量分布受其影响不大。因此，建议在使用该方法时小于，这样可确保Qi-Φi关系曲线中有效数据的数目，以更加真实地反映曲线的分布规律。
图6 花园口1986—2015年不同ΔQ2下Qi-Φi关系曲线
Fig. 6 Relationship curve between Qi and Φi with different ΔQ2 based on data from the Huayuankou hydrological station from 1986 to 2015
2.2 不同代表年段花园口造床流量计算
根据图5和图6计算结果可知：在200~600 m3·s-1范围内、在1~600 m3·s-1范围内对造床流量的确定没有影响。考虑黄河造床流量计算时经常采用的马氏流量级组数以及黄河日均流量分布，计算时选定为500 m3·s-1，选定为 1 m3·s-1。近年黄河水沙变化较典型年份有1986年和1999年：1986年，黄河上游龙羊峡水库和刘家峡水库得到联合运用；1999年，黄河下游小浪底水库得到运用。本文将研究时段分为1950—1985年、1986—1999年、2000—2015年以及1950—1999年、1986—2015年，相应的Qi-Φi关系曲线如图7所示。在图7中分别找到Φi的最大峰值和次大峰值对应的特征流量，其中，较大流量为第一造床流量，较小流量为第二造床流量，如表1所示。
由图7和表1可知：1950—1985年花园口第一造床流量的造床能力大于第二造床流量的造床能力；1986—1999年水沙条件发生变化，花园口第一和第二造床流量均显著减小，且第一造床流量的造床能力小于第二造床流量的造床能力。可以看出，1986年以后黄河下游水沙关系发生了显著变化，第一造床流量明显降低，意味着洪水频次和洪峰流量显著减少；第二造床流量也显著降低，表明来水流量减少，这是造成下游河道淤积萎缩的主要原因。
另一方面，相比于1986—1999年，花园口2000—2015年第一造床流量降低，第二造床流量增加，且第一造床流量的造床能力小于第二造床流量的造床能力。通过分析可知，自1999年小浪底水库运用及调水调沙以来，下游河道发生持续冲刷，主槽过流输沙能力提高，第二造床流量增加。为了更加全面分析目前黄河下游河道造床流量分布现状，本文还计算了黄河下游高村水文站(距上游的花园口水文站近189 km，地理位置如图2所示)和泺口水文站(距上游的花园口水文站近469 km，地理位置图2所示)1986—2015年第一和第二造床流量。计算时，高村水文站比降为实测比降的平均值即0.000 15，拟合Gs与Q得到m为2.0；泺口水文站比降为实测比降的平均值即0.000 1，拟合Gs与Q得到m为1.8，计算结果如图8和表1所示。可以看出，两者造床能力现状均为第二造床流量大于第一造床流量。
2.3 滑动分析法与马氏方法对比
采用滑动分析法(和分别为图5和图6所示范围)和马氏方法(流量级为100和800 m3·s-1)对花园口1950—1985年、1986—2015年造床流量进行计算，所得造床流量如表2所示。由表2可知：马氏方法中流量级对造床流量影响很大，流量级取值不同，计算得到的造床流量确定值差异很大；而当在200~600 m3·s-1范围内、在1~600 m3·s-1范围内时，采用滑动分析法均可准确计算出第一和第二造床流量。
表1 不同研究时段花园口第一、第二造床流量
Table 1 The first and second dominant discharge based on data from the Huayuankou hydrological station under different research periods
表2 基于滑动分析法和马氏方法计算得到的花园口造床流量
Table 2 Dominant discharge based on data from the Huayuankou hydrological station according to moving analysis method and the Makayev method
图7 花园口不同研究时段Qi-Φi关系曲线
Fig. 7 Relationship curves between Qi and Φi based on data from the Huayuankou hydrological station under different research periods
图8 高村和泺口水文站1986—2015年Qi-Φi关系曲线
Fig. 8 Relationship curves between Qi and Φi based on data from the Gaocun and Luokou hydrological stations from 1986 to 2015
3 结论
1) 滑动分析法的引入有效避免了马氏造床流量计算方法中人为划分流量级差带来的误差。当流量级长度参数在200~600 m3·s-1范围内、流量级滑动距离参数在1~600 m3·s-1范围内时，改进后的造床流量计算方法对第一和第二造床流量的确定没有影响，均可准确计算出造床流量。
2) 1986—1999年，黄河下游水沙条件发生了显著的变化，第一造床流量和第二造床流量显著减少，下游河道淤积萎缩严重；自1999年小浪底水库运用及调水调沙以来，第二造床流量增加，下游河道主槽过流输沙能力提高。
3) 目前黄河下游造床能力已由第一造床流量大于第二造床流量转变为第二造床流量大于第一造床流量。
参考文献：
[1] 谢鉴衡. 河流模拟[M]. 北京: 水利电力出版社, 1990: 1-36.
XIE Jianheng. River simulation[M]. Beijing: Water Resources and Electric Power Press, 1990: 1-36.
[2] 钱宁, 张仁, 周志德. 河床演变学[M]. 北京:科学出版社, 1987: 1-28.
QIAN Ning, ZHANG Ren, ZHOU Zhide. Riverbed evolution[M]. Beijing: Science Press, 1987: 1-28.
[3] 谢鉴衡. 河床演变及整治[M]. 2版. 北京: 中国水利水电出版社, 1997: 23-24.
XIE Jianheng. Riverbed evolution and regulation[M]. 2nd ed. Beijing: China Water & Power Press, 1997: 23-24.
[4] 张红武, 江恩惠, 吕昕. 黄河高含沙洪水异常现象成因分析[J]. 人民黄河, 1993(9): 8-9.
ZHANG Hongwu, JIANG Enhui, LU Xin. Analysis of the causes of the abnormal phenomena of the high sediment concentration flood in the Yellow River[J]. Yellow River, 1993(9): 8-9.
[5] 张红武,张清,江恩惠.黄河下游河道造床流量的计算方法[J].泥沙研究,1994(4): 50-55.
ZHANG Hongwu, ZHANG Qing, JIANG Enhui. Calculation of dominated discharge in the lower Yellow River[J]. Journal of Sediment Research, 1994(4): 50-55.
[6] 韩其为. 黄河下游输沙及冲淤的若干规律[J]. 泥沙研究, 2004(3): 1-13.
HAN Qiwei. Some rules of sediment transportation and deposition-scouring in the lower Yellow River[J]. Journal of Sediment Research, 2004(3): 1-13.
[7] 韩其为. 小浪底水库初期运用及黄河调水调沙研究[J]. 泥沙研究, 2008(3): 1-18.
HAN Qiwei. Study on preliminary operation of Xiaolangdi Reservoir and flow-sediment regulation of the Yellow River[J]. Journal of Sediment Research, 2008(3): 1-18.
[8] 韩其为. 第一造床流量及输沙能力的理论分析: “黄河调水调沙的根据、效益与巨大潜力”之三[J]. 人民黄河, 2009, 31(1): 1-4, 7.
HAN Qiwei. Theoretical analysis on first bed-forming discharge and sediment transporting capacity[J]. Yellow River, 2009, 31(1): 1-4, 7.
[9] 韩其为, 江恩惠, 陈绪坚. 黄河下游第二造床流量研究: “黄河下游调水调沙的根据、效益和巨大潜力”之四[J]. 人民黄河, 2009, 31(2): 1-4.
HAN Qiwei, JIANG Enhui, CHEN Xujian. Study on second bed-forming discharge of the lower Yellow River[J]. Yellow River, 2009, 31(2): 1-4.
[10] 江恩惠, 曹永涛, 张林忠, 等.黄河下游游荡型河段河势演变规律及机理研究[M]. 北京: 中国水利水电出版社, 2006: 1-30.
JIANG Enhui, CAO Yongtao, ZHANG Linzhong, et al. Study on evolution law and mechanism of river regime in wandering reach of the lower Yellow River[M]. Beijing: China Water & Power Press, 2006: 1-30.
[11] 陈建国, 胡春宏, 董占地, 等. 黄河下游河道平滩流量与造床流量的变化过程研究[J]. 泥沙研究, 2006(5): 10-16.
CHEN Jianguo, HU Chunhong, DONG Zhandi, et al. Change of bankfull and bed-forming discharges in the Lower Yellow River[J]. Journal of Sediment Research, 2006(5): 10-16.
[12] 陈绪坚, 韩其为, 方春明. 黄河下游造床流量的变化及其对河槽的影响[J]. 水利学报, 2007, 38(1): 15-22.
CHEN Xujian, HAN Qiwei, FANG Chunming. Variation of dominant discharge in Lower Yellow River and its influence on river channel[J]. Journal of Hydraulic Engineering, 2007, 38(1): 15-22.
[13] 胡春宏, 陈建国, 郭庆超. 黄河水沙过程调控与下游河道中水河槽塑造[J]. 天津大学学报(自然科学与工程技术版), 2008, 41(9): 1035-1040.
HU Chunhong, CHEN Jianguo, GUO Qingchao. Regulation of flow and sediment processes and reform of medium-size channel in the Lower Yellow River[J]. Journal of Tianjin University(Science and Technology), 2008, 41(9): 1035-1040.
[14] 孙东坡, 王勤香, 王鹏涛, 等. 基于水沙关系系数法确定黄河下游造床流量[J]. 水力发电学报, 2013, 32(1): 150-155.
SUN Dongpo, WANG Qinxiang, WANG Pengtao, et al. Determination of dominant discharge of the lower Yellow River based on coefficient frequencies of flow-sediment relationship[J]. Journal of Hydroelectric Engineering, 2013, 32(1): 150-155.
[15] 马良, 张红武, 马睿, 等. 多沙河流治导线流路方程研究[J]. 水利学报, 2017, 48(3): 279-290.
MA Liang, ZHANG Hongwu, MA Rui, et al. Study on the regulation line flow path equation of the sandy river[J]. Journal of Hydraulic Engineering, 2017, 48(3): 279-290.
[16] COPELAND R, SOAR P, THORNE C. Channel-forming discharge and hydraulic geometry width predictors in meandering sand-bed Rivers[C]// Americam Society of Civil Engineers: Impacts of Global Climate Change. Anchorage, Alaska, United States: World Water & Environmental Resources Congress, 2005: 1-12.
[17] FERRO V, PORTO P. Identifying a dominant discharge for natural Rivers in southern Italy[J]. Geomorphology, 2012, 139/140: 313-321.
[18] KEYLOCK C J, SINGH A, FOUFOULA-GEORGIOU E. The complexity of gravel bed river topography examined with gradual wavelet reconstruction[J]. Journal of Geophysical Research: Earth Surface, 2014, 119(3): 682-700.
[19] SARKAR K, MAZUMDER B S. Space time evolution of sand bed topography and associated flow turbulence: experiments with statistical analysis[J]. Stochastic Environmental Research and Risk Assessment, 2018, 32(2): 501-525.
[20] 董冠营, 赵文宾. 淮河盱眙段河道造床流量的计算[J].中国水运, 2013, 13(8): 164-165.
DONG Guanying, ZHAO Wenbin. Calculation of bed forming discharge in Xuyi section of Huaihe River[J]. China Water Transport, 2013, 13(8): 164-165.
[21] 楚万强, 李永丽, 谢龙. 长江干流重庆主城河段造床流量计算[J]. 重庆交通大学学报(自然科学版), 2015, 34(2): 69-71, 117.
CHU Wanqiang, LI Yongli, XIE Long. Calculation of dominant discharge of Chongqing urban reach in Yangtze River main stream[J]. Journal of Chongqing Jiaotong University(Natural Science), 2015, 34(2): 69-71, 117.
[22] 张玮, 徐晶鑫. 长江下游感潮河段造床流量计算研究[J]. 水运工程, 2017(4): 8-14.
ZHANG Wei, XU Jingxin. Study on dominant discharges in the lower tidal reach of the Yangtze River[J]. Port & Waterway Engineering, 2017(4): 8-14.
[23] 张为, 高宇, 许全喜, 等. 三峡水库运用后长江中下游造床流量变化及其影响因素[J]. 水科学进展, 2018, 29(3): 331-338.
ZHANG Wei, GAO Yu, XU Quanxi, et al. Changes in dominant discharge and their influential factors in the middle and lower reaches of Yangtze River after the Three Gorges Dam impoundment[J]. Advances in Water Science, 2018, 29(3): 331-338.
[24] 宝婉宁, 郭维东, 李翔, 等. 辽河干流中下游河段造床流量和平滩流量及有效流量分析[J]. 泥沙研究, 2018, 43(2): 55-60.
BAO Wanning, GUO Weidong, LI Xiang, et al. Study on dominant discharge, bankfull discharge and effective discharge of Liaohe River[J]. Journal of Sediment Research, 2018, 43(2): 55-60.
[25] 贾艳红,伍悦滨,王励,等. 浅析马卡维耶夫法在松花江中下游河道造床流量分析中的应用[J]. 东北水利水电, 2006, 24(2): 35-36, 40.
JIA Yanhong, WU Ruibin, WANG Li, et al. Application of Makayev in analysis of bed-building discharge in middle and downstream of Songhuajiang river[J]. Water Resources & Hydropower of Northeast China, 2006, 24(2): 35-36, 40.
[26] 肖立敏, 朱立俊, 范红霞. 造床流量法中变指数m值的确定[J]. 东北水利水电, 2007, 25(7): 1-3.
XIAO Limin, ZHU Lijun, FAN Hongxia. Determination of parameter m in method of bed-building discharge[J]. Water Resources & Hydropower of Northeast China, 2007, 25(7): 1-3.
[27] 虞邦义, 郁玉锁, 赵凯. 淮河中游造床流量计算[J]. 河海大学学报(自然科学版), 2010, 38(2): 210-214.
YU Bangyi, YU Yusuo, ZHAO Kai. Calculation of dominant discharge in the middle reaches of Huaihe River[J]. Journal of Hohai University(Natural Sciences), 2010, 38(2): 210-214.
[28] 钱胜, 马翠丽, 崔振华. 维持黄河下游中水河槽的小浪底调控指标分析[J]. 人民黄河, 2017, 39(11): 15-18.
QIAN Sheng, MA Cuili, CUI Zhenhua. Analysis of regulation indicators of the Xiaolangdi reservoir based on maintaining the lower yellow river for normal channel[J]. Yellow River, 2017, 39(11): 15-18.
[29] 陈琳, 胡春宏, 陈绪坚.黄河下游河道平滩流量与水沙过程响应关系研究[J]. 泥沙研究, 2018, 43(4): 1-7.
CHEN Lin, HU Chunhong, CHEN Xujian. Relationship between bank-full discharges and processes of flow-sediment in lower Yellow River[J]. Journal of Sediment Research, 2018, 43(4): 1-7.
(编辑 杨幼平)
收稿日期： 2019 -11 -19; 修回日期： 2020 -04 -23
基金项目(Foundation item)：国家自然科学基金资助项目(51539004，51909102)；中央级公益性科研院所基本科研业务费专项项目(HKY-JBYW-2019-02)；第四届中国科协青年人才托举工程全额资助项目(水学合同201938)；黄河水利科学研究院科技发展基金资助项目(黄科发201907) (Projects(51539004, 51909102) supported by the National Natural Science Foundation of China; Project(HKY-JBYW-2019-02) supported by the Fundamental Research Fund of Central Public-Interest Scientific Institution; Project(201938) supported by the Young Elite Scientists Sponsorship Program by China Association for Science and Technology (CAST); Project(YRIHR201907) supported by the Science and Technology Development Fund of YRIHR)
通信作者：赵连军，博士，教授级高级工程师，从事河流泥沙动力学研究；E-mail：zhaolianjun88@163.com
