張治昊，戴 清，李敬義，聶莉莉

（1.中國水利水電科學研究院，北京100044；2.黃河河口管理局，東營257091）

黃河口是黃河水沙的歸宿區(qū)，是河海動力的交匯區(qū)，是黃河流域系統(tǒng)的重要組成部分。人民治黃以來，黃河口的研究和治理工作得到了國家的高度重視，取得了豐碩的成果［1］。黃河口是世界公認的最復雜的河口，老問題未完全解決，新問題又不斷涌現(xiàn)。1986年后，黃河口水沙過程發(fā)生顯著變異，對河口演變造成巨大影響。黃河口尾閭河道作為河海交匯作用最劇烈的地帶，為適應水沙過程變異，經歷了復雜的響應過程；其間，尾閭河道萎縮的特征最為突出［2］。尾閭河道的萎縮，對黃河防洪安全構成極大威脅，全面深入地研究水沙過程變異造成的尾閭河道萎縮迫在眉睫。本文在深入分析黃河口水沙過程變異的基礎上，揭示了黃河口尾閭河道斷面形態(tài)的萎縮特征，探求了尾閭河道斷面形態(tài)調整與黃河口水沙過程的響應關系，引入了尾閭河道斷面形態(tài)萎縮的判別指標，試圖為科學治理黃河口提供技術支持。

1 黃河口水沙過程變異特征

黃河口行水清水溝流路以來，水沙過程發(fā)生顯著變異，除了自然水文情況影響外，主要影響因素還是人類的活動，如大型水利工程的建設和運用會劇烈改變河流下游的來水來沙過程［3］。1986年，黃河上游龍羊峽水庫建成并投入運用，對下游來水來沙過程影響巨大，是黃河口水沙過程發(fā)生變異的根本原因；2002年，小浪底水庫開始調水調沙運用，也使進入下游河道的水沙過程變化較大，綜合上述分析，以1986年、2002年為界，將黃河口清水溝流路水沙過程分為1976～1985年、1986～2001年、2002～2006年3個時段，以I時段為基準階段，用II、III時段與之進行對比，分析黃河口水沙過程變異特征。

1.1 年際間變異

由黃河口利津站水文特征值統(tǒng)計（表1）可見，黃河口水沙過程變異特征主要表現(xiàn)為：（1）來水來沙大幅度減少，II時段河口年均來水量為137.7億m3，僅為I時段的40.6%，年均來沙量也由I時段的8.31億t減少到3.51億t；（2）水沙搭配關系惡化，I時段河口年均來沙系數為0.027，II時段高達0.061，是I時段的2.26倍，表明II時段黃河口水沙搭配關系急劇惡化，為小水帶大沙時段；（3）來沙組成有所變細，但變化幅度較小，河口年均中值粒徑由I時段0.02 mm變?yōu)镮I時段0.018 mm；（4）利津以下河道水沙量沿程大幅減少，利津以下有王莊等9處引水口，設計引水流量310 m3/s。1986年后，由于黃河口來水來沙減少，而工農業(yè)發(fā)展的需水量增大，利津以下河道引水量明顯增多，所以即便在利津站不斷流期間，利津至口門河段水沙量沿程減少也十分嚴重。（5）III時段黃河口水沙過程變異有所緩解，主要表現(xiàn)為來水量增大，來沙量減少，水沙搭配關系明顯改善，由于小浪底水庫調水調沙的持續(xù)運用，造成下游河道沖刷，河床粗化，進入河口的泥沙顆粒變粗。

表1 利津站水文特征統(tǒng)計Tab.1 Characteristic values of hydrology at Lijin hydrological station

1.2 年內分配變異

由黃河口來水來沙年內分配過程（圖1）可見：與I時段相比，II時段黃河口水沙年內分配的豐枯歷時、量值大小、變化趨勢、極值、變化幅度均發(fā)生了變化；最突出的變化是II時段水沙量從8月峰值迅速減至10月枯值，改變了I時段水沙量從汛期峰值緩慢減至12月枯值，表明II時段水沙量8月到10月減小幅度最大。與前兩時段相比，III時段黃河口水沙年內分配明顯不同，最大的不同之處在于從6月份開始，水沙過程幾乎沒有峰值，來水來沙量值分別維持在25億m3、0.35億t左右，一直持續(xù)到10月汛期結束，充分反映出小浪底水庫調水調沙運用對黃河口汛期水沙過程有較強的控制作用。

