蔡蓉蓉 張紅武 卜海磊

摘 要：為避免歷史調(diào)查洪水誤差較大對設(shè)計洪水計算的影響和我國過于保守的分析思路導(dǎo)致設(shè)計洪水值普遍偏大而影響黃河下游防洪體系合理運(yùn)行的現(xiàn)狀，利用黃河下游花園口站百年實(shí)測洪峰流量還原資料，通過洪水頻率分析和調(diào)洪計算對花園口站在小浪底等水庫運(yùn)用前后的設(shè)計洪水頻率及洪水過程進(jìn)行研究，并結(jié)合洪峰流量相關(guān)圖確定了黃河下游各站的設(shè)計洪峰流量。結(jié)果表明：花園口站天然狀態(tài)和考慮小浪底等水庫運(yùn)用的100 a一遇設(shè)計洪峰流量分別為22 949、12 948 m3/s，黃河下游夾河灘、高村、孫口站在小浪底等水庫運(yùn)用后100 a一遇設(shè)計洪峰流量分別為11 730、11 050、10 100 m3/s;花園口站天然狀態(tài)和考慮小浪底等水庫運(yùn)用的1 000 a一遇設(shè)計洪峰流量分別為32 301、20 865 m3/s，小浪底等水庫運(yùn)用后相應(yīng)的夾河灘、高村、孫口站1 000 a一遇設(shè)計洪峰流量分別為19 125、17 100、15 350 m3/s。對比其他學(xué)者的成果和目前采用的設(shè)計洪水標(biāo)準(zhǔn)，所得結(jié)果較小，現(xiàn)有堤防設(shè)防標(biāo)準(zhǔn)自然相應(yīng)提高，表明黃河下游防洪安全客觀上存在一定余地，為解放北金堤滯洪區(qū)與明顯減小東平湖滯洪區(qū)的分洪概率和分洪量提供了良好條件，有利于采用“兩道防線”與生態(tài)治理、因?yàn)┦┲蜗嘟Y(jié)合的“三灘分治”模式治理黃河，進(jìn)而使黃河下游成為造福兩岸人民的幸福河。

關(guān)鍵詞：洪水頻率分析;調(diào)洪計算;設(shè)計洪水;花園口站;黃河下游

Abstract：In order to avoid the influence of the potential errors of the estimated historical floods to the design flood estimation and the potential influence of the overly conservative results of the design flood estimation to the current operation of the flood control system in the lower Yellow River （LYR）， this study conducted the flood frequency analysis of 100-year restored flood peak flow series of the Huayuankou Station （HYK） to quantify the flood frequency of the HYK before and after the operation of the Xiaolangdi Reservoir （XLDR） and estimate the corresponding design flood discharges of several stations （Jiahetan （JHT）， Gaocun （GC） and Sunkou （SK）） over the LYR by using the flood flow correlation diagrams. The results show that the calculated natural （before the XLDR） and current （after the XLDR） 100-year flood peak discharges over the HYK are 22 949 and 12 948 m3/s， while the corresponding values for the JHT， GC， and SK after the XLDR are 11 730， 11 050 and 10 100 m3/s respectively. Besides， the natural and current 1 000-year flood peak discharges of the HYK are 32 301 and 20 865 m3/s， while those values of the JHT， GC and SK after the XLDR are 19 125， 17 100 and 15 350 m3/s respectively. In comparison with the estimations of many other studies and the current design flood standards， the results shown in this study are smaller， indicating that the current levee protection standard over the LYR can provide sufficient security margins. The potential sufficient margins can provide proper conditions for the relief of the flood detention zones over the Beijindi and the reduction of the usage of the flood detention zones over the Dongping Lake and is conducive to the adoption of the “Two Lines of Defense” river management strategies combined with the ecological management and the basin specific management strategy “Three basins， Three Strategies”， thus making the LYR a “happiness river” for the benefit of the people.

