坎門(mén)漁港三期工程潮流數(shù)模匯報(bào)

1、河海大學(xué)河海大學(xué)2012.08.16 2一、前言二、自然條件三、坎門(mén)漁港潮流數(shù)學(xué)模型 四、防波堤工程潮流數(shù)值模擬五、港內(nèi)泥沙沖淤分析六、結(jié)論與建議主要內(nèi)容一、前言1.1 研究目的研究目的 玉環(huán)縣擬在坎門(mén)灣口的南排山與黃門(mén)山之間進(jìn)行三期防波堤工程建設(shè)，將坎門(mén)漁港建成深水港。為了從潮流變化和港內(nèi)回淤方面探討防波堤的合理布置，開(kāi)展本項(xiàng)研究，成果作為坎門(mén)漁港三期（南排山外黃門(mén)山）防波堤工程可行性研究的依據(jù)。 一、前言1.2 工作依據(jù)工作依據(jù) “浙江省玉環(huán)縣坎門(mén)漁港三期（南排山外黃門(mén)山）防波堤工程潮流數(shù)學(xué)模型及泥沙分析” 合同書(shū)，2010年6月。 中華人民共和國(guó)交通運(yùn)輸部發(fā)布的海港水文規(guī)范（JTJ213

2、-98）。 中華人民共和國(guó)交通運(yùn)輸部發(fā)布的海岸與河口潮流泥沙模擬技術(shù)規(guī)程（JTS/T231-2-2010）。 一、前言溫州灣及附近（圖號(hào)13710，比例尺1:100000，2004年5月出版，中國(guó)人民解放軍海軍司令部航海保證部，理論深度基準(zhǔn)面）。漩門(mén)灣及附近（圖號(hào)13711，比例尺1:25000，2005年10月出版，中國(guó)人民解放軍海軍司令部航海保證部，理論深度基準(zhǔn)面）?？查T(mén)漁港水下地形測(cè)量圖（比例尺1:5000，2001年11月施測(cè)，溫州綜合測(cè)繪院，1985國(guó)家高程基準(zhǔn)面）。長(zhǎng)江委水文局長(zhǎng)江口水文水資源勘測(cè)局編制的玉環(huán)縣坎門(mén)南排山避風(fēng)港防波堤工程水文測(cè)驗(yàn)技術(shù)報(bào)告，2009年11月。 一、前言

3、1.3 研究?jī)?nèi)容和技術(shù)路線研究?jī)?nèi)容和技術(shù)路線 在分析地形、水沙因素等自然條件的基礎(chǔ)上，建立平面二維潮流數(shù)學(xué)模型，并進(jìn)行模型驗(yàn)證計(jì)算；然后，進(jìn)行工程方案的潮流模擬，比較工程前、后流場(chǎng)的變化（著重比較工程前、后流速的變化）；最后，計(jì)算分析泥沙沖淤特征，比較各工程方案的優(yōu)劣。 二、自然條件 2.1 地形地形 二期防波堤深度：m，理論最低潮面門(mén)山黃南排山黃門(mén)山內(nèi)坎門(mén)媽祖印礁10.010.010.010.010.020.00.00.02.04.00.0 內(nèi)港水域的等深線呈現(xiàn)以二期防波堤口門(mén)為中心的弧形，口門(mén)附近水深最深，離口門(mén)愈遠(yuǎn)，水深愈淺，西北部水深大多不足0m。 坎門(mén)灣口有南排山和黃門(mén)山，形成三條水

4、道與外海相通：（1）東口門(mén)，介于南排山與北側(cè)坎門(mén)頭之間，寬約0.7km，最小水深6m，其間有寬約0.3km、水深大于10m的深槽；（2）南口門(mén)，介于黃門(mén)山與南排山之間，寬約2.5km，水深1014m；（3）西口門(mén)，介于黃門(mén)山與小門(mén)仔之間，寬約0.5km，水深多大于10m。 2.2 潮位潮位 根據(jù)坎門(mén)海洋站多年實(shí)測(cè)潮位資料分析，坎門(mén)灣的潮汐屬正規(guī)半日潮；潮差較大，多年平均潮差為4.05m；漲潮歷時(shí)比落潮歷時(shí)略長(zhǎng)，多年平均漲潮歷時(shí)為6h18min，多年平均落潮歷時(shí)為6h7min。 2.3 潮流潮流 坎門(mén)灣的潮流有往復(fù)流、旋轉(zhuǎn)流和混合流三種類(lèi)型。 坎門(mén)灣的潮流屬駐波型。漲、落急出現(xiàn)在中潮位附近，漲、

