典型海水溫度年計算方法及其應用（續）

端木琳 賈 欣 馬偉強 王 蕓

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

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

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

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

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

1 數據來源及處理

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

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

2 典型年計算

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

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

續表1

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

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

（1）

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

3 計算結果分析

3.1 數據代表性分析

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

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

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

3.2 海水溫度極端年

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

3.3 海水溫度逐年變化分析

為研究近17年來海水溫度的變化趨勢，本文將2005—2021年的海水溫度數據進行對比。圖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 ℃；經過同樣的方法分析得出，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年典型年進行對比發現，除了2月、5月和6月這三個月份的典型月沒有改變，其他月份的典型月均發生了變化，且溫度值均增加，增幅最大的為1月，由原來的4.07 ℃變為5.05 ℃增加了0.98 ℃，具體數據如表4所示。

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

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

圖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

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

4 結論

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

2）通過分析計算結果發現，近17年來，大連地區黃海近海海域的海水溫度均處于上升趨勢，黃海近海海域8月份近五年月平均溫度較2005—2009這五年上升約0.67 ℃，2月份上升約0.61 ℃，與2014年典型年相比，本文典型年月平均溫度除3個月份未發生改變之外，其他月份平均溫度均有不同程度的增加，其中1月增值最大，增加了0.98 ℃。大連地區渤海近海海域溫度變化趨勢與黃海近海海域不同，冬夏季海水溫差變大，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）根據大連市附近黃海、渤海海域海水溫度的特點，黃海海域夏季溫度更低，有利于熱泵夏季工況運行效率的提升，冬季溫度更高，有利于熱泵冬季工況運行效率的提升，因此黃海近海海域更適合作為海水源熱泵系統的低位冷熱源。但渤海近海海域冬季海水溫度較低，最低月平均溫度為-0.4 ℃，海水極端最低溫度為-1.8 ℃，海水不適合直接或單獨作為低位熱源，需要采用岸邊打井取水方式或具有其他輔助熱源的海水源熱泵形式。

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

致謝

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