武賀，趙世明，張松，王鑫，馬治忠

（國家海洋技術(shù)中心，天津 300111）

老鐵山水道潮流能初步估算

武賀，趙世明，張松，王鑫，馬治忠

（國家海洋技術(shù)中心，天津 300111）

采用ECOMSED海洋數(shù)值模型較好地模擬了渤海的潮流運動狀況，估算了老鐵山水道的潮流能的理論蘊藏量，并采用FLUX方法對該水道的技術(shù)可開發(fā)量進行了評估。結(jié)果表明，老鐵山水道北側(cè)近岸海域最大可能流速約2.5 m/s，平均功率密度超過500 W/m2，大潮平均最大功率密度為3 700 W/m2，小潮平均最大功率密度為1 400 W/m2，水道潮流能資源理論蘊藏量為26.9×104kW，技術(shù)可開發(fā)量為2.69×104kW。

ECOMSED；老鐵山水道；潮流能

老鐵山水道位于渤海、黃海的分界線上，北起老鐵山燈塔，南至北隍城，全長約為13 km，水道北側(cè)水深流急，平均水深約為60 m，水道南側(cè)水深較淺，平均水深約為30 m，該水道實測最大流速超過2.3 m·s-1，潮流能源蘊藏量豐富。我國早期的潮流能資源評估工作[1]曾利用鄭志南[2]提出的評估方法計算了該水道蘊藏量并得到了初步結(jié)論，這對于潮流數(shù)據(jù)匱乏的評估研究初期具有重要的開創(chuàng)性意義。然而，由于受到水深地形、海水密度等因素的影響，水道截面的水平方向和垂直方向上的流速不可能相同，因此以單點流速代表整個水道的截面流速的做法值得商榷，加之能量與流速存在著3次方的關(guān)系，使得評估結(jié)果有可能存在不小的誤差。隨著計算機和數(shù)值計算技術(shù)的發(fā)展，已經(jīng)可以實現(xiàn)準(zhǔn)確的潮波模擬和潮流預(yù)報，獲取豐富的四維（三維空間和時間）潮流數(shù)據(jù)，從而徹底解決潮流能資源評估工作中的數(shù)據(jù)短缺問題，事實上，國外不少學(xué)者已經(jīng)成功地采用此類方法進行了一些水道的潮流資源評估試驗[3-8]，相比之下，國內(nèi)的潮流數(shù)值模擬研究雖然也較為成熟[10-12]，卻很少應(yīng)用于潮流資源評估工作[9]。因此，本文將利用國際上廣泛使用的ECOMSED海洋環(huán)流數(shù)值模式對渤海、北黃海進行潮波模擬計算，在較好地完成潮流模擬的基礎(chǔ)上，采用FLUX計算方法分析計算老鐵山水道的潮流能蘊藏量和可開發(fā)量，該項研究工作對進一步了解老鐵山水道潮流能資源的時空分布和日后的潮流能資源開發(fā)利用具有重要的實踐意義。

2 數(shù)值模式及其設(shè)置

ECOMSED海洋環(huán)流數(shù)值模式是一個模擬水動力、波浪和沉積物輸運的三維數(shù)值模型，該模式具有熱力學(xué)方程完整、擬合水深地形較好、收斂速度快等特點，被廣泛地應(yīng)用于近海動力科學(xué)研究和海洋工程。本文考慮M2、S2、K1、O1 4個分潮。模式計算區(qū)域包括渤海和北黃海，水深數(shù)據(jù)采用1︰100萬中國海圖，開邊界潮汐調(diào)和常數(shù)采用T/P反演數(shù)據(jù)。水平方向采用seagrid軟件生成的正交曲線網(wǎng)格（見圖1），并在黃渤海分界線海域進行網(wǎng)格加密以便更好地了解老鐵山水道的潮流運動狀況，其最小分辨率為1 km×1 km，垂向上均勻分5層。

圖 1 模擬海域網(wǎng)格設(shè)置及水深分布（m）Fig.1 Depth distribution and the grid settings of the simulative area (m)

3 數(shù)值模擬結(jié)果

潮位方面，本文分別選取具有代表性鲅魚圈、秦皇島、塘沽、龍口和蓬萊5個驗潮站潮位數(shù)據(jù)進行調(diào)和分析，并與計算結(jié)果進行比較。表1顯示，渤海灣內(nèi)5個海洋觀測站的實測潮位數(shù)據(jù)與計算結(jié)果吻合良好， 此外，4個分潮（M2、S2、K1、O1）的同潮圖（圖2a-圖2d）也很好地再現(xiàn)了渤海海域各潮波系統(tǒng)的駐波性質(zhì)。

潮流方面，采用2009年3月1個ADCP定點站位（121°05.8′E，38°41.8′N）連續(xù)25 h的實測資料對中海域模擬結(jié)果進行驗證（見圖3a-3b），由各站位流速、流向?qū)Ρ葓D可知，計算值與實測結(jié)果吻合較好。

