石浦港水動(dòng)力泥沙問題研究

李孟國，麥 苗，李文丹，張義豐

(交通運(yùn)輸部天津水運(yùn)工程科學(xué)研究所 港口水工建筑技術(shù)國家工程實(shí)驗(yàn)室 工程泥沙交通行業(yè)重點(diǎn)實(shí)驗(yàn)室,天津 300456)

石浦港水動(dòng)力泥沙問題研究

李孟國，麥 苗，李文丹，張義豐

(交通運(yùn)輸部天津水運(yùn)工程科學(xué)研究所 港口水工建筑技術(shù)國家工程實(shí)驗(yàn)室 工程泥沙交通行業(yè)重點(diǎn)實(shí)驗(yàn)室,天津 300456)

石浦港是一個(gè)具有五個(gè)口門的狹長(zhǎng)港灣。作為一個(gè)港口，現(xiàn)有航道的通航能力明顯不足，亟需升級(jí)現(xiàn)有航道至五萬噸級(jí)航道以適應(yīng)發(fā)展需要。另外，近年來石浦港大部存在泥沙淤積明顯加重的趨勢(shì)。本文通過地貌特征分析、遙感分析、現(xiàn)場(chǎng)實(shí)測(cè)水文泥沙資料分析、數(shù)學(xué)模型試驗(yàn)和物理模型試驗(yàn)對(duì)石浦港的潮流泥沙特征、泥沙淤積原因、石浦港航道開挖及下灣門炸礁通航等水動(dòng)力泥沙問題進(jìn)行了研究。研究結(jié)果表明：(1)石浦港是三門灣的外海通道之一，潮流是維持水道水深的主要?jiǎng)恿Γ?2)石浦港及附近海域?qū)儆谟倌噘|(zhì)海岸，含沙量較大；(3)五個(gè)口門中，三門和下灣門潮量較大，是主要口門；(4)從水動(dòng)力泥沙角度考慮，石浦港五萬噸級(jí)航道開挖及下灣門炸礁通航是可行的；(5)下洋涂圍墾造成了石浦港潮動(dòng)力減弱，是石浦港西、中部泥沙淤積加重的主要原因。

石浦港；石浦水道；潮汐汊道；水動(dòng)力；泥沙；三門灣；數(shù)學(xué)模型；物理模型

1 研究背景

圖1 石浦水道形勢(shì)圖Fig.1 Sketch of Shipu Waterway

石浦港位于浙江中部沿海的石浦水道內(nèi)，是典型的多口門狹長(zhǎng)潮汐汊道型港灣(圖1)，其是由東門島、對(duì)面山、南田島、高塘島、坦塘島圍列而成的“月牙”狀封閉型港灣，有銅瓦門、東門、下灣門、蜊門及三門共五個(gè)口門與外海相通，灣內(nèi)分布著打鼓峙、中界山、汏網(wǎng)嶼、飯甑山、庵山等島嶼及深槽和淺灘、中心線全長(zhǎng)約20 km，寬0.4～3 km，面積約27 km2。石浦港水深4～61 m(理論基面，下同)，大部分水深在10 m以上，可泊萬艘漁船，可行數(shù)千噸級(jí)海輪；港內(nèi)風(fēng)輕浪平，可避10級(jí)以上大風(fēng)，是東南沿海著名的避風(fēng)良港。目前石浦港的五個(gè)口門均具備一定的通航能力，船舶通行以下灣門、東門及銅瓦門為主，大型船舶通航以下灣門為主，最大通航能力5 000 t[1-2]。

目前石浦港遇到兩個(gè)問題：其一是近幾年來石浦港淤積有明顯加重的趨勢(shì)，原因不明；其二是提高通航能力問題：隨著臨港工業(yè)發(fā)展，現(xiàn)有千噸級(jí)航道的通航能力明顯不足，亟需升級(jí)現(xiàn)有航道，建設(shè)萬噸級(jí)航道以適應(yīng)發(fā)展需要。建設(shè)萬噸級(jí)航道除了在石浦港內(nèi)進(jìn)行局部航道段浚深外，還要在五個(gè)口門中選擇出海口。在五個(gè)口門中，下灣門適合開發(fā)萬噸級(jí)出海航道[1-2]。下灣門水道全長(zhǎng)近3 km，寬0.4～0.8 km，水深20 m以上，口門處水深可達(dá)80 m左右，萬噸級(jí)航道無需開挖，但在口門附近存在金龍礁和石欄礁等礁石群(圖2)，需要炸礁通航。航道建設(shè)的浚深和炸礁對(duì)具有多口門的石浦港水沙環(huán)境有何影響是讓人關(guān)心和需要論證的。