Key words： flood frequency analysis; flood regulation calculation; design flood; Huayuankou Station; Lower Yellow River

1 引 言

設(shè)計洪水量級和河流治理對策的確定與具體防洪工程的布局密切相關(guān)[1-2]。為解決水文實(shí)測資料短缺問題，我國自20世紀(jì)50年代初期以來，借鑒日偽時期修建豐滿、水豐、二龍山等大型水庫工程確定設(shè)計洪水時，將發(fā)掘的遠(yuǎn)年特大歷史洪水資料應(yīng)用到設(shè)計洪水頻率分析計算中的做法[3]，采用歷史調(diào)查洪水和實(shí)測洪水資料相結(jié)合的方式，通過洪水頻率分析方法確定設(shè)計洪水[4-7]。具體到黃河流域，筆者通過查閱《黃河流域歷史洪水調(diào)查考證及計算工作的體會》[8]，發(fā)現(xiàn)了進(jìn)行黃河流域洪水頻率分析時計入歷史調(diào)查洪水的理由。該報告指出：黃河流域干支流水文資料長度較短，新中國成立以前資料殘缺不全，采用數(shù)理統(tǒng)計方法進(jìn)行頻率計算時，通過長度較短的水文資料進(jìn)行外延以分析百年、千年、萬年一遇設(shè)計洪水會產(chǎn)生較大的抽樣誤差;進(jìn)行洪水頻率分析時，若只采用實(shí)測（觀測）資料，統(tǒng)計參數(shù)和設(shè)計結(jié)果會隨資料長度的增加產(chǎn)生顯著變化，而在計算中加入歷史調(diào)查洪水后，即使資料長度不斷增加，成果也不會產(chǎn)生大的變化。由此看出，該理由前部分是合理的，“通過長度較短的水文資料進(jìn)行外延以分析百年、千年、萬年一遇設(shè)計洪水”，由于所選擇資料系列相對于洪水系列總體不具代表性，因此會產(chǎn)生較大的誤差而缺乏可靠性。不過，該理由后部分則暴露出一些問題，“采用實(shí)測（觀測）資料，統(tǒng)計參數(shù)和設(shè)計結(jié)果會隨資料長度的增加產(chǎn)生顯著變化”實(shí)際是符合自然規(guī)律的（只是所用“顯著”一詞不那么確切，除非實(shí)際出現(xiàn)非常洪水，才可能產(chǎn)生“顯著變化”），不能作為加入歷史調(diào)查洪水資料的理由，至于“在計算中加入歷史調(diào)查洪水后，即使資料長度不斷增加，成果也不會產(chǎn)生大的變化”，卻暴露出加入歷史調(diào)查洪水資料的最大缺陷[9]，無論如何，都不能排除實(shí)測（觀測）資料對設(shè)計洪水成果的影響。

具體而言，進(jìn)行黃河流域歷史洪水洪峰流量計算時（即確定歷史調(diào)查洪水時）多采用比降法[3，10-13]?？紤]到水深與糙率n之間的關(guān)系通常并不十分明確[14]，比降J不易測量[15]，采用綜合系數(shù)Kf（Kf= J/n=V/R2/3）代替比降計算洪峰流量[16]。以1843年洪水事件時三門峽站洪峰流量計算為例，采用三門峽上口站實(shí)測資料繪制水力半徑R與綜合系數(shù)Kf的關(guān)系曲線，根據(jù)洪水時的水力半徑求出對應(yīng)的綜合系數(shù)Kf，再利用綜合系數(shù)、斷面面積、水力半徑計算出洪峰流量[17]，原始數(shù)據(jù)見表1。雖然采用綜合系數(shù)法減小了直接使用糙率n和比降J計算帶來的誤差，但仍存在一定的問題，如計算使用的斷面面積要素等實(shí)際難以確切量化，且計算假定特大洪水時的R—Kf關(guān)系與中常洪水時的R—Kf關(guān)系相同，也不符合實(shí)際。鑒于特大洪水發(fā)生年代久遠(yuǎn)，調(diào)查結(jié)果與真實(shí)情況之間難免存在偏差，再考慮到高含沙水流的特殊性質(zhì)，更不能直接將適用于中常洪水的水文指標(biāo)關(guān)系式應(yīng)用于特殊洪水計算，歷史洪水調(diào)查成果確實(shí)需要進(jìn)一步修正[9]。