綜上所述，本文較好地模擬了渤海的潮波運動狀況，成功地再現(xiàn)了渤海海域各主要分潮潮波系統(tǒng)的駐波性質(zhì)，為進行準(zhǔn)確的潮流能計算奠定了基礎(chǔ)。

表 1 渤海沿岸5個驗潮站調(diào)和常數(shù)實測與計算誤差Tab.1 Errors of harmonic constants between measurements and calculations of five tide stations around the Bohai Sea

圖 2a M2分潮同潮圖Fig.2a Co-tidal line of M2 tide

圖 2b S2分潮同潮圖 Fig.2b Co-tidal line of S2 tide

圖2c K1分潮同潮圖Fig.2c Co-tidal line of K1 tide

圖2d O1分潮同潮圖 Fig.2d Co-tidal line of O1 tide

圖 3a 實測與預(yù)報流速曲線圖(3月10日-3月11日)Fig.3a Speed of the tidal current between measurement forecast (Mar 10 to Mar 11)

圖 3b 實測與預(yù)報流向曲線圖(3月10日-3月11日) Fig.3b Direction of the tidal current between measurement and and forecast (Mar 10 to Mar 11)

4 老鐵山水道潮流能分析

4.1 估算方法

潮流能估算的基本方法有兩類[9]，一類以能通量為基礎(chǔ)，一類基于動力分析。前者包含了FLUX方法、FARM方法等，后者則是Garrett[3]提出的一種新方法，應(yīng)用上，F(xiàn)ARM方法最為廣泛，F(xiàn)LUX方法次之，Garrett方法還很少，由于本文僅對老鐵山水道的潮流能蘊藏量及技術(shù)可開發(fā)量進行初步評估，并未涉及潮流能開發(fā)裝置類型及轉(zhuǎn)換效率參數(shù)等，因此，本文采用使用廣泛的FLUX方法。平均功率的理論公式及實際計算公式分別為：

式中：

P是指周期T內(nèi)潮流能的平均功率(W)；

t初始時刻，T為評估周期，一般取一個朔望月即可(s)；

L水道寬度(m)；

H水深(m)；

ρ為海水密度，取1025kg·m-3。

實際計算公式

SIF為有效影響因子。文獻[15]顯示，全部潮流能中可開發(fā)部分的比例與地點有關(guān)。在水道或海峽、海岬、海湖、共振河口等各海域中SIF的取值各有不同。由于老鐵山水道屬于海峽范疇，因此本文取10%。

4.2 最大可能流速分布

采用文獻[14]中所記載的最大可能流速計算公式，根據(jù)潮流類型的不同，分別計算各網(wǎng)格點最大可能流速。圖4顯示，最大可能流速較強海域主要為渤海海峽、黃河口附近、遼東灣兩處海域以及河北曹妃甸附近，其最大可能流速都超過1 m·s-1；強流速區(qū)位于老鐵山水道北側(cè)，最大可能流速值大都超過1.5 m·s-1，其中有2個網(wǎng)格點超過2.5 m·s-1，其中一點的最大功率密度可達12.5 kW·m-2，且離岸較近，具有較大的潮流能開發(fā)潛力。此外，北隍城至蓬萊的廟島群島諸水道最大可能流速也較強，大都超過1.5 m·s-1。

圖 4 渤海海域最大可能潮流速度分布（m·s-1）Fig.4 Most possibility speed of tidal current in Bohai Sea (m·s-1)

4.3 潮流能蘊藏量

圖5、圖6、圖7分別為老鐵山水道平均能流密度分布、大潮平均最大功率密度分布和小潮平均最大功率密度分布，圖中顯示，此3個重要參數(shù)分布基本相同，且與最大可能流速分布相似。其中有2個網(wǎng)格點的平均能流密度超過500 W·m-2，大潮最大功率密度可達3 700 W·m-2，小潮最大功率為密度為1 400 W·m-2，涉及海域面積約為10 km2。過最大平均功率密度點，繪制渤海海峽潮流漲急時刻的法向線（圖5中虛線），以此虛線為截面，采用公式（2），計算出老鐵山水道潮流能理論蘊藏量為26.9 ×104kW，技術(shù)可開發(fā)量為2.69×104kW。

圖 5 老鐵山水道平均能流密度水平分布（W·m-2）Fig.5 Distribution of average tidal energy density at LaoTieshan channel（W·m-2）

圖 6 老鐵山水道大潮平均最大能流密度水平分布（W·m-2）Fig.6 Distribution of average tidal energy density during the spring tide period at LaoTieshan channel（W·m-2）

圖 7 老鐵山水道小潮平均最大能流密度水平分布（W·m-2）Fig.7 Distribution of average tidal energy density during the neap period at LaoTieshan channel（W·m-2）

眾所周知，潮流能隨著潮流流速的變化也進行著周期性的變化。由于海峽（水道）海域的年平均潮流能與月平均潮流能差別不大[9],本文僅給出老鐵山水道一個朔望月的潮流能時間分布（圖8）。經(jīng)統(tǒng)計計算，該水道小于50×104kW的時間比例約為81.4 %，介于50×104kW和100×104kW的時間比例約為13.4 %，大于100×104kW的時間比例約為6 %。

