近十年來(lái)長(zhǎng)江口北支動(dòng)力地貌演化過(guò)程

朱 林，陳曉紅，郭興杰，徐 健

（1. 海洋出版社有限公司，北京 100161；2. 長(zhǎng)江上海航道處，上海 200010；3. 上海市地質(zhì)調(diào)查研究院，上海 200072；4. 上海市水務(wù)規(guī)劃設(shè)計(jì)研究院，上海 200233）

河口位于陸域、海洋和河流的過(guò)渡地帶，擁有豐富的河道、灘涂、濕地等獨(dú)特的河口生態(tài)資源，也是重要的水運(yùn)航道。長(zhǎng)江口呈“三級(jí)分汊，四口入?！被靖窬諿1]。北支位于長(zhǎng)江口第一級(jí)分汊口北汊道，西起崇明島頭，東至入海口，全長(zhǎng)約79 km，平面呈喇叭形，水深較淺，發(fā)育大量的心灘和邊灘[2-3]，也是長(zhǎng)江口重要的支航道之一。歷史上長(zhǎng)江口北支曾是長(zhǎng)江口主要的分流汊道之一，但受河道兩側(cè)人類活動(dòng)和長(zhǎng)期地轉(zhuǎn)偏向力的影響，北支處于不斷淤漲衰退的演化趨勢(shì)[2]，落潮分流比從1915年的15%轉(zhuǎn)變?yōu)?959年的0.7%，并成為潮控型漲潮優(yōu)勢(shì)河道[4-5]，1999年以來(lái)的分流量穩(wěn)定在1%以下[6]。由于長(zhǎng)三角經(jīng)濟(jì)發(fā)展的需求，北支兩岸經(jīng)歷大量圍墾工程[7-8]。許多學(xué)者從北支水沙輸運(yùn)機(jī)制[9-11]、北支沉積物變化特征[8,12-13]、地形地貌演化[14-19]、人類活動(dòng)影響[2,7-8]等方面進(jìn)行了深入研究，但對(duì)近期北支演化特征研究較少，且相關(guān)分析在不同時(shí)期驅(qū)動(dòng)因素也不盡相同，鑒于此，研究近年來(lái)北支河道演化特征具有重要的學(xué)術(shù)意義。

本文利用高分辨率單波束測(cè)深數(shù)據(jù)，并結(jié)合GIS分析方法對(duì)崇明東灘演化過(guò)程進(jìn)行研究，希望解決以下問(wèn)題：近十年來(lái)長(zhǎng)江口北支演化特征；大型工程對(duì)地形演變的影響；長(zhǎng)江口北支的灘涂濕地保護(hù)策略。

1 數(shù)據(jù)來(lái)源與研究方法

本文地形數(shù)據(jù)為長(zhǎng)江上海航道處實(shí)測(cè)，定位系統(tǒng)采用SHCORS系統(tǒng)，利用單波束測(cè)深儀，儀器精度約（1～10 cm）±1%測(cè)深，并結(jié)合GPS接收機(jī)進(jìn)行間距約1 km的斷面測(cè)量，測(cè)深數(shù)據(jù)經(jīng)過(guò)多潮位站修正。坐標(biāo)采用理論最低潮面，北京54坐標(biāo)。平面定位最大和最小誤差分別約±1.99 m和±0.95 m。

采集數(shù)據(jù)通過(guò)ArcGIS軟件中Kriging法進(jìn)行插值和網(wǎng)格化處理，并建立DEM地形模型，在地形模型中提取斷面，地形模型相減形成沖淤模型。數(shù)據(jù)選擇年份為2007年和2017年，為更好地定量對(duì)比長(zhǎng)江口北支的地形變化，本文研究區(qū)域選取固定面積的區(qū)域進(jìn)行對(duì)比研究，具體邊界位置見(jiàn)圖1，選擇固定邊界內(nèi)研究區(qū)定量計(jì)算體積，具體計(jì)算邊界見(jiàn)圖1b，相鄰兩年體積相減，得到體積差，與研究區(qū)面積相除，得到?jīng)_淤厚度，研究區(qū)泥沙容重選用1.4 t/m3[20]，計(jì)算北支沖淤厚度變化。