綜上，過去進(jìn)行洪水頻率分析計算時資料年限較短，不具有代表性，為了穩(wěn)定設(shè)計洪水結(jié)果，將歷史調(diào)查洪水計入資料系列內(nèi)?，F(xiàn)如今，已積累了大量實(shí)測水文資料，同時考慮到當(dāng)年歷史調(diào)查洪水成果的誤差較大，且經(jīng)計算發(fā)現(xiàn)進(jìn)行洪水頻率計算時即使不計入歷史調(diào)查洪水，只要洪水資料連續(xù)系列較長且包含典型豐中枯水的資料，近30多a沒有出現(xiàn)大洪水，對設(shè)計洪水結(jié)果影響有限，洪水頻率也相對穩(wěn)定，且分析成果更為合理[9，18]。已有學(xué)者指出，我國進(jìn)行設(shè)計洪水計算時在各個環(huán)節(jié)偏于保守，使得設(shè)計洪水值普遍偏大，造成不必要的投資浪費(fèi)[6]，甚至影響大江大河防洪體系的合理運(yùn)行[9]。在當(dāng)今的水沙情勢下，應(yīng)科學(xué)修正設(shè)計洪水成果，使得防洪減災(zāi)體系合理科學(xué)，與新時代的發(fā)展情勢相適應(yīng)[9]。在如今可靠的實(shí)測洪水資料已較為豐富的背景下，根據(jù)國家重點(diǎn)科技研發(fā)計劃專項(xiàng)“黃河下游河道與灘區(qū)治理研究”項(xiàng)目第一課題“黃河下游未來洪水條件及其災(zāi)害情景”任務(wù)書要求，為克服傳統(tǒng)的先利用稀遇洪水的歷史調(diào)查洪水點(diǎn)據(jù)與部分實(shí)測資料將頻率曲線走向確定后，再調(diào)整CS/CV進(jìn)行經(jīng)驗(yàn)適線所導(dǎo)致的“即使資料長度不斷增加，成果也不會產(chǎn)生大的變化”的缺陷，盡可能利用實(shí)測資料而不采用歷史調(diào)查洪水點(diǎn)據(jù)，對黃河下游洪水頻率展開分析計算，以給出適用于當(dāng)前與未來水沙情勢的設(shè)計洪水成果。

2 研究數(shù)據(jù)與方法

2.1 研究數(shù)據(jù)

收集花園口站1919—1949年歷年最大洪峰流量資料、1950—2010年實(shí)測洪峰流量還原資料[7]、三門峽站與花園口站2011—2020年年最大流量資料，其中花園口站歷年最大洪峰流量如圖1所示。2000年后由于缺少大洪水流量過程，洪水資料還原與否對設(shè)計洪水結(jié)果影響基本在分析的誤差范圍之內(nèi)，因此2011—2020年資料基本可以直接采用實(shí)測資料。為偏于安全，對這10 a資料，如果某年份花園口站年最大流量小于三門峽站的，則該年采用三門峽站最大流量資料，相當(dāng)于小浪底水庫沒有起到攔蓄作用。

2.3 水庫運(yùn)用方式

