長江口青草沙水庫前沿河床演變與失穩風險研究

盛 皓，戴志軍，梅雪菲，葛振鵬，黎樹式，高近娟

(華東師范大學 河口海岸學國家重點實驗室，上海 200062)

長江口青草沙水庫前沿河床演變與失穩風險研究

盛 皓，戴志軍，梅雪菲，葛振鵬，黎樹式，高近娟

(華東師范大學 河口海岸學國家重點實驗室，上海 200062)

位于長江口南北港分流口的青草沙水庫是上海市最大的城市供水水源地，開展該水庫前沿河床穩定性和失穩風險研究具有重要的理論意義和工程指導價值?；诖?，利用長江口南北港分流口2000年至2013年期間實測的高精度地形資料，分析青草沙水庫前沿河床近期地貌變化特征，為水庫庫堤安全預警提供相關理論支撐。結果表明：青草沙水庫前沿河床沖淤變化在建庫前后有明顯差異，建庫前總體表現為河槽沖刷、沙洲淤積，呈現準沖-淤振蕩的性質；建庫后變為幅度逐漸減弱的持續沖刷；沿河床河槽形態由U型向U型與V型河槽疊加的復式河槽轉變。青草沙水庫前沿沙體沿落潮主流方向下移，成形沙體呈先增長、再減少，最后趨于平衡的態勢，隱形沙體大體呈減少趨勢。青草沙水庫庫堤前沿近600 m位置是河床失穩的重點風險區域，尤其是水庫庫堤中上部河床處于不穩定狀態，如前沿邊坡所在河槽進一步逼近水庫前沿，則河床面臨進一步沖刷的可能。

河床演變；青草沙水庫；長江口；水庫前沿；分流口；河槽沖刷

Abstract: The Qingcaosha Reservoir (QCSR), located in the South and North Channel bifurcation, is the largest urban water supply source in Shanghai. It is of great significance for theory and practices to conduct the stability assessment of the frontal river bed in the QCSR. Thereafter, based on high resolution topography and sediment data measured during 2000-2013, the geomorphology characteristics of the frontal river bed of QCSR and associated stability of the river bed were analyzed to provide theoretical support for safety pre-warning of QCSR. The results showed that: 1) The variations in erosion and accretion of the river bed present clearly differences before and after QCSR operation. Before the reservoir operation, the river bed suffered channel erosion and shoal deposition with oscillations, which depended on the diversion of flow and sediment of the South and North Channel bifurcation. After the reservoir operation, the changes in the frontal river bed had been shifted from the fluctuations of erosion/accretion to the uninterrupted erosion with weak scouring trends. Moreover, the corresponding channel configuration appeared changes from the U-shaped to the combined U-shaped with V-shaped formation. 2) The frontal sand bodies along QCSR moved downward to follow the mainly ebb-tidal current direction. During the period of movement, the developmental sand bodies firstly expanded, and then shrinked, to the end with tendency to reach equilibrium, while the clocked sand body occurred reduced trends. 3) It is the key risked region that the location of about 600 m distance to the frontal area of the reservoir is seriously unstable, especially to the upper part of the river bed along QCSR. It can be expected that there is further scouring possibility for the river bed if the channel located at the front river bed slope will be moved approximately to the dyke of QCSR.

Keywords: river bed evolution; Qingcaosha Reservoir; Changjiang Estuary；the frontal river bed of the reservoir; bifurcation mouth; channel scouring

河口地貌的發育和演變是一個極其復雜的過程，它不僅受控于上游流域來水來沙條件，同時也受到海域以及局部人類活動作用的影響[1]。河口上游以及分流口來水來沙的變化，將會引起下游河道河床形態的變化，而河床形態的變化直接影響河道水流結構和水流阻力，進而控制進入下游河道的能量分配和泥沙輸移量，由此成為控制河口三角洲沖淤進退的重要因素之一[2]。水流和泥沙二者相互配合，共同塑造河床的邊界形態，河道的邊界狀態反過來又進一步控制河流的水流和泥沙過程，進而影響局部水利工程的環境?；诖?，不少研究指出流域和河口水利工程的建設往往會不同程度地改變工程區域的邊界形態，進而對河床形態產生影響[3-6]。因此河口以及流域水利工程對河床或河床邊界地貌形態的影響，成為當前分析和研究河口局部工程河勢穩定與否的關鍵。作為控制河道水沙運動的邊界，河床形態變化成為河口地貌研究的重要組成部分。

