近50年珠江三大支流水沙關系變化研究

吳 堯，官明開，曹 煜，龔偉杰

（河海大學港口海岸與近海工程學院，南京210098）

近50年珠江三大支流水沙關系變化研究

吳 堯，官明開，曹 煜，龔偉杰

（河海大學港口海岸與近海工程學院，南京210098）

基于西江高要站、北江石角站和東江博羅站的水沙數據，分析近50年來珠江三大支流水沙關系的變化。利用交叉小波變換分析水沙關系時發現，西江和北江水沙波動能量的穩定性及水沙的正相關性強于東江。水沙冪指數曲線擬合反映了水沙在不同年代的變化關系及相關程度，其中西江水沙擬合斜率和相關系數整體較大。此外，月均水沙的年內分布曲線揭示了三大支流間水沙年內分布的差異。水沙關系的變化不僅與流域風化沉積物的來源有關，而且更是受到水土保持工程和水庫建設等人類活動的影響。

水沙關系；交叉小波變換；水沙冪指數；珠江

由于豐富的自然資源和優越的地理條件，近幾十年來珠江三角洲的經濟飛速發展，已經成為華南對外經貿的門戶。泥沙等陸相物質通過徑流的運輸，在河口區淤積形成珠江三角洲，形成了縱橫交錯的水網系統。河口三角洲地區水沙動力特性復雜，影響到河道的泄洪、通航［1］。大規模水利工程的興建為經濟的發展奠定了堅實的基礎，與此同時，這些人類活動一定程度上也引起了珠江水沙關系發生改變。河流的水沙關系是河流系統中關鍵部分，關系到河流泥沙的輸運、沉積等特性。更進一步地，還影響到水道通航、三角洲地貌和河口生態平衡［2］。水沙關系的變化是一個十分復雜的過程，受到多方因素的共同調控作用。包括河道形態、河流水情、降雨量等自然因素，以及土地利用，水庫建設和人工采砂［3］等人為影響。以往對珠江水沙關系的研究更多的是集中在對水沙通量的變化分析和對水沙過程的統計分析。本文依據珠江三大支流上主要控制水文站長時間尺度的實測水沙資料，利用交叉小波變換，水沙冪指數關系和繩套關系等角度分析近50年各支流水沙關系的變化。對珠江流域長時間尺度水沙關系的研究，不僅對河道整治、河口生態治理等有重要意義，也能為各流域水資源的統籌開發和優化配置提供依據。

1區域概況及數據

珠江是中國流量第二大的河流，其水系主要由西江、北江和東江三大支流組成。徑流流經江西、廣西、湖南和廣東等省份。流域總面積約為45.4×104km2，上游的徑流通過復雜的珠江三角洲河網區匯入南海。珠江多年平均的水沙通量分別大致為2.85×1011m3/a和0.75×108kg/a。其中西江是珠江最大的支流，根據高要站的實測資料，其多年平均的水沙通量占珠江水沙總量的77％和89％。而北江和東江的水沙通量相對較小。流域受到熱帶和亞熱帶季風氣候的影響，大部分的降雨主要集中在夏季，河流水情具有明顯的季節性特征。

高要站、石角站和博羅站分別為西江、北江和東江干流的主要水文控制站，本文研究資料為西江高要站，北江石角站和東江博羅站3個站點的月平均徑流量和含沙量，研究的時間序列跨度為1961～2012年。為了分析長時間尺度下水沙關系的變化，將研究數據的時間跨度依據年代不同劃分成5個部分。資料來源于水利部發布的水文年鑒，文中對泥沙的研究均指懸移質泥沙。

2研究方法介紹

2.1交叉小波分析

交叉小波變換能反映兩組離散時間序列的相關程度，而且能反映兩者在各時頻域的位相特征和細部特征，被廣泛的應用于水文、氣象等資料間的相關分析［4］。其中，交叉小波功率譜可以定義為

式中：S為平滑算子。

2.2水沙冪函數曲線

水沙的冪函數關系為

公式兩邊經過對數變換后轉化為線性關系

式中：Cs為含沙量，Q為流量；a和b都為擬合的系數。a越大意味著易受徑流侵蝕和輸運的風化沉積物來源越充足，b越大意味著徑流的挾沙能力和侵蝕能力越強［5］。

3結果分析和討論

3.1水沙通量的交叉小波分析

交叉小波分析能夠揭示徑流量和含沙量在不同時頻域的相關關系及位相特征。因此對西江高要站、北江石角站和東江博羅站的月均水沙通量運用交叉小波分析，來研究各支流水沙關系的變化。圖1中粗實線區域表示相關性通過了紅噪聲標準譜的顯著性檢驗。細實線圍成的“U”型區域內表示交叉小波譜的有效范圍，避免了由于數據長度有限性導致的邊界效應。矢量箭頭表示水沙的位相關系，箭頭向右為正相關，箭頭向左為負相關。