圖1 黃河口來水來沙年內分配過程Fig.1 Processes of monthly distribution of flow and sediment in the Yellow River Estuary

1.3 典型水文特征變異

由黃河口典型水文特征統(tǒng)計（表2）可見：與I時段相比，II時段黃河口汛期洪峰流量減少劇烈，且流量級越大，減少的越多；同時，非汛期斷流次數不斷增多，斷流天數不斷增加，斷流長度不斷增長，III時段黃河口 Q≥3 000 m3/s、3 000＞Q≥2 000 m3/s、2 000＞Q≥1500 m3/s三個流量級天數均有增加，說明 III時段黃河口汛期洪水量級有所恢復。1999年后，黃河水利委員會采取了全流域水量統(tǒng)一調度的人工干預措施，利津站斷流的繼續(xù)惡化才得到遏制，但黃河口功能性斷流問題仍沒有得到徹底解決。

表2 黃河口典型水文特征統(tǒng)計Tab.2 Typical characteristic of hydrology in the Yellow River Estuary

2 黃河口尾閭河道斷面形態(tài)的萎縮調整

黃河口水沙過程變異后，為適應持續(xù)的枯水枯沙條件，尾閭河道萎縮十分明顯。表3為尾閭河道相應的3個典型斷面不同時段主河槽平灘特征均值統(tǒng)計，由表3可見，水沙過程變異，黃河口尾閭河道橫斷面形態(tài)變化特征主要表現(xiàn)為：（1）平灘寬度縮窄，水深變淺，面積大幅度減小。與I時段相比，利津、漁洼、清6三個斷面II時段主河槽平灘河寬均值分別減少了92 m、84 m、285 m；主河槽平灘河寬明顯縮窄的同時，3個斷面主河槽平灘水深均值也呈減小的趨勢，減小的幅度為0.3～1.05 m，由于主河槽平灘河寬和水深同時減少，尾閭河道3個斷面平灘面積均值大幅度減少，分別由I時段的2 047 m2、1 984 m2、2 137 m2減少至II時段的1 230 m2、1 644 m2、1 361 m2，減少幅度達20～39%。（2）平灘寬深比有所增大，主河槽更加寬淺。利津、漁洼、清6三個斷面主河槽平灘寬深比均值分別由I時段的5.53、10.2、19.1增加至II時段的6.67、11.63、20.1；主河槽平灘寬深比均值增大表明，尾閭河道萎縮過程中，橫斷面淤積的部位集中在斷面主河槽的窄深處，尾閭河道橫斷面呈寬淺碟型。（3）受黃河口尾閭河道清8出汊工程的影響，近海段橫斷面形態(tài)調整與尾閭上段有所不同。1996年黃河口尾閭河道實施了清8出汊工程，尾閭河道近海段長度大大縮短，比降相應增大，當年汛期適逢黃河“96.8”洪水的大流量水沙過程，尾閭河道清8斷面以上河段發(fā)生了長距離的溯源沖刷，由距離清8斷面較近的清6斷面統(tǒng)計資料明顯看出，1996年后，清6斷面河道主河槽有大幅度的沖深，1996～2001年時段其主河槽平灘水深均值比1986～1995年時段沖深1.01 m，其主河槽平灘面積均值增加148 m2，增大幅度為11%，雖然主河槽平灘面積增大幅度較小，但其主河槽斷面形態(tài)卻有較大程度的調整，其主河槽平灘寬深比均值由1996年前的20.1調整至9.57，說明1996年后近海段主河槽橫斷面形態(tài)變得十分窄深，從而反映出，通過調整口門，增大比降，造成溯源沖刷，對塑造近海段河道橫斷面形態(tài)效果明顯。（4）III時段尾閭河道主河槽萎縮程度有所恢復。3個典型斷面主河槽沖深十分明顯，主河槽平灘水深均值增大0.78～1.03 m，主河槽平灘面積均值增大137～267 m2，由于主河槽平灘河寬均值變化不大，而平灘水深明顯增大，所以3個典型斷面的主河槽平灘寬深比均值大幅度減小，說明該時段尾閭河道橫斷面形態(tài)均朝著窄深方向調整，小浪底水庫調水調沙的成功運用，對于尾閭河道主河槽過水斷面的恢復發(fā)揮了良好的作用［4］。