圖1 研究區域和斷面位置Fig. 1 Study area and section location

1 資料來源與研究方法

1.1資料來源與處理

數據及其處理主要分為以下兩部分: 1) 利用雙頻測深儀和GPS設備在長江口南北港分流口區域采集2000年、2002年、2004年、2007年以及2009～2013年共9年實測水深資料(資料為上海地質調查研究院和河口海岸學國家重點實驗室每年4～5月共同所測)?；谝陨蠈崪y水深數據，用ArcGIS軟件通過Kriging插值方法生成規則矩形網格數字高程模型(DEM)，進而計算出該區域不同年份的沖淤值。沿南北港分流口NE方向提取七個河槽斷面，以分析青草沙水庫前沿河床的變化。所有斷面南側起點記為固定樁位A點，北側終點記為B點樁位，見圖1(b)。2)基于DEM進行柵格計算，提取了研究區域相關年份的河槽容積、-2 m等深線以及深泓線，以分析一系列大型水利工程前后水庫前沿河床演變規律。

1.2研究方法

1.2.1 成形沙體與隱形沙包判別法

諸多研究表明，成形河槽的活動沙體以及隱形沙包淤積體的輸移是影響河勢變化和沙島成陸的重要因素[15-18]。根據前人對活動以及隱形沙體的定義，本文以-5 m等深線閉合圈周的暗沙作為成形沙體識別標志，以不同時間尺度沖淤圖上淤積厚度大于3 m的環形沙包作為隱形沙包識別標志[15-17]。基于此，通過ArcGIS平臺對研究區域的成形沙體和隱形沙包進行鑒別，進而利用淺水區成形沙體自上而下的遷移、沙洲洲頭后退的速度以及洲尾向下延伸的速度作為成形沙體廓線的移動速度；利用深水區沿程沖淤形成的隱形沙包的規模變化以及移動規律對隱形沙體廓線移動速率進行估算[16]。

1.2.2Mann-Kendall法

Mann-Kendall法是一種檢驗時間序列發生突變的方法，其也是在水文、氣象領域應用較廣的非參數檢驗方法，用于檢驗時間序列的變化趨勢，相關計算如下[19-20]。

對于有n個樣本的序列x，假設序列隨機獨立且具有相同的概率分布，構造一個秩序列：

式中：UF1=0，E(sk)為均值，var(sk)為方差。

UFi為標準正態分布，它是按x順序x1，x2，...，xn計算出的統計量序列，給定顯著性水平α，查正態分布表臨界值u0.05=±1.96，若|UFi|>Uα，則表明序列存在明顯的趨勢變化。按序列x逆序xn，xn-1，...，x1，再重復以上過程，同時使UBk=-UFk(k=n，n-1，…，1)，UB1=0。將UFk和UBk繪在圖上，若UFk和UBk兩條曲線在置信區間存在交點，那么交點處對應的就是突變點。

本文采用Mann-Kendall法對青草沙水庫前沿河床斷面進行坡度突變檢驗以尋找河床邊坡坡度突變點，為進一步探討水庫前沿河床的穩定性提供依據。

2 結果分析

2.1青草沙水庫前沿河床地貌沖淤變化

2007～2013年是長江口南支河段大型水利工程實施的集中期，其中在南北港分流口附近實施的水利工程主要有青草沙水庫工程(2007～2011年)、新瀏河沙護灘工程和南沙頭通道潛堤工程(2007～2009年)[21-22]。選擇的地形資料可反映工程實施前后水庫前沿河床的變化，進而將區域地貌變化分為如下三個階段分析：1)2000～2007年期間；2)2007～2011年期間；3)2011～2013年期間。由此以不同時間尺度下的河床沖淤和河槽斷面形態變化來反映水庫建設前后河床地貌變化，見圖2、圖3。

1) 2000～2007年期間

2000～2002年期間，新橋通道、新橋水道以及北港主河槽均出現河槽泥沙沖刷，總體沖刷幅度大致為2～4 m，局部沖刷高達4 m以上，見圖2(b)中a。然而，相應的分流口沙洲如新橋沙、扁擔沙則處于淤積狀態，新橋沙平均淤積幅度由沙洲軸心2～4 m向外延降低到0～2 m的幅度，河床總體年平均沖刷幅度為0.5 m。同時，水庫庫堤前沿中部出現小幅淤積的態勢，淤積幅度為0.27 m，下部呈現小幅沖刷態勢，幅度為0～2 m，見圖2(c)中a。河槽形態變化主要表現上部和下部近岸淺V型河槽刷深、邊坡陡化，中部單一U型河槽向雙U型河槽演變，見圖3。相應的河槽容積由25.81億m3增大至30.28億m3，較2000年增加17.32%，見圖4。