從圖1中的交叉小波功率譜可以明顯的看出，3條支流的水沙波動能量都主要集中在1年左右的主周期上。而在其他頻率尺度上只檢測到分布不規律且持續時間極短的顯著性特征。因此忽略高頻率尺度的波動，主要分析3條支流在1年主周期上的水沙相關性及其在近50年的變化規律。其中，高要站的交叉小波功率譜在1年主周期上的波動能量在研究周期內均較穩定，只在2005年后才出現波動能量的消失，表明此時的水沙相關性發生改變（圖1（a））。石角站的波動能量同樣表現穩定，只有在2002～2003年出現了波動能量的缺失（圖1（c））。博羅站的水沙相關性的變化相對較明顯。在1990年后，東江水沙關系1年主周期的波動能量出現明顯的波動和變化（圖1（e））。交叉小波凝聚譜比交叉小波功率譜揭示了范圍更大的顯著性時頻域，其主周期同樣集中在8～16個月。三大支流的水沙交叉小波凝聚譜在周期小于1年的頻段內均檢測到不連續且分布不均勻的波動能量。而在更長周期的頻段上，石角站的交叉小波凝聚譜檢測的水沙關系的顯著性覆蓋范圍最廣，除了20世紀70年代和20世紀80年代在16個月到64個月的部分區域，其余時頻段均檢測到顯著的波動能量（圖1（d））。高要站在20世紀90年代后，在16個月到64個月的頻段內存在顯著的水沙相關性，覆蓋范圍相對較小（圖1（b））。而博羅站從20世紀90年代后，在所有的頻域范圍內波動能量均較顯著（圖1（f））。從位相關系來看，3個測站的水沙關系都呈現顯著的正相關。高要站和石角站水沙位相關系的正相關性要強于博羅站。

圖1 高要站（a/b）、石角站（c/d）和博羅站（e/f）水沙的交叉小波功率譜（第一列）和交叉小波凝聚譜（第二列）Fig.1 The cross wavelet transform（first line）and wavelet coherence（second line）between monthly runoff and sediment concentration of Gaoyao（a/b），Shijiao（c/d）and Boluo（e/f）station

3.2水沙的冪函數關系分析

冪函數關系曲線被廣泛的應用于水沙關系的研究。從圖2可以直觀地看出珠江3個支流的水沙冪函數關系在各個年代的變化。其中高要站在各個年代水沙關系的冪函數擬合直線的斜率（b）最大，明顯大于石角站和博羅站擬合直線的斜率。這可能是由于西江徑流量較大使得徑流挾沙能力較強。而以石角站和博羅站為代表的北江和東江徑流則具有較相近的挾沙能力。另外，3個測站在20世紀70年代和80年代擬合直線的斜率均明顯大于各自在近20 a擬合的斜率，相應的直線的截距（ln（a））均較小，說明珠江三大支流挾沙能力在近20 a都明顯減弱。風化沉積物的來源與珠江流域的水土保持措施密切相關。研究發現，20世紀50年代到80年代，由于土地改革等農業、經濟改革，流域內出現了大規模的森林采伐，水土流失加劇。尤其從1975～1984年，珠江流域的森林覆蓋率從38％猛降到27％［6］。而從1990年后，隨著人們在資源開發過程中環保意識的增強，珠江流域的《水土保持法》得到貫徹實施。大量水土保持工程的建設改善了流域中水土流失嚴重的現象，流域的水土治理面積逐漸增大，水土流失匯入徑流的泥沙來源顯著減少［7］。

冪函數擬合的參數ln（a）和b以及擬合的相關系數R2見表1。3個測站的水沙關系在不同的年代有不同的相關關系。其中，高要站的水沙的相關系數R2最大（R2∈(0.82,0.90)），說明高要水沙關系相關性最好。其最大的相關系數出現在20世紀70年代和80年代，在2000年后，水沙相關性減弱。博羅站的R2整體較小（R2∈(0.53,0.73)），表明東江的水沙相關性在3條支流中最弱，這與前文位相關系分析結果相吻合。而且在近20年東江的水沙相關性還有進一步減弱的趨勢。石角站的水沙相關系數R2適中，最小的水沙相關系數出現在20世紀90年代（R2=0.57）。

圖2 高要站、石角站和博羅站在各個時期的水沙冪函數擬合Fig.2 The power function of water and sediment of Gaoyao，Shijiao and Boluo station at different periods

