黃河口潮灘泥沙絮凝研究

于上，謝衛明*，何青，王憲業，趙中豪，郭磊城，徐凡

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

1 引言

絮凝是泥沙的重要特性，小顆粒泥沙會團聚成為大的絮凝體，進而改變沉降速度等特性[1–6]。在河口海岸區域，超過90%的泥沙都是以絮凝體的形式運輸[7]。研究泥沙絮凝特性及影響因素，對深化認識河口海岸泥沙輸運規律有重要意義。已有諸多研究對天然水體中的絮凝特性與影響因素進行了探討：陳錦山等[8]、程江等[9–10]和唐建華[11]利用LISST 觀測了長江口細顆粒泥沙的絮凝特性，發現長江口絮團平均粒徑為86.3 μm，受水動力影響明顯；Xia 等[12]、鄧智瑞等[13]、林建良等[14]對珠江口絮凝體特性也進行了研究，發現珠江口絮團粒徑略大于長江口，絮團粒徑受鹽度影響較小但受動力影響較大；長江中下游淡水環境中也發現有絮凝現象的存在，淡水絮團的粒徑略低于河口絮團[15–16]。綜合現場觀測與室內實驗發現，影響絮凝的因素包括水體紊動剪切強度、懸沙濃度、鹽度、泥沙粒徑、泥沙級配以及生物作用等[8,13,17–20]。

過往對絮凝的現場觀測研究主要集中在河流與河口區域，對潮灘區域絮凝特性的研究相對較少。但是另一方面，潮灘上存在明顯的絮凝現象，且絮凝對潮灘泥沙運移與沖淤變化有明顯影響[21–23]。因此，研究潮灘區域泥沙絮凝對進一步揭示河口海岸區域泥沙運動規律有重要意義。近年來，已經有一些研究者在美國、荷蘭的潮灘上進行了絮凝觀測研究，但是這些研究普遍集中于中潮（潮差在2～4 m 之間）或強潮河口（潮差大于4 m），鮮見在弱潮河口潮灘進行的泥沙絮凝觀測[21,23–24]。

基于以上背景，本文選取黃河口南部潮灘開展絮凝研究。黃河是著名的高含沙量河流，2019 年黃河口利津站實測輸沙量為2.71×108t，大量泥沙在河口輸移沉降。黃河口屬于弱潮河口，潮差在0.7～1.7 m 之間[25]。因此，黃河口區域適于開展泥沙絮凝研究。通過現場觀測獲取黃河口絮團粒徑、懸沙濃度、波浪、潮流動力等數據，可以研究弱潮河口波流共同作用下潮灘泥沙的絮凝規律與控制因素，豐富對我國河口海岸地區水沙運動規律的認識。

2 數據和方法

2.1 研究區域概況

研究區域選定在黃河口南部、萊州灣西岸的潮灘上，觀測點位于墾東大堤外側，屬于中國山東省東營市（圖1a）[26]。該區域潮灘坡度較小，退潮時有大片潮間帶露出水面。研究區域的波浪主要受季風控制，全海區的風浪以波浪為主[25]。潮汐類型屬于不正規混合半日潮，平均高潮位為1.6 m，平均低潮位為0.6 m（理論深度基準面）。該區域沉積物來源既包括黃河攜帶的陸相泥沙，也包括漲潮流攜帶的海相泥沙[27–28]。對采集觀測架附近表層沉積物于實驗室內進行粒度測量，測得中值粒徑為64.66 μm。

圖1 研究區域及儀器布置示意圖Fig. 1 Sketch of study area and instrument layout

2.2 觀測方法與過程

實施現場觀測的時間為2019 年11 月18–23 日。觀測點位于黃河口南部、萊州灣西岸潮灘上，距墾東大堤垂線距離約2 400 m，經緯度坐標為37°34′40″N，118°57′01″E。觀測儀器被固定于觀測架上，安置儀器前將觀測架錘入泥灘約1 m 以盡可能減小振動，避免對儀器觀測產生影響（圖1b）。

