黃河下游近期來水來沙年內年際波動特征研究

關鍵詞：水沙條件；趨勢性；水沙搭配；游蕩段；黃河下游

0引言

21世紀以來，受黃河水沙變化與小浪底水庫調節的影響，進入黃河下游的水沙特征發生了新的變化，由此引發了下游河道尤其是游蕩河段治理的新問題，例如：局部河段河勢大幅度調整、多河段出現畸形河勢等，成了黃河下游防洪安全的潛在風險點。在深入推動黃河流域生態保護和高質量發展的關鍵時期，氣候變化與人類活動影響雙重作用下黃河下游游蕩段水沙變化新情勢亟待明晰。

按照不同的河床演變特點，黃河下游通常可分為游蕩、過渡和彎曲3個河段，其中游蕩段有花園口、夾河灘與高村3個水文站。黃河下游來水來沙量不僅受自然條件影響，而且受人類活動的影響也越發強烈。截至2004年，黃河流域共有超過3380個水庫建成運營，其中小浪底水庫1999年下閘蓄水后，2002年開始實施水庫群聯合調度，這在很大程度上改變了進入黃河下游的水沙條件。近年來，國內外許多學者基于不同水文站實測水沙數據，對進入黃河下游的水沙過程開展了大量研究。例如：王鴻翔等利用Mann-Kendal檢驗法和小波分析法對花園口站1960-2016年徑流量和輸沙量進行分析，發現水、沙量總體都呈下降趨勢，其突變點分別出現在1985年和1998年；申紅彬等基于黃河下游花園口、高村、艾山、利津4個水文站建立下游水沙關系模型，采用理論分析法、線性回歸法得出符合實際的河道沿程水沙搭配參數調整方程，揭示了水沙關系由不平衡向平衡的過渡過程：韓沙沙等基于1950-2018年利津站實測水沙資料，分析了利津站水沙變化特征，認為黃河口來水來沙過程逐漸由自然模式轉變為人為調控模式。

為進一步探究近期黃河下游游蕩段來水來沙年內年際波動特征，重點關注1997年小浪底壩址大河截流后強人類活動影響階段，本文基于花園口水文站1997-2020年實測水沙數據，并對比同時段潼關站水沙數據，采用Mann-Kendall（M-K）檢驗法、滑動t檢驗法和小波分析法分析兩站水沙序列的趨勢性、變異性、周期性、年內分布及水沙搭配特點，以期為黃河下游河道演變等相關科學研究提供支撐，也為未來水沙調控等工程實踐提供決策參考。

1趨勢性與變異性

本文以花園口水文站實測水沙數據代表黃河下游水沙條件，其水沙過程直接受小浪底水庫運用的影響。相對而言，潼關水文站雖然也受黃河上中游水庫及水保工程等影響，但其水沙過程更能反映黃河天然水沙變化特別是泥沙的變化，因此本文選用潼關站水沙數據代表黃河天然水沙過程，對比花園口水沙過程，分析人工調控作用下進入黃河下游的水沙特征。

潼關、花園口站1997-2020年年來水量、來沙量變化如圖1所示。由圖1（a）可知，兩站來水量稍有差異，但其整體變化趨勢基本一致，均呈現階段性增減交替的過程，這是因為小浪底水庫為年調節水庫，年內水量基本平衡，所以黃河下游來水量年際變化與天然來水量過程基本一致。

潼關、花園口站來沙量變化差異明顯。由圖1（b）可知，潼關站來沙量除2000-2001年較小外，1997-2009年總體波動下降，2009-2013年呈波動上升趨勢，2013-2015年呈下降趨勢，隨后除2018年來沙量較大外，2015-2020年總體呈上升趨勢。花園口站來沙量1997-1999年較大，平均為4.09億t，但仍較潼關站來沙量偏少28%：1999年小浪底水庫投入運用后，進入下游的沙量銳減，潼花（潼關一花園口）間的來沙量差距加大：2002年調水調沙開始實施后至2009年，花園口站來沙量不斷減少，這一方面是因為黃河天然來沙量呈下降趨勢，另一方面也得益于小浪底水庫的攔沙作用，但兩站來沙量差距逐漸減小；2009-2013年受天然來沙量增加趨勢的影響，雖然有小浪底水庫的攔截，但是花園口來沙量仍呈小幅度增加趨勢：2014-2017年雖然天然來沙量呈增加趨勢，但是來沙量相對較少，在小浪底水庫攔沙作用下，花園口來沙量仍呈下降趨勢：2018-2020年花園口站年均來沙量較大，為1.49t，占潼關站來沙量的57%。

