黃驊港綜合港區單、雙堤延伸防沙堤方案減淤效果分析

張文杰

（滄州港務集團有限公司，滄州061113）

黃驊港綜合港區單、雙堤延伸防沙堤方案減淤效果分析

張文杰

（滄州港務集團有限公司，滄州061113）

在分析工程海區水沙條件基礎上，通過數學模型試驗的方法對黃驊港綜合港區單、雙堤延伸防沙堤方案減淤效果進行研究，論證了單堤延伸防沙堤的合理性、可行性，單獨延伸綜合港區北防沙堤至-8 m等深線位置能夠起到較好的減淤效果，先期單堤延伸防沙堤不失為一種好的起步方式。

數學模型試驗；單堤；雙堤；減淤效果

黃驊海域為中國典型的粉砂質海岸，泥沙活躍，航道回淤嚴重。根據神華煤碳港區的經驗，目前綜合港區航道已達20萬噸級，且正考慮雙向航道建設的可能性。隨著航道規模的持續增大，其泥沙年回淤量及維護疏浚量也隨之增大，每年的疏浚土方都在2 000萬m3以上，對港口運營來說負擔較重，需要采取防淤減淤措施。以往相關研究成果［文獻3-6］表明，對于20萬噸級航道而言，防沙堤堤頭宜設于-8 m水深附近是較為適宜的，延堤長度8.8 km，堤頂高程由-1 m漸變至-3 m（當地理論最低潮面，下同），本研究以此為基本前提，對黃驊港綜合港區單、雙堤延伸方案減淤效果采用數學模型試驗的方法進行論證。

1 黃驊海域基本水動力條件

本海區潮汐屬不正規半日潮，平均高潮位為3.58 m，平均低潮位為1.28 m，平均潮差為2.30 m。本海區潮流屬規則半日潮型，從總的平面分布趨勢看，外海流速大于近岸流速，漲潮流速大于落潮流速。漲潮潮段平均流速在0.29～0.42 m/s之間，流向為240o～300o；落潮潮段平均流速在0.25～0.37 m/s之間，流向為46o～97o。

本海區的波浪是以風浪為主，涌浪為輔。純風浪頻率為66.81%，涌浪為主的混合浪頻率為27.1%，風浪為主的混合浪頻率為4.64%，風涌混合浪頻率為0.12%。全年常浪向為E向，頻率為10.06%，次常浪向為ENE向，頻率為9.38%，強浪向為NE～E向。ENE為平均最大波高向，該向累年平均波高為0.97 m。

2 黃驊海域泥沙環境

黃驊海域為典型的粉砂質海岸，-10 m以內泥沙平均中值粒徑在0.03 mm左右，粘土含量在20%以內，泥沙活躍性很強。

黃驊港海域泥沙運動主要受風浪控制，風浪越大含沙量越大。在含沙量平面分布上，近岸含沙量明顯大于遠岸，同時近岸泥沙在落潮流帶動下向外海輸移，在口門處形成懸沙濃度較高區域，對航道局部淤積起著重要作用。從實測值來看：正常天氣下，-6 m、-10 m、-14 m、-16 m、-19 m水深處，漲潮潮段平均含沙量分別為0.107 kg/m3、0.110 kg/m3、0.098 kg/m3、0.065 kg/m3、0.064 kg/m3；落潮潮段平均含沙量分別為0.165 kg/m3、0.157 kg/m3、0.120 kg/m3、0.091 kg/m3、0.108 kg/m3；有明顯的隨水深增加而減小的趨勢，并且落潮含沙量大于漲潮含沙量。大風浪作用下，海區含沙量隨風浪變化而變化，對黃驊港口門附近進行底部含沙量觀測時發現，灘面以上0.5 m水層處曾觀測到10 kg/m3以上含沙濃度。含沙量的變化與波浪的變化過程有密切關系，當波高達到1.0 m以上后底層含沙量的增長迅速；當波高達到最大時，含沙量并未達到峰值，含沙量峰值置后于波高峰值約3 h；當波高開始衰減以后，底層含沙量衰減不明顯，仍維持一段時間。

造成黃驊港外航道泥沙淤積物的來源主要是灘面泥沙運動。波浪是本海區泥沙起動的主要動力，潮流為灘面物質的運動提供了條件。大風浪天氣在破波帶可形成高含沙區，高濃度泥沙隨潮流運動時在穿越航道的過程中，泥沙落淤，對外航道局部區域回淤帶來較大影響。

3 數學模型的建立和驗證

本模型風浪場模型采用SWAN模型，其基本控制方程可參考SWAN用戶手冊。水動力模型合理考慮了波浪對水流的影響，這些影響主要包括波浪輻射應力、波浪破碎引起的表面水滾應力、波浪存在引起水流底部摩阻的變化（底部剪切應力隨著底摩阻變化而變化）等。

3.1 水動力模型

式中：x,y,z為笛卡爾坐標系下的三維坐標，分別為東西方向、南北方向以及垂向的坐標；u,v分別為水平方向的東分量速度和北分量速度；w為垂向速度；g為重力加速度；t為時間；ρ0為平均密度；f為科氏參數；Km為垂向渦粘系數，Fu和Fv為水平動量擴散系數。