觀測架上儀器布置情況如下：ADV，測定觀測期間的水體紊動數據及灘面高程變化；OBS，測定濁度與水深數據；LISST-200（Laser In Situ Scattering and Transmissometry），采集現場懸浮體粒徑及體積濃度數 據；RBRsolo3 D|wave16（示 意 圖 中 簡 記 為RBR），采集波浪數據。以上儀器數量均為1 臺。OBS、RBRsolo3 D|wave16 及LISST-200 探 頭 距 床 面 高 程 均為5 cm；ADV 探頭距床面高程為15 cm，該儀器的水動力采樣點被設置為探頭下方10 cm，即測量距床面高程為5 cm 處的水體紊動數據。4 臺儀器時間同步，采樣間隔均為5 min。各儀器探頭位置接近，可以認為反映同一點的水沙特征。

安置及回收儀器時，均同步采集了觀測站附近的表層沉積物，用于分析床沙粒徑組成及對OBS 進行標定。研究時間內風速、風向收集自山東東營墾利氣象站（37°34.8′N ，118°33′E）實測數據。

2.3 數據處理與計算方法

2.3.1 懸沙濃度

觀測期間的懸沙濃度（SSC）根據OBS 測定的濁度（NTU）數據反演后獲得。在實驗室內使用現場獲取的表層沉積物配置不同濃度的懸沙濁液并測定對應的濁度，獲得濁度與含沙量間的相關關系[29]。標定結果顯示懸沙濃度與濁度間關系可以用關系式SSC=0.006 9 NTU–0.240 8 表 示， 相 關 系 數 的 平 方R2=0.89（數據點數量N=13）表示。

2.3.2 表層沉積物粒度

安置及回收儀器時各采集觀測站附近的表層沉積物樣品2 份，在實驗室中每份樣品測驗2 次，共測驗8 次。測驗時隨機從樣品袋中取出5 g 左右泥沙，加入分散劑六偏磷酸鈉之后超聲波震蕩30 min，使絮團充分破碎之后，使用Coulter LS-100Q 激光粒度儀測定[30–31]。計算8 次測驗所得的中值粒徑平均值，作為觀測站處表層沉積物中值粒徑，計算結果為64.66 μm。

2.3.3體積濃度及中值粒徑

本研究使用美國Sequoia 公司的現場激光粒度儀LISST-200來觀測天然水體中的絮團。射入水中的激光束會被懸浮顆粒所散射，顆粒物粒徑不同，對激光的散射角度也不同。大顆粒以小角度散射、小顆粒以大角度散射[1]。LISST-200 將向前散射的角度分為36 級，分別對應1～500 μm 之間的36 個粒級。運用Mie 散射理論，可以反推得到36 個粒級的懸浮體的體積濃度，并計算得到懸浮體中值粒徑[1,9]。綜合前人研究發現，潮灘環境中約80%的泥沙以絮團的形式輸運[1,9]。綜合郭超和何青[16]利用LISST-100 在荷蘭Scheldt 河口Kapellebank 潮灘、鄧智瑞等[13]在珠江口磨刀門海域、程江等[10]在長江口海域、Guo 和He[15]在長江流域的觀測實踐方法，以上在不同水體環境中開展的絮凝研究將LISST 觀測到的懸浮體視為絮團。因此在本研究中，也將LISST-200 的觀測數據整體上視為絮團進行處理。

2.3.4有效密度

有效密度是反映絮團結構的重要參數，且對沉降速度有明顯影響[4,17]。本研究采用Fettweis 方法計算有效密度[3]，公式為

式中，ρw與ρs分別為水與組成絮團的分散態泥沙顆粒密度（單位：kg/m3），分別取為1 030 kg/m3與2 650 kg/m3；SSC 為懸沙濃度（單位：kg/m3）；TVC 為顆粒總體積濃度，即36 個粒級組的絮團體積濃度之和（單位：μL/L）。

2.3.5沉速

沉速可以認為是絮團重力W與下沉過程中受到的阻力F達到平衡時等速沉降的值，采用Stokes 公式計算[5]：

式中，Dm為絮團平均粒徑（單位：μm）；μ為水的動力黏滯系數，根據水溫選取相應的值（單位：Pa·s）；g為重力加速度，取為9.8 m/s2。

2.3.6紊動剪切率

在泥沙絮凝研究中，通常使用紊動剪切率G（s-1）來描述水體的紊動強度，其計算方法為[17,32]