總體來看，近20a來黃河下游水沙通量變化顯著，人工調控作用下的黃河下游來水量與天然來水量變化過程基本一致，而來沙量則受天然來沙與小浪底水庫運用的共同影響。

為進一步科學描述黃河下游水沙變化規律，采用M-K檢驗法與滑動t檢驗法檢驗潼關、花園口站年均水、沙量的趨勢性與突變點。M-K檢驗法是一種非參數檢驗法，對所有的分布均適用，部分數據缺失也不會影響結果，常用于對降水、徑流、氣溫和水質等要素時間序列變化趨勢和突變點分析：滑動t檢驗法則是通過兩組樣本平均值差異是否顯著來檢驗突變，要求變量連續且有一定的正態性。檢驗結果表明，兩站在1997-2020年來水量呈顯著增加趨勢、來沙量呈顯著減少趨勢。潼關站水量在2018年附近發生突變，花園口站水量在2017年發生顯著突變，兩站突變年份相近；潼關站沙量在2004年附近發生突變，花園口站沙量在2001-2002年附近發生顯著突變。這說明小浪底水庫的運用方式對黃河下游來水量變化過程影響不大，但直接影響來沙量的變化過程。

2周期性

河流徑流與輸沙是一類受不同因素影響、十分復雜的多時間尺度變化過程，小波分析法常被用于識別此類時間序列中的多種變化周期，例如河川徑流、地震波、暴雨、洪水等。本文采用Morlet連續復小波變換的小波分析法識別1997-2020年黃河下游來水來沙量的周期性變化特征。

結合Matlab、Suffer等軟件計算并繪制潼關、花園口站來水量、來沙量時間序列的小波系數實部等值線圖、小波方差圖及主周期變化趨勢圖。小波系數實部等值線圖可反映序列在不同時間尺度的周期變化與在時間域中的分布，小波方差圖可反映序列波動能量隨時間尺度的變化，主周期變化趨勢圖可反映序列在不同時間尺度存在的平均周期。

小波系數實部等值線如圖2所示。小波系數實部為正代表水量或沙量充沛，圖中以實線、紅底表示，“H”為正值中心：小波系數實部為負則代表水量或沙量處于枯期，以虛線、藍底表示，“L”為負值中心。由圖2可識別出潼關、花園口站來水量在演變過程中不同時間尺度的變化周期，潼關站為3～7、7～11、12～16a，花園口站為3～8、9～16a，相對于潼關站，花園口站來水量兩個時間尺度的周期變化在整個分析時段更為穩定，具有全域性；就來沙量而言，潼關站來沙量在4～9、10～16a兩個時間尺度的周期變化較為穩定，具有全域性，而花園口站來沙量識別出的3～8、9～16a兩個時間尺度的周期變化在2002年調水調沙實施前相對穩定。

計算潼關、花園口站水、沙量小波方差，如圖3所示。由圖3可知，潼關、花園口站來水量最大峰值即第一主周期分別為14、13a，來沙量第一主周期分別為13、12a。根據潼關、花園口站水、沙量小波方差最大峰值可繪制主周期變化趨勢圖，如圖4所示。潼關站水、沙量變化平均周期為9.5a，研究時段內經歷了2.5個周期的豐一枯變化，花園口站水、沙量變化平均周期為8a，經歷了約3個周期的豐一枯變化。

總體來說，相對于潼關站水沙過程，花園口站來水量周期性變化在整個分析時段較為穩定，具有全域性，來沙量在2002年調水調沙后在小時間尺度的周期性變化相對不穩定，但兩站水沙的主周期變化趨勢仍較為一致。這說明人工調控作用下的水沙過程仍能反映天然水沙的周期性變化，但小時間尺度來沙過程受水庫運用方式影響較大。

3年內分布

基于潼關、花園口站年均水沙變化過程，結合趨勢性、變異性及周期性分析，將研究時段劃分為1997-1999年、2000-2001年、2002-2012年、2013-2017年、2018-2020年共5個階段。為量化黃河下游水沙年內分布特征，用變差系數CV表示兩站水沙年內分布不均勻度，CV越大表明水沙年內分布越不0MEP1NNkeBHfG4JuFXo0hw==均勻，波動幅度越大。

潼關、花園口站月均流量、含沙量年內分布變差系數計算結果見表1。由表1可知兩站流量變差系數CV相近，5個階段均經歷了減小一增大的交替過程。兩站含沙量變差系數有所差別，潼關站5個階段的含沙量變差系數均較為接近，花園口站各階段含沙量則直接受小浪底水庫運用方式的影響，其年內分布不均勻度變化較大。此外，兩站含沙量變差系數大于流量變差系數，含沙量年內分布的不均勻性更為明顯。