本文通過地貌特征分析、懸沙分布遙感分析、現(xiàn)場(chǎng)實(shí)測(cè)水文泥沙資料分析、波浪潮流數(shù)學(xué)模型和波浪潮流泥沙物理模型對(duì)石浦港的潮流泥沙特征、泥沙淤積原因、五萬噸級(jí)航道開挖及下灣門炸礁通航(圖2)的水動(dòng)力泥沙問題進(jìn)行研究。

圖2 石浦港航道軸線示意圖Fig.2 Sketch of navigational channel axis in Shipu Harbor

2 地貌特征

從地貌上看，石浦水道是三門灣的一部分(圖1)。三門灣是浙江省主要海灣之一，海域面積775 km2，貓頭水道、滿山水道、珠門港水道及石浦水道是其四個(gè)口門段汊道(潮汐汊道)，其中石浦水道與珠門港水道一起控制著三門灣北部的白礁水道、下洋涂東部潮灘區(qū)域[3]。從地貌學(xué)角度看，四個(gè)口門段水道既是三門灣納潮水體的進(jìn)出通道，又要靠三門灣的納潮來維持其水道水深。

3 水動(dòng)力泥沙特征

3.1潮汐與潮流[4-8]

2011年7月24～25日和7月31日～8月1日在石浦港內(nèi)外進(jìn)行了6個(gè)潮位測(cè)站、13條水文垂線的大、小潮全潮同步觀測(cè)(站位見圖1，H6潮位站位于下洋涂西側(cè)的鵓鴣頭漁村，圖1中未能顯示)，根據(jù)實(shí)測(cè)資料，石浦港的潮汐和潮流特征如下：

(1)各站潮汐類型(HK1+HO1)HM2的值均小于0.5，表明潮汐屬正規(guī)半日潮型，潮汐漲潮平均歷時(shí)略短于落潮平均歷時(shí)。

(2)水文全潮測(cè)驗(yàn)期間H1、H2 、H3站的最大潮差分別為5.90 m、5.61 m和5.21 m，最小潮差分別為1.43 m、1.36 m和1.27 m，平均潮差分別為3.21 m、3.05 m和2.83 m。石浦港從西向東潮差逐漸減小。

(3)各測(cè)站潮流類型F=(WK1+WO1)WM2介于0.08～0.26之間，平均為0.16，表明潮流類型屬于規(guī)則半日潮流性質(zhì)，潮流基本上沿岸線或者深槽呈往復(fù)運(yùn)動(dòng)。

(4)總體上講，石浦港漲潮時(shí)潮流從東、南四個(gè)口門(蜊門、下灣門、東門和銅瓦門)流入，從西部三門流出，落潮則相反(與文獻(xiàn)[9]報(bào)道結(jié)論一致)，但除三門和下灣門漲落潮流無位相差外，其與其它各口門有1～2 h的位相差。

(5)根據(jù)V3和V5站資料，石浦港內(nèi)大潮漲落潮平均流速介于0.66～0.71 ms之間，最大流速介于1.14～1.27 ms之間；三門口(V2站)大潮漲落潮平均流速分別為0.50 ms和0.56 ms，最大漲落潮平均流速分別為0.84 ms和1.15 ms；蜊門(V4站)大潮漲落潮平均流速分別為0.69 ms和0.53 ms，最大漲落潮平均流速分別為1.16 ms和0.87 ms；銅瓦門(V6站)大潮漲落潮平均流速分別為0.22 ms和0.62 ms，最大漲落潮平均流速分別為0.42 ms和1.20 ms；東門(V7站)大潮漲落潮平均流速分別為0.54ms和0.58 ms，最大漲落潮平均流速分別為0.85 ms和0.92 ms；下灣門(V8站)大潮漲落潮平均流速分別為0.74 ms和0.61 ms，最大漲落潮流速分別為1.29 ms和1.12 ms。

