張俊彪，繆 斌 ，許雪峰 *,3,楊萬康 ，張 峰 ，潘 沖 ，施偉勇

（1.國家海洋局第二海洋研究所，浙江 杭州 310012；2.國家海洋局 閩東海洋環(huán)境監(jiān)測中心站，福建 寧德 352000；3.浙江大學(xué)，浙江 杭州 310006)

海灣內(nèi)外的動態(tài)水位差模擬研究
——一種全新的潮汐能開發(fā)方式的探索

張俊彪1，繆 斌2，許雪峰*1,3,楊萬康1，張 峰1，潘 沖1，施偉勇1

（1.國家海洋局第二海洋研究所，浙江 杭州 310012；2.國家海洋局 閩東海洋環(huán)境監(jiān)測中心站，福建 寧德 352000；3.浙江大學(xué)，浙江 杭州 310006)

首次提出了一種全新的海灣潮汐能開發(fā)方式——利用海灣內(nèi)外的動態(tài)水位差進行潮汐能發(fā)電。原理是利用潮波到達海灣底部和海灣外部的相位不一致性，形成了海灣內(nèi)外動態(tài)水位差，可作為潮汐能發(fā)電水頭。并構(gòu)建了一個理想潮汐能海灣進行水動力數(shù)值模擬研究，模擬結(jié)果顯示該新型的潮汐能是存在的，且相比傳統(tǒng)潮汐能利用方式對海灣水動力環(huán)境和通航影響更小。

數(shù)值模擬；潮汐能；海灣；動態(tài)

潮汐能的最早利用方式是11世紀出現(xiàn)的潮汐磨坊。潮汐發(fā)電的實際應(yīng)用開始于1912年德國胡蘇姆興建的一座小型潮汐電站。1966年法國建成的朗斯潮汐電站，裝機容量為24萬kW，年均發(fā)電量為5.44億kW·h，是目前最大的潮汐電站[1-2]。我國潮汐能發(fā)電始于20世紀50年代后期[3]，迄今建成潮汐電站8座，總裝6 120 kW，其中最大的是浙江江廈潮汐試驗電站，為3 900 kW。

但潮汐能在商業(yè)上未能取得公認的常規(guī)能源地位。朗斯電站總結(jié)為“在技術(shù)上可行，經(jīng)濟上不合算”的論點仍是其后許多能源決策者對潮汐能評估的依據(jù)。2 0余年來世界上大型潮汐電站仍處于前期論證和試驗研究階段，進展較緩，迄今未見規(guī)模大于朗斯的潮汐電站開工建造[4]。

這主要是由于傳統(tǒng)的潮汐電站存在以下的問題：（1）需要在海灣內(nèi)建設(shè)大壩，投資巨大（圖1）；（2）受到庫容的限制，電站的裝機總量有限，如我國最大的江廈潮汐實驗電站[5]裝機總量僅為3 900 kW；（3）建壩對海灣進行截流，會對海灣和庫區(qū)內(nèi)造成較大淤積，影響電站使用壽命和增加維護成本；（4）海灣截流會對魚類洄游、海灣水交換造成影響，破壞海灣生態(tài)；（5）需要納潮發(fā)電，導(dǎo)致實際發(fā)電時間較短。

由于以上缺陷，因此目前潮汐能電站遭遇了發(fā)展的瓶頸，中國僅存江廈潮汐實驗電站正常發(fā)電。國際上對潮汐能資源的開發(fā)也遠遠落后于風(fēng)能、太陽能等新能源。這對于蘊藏豐富的潮汐能資源是極大的浪費。

本文提出的全新潮汐能發(fā)電方法不同于傳統(tǒng)的潮汐電站，是潮汐能開發(fā)利用的新思路和新探索。

1 傳統(tǒng)潮汐發(fā)電方法及全新方法的提出

1.1 傳統(tǒng)潮汐發(fā)電方法

傳統(tǒng)潮汐電站主要由5個基本部分組成：潮汐水庫、堤壩、閘門水道、發(fā)電機組和廠房、輸電設(shè)施等[6]。

傳統(tǒng)潮汐電站主要有兩類[7]：單庫型（單向型、雙向型）和雙庫型。由于雙庫型潮汐電站發(fā)電效率較低且投資較高，因此目前的潮汐電站主要為單庫型。

圖1 傳統(tǒng)單庫型潮汐電站發(fā)電工況的平面簡面

