長江中下游干流河道，宜昌至湖口（955 km）為中游，湖口至大通（338km）為下游，大通以下為河口段（600km）。據(jù)大通站多年實測資料統(tǒng)計，流域平均每年匯集河道的徑流總量達9.051×1011m3（1950～2000年），并挾帶約4.33×108t（1950～2000年 ）泥沙入海［1］。巨量的水沙下泄不僅對長江河口，而且對臨近陸架海域的水文、泥沙、沉積和生態(tài)環(huán)境等都 有重大影響，多年來許多研究工作者已經(jīng)從不同的角度對這些影響進行了討論［2～5］。本文在前人的研究成果和大量實測資料基礎(chǔ)上，采用水文學(xué)、泥沙運動力學(xué)和數(shù)理統(tǒng) 計學(xué)相結(jié)合的分析方法，著重研究長江中下游河道的水沙通量的變化規(guī)律，以及對河口水沙通 量的貢獻，這對于進一步探討河口發(fā)育及演變的成因?qū)⒕哂猩钸h的意義。

　　1 水沙來源

　　長江中下游河道水量及其變化規(guī)律主要取決于流域降雨的時空分布，前人曾做過長江中下游各站年平均流量與上游流域年降雨量的比較，結(jié)果表明峰谷起伏基本相應(yīng)［4］。根據(jù)統(tǒng)計資料，宜昌站多年平均徑流量為4.382×1011m3（1950～2000年），與 入河口區(qū)（大通站） 徑流總量9.051×1011m3（1950～2000年）相比［1］，只占其48.4%，由此可見，進入河口區(qū)徑流量有一半來自上游，一半來自中下游流域。

　　長江中下游流域?qū)儆跊_積平原河流，湖泊對于水沙的調(diào)蓄和沉積作用不容忽視，這對于河道輸沙量和含沙量的沿程變化具有重要影響，洞庭湖多年平均（1956～1995）入湖沙量為1.67×108t，其中有1.32×108t泥沙通過荊江河段的松滋口、太平口和藕池口（其中 調(diào)弦口已于1959建閘）輸入洞庭湖［6］；含沙量在宜昌以下呈單向遞減的趨勢，中游宜昌站多年平均含沙量為1.14 kg/m3（1950～2000年），下游大通站的多年平均含沙量僅為0.486kg/m3（1950～2000年）［1］，只相當于宜昌站的42.6%（表1）；同時長江上游來沙量的變化對于中下游 輸沙規(guī)律的變化也起著重要的作用，上的金沙江和嘉陵江是長江流域的主要產(chǎn)沙區(qū)，兩江多年（1954～1998年）平均輸沙量分別占宜昌站輸沙總量的50.6%和23.5%，而其他河流和地區(qū)來沙約共占1/4［7］。此外長江來沙又以懸沙為主，占總輸沙量的98%左右［8］，這就為河口三角洲的形成提供了豐富的泥沙來源。

　　2 水沙通量時間系列的特征分析

　　在數(shù)理統(tǒng)計學(xué)中對時間系列進行特征分析時可以用不同的參數(shù)來說明數(shù)據(jù)量的特征，其中用來反映集中趨勢的綜合特征數(shù)有：算術(shù)平均數(shù)、中位數(shù)和眾數(shù)等；描述分配離散情況的特征數(shù)有：方差、標準差和全距等；另外還可以通過偏態(tài)系數(shù)和峰度系數(shù)來描述系列曲線的分布特征，宜昌站和大通站多年平均流量和輸沙量時間系列的相關(guān)表征參數(shù)見表2.

　　從對上述參數(shù)的分析中可知，兩個測站的流量和輸沙量時間系列的特征參數(shù)均存在眾數(shù)＜中位數(shù)＜算術(shù)平均數(shù)的現(xiàn)象，故中下游水沙通量時間系列的分布都呈正偏趨勢，這與我國大多數(shù)河流的水文序列是一致的；由于長江中下游河段的水沙來源不同，造成了宜昌站與大通站的年均流量和輸沙量時間系列之間存在著相異的分布特征，長江中游宜昌站的水沙來量均以上游為主，下游大通站泥沙主要來源于上游，而水量則為上、中、下游的匯集量，故宜昌站的多年平均流量只相當于大通站的48.4%，但輸沙量卻是大通站的1.16倍；宜昌站多年平均流量時間系列曲線的峰度系數(shù)為0.026，大通站為2.602，這表明大通站的年均流量分布曲線較宜昌站的起伏要大；而兩個測站的年均輸沙量時間系列曲線的峰度系數(shù)絕對值相近，即二者的分布曲線形式類似。