3.2 泥沙運動模型

（1）控制方程

非粘性泥沙運動以懸移質為主，控制方程為對流-擴散方程

式中：AH和Kh分別為水平和垂向泥沙質量擴散系數，AH由Mellor-Yamada 2.5階模式［2］計算得到，Kh采用van Rijn的泥沙摻混系數公式計算得到；u,v,w分別為x、y、z向水流速度；C為懸浮泥沙濃度；ws為泥沙沉速。

由于黃驊港區域泥沙運動同時受波浪和潮流兩種動力因素的影響，本研究在對原有模型進行了改進，即采用van Rijn的模式計算波流共同作用下懸沙垂向擴散系數，將波流共同作用下垂向擴散系數表示為潮流和波浪單獨作用時擴散系數的非線性疊加。

（2）邊界條件

自由水面要求含沙量的凈通量為零，即在z=ζ時，

底部邊界條件表示為

式中：Ei和Di分別為淤積率和沖刷率。當底部剪應力超過臨界沖刷應力時，則產生沖刷。沖刷率表示為

式中：E是沖刷率；E0是河床侵蝕強度；Fi是第i組份泥沙在所有組份中所占的百分比；單組份時，Fi取1；Pb是表層泥沙的孔隙率；τb是床面剪切應力；τe是臨界沖刷應力。

泥沙的沉積作用由以下方程控制

計算沉積通量和沉積物濃度更新需要滿足SLIP限制。它具有二階計算精度，對濃度的局部劇烈變化能夠比較準確的進行計算。

3.3 模型網格剖分

為準確掌握工程所在海域水動力特性，波浪和潮流均采用二重嵌套的方法進行計算。大模型計算范圍北起40°55′N，南至37°07′N，東至122°29′E，采用非結構化網格。小模型范圍北起39°12′N，南至38°00′N，東至119°00′E，采用局部加密的非結構化網格，網格空間步長最大為4 500 m，最小為30 m。大、小模型計算網格如圖1所示。

圖1 模型計算網格Fig.1 Computational scope and grid

3.4 數學模型的驗證

圖2為各測站位置示意圖。

圖2 各測站位置示意圖Fig.2 Sketch of measuring location

圖3 2003年10月溫帶風暴過程風速、風向驗證結果Fig.3 Verification of wind magnitude and direction in 2003-10

圖4 2009年11月溫帶風暴過程風速、風向驗證結果Fig.4 Verification of wind magnitude and direction in 2009-11

（1）風場模型的驗證。

本研究通過2003年10月份和2009年11月份溫帶風暴對風場模型進行驗證，驗證結果見圖3、圖4。可見本風場模型能很好的反應大風的的發展、成長、衰減過程，可用于為波浪模擬提供風場動力條件。

（2）波浪模型的驗證。

本研究通過2003年10月和2009年11月溫帶風暴實測波浪過程進行驗證，驗證結果見圖5。由圖可見，本風浪模型能很好的反應工程海域波浪的發展、成長、衰減過程，可用于本研究波浪模擬，并為后續泥沙模型提供動力條件。

（3）潮流模型的驗證。

本研究采用黃驊港海域2003年10月10日～14日風暴潮觀測潮位資料和2007年8月29日～30日大潮水文資料對潮流模型進行驗證，限于篇幅只給出部分測點數據，驗證結果見圖6、圖7。由圖可知，計算結果與實測符合良好，可為后續泥沙模型提供動力場。

（4）泥沙運動模型的驗證。

圖5 實測與計算有效波高、周期的比較Fig.5 Verification of wave height and period

圖6 潮位驗證結果Fig.6 Verification of tidal level

圖7 2007年8月29～30日實測與計算平均流速、流向比較Fig.7 Verification of flow magnitude and direction in 2007-08-29～2007-08-30

采用2007年8月29日～30日水文全潮含沙量資料進行驗證，驗證結果見圖8。由驗證結果可見，本模型能夠良好的闡述工程海域泥沙運動情況，可用于工程海域泥沙運動模擬。

（5）大風驟淤強度的驗證。

采用2003年10月10日和2009年11月黃驊港航道實測回淤驗證，驗證結果符合良好，見圖9，本泥沙回淤模型可用于本航道回淤計算。

圖8 2007年8月29日～30日垂向平均含沙量驗證結果Fig.8 Verification of sediment concentration in 2007-08-29～2007-08-30

圖9 溫帶風暴潮過程煤港航道回淤計算與實測結果Fig.9 Verification of channel deposition

4 論證方案

4.1 計算方案

綜合港區防沙堤分雙堤延伸或單堤（北堤）延伸兩種方案，延伸至－8 m（從當地理論最低潮面起算），防沙堤以潛堤形式平行灘面向外海延伸，潛堤高度由-1 m漸變至-3 m，進出港航道按20萬t級考慮（航道底寬250 m、-19 m深）。綜合港建港時航道設計水深以10 a一遇驟淤不礙航作為設計標準，本研究延續此標準，以10 a一遇大風驟淤作為單雙堤減淤效果對比條件。計算方案見圖10。