表1 高要、石角和博羅站在各個時期的水沙冪函數擬合參數（ln（a）和b）和相關系數R2Tab.1 Fitting parameters（ln（a）and b）and coefficient of determination（R2）of Gaoyao，Shijiao and Boluo station for different decades

圖3 高要站、石角站和博羅站在各個時期的月平均水沙年內分布曲線Fig.3 Distribution curves within year of monthly average discharge and sediment concentration for different decades at the Gaoyao，Shijiao and Boluo station

3.3月均水沙關系年內分布曲線

正常情況下，月平均徑流量從枯季到洪季有一個明顯的增加過程，再從洪季到枯季有一個減小過程，可以分別稱之為漲水期和落水期。與徑流密切相關的月平均含沙量隨之變化，但它與徑流的變化過程存在相位差，使得水沙關系在一年的12個月份里形成一個具有上升過程和下降過程的繩套曲線，即水沙關系年內分布曲線。通常在漲水期階段，年內首次降雨形成的洪水使得地表泥沙易受侵蝕，流域產沙量大。當洪峰落后于沙峰時，則呈現順時針的繩套曲線，當沙峰落后于洪峰，則為逆時針的繩套曲線。同樣地，將研究周期依據年代分成5個研究階段，利用月均水沙年內分布曲線來探究近50年水沙關系在不同年代的年內分布變化。

從圖3可以看出，除了20世紀80年代，在其余4個研究階段西江高要站的水沙年內分布曲線的漲水期和落水期都接近重合。20世紀80年代的水沙曲線呈現明顯的順時針環形，漲水期的含沙量大于落水期的含沙量，說明年內首次洪水沖刷侵蝕大量的風化沉積物后匯入徑流，沙峰領先洪峰。且相比其余4個階段，20世紀80年代的洪季平均含沙量相對較大，最大的月均含沙量（7月）達到0.7 kg/m3左右。前面分析已經提到過，20世紀80年代珠江流域的森林采伐現象最為嚴重。而西江作為珠江最大的支流，其上游流經云南、貴州和廣西等省份，這部分流域屬于巖溶地貌，土壤貧瘠且脆弱。植被的破壞將造成難以修復的土地荒漠化現象［8］。在降雨形成地表洪水侵蝕作用下，風化物被水流帶走使得徑流中含沙量顯著增大。北江石角站的水沙年內分布曲線在前4個研究階段均呈現順時針環形且差異不大，表明其水沙關系的年內分布特征較穩定。而對于東江博羅站，前4個階段的曲線也都呈現順時針環形。與北江不同的是，東江前4個階段的水沙年內分布繩套曲線不斷收縮，沙峰和洪峰的峰值均不斷減小，尤其是含沙量減小程度更大。泥沙的減少與水庫攔沙密切相關。東江在近50內修建的許多大中型水庫，尤其是新豐江（1962）、楓樹壩（1973）和白盆珠（1985）3座大型水庫，其庫容量占廣東省大型水庫庫容總量的63.4％［9］。水庫建設已經被證明是入海泥沙通量減少的主要原因。值得注意的是，三大支流在2000年后繩套曲線在漲水期和落水期趨于重合，表明洪峰和沙峰的相位趨于一致。這種現象可以歸結為2000年后大量水土保持工程的建設改善了流域水土流失的狀況，受到首次降雨沖刷侵蝕的泥沙來源顯著減少。

4結論

珠江三大支流在近50年的水沙關系都發生了不同程度的變化。對水沙位相關系分析和相關系數R2的計算結果都發現，3條支流中東江的水沙相關性最弱。尤其是在1990年后，東江水沙交叉小波譜在1年主周期的波動能量變化明顯。此外，從20世紀60年代起東江的沙峰和洪峰峰值不斷減小，反映在繩套曲線上是前4個年代的曲線形態不斷收縮，尤其是含沙量減小顯著。這可以歸結為東江上大型水庫對徑流的調控和攔沙作用。北江的水沙關系在研究周期內相對較穩定，除了20世紀90年代的水沙相關系數較小（R2= 0.57）。西江的水沙相關系數整體較大，說明其水沙相關性最顯著。此外，冪指數擬合結果發現西江水沙的擬合直線斜率整體均較大，揭示了西江徑流較強的侵蝕能力和挾沙能力。