4水沙搭配

采用來沙系數t與水沙系數K兩個經驗參數表示水沙搭配狀況。來沙系數t是指單位流量的含沙量大小，一般用于描述河道輸沙能力與水沙關系的變化：水沙系數K是一種時段水沙搭配判定指標，可定性判斷研究年份的水沙條件相對于多年平均水沙條件的優劣。

潼關、花園口站來沙系數t、水沙系數K變化如圖5所示，各水沙參數階段性變化見表2。由圖5和表2可知，潼關站來沙系數t第1階段減小，第Ⅱ階段增大，均在0.02以上，水沙系數K較小，前兩階段水少沙多，水沙搭配惡劣；第Ⅲ階段t穩步減小，較上一階段顯著減小約53%，K不斷增大，且階段后期增大速度較快，階段平均值為1.48，該階段來沙量大幅減少，來水量增加，水沙搭配不斷優化；第Ⅳ階段t呈增大趨勢，K呈減小趨勢，但該階段水沙搭配仍持續優化，減小至0.0055，K增大至2.21;第V階段t繼續減小，K快速增大，該階段來沙量與上一階段持平，來水量大幅增加，水沙搭配顯著優化。

相對于黃河天然水沙過程，黃河下游的水沙搭配較潼關站稍有差異。花園口站第1階段來沙系數t最大，為0.0220，水沙系數K較小，為0.12，該階段黃河下游水沙條件與潼關站相似，水少沙多，水沙搭配惡劣；第Ⅱ階段t顯著減小約62%，K仍較小，為0.09，該階段小浪底水庫開始運用，進入黃河下游的沙量顯著減少，水量雖仍與潼關站持平，但由于汛期水庫蓄水，最大流量被削峰，因此雖然水沙搭配優化，但水沙系數仍維持在較小值；隨后的第Ⅲ、Ⅳ階段，t繼續減小，K快速增大，調水調沙成效顯著，使黃河下游的水沙搭配得到了明顯改善；第V階段t明顯回升，K有所減小，該階段水沙搭配仍較優。

對比分析花園口站與潼關站兩個經驗參數的變化過程可知，相對于潼關站天然水沙過程來沙系數t逐階段減小、水沙系數K逐階段增大的變化特征，花園口站來沙系數t前4個階段逐漸減小，第V階段有所回升，水沙系數K前2個階段維持在較小值，Ⅲ一Ⅳ階段逐漸增大，第V階段有所減小。這說明，花園口水沙搭配變化一方面直接受小浪底水庫運用方式的影響，另一方面又反映了天然來水來沙過程。總體來看，在天然水沙條件變化與人工水沙調控綜合作用下，近20a來黃河下游水沙條件整體相對較優。

5結論

以同時段潼關站水沙過程代表黃河天然來水來沙條件，對比分析1997-2020年黃河下游花園口站水沙年際年內變化特征可知：

1）兩站來水來沙量變化的趨勢性基本一致，來水量均呈顯著增加趨勢、來沙量均呈顯著減少趨勢，但變異性稍有差別，潼關站來水量在2018年附近發生突變，來沙量在2004年附近發生突變，花園口站水、沙量發生顯著突變的時間分別為2017年、2001-2002年。

2）兩站水、沙量存在不同時間尺度的周期變化，但其水沙的主周期變化較為一致，潼關站水、沙量第一主周期分別為14、13a，平均周期為9.5a，花園口站水、沙量第一主周期分別為13、12a，平均周期為8a。

3）基于兩站來水來沙量變化過程，結合趨勢性、變異性及周期性分析，可將研究時段劃分為5個階段。5個階段兩站流量變差系數CV相近，均經歷了減小一增大的交替過程，含沙量變差系數則有所差別，潼關站5個階段的含沙量變差系數均較為接近，花園口站各階段含沙量則直接受小浪底水庫運用方式的影響，其年內分布不均勻度變化較大。

4）相對于潼關站天然水沙過程來沙系數t逐階段減小、水沙系數K逐階段增大的變化特征，花園口站來沙系數t前4個階段段漸減小，第V階段有所回升，水沙系數K前2個階段維持在較小值，Ⅲ一Ⅳ階段逐漸增大，第V階段有所減小。

5）小浪底水庫的運用方式對黃河下游來水過程影響不大，但直接影響來沙過程。在天然水沙條件變化與人工水沙調控綜合作用下，近20a來黃河下游水沙過程呈明顯的階段性變化特征，水沙搭配總體相對較優。