2 北支地形演化特征

2.1 歷史地形演化

長(zhǎng)江口北支河道的形成和演化與崇明島的發(fā)育有較大關(guān)系，歷史上北支大部分時(shí)期處于自然演化態(tài)勢(shì)，人類活動(dòng)影響較小，1915年北支上口0 m等深線河寬約5.8 km，-10 m等深線基本覆蓋北支河段，分流比約占長(zhǎng)江口整體徑流量的25%；1954年由于長(zhǎng)江大洪水的造床作用，南支河道被沖深加寬，南支分流比增加，而北支分流角繼續(xù)增大，徑流分流比減少，20世紀(jì)60年代后北支兩側(cè)圍墾工程增加，加速了北支的萎縮程度，北支逐漸由河控型河道變化為潮控型河道，徑流分流比逐漸下降至1%左右，潮流帶來(lái)的泥沙逐漸淤堵河道，河道逐年淤淺[18,21]。

圖1 研究區(qū)域示意圖Fig.1 Schematic diagram of the study area

2.2 地形沖淤變化

圖2 2007～2017年長(zhǎng)江口北支沖淤變化Fig.2 Erosion and deposition changes in the North Branch of Yangtze River Estuary from 2007 to 2017

從圖2長(zhǎng)江口北支沖淤變化結(jié)果可見(jiàn)，2007～2017年地形變化可分為3個(gè)不同河段，分別為中上段、下段、口門(mén)段。北支中上段河道較窄，基本以侵蝕為主，局部河段淤積，近十年間沖淤變化在0.5～1 m之間，其中侵蝕范圍明顯大于淤積范圍，最大侵蝕或者淤積區(qū)發(fā)生于靠岸一側(cè)，侵蝕深度在1～2 m之間；中下段為淤積河段，此河段河寬逐漸加寬，發(fā)生明顯的淤積，北側(cè)淤積大于南側(cè)，淤積厚度在0.5～1 m之間，河道中間形成了一個(gè)窄的侵蝕條帶，條帶侵蝕深度在1 m以上；口門(mén)河段呈“北沖南淤”態(tài)勢(shì)，存在明顯的沖淤分界線，侵蝕區(qū)范圍略大于淤積段范圍，北側(cè)侵蝕深度在0.5～1 m之間，而南側(cè)淤積厚度大部分大于1 m。

2007年長(zhǎng)江口北支河道容積約為0.66 km3，2017年河道容積增大，約為0.67 km3，研究區(qū)十年內(nèi)侵蝕量為13.97×106t，研究區(qū)面積為218.61 km2，年均侵蝕速率為0.46 cm。

2.3 典型斷面變化特征

在研究區(qū)內(nèi)選擇5個(gè)典型斷面進(jìn)行分析，其中AA’斷面在北支上口門(mén)，分析河道徑流入口斷面變化，BB’和CC’斷面都在北支河道轉(zhuǎn)折位置，在轉(zhuǎn)折角兩側(cè)分別進(jìn)行斷面分析，DD’斷面位于淤積區(qū)，位于崇啟大橋西側(cè)位置，EE’斷面位于北支下口門(mén)位置，監(jiān)測(cè)口門(mén)斷面變化。

