基于ADCP的長江河口推移質運動特性研究

陸雪駿，程和琴，胡浩，郭興杰，劉高偉，吳帥虎

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

推移質運動是河流中泥沙運動的普遍形式之一，對地貌演變有重要影響，所以理解推移質運動對實際生產有著重要的理論指導意義。由于推移質運動處于水底床面附近，以往受限于觀測方法和儀器條件，觀測難度較大、耗資巨大，成為人類工程建設的一個技術難點。

傳統推移質觀測由于采樣器的阻水作用，改變了床面運行的推移泥沙的水力條件，造成偏差（張瑞瑾等，1998），隨著測量技術的發展，利用新型的測量儀器嘗試觀測推移質運動成為可能。國外Rennie 等（2002） 首先利用ADCP（聲學多普勒流速剖面儀） 的底跟蹤（Bottom Track） 功能觀測河床推移質視速度（va），發現與傳統采樣器測的推移質輸沙率有高度相關性，確認ADCP 觀測推移質運動的可行性（Colin，2002；Colin，2004）；國內吳中等（2002） 建立ADCP 的底跟蹤和GPS 信號的聯系，證實了利用ADCP 底跟蹤與GPS 定位差距分析底沙運動的理論正確、方法可行（吳中等，2002）。ADCP 觀測推移質方法優于傳統采砂器測量方法，采用非侵入式手段，不破壞現場水流結構。在實際測量應用中，把ADCP 測得推移質運動視速度近似看作推移質運動速度（Kostaschuk，2005），對了解和掌握長江口推移質運動規律有重要意義。

本文采用定點測量方法，利用ADCP 對長江口砂質河床推移質運動進行觀測，分析了推移質運動在漲、落潮周期內的變化情況。

1 數據采集

ADCP 的實測資料來自筆者所在課題組在長江口S1（長興島北側）、S2（南北槽分流口）、S3（崇明東灘南側）、S4（橫沙島北側）、S5（南匯南灘） 和S6（橫沙通道北側） 共6 個測點的定點水文觀測（圖1） 資料，時間尺度為二個潮周期。S1、S3、S4、S5 測點于2012年6月6-7日測量，S2 測點于2013年7月10-11日測量，S6 測點于2013年7月1-2日測量。

此次數據采集使用的是美國RDI 公司生產的600 kHz 型ADCP，具體由WinRiver 軟件對ADCP的換能器實現操控和采樣，采樣時間間隔設置為5 s，換能器入水深度為1 m。同時，采用高精度差分GPS 和外部GPS 羅經進行定位。

底層泥樣采集采用帽式采集器，并用聚乙烯塑料袋密封存放。在實驗室內，用超聲波震蕩分散，使用MASTER SIZER 2000 型激光粒度分析儀經行粒度分析。

圖1 長江口自然地理概況及測站位置

2 觀測原理

ADCP 是目前世界上先進的流速流量測量設備，其利用多普勒原理探測水體中微顆粒的運動速度而推測水體流速。ADCP 提供兩種方法測量安裝平臺的速度：一是“底跟蹤”技術，ADCP 通過接收和處理來自水底的回波信號而計算得到ADCP 的安裝平臺與底床的相對速度（vBT）；二是利用GPS全球衛星定位系統測量安裝平臺的速度（vDGPS）（趙勝凱 等，2007；楊俊輝，2009）。在測量時，如果床面沒有推移質運動，那么ADCP 的安裝平臺相對于底床的速度等于安裝平臺參照GPS 得到的速度，即：

圖2 各測點漲落急階段航跡圖（底跟蹤航跡為三角形，GPS航跡為正方形）

如果是存在推移質運動的情況下，那么ADCP的安裝平臺相對于水底的速度和安裝平臺參照GPS得到的速度之間存在偏差，差值就是推移質運動速度（va） （Colin D，2002），即：

ADCP 所提供的測量速度有東西（ve）、南北（vn） 方向上的分量，利用公式計算出速度（Colin，2007）：

3 結果

3.1 各測點航跡圖