式中，u*為摩阻流速（單位：m/s-1）；z和H分別為距離床面高程與水深（單位：m）；ν為水體運動黏度，根據水溫選取相應的值（單位：m2/s）；κ為卡門常數，取為0.4。水流與波浪共同引起的床面切應力τcw（單位：N/m2）為[21,33]

式中，τc為水流引起的切應力，使用改進的湍流動能法（Turbulent Kinetic Energy，TKE）計算[34]：

τw為風浪引起的底部剪切力，計算方法為[35–36]

式中Aδ為波質點軌道流速（單位：m/s）；Hs為有效波高；k為波數；Ts為波浪有效周期（單位：s）；L為波長（單位：m）；fw為波浪摩擦系數；Rew與r分別為波浪雷諾數（無量綱）與相對糙率（無量綱），分別由下式求得[35–36]：

式中，ks為床面粗糙高度，取2.5d，d為潮灘表層沉積物中值粒徑。根據實驗室測驗結果，觀測期間天氣較為平靜，潮灘表層沉積物粒度沒有明顯變化。d取觀測前、后的表層沉積物中值粒徑平均值，為64.66 μm。根據式（4）至式（13）計算得到床面剪切力τcw之后，可由式（14）求得摩阻流速，并代入式（3）得到紊動剪切率[5]。

2.3.7臨界侵蝕切應力及臨界沉降剪切應力

根據前人在黃河口潮灘的研究，李華國等[37]提出的模型能夠較好地預測潮灘表面沉積物的臨界剪切應力：

式中，c為黏結力系數，取2.9×10-5kg/m；γd、γdc分別為表層沉積物的干容重及穩定干容重（單位：kg/m3）；其余符號含義同前[37]。γd的值通過“環刀法”采樣之后在實驗室中測量，為1 490 kg/m3；穩定干容重γdc指經過長期的一般壓密后的沉積物干容重，根據張耀哲和王敬昌[38]提出的經驗公式計算為

式中，γs為泥沙顆粒干容重（單位：kg/m3），取2 650 kg/m3；其余符號含義同前。計算得γdc的值為1 373.85 kg/m3，計算結果接近前人的實測值，結果可信[38–39]。根據式（15）、式（16），計算得觀測架附近表層沉積物臨界侵蝕切應力τe為0.086 N/m2。

實驗室分析結果顯示，研究區表層沉積物成分為黏土（8.13%）、粉砂（46.95%）和砂（44.92%）。根據Lumborg[40]實驗結果，研究區表層沉積物懸沙臨界沉降剪切應力τd取0.05 N/m2。也有研究表明，τd約為臨界侵蝕切應力τe的0.5 倍[41]。兩種方法得到的值接近，證明結果可信。

2.3.8相關性檢驗

本文選用適用性較強的Spearman 秩相關系數ρ來檢測水文變量之間的相關關系，這里的“秩”指的是順序或者排序。所以，計算時先對由n個樣本組成的時間序列x1，x2，…，xn，和由n個樣本組成的時間序列y1，y2，…，yn進行排序（同時為升序或者降序），得到排序后的集合X和Y，Xi、Yi分別是xi、yi在X、Y中的排行，i=1，2，3…，n，則ρ的計算公式為

計算所得的相關系數ρ值，一般取小數點后兩位來表示，取值范圍在–1.00～1.00 之間，負數表示負相關，正數表示正相關。ρ值對應的兩變量之間的相關程度為：

ρ= ±0～±0.30 時，不相關；

ρ= ±0.30～±0.50 時，微相關；

ρ= ±0.50～±0.80 時，顯著相關；

ρ= ±0.80～±1.00 時，高度相關。

3 水沙動力特性與變化過程

3.1 潮流變化過程

研究區域潮型為不正規混合半日潮，只有高潮情況下才會淹沒觀測架，觀測時段內觀測架一共被淹沒4 次，漲落潮歷時大致相等。11 月19–22 日，分別記作T1–T4，最大水深分別為0.19 m、0.62 m、0.53 m與0.72 m（圖2a）。T2–T4 潮位接近，且均顯著大于T1。觀測期間的潮流呈現明顯的旋轉流特性。漲潮階段由約N100°順時針旋轉至約N50°，落潮階段由約N50°順時針旋轉至約N200°（圖2b）。