5、落憩出現(xiàn)在高、低潮位前后。 二、自然條件 2.4 泥沙泥沙坎門(mén)灣為淤泥質(zhì)海灣，懸移質(zhì)造床。東口門(mén)和西口門(mén)的底質(zhì)相對(duì)較粗，中值粒徑分別為0.018mm和0.022mm，其他水域的底質(zhì)中值粒徑為0.0070.008mm。懸移質(zhì)中值粒徑為0.008mm左右?？查T(mén)灣含沙量大小按大、中、小潮的順序依次遞減，漲潮含沙量大于落潮含沙量；在灣口的三個(gè)口門(mén)中，南口門(mén)的含沙量最大，西口門(mén)次之，東口門(mén)最小；坎門(mén)灣灣口的平均含沙量為0.3kg/m3左右。 二、自然條件 10三、坎門(mén)漁港潮流數(shù)學(xué)模型3.1 數(shù)學(xué)模型數(shù)學(xué)模型 水流運(yùn)動(dòng)的基本方程包括水體連續(xù)性方程和動(dòng)量守恒方程，即 0)()(yHVxHUtH)()()()

6、()(2222yUxUHNfHVxZHgHyHVUxHUUtHUxbbxsx)()()()()(2222yVxVHNfHUyZHgHyHVVxHUVtHVybbysy3.2 計(jì)算模式計(jì)算模式 采用有限體積法離散上述基本方程；通過(guò)離散得到的線性代數(shù)方程組，用高斯消去法求解。 11三、坎門(mén)漁港潮流數(shù)學(xué)模型3.3 模型范圍及網(wǎng)格尺寸模型范圍及網(wǎng)格尺寸圖3.1 潮流數(shù)模計(jì)算范圍及網(wǎng)格剖分 模型范圍，南起北緯275822，北至北緯281145，南北間距24.8km；西起東經(jīng)1211046，東至東經(jīng)1212315，東西間距20.5km；其面積為508.4km2。 采用正方形網(wǎng)格剖分計(jì)算域，網(wǎng)格尺寸為20m

7、20m。 123.4 邊界條件邊界條件3.5 參數(shù)確定參數(shù)確定 曼寧系數(shù)：n0.020三、坎門(mén)漁港潮流數(shù)學(xué)模型+ + + + + + + + + + + + 長(zhǎng)海大青蘇爺島大連秦皇島連云港大戢山濟(jì)州島福江寶島儀間大陳島平壇臺(tái)灣朝鮮半島九州圖3.2 中國(guó)東部海區(qū)潮波模擬計(jì)算網(wǎng)格圖南澳布袋蘇澳24N 26N 28N 30N 32N 34N 36N 38N 40N 117E 119E 121E 123E 125E 127E 129E Scale0 100 200km40 N 24N 26N 28N 30N 32N 34N 36N 38N 40N 117E 119E 121E 123E 125E 12

8、7E 129E Scale0 100 200km40 N 35 N 30 N 0 100 200km25 N 35 N 30 N 0 100 200km25 N 模型的水邊界采用潮位控制，潮位過(guò)程由東中國(guó)海潮波數(shù)學(xué)模型提供。 對(duì)于隨潮位變化的淺灘動(dòng)邊界，采用網(wǎng)格干、濕判斷法進(jìn)行處理。 13三、坎門(mén)漁港潮流數(shù)學(xué)模型3.6 驗(yàn)證計(jì)算驗(yàn)證計(jì)算 模型驗(yàn)證取用2009年11月3日11月10日的實(shí)測(cè)大、中、小潮水文資料，測(cè)站布設(shè)見(jiàn)圖2.3，有2個(gè)潮位測(cè)站（W1、W2）和6條測(cè)流垂線（S1 S6）。 潮位站測(cè)流點(diǎn)圖2.3 現(xiàn)場(chǎng)水文泥沙測(cè)驗(yàn)站位圖14三、坎門(mén)漁港潮流數(shù)學(xué)模型圖3.3(a) 大潮潮位驗(yàn)證曲線黃門(mén)

9、山水產(chǎn)碼頭 *：實(shí)測(cè)值 -：計(jì)算值圖3.4(a) 大潮流速驗(yàn)證曲線(2009年11月3日) *：實(shí)測(cè)值 -：計(jì)算值大潮驗(yàn)證曲線15三、坎門(mén)漁港潮流數(shù)學(xué)模型圖3.3(b) 中潮潮位驗(yàn)證曲線黃門(mén)山水產(chǎn)碼頭 *：實(shí)測(cè)值 -：計(jì)算值 *：實(shí)測(cè)值 -：計(jì)算值(2009年11月6日)圖3.4(b) 中潮流速驗(yàn)證曲線中潮驗(yàn)證曲線16三、坎門(mén)漁港潮流數(shù)學(xué)模型圖3.3(c) 小潮潮位驗(yàn)證曲線黃門(mén)山水產(chǎn)碼頭 *：實(shí)測(cè)值 -：計(jì)算值(2009年11月9日)圖3.4(c) 小潮流速驗(yàn)證曲線 *：實(shí)測(cè)值 -：計(jì)算值小潮驗(yàn)證曲線17三、坎門(mén)漁港潮流數(shù)學(xué)模型流場(chǎng)特征流場(chǎng)特征圖3.5(a) 大潮流場(chǎng)圖（工程前 漲急）水產(chǎn)碼頭

10、潮位相位示意圖 （1）計(jì)算水域的潮汐屬正規(guī)半日潮型，潮波為駐潮波。 （2）漲潮時(shí)，潮水由東口門(mén)和南口門(mén)進(jìn)入坎門(mén)灣內(nèi)，其中，由南口門(mén)進(jìn)入灣內(nèi)的水流大部分通過(guò)外黃門(mén)（西口門(mén)）流出，然后沿西南方向流往甌江口；由東口門(mén)進(jìn)入灣內(nèi)的水流主要經(jīng)過(guò)二期防波堤口門(mén)進(jìn)入坎門(mén)漁港現(xiàn)有港區(qū)（內(nèi)港區(qū)）。三、坎門(mén)漁港潮流數(shù)學(xué)模型圖3.5(b) 大潮流場(chǎng)圖（工程前 落急）水產(chǎn)碼頭潮位相位示意圖 （3）落潮時(shí)，由西口門(mén)和南口門(mén)西側(cè)（馬祖印礁至黃門(mén)山之間）流來(lái)的潮水由南口門(mén)東側(cè)（馬祖印礁至南排山之間）退出；坎門(mén)漁港現(xiàn)有港區(qū)的潮水經(jīng)過(guò)二期防波堤口門(mén)由東口門(mén)退出。 （方案1：南防波堤，東偏北25走向，2310m長(zhǎng)）圖4.1(a)

11、三期防波堤工程平面布置圖 圖4.1(b) 三期防波堤工程平面布置圖 （方案2：南防波堤，東偏北25走向，2310m長(zhǎng)；東防波堤，南北走向，200m長(zhǎng)） 玉環(huán)縣坎門(mén)漁港三期（南排山-黃門(mén)山）防波堤工程擬建南防波堤和東防波堤，南防波堤位于南排山與黃門(mén)山之間的南口門(mén)，東防波堤南排山與坎門(mén)頭之間的東口門(mén)。根據(jù)當(dāng)?shù)氐牡匦蔚孛布安ɡ颂攸c(diǎn)，防波堤有六種平面布置方案。4.1 三期防波堤工程布置方案三期防波堤工程布置方案四、防波堤工程潮流數(shù)值模擬圖4.1(c) 三期防波堤工程平面布置圖 （方案3：南防波堤，東偏北25走向，2310m長(zhǎng)；東防波堤，南北走向，250m長(zhǎng)）四、防波堤工程潮流數(shù)值模擬圖4.1(d) 三