通過ADCP 數據處理軟件WinRiver 查看航跡圖，筆者選取了各站點觀測期間漲、落急時段的航跡圖（圖3），三角形代表測量船參考底跟蹤的航跡，正方形代表測量船參考GPS 定位的航跡。通過航跡圖，能夠觀察到，底跟蹤的航跡與GPS的航跡有明顯的偏移，且各測點落潮與漲潮時刻的偏移方向是相反的。S1 測點落急時段，GPS 航跡基本不變，符合定點測量的事實；底跟蹤航跡持續向西北方向延伸了一段距離，起點為（0，0），終點為（-256，230），航跡向西北方向移動了344 m，歷時為2 400 s，由此可算出船的平均相對運動速度為0.51 m/s，即底床相對于測量船的運動平均速度為0.51 m/s。由此說明底床存在推移質運動。同樣，S1 漲急時段和S2、S3、S4、S5、S6 漲、落急時段都存在推移質運動。

可見，在漲、落急時段流速較強、推移質運動速度較大時，ADCP 可以直觀地觀測到推移質運動速度變化，且通過觀察水流流向（表1） 和底跟蹤航跡偏移方向，可得出推移質運動方向和水流流向基本一致。

3.2 各測點潮流特征

S1、S2、S3、S4、S5、S6 測點處的潮流性質屬于非正規半日淺海潮流，水流運動形式為往復流，漲落潮流向相反（沈煥庭等，1979）。觀測數據顯示（表1），受長江徑流量及測點位置的影響，S1、S2、S3、S4 和S6 測點落潮平均流速一般比漲潮平均流速大；因S5 測點位于最大渾濁帶的攔門沙區域，潮流和徑流動力條件相當，故漲、落潮平均流速差別不大。S1、S2、S4 和S6 測點落潮最大流速與漲潮最大流速，分別相差0.56 m/s、0.12 m/s、0.29 m/s、0.16 m/s；S3 和S5 測點情況相反，漲潮最大流速比落潮分別大了0.11 m/s 和0.82 m/s。

表1 各測點落潮、漲潮流速對比單位：m/s

S1、S2、S3 和S4 的流速過程線存在著明顯的不對稱現象，落潮流速過程線相對平緩且持續時間較長；漲潮流速過程線的頂峰明顯且持續時間短，落潮流的強度超過漲潮流，落潮歷時長于漲潮歷時，強流速時間持續3～4 小時。漲、落潮流流向基本上與河槽主軸線平行。S5 和S6 測點漲、落潮流速過程線基本上對稱，漲落潮峰值和持續時間都相近。

3.3 各測點床沙粒徑

各測點床沙粒徑統計如表2，S1 測點的中值粒徑最大，D50=0.094 mm，為極細砂質。S6 測點的中值粒徑最小，D50=0.010 mm，為中黏土質。所以此次研究討論范圍在粗粉砂至極細砂質河床。

表2 各測點河床粒徑統計表

4 討論

4.1 推移質運動速度隨流速變化過程

圖3 各測點垂線平均流速（實線）、平均流向（虛線） 和推移質運動速度

通過公式（2） 和公式（3） 可計算出推移質運動速度，可以直觀地看到推移質運動速度在潮周期內的變化過程（圖3）。觀察推移質運動速度變化過程曲線和潮流流速變化過程曲線，我們可以發現兩者呈現“對應、不對稱”的特征。

“對應”指在觀測的二個潮周期內，各測點的推移質運動速度變化都與垂線平均流速變化一致，即在落急和漲急時段，推移質運動速度會達到最大值。但在S1（10 ∶25-15 ∶25）、S2（10 ∶25-15 ∶25） 測點的一次漲潮過程中和S3（7 ∶45-22 ∶45）測點的一次漲潮至落潮的變化過程中，推移質運動速度并沒有明顯的增大后又減小的變化過程，這與潮流流速和底質粒徑大小有關，過小的流速不能使該測點底質發生推移質運動（錢寧等，1983）。

“不對稱”現象指推移質運動速度變化歷時比潮流流速變化歷時短（表3），推移質運動速度變化過程曲線與潮流流速變化曲線比較，其峰值較尖。以泥沙起動流速為臨界劃分推移質運動狀態（胡浩等，2014），當流速達到臨界流速時，明顯的推移質運動才會發生，這與測點的底質粒徑大小有關（圖4）。

各測點底質粒徑與推移質運動起動臨界速度形成反比趨勢，受細顆粒泥沙之間的粘結力的影響，隨著粒徑的減小，也變得越來越不容易起動。

表3 各測點潮流、推移質運動速度變化歷時單位：h

圖4 各測點臨界流速和中值粒徑關系

4.2 推移質運動速度變化特征

表4 為各測點推移質運動速度的統計，S1、S2、S3、S4 和S6 測點落急的推移質運動速度均大于漲急的推移質運動速度，且落潮階段推移質運動持續時間大于漲潮階段（圖3）；S5 測點的推移質運動漲、落潮時段沒有明顯的差別，這與當地的潮流特征保持一致。