圖10 計算方案布置示意圖Fig.10 Sketch of calculation scheme

4.2 計算結果及分析

（1）流場情況。

計算結果表明：單堤（北堤）或雙堤延伸后，港池及航道內水流平順，未出現明顯回流情況，水流延堤向沒有明顯增大趨勢。防沙堤延伸前后，漲、落急時刻平均流速大部分觀測點基本相同，僅新、老口門處有所不同。

圖11 單堤延伸方案流速矢量圖Fig.11 Flow vector diagram of single dike project

圖12 延堤前后航道內橫流分布Fig.12 Cross flow distribution of each program

圖13 單、雙堤延伸方案10 a一遇大風淤積比較Fig.13 Comparison of two extension schemes

原口門附近Z4點最大橫流由延伸潛堤前的0.45 m/s減小到單堤延伸方案的0.39 m/s和雙堤延伸方案的0.35 m/s；新口門附近Z22點最大橫流由延伸潛堤前的0.24 m/s增大到單堤延伸方案的0.34 m/s和雙堤延伸方案的0.37 m/s，呈現新口門外側附近橫流有所增大趨勢；新掩護段內其它觀測點最大橫流相差不多。圖11為單堤延伸方案流速矢量圖，圖12為延堤前后航道內橫流分布。

（2）泥沙回淤情況。

10年一遇大風作用下，綜合港區防沙堤雙堤延伸至-8 m情況時，20萬t級航道內最大淤積強度1.04 m；全航道淤積總量為585萬m3，其中掩護段內為228萬m3，掩護段外357萬m3；綜合港區北防沙堤單堤伸至-8 m情況下，20萬t級航道內最大淤積強度1.05 m，全航道淤積總量為621萬m3，其中掩護段內266萬m3，掩護段外355萬m3。大風淤積分布見圖13。

對比雙堤掩護及北側單堤掩護效果，掩護段內延伸單堤減淤率能到達延伸雙減淤率的70%以上，可見在已形成航道相互掩護的基礎上，單獨延伸北堤單堤掩護效果良好，可作為20萬t級航道減淤整治工程的起步。

若神華南側防沙堤也延伸至-8 m，與綜合港區北堤一起對該港兩航道形成掩護，其減淤效果會更好。

5 結語

（1）單獨延伸綜合港區北防沙堤至-8 m等深線位置能夠起到較好的減淤效果。考慮到-6 m水深以往建設防沙堤的高額造價，先期單堤延伸防沙堤不失為一種好的起步方式。

（2）在綜合港區北堤延伸情況下，若神華煤碳港區南防沙堤也延伸至-8 m，其減淤效果應會更佳。

［1］JTS145-2015，《港口與航道水文規范》.

［2］TJS/T231-2-2010，《海岸與河口潮流泥沙模擬技術規程》.

［3］劉國亭.黃驊港海岸泥沙特性及防沙減淤工程措施的研究［D］.天津：天津大學，2005.

［4］季則舟.粉砂質海岸港口水域平面布局特點［J］.海洋工程，2006，24（4）81-85. JI Z Z.The layout arrangement features for water areas of port on silty beach［J］.The Ocean Engineering，2006，24（4）81-85.

［5］曹祖德，楊樹森，楊華.粉砂質海岸的界定及其泥沙運動特點［J］.水運工程，2009（1）：108-111. CAO Z D，YANG S S，YANG H.Definition of silt?sandy beach and its characteristics of sediment movement［J］.Port&Waterway Engineering，2009（1）：108-111.

［6］嚴冰，劉濤.黃驊港綜合港區單側延伸防沙堤潮流泥沙數模研究報告［R］.天津：交通運輸部天津水運工程科學研究所，2012.

［7］侯志強，劉濤.黃驊港綜合港區10萬噸級航道工程波浪潮流泥沙數學模型試驗研究報告［R］.天津：交通運輸部天津水運工程科學研究所，2010.

［8］侯志強，劉濤，張書莊.黃驊港綜合港區20萬噸級航道工程防沙堤減淤效果分析研究報告［R］.天津：交通運輸部天津水運工程科學研究所，2009.

［9］侯志強.黃驊港綜合港區港池、航道泥沙淤積數學模型試驗研究報告［R］.天津：交通運輸部天津水運工程科學研究所，2008.

Deposition reduction effect study of single?double sediment?protecting dyke with numerical model test

ZHANG Wen?jie
（Cangzhou Port Group Co.,Ltd.,Cangzhou 061113,China）

Based on the analysis of water and sediment conditions in the engineering sea area,the rationality and feasibility of extending single sediment?protecting dyke were demonstrated through deposition reduction effect study of single?double sediment?protecting dyke with numerical model method.This will save the construction cost and improve the comprehensive social benefits.

mathematical model test;single dike;double dike;deposition reduction effect

TV 142；O 242.1

A

1005-8443（2016）03-0247-08

2016-04-21；

2016-05-12

張文杰（1982-），男，河北省南皮人，工程師，主要從事港口與航道工程建設管理工作。

Biography：ZHANG Wen?jie（1982-），male，engineer.