長江口徐六涇斷面近底懸沙測驗及影響分析

劉 傳 杰，朱 巧 云

(長江水利委員會水文局 長江口水文水資源勘測局，上海 200136)

0 引 言

長江口為中等強度的潮汐河口，在徑流、潮流、波浪、鹽水楔以及增減水等多種變化的動力因素作用下，河口的泥沙運動規(guī)律及其分布較為復雜[1]。而近底懸沙濃度主要反映水體懸沙和河床床面泥沙的交換及與水動力的關系，研究其較為真實的分布及運動規(guī)律，不僅直接關系到取、排水，防、排沙工程的規(guī)劃設計，而且也是泥沙研究的重要方面。常規(guī)水文觀測方法測得的底部含沙量一般只相當于距底70～120 cm高度含沙量，為了有效解決這個問題，多年來許多專家和學者對此做了大量的工作，如方鐸和王貴道[2-3]研制了NBS-84型臨底懸沙采樣器，主要應用于5 m以下水深的河道測驗。近年來，以華東師范大學何青等為代表的專家學者利用先進的聲學和光學儀器組建四角架近底觀測系統(tǒng)，在長江口南、北槽河段進行了應用研究[4-8]。此系統(tǒng)雖然可獲取常規(guī)水文觀測難以獲取的近底溫鹽沙資料，但沒有自動采水系統(tǒng)，不能同步采集同高度的懸沙水樣，近底懸沙濃度是采用上層水樣與濁度計(OBS)獲取的水體中懸浮物濁度率定推求，可OBS的率定曲線很大程度上依賴采樣的環(huán)境，如懸浮水體的粒徑、水溫等，尤其是懸浮水體的粒徑[9-10]，所以采用上層水樣與近底層濁度率定曲線推求出的近底懸沙濃度存在一定的誤差。

本次研究是在洪、枯季同步水文測驗期間，應用一種新型的近底懸沙采樣系統(tǒng)——雙尾翼近底懸沙采樣系統(tǒng)，在長江口徐六涇斷面的淺灘和深槽的代表垂線上同步采集距河底50 cm和10 cm的近底水樣；并根據(jù)現(xiàn)場實測資料分析含沙量垂向分布規(guī)律，以及近底懸沙濃度對測驗成果精度的影響程度。

1 研究區(qū)域概況

徐六涇斷面位于長江口徐六涇節(jié)點段，位置如圖1 所示。所在河段上承澄通河段通州沙汊道段，下連長江口北支和白茆沙汊道段。通州沙汊道段上起十三圩，下至徐六涇，全長約39.0 km，為多灘多分汊河道。白茆沙汊道段上起白茆河口，下至七丫口，全長35.5 km，大部分徑流從白茆沙水道下泄至南支主槽，小部分徑流分流北支，該河段同樣是復雜的多汊分流河段。

徐六涇節(jié)點段進口河寬約4.6 km(2011年新通海沙圍墾前為5.7 km)，斷面形狀呈不對稱“W”形(見圖2)，主槽寬約2.35 km(-10 m等高線間距，1985國家高程基準)[11]，主流偏南，斷面主槽多年來基本保持穩(wěn)定。根據(jù)以往的實測數(shù)據(jù)分析：徐六涇斷面主槽以南底質(zhì)為淤泥質(zhì)亞黏土，抗沖擊性較強；主槽以北為粉砂，易起動[12]。懸沙以細粒組分為主，多為極細粉砂和黏土；枯季懸沙中值粒徑粗于洪季；洪季懸移質(zhì)含沙量高于枯季[13]。

2 現(xiàn)場測驗及方法

2.1 近底懸沙采樣系統(tǒng)結構及工作原理

為了能在徐六涇斷面獲取近底懸沙含沙量，采用了新型雙尾翼近底懸沙采樣系統(tǒng)(見圖3)。受船舶中的雙舵設計啟發(fā)，該采樣系統(tǒng)采用雙尾翼設計，在高流速的紊流中，采樣系統(tǒng)縱向穩(wěn)定性極好，不左右擺動，原地不旋轉(zhuǎn)、不轉(zhuǎn)圈，采用流線型雙尾翼設計，由上、下2個橫式采樣器組成，上層采樣器距離底部垂直距離為50 cm，下層采樣器距離底部10 cm，利用自重為動力，在底部裝有觸底關閉裝置，采用精密軸承傳動將水樣容積倉自動關閉。

現(xiàn)場取樣時，以開艙的狀態(tài)入水，當采樣器下放到河底時，底部托盤受到重力的反作用力，向上頂開容積倉門的鎖緊裝置，在彈簧的拉力作用下，同步一次性關閉上下兩層橫式采樣器。此采樣系統(tǒng)優(yōu)點是：提高倉門關閉的成功率，同時確保了水樣采集位置的準確性，保證了水樣的代表性；缺點是系統(tǒng)自身重量較重，測船必須安裝有專用絞車才能使用。