表3 黃河口尾閭河道典型橫斷面主河槽平灘特征值統(tǒng)計表Tab.3 Statistics of bankfull value of tail channels in the Yellow River Estuary

3 尾閭河道斷面形態(tài)調整與水沙過程的響應關系

3.1 尾閭河道平灘面積與水沙過程的關系

與其他特征參數相比，在表征尾閭河道橫斷面形態(tài)的特征參數中，尾閭河道平灘面積A隨黃河口水沙過程的變化最為靈敏，跟隨性良好。圖2點繪了尾閭河道利津斷面平灘面積A與黃河口年際水沙過程的關系，利用回歸分析法分別建立了尾閭河道利津、漁洼、清6三個斷面平灘面積A與黃河口來水量w、汛期來沙系數ρ的相關關系式（表4）。

圖2-a中關系線的變化趨勢表明，尾閭河道平灘面積A隨黃河口來水量w的增大而增大。黃河口來水量越大，其造床能力越強，尾閭河道淤積減輕，沖刷增強，平灘面積必然擴大。反之，黃河口來水量小，造床能力弱，尾閭河道淤積加重，平灘面積必然減小。圖2-b中關系線的變化趨勢表明，尾閭河道平灘面積A隨黃河口來沙系數ρ的增大而減小。黃河口來沙系數越大，說明相同來水過程，其水流含沙量越大，水沙搭配關系越差，河床淤積必然加重，尾閭河道平灘面積受到損失。

表4中式（1）～式（3）的相關系數 R 值分別為0.83、0.78、0.76，反映出尾閭河道平灘面積A與黃河口來水量w的關系非常密切；式（4）～式（6）的相關系數 R 值分別為 0.72、0.69、0.65，反映出尾閭河道平灘面積A與黃河口來沙系數ρ的關系比較密切，但與來水量w相比，兩者的關系明顯偏弱。仔細觀察表4中式（1）～式（6）的相關系數R值還可發(fā)現(xiàn)，3個斷面平灘面積A與黃河口水沙過程2個參數的關系R值，都是呈沿程減小的趨勢，說明越靠近河口，受到海洋動力、邊界條件等其他因素影響越大，尾閭河道演變與黃河口來水來沙過程的關系沿程逐漸減弱。

圖2 黃河口尾閭河道利津斷面平灘面積與水沙過程關系圖Fig.2 Relationship between bankfull area and flow-sediment process of Lijin section

表4 黃河口尾閭河道典型橫斷面平灘面積與汛期水沙過程關系統(tǒng)計Tab.4 Statistics of the relationship between bankfull area and flow-sediment process of tail channels

3.2 尾閭河道平灘寬深比與水沙過程的關系

式（7）～式（8）的相關系數R值分別為0.61、0.63，反映出尾閭河道平灘寬深比與黃河口水沙過程的關系比尾閭河道平灘面積與黃河口水沙過程的關系要明顯偏弱。相比而言，尾閭河道平灘寬深比與黃河口來沙系數ρ的關系比尾閭河道平灘寬深比與來水量w的關系更密切，說明黃河口水沙搭配關系對尾閭河道橫斷面形態(tài)的塑造非常重要，同時來水量的大小對尾閭河道橫斷面形態(tài)的塑造也不容忽視。

圖3 黃河口尾閭河道利津斷面平灘寬深比與水沙過程關系圖Fig.3 Relationship between bankfull width-depth ratio and flowsediment process of Lijin section

4 尾閭河道斷面形態(tài)萎縮的判別指標