珠江三大支流的水沙關系變化不僅各有特點，而且還存在一些共性。整體來看，三大支流在近20年擬合直線的斜率均明顯小于各自在20世紀70年代和80年代的斜率，表明在近20年珠江流域整體的風化沉積物等泥沙來源減少。這與流域的水土保持工程的實施有關，近20年流域的森林覆蓋率增加且水土流失減少。此外，水沙關系的年內分布曲線結果還發現，2000年后三大支流的繩套曲線形態在漲水期和落水期均趨于重合，意味著洪峰和沙峰的相位趨于一致，受到年內首次降雨沖刷侵蝕的泥沙來源顯著減少。

［1］張華慶，呂忠華.珠江河口水沙數值模擬系統［J］.水道港口，2002，23（2）：51-53. ZHANG Q H，LV Z H.A Numerical Model of Current and Sediment in Pearl River Estuary［J］.Journal of Waterway and Harbor，2002，23（2）：51-53.

［2］趙玉，穆興民，何毅，等.1950-2011年黃河干流水沙關系變化研究［J］.泥沙研究，2014（4）：32-38. ZHAO Y，MU X M，HE Y，et al.Relationship between runoff and sediment discharge in the main channel of Yellow River from 1950 to 2011［J］.Journal of Sediment Research，2014（4）：32-38.

［3］張蔚，嚴以新，諸裕良，等.人工采沙及航道整治對珠江三角洲水流動力條件的影響［J］.水利學報，2008，39（9）：1 098-1 104. ZHANG W，YAN Y X，ZHU Y L，et al.Impact of sand excavation and waterway regulation on hydrodynamics of Pearl River net?works［J］.Journal of Hydraulic Engineering，2008，39（9）：1 098-1 104.

［4］Torrence C，Compo G P.A practical guide to wavelet analysis［J］.Bulletin of the American Meteorological Society，1998，79（1）：61-78.

［5］Zhang W，Wei X Y，Zheng J H，et al.Estimating suspended sediment loads in the Pearl River Delta region using sediment rating curves［J］.Continental Shelf Research，2012，38：35-46.

［6］夏漢平.論長江與珠江流域的水災、水土流失及植被生態恢復工程［J］.熱帶地理，1999，19（2）：124-129. XIA H P.Flood disasters，soil erosion，and eco?restoration of vegetation in the Yangtze and Pearl River valleys［J］.Tropical Geogra?phy，1999，19（2）：124-129.

［7］吳創收，楊世倫，黃世昌，等.1954-2011年間珠江入海水沙通量變化的多尺度分析［J］.地理學報，2014，69（3）：422-432. WU C S，YANG S L，HUANG S C，et al.Multi?scale variability of water discharge and sediment load in the Pearl River during 1954-2011［J］.Acta Geographica Sinica，2014，69（3）：422-432.

［8］賴天锃，張強，陳永勤.1960～2010年西江流域水沙變化特征及其成因［J］.武漢大學學報：理學版，2015（3）：271-278. LAI T Z，ZHANG Q，CHEN Y Q.Sediment and runoff changes in the West River basin：changing properties，causes and implica?tions［J］.Journal of Wuhan University：Natural Science Edition，2015（3）：271-278.

［9］Du H Q，Wang Y H，Gao L H，et al.Effects of reservoirs on hydrological regimes in several reaches of East River［J］.Engineering Journal of Wuhan University，2011，4：012.

Detecting the changing runoff?sediment relationship among three main tributaries of the Pearl River in the last 50 years

WU Yao,GUAN Ming?kai,CAO Yu,GONG Wei?jie
(College of Harbor,Coastal and Offshore Engineering,Hohai University,Nanjing 210098,China)

Based on the runoff and sediment dataset from Gaoyao station of Xijiang,Shijiao station of Beijiang and Boluo station of Dongjiang during the last 50 years,the changing runoff?sediment relationship among three main tributaries of the Pearl River was detected.The cross wavelet analysis reveals that the stability of runoff?sedi?ment fluctuant energy and the positive correlation property of Xijiang and Beijiang are more significant than Dongji?ang.Additionally,the power function of water and sediment was applied to study the changing runoff?sediment rela?tionship and their correlation characteristic during different decades.It finds that the slope of runoff?sediment power function is steeper and the correlation coefficient of determination is larger in Xijiang.Furthermore,the distribution curves within year of monthly average runoff and sediment illustrate different characteristics for three tributaries. The changing runoff?sediment relationship is affected by the source of weather sediment and the human impact via soil-water conservation projects and the construction of reservoirs.

runoff?sediment relationship;cross wavelet transform;power function of water and sediment;the Pearl River

P 333；TV 11

A

1005-8443（2016）04-0405-06

2016-01-19；

2016-03-30

吳堯(1992-)，男，福建省漳州人，碩士研究生，主要從事河口水動力學研究。Biography：WU Yao（1992-）,male,master student.