圖1為單庫型潮汐電站發(fā)電工況的平面簡面。傳統(tǒng)的單庫單向潮汐電站工作原理如下：高潮時段，打開水閘納水，對庫區(qū)進行蓄水，當庫區(qū)水位達到最高時關(guān)閉水閘。等到低潮位時段，水閘保持關(guān)閉狀態(tài)，庫區(qū)內(nèi)與庫區(qū)外達到一定水位差時，打開電站機組進行發(fā)電。而單庫雙向潮汐電站與單庫單向型的主要區(qū)別是采用雙向發(fā)電機組，在漲潮納水階段仍然可以進行發(fā)電。傳統(tǒng)潮汐能電站需在海灣內(nèi)建設(shè)截流大壩，因此對生態(tài)環(huán)境影響較大。

圖2為傳統(tǒng)潮汐發(fā)電過程示意圖[7]，由圖可見在水庫水位在高水位和低水位分別有較長時間是等待期（蓄水等待大壩兩側(cè)形成較大水位差），在此期間發(fā)電機不能發(fā)電。從江廈電站運行情況看，裝機年利用小時為2 300～2 500 h之間，發(fā)電小時數(shù)為5 000～6 000 h。

圖2 傳統(tǒng)單庫型潮汐發(fā)電示意圖[7]

1.2 全新的潮汐能發(fā)電方法

發(fā)電原理：潮波由外海向海岸傳播，高潮等時線到達海灣外要早于到達海灣內(nèi)，導(dǎo)致海灣外與海灣內(nèi)的潮汐相位差異，灣內(nèi)和灣外的水位存在動態(tài)的水位差。因此選擇合適的地點，開挖流道，連通灣內(nèi)與灣外的水體，并在流道上安裝新型雙向潮汐發(fā)電燈泡貫流式水輪發(fā)電機組[2]，就可以利用灣內(nèi)和灣外的動態(tài)水位差進行雙向發(fā)電。

方法介紹：潮波從外海傳到海灣底部需要一定的時間，如圖3所示，結(jié)合當?shù)氐牡匦蝺?yōu)勢，電站位置設(shè)置在海灣底部，開挖流道即可打通海灣底部與外海，在流道上安裝發(fā)電機組即完成電站主體工程。發(fā)電各階段如下：

發(fā)電階段一（圖3）：外海潮波向海岸和海灣底部傳播，當中低潮位線到達海灣內(nèi)部（即電站的內(nèi)側(cè)）時，中高潮位線已經(jīng)到達海灣外側(cè)（即電站外側(cè)）。此時段，電站內(nèi)側(cè)是中低潮位，電站外側(cè)是中高潮位，存在水位差，可開啟發(fā)電機組進行發(fā)電。此階段是從灣外向灣內(nèi)排水發(fā)電。

發(fā)電階段二（圖4）：當中高潮位線到達海灣內(nèi)部（即電站的內(nèi)側(cè)）時，中低潮位線到達海灣外側(cè)（即電站外側(cè)）。此時段，電站內(nèi)側(cè)是中高潮位，電站外側(cè)是中低潮位，存在水位差，發(fā)電機組可向外排水發(fā)電。

發(fā)電階段一和階段二是交替循環(huán)工作的，以保證在大多數(shù)的時間內(nèi)電站保持發(fā)電。

圖3 新型潮汐能發(fā)電示意圖（階段一）

圖4 新型潮汐能發(fā)電示意圖（階段二）

1.3 新方法與傳統(tǒng)方法相比的優(yōu)點

根據(jù)上述的新型潮汐能發(fā)電的方法，有以下優(yōu)點：（1）不在海灣內(nèi)建設(shè)大壩，不建設(shè)封閉的水庫，僅需在海灣底部開挖通道，因此在發(fā)電期間不會減少海灣的納潮量，且有利于海灣底部水交換；（2）與海灣的通航和其他資源利用的矛盾較小；（3）海灣內(nèi)的水體均可作為發(fā)電水體，而發(fā)電量主要受流量及內(nèi)外水位差的限制，因此只要開挖合適的流道、設(shè)計合適的流量便可獲得較大的裝機容量。

2 理想潮汐能海灣的水動力模型

為了驗證該新型的潮汐能開發(fā)方式是否可行，本文構(gòu)建了理想的潮汐能海灣水動力模型，進行數(shù)值模擬試驗。

2.1 理想海灣網(wǎng)格構(gòu)建

構(gòu)建如圖5所示的理想海灣，水深均為10 m，海灣尺度為南北跨度95 km，東西寬60 km。為了便于計算和給定邊界條件，本模型取均勻的正方形網(wǎng)格，網(wǎng)格邊長取3.5 km。采用較大且正規(guī)的計算網(wǎng)格具有計算速度較快、模型更穩(wěn)定的特點，且也能滿足本研究的定性定量的計算要求。