根據(jù)胡一三等[22]的研究成果查得小浪底水庫2002年汛前的庫容及泄流量（見表2），小浪底水庫防洪限制水位為275 m，允許防洪運(yùn)用庫容為40.5億m3。根據(jù)水位庫容曲線與40.5億m3防洪庫容，確定起始計算庫容為71.5億m3，通過線性插值求得對應(yīng)水位為258 m，根據(jù)水位通過線性插值得到對應(yīng)的泄流能力。擬定小浪底水庫100 a一遇洪水時的運(yùn)行方式：①當(dāng)入庫流量小于10 000 m3/s時，出庫流量等于入庫流量;②當(dāng)入庫流量大于等于10 000 m3/s時，按湊泄花園口站流量為10 000 m3/s（當(dāng)入庫流量小于12 000 m3/s時，陸渾、故縣、河口村水庫按合計1 000 m3/s下泄，小花區(qū)間來水按2 000 m3/s計算;當(dāng)入庫流量大于等于12 000 m3/s時，陸渾、故縣、河口村水庫不下泄，小花區(qū)間來水按2 000 m3/s計算）;③當(dāng)入庫流量大于等于10 000 m3/s，但湊泄花園口流量為10 000 m3/s無法滿足時，出庫流量按對應(yīng)水位的泄流能力計算。

小浪底水庫300 a一遇洪水時的運(yùn)行方式：①當(dāng)入庫流量小于10 000 m3/s時，出庫流量等于入庫流量;②當(dāng)入庫流量大于等于10 000 m3/s時，按湊泄花園口流量為10 000 m3/s（當(dāng)入庫流量小于12 000 m3/s時，陸渾、故縣、河口村水庫按合計1 250 m3/s下泄，小花區(qū)間來水按2 500 m3/s計算;當(dāng)入庫流量大于等于12 000 m3/s時，陸渾、故縣、河口村水庫不下泄，小花區(qū)間來水按2 500 m3/s計算）;③當(dāng)入庫流量大于等于10 000 m3/s，但湊泄花園口流量為10 000 m3/s無法滿足時，出庫流量按對應(yīng)水位的泄流能力計算。

小浪底水庫1 000 a一遇洪水時的運(yùn)行方式：①當(dāng)入庫流量小于10 000 m3/s時，出庫流量等于入庫流量;②當(dāng)入庫流量大于等于10 000 m3/s時，按湊泄花園口流量為12 000 m3/s（當(dāng)入庫流量小于12 000 m3/s時，陸渾、故縣、河口村水庫按合計2 000 m3/s下泄，小花區(qū)間來水按4 000 m3/s計算;當(dāng)入庫流量大于等于12 000 m3/s時，陸渾、故縣、河口村水庫不下泄，小花區(qū)間來水按4 000 m3/s計算）;③當(dāng)入庫流量大于等于10 000 m3/s，但湊泄花園口流量為12 000 m3/s無法滿足時，出庫流量按對應(yīng)水位的泄流能力計算。

3 花園口站設(shè)計洪水分析

3.1 天然狀態(tài)下花園口站設(shè)計洪水成果

采用2.1節(jié)中的研究數(shù)據(jù)，即nf=100，根據(jù)2.2節(jié)中的分析方法，得到花園口站洪水頻率分析結(jié)果，見圖2和表3。由此得到考慮小浪底等水庫影響前，天然狀態(tài)下花園口站的設(shè)計洪水成果，100 a一遇設(shè)計洪峰流量為22 949 m3/s，300 a一遇設(shè)計洪峰流量為27 424 m3/s，1 000 a一遇設(shè)計洪峰流量為32 301 m3/s。將張紅武等修正的較為符合實(shí)際的1843年洪水花園口洪峰流量為33 000 m3/s、重現(xiàn)期為1 200 a[9]，點(diǎn)繪在圖2中，可看出其點(diǎn)據(jù)恰好在設(shè)計頻率曲線上，表明上述成果較為合理。

3.2 考慮水庫影響時花園口站設(shè)計洪水成果