(6)在五個(gè)口門中，只有蜊門和下灣門的漲潮流速大于落潮流速(其他三個(gè)口門則反之)，只有下灣門大潮漲落潮流速均大于1.0 ms。

3.2波浪

由于受四周島嶼掩護(hù)，外海波浪對(duì)石浦港影響很小。港內(nèi)波浪主要是風(fēng)成浪(本海區(qū)春夏季盛行SW風(fēng)，夏末至秋冬盛行偏N風(fēng)，常風(fēng)向?yàn)镾W和N向；強(qiáng)風(fēng)向?yàn)镹E向，全年>8級(jí)以上的大風(fēng)日數(shù)為98.9 d，大風(fēng)主要由臺(tái)風(fēng)侵襲形成)，波高一般小于0.8 m，即使臺(tái)風(fēng)過境時(shí)，最大波高僅為1.5 m左右[4,7]。

3.3含沙量與懸沙粒徑

根據(jù)2011年7～8月大、小潮水文全潮測(cè)驗(yàn)資料分析，大潮含沙量相對(duì)較大，各測(cè)站大潮漲落潮平均含沙量分別介于0.135～0.343 kgm3和0.190～0.393 kgm3之間，大潮漲落潮最大含沙量分別介于0.272～0.605 kgm3和0.257～0.596 kgm3之間；小潮含沙量明顯小于大潮，除個(gè)別點(diǎn)外，漲落潮平均含沙量和最大含沙量一般均小于0.1 kgm3。

含沙量遙感分析表明[4]：(1)石浦港及附近海域表層含沙量基本介于0.1～0.3 kgm3之間(石浦港內(nèi)含沙量略高)，屬于含沙量相對(duì)較大的海域，含沙量大致呈現(xiàn)從近岸至外海遞減趨勢(shì)；(2)冬季含沙量大，夏季含沙量小，冬季含沙量是夏季含沙量的2～3倍，與文獻(xiàn)[3]結(jié)論一致；(3)含沙量隨潮差大小而變化，潮差大則含沙量大，潮差小則含沙量小，與實(shí)測(cè)資料結(jié)果一致；(4)近岸水體含沙量受風(fēng)影響較大，大風(fēng)浪天水體含沙量相對(duì)增高；(5)綜合以上4點(diǎn)，從遙感反演結(jié)果看，泥沙來源主要為近岸淺灘泥沙的起動(dòng)和搬運(yùn)。

根據(jù)2011年7～8月大、小潮水文全潮測(cè)驗(yàn)期間各站漲憩、落憩、漲急、落急時(shí)段取懸沙水體中的沙樣分析，懸沙中值粒徑變化范圍為0.003 7～0.013 7 mm，平均d50為0.008 2 mm。

3.4底質(zhì)

根據(jù)2011在石浦港內(nèi)外85個(gè)底質(zhì)取樣分析，石浦港內(nèi)底質(zhì)以粘土質(zhì)粉砂為主，其中值粒徑基本介于0.005～0.012 5 mm之間，分選系數(shù)基本介于1.4～2.0之間，分選屬中常；粘土百分含量為20%～40%之間，為典型的淤泥質(zhì)海岸。

3.5單寬輸沙率

根據(jù)2011年7～8月在石浦港內(nèi)外進(jìn)行大、小潮水文全潮測(cè)驗(yàn)，各口門及附近測(cè)點(diǎn)的單寬輸沙量統(tǒng)計(jì)計(jì)算結(jié)果見表1。

表1 水文全潮測(cè)站單寬輸沙量(t全潮)Tab.1 Unit width sediment transport in a complete tidal cycle at observation stations t

表1 水文全潮測(cè)站單寬輸沙量(t全潮)Tab.1 Unit width sediment transport in a complete tidal cycle at observation stations t

站位漲潮大潮小潮平均落潮大潮小潮平均凈輸沙量漲潮-落潮V139.5600.81720.18922.5630.49411.5298.660V2129.5289.04369.286178.8518.03393.442-24.156V396.1006.95451.527119.4378.03863.738-12.211V4113.6016.10759.85487.573.89045.7314.124V635.69529.54032.618172.39224.12298.257-65.639V739.9706.56123.26656.4486.80731.628-8.362V8333.11733.820183.469353.89723.656188.777-5.308V1179.9913.51841.75543.9972.31423.15618.599