研究區域流速存在兩種不同的隨時間變化模式。T1、T2 流速呈現“單峰型”，潮周期內只出現一次流速極值，漲潮平均流速約為落潮平均流速的4 倍，最大流速出現在漲潮中期，分別為0.52 m/s 與0.46 m/s。T3、T4 流速呈現“雙峰型”，漲潮、落潮各出現一次流速極值，漲落潮平均流速接近，漲落潮中期均存在明顯的流速峰值。T3 漲潮最大流速為0.49 m/s，落潮最大流速為0.50 m/s；T4 漲潮最大流速為0.60 m/s，落潮最大流速為0.53 m/s（圖2c）。

3.2 風浪變化過程

T1 期間波浪動力較弱，有效波高Hs低于0.05 m。T2–T4 期間波浪動力較強，有效波高Hs最大值分別為0.24 m、0.21 m 及0.27 m。T2–T4 期間，漲落潮初期流速、流向波動較大，這是因為波浪動力較強，在水深較淺的情況下，波浪對流速流向的影響較顯著（圖2d）。T1 期間有效波浪周期Ts在1.24～3.21 s 之間，平均為2.32 s。T2–T4 期間有效波浪周期Ts在1.16～2.11 s 之間，平均為1.68 s，高于T1 期間。

T1 期間平均風速為1.33 m/s，最低為0.5 m/s，最高為2.1 m/s，風向由向岸（328°）轉變為離岸（73°）。T2–T4 期間平均風速為2.67 m/s，最高為4.60 m/s，最低為0.10 m/s，風向為離岸風，平均風向為153.94°（圖2e）。

圖2 觀測期間水動力、波浪與風速隨時間變化序列Fig. 2 Time series of hydrodynamic force, wave and wind during the observation period

3.3 床面剪切力及紊動剪切率變化過程

T1 期間，波流共同引起的床面剪切力（τcw）平均為0.09 N/m2，τcw最大達到0.34 N/m2，中潮期間τcw平均為0.23 N/m2，最大可以達到0.99 N/m2。T1 期間潮流引起的床面切應力τc平均為0.05 N/m2，最大達到0.21 N/m2；T2–T4 期間τc平均為0.14 N/m2，最大達到0.68 N/m2。觀測期的4 d 中，總床面剪切力τcw與潮流引起的床面剪切力τc的變化趨勢接近，峰值都出現在落潮中期（圖3a）。

黃河口為弱潮河口，研究區屬于開敞性海岸，波浪沒有遮蔽，波浪引起的床面剪切力τw較高[42]。T1 期間τw在0.90～1.19 N/m2之間，平均為0.98 N/m2；T2–T4 期間τw在0.99～1.30 N/m2之間，平均為1.11 N/m2，T2–T4 期間τw高于T1。波浪引起的床面剪切力變化圖線呈“U”形，漲潮初期與落潮后期τw較高，水深較大時的τw反而較小。說明水深較淺時，波浪對床面切應力的影響會明顯增加（圖3a）。

紊動剪切率G的變化規律與床面切應力大致相同，在落潮中期出現峰值。T1 期間紊動剪切率在1.88～9.63 s-1之間，平均為3.33 s-1；T2–T4 期間紊動剪切率在2.51～28.48 s-1之間，平均為9.00 s-1，T2–T4期間水體紊動剪切率大于T1（圖3b）。

3.4 懸沙濃度變化過程

研究區域水體懸沙濃度較低，T1 期間水體平均懸沙濃度為0.08 kg/m3，最大可以達到0.26 kg/m3。T2–T4 期間，落潮中后期水體懸沙濃度會明顯增加，出現一個明顯的懸沙濃度峰。T3–T4 期間水體平均懸沙濃度為0.10 kg/m3，最大可以達到0.40 kg/m3。懸沙濃度峰出現的時間與床面切應力峰值出現的時間接近，且出現懸沙濃度峰時，床面切應力顯著高于臨界剪切應力τe（0.086 N/m2），說明這一懸沙濃度峰的出現源于泥沙再懸浮（圖3c）。

圖3 觀測期間床面切應力與懸沙濃度（SSC）隨時間變化序列Fig. 3 Time series of bed shear stress, bottom turbulent shear rate and suspended sediment concentration (SSC) during the observation period