李保國等[7]通過洪水頻率分析得到天然狀態(tài)下的設(shè)計洪水成果后，還根據(jù)不同量級設(shè)計洪水成果與典型洪水過程線放大得到對應(yīng)量級的設(shè)計洪水過程線，再通過調(diào)洪計算得到考慮水庫等影響后的設(shè)計洪水成果。本文選取1958年7月實(shí)測洪水過程為典型洪水過程（見圖3），以千年一遇分析結(jié)果為例，可根據(jù)32 301 m3/s與22 300 m3/s的比值同倍比放大得到設(shè)計洪水過程線（見圖3）。

使用放大后的洪水過程線，根據(jù)2.3節(jié)介紹的水庫運(yùn)用方式進(jìn)行調(diào)洪計算，得到小浪底水庫的出庫流量。將此流量過程與對應(yīng)陸渾、故縣、河口村水庫及小花區(qū)間來水的總和作為花園口站的洪水過程，該洪水過程的最大值即為考慮小浪底等水庫影響后對應(yīng)量級設(shè)計洪水花園口站的洪峰流量。最終，計算得到考慮小浪底等水庫運(yùn)用后花園口站100 a一遇設(shè)計洪峰流量為12 948 m3/s，使用相同的方法計算得到考慮小浪底等水庫運(yùn)用后花園口站300 a一遇、1 000 a一遇設(shè)計洪峰流量分別為13 611、20 865 m3/s。將本文計算結(jié)果與已有文獻(xiàn)中的計算結(jié)果[7，23]進(jìn)行對比，見表4，發(fā)現(xiàn)本文結(jié)果比李保國等[7]計算得到的結(jié)果還小。

4 黃河下游各站設(shè)計洪水分析

采用陳贊廷等[24]編制的“黃河流域?qū)嵱盟念A(yù)報方案”中花園口站—夾河灘站、夾河灘站—高村站、高村站—孫口站的洪峰流量相關(guān)圖，根據(jù)花園口站的100 a一遇、300 a一遇和1 000 a一遇設(shè)計洪峰流量得到夾河灘站、高村站、孫口站的對應(yīng)設(shè)計洪峰流量，如圖4～圖6所示。

將上述結(jié)果與已有文獻(xiàn)中的計算結(jié)果[7，23]對比，見表5。本文方法所得孫口站計算結(jié)果比李保國等[7]計算結(jié)果還小，本文方法所得夾河灘站、高村站、孫口站計算結(jié)果明顯小于黃委[23]的計算結(jié)果。此外，張紅武等使用1919—1999年實(shí)測洪水還原資料，結(jié)合調(diào)洪計算，得到小浪底等水庫影響前后花園口站的設(shè)計洪水成果，再進(jìn)一步利用下游各站洪峰流量相關(guān)關(guān)系與以往所開展的黃河下游模型試驗(yàn)經(jīng)驗(yàn)，給出黃河下游各站設(shè)計洪水成果[9]，同本文方法所得結(jié)果頗為接近，只是本文方法所得夾河灘、高村等站的洪峰流量更小一些。

著名治黃專家陳先德在指導(dǎo)本文研究工作時指出，大型工程對黃河水沙變化的影響具有階躍特性，例如龍羊峽與劉家峽水利樞紐使托克托站汛期流量由2 000 m3/s降至500 m3/s左右，基流減少使北干流的洪水峰高量小，洪水進(jìn)入干流迅速坦化，龍門很難出現(xiàn)20 000 m3/s以上洪水流量，洪水流量基本在10 000 m3/s以下?；▓@口站流量超過10 000 m3/s的萬年一遇最大洪量為26億m3，完全能被現(xiàn)有工程和桃花峪工程攔蓄。

