陳 磊，朱曉曉，李國成，王 怡，葛佳艷

（漣水縣氣象局，江蘇漣水 223400）

漣水縣地處江蘇省北部、淮河下游、淮安境內(nèi)東北部，有“自古漣漪佳絕地”的美譽(yù)。境內(nèi)河網(wǎng)密布，水鄉(xiāng)特色濃郁，自然條件良好，適合水產(chǎn)品的種植、養(yǎng)殖。2020 年漣水縣全縣水產(chǎn)品種植養(yǎng)殖面積3 666.6 hm2，淡水產(chǎn)品產(chǎn)出1.8 萬t，實現(xiàn)產(chǎn)值4.02 億元。在水產(chǎn)種植養(yǎng)殖產(chǎn)業(yè)中，仍以露天作業(yè)環(huán)境為主，導(dǎo)致了該產(chǎn)業(yè)極易受到天氣因素的影響和制約，天氣條件對水產(chǎn)種植養(yǎng)殖業(yè)的影響巨大。

池塘水溫是水產(chǎn)種植、養(yǎng)殖最重要的水環(huán)境要素，不僅水體中的溶解氧、pH、氨氮含量等多項要素和水溫關(guān)系密切，而且水體中的各類魚、蝦、蟹和水生植物對水溫都非常敏感，影響水生生物的生長和繁殖。為了準(zhǔn)確預(yù)測水體溫度的變化，科技工作者進(jìn)行了大量的研究。倪玉紅等［1］應(yīng)用BP 神經(jīng)網(wǎng)絡(luò)對盱眙龍蝦池塘夏季水溫建立預(yù)報模型；張志勇等［2］利用數(shù)學(xué)統(tǒng)計方法，詳細(xì)探討了江蘇省南部沿海養(yǎng)殖池塘水溫變化規(guī)律；鄧愛娟等［3］利用逐步回歸法對洪湖地區(qū)水溫進(jìn)行了預(yù)報；楊文剛等［4］對湖北省池塘水溫建立了預(yù)報模型；還有許多研究成果也運(yùn)用在水溫預(yù)報中［5-28］。本研究擬從漣水觀測設(shè)備出發(fā)，監(jiān)測池塘中水體溫度的變化特征，挖掘其與氣溫的關(guān)系，為漣水縣池塘水溫預(yù)報提供可靠結(jié)論。

1 數(shù)據(jù)來源和處理方法

1.1 數(shù)據(jù)來源

水溫資料包括50 cm 水溫和80 cm 水溫，來源于漣水縣氣象局建設(shè)的水體監(jiān)測站（119°17＇54″E，33°46＇04″N），為了獲得真實試驗數(shù)據(jù)，該池塘不進(jìn)行抽水、灌水、增氧等人工調(diào)節(jié)活動，水體監(jiān)測站設(shè)置在面積為0.2 hm2的池塘內(nèi)，池塘水深1.8～2.0 m，傳感器距池塘岸3 m，水溫監(jiān)測設(shè)備為上海氣象儀器廠有限公司生產(chǎn)的DZZ3 型水產(chǎn)養(yǎng)殖自動氣象站。氣溫數(shù)據(jù)來源于漣水國家氣象觀測站，水體監(jiān)測站和漣水國家氣象觀測站相距150 m。水溫、氣溫資料的選用時間為2020 年3 月至2021 年2 月。

1.2 數(shù)據(jù)質(zhì)控

氣溫數(shù)據(jù)資料完整準(zhǔn)確，數(shù)據(jù)可用率100%。水溫資料因設(shè)備斷電維護(hù)，部分時段未采集數(shù)據(jù)，數(shù)據(jù)完整率為96.5%，缺測數(shù)據(jù)采用拉格朗日差值法補(bǔ)齊，處理后數(shù)據(jù)可用率為100%。

1.3 統(tǒng)計分析方法

水溫、氣溫日平均值為當(dāng)日20：00、08：00、14：00、02：00 的實測值的平均值。應(yīng)用SPSS 22 軟件進(jìn)行Pearson 相關(guān)性分析。Pearson 相關(guān)系數(shù)計算公式如下。

用逐步回歸法分別建立日平均水溫、日最高水溫、日最低水溫與氣溫要素的預(yù)報模型。逐步回歸法的具體步驟如下。