AA’斷面兩側(cè)變化相反，江蘇省滸通港一側(cè)發(fā)生明顯淤積，岸邊最大淤積厚度約8 m，而崇明島一側(cè)發(fā)生侵蝕，最大侵蝕厚度在3 m左右（圖3a）；BB’斷面在江蘇省中興村一側(cè)發(fā)生侵蝕，侵蝕深度在3 m左右，而崇明島一側(cè)發(fā)生淤積，淤積厚度在1.5 m左右（圖3b）；CC’在靠近江蘇省崇海鎮(zhèn)一側(cè)發(fā)育有一條深泓，在2007～2017年間深泓發(fā)生侵蝕增加，最大侵蝕深度在3 m左右，深泓寬度也增寬，而崇明島一側(cè)沒(méi)有明顯的沖淤變化（圖3c）；DD’斷面位于淤積區(qū)，在2007年河道中間發(fā)育有一個(gè)心灘，但是在2017年，心灘基本消失，周圍發(fā)生明顯淤積，淤積厚度在1 m左右（圖3d）；EE’斷面位于北支下口門(mén)，靠近江蘇省江海村一側(cè)發(fā)生明顯侵蝕，侵蝕深度在1.5 m左右，而崇明島一側(cè)發(fā)生淤積，淤積厚度在0.5 m左右，整體侵蝕略大于淤積（圖3e）。

3 北支演化驅(qū)動(dòng)因素

從上述結(jié)果可知，2007～2017年間長(zhǎng)江口北支近年來(lái)整體呈侵蝕態(tài)勢(shì)，局部中上段發(fā)生侵蝕，下段發(fā)生淤積，口門(mén)段發(fā)生“北沖南淤”變化特征。一般而言，河勢(shì)的穩(wěn)定是由于水流、泥沙和地形長(zhǎng)期相互作用的結(jié)果，任何一個(gè)因素的改變都會(huì)打破河流原有的態(tài)勢(shì)。近年來(lái)長(zhǎng)江口演化受人類活動(dòng)影響明顯，人類活動(dòng)打破原有的河流演化平衡，后河道經(jīng)過(guò)自適應(yīng)演化再次達(dá)到新的平衡，整個(gè)演化過(guò)程為平衡—不平衡—平衡[3,22]。

3.1 泥沙來(lái)源及與地形沖淤的關(guān)系

圖3 2007～2017年長(zhǎng)江口北支典型斷面水深變化Fig.3 Water depth changes of typical sections in the North Branch of Yangtze River Estuary from 2007 to 2017

動(dòng)力和泥沙來(lái)源是決定河道演化的主要因素，北支的水沙來(lái)源一般為長(zhǎng)江徑流輸入和近海沿岸流輸送。由于1954年大洪水對(duì)南支的沖刷作用，南支河道擴(kuò)大，基本成為長(zhǎng)江徑流的主要分流口，加之20世紀(jì)60年代后，北支兩側(cè)大量圍墾工程的建設(shè)，北支的徑流分流比減少至1%左右，北支河道與長(zhǎng)江主河道夾角基本接近90°（圖4），所以減少了北支中徑流攜帶的泥沙量，但流域輸沙仍然是北支重要的泥沙來(lái)源之一。由于長(zhǎng)江上游攔水大壩工程、南水北調(diào)工程和一些水土保持工程的實(shí)施，長(zhǎng)江泥沙輸運(yùn)量減少至過(guò)去的28%[3,23]，長(zhǎng)江口整體懸沙含量較低，而根據(jù)趙捷等[4]分別于2009年和2010年在長(zhǎng)江口徐六徑、南支崇西、北支永隆沙懸沙觀測(cè)可知，長(zhǎng)江口徐六徑懸沙含量最低，北支永隆沙懸沙含量遠(yuǎn)大于南支崇西，最大懸殊在10倍以上，且北支口外存在蘇北沿岸流，蘇北沿岸流搬運(yùn)廢黃河三角洲泥沙，輸送至長(zhǎng)江口[24]，在漲潮作用的輸送下，進(jìn)入北支，導(dǎo)致北支淤積變淺，是北支淤積的另一個(gè)重要物源之一。

圖4 長(zhǎng)江口北支上口分流角特征Fig.4 Diversion angle’s characteristics of upper entrance in the North Branch of Yangtze River Estuary

北支中上段以侵蝕為主，下段以淤積為主，根據(jù)北支的兩個(gè)泥沙來(lái)源，徑流輸沙減少可能是導(dǎo)致北支中上段侵蝕的主要原因之一，而沿岸流輸沙為北支下段淤積提供了物質(zhì)基礎(chǔ)。