2.2 測驗方法

雙尾翼近底懸沙采樣系統(tǒng)于2019年6月開始在徐六涇斷面進行現(xiàn)場測試(見圖4)。首次使用時，由于近底采樣系統(tǒng)底盤較小，配重不夠，河床底質(zhì)較松軟等原因，在斷面深槽處出現(xiàn)多次觸底不能閉合的情況[14]。經(jīng)過多次改進，新型雙尾翼近底懸沙采樣系統(tǒng)可成功運用于徐六涇斷面的近底懸沙取樣。

為了研究徐六涇斷面淺灘和深槽不同位置懸移質(zhì)含沙量沿垂線分布規(guī)律，現(xiàn)場測驗時在斷面上選取1號(平均水深約15 m)和4號(平均水深約45 m)，這兩條垂線分別代表淺灘和深槽，具體位置見圖2?，F(xiàn)場測驗時，采用常規(guī)六點法與近底兩點相結合采樣，每條垂線沿垂線取8層水樣，分別為水下0.5 m、相對水深0.2H層、0.4H層、0.6H層、0.8H層、水深H-0.7 m層(常規(guī)法的底層)、水深H-0.5 m層和水深H-0.1 m層。采用焙干稱重法直接得出各分層含沙量。

2.3 測驗實施

徐六涇斷面具有雙向水沙運動特性，懸移質(zhì)含沙量有明顯的漲落潮、大小潮等的潮相變化，還有洪季與枯季季節(jié)性的變化，考慮到分析樣本的代表性，安排在洪季、枯季各進行一次全潮測驗?，F(xiàn)場具體測驗實施如下：在淺灘1號垂線處，利用常規(guī)法與近底懸沙采樣系統(tǒng)進行了洪季和枯季2次測驗；在深槽4號垂線處，利用常規(guī)法進行了洪、枯季2次測驗，近底懸沙采樣系統(tǒng)只在枯季進行了1次測驗。洪、枯季均開展大、中、小潮完整潮周期的測驗。具體測驗時間見表1。

表1 測驗時間Tab.1 The timetable of the observations

3 測驗成果分析

3.1 懸移質(zhì)顆粒級配分析

3.1.1洪、枯季

洪、枯季期間，淺灘1號和深槽4號垂線懸移質(zhì)顆粒特征值見表2～3。由表2～3可見：1號垂線懸移質(zhì)顆粒粒徑沿垂線分布，枯季較洪季存在明顯的粗細變化，洪季4種流態(tài)下，顆粒粒徑沿垂線分布變化不明顯，中值粒徑基本在8～12 μm之間；枯季急流時，中上層的泥沙中值粒徑變化不明顯，而越靠近河床泥沙粒徑越粗，近底層粗顆粒泥沙(中值粒徑0.062 mm以上)比重達70%以上，中值粒徑為0.121 mm，約是表層中值粒徑的8倍多；憩流時由于水流動力減弱，粒徑粗細變化較急流時減小，近底層中值粒徑約為表層的1.5～3.5倍。深槽處4號垂線與淺灘1號垂線不同，洪、枯季粒徑粗細差異不大，整個剖面分布較均勻，洪季中值粒徑基本在10～12 μm之間，沿垂線分布粗細變化不明顯；枯季粒徑較洪季粗一些，中值粒徑基本在12～19 μm之間，近底層顆粒較粗。

表2 淺灘1號垂線懸移質(zhì)顆粒特征值Tab.2 Characteristic values of the suspended sediment in profile 1 on the shoal

表3 深槽4號垂線懸移質(zhì)顆粒特征值Tab.3 Characteristic values of the suspended sediment in profile 4 in the deep channel

3.1.2不同潮型

表4和表5分別為1號和4號垂線懸移質(zhì)顆粒大、中、小潮期表層和近底層(4種流態(tài))的顆粒特征值。由表可見：1號垂線大潮和中潮期，顆粒粒徑沿垂線分布存在明顯的粗細變化，即越靠近河床泥沙粒徑越粗；而小潮期顆粒粒徑沿垂線分布粗細比較均勻，中值粒徑基本在9～20 μm之間。深槽處4號垂線，大、中、小潮期，顆粒粒徑沿垂線分布粗細均比較均勻，中值粒徑基本在10～19 μm之間。淺灘與深槽處對于不同潮型，總體上都是小潮期粒徑最細，大潮期粒徑較粗；4種流態(tài)比較，急流時粒徑較憩流時粒徑粗。

表4 不同潮型下淺灘1號垂線懸移質(zhì)顆粒特征值Tab.4 Characteristic values of the suspended sediment in profile 1 on the shoal under different tidal types