2002～2004年期間，新橋通道仍以沖刷為主，沖刷幅度在3 m左右，而新橋水道和北港河槽由之前的沖刷態勢變為大幅淤積，總體平均淤積厚度為3 m左右。扁擔沙沙尾北翼由淤積態勢變為大幅沖刷，局部高達4 m左右，前沿河床總體呈現淤積態勢，年平均淤積厚度為0.75 m，見圖2(b)中b。水庫庫堤前沿上部呈現1 m左右小幅沖刷態勢，中下部呈2 m左右大幅淤積態勢，見圖2(c)中b。相應的河槽形態變化表現為由U型與V型疊加的復式河槽向單一型U型河槽轉變，呈展寬拉伸態勢,見圖3。由于大幅淤積，河槽容積大幅減小，由先前的30.28億m3降低至24.13億m3，下降幅度高達20.31%，見圖4。

2004～2007年期間，水庫前沿河床仍以沖刷為主，但是較2000～2002年沖刷略有減弱，年平均沖刷幅度為0.45 m，其中中部河床出現4 m以上大幅沖刷，見圖2(b)中c。水庫庫堤前沿河床呈現出與2000～2002年相近的沖淤態勢，中部呈現1.5 m左右的淤積態勢，上部和下部分別呈現2.5 m和1 m的沖刷態勢，見圖2(c)中c。河槽形態發生劇烈變化，呈刷深束窄態勢，由單一型U型河槽向U型與V型疊加的復式河槽轉變，如圖3所示。河槽容積增幅度高達25.16%，恢復至30.2 億m3，如圖4所示?？傮w來說，2000～2007年青草沙水庫周圍河床總體表現為河槽沖刷、沙洲淤積、水庫庫堤前沿河床沖淤變化呈現準沖-淤的年際振蕩規律。

圖3 青草沙水庫前沿河槽斷面形態變化Fig. 3 Changes of the river channel in front of the Qingcaosha Reservoir

圖4 河槽容積和平均水深變化Fig. 4 Volume and average depth variations of the channel in front of the Qingcaosha Reservoir

2) 2007～2011年期間

2007～2009年期間，青草沙河段的河槽和沙洲都出現不同程度的淤積，扁擔沙沙尾出現3～4 m左右強淤積，見圖2(b)中d。水庫庫堤前沿河床表現為上部強沖刷，中下部弱沖刷，年平均沖刷幅度為0.63 m左右，中部前沿河床由幅度2 m左右的強淤積轉為0.3 m的弱沖刷態勢，見圖2(c)中d。期間河槽形態大體不變，水庫庫堤前沿邊灘刷深，向庫堤逼近，邊坡形態較之前變陡，見圖3。區域河槽容積略有增加，變為31.17億m3，見圖4。

2009～2010年期間，新橋通道、新橋水道以及分流口沙洲扁擔沙和新橋沙等出現不同程度的淤積，見圖2(b)中e。水庫庫堤中部河床呈現幅度為2.7 m的沖刷態勢，水庫庫堤前沿上部河床的侵蝕依舊存在，見圖2(c)中e。水庫尾閭的前緣則出現淤積，可能是水庫中部河床出現沖刷而侵蝕下的泥沙堆積于此所致。除青草沙水庫上部河槽略展寬外，整體河槽形態與2009年基本保持一致，如圖3所示。在此期間，河槽容積小幅增加至32.26億m3，見圖4。

2010～2011年期間，扁擔沙軸心出現淤積態勢，新新橋通道則出現局部沖刷接近4 m的情勢。此外，新橋沙和新橋通道都出現淤積接近2 m的狀態，水庫前沿河床中部雖有局部沖刷，但沖刷幅度比先前2009～2010年減弱，見圖2(b)中f。水庫庫堤前沿河床中部由之前的大幅沖刷變為1 m左右的淤積態勢，而整體由2 m左右沖刷態勢變為0.6 m左右的淤積態勢，見圖2(c)中f。上部河槽近水庫邊緣刷深，由呈淺U型河槽變為深V型河槽；中下部河槽形態變化較小，略有刷深，如圖3所示。相應的河槽容積因此出現小幅下降，變為31.51億m3，如圖4所示?？傮w而言，2010～2011年青草沙水庫前沿河床由沖刷向淤積轉變。

