端木琳 賈 欣 馬偉強 王 蕓

（1 大連理工大學(xué)土木工程學(xué)院 大連 116024；2 大連大學(xué)建筑工程學(xué)院 大連 116622）

目前，我國傳統(tǒng)的區(qū)域供熱和供冷方式主要為鍋爐和冷凍機，該方式通常需要消耗大量化石能源且污染環(huán)境。在“雙碳”目標之下，區(qū)域供熱和供冷方式的改革勢在必行，這使海水源熱泵的普及成為可能。在沿海地區(qū)，海水源熱泵系統(tǒng)雖然可以代替?zhèn)鹘y(tǒng)的鍋爐房和冷凍機為建筑物供熱、制冷，但并非任意從海洋中取出的海水都可作為熱泵系統(tǒng)的冷、熱源。

海水溫度對熱泵系統(tǒng)的設(shè)計和運行起著重要作用，是海水源熱泵技術(shù)應(yīng)用成敗的關(guān)鍵，主要考慮應(yīng)用海水源熱泵地區(qū)的海水的極端溫度[1-2]，極端溫度主要涉及海水最冷月和最熱月海水各層的溫度。我國大部分海域有很好的水溫條件，符合機組制冷和制熱對冷熱源溫度的要求[3]。海水源熱泵供熱設(shè)計時要求海水溫度不得低于2 ℃[4]，張莉等[5]針對青島市一海水源熱泵試點工程，通過實驗數(shù)據(jù)分析青島市附近黃海海域海水作為熱泵系統(tǒng)熱源時系統(tǒng)運行的可靠性，由實驗結(jié)果分析可知，1、2月份氣溫最低的冬季，系統(tǒng)運行可靠。陳高峰[6]對天津港附近渤海海域海水冬季水溫進行了測試，冬季海水作為低溫熱源，水溫較低，可能會出現(xiàn)結(jié)冰現(xiàn)象，設(shè)計中要考慮海水溫度過低系統(tǒng)無法運行時的措施。因此，渤海海域海水源熱泵實際工程需要考慮不同的設(shè)計方案，主要有三種：1）熱泵機組+冬季板式換熱器市政熱網(wǎng)供暖；2）熱泵機組+冬季備用熱源（鍋爐、電加熱器等）；3）采用海岸井滲濾取水的海水源熱泵系統(tǒng)[7-11]。

海水溫度全年變化，而海水溫度條件的改變將對熱泵循環(huán)特性乃至機組性能產(chǎn)生影響。季節(jié)的變化對陸地和海洋溫度變化起著決定性的作用，因此海水溫度的變化具有顯著的季節(jié)特征[12]。胡瑩英[13]分析了東海不同區(qū)域、不同斷面和不同層次的海水溫度年際變化規(guī)律。為滿足工程設(shè)計需要，得到具有代表性的海水溫度數(shù)據(jù)，端木琳等[14]基于2005—2014年海水溫度數(shù)據(jù)提出了典型海水溫度年（簡稱:2014年典型年），并給出了典型海水溫度年的計算方法，為海水源熱泵系統(tǒng)的設(shè)計工作奠定了基礎(chǔ)。但典型海水溫度年計算并非一勞永逸之事，全球氣候變暖，海洋溫度也隨之變化，根據(jù)其他國內(nèi)外多位學(xué)者的調(diào)查研究可知，近年來海水溫度一直處于上升趨勢[15-16]，我國的黃海、渤海海域也不例外，因此隨著全球氣候的變化，黃海、渤海海域典型海水溫度年需要進行更新。

綜上所述，本文在前期工作[14]的基礎(chǔ)上，統(tǒng)計了2005—2021年大連市黃海、渤海近海海域海水溫度數(shù)據(jù)，計算出大連市黃海、渤海海域（本文提到的“黃海、渤海海域”，無特殊說明均指大連地區(qū)近海海域）最新的典型海水溫度年并進行數(shù)據(jù)分析。