T2 期間水體懸沙濃度顯著高于T3 和T4，水體平均懸沙濃度為0.38 kg/m3，最大可以達到0.62 kg/m3。T2 期間，懸沙濃度隨漲潮不斷增加，在最高潮位前后保持峰值，后隨水深降低而不斷減小，在落潮中后期出現第2 個懸沙濃度峰。T2 期間兩個懸沙濃度峰值出現的時間與床面切應力峰值出現的時間不一致，懸沙濃度變化過程與床面切應力變化過程不同步。但是T2 期間懸沙濃度增加的階段，與流速峰值出現的時間基本吻合。推測T2 期間的懸沙濃度峰值并不是來源于泥沙再懸浮?？赡苁怯捎跐q潮水流從附近攜帶大量泥沙，并在高水位期間沉降造成的[43]（圖3c）。

觀測中發現，研究區潮灘存在大量的“泥丸”型結構，即大型的泥沙聚集體（圖4）。這種聚集體呈球狀，可以隨潮汐運動在灘面遷移滾動，但又容易破壞，一旦分散會導致所在局部懸沙濃度大幅上升。T2 期間出現懸沙濃度峰值，可能是發生了“泥丸”的遷移與分散。

圖4 研究區域的“泥丸”型泥沙聚合體Fig. 4 The “pill-like” sediment aggregates in the study area

4 絮凝特性與變化

4.1 絮團粒徑

T1 期間絮團粒徑先增大后減小，在最高潮位階段達到最大，為264.44 μm。T2–T4 期間絮團粒徑呈現緩慢遞減趨勢，中潮超過60%的時間內，絮團中值粒徑都低于100 μm。但是在T2–T4 每個落潮的中后期絮團中值粒徑會出現一個峰值，最大可達到203.83 μm。小潮期間絮團平均粒徑為123.90 μm，大于中潮期間的91.81 μm（圖5a）。

4.2 絮團體積濃度

T1 期間絮團總體積濃度平均為221.05 μL/L，最大為810.22 μL/L；T2–T4 期間絮團總體積濃度平均為390.77 μL/L，最大為1 100.29 μL/L。在T1 最高潮位及T2、T3 的落潮中期，絮團總體積濃度均出現一個峰值（圖5b）。

圖5 觀測期間絮團粒徑、體積濃度、有效密度與沉速隨時間變化序列Fig. 5 Time series of floc size, volume concentration, effective density and settling velocity during the observation period

4.3 絮團有效密度與沉速

T1 期間黃河口潮灘絮團有效密度在56.55～452.44 kg/m3之間，平均為257.03 kg/m3。T2–T4 期間絮團有效密度逐日遞減，T2 期間平均有效密度為427.77 kg/m3，T3 期 間 絮 團 平 均 密 度 為319.38 kg/m3，T4 期間平均有效密度為166.32 kg/m3（圖5c）。T1 期間，絮團沉速在0.99～1.86 mm/s 之間，平均為1.41 mm/s。T2–T4 期間，絮團沉速在0.17～4.07 mm/s之間，平均為0.94 mm/s?？傮w來看，中潮期T2–T4 期間絮團的沉速大于小潮期T1（圖5d）。

5 分析與討論

5.1 絮團粒徑變化的影響因素

從泥沙運動學的角度，天然水體中的絮凝主要受到懸浮泥沙粒徑、水動力環境、懸沙濃度的控制[17]。本研究期間受交通與天氣條件限制，未采集懸沙樣品，無法分析懸浮泥沙粒徑對絮凝的影響，但可分析水動力與懸沙濃度與絮凝的相關性。采用相關性分析方法，分別計算各個潮周期內水體紊動剪切、流速及懸沙濃度與絮團粒徑之間的Spearman 系數，計算結果及對應的相關性如表1。絮團粒徑與流速的關系主要為不相關，與懸沙濃度及底部紊動剪切率則存在著較明顯的相關關系，說明黃河口潮灘絮團粒徑主要受到底部紊動剪切率與懸沙濃度的影響，與流速關系不大。

表1 各個潮周期絮團粒徑與紊動剪切率、懸沙濃度及流速的Spearman 系數Table 1 Spearman correlation coefficient of floc particle size and bottom turbulent shear rate, suspended sediment concentration and flow velocity