與2000～2007年相比，2007～2011年期間青草沙水庫前沿河床沖刷減弱，由沖刷轉為弱淤積。具體表現為河槽整體沖刷，局部淤積，沙洲弱淤且呈縮窄拉長態勢。庫堤前沿河床表現為不同程度的沖刷，年平均沖刷幅度為0.6 m，中部河床變化最大，由之前的持續淤積變為持續沖刷，沖刷幅度在1.1 m左右。河槽容積基本穩定，出現小幅波動上升，如圖4所示。

3)2011～2013年期間

2011～2012年期間，青草沙水庫前沿河床整體呈現出弱淤積態勢，平均淤積厚度為0.1 m，見圖2(b)中g，而水庫庫堤前沿上部和中部河床表現為較大幅度的淤積，平均淤積厚度為0.9 m，局部高達4 m，下部呈弱沖態勢，見圖2(c)中g。然而，在2012～2013年期間，水庫前沿河床轉變為以沖刷為主，平均沖刷幅度在0.5 m左右，見圖2(b)中h，水庫庫堤前沿河床表現為中上部強沖刷，下部弱沖刷，平均沖刷厚度為1.2 m，其中水庫頭緣局部沖刷超過2 m，見圖2(c)中h。2011～2013年期間，青草沙水庫前沿河床表現為先淤后沖，年平均沖刷0.2 m。河槽形態整體較為穩定，上部河槽由于淤積水深變淺，近水庫邊緣深V型河槽衰退，如圖3所示。由于河床由弱淤轉為弱沖刷，相應的河槽容積從先前的31.51億m3減至30.51億m3，然后又增至32.79億m3，見圖4。

整體而言，2000～2013年青草沙水庫前沿河床沖淤變化在大型水利工程建設前后有明顯差異：建庫前總體表現為河槽沖刷、沙洲淤積，呈現準沖-淤振蕩性質；建庫后表現為持續沖刷。水庫庫堤前沿河床也由之前的周期性沖-淤變化轉為幅度逐漸減弱的持續沖刷。其中庫堤上部河床一直處于沖刷狀態，沖刷強度先增強再減弱；庫堤中部由淤積變為小幅持續沖刷；庫堤下部除2002～2004年外均呈現沖刷態勢，建庫后有減弱趨勢，其沖淤振蕩幅度較小，相對穩定。由此可見，庫堤上部和中部河床受工程影響較大，應重點監測。

2.2青草沙水庫前沿沙體輸移規律

河勢演變對于底沙輸移有重要影響，為直觀地反映北港主槽河勢演變，在此進一步提取了2000～2013年期間相關年份主槽深泓線和-2 m等深線數據，見圖5、圖6。結果發現，水庫前沿區域深泓線北偏，下段大幅南移，并且深槽逐步趨于穩定，如圖5所示。2000～2010年間，-2 m等深線在扁擔沙沙尾、新橋沙以及堡鎮沙附近發生劇烈發變化，沙體先后歷經切割重組、再切割的過程。2010～2013年期間，除新橋沙附近線圈略有下移，水庫庫堤中上部有切割重組跡象外，-2 m等深線圈基本穩定，見圖6。由此可見，中央沙圈圍和青草沙水庫工程以及同時期施工的南北港分流口工程對穩定南北港分流口河勢和增強北港河勢的穩定性起到了積極的作用。

圖5 北港深泓線變化Fig. 5 Changes of the thalweg of the North Channel

圖6 青草沙水庫前沿-2 m等深線變化Fig. 6 Changes of the -2 m isobaths in front of the Qingcaosha Reservoir

北港底沙輸移呈現成形沙體和隱形沙體兩種形式。在此，以-5 m等深線圈圍暗沙作為成形沙體識別標志，以連續沖淤圖上淤積厚度大于3 m的環形沙包作為隱形沙體識別標志，繪制出2000年至2013年期間水庫前沿沙體的分布圖，見圖7、圖8。成形沙體主要有扁擔沙，新橋沙，堡鎮沙，河槽幾乎無沙體分布。從成形沙體的沙體廓線移動來看，自2000年到2013年總體表現為沙體整體拉伸，從廓線移動路徑可以計算出：沙頭下移約2 500 m，沙尾下移約4 200 m，輸移路線與落潮主流方向一致，平均輸移速率分別為192與323 m/a，見圖7。