3.2 河道形態(tài)影響

河道形態(tài)也是決定地形演化的主要因素，根據(jù)前述結(jié)果可知，北支中上段沖刷，下段淤積，中上段河寬較窄，下段逐漸變寬，動(dòng)力減弱，漲潮流泥沙在中下段大量淤積，上溯的過(guò)程中，隨著河寬變窄，流速增大，動(dòng)力增強(qiáng)，泥沙不容易落淤，所以導(dǎo)致中上段侵蝕大于淤積。

長(zhǎng)江口北支上段基本為侵蝕態(tài)勢(shì)（圖5a），在北支河道轉(zhuǎn)角位置也出現(xiàn)典型的“凹岸沖刷、凸岸淤積”現(xiàn)象（圖5b），局部河道由于從窄河道進(jìn)入寬河道后發(fā)生淤積現(xiàn)象。

圖5 長(zhǎng)江口北支重要河段沖淤變化Fig.5 Erosion and deposition changes of the important sections in the North Branch of Yangtze River Estuary

3.3 地轉(zhuǎn)偏向力影響

北支河道主要由漲潮流控制[7,19]，尤其是口門(mén)位置，漲潮流動(dòng)力較強(qiáng)，在地轉(zhuǎn)偏向力的作用下，漲潮流向右偏轉(zhuǎn)，動(dòng)力較強(qiáng)，導(dǎo)致江蘇省一側(cè)發(fā)生侵蝕，而崇明島一側(cè)主要由落潮流影響，北支落潮動(dòng)力較弱，發(fā)生淤積現(xiàn)象，而南側(cè)為崇明東灘位置，近年來(lái)崇明東灘基本為淤積態(tài)勢(shì)[3]，河道寬闊、地轉(zhuǎn)偏向力和多汊道交匯共同作用導(dǎo)致南側(cè)發(fā)生淤積現(xiàn)象。

3.4 人類活動(dòng)影響

人類活動(dòng)已成為影響長(zhǎng)江口地形演化的主要驅(qū)動(dòng)因素之一，灘涂圍墾是典型的長(zhǎng)江口人類活動(dòng)，也是河口演化的主要驅(qū)動(dòng)力。崇明北側(cè)邊灘是長(zhǎng)江口重要的圍墾區(qū)之一，也是上海市主要的土地來(lái)源之一，從1958～2005年間，北支兩側(cè)進(jìn)行了大量的圍墾工程（表1），共計(jì)圍墾491.4 km2。本文研究時(shí)間為2007～2017年，通過(guò)文獻(xiàn)[3]遙感影像可知，崇明北沿在2011～2017年間北支兩側(cè)仍然進(jìn)行了大量圍墾工程。

表1 長(zhǎng)江口北支灘涂圍墾統(tǒng)計(jì)Table 1 Statistics of tidal flat reclamation within the North Branch of Yangtze River Estuary

北支兩側(cè)的圍墾工程極大的束窄了河道河寬，也減弱了河道的擺動(dòng)幅度，加速了河道的淤積萎縮程度，其次口門(mén)位置河寬較大，2007年存在一條明顯的潮流沙脊，此區(qū)域?yàn)槊黠@的漲落潮主流路不同的復(fù)式河道，由于北支口門(mén)南側(cè)圍墾工程的實(shí)施[3]，束窄了河道，此區(qū)域發(fā)生明顯的淤積現(xiàn)象，但是河道中間的潮流沙脊基本消失（圖6）。

圖6 2007～2017年間長(zhǎng)江口北支下段潮流沙脊消失特征Fig.6 Disappearing characteristics of tidal sand ridge at the lower section in the North Branch of Yangtze River Estuary from 2007 to 2017