黃河口尾閭河道萎縮導致相同洪水流量尾閭河道水位急劇升高，黃河防洪安全受到極大威脅，已引起政府部門的高度重視［6］。1988年至今，國家多次對黃河口尾閭河道實施了大型拖淤、挖河疏浚等工程措施，對減緩尾閭河道萎縮起到了一定作用，但工程實施過程中，急需引入一個判別指標，作為實施尾閭河道萎縮治理工程的控制標準。

利津斷面是黃河口尾閭河道演變的進口控制斷面，水文資料較齊全，可作為黃河口尾閭河道萎縮的控制斷面。為了對尾閭河道萎縮進行量化控制，首先要提出一個能確切表征斷面形態(tài)萎縮的變量。尾閭河道斷面形態(tài)萎縮最突出的特征是過水斷面面積損失嚴重，河道過流能力減小，所以選用利津斷面的變量是利津斷面第i年過水面積損失占多年平均過水面積的百分比，其中ΔAj為利津斷面第i-1年的過水面積減去第i年過水面積的差值；為利津斷面多年過能比較確切的表征尾閭河道第i年的萎縮程度。變量水面積的平均值。當變量為正值時越大，利津斷面第i年過水面積損失程度越大；變量越小，利津斷面第i年過水面積損失程度越小。

影響黃河口尾閭河道萎縮的因素主要有2方面：（1）來水來沙過程：黃河口水沙過程變異是尾閭河道萎縮的主要原因；（2）前期河床邊界條件：尾閭河道前期斷面形態(tài)與河道萎縮關系密切，前期萎縮的斷面形態(tài)將加劇尾閭河道的萎縮。根據實測資料，運用多元回歸方法，建立第i年黃河口尾閭河道斷面形態(tài)萎縮變量與第i年水沙過程變異參數ξ和利津斷面第i-1年寬深比的關系式

其綜合相關系數R=0.82，表明相關性良好，將水沙過程變異參數ξ的表達式［2］代入式（9）可得細化后的表達式

由式（11）可知，尾閭河道斷面形態(tài)萎縮判別指標Tn是變量自1986年后至第n年的累加值，其物理意義十分明確：1986年后，尾閭河道萎縮發(fā)展到第n年過水斷面面積累計損失的百分比。圖4為尾閭河道斷面形態(tài)萎縮判別指標Tn隨年份變化過程，由圖4可見，1986～1995年，除1990年外，Tn持續(xù)增大，由1986年的0.21增大到1995年的0.5；1996～2001年，Tn有減有增，出現(xiàn)一定波動；2002～2006年，Tn明顯呈減小的趨勢。Tn的變化反映出，1986～1995年尾閭河道萎縮隨時間不斷發(fā)展加重，最嚴重的1995年過水斷面面積累計損失達50%；1996～2001年，受清8出汊溯源沖刷的影響，尾閭河道過水斷面面積有所恢復，但萎縮的趨勢并未改變；2002～2006年，小浪底水庫調水調沙的成功運用，明顯減輕了尾閭河道萎縮程度。

由圖4可見，Tn分布在直線Tn=0.2之下的為0 a；分布在直線Tn=0.2和直線Tn=0.3之間的有4a，約占圖4中總年數20%；分布在直線Tn=0.3和直線Tn=0.45之間的有13 a，約占圖4中總年數60%；分布在直線Tn=0.45之上有4 a，占圖4中總年數20%。根據圖4中點群上述分布特征，可取Tn值為0.2、0.3、0.45，對河道萎縮發(fā)展過程進行階段劃分：Tn＜0.2，尾閭河道尚未出現(xiàn)萎縮階段；0.2≤Tn＜0.3，尾閭河道萎縮出現(xiàn)和發(fā)育階段；0.3≤Tn＜0.45，尾閭河道萎縮形成和發(fā)展階段；Tn≥0.45，尾閭河道萎縮十分嚴重階段。

判別指標Tn對尾閭河道萎縮發(fā)展過程的階段劃分，為治理尾閭河道萎縮工程的實施提供了量的控制。當萎縮判別指標Tn達到0.2時，標志著尾閭河道已開始出現(xiàn)萎縮，應加強觀察，密切注視河道萎縮是否繼續(xù)發(fā)展；Tn升高到0.3時，有條件的前提下，應采取一定的工程措施減小河道萎縮程度；Tn升高到0.45時，必須采取有力的工程措施或綜合性工程措施大幅度減小河道萎縮程度。