從表1的單寬凈輸沙量看，泥沙從三門和蜊門向石浦港內(nèi)輸送，從下灣門、東門和銅瓦門向港外輸沙，與有關(guān)文獻(xiàn)[2,9,10]結(jié)論相同。

3.6泥沙來源

石浦港位于三門灣北岸，附近無大河注入，因此受河流來沙的直接影響很小。從宏觀上來說，長(zhǎng)江口外泄泥沙隨江浙沿岸流南下擴(kuò)散并在本海域近岸淺灘不斷淤積，是形成石浦港附近淤泥質(zhì)海岸的主要泥沙來源。石浦港附近近岸淺灘的細(xì)顆粒泥沙在波浪和潮流共同作用下懸浮再搬運(yùn)并隨潮流進(jìn)入石浦港內(nèi)及石浦港內(nèi)底床泥沙在波流作用下懸浮再搬運(yùn)是當(dāng)前石浦港內(nèi)淤積的主要泥沙來源。

3.7沖淤演變

根據(jù)文獻(xiàn)[2],石浦港港域在1981～1995年的14 a中，大面積處于淤積狀態(tài)，淤積厚度一般在1～1.5 m以內(nèi)。年淤積強(qiáng)度為0.053 ma的區(qū)域占總淤積面積的13，淤積強(qiáng)度大于或等于0.22 ma的占13弱，其余13強(qiáng)的面積淤積強(qiáng)度均小于或等于0.12 ma，整個(gè)港域14 a中的淤積強(qiáng)度為0.13 ma。

根據(jù)2006年6～7月1:1萬和2011年7～8月1:1萬水深測(cè)圖對(duì)比分析[4]，石浦港呈大部分淤積、局部區(qū)域沖刷的狀態(tài)。5 a內(nèi)平均淤厚0.60 m，年淤積強(qiáng)度0.12 ma，其中三門口—打鼓峙和中界山—?dú)W(wǎng)嶼淤積較重，年均淤強(qiáng)為0.15～0.18 ma，打鼓峙—中界山淤積較輕，年均淤強(qiáng)約0.06 m。銅瓦門淤強(qiáng)0.18 ma，東門和下灣門沖刷強(qiáng)度分別為0.29 ma和0.47 ma。

根據(jù)2012～2014年石浦港內(nèi)3個(gè)固定斷面(圖2中的D1、D2和D3)共7次1:1 000水深測(cè)圖并結(jié)合2006年和2011年的水深測(cè)圖分析，2012～2014年期間，石浦港東部的D1斷面呈灘淤槽沖、整體略有沖刷的趨勢(shì)，整體沖刷了0.33 m，淤強(qiáng)為-0.13 ma；石浦港中部的D2斷面整體淤積了0.48 m，淤強(qiáng)為0.21 ma；石浦港西部的D3斷面整體呈淤積趨勢(shì)，整體淤積了1.85 m，淤強(qiáng)為0.82 ma。

固定斷面地形在2011年前后兩個(gè)時(shí)段沖淤變化趨勢(shì)存在明顯差異，2011年之前整體呈基本穩(wěn)定狀態(tài)；2011年后則出現(xiàn)淤積加重趨勢(shì)，有西部淤積較重、中部次之、東部深槽略呈沖刷的變化特點(diǎn)。

4 模型試驗(yàn)

通過數(shù)學(xué)模型和物理模型試驗(yàn)，對(duì)石浦港5萬噸級(jí)航道開挖和下灣門口門炸礁工程的水沙環(huán)境問題進(jìn)行了研究。五萬噸級(jí)航道軸線如圖2所示，航道底寬165 m，底高程-14.8 m(國家85基面)。航道有2處開挖段，打鼓峙西段長(zhǎng)約3.55 km，最大挖深4.39 m，中界山東段長(zhǎng)約3.76 km，最大挖深6.84 m。下灣門水道無需開挖，但要對(duì)石欄礁炸礁和園山突嘴礙航部分炸除。

4.1數(shù)學(xué)模型和物理模型的建立