為了確定海灣底部與灣外的動態(tài)水位差，分析潮汐能發(fā)電的可行性，在模型中設(shè)置了兩個點（out、in），如圖5所示。記錄該兩點的全程模擬數(shù)據(jù)，in、out兩點的水上距離最近約150 km，而in點到灣口的水上距離約100 km。

圖5 理想潮汐能海灣計算網(wǎng)格

2.2 控制方程

控制方程采用N-S的淺水方程[8-9]：

式中：ζ為潮位；g為重力加速度；Aζ，Aη為渦動粘滯系數(shù)；u，v為 ξ，η 方向上的流速分量；h 為水深；H=h+ζ為總水深；t為時間；C 為謝才系數(shù)，C=1/n·H1/6，n 為糙率系數(shù)；f為柯氏系數(shù)，f=2ωsinφ，ω為地轉(zhuǎn)角速度，φ 為緯度。

對水動力方程的離散引用傳統(tǒng)的ADI法數(shù)值求解，計算時間步長為1 min。

2.3 邊界條件

（1）邊界條件：開邊界采用水位控制，邊界水位采用潮位預(yù)報的方法得到[10-11]：

式中：A0為平均海面；Fi，(v0+u)i為天文要素；Hi，gi為調(diào)和常數(shù)。為了簡化模型，且更清楚地顯現(xiàn)理想海灣的動態(tài)水位差，這里調(diào)和常數(shù)僅選用2個分潮M2和S2，振幅取值分別為1.8 m和0.5 m，相位取值分別為260°和290°，大致和中國東海強潮區(qū)域調(diào)和常數(shù)相符，計算潮型取大潮。

（2）初始條件：

2.4 參數(shù)選取

柯氏力 f=2ω*sin(29.1°)；粘滯系數(shù) A 取 20～60 m2/s；水容重 ρ=1 020 kg/m3；糙率系數(shù) n=0.010。

3 理想潮汐能海灣的模擬結(jié)果分析

3.1 潮流場分析

海域潮流漲落潮流場分布見圖6。計算結(jié)果顯示，理想海灣內(nèi)漲潮流速最大為2.8 m/s，落潮流速最大為2.2 m/s；平均漲潮流速為1.5 m/s,平均落潮流速為1.3 m/s；由此可見，流速強度和我國的強潮海灣大致相當。而流速較強的位置一般出現(xiàn)在海灣岸線的岬角前沿，海灣底部的流速相對較弱，水交換強度也較弱，這和我國大多數(shù)海灣的水動力特征一致。

圖6 圍墾工程實施前雙盤涂附近漲、落急流矢圖

3.2 海灣內(nèi)外動態(tài)水位差發(fā)電可行性分析

圖7 動態(tài)水位差平面分布圖

圖8 電站內(nèi)外兩側(cè)水位變化過程曲線

圖9 電站內(nèi)外兩側(cè)的動態(tài)水位差過程曲線

模型計算結(jié)果中，每2 h輸出水位平面分布，如圖7所示。可見，海灣內(nèi)外動態(tài)水位差最大的位置出現(xiàn)海灣底部附近，暨圖示的電站位置。

模型試驗中，在海灣內(nèi)外（暨電站的內(nèi)外兩側(cè)）分別設(shè)置兩個點（圖5），記錄水位變化過程曲線。海灣底部（電站內(nèi)側(cè)）和海灣外（電站外側(cè)）的水位變化過程曲線見圖8所示，可見電站內(nèi)外兩側(cè)的潮時相差約4.5 h。電站外側(cè)水位表現(xiàn)較為正規(guī)；而電站內(nèi)側(cè)低水位時間長于高水位時段，水位抬升過程較快，略有涌潮特征。因此，電站內(nèi)外兩側(cè)的潮時和潮汐特征均表現(xiàn)不一致，使得大多數(shù)時間內(nèi)兩側(cè)存在動態(tài)水位差。

圖9所示為電站內(nèi)外兩側(cè)的動態(tài)水位過程曲線。由圖可見，電站兩側(cè)的動態(tài)水位差最大可達6.0 m，幾乎接近最大潮差；兩側(cè)平均水位差在3.0 m左右，因此可達到較好的發(fā)電效果。從過程曲線來看，海灣內(nèi)側(cè)向外側(cè)排水（正向）發(fā)電的時間更長，而外側(cè)向內(nèi)側(cè)排水（逆向）發(fā)電的水頭更大。