1 數(shù)據(jù)來源及處理

本文收集記錄了近17年（2005—2021年）大連市棒棰島海域、星海海域、付家莊海域和夏家河子海域近海岸表層海水溫度，海水溫度數(shù)據(jù)源于大連市海洋漁業(yè)局[17]。其中，以棒棰島、星海和付家莊海域近海岸表層海水溫度三者的平均值代表黃海海域海水溫度，以夏家河子海域近海岸表層海水溫度代表渤海海域海水溫度。

由于實測人員節(jié)假日休息或其他原因無法進行實測記錄，導(dǎo)致個別天數(shù)溫度數(shù)據(jù)缺失，則利用線性內(nèi)插法對缺失數(shù)據(jù)進行補充。夏家河子浴場個別年份1、2月的海水溫度數(shù)據(jù)缺失過多，采用差值法補充后導(dǎo)致計算誤差過大，計算結(jié)果不具有代表性，故未將這些年份的數(shù)據(jù)包含于累計分布函數(shù)。

2 典型年計算

本文沿用2014年典型年采用的FS（Finkelstein-Schafe）統(tǒng)計方法[18-20]進行計算。以黃海為例，通過FS統(tǒng)計方法計算每個月份的Fs值，計算結(jié)果如表1所示。

表1 海水溫度Fs計算值Tab.1 Fs statistics of seawater temperature

續(xù)表1

由計算過程可知，F(xiàn)s值的大小反映當月數(shù)據(jù)分布規(guī)律與多年累計分布規(guī)律的相似程度，F(xiàn)s值越小，其所在月份溫度分布規(guī)律與該月份過去多年長期分布越接近，二者分布規(guī)律越相似。因此，為保證典型年的溫度分布規(guī)律與過去多年海水溫度的長期分布規(guī)律相似，選取具有最小Fs值的月份作為典型月，由12個不同的典型月組成典型年[14]，如表2所示。同理計算出渤海海域典型海水溫度年，如表3所示。

基于DataFit軟件對得到的典型年黃海、渤海海水逐日溫度進行計算分析，得到黃海、渤海典型年海水溫度表達式：

（1）

式中：（1）對于黃海海域海水溫度，τ為天數(shù)（1月1日為起始）；a=12.52；b=-10.27；T=367.65；τ0=15.68；e=8.70，均為擬合系數(shù)。（2）對于渤海海域海水溫度，τ為天數(shù)（1月1日為起始）；a=12.04；b=-13.05；T=362.02；τ0=4.25；e=10.27，均為擬合系數(shù)。圖1所示為黃海、渤海海域典型年海水逐日溫度，其中散點表示典型年海水逐日溫度，實線表示擬合曲線。該公式能夠便于工程師根據(jù)實際工程設(shè)計需求計算海水溫度，作為工程設(shè)計基礎(chǔ)。

3 計算結(jié)果分析

3.1 數(shù)據(jù)代表性分析

表2 黃海典型海水溫度年組成Tab.2 Typical seawater temperature year in the Yellow Sea

表3 渤海典型海水溫度年組成Tab.3 Typical seawater temperature year in the Bohai Sea

對典型年月平均溫度和累年月平均溫度進行對比，可以觀察典型年溫度分布規(guī)律與溫度累年分布規(guī)律的相似性。圖2所示為黃海、渤海海域海水的月平均溫度對比，由圖2可知，典型年各月平均溫度與該月累年平均溫度很接近，兩條曲線幾乎重合。黃海、渤海海域典型年月平均溫度和累年月平均溫度最大差值分別出現(xiàn)在2月和12月，為0.9 ℃和0.63 ℃。因此，本文計算確定的典型海水溫度年溫度全年變化趨勢和分布規(guī)律與過去17年里的溫度變化趨勢和分布相似，能夠反映過去多年海水溫度的變化情況，具有較好的代表性。

3.2 海水溫度極端年