圖8 青草沙水庫前沿隱形沙體分布Fig. 8 Distribution of the clocked sand bodies in front of the Qingcaosha Reservoir

2001年扁擔沙尾南翼-10 m線切穿，形成新新橋通道[14]，新新橋通道下部沙體出現淤積擴張，到2007年則南延300 m左右。在此期間，新新橋通道先淤后沖，其中部在2004年出現大量淤積，到2007年沙體重新被切穿，新新橋通道較之前略微北移；2007年到2010年，沙體與扁擔沙沙尾相連，新新橋通道受阻；2010年以后，沙體出現明顯分化，扁擔沙沙尾南側再次被切穿；到2013年較之前顯著下移，且出現小幅萎縮，見圖7。

2000年到2007年青草沙水庫庫堤前沿上部沙體沿庫堤下移，在中下部快速堆積。2007年至2010年，沙體較2007年略有萎縮，與沖淤圖結合可以發現，總體表現為上部弱沖刷，中下部先淤后沖。2010年以后，沙體整體較穩定，見圖7。

結合-2 m等深線和成形沙體變化來看，2000～2004年期間，堡鎮沙區域-2 m等深線圈收縮，由條帶狀變為波狀分布，堡鎮沙呈淤長態勢。2004～2007年期間，受北港中段深泓線北偏影響，-2 m等深線圈向北擴張，堡鎮沙略有減少。2007～2011年期間，受青草沙工程影響，中段深泓線繼續北偏，堡鎮沙受到沖刷，堡鎮頭部條帶狀等深線圈被切割成若干個點狀線圈，堡鎮沙逐漸萎縮。2011年以后，青草沙工程完工，北港主流逐漸穩定，堡鎮沙趨于穩定，見圖6、圖7。

總體而言，水庫前沿河床成形沙體先增長、再減少，最后趨于平衡。沙體整體沿落潮主流方向下移，從而導致青草沙水庫前沿河床呈河槽沖刷，沙洲弱淤積態勢，而庫堤前沿河床出現上部先沖再淤，中下部由強到弱的淤積過程。

此外，隱形沙體主要分布在南北港分流口以及水庫庫堤前沿，堡鎮沙附近有少量分布，說明水庫庫堤前沿邊灘以及沙尾部位是重點淤積區域，北港主槽淤積不明顯，見圖8。2000年至2013年，隱形沙體整體呈減少趨勢。隱形沙體沿落潮主流方向向下輸移7 800 m，移動速率為1 560 m/a。

值得提及的是自三峽水庫蓄水之后，長江流域徑流量相對比較穩定，仍呈現洪季大、枯季小的特征[23]，而入海的泥沙通量卻急遽減少，以大通站為例，1953～2002年平均年輸沙量為4.2億t，而2003～2013年平均年輸沙量僅有1.4億t，下降幅度約高70%，使得部分泥沙被沖刷向下游輸送[2, 24]，北港2000～2004年期間隱形沙體大幅增加可能與此有關。爾后隨著南北港分流口一系列護灘工程實施后，分流口和北港的河勢穩定性增強，進入河口水體的含沙量有所下降，導致隱形沙體呈現大幅減少的趨勢。

3 青草沙水庫河床穩定狀態分析

利用Matlab編程計算歷年實測斷面不同等深線，并運用傳統的剖面沖淤對比統計方法分析不同年份在同一位置的沖淤情況，探討不同位置河床的穩定狀態。

斷面1所處位置為南北港分流口，除2010年斷面水深明顯變淺外，河床斷面形態基本維持不變，而僅在相同的位置出現泥沙淤積或泥沙侵蝕，見圖3(a)。

斷面2恰位于青草沙水庫的前沿位置，從2000～2013年，青草沙水庫前沿河床地形變化較大。斷面2中距青草沙水庫最近的(a點固定樁位)位置由2007年的V型斷面(V型尖端距1 200 m位置)，以每年150 m的速度向青草沙水庫庫堤不斷逼退，其中V型尖端到2013年已離固定樁位a點不到300 m，同時V型河槽近岸坡度逐漸變陡，在2013年邊坡已經完全沖刷，可能將對水庫的庫堤泥沙沖刷產生較大影響，見圖3(b)。