圖4 黃河口尾閭河道斷面形態(tài)萎縮判別指標Tn變化過程Fig.4 Variational processes of distinguishing target Tnof atrophy of tail channels in the Yellow River Estuary

5 結論

（1）1986年后，黃河口年際間水沙劇烈減少，水沙搭配關系變壞；年內分配以汛期8～10月減少的幅度最大；汛期洪水與非汛期小水的對比發(fā)生質的變化，典型水文條件兩極分化；河道水沙沿程衰減嚴重，入海水量很少，甚至為0；黃河口水沙過程變異是上述各種水沙條件變異綜合作用的結果。2002年后，由于小浪底水庫調水調沙的持續(xù)運用，黃河口水沙過程變異有所緩解，來水量增大，來沙量減少，水沙搭配關系明顯改善。

（2）黃河口水沙過程變異，尾閭河道斷面形態(tài)萎縮十分明顯，主要表現(xiàn)為平灘寬度縮窄，平灘水深變淺，平灘面積大幅度減小，平灘寬深比有所增大，斷面形態(tài)變的更加寬淺。2002年后，尾閭河道斷面萎縮有所恢復，平灘水深與平灘面積明顯增大，平灘寬深比大幅度減小，斷面形態(tài)朝著窄深方向調整，說明小浪底水庫調水調沙運用對于尾閭河道斷面形態(tài)的恢復發(fā)揮了良好的作用?；貧w分析表明，尾閭河道斷面形態(tài)調整與水沙過程有著良好的響應關系，尾閭河道平灘面積隨黃河口來水量的增大而增大；隨黃河口來沙系數的增大而減?。黄綖捝畋入S黃河口來水量的增大而減?。浑S黃河口來沙系數的增大而增大。

（3）黃河口造床流量的減小、河床阻力的增大、河海耦合動力減小均降低了尾閭河道的輸水輸沙能力，其持續(xù)作用的結果是塑造了萎縮的河槽水力幾何形態(tài)；而萎縮后的河床形態(tài)反影響于尾閭河道的輸水輸沙能力，其結果是加劇了河道萎縮的發(fā)展。為解決治理尾閭河道萎縮缺乏控制標準的問題，引入尾閭河道斷面形態(tài)萎縮的判別指標Tn，Tn的物理意義是1986年后尾閭河道橫斷面累計損失的過水面積百分比，依據Tn點群分布特征，取Tn值0.2、0.3、0.45作為實施尾閭河道萎縮治理工程的控制界限點。

［1］胡春宏.黃河水沙過程變異及河道的復雜響應［M］.北京：科學出版社，2005.

［2］胡春宏，陳建國，郭慶超，等.黃河水沙調控與下游河道中水河槽塑造［M］.北京：科學出版社，2007.

［3］胡春宏，張治昊.黃河口尾閭河道平灘流量與水沙過程響應關系［J］.水科學進展，2009，20（2）：209-214.HU C H，ZHANG Z H.Response relationship between the bank-full discharge in tail channel and the process of flow-sediment in the Yellow River Estuary［J］.Advances in Water Science，2009，20（2）：209-214.

［4］王光謙，劉家宏，李鐵鍵.黃河數字流域模型原理［J］.應用基礎與工程科學學報，2005，13（1）：1-8.WANG G Q，LIU J H，LI T J.Digital Watershed Model of Yellow River［J］.Journal of Basic Science and Engneering，2005，13（1）：1-8.

［5］錢寧，張仁，周志德.河床演變學［M］.北京：科學出版社，1989.

［6］張治昊，胡春宏.黃河口水沙過程變異及其對河口海岸造陸的影響［J］.水科學進展，2007，18（3）：336-341.ZHANG Z H，HU C H.Variation of the processes of flow-sediment and its effect on epeirogenesis of seacoast in the Yellow River Estuary［J］.Advances in water science，2007，18（3）：336-341.