綜合本文結(jié)果，小浪底等水庫運(yùn)用后黃河下游各水文站設(shè)計洪峰流量都有所減小，自然意味著現(xiàn)有堤防設(shè)防標(biāo)準(zhǔn)都相應(yīng)提高，例如，由表5數(shù)據(jù)估計出花園口現(xiàn)狀設(shè)防流量22 000 m3/s相應(yīng)的洪水頻率可達(dá)1 560 a一遇，而按原設(shè)計該流量洪水頻率為“近千年一遇”，表明黃河下游防洪安全確實(shí)存在一定余地。若再進(jìn)一步修建桃花峪水庫，可在黃河下游發(fā)生千年一遇洪水時不使用北金堤滯洪區(qū)，徹底解放北金堤滯洪區(qū)[9]，且可在黃河發(fā)生300 a一遇以上洪水時，明顯減小東平湖滯洪區(qū)的分洪概率和分洪量，從而為黃河下游成為造福山東、河南人民的幸福河提供前提條件。與此同時，在黃河下游10 000 m3/s流量實(shí)際出現(xiàn)頻率已經(jīng)較小的條件下，清華大學(xué)與黃河勘測規(guī)劃設(shè)計研究院有限公司等單位共同承擔(dān)的國家重大科技研發(fā)計劃專項(xiàng)“黃河下游河道與灘區(qū)治理研究”提出了“兩道防線”布局[25]與生態(tài)治理、因?yàn)┦┲蝃26]相結(jié)合的“三灘分治”方案，即防護(hù)堤按高標(biāo)準(zhǔn)修建，作為保障“二灘”和“嫩灘”區(qū)域的相對穩(wěn)定的“第一道防線”，從而將從黃河大堤到主槽的灘地依次分區(qū)改造為“高灘”“二灘”和“嫩灘”[27]。以該方案實(shí)施黃河下游河道和灘區(qū)綜合提升治理工程，可實(shí)現(xiàn)黃委的“寬河固堤、穩(wěn)定主槽、因?yàn)┦┲?、綜合治理”思路和目標(biāo)[28]，使黃河下游形成生態(tài)保護(hù)和高質(zhì)量發(fā)展的嶄新局面。

5 結(jié) 論

以往分析設(shè)計洪水頻率曲線時，多先利用稀遇洪水的歷史調(diào)查洪水點(diǎn)據(jù)與部分實(shí)測資料將頻率曲線走向確定后，再調(diào)整CS/CV進(jìn)行經(jīng)驗(yàn)適線，不僅人為性強(qiáng)，而且存在“即使資料長度不斷增加，成果也不會產(chǎn)生大的變化”的缺陷。為此，本文不采用歷史調(diào)查洪水點(diǎn)據(jù)，利用花園口站1919—1949年歷年最大洪峰流量資料、1950—2010年實(shí)測洪峰流量還原資料與2011—2020年最大流量資料，對黃河下游花園口站的設(shè)計洪水頻率及洪水過程進(jìn)行計算分析。結(jié)果表明，天然狀態(tài)下花園口站100 a一遇設(shè)計洪峰流量為22 949 m3/s，考慮小浪底等水庫影響后花園口站100 a一遇設(shè)計洪峰流量為12 948 m3/s，結(jié)合洪峰流量相關(guān)圖計算出黃河下游孫口站對應(yīng)量級洪峰流量為10 100 m3/s。本文所得設(shè)計洪水成果比其他學(xué)者推薦修訂成果小，人為性小且適用于當(dāng)前與未來水沙情勢。小浪底等水庫運(yùn)用后黃河下游各水文站設(shè)計洪峰流量都有所減小，花園口現(xiàn)狀設(shè)防流量22 000 m3/s相應(yīng)的洪水頻率由“近千年一遇”減小到1 560 a一遇，意味著現(xiàn)有堤防設(shè)防標(biāo)準(zhǔn)都相應(yīng)提高，有利于按照“兩道防線”與生態(tài)治理、因?yàn)┦┲蜗嘟Y(jié)合的“三灘分治”模式，實(shí)施黃河下游河道和灘區(qū)綜合提升治理工程。

參考文獻(xiàn)：