斷面3經歷2007年的大沖刷后，自2009年到2013年，斷面整體維持形態不變，與斷面2類似的是，斷面3距離青草沙水庫最近的河槽由U型不斷向V型變形，且邊坡角度越來越陡，特別是2010年河槽V型尖端距離水庫前沿樁點a不到100 m，其近岸邊坡接近于直角，盡管在2011～2012年河槽出現少量淤積，但邊坡的形態不變。如斷面3所處位置進一步發展，將可能對水庫的河床失穩產生嚴重影響，見圖3(c)。

斷面4基本可表征水庫中部庫堤所在河床的位置。與斷面2、3比較，斷面4的形態較為簡單，中部發育U型河槽，斷面形態較為平緩，平均水深小于10 m，但在近岸1 000 m以內，河槽形態變化較大，呈現不穩定狀態，見圖3(d)。

斷面5和斷面6基本反映了寬而深的U型槽上疊加有典型的V型復式河槽，見圖3(e)、3(f)。斷面5距離庫堤較近的V型河槽盡管邊坡沒有斷面2陡，但邊坡處在不斷沖刷的狀態，V槽有向庫堤靠攏的趨勢，不容忽視。另外從斷面6的剖面變化來看，V型河槽尖端亦在近岸300 m以淺位置，雖然V型河槽最深處位置基本穩定，但近岸邊坡相對較陡，同時斷面中部U型河槽有向庫堤遷移的趨勢，V型河槽受擠壓，離岸邊坡可能變陡，進一步加劇沖刷，對于水庫河床穩定可能會有一定影響。

斷面7位于青草沙水庫的尾閭位置，其剖面的變化相對5、6較為強烈，但總體來看，相對平穩，見圖3(g)。從圖中可以看到距岸不到400 m的水深以淺，近岸邊坡雖然較陡，但基本穩定，水下亦有平臺，為青草沙水庫前沿可能發生的侵蝕提供較大緩沖余地，有利于河床穩定。

整體而言，2000～2013年期間，青草沙水庫前沿河床河槽形態由U型向U型和V型河槽疊加的復式河槽轉變?？拷畮焐喜康暮哟策厼┎粩嗨⑸睿硬坌螒B變化劇烈；中部河床沖淤平衡，河槽形態變化較小，U型河槽持續發育；下部河床表現為U型槽上疊加V型復式河槽，兩翼邊灘被切深，V型河槽向近岸收縮，中部U型河槽向庫堤遷移。

利用斷面數據，進行斷面坡度序列檢驗，得出2000年至2013年7個斷面坡度突變點的位置坐標(表1)，突變點位越多，說明斷面穩定性越差。從表1中可以發現，突變點大多集中在斷面兩側靠近岸灘處。斷面2產生突變點位最多，自2000年至2013年每年都存在，突變位置集中在近岸1 000 m以內，靠近水庫庫堤600 m以內，多達3次突變，最近點離岸140 m，整個斷面形態溝槽交錯。其次是斷面4與斷面5，除個別年份外其余年份均有突變，靠近庫堤400 m以內分別有5次和4次突變。說明水庫庫堤中上部靠近庫堤600 m以淺穩定性較差，有失穩風險，屬于重點監測區域。

表1 斷面坡度突變檢驗Tab. 1 Mutation test of section slope

4 結 語

基于南北港分流口區域2000～2013年期間實測水深數據，分析水庫前沿河床沖淤變化、沙體輸移規律和斷面時空變化，進而探討水庫前沿河床演變以及穩定性，主要結論為：

1)青草沙水庫前沿河床沖淤變化在建庫前后有明顯差異，建庫前總體表現為河槽沖刷、沙洲淤積，呈現準沖-淤振蕩的性質；建庫后變為幅度逐漸減弱的持續沖刷；沿河床河槽形態由U型向U型與V型河槽疊加的復式河槽轉變。

2)青草沙水庫前沿沙體沿落潮主流方向下移，成形沙體呈先增長、再減少，最后趨于平衡，隱形沙體大體上呈減少趨勢。

3)青草沙水庫庫堤前沿近600 m位置是河床失穩的重點風險區域，尤其是水庫庫堤中上部河床處于不穩定狀態，如前沿邊坡所在河槽進一步逼近水庫前沿，則河床面臨進一步沖刷的可能。

[1] SYVITSKI J P M, KETTNER A J, CORREGGIARI A, et al. Distributary channels and their impact on sediment dispersal [J]. Marine Geology, 2005, 222-223: 75-94.