S1、S2、S3、S4、S5、S6 各測點受淺海、河口水下地形、徑流等影響，潮流變化均有顯著的潮汐不等特征，相鄰的二次落潮或（漲潮） 的流速不等，導致了各測點相鄰的二次落急（或漲急） 的推移質運動速度不等。這種不等在S1 和S3 測點體現得尤為明顯，S1 測點第一次落潮推移質運動速度比第二次落潮推移質運動速度大了一個數量級（表4）；在S3 測點，第一次漲落潮推移質運動速度明顯大于第二次落潮推移質運動速度，甚至在第二次落潮和漲潮之間的推移質運動速度始終保持在一個較小的值，沒有發生明顯起伏變化（圖3）；同樣，在S2 的兩次漲潮時段和S4、S6 的兩次落潮時段都出現了推移質運動速度日變化不相等的情況。在各測點，潮汐不等現象導致了推移質運動速度的日變化不等，又因為臨界起動速度的存在，使推移質運動速度的變化不等現象比潮汐日不等現象更加顯著，說明在較長時間尺度下，長江河口凈推移質輸沙方向與落潮流方向一致。

表4 各測點漲、落急時段推移質運動速度單位：m/s

4.3 誤差分析

基于ADCP 測量推移質運動的誤差主要由DGPS 系統誤差、船體晃動產生的傾斜誤差和ADCP 自身噪聲組成。

在潮流轉流時段，因為流速變小，甚至為零，所以實際推移質運動速度為零。但在圖3 中我們可以看到，在各測點轉流時段，推移質運動速度雖然小但不為零，這是因為GPS 存在誤差，無法與底跟蹤完全同步。可以認為在潮流轉流時段，流速條件不足以發生推移質運動，而此時推算出的推移質運動速度就是絕對誤差的下限，對實際推移質運動速度存在高估，通過計算高估值范圍在0.016-0.038 m/s 之間（表4）。可見，S1 落急1，S2 落急1、漲急1、落急2 和S6 落急2 計算得出的推移質運動速度比高估值大了一個量級，具有較高的可信度。

5 結語

（1） 通過ADCP 多普勒流速剖面儀觀測，在漲急、落急時段有明顯的推移質運動，能直觀體現半日潮對底沙運動的周期性影響。

（2） 推移質運動速度變化過程曲線和潮流流速變化過程曲線呈現“對應、不對稱”的特征。

（3） 推移質運動速度的日不等現象比潮流流速日不等現象更加顯著。

（4） 通過轉流時段推移質運動不存在的特征，計算出使用ADCP 觀測推移質運動速度的高估范圍為0.016-0.038 m/s。

C.D.Rennie,F.Rainville,S.Kashyap,2007.Improved estimation of ADCP Apparent bed-load velocity using a real-time Kalman filter.Journal of Hydraulic Engineering. (12) :1337-1344.

Colin D.Rennie,Robert G.Millarand, Michael A, 2002.Measurement of Bed Load Velocity using an Acoustic Doppler Current Profiler Church.Journal of Hydraulic Engineering,128:473-483.

Colin D. Rennie, Paul V. Villard, 2004. Site specificity of bed load measurement using an acoustic Doppler current profiler. Journal of Geophysical Research,109.

Kostaschuk R,Best J,Villard P,et al, 2005.Measuring flow velocity and sediment transport with an acoustic Doppler current profiler.Geomorphology.May 15,68 (1/2) :25-37.

胡浩，程和琴，楊忠勇，2014.基于ADCP 的長江口推移質運動遙測技術研究.人民長江，45（14）：6-9.

錢寧，萬兆惠，1983.泥沙運動力學，北京：科學出版社.

沈煥庭，潘定安，1979.長江河口潮流特性及其對河槽演變的影響.華東師范大學學報（自然科學版），01：131-144.

吳中，陳力平，游目林，等，2002.底部浮泥表層推移速度分布的ADCP-GPS 估測方法.海洋工程，20（4）：85-88.

楊俊輝，2009.ADCP、OBS 在底部泥沙運動觀測中的應用探討.港工技術，46（z1）：90-93.

張瑞瑾，謝鑒衡.河流泥沙動力學.北京：中國水利水電出版社，1998.

趙勝凱，王志芳，2007.ADCP 基本原理及應用.河北水利， （11）：25.