5.2 水體紊動剪切促進絮凝的上限

研究發現，不同潮周期中，黃河口潮灘絮團粒徑與紊動剪切率的關系不盡相同。T1 期間，絮團粒徑與紊動剪切率變化趨勢基本一致，粒徑隨紊動剪切率增大而增大，紊動剪切率減小時粒徑也隨之減小。但是T2–T4 各個潮周期內，水體紊動剪切率與絮團粒徑呈現負相關趨勢，紊動強度較大時對應粒徑較小。

在影響絮凝的外部因素中，水體紊動是最重要的因素之一，水體紊動強度的增大可以增加顆粒的碰撞頻率，促進顆粒聚合；但當紊動強度增大到一定程度時，反而會增加絮團破碎頻率抑制絮凝[2,44]。T1 期間為小潮，水動力較弱，底部紊動剪切率平均僅為3.33 s–1。T2–T4 期間底部紊動剪切率較高，平均為9.00 s–1，最高可達到28.48 s–1。在T1 期間紊動剪切較弱的情況下，水體紊動增加對絮凝起促進作用；T2–T4 期間紊動強度較高，水體紊動增加抑制絮凝。

為了進一步探究黃河口潮灘水體紊動影響絮凝的上限，首先對實測數據進行等距分箱以體現規律性，分箱間距為4[22–23]。發現絮團粒徑隨紊動強度增加先增后減。選擇不同的臨界紊動剪切率值（Gl），將觀測結果分為“絮團粒徑隨紊動強度增加而增大”與“絮團粒徑隨紊動強度增加而減小”階段，并對兩個階段分別使用冪函數擬合（圖6）。

圖6 絮團粒徑與底部紊動剪切率關系Fig. 6 Variation of floc size with bottom turbulent shear rate

結果發現，當Gl=3.76 s-1時，取得的擬合結果最為理想（表2）。T1 期間，水體紊動剪切率低于3.76 s-1，T2–T4 期間則高于這一值，統計分析得到的結果與實測數據能夠相互驗證。因此可以認為，黃河口當地水體紊動促進絮凝的上限為Gl=3.76 s-1，水體紊動強度低于這一數值時水體紊動促進絮凝，高于這一數值時水體紊動破壞絮凝。

表2 不同Gl 對應擬合結果對比Table 2 Fitting results of different Gl

5.3 懸沙濃度對絮凝影響

懸沙濃度對絮凝的影響尚未有確切定論，一些實驗發現，懸沙濃度增加增大了顆粒碰撞概率，對絮凝起促進作用[2,20,45]；但也有研究認為，懸沙濃度增加導致絮團碰撞更頻繁，造成絮團破壞[17,20,45–46]。參考郭超[17]的研究方法，對黃河口潮灘絮團粒徑與懸沙濃度之間的關系整體采用冪函數擬合，發現二者呈負相關關系，相關系數為R2=0.778 8。說明在黃河口潮灘，懸沙濃度增加對絮凝起抑制作用（圖7）。

圖7 絮團粒徑與懸沙濃度（SSC）關系Fig. 7 Variation of floc size with suspended sediment concentration

在5.1 節的討論中，發現T2–T4 期間懸沙濃度與絮團粒徑呈負相關，但是T1 期間懸沙濃度與絮團粒徑沒有明顯相關性。這是因為絮凝同時受到懸沙濃度與水動力共同控制，T1 期間水體紊動強度先增后減，懸沙含量則緩慢降低。絮團粒徑的變化趨勢與水體紊動強度類似，也是先增后減。這說明相比于懸沙含量變化造成的影響，水體紊動剪切對絮凝的影響更大。懸沙含量降低對絮凝的促進作用，不足以抵償水體紊動強度變化對絮凝的影響，水體紊動剪切才是絮團粒徑變化的主導因素。因此在T1 期間，絮團粒徑與水體紊動剪切強度顯著相關，但是與懸沙濃度幾乎沒有相關性。

5.4 絮團有效密度與沉速變化規律