表5 不同潮型下深槽4號垂線懸移質(zhì)顆粒特征值Tab.5 Characteristic values of the suspended sediment in profile 4 in the deep channel under different tidal types

3.2 含沙量剖面分析

3.2.1洪、枯季

圖5為1號垂線洪季(2020年7月)、枯季(2021年3月)含沙量剖面圖(z為距河底距離，下同)。從圖5可知：含沙量沿垂線分布基本上與顆粒粒徑一致，中上層區(qū)域含沙量變化比較小，越靠近河床含沙量越大，使得整個含沙量剖面形成“L”形分布，漲、落急時更顯著。通過懸沙顆粒粒徑分析可知，急流時刻的泥沙粒徑要明顯粗于憩流時刻，急流時粗顆粒泥沙集中于近河底部分，但是由于水體紊動強度不足，部分粗粒徑級泥沙只能懸浮到某一高度，根本到不了上層區(qū)域，從而造成了含沙量剖面形成“L”形分布，同時當水流動力逐步減弱至憩流時刻，底層較粗、沉降速度大的泥沙重新落淤到床面，造成近底層含沙量迅速減少，上層區(qū)域由于泥沙粒徑細沉降速度小含沙量變化不大。洪、枯季規(guī)律相似。

3.2.2不同潮型

圖6和圖7分別為1號和4號垂線枯季(2021年3月)大、中、小潮期4種不同流態(tài)下含沙量剖面圖。由圖6～7可知：不同潮型含沙量剖面分布規(guī)律與上述分析一致。大、中潮期，1號垂線漲、落急含沙量沿垂線分布存在明顯的梯度，憩流時含沙量沿垂線變化比較小；小潮期，由于動力較弱，泥沙粒徑進一步細化，含沙量比較小，使得含沙量沿垂線分布更加趨于均勻。4號垂線含沙量明顯小于1號，尤其在近底層，無論是急流時刻還是憩流時刻都遠小于1號垂線，且無論大、中、小潮期含沙量沿垂線分布均比較均勻。

3.3 近底懸沙影響分析

3.3.1垂線平均含沙量計算方法

(1) 常規(guī)六點法。垂線平均含沙量通常采用垂線上各分層流速加權計算[15]，計算公式如下：

(1)

(2) 加入近底2點(H-0.5 m、H-0.1 m)懸沙的垂線平均含沙量按平均流量法計算。沿水深方向?qū)⒋咕€離散為m個微小單寬面積，根據(jù)垂線流速分布曲線、含沙量分布曲線，分別計算出q和qs。計算公式如下：

(2)

式中：q為垂線單寬流量；qs為垂線單寬輸沙率。

3.3.2誤差統(tǒng)計方法

相對誤差：

(3)

系統(tǒng)誤差：

(4)

標準差：

(5)

3.3.3垂線平均含沙量誤差分析

根據(jù)現(xiàn)場實測資料，用常規(guī)六點法計算值與計入近底懸沙后計算的垂線平均含沙量擬合關系曲線，如圖8～9所示。淺灘1號垂線關系點成散射狀，含沙量越大，偏離相對越大(見圖8)，但總體上兩者呈高度顯著相關。相比較而言，在含沙量及含沙量垂向分布梯度均較小的深槽4號垂線(見圖9)，無論漲、落潮，相關系數(shù)均達0.99以上，即計入近底懸沙后計算的垂線平均含沙量與常規(guī)六點法的計算值基本相當。

根據(jù)數(shù)據(jù)統(tǒng)計分析可知：常規(guī)法計算的垂線平均含沙量偏小，尤其在淺灘處，相對誤差大于10%的樣本數(shù)中漲潮為26.3%、落潮為28.6%，漲、落系統(tǒng)誤差潮分別為-6.50%，-5.63%，標準差在10%左右，即常規(guī)法計算的垂線平均含沙量無論漲、落潮時均偏小。深槽處，由于泥沙粒徑比較細，含沙量垂向分布比較均勻，80%以上的樣本數(shù)相對誤差小于5%，系統(tǒng)誤差基本為0，標準差在3%以內(nèi)，即深槽處近底懸沙對垂線平均含沙量影響比較小。

3.3.4潮平均含沙量誤差分析

表6～7分別為淺灘1號和深槽4號垂線潮平均含沙量及誤差統(tǒng)計，表中以計入近底懸沙后的值為“真值”。根據(jù)表中的結果分析可知：計入近底懸沙后，垂線潮平均含沙量將增大，增大幅度與含沙量垂向分布密切相關，根據(jù)前述章節(jié)分析，淺灘1號垂線處含沙量垂向梯度大，且枯季比洪季大，則枯季相應的增大幅度較洪季大，洪、枯季增幅平均分別為5.65%，9.08%；落潮增幅大于漲潮；大、中潮比小潮垂向梯度大，則相應的增幅也大。深槽4號垂線處無論大、中、小潮期，含沙量沿水深的分布均相對較均勻，計入近底懸沙后，漲、落潮潮平均含沙量增幅基本在1%以內(nèi)，即深槽處近底懸沙對潮平均含沙量影響也比較小。