北支在2007～2017年間發(fā)生整體侵蝕，與過(guò)去淤積[7,17-19]演化趨勢(shì)不同，北支口門(mén)位置崇明東灘發(fā)生淤積，且其主要沙源來(lái)自海域來(lái)沙，且北支口內(nèi)懸沙量并未發(fā)生明顯減少，但是整體侵蝕卻大于淤積，可能與近期周圍大量圍墾工程相關(guān)，圍墾工程就近河道采集泥沙進(jìn)行促淤，導(dǎo)致長(zhǎng)江口北支發(fā)生持續(xù)侵蝕，在長(zhǎng)江大保護(hù)的號(hào)召下，北支兩側(cè)圍墾工程將逐漸減少，未來(lái)北支受海域來(lái)沙的長(zhǎng)期影響，淤積衰退仍然是北支的主要演化趨勢(shì)。

4 北支灘涂濕地及其保護(hù)建議

過(guò)去長(zhǎng)江流域輸沙量豐沛[1]，且自南北支分流口工程的建設(shè)和北支的圍墾工程的開(kāi)展，北支過(guò)去一直為淤淺的發(fā)育態(tài)勢(shì)，所以發(fā)育有大量的灘涂濕地。20世紀(jì)70年代，特別是2003年三峽工程后長(zhǎng)江入海泥沙通量急劇下降[23]，2007～2017年長(zhǎng)江輸沙率為128×106t/a左右，較三峽工程前減少約70%[25]，長(zhǎng)江口北支發(fā)生“中上段侵蝕，下段淤積”的發(fā)育態(tài)勢(shì)，整體侵蝕大于淤積，極大地影響到北支灘涂濕地的發(fā)育[26]。

北支的灘槽地形受到流域輸沙減少和工程采砂的雙重脅迫，近期發(fā)生侵蝕，為保護(hù)長(zhǎng)江口灘涂濕地，首先，沿江省市需堅(jiān)持“長(zhǎng)江大保護(hù)”國(guó)家戰(zhàn)略，減少北支河槽的工程采砂；其次，由于北支兩側(cè)分別由江蘇和上海分別管理，在北支的保護(hù)中需上海和江蘇協(xié)同進(jìn)行，借助“長(zhǎng)三角一體化發(fā)展戰(zhàn)略”，兩省市協(xié)同管理和保護(hù)北支的灘涂濕地。

5 結(jié)論

本文收集2007～2017年間高精度實(shí)測(cè)地形數(shù)據(jù)，并利用GIS分析方法研究長(zhǎng)江口北支的演化特征，結(jié)合北支周圍水沙輸送、河道形態(tài)和圍墾工程等影響因素進(jìn)行分析，得出以下結(jié)論：

（1）長(zhǎng)江口北支在2007～2017年間整體發(fā)生侵蝕，其中，中上段侵蝕明顯，下段輕微淤積，入?？陂T(mén)呈“北沖南淤”的特征，期間共侵蝕泥沙量約13.97×106t，年均侵蝕厚度約0.46 cm，侵蝕峰值發(fā)生于北支上段，年均最大侵蝕為10～20 cm。

（2）影響北支演化的主要泥沙物源來(lái)自長(zhǎng)江流域和蘇北沿岸流的泥沙輸送，河道轉(zhuǎn)角位置發(fā)生明顯的“凹岸沖刷、凸岸淤積”現(xiàn)象，豐沛的懸沙和較弱的動(dòng)力條件是導(dǎo)致北支下段淤積的主要原因。

（3）北支兩側(cè)的圍墾工程對(duì)河道的演化影響日趨增大，長(zhǎng)江口北支上口分流角變化是引起徑流動(dòng)力減弱的主要原因，大量圍墾工程的河道采集泥沙可能是引起研究期間河道侵蝕的主要原因。

（4）建議未來(lái)堅(jiān)持“長(zhǎng)江大保護(hù)”和“長(zhǎng)三角一體化”國(guó)家戰(zhàn)略，減少北支周圍的工程采砂，長(zhǎng)三角省市協(xié)同保護(hù)北支的灘涂濕地。