基于Mike系列軟件中的三角形網(wǎng)格波浪SW模型、水動(dòng)力模塊(FM模塊)、平面二維懸沙輸移擴(kuò)散方程和床面沖淤變化方程，建立了考慮波浪作用的潮流泥沙數(shù)學(xué)模型[5,8,11]。模型范圍的邊界北至韭山列島、南到東磯列島、外海至約-30 m等深線，包括了整個(gè)三門灣海域。為擬合復(fù)雜岸線和航道等細(xì)致建筑物邊界，采用了無結(jié)構(gòu)三角形網(wǎng)格對(duì)計(jì)算域進(jìn)行剖分。為模擬石浦港航道工程的細(xì)部流態(tài)，模型網(wǎng)格最小空間步長(zhǎng)應(yīng)保證精確分辨航道及港內(nèi)各水道，網(wǎng)格節(jié)點(diǎn)90 000個(gè)，其中最小網(wǎng)格尺度近10 m。外海潮位邊界由Chinatide提供[12]。

圖3 物理模型試驗(yàn)布置Fig.3 Layout of the physical model

基于幾何相似、潮流運(yùn)動(dòng)和泥沙運(yùn)動(dòng)相似的要求，建立了定床懸沙物理模型[6]。模型按包括下灣門炸礁工程區(qū)及石浦港內(nèi)的航道布置，南北方向包括東門島—南田島的范圍，西側(cè)至三門口，東側(cè)至檀頭山，即東西長(zhǎng)約22.5 km，南北寬約11.25 km的海域，面積約253.1 km2。模型采用水平比尺為300，根據(jù)潮流運(yùn)動(dòng)和泥沙運(yùn)動(dòng)相似以及選擇模型沙的要求，選取模型垂直比尺為70，模型變率4.29。模型占地86 m×42 m，見圖3。本模型各開邊界的潮汐控制均采用可逆泵生潮系統(tǒng)來實(shí)施，在銅瓦門、東門、下灣門、汏網(wǎng)嶼、中界山三門口各布置一條加沙斷面。

工程海域懸沙中值粒徑介于0.003 7～0.013 7 mm之間，平均中值粒徑約為0.008 2 mm左右，屬粘土質(zhì)粉砂。這種懸沙是會(huì)產(chǎn)生絮凝的，對(duì)應(yīng)的靜水沉速約為0.48 mms左右。物理模型選擇密度為1.45 gcm3的寧夏煤作為模型懸沙。在模型平面比尺為300，垂直比尺為70的條件下，計(jì)算沉速比尺為1.95，選取模型沙平均中值粒徑為0.035 mm，模型沙沉速約為0.25 mms左右，換算成實(shí)際沉速比尺為1.92，基本可以滿足泥沙沉降相似要求。

數(shù)學(xué)模型和物理模型均使用2011年7～8月大、小潮水文全潮資料及石浦港內(nèi)2006～2011年的淤積資料進(jìn)行了驗(yàn)證[5,6,8,11]，均符合《海岸與河口潮流泥沙模擬技術(shù)規(guī)程》[13]的要求。

4.2主要試驗(yàn)成果

(1)航道流速變化[5,6,8,11,14]。航道工程實(shí)施后，僅在沿航道軸線兩側(cè)局部區(qū)域流速有所變化。石浦港東西向的三門口至汏網(wǎng)嶼航道段水域，受航道開挖的影響，水流存在歸槽現(xiàn)象，航道流速略有增大，漲落潮平均流速增加0.01～0.06 ms，流速增幅介于5%～6%之間。下灣門水道的汏網(wǎng)嶼至園山航道段，漲落潮流速呈增加的趨勢(shì)，平均流速增加0.02～0.04 ms，流速增幅介于3%～4%之間；下灣門水道的石欄礁附近航道段，由于石欄礁被炸除，水流通暢，漲落潮平均流速分別增加了0.33 ms和0.34 ms，增幅約為73%和45%。

航道線上的水流基本呈現(xiàn)為平行于航道軸線的往復(fù)流，滿足航行安全需要，航道線的布置是合理可行的。

(2)口門內(nèi)外潮位變化。航道工程實(shí)施后，數(shù)學(xué)模型試驗(yàn)的下灣門口門內(nèi)外高低潮位偏差均在±1 cm以內(nèi)，物理模型試驗(yàn)的高低潮位變化幅度在±0.02 m之內(nèi)，處于試驗(yàn)誤差范圍，因此，航道工程實(shí)施后對(duì)潮位基本沒有影響。