表6 淺灘1號垂線潮平均含沙量及誤差統(tǒng)計Tab.6 Vertical tide-averaged suspended sediment concentration and the error statistics in profile 1 on the shoal

表7 深槽4號垂線平均含沙量及誤差統(tǒng)計Tab.7 Vertical tide-averaged suspended sediment concentration and the error statistics in profile 4 in the deep channel

4 結果討論

在懸移質(zhì)運動過程中，貫穿著重力作用與紊動擴散作用的矛盾統(tǒng)一。根據(jù)Rouse方程，懸移質(zhì)泥沙垂向分布的不均勻程度取決于Z值[16](Z=w/(ku*)，w為泥沙沉降速度，u*為摩阻流速，k為卡門常數(shù))。在一定時間內(nèi)水力條件變化不大的情況下，Z值與顆粒的沉速成正比。細顆粒的Z值小，含沙量沿水深的分布相對較均勻；粗顆粒的Z值大，泥沙將集中于近河底部分，部分粗粒徑級泥沙只能懸浮到某一高度，到不了上層區(qū)域，也就是說，在河底各粒徑級的泥沙都存在，且較粗粒徑級含量大，故含沙量大；在上層，只有部分細顆粒徑級的懸浮泥沙，故含沙量比較小。

通過上文比較分析，徐六涇斷面淺灘和深槽含沙量垂向分布規(guī)律不同。淺灘漲落急流時刻的泥沙粒徑要明顯大于憩流時刻，越靠近河床泥沙粒徑越大，其中近底70 cm處的泥沙要明顯比近底50 cm處的細，而中上層的泥沙中值粒徑變化則不再明顯，這一現(xiàn)象說明在徐六涇斷面泥沙從床面泥沙受到急流侵蝕進入水體后，在70 cm以內(nèi)形成比較顯著的含沙量梯度，但是由于水體紊動強度不足，占比較大的粗顆粒泥沙往往無法進入床面1 m以上的水體中，從而造成了整個含沙量剖面呈現(xiàn)出“L”形分布，同時當水流動力逐步減弱至憩流時刻，底層較粗、沉降速度大的泥沙重新落淤到床面，造成近底層70 cm以內(nèi)的含沙量迅速減少，不過床面1 m以上的水體由于泥沙粒徑細，沉降速度小含沙量變化不大。深槽含沙量垂向分布比較均勻，深槽泥沙粒徑在15 μm左右且變化比較小，這與淺灘上下層粒徑存在明顯梯度是不同的，這些都表明深槽的床沙與懸沙交換不多，此處的懸沙主要以過境泥沙為主。含沙量垂向分布規(guī)律不同對測驗成果影響程度也不同，淺灘近底層含沙量大，對測驗成果的影響就比較大；而深槽近底層含沙量與中上層水體相當，對測驗成果的基本沒有影響。洪、枯季含沙量垂向分布均符合于上述規(guī)律和特點。

5 結 論

本文采用雙尾翼近底懸沙采樣系統(tǒng)，觀測得到了徐六涇斷面淺灘和深槽的洪、枯季全潮實測資料，并通過計算、分析得出如下幾點結論。

(1) 雙尾翼近底懸沙采樣系統(tǒng)，采用雙尾翼設計，在長江口河段往復流的復雜水流條件下縱向穩(wěn)定性好，確保了采樣位置的準確性，為測驗成果精度提供了重要支撐。

(2) 徐六涇斷面淺灘和深槽含沙量垂向分布規(guī)律不同。含沙量垂向分布的不均勻程度取決于泥沙顆粒特性，淺灘泥沙粒徑較粗且越靠近河床泥沙粒徑越粗，含沙量沿垂線分布基本上與顆粒粒徑一致，中上層水體含沙量變化較小，越靠近河床含沙量越大，含沙量剖面呈“L”形分布。深槽泥沙粒徑比較細，含沙量垂向變化比較小，含沙量剖面分布比較均勻。

(3) 近底懸沙對測驗成果的影響程度與取樣位置有關。計入近底懸沙后，淺灘處平均增幅為7.37%，深槽處平均增幅小于1%，淺灘處影響程度明顯大于深槽。

(4) 近底懸沙對測驗成果的影響程度在各個時間尺度上有差異。計入近底懸沙后，含沙量一般情況下是增大的，但增幅不同，枯季大于洪季，落潮大于漲潮，大潮大于小潮。如淺灘處平均含沙量在洪、枯季增幅分別為5.65%，9.08%，枯季是洪季的1.6倍；2021年3月枯季淺灘處漲、落潮平均含沙量增幅分別為8.98%，9.19%。