王鑫,武賀,王靜,張巍
(國(guó)家海洋技術(shù)中心,天津 300112)
在目前全球能源危機(jī)與環(huán)境問(wèn)題日益突出的背景下,潮流能這一可再生能源日益顯現(xiàn)出其廣闊的開(kāi)發(fā)利用前景。在各國(guó)清潔能源激勵(lì)政策與新興能源產(chǎn)業(yè)推動(dòng)政策的共同作用下,潮流能開(kāi)發(fā)利用技術(shù)進(jìn)入一個(gè)前所未有的快速發(fā)展階段。潮流能裝置的年發(fā)電量是其運(yùn)行效果的重要評(píng)判依據(jù)之一,但通過(guò)統(tǒng)計(jì)裝置長(zhǎng)時(shí)間運(yùn)行的累計(jì)電量來(lái)確定其年發(fā)電量,需要耗費(fèi)大量的人力物力,且在現(xiàn)場(chǎng)惡劣環(huán)境下監(jiān)測(cè)裝置本身的損壞與丟失也時(shí)有發(fā)生,顯然這種方法的可操作性較差。鑒于此,本文的研究基于統(tǒng)計(jì)分析相關(guān)原理,對(duì)短期現(xiàn)場(chǎng)監(jiān)測(cè)獲得的潮流數(shù)據(jù)和電功率輸出數(shù)據(jù)進(jìn)行分析,進(jìn)而實(shí)現(xiàn)了對(duì)裝置年發(fā)電量的估算,從而為評(píng)判潮流能裝置的運(yùn)行效果提供了更具可操作性的方法。
為獲得用于估算的潮流數(shù)據(jù)和電功率輸出數(shù)據(jù),在潮流能裝置正前方來(lái)流方向設(shè)置潮流觀測(cè)儀,布放的儀器距潮流能裝置的距離既夠保證測(cè)得流速代表裝置布放位置流速,又可避免潮流能裝置受到流速儀器引起的渦旋與渦流影響。根據(jù)裝置運(yùn)行的具體環(huán)境,一般取大于10 倍裝置葉輪直徑的距離布放測(cè)流裝置。對(duì)于雙向流發(fā)電的潮流能裝置,則需要在裝置前后各設(shè)置一臺(tái)潮流觀測(cè)儀,以保證在潮流能裝置正反向發(fā)電時(shí),均能夠獲得來(lái)流方向的流速。在大潮潮差日至小潮潮差日之間連續(xù)觀測(cè)15 d,對(duì)于在正規(guī)半日潮或正規(guī)日潮海域連續(xù)觀測(cè)7~8 d(The European Marine Energy Centre Ltd,2009)。
由潮流觀測(cè)儀測(cè)試潮流能裝置設(shè)置處潮流流速,由潮流能裝置電力輸出端設(shè)置的電功率檢測(cè)儀測(cè)試電功率數(shù)據(jù),潮流測(cè)試與電功率測(cè)試需同步進(jìn)行。每5 min 獲得一次平均潮流流速數(shù)據(jù),同時(shí)每5 min 獲取一次平均電功率數(shù)據(jù),每個(gè)潮流流速數(shù)據(jù)需與同時(shí)間段采集的電功率數(shù)據(jù)一一對(duì)應(yīng)。相對(duì)應(yīng)的每一個(gè)潮流流速數(shù)據(jù)與電功率數(shù)據(jù)組成一個(gè)數(shù)據(jù)點(diǎn)。
潮流流速低于潮流能裝置的最低設(shè)計(jì)運(yùn)行流速或高于最高設(shè)計(jì)運(yùn)行流速時(shí),潮流能裝置會(huì)停止工作,對(duì)于單向發(fā)電的潮流能裝置而言,反向來(lái)流時(shí)裝置也會(huì)停止工作,裝置停止工作期間沒(méi)有電力輸出,在進(jìn)行數(shù)據(jù)分析處理前需要將上述潮流數(shù)據(jù)篩選出來(lái)并剔除掉。
分析裝置的功率特性之前,需要首先將獲得的數(shù)據(jù)點(diǎn)進(jìn)行分組,分組依據(jù)每個(gè)數(shù)據(jù)點(diǎn)所包含流速值的大小來(lái)進(jìn)行,劃分到每組的數(shù)據(jù)點(diǎn),均是滿足其包含的流速值大于等于某一流速值,且小于另一流速值的數(shù)據(jù)點(diǎn),作為下限的流速值定義為組的起點(diǎn),作為上限的流速值定義為組的終點(diǎn)。因?yàn)榉纸M的目的是對(duì)被測(cè)試的潮流能裝置進(jìn)行功率特性分析,所以所分的組就命名為分析區(qū)間,則組的起點(diǎn)和終點(diǎn),即為分析區(qū)間的起點(diǎn)和終點(diǎn),每個(gè)分析區(qū)間終點(diǎn)值和起點(diǎn)值的差值定義為該分析區(qū)間的區(qū)間寬度。劃分的所有分析區(qū)間需首尾相連,區(qū)間與區(qū)間之間無(wú)重合,且保證所有的數(shù)據(jù)點(diǎn)均能夠且僅能夠分入一個(gè)所劃分的分析區(qū)間中。
進(jìn)行潮流能功率特性分析時(shí),分析區(qū)間的區(qū)間寬度取0.2 m,保證有不少于10 個(gè)數(shù)據(jù)點(diǎn)落入每個(gè)分析區(qū)間中,如落入某個(gè)分析區(qū)間的數(shù)據(jù)點(diǎn)少于10 個(gè),則可擴(kuò)大分析區(qū)間寬度,以增加落入?yún)^(qū)間內(nèi)的數(shù)據(jù)點(diǎn),但分析區(qū)間寬度不大于0.3 m。分析區(qū)間確定后,考查各數(shù)據(jù)點(diǎn)所包含的流速值,則根據(jù)流速值的大小,所得到的數(shù)據(jù)點(diǎn)會(huì)落入不同的分析區(qū)間,將每個(gè)區(qū)間中數(shù)據(jù)點(diǎn)的流速值與電功率值分別作算術(shù)平均,每個(gè)分析區(qū)間的流速與電功率算術(shù)平均值構(gòu)成該分析區(qū)間的代表數(shù)據(jù)點(diǎn)。利用計(jì)算得到的各分析區(qū)間代表數(shù)據(jù)點(diǎn),利用三次最小二乘法對(duì)代表數(shù)據(jù)點(diǎn)進(jìn)行數(shù)值擬合,獲得裝置的功率輸出功率特性曲線。最后,統(tǒng)計(jì)落入每個(gè)分析區(qū)間的數(shù)據(jù)點(diǎn)個(gè)數(shù)在數(shù)據(jù)點(diǎn)總數(shù)中所占的比例,全年按8 760 小時(shí)計(jì)算,利用此前計(jì)算得到的各分析區(qū)間的代表數(shù)據(jù)點(diǎn)中包含的電功率算數(shù)平均值,利用如下公式估算出裝置的年發(fā)電量(盛驟等,1999)。
式中:AEP 為年發(fā)電量;hAE為全年小時(shí)數(shù);n 為參與計(jì)算的分析區(qū)間個(gè)數(shù);ai為落入第i 個(gè)分析區(qū)間的數(shù)據(jù)點(diǎn)個(gè)數(shù);pi為第i 個(gè)分析區(qū)間的代表輸出功率;N 為數(shù)據(jù)點(diǎn)總數(shù);
本文以2011年12月24日-12月29日(農(nóng)歷初一朔望大潮期),依照前文所述數(shù)據(jù)獲取方法,對(duì)岱山10 kW 潮流能發(fā)電裝置開(kāi)展測(cè)試獲得數(shù)據(jù)為例說(shuō)明方法的具體分析步驟,同時(shí)利用該潮流能裝置長(zhǎng)期運(yùn)行獲得的歷史數(shù)據(jù),對(duì)本方法科學(xué)性進(jìn)行比對(duì)說(shuō)明。本次測(cè)試選用浪龍作為測(cè)流設(shè)備,得到的垂向剖面上的流速,分為4 層,1-4 層距底的距離分別為1.5 m、2.5 m、3.5 m 和4.5 m,對(duì)剖面流速進(jìn)行垂向平均,作為該時(shí)刻的平均流速進(jìn)行分析。
電能輸出功率計(jì)算公式:
式中:P 為輸出功率,單位瓦(W);UDC為蓄電池電壓,單位伏特(V);ID為充電電流,單位安培(A)
此次觀測(cè)每5 min 獲得1 組潮流數(shù)據(jù),共獲得1 331 組,獲得數(shù)據(jù)范圍為0~2.48 m/s,分成若干區(qū)間:小于0.6 m/s 的流速為范圍Ⅰ;0.6~2.6 m/s流速為范圍Ⅱ,該范圍內(nèi)區(qū)間寬度取0.2 m/s。檢測(cè)試驗(yàn)獲得的全部1 331 個(gè)數(shù)據(jù),落入范圍Ⅰ的基本數(shù)據(jù)記錄點(diǎn)數(shù)共為259 個(gè),落入范圍Ⅱ的基本數(shù)據(jù)記錄點(diǎn)共為1 072 個(gè)(圖1 和表1)。統(tǒng)計(jì)落入各區(qū)間的流速值,取區(qū)間內(nèi)的所有流速值的算術(shù)平均值,將獲得的流速值,作為該區(qū)間的代表流速值。
圖1 各區(qū)間流速數(shù)據(jù)個(gè)數(shù)統(tǒng)計(jì)圖
表1 各區(qū)間代表流速和代表輸出功率數(shù)據(jù)表
每5 min 取一次電平均輸出功率,此次檢測(cè)共獲得506 個(gè)數(shù)據(jù)。獲得電平均輸出功率數(shù)與同時(shí)段的各流速值對(duì)應(yīng),將各流速區(qū)間內(nèi)與流速對(duì)應(yīng)的電輸出功率值取算術(shù)平均,得到值作為該區(qū)間代表電輸出功率值。
代表流速與代表電輸出功率作為該區(qū)間的代表點(diǎn),該點(diǎn)可以較準(zhǔn)確反映對(duì)應(yīng)流速下的電功率輸出特性,將各代表點(diǎn)做擬合處理,得出潮流能系統(tǒng)功率輸出特性曲線(圖2)。
可依式(1)、(2),全年小時(shí)數(shù)hAE取8 760,計(jì)算得到數(shù)據(jù)點(diǎn)總數(shù)N 取1 688,參與計(jì)算的分析區(qū)間個(gè)數(shù)n 取11,計(jì)算得到裝置全年連續(xù)正常運(yùn)行的情況下,其年發(fā)電量約為19 896 千瓦時(shí)。
圖2 潮流能系統(tǒng)輸出功率特性曲線
本文的研究成果基于短期的實(shí)海況觀測(cè),對(duì)潮流能裝置的年發(fā)電量進(jìn)行評(píng)價(jià),評(píng)價(jià)周期短,評(píng)價(jià)成本低,且評(píng)價(jià)方法具有可靠的理論依據(jù),可為政府科技主管部門(mén)的項(xiàng)目驗(yàn)收,科研機(jī)構(gòu)的技術(shù)研發(fā)與改進(jìn),以及企業(yè)的投資技術(shù)選型,提供科學(xué)參考依據(jù),該方法已在相關(guān)項(xiàng)目成果驗(yàn)收的過(guò)程中得到了很好的應(yīng)用。當(dāng)然,應(yīng)該看到,不同月份大潮的潮位流速變化通常是不完全相同的,利用短期的觀測(cè)結(jié)果估算年發(fā)電量會(huì)產(chǎn)生一定的誤差,如可獲得當(dāng)?shù)馗L(zhǎng)時(shí)間的潮流歷史資料,結(jié)合該方法,則可明顯提高估算結(jié)果的質(zhì)量。此外,本文所述方法對(duì)于波浪能裝置年發(fā)電量的估算同樣適用。
The European Marine Energy Centre Ltd, 2009. Assessment of Performance of Wave Energy Conversion Systems,Marine Renewable Energy Guides.
The European Marine Energy Centre Ltd, 2009. Assessment of Performance of Tidal Energy Conversion Systems,Marine Renewable Energy Guides.
The European Marine Energy Centre Ltd,2009.Assessment of Tidal Energy Resource,Marine Renewable Energy Guides.
The European Marine Energy Centre Ltd,2009.Guidelines for Manufacturing,Assembly and Testing of Marine Energy Conversion Systems,Marine Renewable Energy Guides.
盛驟,謝世千,潘承毅,1999年3月,《概率論與數(shù)理統(tǒng)計(jì)》,高等教育出版社.