圖 8 老鐵山水道潮流能時間變化曲線（一個朔望月）Fig.8 Change of the tidal energy during one month at LaoTieshan channel

5 結(jié) 論

本文采用ECOMSED海洋數(shù)值模型模擬了渤海的潮流運動狀況，繪制了K1、O1、M2、S2 4個分潮的同潮圖，并與實測潮位、潮流數(shù)據(jù)進行了對比分析。在此基礎(chǔ)上估算了老鐵山水道潮流能的理論蘊藏量，并采用FLUX方法對該水道的技術(shù)可開發(fā)量進行了評估。結(jié)果表明，老鐵山水道北側(cè)近岸海域最大可能流速大于2.5 m/s，平均功率密度超過500 W·m-2，大潮平均最大功率密度為3 700 W·m-2，小潮最大功率密度為1 400 W·m-2，水道潮流能資源理論蘊藏量為26.9×104kW，技術(shù)可開發(fā)量為2.69 ×104kW，潮流能資源隨時間不斷變化，在一個朔望月的周期內(nèi)，能量大于50×104kW的時間比例約為81.4 %，介于50×104kW和100×104kW之間的時間比例約為13.4 %，大于100×104kW的時間比例約為6 %。

[1]王傳崑, 陸德超.中國沿海農(nóng)村海洋能資源區(qū)劃 [M].1989: 1-127.

[2]鄭志南.海洋潮流能的估算 [J].海洋通報, 1987, 6(4): 70-75.

[3]Garrett C, Cummins P.Generating power from tidal currents [J].Waterway Port Coastal Ocean Eng, 2004, 130: 114-118.

[4]Karsten R H, McMillan J M, Lickley M J, et al.Assessment of tidal current energy in the Minas Passage, Bay of Fundy [J].Power and Energy, 2008, 22(A): 493-507.

[5]David A B.The tidal stream-energy resource in Passamaquoddy-Cobscook Bay: A fresh look at an old story [J].Renewable Energy, 2006, 31(14): 2284-2295.

[6]Ian G B, Scott J C.How much energy can extracted from moving water with a free surface: A question of importance in the field of tidal current energy [J].Renewable Energy, 2007, 32: 1961-1966.

[7]Sutherland G, Foreman M, Garrett C.Tidal current energy assessment for Johnstrone strait Vancouver Island [J].Power and Energy, 2007, 221: 147-157.

[8]Carballo R, Iglesias G, Castro A.Numerical model evaluation of tidal stream energy resources in the Ria de Muros(NW Spain) [J].Renerable Energy, 2009, 34: 1517-1524.

[9]呂新剛, 喬方利, 趙昌, 等.海洋潮流能資源的數(shù)值估算——以膠州灣口為例 [J].太陽能學(xué)報, 2010, 31(2): 137-143.

[10]趙保仁, 方國洪.渤、黃、東海潮汐潮流的數(shù)值模擬 [J].海洋學(xué)報, 1994, 16(5): 1-10.

[11]王凱, 馮士筰.渤海、黃海、東海M2潮汐潮流的三維數(shù)值模擬[J].海洋學(xué)報, 1999, 21(4): 1-13.

[12]萬振文, 喬方利, 袁業(yè)立.渤、黃、東海三維潮波運動數(shù)值模擬[J].海洋與湖沼, 1998, (29)6: 611-616.

[13]呂新剛, 喬方利.海洋潮流能資源估算方法研究進展 [J].海洋科學(xué)進展, 2008, 26(1): 98-108.

[14]中華人民共和國行業(yè)標(biāo)準(zhǔn) [S].海港水文規(guī)范, 1999.

[15]Black & Veatch Consulting, Ltd.Phase II, UK tidal stream energy resource assessment [R].Marine Energy Challenge Report, 2005, No.107799/D/2200/03.

Preliminary assessment of tidal energy in Lao Tieshan channel

WU He, ZHAO Shi-ming, ZHANG Song, WANG Xin, MA Zhi-zhong

（National Ocean Technology Center, Tianjin 300111, China）

ECOMSED ocean model is applied to simulate the tidal wave system on Bohai sea, and the error between investigation and simulation can be acceptable.Based on this result, we assess the tidal energy of LaoTieshan channel.By the method of FLUX, a conclusion can be extracted like this: the probable velocity of maximum of north of this channel can reach to 2.5 m/s, and the average density of power is above 500W/m2, meanwhile, the maximal densities of power of spring tide period and neap tide period are 3 700 W/m2and 1 400 W/m2respectively, and theoretical tidal energy and the available tidal energy of LaoTieShan Channel are 26.9×104kW and 2.69×104kW.

ECOMSED; LaoTieshan channel; tidal energy

P743.1

A

1001-6932(2011)03-0310-05

2010-11-19；收修改稿日期：2011-01-20

我國近海海洋綜合調(diào)查評價專項 (908-01-NY, 908-02-05-01, 908-ZC-I-21 )。

武賀 ( 1981－ )，男，碩士，主要從事近海環(huán)境動力學(xué)研究。電子郵箱：wh-crane@sohu.com。