[2] DAI Z J, LIU J T, WEI W, et al. Detection of the Three Gorges Dam influence on the Changjiang (Yangtze River) submerged delta [J]. Scientific Reports, 2014, 4: 6 600.

[3] YANG S L, MILLIMAN J D, XU K H, et al. Downstream sedimentary and geomorphic impacts of the Three Gorges Dam on the Yangtze River [J]. Earth-Science Reviews, 2014, 138: 469-486.

[4] XU K H, MILLIMAN J D. Seasonal variations of sediment discharge from the Yangtze River before and after impoundment of the Three Gorges Dam [J]. Geomorphology, 2009, 104: 276-283.

[5] NAKAYAMA T, SHANKMAN D. Impact of the Three-Gorges Dam and water transfer project on Changjiang floods [J]. Global and Planetary Change, 2013, 100: 38-50.

[6] CHEN X Q, YAN Y X, FU R S, et al. Sediment transport from the Yangtze River, China, into the sea over the Post-Three Gorge Dam Period: A discussion [J]. Quaternary International, 2008, 186: 55-64.

[7] 朱建榮, 顧玉亮, 吳輝. 長江河口青草沙水庫最長連續不宜取水天數[J]. 海洋與湖沼, 2013, 44(5): 1 138-1 145. (ZHU Jianrong, GU Yuliang, WU Hui. Determination of the period not suitable for taking demestic water supply to the Qingcaosha Reservoir near Changjiang River estuary [J]. Oceanologia Et Limnologia Sinica, 2013, 44(5):1 138-1 145. (in Chinese) )

[8] 王曉鵬. 青草沙水庫工程規劃布局節能分析[J]. 城市道橋與防洪. 2012, (8): 216-219. (WANG Xiaopeng. Analysis on energy saving of planning layout of Qingcaosha Reservoir project [J]. Urban Roads Bridges & Flood Control, 2012, (8): 216-219. (in Chinese) )

[9] 宋少紅. 淺析青草沙水庫工程風險因素及對策[J]. 城市道橋與防洪. 2007, (4): 44-46. (SONG Shaohong. Elementary analyses to engineering risk factors of Qingcaosha Reservoir and countermeasures [J]. Urban Roads Bridges & Flood Control, 2007, (4): 44-46. (in Chinese) )

[10] 阮偉，曹慧江，龔鴻鋒. 長江口南北港分汊口河勢控制工程及實施效果研究[J]. 海洋工程, 2011, 29(3): 76-81. (RUAN Wei, CAO Huijiang, GONG Hongfeng. The study of the north-south channel bifurcation control system and its effect [J]. The Ocean Engineering, 2011, 29(3): 76-81. (in Chinese) )

[11] 劉瑋祎，唐建華，繆世強. 長江口北港河勢演變趨勢及工程影響分析[J]. 人民長江, 2011, 42(11): 39-43. (LIU Weiyi，TANG Jianhua，MIU Shiqiang. Study on evolution trend of North Channel of Yangtze River Estuary and influences of related projects [J]. Yangtze River, 2011, 42(11): 39-43. (in Chinese) )

[12] 陸忠民，盧永金，宋少紅，等. 青草沙水庫建設與長江口綜合整治關系研究[J]. 給水排水, 2009, 35(1): 55-58. (LU Zhongmin, LU Yongjin, SONG Shaohong, et al. The relationship between the construction of Qingcaosha Reservoir and the comprehensive regulation of the Yangtze Estuary [J]. Water & Wastewater Engineering, 2009, 35(1): 55-58. (in Chinese) )

[13] 郭超碩, 朱建榮. 青草沙水庫工程對附近水域河床沖淤的影響[J]. 海洋學研究, 2015, 33(3): 34-41. (GUO Chaoshuo, ZHU Jianrong. Impact of Qingcaosha reservoir project on the bed erosion and deposition nearby the water area [J]. Journal of Marine Sciences, 2015, 33(3): 34-41. (in Chinese) )

[14] 黃國玲, 吳焱. 青草沙水源地河勢演變分析[C]//中國水利學會2007年學術年會人類活動與河口分會場論文. 北京:中國水利水電出版社, 2007:180-186. (HUANG Guoling, WU Yan. Analysis of river regime change in Qingcaosha water source [C]//China Hydraulic Engineering Society Annual Meeting, Paper Collection of Human Activity and Estuary Group. Beijing: China Water & Power Press, 2007: 180-186. (in Chinese) )