黃河口潮灘絮團粒徑與有效密度呈負相關關系，粒徑更大的絮團對應的有效密度更小。過往研究中，通常使用 Δp=·Dam的關系來描述絮團粒徑與有效密度的關系，a值在–0.46（Sternberg 等[6]在美國加利福尼亞海域測得）到–1.58（Xia 等[12]在中國珠江口測得）之間[5,10,19]。郭超[17]、Guo 等[21]在荷蘭Kapellebank 潮灘常態天氣下測得a=–0.53，風暴天氣下測得a=–1.28。

分別對黃河口潮灘小潮（T1）與中潮（T2–T4）期間的絮凝特性進行研究，計算得小潮期間a=–1.52，中潮期間a=–0.79，屬于上述文獻中所列的取值范圍（圖8）。無論在何種動力條件下，有效密度都與絮團粒徑成反比，粒徑越大的絮團對應的有效密度越低。在水體紊動剪切促進絮凝的情況下，泥沙聚合體在由小變大的過程中，先由易絮凝的細顆粒泥沙發生絮凝，而后捕獲粒度相對更大的泥沙，形成的絮凝體越高階，其孔隙率越大，有效密度越低[10]。在水體紊動剪切抑制絮凝的情況下，結構疏松、有效密度低的大絮團更容易受到紊動影響破碎，形成粒徑更小但是有效密度更高的小泥沙聚合體[10,17,19]。

圖8 絮團粒徑與有效密度關系Fig. 8 Variation of floc size with effective density

由于絮凝受到泥沙特性、當地生物化學特性等多方面影響，a值與水動力強弱的關系目前尚無定論。Guo 等[21]在荷蘭Kapellebank 潮灘的研究中，水動力較強的風暴天氣下計算得的a值低于水動力相對弱的常態天氣；但是本研究與其他觀測結果中，也發現了相反的結論[9,16]。本研究在黃河口觀測的結果顯示，小潮期間a值低于中潮；程江等[9]在長江口徐六涇的觀測結果也類似，小潮期間的a值（–0.646）低于大潮期間（–0.368）。

根據式（2），絮團的沉速與絮團粒徑及有效密度有關。過往研究發現，沉速與粒徑之間的關系與“有效密度-粒徑”關系類似，可以用 ωs=c·Ddm的形式來表示。c取值在0.290～1.632 之間；d取值在0.000 022 8～0.008 20 之 間[5–6,10,12]。黃 河 口 潮 灘 絮 團沉速與粒徑之間相關性較好，沉速與粒徑間關系可以用ωs=0.0028·D1m.2679表示，屬于過往文獻中所列取值范圍（圖9）。但是對有效密度與沉速關系進行分析時發現，二者之間不存在明顯相關性。說明黃河口潮灘絮團沉速主要受到絮團粒徑的影響，與有效密度的關系不明顯。

圖9 絮團粒徑與沉速關系Fig. 9 Variation of floc size with settling velocity

6 結論

本文通過研究黃河口南部潮灘絮團特征的變化，得到如下幾點認識：

（1）黃河口潮灘絮團粒徑在25.42～264.44 μm 之間，平均為95.20 μm。小潮潮周期T1 期間絮團總體積濃度平均為221.05 μL/L，平均有效密度為257.03 kg/m3，平均沉速為1.41 mm/s；中潮潮周期T2–T4期間絮團總體積濃度平均為390.77μL/L，平均有效密度為304.49 kg/m3，平均沉速為0.94 mm/s。

（2）水體紊動較低時對黃河口潮灘絮凝起促進作用，較高時則抑制絮凝。隨紊動強度增加，黃河口潮灘絮團粒徑先增后減。分析得底部紊動剪切率促進絮凝的上限值為Gl=3.76 s-1。

（3）懸沙濃度增加對黃河口潮灘泥沙絮凝起抑制作用，但是懸沙濃度變化對絮凝的影響低于水體紊動剪切變化對絮凝的影響。

（4）黃河口潮灘絮團有效密度與粒徑呈現Δp=b·Dam的關系。T1 潮周期a值為–1.52，T2–T4 潮周期a值為–0.79。

（5）黃河口潮灘絮團沉速主要受到泥沙粒徑影響。絮團沉速與粒徑之間相關性較好，但與有效密度之間不存在明顯相關關系。

致謝：本文野外觀測工作得到水利部黃河水利委員會黃河水利科學研究院黃河河口研究所、山東黃河三角洲國家級自然保護區等單位的大力支持，在此表示誠摯的感謝！