參考《中國建筑熱環(huán)境分析專用氣象數(shù)據(jù)集》[21]中獲得極端氣象數(shù)據(jù)的方法，本文采用歷年海水月平均溫度最大的月份組成海水溫度極高年；歷年海水月平均溫度最小的月份組成海水溫度極低年；采用累年海水日最高溫度為極端最高溫度；采用累年海水日最低溫度為極端最低溫度。黃海海水溫度極高年和溫度極低年月平均海水溫度如圖2（a）所示。渤海海水溫度極高年和溫度極低年月平均海水溫度如圖2（b）所示。黃海海水極端最高溫度為26.2 ℃；極端最低溫度為-0.4 ℃。渤海海水極端最高溫度為28.6 ℃；極端最低溫度為-1.8 ℃。

3.3 海水溫度逐年變化分析

為研究近17年來海水溫度的變化趨勢，本文將2005—2021年的海水溫度數(shù)據(jù)進行對比。圖3（a）所示為黃海海域2月海水平均溫度的逐年變化，2005—2009這五年的海水平均溫度為3.22 ℃，而2017—2021近五年的海水平均溫度為3.83 ℃，上升了0.61 ℃；圖3（b）所示為黃海海域8月平均溫度的逐年變化，2005—2009這五年的海水平均溫度為22.82 ℃，而2017—2021近五年的海水平均溫度為23.5 ℃，上升了0.67 ℃；經(jīng)過同樣的方法分析得出，2005—2009這五年的渤海2月海水平均溫度為1.14 ℃，而2017—2021近五年的渤海海水平均溫度為-0.72 ℃，2005—2009這五年的渤海8月海水平均溫度為25.16 ℃，而2017—2021近五年的渤海海水平均溫度為26.27 ℃。渤海海域2月份海水平均溫度降低了1.86 ℃，8月份的海水平均溫度上升了1.11 ℃。

圖1 典型年海水逐日溫度Fig.1 Seawater temperature of the typical year

圖2 黃海、渤海海域海水月平均溫度對比Fig.2 Comparison of mean monthly seawater temperature in the Yellow Sea and Bohai Sea

本文通過將計算確定的典型年與2014年典型年進行對比發(fā)現(xiàn)，除了2月、5月和6月這三個月份的典型月沒有改變，其他月份的典型月均發(fā)生了變化，且溫度值均增加，增幅最大的為1月，由原來的4.07 ℃變?yōu)?.05 ℃增加了0.98 ℃，具體數(shù)據(jù)如表4所示。

3.4 黃海、渤海海水溫度對比分析

受地理位置的影響，不同海域的海水溫度也存在差異。黃海海域海水全年月平均溫度約在2.8～22.9 ℃范圍內(nèi)變化，渤海海域海水全年月平均溫度約在-0.4～25.5 ℃范圍內(nèi)變化，渤海海域海水全年月平均溫度變化范圍大于黃海海域海水全年月平均溫變化范圍。黃海、渤海海域海水月平均溫度最大差值為5.8 ℃，出現(xiàn)在12月；最小差值為-0.7 ℃，出現(xiàn)在9月，具體數(shù)據(jù)如表5所示。圖4所示為黃海、渤海海域月平均溫度對比。由圖4可知，4—9月渤海海域海水月平均溫度高于黃海海域，10月至來年3月黃海海域海水月平均溫度高于渤海海域，即黃海海域海水夏季溫度更低，有利于熱泵夏季工況運行效率的提升（海水作為冷卻水，降低了機組冷凝溫度，提高了機組效率），冬季溫度更高，有利于熱泵冬季工況運行效率的提升（海水作為冷凍水，提高了機組蒸發(fā)溫度和效率），因此黃海近海海域更適合作為海水源熱泵系統(tǒng)的低位冷熱源。

圖3 黃海海域海水平均溫度逐年變化Fig.3 Changes in the mean seawater temperature in the Yellow Sea

表4 黃海海域典型海水溫度年月平均氣溫Tab.4 Mean monthly seawater temperature of the typical years in the Yellow Sea