1）對p個回歸自變量X1，X2，…，Xp，分別同因變量Y建立一元回歸模型，

Y=β0+βiXi+?,i=1,…,p

計算變量Xi，相應(yīng)的回歸系數(shù)的F檢驗統(tǒng)計量的值，記為，取其中的最大值，即

依此方法重復(fù)進(jìn)行，每次從未引入回歸模型的自變量中選取1 個，直到經(jīng)檢驗沒有變量引入為止。

2 結(jié)果與分析

2.1 水溫與氣溫的變化規(guī)律

2.1.1 水溫與氣溫的月變化規(guī)律 50 cm 和80 cm的水溫的最大值均出現(xiàn)在8 月，分別為30.5 ℃和31.3 ℃；50 cm 和80 cm 水溫最小值均出現(xiàn)在1 月，分別為8.9 ℃和4.2 ℃。2—8 月水溫持續(xù)上升，9 月至次年1 月水溫持續(xù)下降，與氣溫變化趨勢一致。全年氣溫均低于水溫，除7—9 月外，50 cm 水溫均高于80 cm 水溫（圖1）。

圖1 漣水縣月平均水溫和氣溫變化

2.1.2 水溫與氣溫的日變化規(guī)律 白天空氣和水受熱輻射升溫，由于水的比熱容遠(yuǎn)大于空氣，水體的升溫速度小于空氣，水溫的變化幅度也小于空氣，夜間空氣輻射降溫的幅度也大于水溫變化的幅度，且日最高最低水溫出現(xiàn)時間均落后于氣溫。氣溫日最大值出現(xiàn)在14∶00，為19.8 ℃，最小值出現(xiàn)在06∶00，為12.4 ℃；50 cm 水溫日最大值出現(xiàn)在16∶00，為20.6 ℃，最小值出現(xiàn)在09∶00，為19.4 ℃；80 cm 水溫日最大值出現(xiàn)在17∶00，為18.9 ℃，最小值出現(xiàn)在07∶00，為17.2 ℃。50 cm 水溫平均值均高于80 cm水溫平均值和氣溫平均值，80 cm 水溫除11∶00 至16∶00 外，均高于氣溫（圖2）。

圖2 漣水縣小時平均水溫和氣溫變化

2.2 水溫預(yù)報模型的建立與檢驗

2.2.1 水溫預(yù)報模型 由于水產(chǎn)種植、養(yǎng)殖主要發(fā)生在春、夏、秋3 季，冬季水溫預(yù)報的意義不大，本研究僅探討研究了春、夏、秋3 季的預(yù)報模型。劃定3—5 月為春季，6—8 月為夏季，9—11 月為秋季，12月至次年2 月為冬季。

60～90 cm 是一般養(yǎng)殖池塘淡水魚生長發(fā)育、產(chǎn)卵繁殖的主要活動范圍，也可以代表池塘水溫的整體狀況，在2 層水溫中采用80 cm 水溫來建立預(yù)報模型。根據(jù)水溫與氣溫的變化規(guī)律發(fā)現(xiàn)，水溫與氣溫的變化趨勢相近，但水溫的變化相對滯后于氣溫，因此選取當(dāng)日和前1～5 d 的平均氣溫、最高氣溫、最低氣溫作為平均水溫、最高水溫、最低水溫進(jìn)行相關(guān)性分析。為方便表達(dá)，分別用T、TMAX、TMIN表示當(dāng)日平均水溫、最高水溫和最低水溫，T0、TMAX0、TMIN0分別表示當(dāng)日平均氣溫、最高氣溫和最低氣溫，T1、TMAX1、TMIN1分別表示前1 d 平均氣溫、最高氣溫和最低氣溫，T2、TMAX2、TMIN2分別表示前2 d 平均氣溫、最高氣溫和最低氣溫，T5、TMAX5、TMIN5分別表示前5 d 平均氣溫、最高氣溫和最低氣溫。

水溫與氣溫的相關(guān)系數(shù)見表1。由表1 可知，除夏季最低水溫外，T0與其他水溫的相關(guān)系數(shù)均大于0.831；T1與各季節(jié)水溫的相關(guān)系數(shù)均大于0.884；TMIN0與水溫相關(guān)系數(shù)均大于0.816；T、TMAX、TMIN與T0、TMIN0和T1的關(guān)系最為緊密。