(3)各口門潮量、分流比及其變化。數(shù)學(xué)模型計(jì)算和物理模型試驗(yàn)得到的大潮過程各口門現(xiàn)狀漲落潮潮量及其工程引起的變化見表2。從表2可見：①下灣門和三門是石浦港主要的口門，無論是下灣門還是三門，其漲落潮潮量均與其他三個(gè)口門的漲落潮潮量之和相當(dāng)；②對(duì)于石浦港多個(gè)口門來說，除三門外，其他四個(gè)口門基本上呈現(xiàn)出漲落潮同步，數(shù)模計(jì)算的銅瓦門、東門、下灣門和蜊門四個(gè)口門的漲潮分流比為0.229:0.164:0.516:0.091，落潮分流比為0.284:0.169:0.481:0.067，下灣門是優(yōu)勢(shì)口門；③航道工程實(shí)施后，三門和下灣門漲落潮潮量呈增加趨勢(shì)，增幅2%～5%(三門增加2%～5%，下灣門增加2%～4%)，其他三個(gè)口門漲落潮潮量呈減小趨勢(shì)，減幅在3%以內(nèi)，不會(huì)引起口門萎縮[5-6]。

表2 各口門大潮過程漲落潮潮量及工程引起的變化Tab.2 Tidal prisms and their changes caused by the project at five entrances during a spring tide process

(5)中界山南北分流比及其變化[6]。天然狀態(tài)下，中界山北側(cè)和南側(cè)漲潮分流比約為63：37，落潮分流比約為59：41；航道工程實(shí)施后，漲潮分流比調(diào)整為67：33，落潮分流比調(diào)整為62：38。分流比及變化說明，一方面，中界山北側(cè)是航道所在水域，流量增加意味著對(duì)航道維護(hù)有利，另一方面，航道工程對(duì)石浦港影響是很小的。

表3 泥沙淤積試驗(yàn)結(jié)果Tab.3 Results of sedimentation test

(6)航道泥沙淤積。航道工程實(shí)施后，數(shù)學(xué)模型和物理模型的泥沙淤積試驗(yàn)結(jié)果見表3。從表3可見，物理模型和數(shù)學(xué)模型的結(jié)果基本一致，石浦港內(nèi)的航道開挖淤積不大，淤積量約90萬m3。

(7)下洋涂圍墾對(duì)石浦港泥沙淤積的影響。下洋涂圍墾工程位于石浦港三門口外部，圍墾面積35.78 km2。利用建立的數(shù)學(xué)模型模擬了下洋涂圍墾工程對(duì)石浦港淤積的影響。模擬結(jié)果表明，①圍墾工程實(shí)施后，石浦港內(nèi)全潮平均流速減小約0.05～0.12 ms，最大流速減小約0.07～0.21 ms；②圍墾工程實(shí)施后，石浦港內(nèi)發(fā)生泥沙淤積，平衡后的淤積分布見圖4，三門口最重，向西淤積厚度逐漸降低。泥沙淤積厚度與前述的2012～2014石浦港內(nèi)固定斷面D3和D2的實(shí)際淤積結(jié)果趨勢(shì)和量級(jí)吻合。

圖4 石浦港內(nèi)淤積分布Fig.4 Distribution of sedimentation in Shipu Harbor

如前所述，石浦水道與珠門港水道一起控制著三門灣北部的白礁水道、下洋涂東部潮灘區(qū)域，下洋涂圍墾工程實(shí)施后，造成納潮面積減小(有關(guān)文獻(xiàn)[15-16]結(jié)論一致)，因此石浦港流速肯定減小，流速減小造成泥沙淤積。

5 結(jié)論