從可發(fā)電時間來看，兩側(cè)水位差在90%的時段內(nèi)大于1.2 m，而現(xiàn)有的潮汐能發(fā)電技術(shù)水位差大于0.8 m即可進行發(fā)電，因此該電站可在大潮90%時段內(nèi)進行發(fā)電，要高于傳統(tǒng)潮汐電站一般大潮的發(fā)電時間。

4 小結(jié)

本文提出了海灣潮汐能開發(fā)的新思路，并構(gòu)建了理想海灣水動力數(shù)值模型進行試驗。試驗結(jié)果顯示，在理想的大潮條件下，采用該新型的潮汐能發(fā)電方式可使得潮汐電站發(fā)電時間提高到90%左右，最大發(fā)電水頭幾乎接近最大潮差，相對傳統(tǒng)潮汐電站提高了潮汐能的利用時間。

研究表明，利用海灣內(nèi)外的動態(tài)水位差進行潮汐能發(fā)電的新思路，在理論上是可行的。與傳統(tǒng)潮汐能電站的發(fā)電模式比較，該全新潮汐能開發(fā)方式還有以下優(yōu)點：不在海灣內(nèi)建設(shè)大壩，因此對海灣的通航等其他資源利用的矛盾較小，也可能對海灣水動力環(huán)境影響更小。

而中國沿海有眾多海灣形態(tài)都類似于文中的理想海灣，如浙江樂清灣和福建三都澳，因此該潮汐能的全新開發(fā)方式可能有較好的發(fā)展前景。

[1]劉全根.世界海洋能開發(fā)利用狀況及發(fā)展趨勢[J].能源工程，1999，2∶5-8.

[2]朱成章.關(guān)于我國潮汐能資源和潮汐電站建設(shè)情況[J].新能源，1995，17(2)∶5-8.

[3]余志.海洋能利用技術(shù)進展與展望[J].太陽能學(xué)報，1999，特刊∶214-226.

[4]郭成濤.潮汐能利用的新概念[J].海洋學(xué)報，1994，16(1)∶142-148.

[5]張發(fā)華.綜合開發(fā)我國潮汐能的探討[J].水力發(fā)電學(xué)報，1996，3∶33-42.

[6]徐錫華，阮世銳.潮汐電站設(shè)計導(dǎo)則[J].國外水電技術(shù)，1991，6(2)∶1-2.

[7]李書恒，郭偉，朱大奎.潮汐發(fā)電技術(shù)的現(xiàn)狀與前景[J].海洋科學(xué)，2006，12∶82-86.

[8]李孟國，曹祖德.海岸河口潮流數(shù)值模擬的研究與進展[J].海洋學(xué)報，1999，21(1)∶61-69.

[9]馮士笮，孫文心，等.物理海洋數(shù)值計算[M].河南科學(xué)技術(shù)出版社，1990.

[10]李孟國，王正林.甌江口潮流數(shù)值模擬[J].長江科學(xué)院院報,2002，19(2)∶8-17.

[11]湯立群，金忠青.正交曲線坐標下河口二維潮流過程計算[J].水動力學(xué)研究與進展，2003，18(2)∶31-41.

Numerical Research on Dynamic Water-head Between Sea and Gulf——A New Exploration on Development Way of Tidal Energy

ZHANG Jun-biao1,MIU Bin2,XU Xue-feng1,3,YANG Wan-kang1,ZHANG Feng1,PAN Chong1,SHI Wei-yong1
(1.The Second Institute of Oceanography,SOA,Hangzhou Zhejiang 310012,China;2.The Marine Environment Monitoring Center Station of Mindong,SOA,Ningde Fujian 352000,China;3.Zhejiang University,Hangzhou Zhejiang 310006,China)

A new way of the Gulf tidal energy development was proposed using the dynamic water-head inside and outside the Gulf for tidal power generation.Tide inconsistency was used when it reaches the bottom of the Gulf and outside the Gulf,leading to delay in the tide on Gulf bottom.A dynamic water-level difference was formed between inside and outside of the Gulf,which can be used for tidal power water-head.An ideal tidal energy gulf was built up for the hydrodynamic numerical simulation.Results show that the new tidal energy is existent.Compared to the traditional tidal energy,its impaction is smaller to the hydrodynamic environment and navigation.

numerical simulation;tidal energy;gulf;dynamic

P743.3，P753

A

1003-2029（2012）04-0083-05

2012-05-11

海洋可再生能源專項資金資助項目（GHME2011CX01，GHME2011ZC05）

張俊彪（1971－），男，高級工程師，主要從事河口近岸動力學(xué)研究。

通訊記者：許雪峰（1981-），男，博士研究生，主要從事河口近岸動力學(xué)研究。Email：xuxuefeng1981@163.com