[1] 李國英.黃河治理的終極目標(biāo)是“維持黃河健康生命”[J].人民黃河，2004，26（1）：1-2，46.

[2] 張紅武.未來黃河下游治理的主要對策[J].人民黃河，2004，26（11）：5-7，46.

[3] 史輔成，易元俊，慕平.黃河歷史洪水調(diào)查、考證和研究[M].鄭州：黃河水利出版社，2002：95-115，211-225.

[4] 謝悅波，李致家.頻率計算加入古洪水資料后對設(shè)計洪水的作用[J].河海大學(xué)學(xué)報，1995，23（6）：99-103.

[5] 金光炎.水文頻率分析述評[J].水科學(xué)進(jìn)展，1999，10（3）：319-327.

[6] 郭生練，劉章君，熊立華.設(shè)計洪水計算方法研究進(jìn)展與評價[J].水利學(xué)報，2016，47（3）：302-314.

[7] 李保國，宋偉華，畢黎明，等.黃河中下游設(shè)計洪水成果修訂及其影響研究[J].人民黃河，2019，41（12）：30-34.

[8] 黃河水利委員會革委會規(guī)劃大隊(duì).黃河流域歷史洪水調(diào)查考證及計算工作的體會[R].鄭州：黃河水利委員會，1974：1-2.

[9] 張紅武，張羅號，景喚，等.山東對黃河流域生態(tài)保護(hù)和高質(zhì)量發(fā)展的作用不可替代[J].水利水電技術(shù)（中英文），2021，52（1）：1-21.

[10] 黃春長，李曉剛，龐獎勵，等.黃河永和關(guān)段全新世古洪水研究[J].地理學(xué)報，2012，67（11）：1493-1504.

[11] 朱向鋒，黃春長，萬紅蓮，等.黃河中游馬蓮河歷史與現(xiàn)代洪痕沉積與水文學(xué)研究[J].干旱區(qū)地理，2012，35（5）：764-771.

[12] 胡貴明，黃春長，周亞利，等.伊河龍門峽段全新世古洪水和歷史洪水水文學(xué)重建[J].地理學(xué)報，2015，70（7）：1165-1176.

[13] HU G， HUANG C C， ZHOU Y， et al. Hydrological Studies of the Historical and Palaeoflood Events on the Middle Yihe River， China[J]. Geomorphology， 2016， 274：152-161.

[14] 張紅武，張羅號，彭昊，等.沖積河流糙率由來與計算方法研究[J].水利學(xué)報，2020，51（7）：774-787.

[15] 張紅武，江恩惠，趙新建，等.黃河高含沙洪水模型的相似條件[J].人民黃河，1995，17（4）：1-3，61.

[16] 張紅武，張清.黃河水流挾沙力的計算公式[J].人民黃河，1992，14（11）：7-9，61.

[17] 黃河水利委員會勘測規(guī)劃設(shè)計院.黃河流域洪水調(diào)查資料[R].鄭州：黃河水利委員會，1983：128-129.

[18] 蔡蓉蓉.黃河水沙特性研究[D].北京：清華大學(xué)，2020：34-46.

[19] 中華人民共和國水利部.水利水電工程設(shè)計洪水計算規(guī)范：SL 44—2006[S].北京：中國水利水電出版社，2006：7-8.

[20] 詹道江，葉守澤.工程水文學(xué)[M].北京：中國水利水電出版社，2000：181-200.

[21] 吳江，翁澤平.用EXCEL繪制PⅢ頻率曲線[J].工程地質(zhì)計算機(jī)應(yīng)用，2005（1）：17-19.

[22] 胡一三，繆鳳舉，曲少軍，等.小浪底水庫運(yùn)用初期三門峽水庫防洪運(yùn)用分析[J].泥沙研究，2001，26（2）：25-28.

[23] 黃河水利委員會.黃河流域防洪規(guī)劃[M].鄭州：黃河水利出版社，2008：44-47.