本文根據(jù)地貌特征、現(xiàn)場(chǎng)實(shí)測(cè)資料和懸沙遙感衛(wèi)片等的分析和數(shù)學(xué)模型、物理模型試驗(yàn)對(duì)石浦港的水動(dòng)力泥沙問題進(jìn)行了研究，有如下主要結(jié)論：(1)擁有五個(gè)口門的石浦港是三門灣的四個(gè)口門潮汐汊道之一，是三門灣納潮水體的進(jìn)出通道，潮動(dòng)力維持著石浦港的水深，是天然避風(fēng)良港；(2)石浦港漲潮時(shí)潮流由東南四個(gè)口門進(jìn)入，由西部三門流出，落潮則相反。但門各口門存在1～2 h的位相差；(3)石浦港潮汐屬于正規(guī)半日潮型，從西向東潮差逐漸減小；潮流類型屬于規(guī)則半日潮流性質(zhì)，潮流基本上沿岸線或者深槽呈往復(fù)運(yùn)動(dòng)；(4)石浦港潮差大，潮流動(dòng)力強(qiáng)。下灣門大潮漲落潮最大流速均在1 ms以上；(5)石浦港海域的含沙量較大，泥沙隨漲落潮向三門和蜊門石浦港內(nèi)輸送，從下灣門、東門和銅瓦門向石浦港外輸送；(6)石浦港海域?qū)儆谟倌噘|(zhì)海岸，懸沙粒徑和底質(zhì)粒徑基本一致；(7)石浦港附近近岸淺灘的細(xì)顆粒泥沙在波浪和潮流共同作用下懸浮再搬運(yùn)并隨潮流進(jìn)入石浦港內(nèi)及石浦港內(nèi)底床泥沙在波流作用下懸浮再搬運(yùn)是石浦港內(nèi)淤積的主要泥沙來源；(8)石浦港航道開挖及下灣門炸礁通航后，基本不改變港內(nèi)水沙環(huán)境，航道水流較平順，各口門分流比及中界山南北分流比均略有調(diào)整，對(duì)潮位無明顯影響，淤積量不大，航道工程是可行的；(9)下洋涂圍墾是造成石浦港西、中部泥沙淤積加重的主要原因。

[1]盧桂營,張義豐,王劍波.石浦港水沙環(huán)境分析[J].浙江水利科技,2013 (4): 39-41.

LU G Y, ZHANG Y F, WANG J B. Analyses of hydrodynamic sediment environment in Shipu Harbor[J]. Zhejiang Hydrotechnics, 2013 (4): 39-41.

[2]王義剛,梁亮,王震.石浦港口門炸礁工程對(duì)港內(nèi)的影響研究[C]中國海洋工程學(xué)會(huì).第十二屆中國海岸工程學(xué)術(shù)研討會(huì)論文集. 北京：海洋出版社,2005.

[3]夏小明. 三門灣潮汐汊道系統(tǒng)的穩(wěn)定性[D].杭州：浙江大學(xué),2011.

[4]麥苗,蔡寅,韓志遠(yuǎn).石浦港航道工程自然條件和岸灘演變分析報(bào)告[R].天津：交通運(yùn)輸部天津水運(yùn)工程科學(xué)研究所,2012.

[5]李文丹,解鳴曉,李孟國.石浦港航道工程潮流泥沙數(shù)學(xué)模型試驗(yàn)研究報(bào)告[R].天津：交通運(yùn)輸部天津水運(yùn)工程科學(xué)研究所,2012.

[6]閆勇,劉國亭,麥苗, 等. 石浦港航道工程潮流泥沙物理模型試驗(yàn)研究報(bào)告[R].天津：交通運(yùn)輸部天津水運(yùn)工程科學(xué)研究所,2012.

[7]麥苗，韓志遠(yuǎn).石浦港地形沖淤分析報(bào)告[R].天津：交通運(yùn)輸部天津水運(yùn)工程科學(xué)研究所,2014.

[8]李文丹,張義豐.寧波舟山港石浦港區(qū)主航道一期工程潮流數(shù)學(xué)模型研究報(bào)告[R].天津：交通運(yùn)輸部天津水運(yùn)工程科學(xué)研究所,2015.

[9]沈和誠. 石浦港淤積的初步探討[J].水運(yùn)工程，1989(9) : 40-43.

SHEN H C. A preliminary investigation of sedimentation in Shipu Harbor[J]. Port & Waterway Engineering, 1989(9) : 40-43.

[10]李金鐸. 石浦港下灣門水道的水沙運(yùn)移特征[J]. 海洋科學(xué),2008, 32(1) :56-59.