[15] 鞏彩蘭，惲才興. 應用地理信息系統研究長江口南港底沙運動規律[J]. 水利學報, 2002, (4): 18-22. (GONG Cailan, YUN Caixing. Investigation on sediment transport in south channel of Yangtze River estuary by means of GIS technology [J]. Journal of Hydraulic Engineering, 2002, (4): 18-22. (in Chinese) )

[16] 鞏彩蘭, 惲才興，虞志英. 長江口南港底沙輸移及其對南、北槽分汊口的影響[J]. 海洋工程, 2003, 21(3): 62-67. (GONG Cailan, YUN Caixing, YU Zhiying. Bed-load movement in the south channel and its effect on the south and north passages in the Yangtze Estuary [J]. The Ocean Engineering, 2003, 21(3): 62-67. (in Chinese) )

[17] 惲才興. 從水沙條件及河床地形變化規律談長江河口綜合治理開發戰略問題[J]. 海洋地質動態, 2004, 20(7): 8-14. (YUN Caixing. Comprehensive control and development strategy on the basis of the changes in water and silt conditions and riverbed topographic features [J]. Marine Geology Letters, 2004, 20(7): 8-14. (in Chinese) )

[18] 惲才興. 長江河口近期演變基本規律[M]. 北京：海洋出版社, 2004. (YUN Caixing. Recent developments of the Changjiang Estuary [M]. Beijing: China Ocean Press, 2004. (in Chinese) )

[19] 魏鳳英. 現代氣候統計診斷與預測技術(第2版)[M]. 北京: 氣象出版社, 2007. (WEI Fengying. Modem climatic statistical diagnosis and forecasting technology [M]. Beijing: China Meteorological Press (second edition), 2007. (in Chinese)

[20] 戴志軍, 惲才興, 韓震, 等. 杭州灣北岸漕涇—拓林岸段的海堤預警分析[J]. 海岸工程，2004, 23(4): 23-31. (DAI Zhijun, YUN Caixing, HAN Zhen, et al. Sea wall warning analysis of Caojing-Zhilin bank segment along north coast of the Hangzhou Bay [J]. Coastal Engineering, 2004, 23(4): 23-31. ( in Chinese))

[21] 鄭文燕, 趙德招. 長江口南北港分汊口河段護灘限流工程效果分析[J].中國港灣建設，2010, (5):10-14. (ZHENG Wenyan, ZHAO Dezhao. Effect analysis on shoal protection and discharge control works of south and north channel bifurcation in Yangtze Estuary [J]. China Harbor Engineering, 2010, (5):10-14. (in Chinese))

[22] 戚定滿，顧峰峰，孔令雙, 等. 長江口深水航道整治工程影響數值研究[J]. 水運工程, 2012, (2):90-96. (QI Dingman, GU Fengfeng, KONG Lingshuang, et al. Numerical study on effect of Yangtze estuary deepwater channel regulation project [J]. Port & Waterway Engineering, 2012, (2): 90-96. (in Chinese))

[23] MEI X F, DAI Z J, VAN GELDER P H A J M, et al. Linking Three Gorges Dam and downstream hydrological regimes along the Yangtze River, China[J]. Earth and Space Science, 2015, 2(4):94-106.

[24] DAI Z J, LIU J T. Impacts of large dams on downstream fluvial sedimentation: an example of the Three Gorges Dam (TGD) on the Changjiang (Yangtze River)[J]. Journal of Hydrology, 2013, 480:10-18.

Research on evolution and instability risk of the frontal river bed along the Qingcaosha Reservoir, Changjiang Estuary

SHENG Hao, DAI Zhijun, MEI Xuefei, GE Zhenpeng, LI Shushi，GAO Jinjuan

(State Key Laboratory of Estuarine and Coastal Research, East China Normal University, Shanghai 200062, China)

1005-9865(2017)02-0105-10

TV148

A

10.16483/j.issn.1005-9865.2017.02.015

2016-04-04

國家自然科學基金資助項目(41576087)；2015廣西高等學?？茖W研究人文社科重點資助項目(KY2015ZD133)

盛 皓(1991-)，女，遼寧喀左縣人，碩士研究生。主要從事流域——河口水文及動力地貌等研究。 E-mail：hs_shenghao@163.com

戴志軍。E-mail：zjdai@sklec.ecnu.edu.cn