表1 水溫與氣溫的相關(guān)系數(shù)

選取相關(guān)系數(shù)最大的前5 項氣溫要素作為平均水溫、最高水溫、最低水溫的預(yù)報因子，通過多元逐步回歸建立水溫預(yù)報模型如下。

春季平均水溫T=0.392×T0+0.231×T1+0.293×T2+4.268

春季最高水溫TMAX=0.595×T0+0.240×T1+0.282×T2-0.196×TMIN0+4.538

春季最低水溫TMIN=0.316×T0+0.277×T1+0.311×T2+3.368

夏季平均水溫T=0.325×T0+0.355×T1+0.325×TMIN2+3.215

夏季最高水溫TMAX=0.363×T0+0.345×T1+0.233×TMIN1+5.588

夏季最低水溫TMIN=0.384×T1+0.354×TMIN0+0.264×TMIN2+3.370

秋季平均水溫T=0.372×T0+0.287×T1+0.400×T3+1.699

秋季最高水溫TMAX=0.409×T0+0.263×T1+0.408×T3+2.126

秋季最低水溫TMIN=0.368×T0+0.281×T1+0.391×T3+1.383

2.2.2 預(yù)報模型檢驗 根據(jù)得出的預(yù)報模型，計算出春、夏、秋3 個季節(jié)平均水溫、最高水溫和最低水溫的預(yù)測值，實測水溫值與模擬水溫值能較好地線性擬合，可直觀地看出散點分布在直線兩側(cè)，春季水溫的平均值、最大值、最低值與實測值的相關(guān)系數(shù)均大于0.94，夏季水溫的平均值、最大值、最低值與實測值的相關(guān)系數(shù)均大于0.87，秋季水溫的平均值、最大值、最低值與實測值的相關(guān)系數(shù)均大于0.97，表明各預(yù)報模型的水溫預(yù)測值均與實測值的相關(guān)性較好（圖3）。

圖3 預(yù)報水溫和實測水溫對比

對實測值和預(yù)報值進(jìn)行誤差分析，結(jié)果見表2。由表2 可知，春、夏、秋3 季的平均水溫、最高水溫、最低水溫的預(yù)報模擬值與實測值的平均絕對誤差在0.62～0.89 ℃，平均相對誤差在2.16%～4.92%。建立的水溫預(yù)報模型對平均水溫、最高水溫、最低水溫的預(yù)測值也可達(dá)到較高的精度，可用于短期和中長期水溫預(yù)報。

表2 預(yù)報水溫和實測水溫誤差對比

3 小結(jié)

1）1 年中池塘水溫最高值出現(xiàn)在8 月，最低值出現(xiàn)在1 月，月變化趨勢穩(wěn)定。池塘水溫與氣溫變化相關(guān)性較高，呈正相關(guān)關(guān)系，但變化幅度小于氣溫，變化時間滯后于氣溫。

2）當(dāng)日平均水溫、最高水溫、最低水溫與當(dāng)日平均氣溫、當(dāng)日最低氣溫和前1 d 平均氣溫的關(guān)系最為緊密。

3）以氣溫為預(yù)報變量，建立春、夏、秋3 季平均水溫、最高水溫、最低水溫預(yù)報模型，得到的模擬值與實測值相關(guān)系數(shù)均大于0.87，具有較好的預(yù)報精度，各季節(jié)預(yù)報模型平均水溫、最高水溫、最低水溫模擬結(jié)果與實測值均低于1 ℃，相對誤差小于5%，模擬效果較好，可以為當(dāng)?shù)厮a(chǎn)種植養(yǎng)殖提供科學(xué)可靠的水溫氣象預(yù)報服務(wù)。

池塘水溫的預(yù)測模擬受池塘大小、深淺和周邊環(huán)境等因素影響較大。本研究結(jié)論僅可作為當(dāng)?shù)丨h(huán)境下的池塘水溫預(yù)測，除相關(guān)定性結(jié)論外，不宜外推直接應(yīng)用到其他地區(qū)環(huán)境下的池塘水溫預(yù)報。