LI J D. Hydrological and sediment features and transportation in Xiawanmen Channel, the Shipu Harbor[J].Marine Sciences, 2008, 32(1) :56-59.

[11]李文丹,李孟國,解鳴曉. 石浦港航道工程潮流泥沙數(shù)值模擬[J].水運(yùn)工程,2014(4): 24-29.

LI W D, LI M G, XIE M X. Tidal current and sediment numerical simulation of channel project in Shipu Harbor[J]. Port & Waterway Engineering, 2014(4): 24-29.

[12]李孟國,鄭敬云. 中國海域潮汐預(yù)報(bào)軟件Chinatide的應(yīng)用[J].水道港口,2007,28(1):65-68.

LI M G, ZHENG J Y. Introduction to Chinatide software for tide prediction in China seas[J]. Journal of Waterway and Harbor, 2007,28(1):65-68.

[14]溫春鵬,李孟國,李文丹. 石浦港下灣門航道工程潮流數(shù)值模擬研究[J].水道港口,2013,34(1):19-25,88.

WEN C P, LI M G, LI W D. Tidal current numerical simulation of Xiwanmen channel project of Shipu Harbor[J]. Journal of Waterway and Harbor, 2013,34(1):19-25,88.

[15]彭姱容,姚炎明,陳琴.石浦港納潮量對(duì)三門灣圍填海的累積響應(yīng)[J].海洋湖沼通報(bào),2014(1):169-178.

PENG K R, YAO Y M, CHEN Q. Cumulative response of tidal prism of Shipu Channel to reclamation projects in Sanmen Bay[J].Transactions of Oceanology and Limnology, 2014(1):169-178.

[16]許雪峰,羊天柱.三門灣多個(gè)圍墾工程的整體影響數(shù)學(xué)模型[J].海洋學(xué)研究,2006, 24(s):49-59.

XU X F, YANG T Z. Mathematical model study on overall impact of Sanwanmen Bay reclamation project[J].Journal of Marine Sciences, 2006, 24(s):49-59.

Study on hydrodynamic sediment problems in Shipu Harbor

LIMeng-guo,MAIMiao,LIWen-dan,ZHANGYi-feng

(TianjinResearchInstituteforWaterTransportEngineering,NationalEngineeringLaboratoryforPortHydraulicConstructionTechnology,KeyLaboratoryofEngineeringSedimentoftheMinistryofTransport,Tianjin300456,China)

Shipu Harbor is a long-narrow harbor with five entrances. As a harbor, its present navigational capacity of navigation channel is insufficient, and it is desiderate to upgrade the existing navigation channel to 50 000 DWT channel to meet the need of development. In addition，it is found that the sedimentation in most part of Shipu Harbor tends to aggravate in recent years. In this paper, by means of analyses of landform features, remote sensing and in-situ hydrology and sediment data，mathematical model and physical model, the hydrodynamic sediment environment problems were studied related to tidal current and sediment characteristics in Shipu Harbor, the sedimentation causes, navigation channel excavation and reef explosion for navigation in Xiawanmen entrances, etc. The study results show that, (1) Shipu Harbor is one of the out-to-sea passages of Sanmen Bay, and the tidal current is the main hydrodynamic force to maintain its water depth. (2) Shipu Harbor and nearby is of muddy coast, and the sediment concentration is higher. (3)Among the five entrances of Shipu Harbor, Sanmen and Xiawanmen are the main entrances with relatively greater tidal capacity. (4) It is feasible to excavate 50 000 DWT navigation channel in Shipu Harbor and to make explosion reef for navigation in Xiawanmen entrances as far as hydrodynamic sediment is concerned. (5) Xiayangtu reclamation weakened the tidal force in Shipu Waterway, and is the main cause to aggravate the sedimentation in the west part and middle part of Shipu Harbor.

Shipu Harbor; Shipu Waterway; tidal inlet; hydrodynamics; sediment; Sanmen Bay；mathematical model; physical model

2017-02-20；

2017-06-21

國家重點(diǎn)研發(fā)計(jì)劃(課題編號(hào)：2016YFC0402605)

李孟國(1964-)，男，天津市人，研究員，主要從事海岸河口水動(dòng)力泥沙研究。

Biography：LI Meng-guo(1964-),male, professor.

TV 142;O 242.1

A

1005-8443(2017)05-0433-07