　　3 水沙通量時間系列的趨勢分析

　　3.1 趨勢分析

　　采用指數(shù)平滑法對時間系列進行修勻，其目的在于通過平均消除個別年度中夾雜其他因素的影響，從而使數(shù)列呈現(xiàn)出更加明顯的長期變化趨勢，已知經(jīng)過指數(shù)平滑后的流量值在較長的時間范圍內(nèi)，更加趨于穩(wěn)定，圍繞平均值上下波動，除了1954年和1998年兩個特大洪水年的年均流量值有較大偏差以外，其余年份波動性均不明顯，大通站平均變幅不超過4%，宜昌站的平均變幅只有0.6%左右；而兩個測站的多年平均輸沙量經(jīng)過指數(shù)平滑后，則呈現(xiàn)出更加明顯的下降趨勢，這主要是由于近年來長江上游干流及主要支流上一些水壩及電站的修建攔截了大量的泥沙所致（據(jù)統(tǒng)計，至80年代末，長江上游已建水庫11.931座，總庫容約2.05×1011m3［1］），宜昌站在葛洲壩水庫蓄水后多年平均枯季（11～4月）含沙量較蓄水前減少了62%以上［9］?？傮w而言，除了特殊 年份以外，流量和輸沙量在較長的時間系列內(nèi)的變化是具有一定的延續(xù)性的，即連續(xù)幾年大于或小于平均值。

　　3.2 相關(guān)性分析

　　相關(guān)性分析主要是確定時間系列自身內(nèi)部的線性相關(guān)程度，并給出相應(yīng)的定量測度，對于給定的時間系列可以得到以下幾種特征量的無偏或漸近無偏估計：均值τ=1，2，…，K（其中K《N），一般要求N≥50，K值約在N/10左右，在本文中取K=5，將宜昌站和大通站流量及輸沙量時間系列的相關(guān)性分。

　　相關(guān)性分析結(jié)果表明，宜昌站流量時間系列自身內(nèi)部線性相依性隨時間τ的增加而降低，而輸 沙量系列則隨著時間τ的增加而降低。

　　4 水沙通量的階段性變化

　　聚類分析中費希爾（Fisher）最優(yōu)分割方法，即對任意指定的分類數(shù)K，總能將N個樣本分為K類，且使各類直徑的總和S達到最小，即

　　將Ｎ個對象分為Ｋ個類的分類數(shù) 目為R（N，K），則（Ｎ－１，Ｋ）＋Ｒ（Ｎ－１，Ｋ－１），且Ｒ（Ｎ，１）＝Ｒ（Ｎ，Ｎ）＝１。應(yīng)用上述方法對宜昌站和大通站的多年平均流量和輸沙量系列進行分析，得到宜昌站年均流量系列在1968和1975年處S出現(xiàn)最小值，大通站年均流量系 列在1955和1986年處S出現(xiàn)最小值；大通站多年平均輸沙量以1966和1980年為界劃分為三個 時段，而宜昌站則以1961和1968年為界劃分為三個時段。

　　長江流域在較長的時段內(nèi)降水情況沒有顯著的變化，故長江中下游流量沒有明顯的階段性變化趨勢。但是輸沙量的變化是徑流變化與人類活動等綜合因素共同作用的結(jié)果，60年代長江中下游地區(qū)全民大辦農(nóng)業(yè)，開墾荒地，表層植被受到了嚴重的破壞，故這一時段內(nèi)長江中下游的輸沙量普遍高于其他時段；自70年代以來，長江全流域大量修建水庫及沿江護岸工程等水利設(shè)施，攔截了部分泥沙，宜昌站該時段年均輸沙量降幅達到21.39%；80年 代以后，各地積極開展植樹造林、綠化荒山，自1988年開始，國家啟動了長江上游水土流失重點防治工程，簡稱“長治工程”，長治工程使當?shù)刂脖桓采w率明顯提高，水土流失減輕，攔沙蓄水能力有所提高，上游嘉陵江（北碚站）1990～1996年時段的年均輸沙量僅為1954～1989年時段的34.2%，減少了0.96×108t［6］。大通站的多年平均輸沙量也由1967～1980年時段的3.72×108t降為1981～2000年時段的3.20×108t，降幅達13.9%.