12、期防波堤工程平面布置圖 （方案4：南防波堤堤根位于南排山，西偏南25走向，1755m長(zhǎng)）圖4.1(e) 三期防波堤工程平面布置圖 （方案5：南防波堤堤根位于黃門(mén)山，東偏北25走向，1473m長(zhǎng)）四、防波堤工程潮流數(shù)值模擬（方案6：南防波堤的一條堤堤根位于黃門(mén)山，東偏北20走向，1443m長(zhǎng)；另一條堤的堤根位于南排山，西偏南50走向，859m長(zhǎng)）圖4.1(f) 三期防波堤工程平面布置圖 （方案5：南防波堤堤根位于黃門(mén)山，東偏北25走向，1473m長(zhǎng)）圖4.6(a) 大潮漲急流場(chǎng)圖水產(chǎn)碼頭潮位相位示意圖四、防波堤工程潮流數(shù)值模擬4.2 成果分析成果分析（方案5：南防波堤堤根位于黃門(mén)山，東偏北25走

13、向，1473m長(zhǎng)）圖4.6(b) 大潮落急流場(chǎng)圖水產(chǎn)碼頭潮位相位示意圖 方案4、方案5及方案6的南防波堤未完全封堵南口門(mén)，坎門(mén)漁港有東口門(mén)、新的南口門(mén)及外黃門(mén)（西口門(mén)）與外海相通。漲潮時(shí)，潮水由東口門(mén)和新的南口門(mén)進(jìn)入港內(nèi)，由外黃門(mén)（西口門(mén)）流出；落潮時(shí)，潮水由外黃門(mén)（西口門(mén)）流入港內(nèi)，由東口門(mén)和新的南口門(mén)退出。這表明，三個(gè)方案的流路與工程前相同。但是，進(jìn)出新的南口門(mén)和外黃門(mén)（西口門(mén) ） 的 潮 量 減 少 了 ， 而 進(jìn) 出 東 口 門(mén) 的 潮 量 增 加 了 。 （方案1：南防波堤，東偏北25走向，2310m長(zhǎng)）水產(chǎn)碼頭潮位相位示意圖圖4.2(a) 大潮漲急流場(chǎng)圖四、防波堤工程潮流數(shù)值模擬（

