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

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

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

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

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

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

6、()(2222yUxUHNfHVxZHgHyHVUxHUUtHUxbbxsx)()()()()(2222yVxVHNfHUyZHgHyHVVxHUVtHVybbysy3.2 計(jì)算模式計(jì)算模式 采用有限體積法離散上述基本方程；通過離散得到的線性代數(shù)方程組，用高斯消去法求解。 11三、坎門漁港潮流數(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三、坎門漁港潮流數(shù)學(xué)模型+ + + + + + + + + + + + 長海大青蘇爺島大連秦皇島連云港大戢山濟(jì)州島福江寶島儀間大陳島平壇臺灣朝鮮半島九州圖3.2 中國東部海區(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 模型的水邊界采用潮位控制，潮位過程由東中國海潮波數(shù)學(xué)模型提供。 對于隨潮位變化的淺灘動邊界，采用網(wǎng)格干、濕判斷法進(jìn)行處理。 13三、坎門漁港潮流數(shù)學(xué)模型3.6 驗(yàn)證計(jì)算驗(yàn)證計(jì)算 模型驗(yàn)證取用2009年11月3日11月10日的實(shí)測大、中、小潮水文資料，測站布設(shè)見圖2.3，有2個(gè)潮位測站（W1、W2）和6條測流垂線（S1 S6）。 潮位站測流點(diǎn)圖2.3 現(xiàn)場水文泥沙測驗(yàn)站位圖14三、坎門漁港潮流數(shù)學(xué)模型圖3.3(a) 大潮潮位驗(yàn)證曲線黃門

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

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

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

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

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

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

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

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

17、化最?。戳魉僮兓手底畲螅?。方案2和方案3中，外黃門（西口門）流速也減小，且流速減小的幅度較方案1大。 （3）各個(gè)方案中，外港區(qū)大部分水域流速減小，由東北端到西南端，流速減小的幅度加大；距南防波堤法向距離愈近，流速減小的幅度愈大。 （4）各個(gè)方案中，內(nèi)港區(qū)流速變化很小。 四、防波堤工程潮流數(shù)值模擬27 工程實(shí)施后的泥沙沖淤強(qiáng)度，可采用文獻(xiàn)4中的公式計(jì)算。 5.1 泥沙沖淤強(qiáng)度計(jì)算泥沙沖淤強(qiáng)度計(jì)算五、港內(nèi)泥沙沖淤分析 坎門灣為淤泥質(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月坎門漁港水文泥沙測驗(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：南防波堤堤根位于黃門山，東偏北25走向，1473m長）圖5.6 工程實(shí)施后泥沙沖淤強(qiáng)度分布圖 S05666 0.00.10.00.00.00.1-0.2-0.20.5-0.2-0.20.1（方案1：南防波堤，東偏北25走向，

19、2310m長）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í)施后泥沙沖淤強(qiáng)度分布圖 S01666 5.1.3 計(jì)算結(jié)果計(jì)算結(jié)果五、港內(nèi)泥沙沖淤分析圖5.4 工程實(shí)施后泥沙沖淤強(qiáng)度分布圖 （方案3：南防波堤，東偏北25走向，2310m長；東防波堤，南北走向，250m長）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長）圖5.8 工程后沖淤平衡時(shí)的水下地形圖 (m)1985國家高程基準(zhǔn)5.2.2 計(jì)算結(jié)果計(jì)算結(jié)果 （1）方案3的東口門，泥沙淤積，但最終淤積厚度較?。?.46m）；其余方案的東口門，海床沖刷；方案2和方案3的東防波堤堤頭前水域海床沖刷分別達(dá)到5.01m和4.37m。 （2）各個(gè)方案的外黃門（西口門），泥沙淤積；方案2的最終淤積厚度最大，達(dá)到7.35 10.70m，但水深仍超過10.0m（85國家高程基準(zhǔn)，下同）；方案5的最終淤積厚度最?。?.081.74m）。 五、港內(nèi)泥沙沖淤分析 （3）各個(gè)方案的外港區(qū)大部分水域，泥沙淤積，由東北端到西南端，最終淤積厚度逐漸增大；距南防波堤法向距離愈近，最終淤積厚度愈大；在南防波堤堤前水域，方案3的最終淤積厚度最大，達(dá)到5.09 6.93m；但南防波堤堤前水深仍超過5.0m。 （4）內(nèi)港區(qū)、二期防波堤所形成的口門及其南側(cè)小范圍水域，海床基本不沖不淤。 五、