表5 黃海、渤海海域典型海水溫度年月平均溫度Tab.5 Mean monthly seawater temperature of the 2021 typical year in the Yellow Sea and the Bohai Sea

圖4 黃海、渤海海域月平均溫度對比Fig.4 Comparison of mean monthly temperature of the Yellow Sea and the Bohai Sea

本文通過計算結(jié)果可知，黃海海域典型海水溫度年最低月平均溫度為2.8 ℃，海水極端最低溫度為-0.4 ℃。相比之下，渤海海域典型海水溫度年最低月平均溫度為-0.4 ℃，海水極端最低溫度為-1.8 ℃，這一點對于海水源熱泵系統(tǒng)是不利的[22]。

4 結(jié)論

1）本文以近17年（2005—2021年）大連市黃海、渤海近海逐日海水溫度數(shù)據(jù)為基礎(chǔ)，沿用2014年典型年的計算方法計算得到了大連市黃海、渤海近海海域最新的典型海水溫度年，經(jīng)過驗證，本文計算得到的典型年能夠反映過去多年海水溫度的變化情況，典型年溫度數(shù)據(jù)具有較好的代表性，并獲得了典型年黃海、渤海逐日海水溫度表達式。

2）通過分析計算結(jié)果發(fā)現(xiàn)，近17年來，大連地區(qū)黃海近海海域的海水溫度均處于上升趨勢，黃海近海海域8月份近五年月平均溫度較2005—2009這五年上升約0.67 ℃，2月份上升約0.61 ℃，與2014年典型年相比，本文典型年月平均溫度除3個月份未發(fā)生改變之外，其他月份平均溫度均有不同程度的增加，其中1月增值最大，增加了0.98 ℃。大連地區(qū)渤海近海海域溫度變化趨勢與黃海近海海域不同，冬夏季海水溫差變大，2005—2009這五年的渤海近海海域2月海水平均溫度為1.14 ℃，而2017—2021近五年的渤海近海海域海水平均溫度為-0.72 ℃，2005—2009這五年的渤海近海海域8月海水平均溫度為25.16 ℃，而2017—2021近五年的渤海近海海域海水平均溫度為26.27 ℃。渤海近海海域2月份海水平均溫度降低了1.86 ℃，8月份的海水平均溫度上升了1.11 ℃。

3）根據(jù)大連市附近黃海、渤海海域海水溫度的特點，黃海海域夏季溫度更低，有利于熱泵夏季工況運行效率的提升，冬季溫度更高，有利于熱泵冬季工況運行效率的提升，因此黃海近海海域更適合作為海水源熱泵系統(tǒng)的低位冷熱源。但渤海近海海域冬季海水溫度較低，最低月平均溫度為-0.4 ℃，海水極端最低溫度為-1.8 ℃，海水不適合直接或單獨作為低位熱源，需要采用岸邊打井取水方式或具有其他輔助熱源的海水源熱泵形式。

典型海水溫度年為海洋能新能源系統(tǒng)、海水源熱泵系統(tǒng)、以海水為冷卻水的發(fā)電系統(tǒng)等相關(guān)設(shè)計和方案論證提供了數(shù)據(jù)基礎(chǔ)。但典型海水溫度年的計算并不是一件簡單的工作，不同海域的典型年不同，同一海域不同區(qū)域的典型年也不同，而且隨著時間的推移，由于氣候變化，典型年還會發(fā)生改變，因此建立我國各海域不同區(qū)域的典型海水溫度年數(shù)據(jù)庫，并不斷更新，為海洋能新能源系統(tǒng)的設(shè)計提供基礎(chǔ)數(shù)據(jù)，助力“雙碳”目標的實現(xiàn)。

致謝

本文作者感謝大連市海洋漁業(yè)局大連市海洋預(yù)報臺提供的大連地區(qū)海水溫度數(shù)據(jù)，向其對海水源熱泵等海洋能利用技術(shù)的貢獻致敬。