14、方案1：南防波堤，東偏北25走向，2310m長(zhǎng)）水產(chǎn)碼頭潮位相位示意圖圖4.2(b) 大潮落急流場(chǎng)圖 方案1的南防波堤完全封堵南口門(mén)，坎門(mén)漁港只有東口門(mén)和外黃門(mén)（西口門(mén)）與外海相通。漲潮時(shí)，潮水由東口門(mén)進(jìn)入港內(nèi)，由外黃門(mén)（西口門(mén)）流出；落潮時(shí)，潮水由外黃門(mén)（西口門(mén)）流入港內(nèi)，由東口門(mén)退出。此時(shí)，進(jìn)出外黃門(mén)（西口門(mén)）的潮量進(jìn)一 步 減 少 ， 而 進(jìn) 出 東 口 門(mén) 的 潮 量 進(jìn) 一 步 增 加 。 （方案3：南防波堤，東偏北25走向，2310m長(zhǎng)；東防波堤，南北走向，250m長(zhǎng)）水產(chǎn)碼頭潮位相位示意圖圖4.4(a) 大潮漲急流場(chǎng)圖（方案3：南防波堤，東偏北25走向，2310m長(zhǎng)；東防波堤，南

15、北走向，250m長(zhǎng)）水產(chǎn)碼頭潮位相位示意圖圖4.4(b) 大潮落急流場(chǎng)圖 方案3為南口門(mén)完全封堵，且在東口門(mén)布置250m長(zhǎng)的東防波堤。它的流路變?yōu)椋簼q潮水流由新的東口門(mén)和外黃門(mén)（西口門(mén)）同時(shí)流入港內(nèi)，落潮水流由新的東口門(mén)和外黃門(mén)（西口門(mén)）同時(shí)退出。此時(shí)，進(jìn)出外黃門(mén)（西口門(mén)）的潮量繼續(xù)減少，但減小幅度??；進(jìn)出新的東口門(mén)的潮量進(jìn)一步減少，且減少幅度大。 四、防波堤工程潮流數(shù)值模擬（方案1：南防波堤，東偏北25走向，2310m長(zhǎng)）圖4.8(a) 大潮漲潮最大流速變化等值線圖工程前后流場(chǎng)變化的特點(diǎn)： （1）方案1、方案4、方案5及方案6中，東口門(mén)流速增大，漲潮平均流速變化率為+7.2%+42.9%，落

16、潮平均流速變化率為+4.1%+71.9%；方案1的流速變化最大，方案5的流速變化最小。東防波堤對(duì)通過(guò)東口門(mén)進(jìn)出灣內(nèi)的水流有阻擋作用，防波堤愈長(zhǎng)，阻擋作用愈大；方案2的東口門(mén)流速增大，而方案3的東口門(mén)大部分水域流速減小，僅在東防波堤堤頭與坎門(mén)頭之間的水域流速增大。 四、防波堤工程潮流數(shù)值模擬（方案1：南防波堤，東偏北25走向，2310m長(zhǎng)）圖4.8(b) 大潮落潮最大流速變化等值線圖 （2）方案1、方案4、方案5及方案6中，外黃門(mén)（西口門(mén)）流速減小，漲潮平均流速變化率為-10.8%-62.2%，落潮平均流速變化率為-6.0%-67.8%；方案1的流速變化最大（即流速變化率值最?。桨?的流速變

17、化最?。戳魉僮兓手底畲螅?。方案2和方案3中，外黃門(mén)（西口門(mén)）流速也減小，且流速減小的幅度較方案1大。 （3）各個(gè)方案中，外港區(qū)大部分水域流速減小，由東北端到西南端，流速減小的幅度加大；距南防波堤法向距離愈近，流速減小的幅度愈大。 （4）各個(gè)方案中，內(nèi)港區(qū)流速變化很小。 四、防波堤工程潮流數(shù)值模擬27 工程實(shí)施后的泥沙沖淤?gòu)?qiáng)度，可采用文獻(xiàn)4中的公式計(jì)算。 5.1 泥沙沖淤?gòu)?qiáng)度計(jì)算泥沙沖淤?gòu)?qiáng)度計(jì)算五、港內(nèi)泥沙沖淤分析 坎門(mén)灣為淤泥質(zhì)海灣，懸移質(zhì)造床。根據(jù)潮流數(shù)學(xué)模擬結(jié)果，采用半理論、半經(jīng)驗(yàn)公式進(jìn)行工程實(shí)施后的港內(nèi)水下地形沖淤變化計(jì)算。 )(11*2*1*SSTSnpd5.1.1 計(jì)算公式計(jì)算公

18、式5.1.2 計(jì)算參數(shù)計(jì)算參數(shù)（1）含沙量 根據(jù)2009年11月坎門(mén)漁港水文泥沙測(cè)驗(yàn)資料，按半潮平均分別統(tǒng)計(jì)水流要素v2/gh和平均含沙量，然后采用最小二乘法進(jìn)行回歸分析，得到相關(guān)關(guān)系如下： 3962.02*629.3ghvS2995.02*245.1ghvS漲潮流落潮流（2）泥沙沉速和干容重 =0.00037m/s d=723.28kg/m3五、港內(nèi)泥沙沖淤分析（方案5：南防波堤堤根位于黃門(mén)山，東偏北25走向，1473m長(zhǎng)）圖5.6 工程實(shí)施后泥沙沖淤?gòu)?qiáng)度分布圖 S05666 0.00.10.00.00.00.1-0.2-0.20.5-0.2-0.20.1（方案1：南防波堤，東偏北25走向，

19、2310m長(zhǎng)）0.10.00.10.30.5-0.20.0-0.20.00.10.30.50.7-0.20.50.3-0.2-0.20.1圖5.2 工程實(shí)施后泥沙沖淤?gòu)?qiáng)度分布圖 S01666 5.1.3 計(jì)算結(jié)果計(jì)算結(jié)果五、港內(nèi)泥沙沖淤分析圖5.4 工程實(shí)施后泥沙沖淤?gòu)?qiáng)度分布圖 （方案3：南防波堤，東偏北25走向，2310m長(zhǎng)；東防波堤，南北走向，250m長(zhǎng)）S03666 0.70.50.30.10.0-0.20.00.0-0.20.70.50.30.10.70.00.10.3 方案5的泥沙淤積最少，港內(nèi)平均淤積強(qiáng)度為0.05m/a；方案3的泥沙淤積最多，港內(nèi)平均淤積強(qiáng)度為0.45m/a；其他

20、方案的港內(nèi)平均淤積強(qiáng)度在0.260.39m/a之間。 5.2 最終沖淤厚度計(jì)算最終沖淤厚度計(jì)算5.2.1 計(jì)算公式計(jì)算公式 假定工程前、后水流挾沙能力不變，可得新沖淤平衡狀態(tài)下的水深： 3/2123/12113)()(qqSShh 床面最終沖淤厚度： 31hhz五、港內(nèi)泥沙沖淤分析（方案1：南防波堤，東偏北25走向，2310m長(zhǎng)）圖5.8 工程后沖淤平衡時(shí)的水下地形圖 (m)1985國(guó)家高程基準(zhǔn)5.2.2 計(jì)算結(jié)果計(jì)算結(jié)果 （1）方案3的東口門(mén)，泥沙淤積，但最終淤積厚度較?。?.46m）；其余方案的東口門(mén)，海床沖刷；方案2和方案3的東防波堤堤頭前水域海床沖刷分別達(dá)到5.01m和4.37m。 （2）各個(gè)方案的外黃門(mén)（西口門(mén)），泥沙淤積；方案2的最終淤積厚度最大，達(dá)到7.35 10.70m，但水深仍超過(guò)10.0m（85國(guó)家高程基準(zhǔn)，下同）；方案5的最終淤積厚度最?。?.081.74m）。 五、港內(nèi)泥沙沖淤分析 （3）各個(gè)方案的外港區(qū)大部分水域，泥沙淤積，由東北端到西南端，最終淤積厚度逐漸增大；距南防波堤法向距離愈近，最終淤積厚度愈大；在南防波堤堤前水域，方案3的最終淤積厚度最大，達(dá)到5.09 6.93m；但南防波堤堤前水深仍超過(guò)5.0m。 （4）內(nèi)港區(qū)、二期防波堤所形成的口門(mén)及其南側(cè)小范圍水域，海床基本不沖不淤。 五、