董雙發(fā), 范 曉, 石海崗, 許莉萍, 章新益
(1.核工業(yè)航測(cè)遙感中心,石家莊 050002; 2.河北省航空探測(cè)與遙感技術(shù)重點(diǎn)實(shí)驗(yàn)室,石家莊 050002; 3.福建福清核電有限公司,福州 350318)
核電運(yùn)行過程中,會(huì)產(chǎn)生大量的余熱,為保證設(shè)備安全,需利用循環(huán)冷卻水將多余熱量排放到周邊環(huán)境中[1]。周邊水體溫度升高到一定程度時(shí)將改變水體質(zhì)量,影響各類水生物的生長(zhǎng)及繁殖活動(dòng)[2]。當(dāng)核電取水區(qū)域受其影響時(shí),可能會(huì)威脅到核電冷卻水取水安全。因此,研究核電附近水域溫排水的分布具有重要意義[3]。
溫排水監(jiān)測(cè)的方法主要有數(shù)值模擬、大面積實(shí)地測(cè)量及遙感監(jiān)測(cè)等方法[4-5]。熱紅外遙感技術(shù)因同步性好、可重復(fù)觀測(cè)和成本相對(duì)較低等優(yōu)勢(shì)成為溫排水監(jiān)測(cè)的重要方法[6]。國內(nèi)外學(xué)者針對(duì)傳感器熱紅外通道設(shè)置情況,提出了諸如單通道法[7]、普適性單窗算法[8]、劈窗算法和多通道法[9]等溫度反演算法。研究者利用這些算法,對(duì)水體進(jìn)行溫度反演取得了較高的反演精度,并將其應(yīng)用于核電附近海域溫排水分布調(diào)查研究[10-14]。無人機(jī)航測(cè)具備機(jī)動(dòng)靈活、分辨率高的優(yōu)勢(shì),隨著無人機(jī)技術(shù)的逐步成熟及熱紅外儀器的小型化,該技術(shù)正在成為核電廠溫排水監(jiān)測(cè)的新趨勢(shì)[15-16]。
目前,基于遙感技術(shù)對(duì)溫排水監(jiān)測(cè)的研究多為衛(wèi)星過境時(shí)刻的結(jié)果,對(duì)不同季節(jié)、不同潮汐狀態(tài)下的研究結(jié)果罕有報(bào)道。本次研究選用Landsat8和無人機(jī)熱紅外數(shù)據(jù),對(duì)福清核電附近海域的海表溫度開展溫度反演,分析不同潮汐、不同季節(jié)福清核電附近海域的溫度場(chǎng)分布特征,利用無人機(jī)監(jiān)測(cè)數(shù)據(jù)對(duì)溫排水影響區(qū)域的水溫分布細(xì)節(jié)進(jìn)行遙感監(jiān)測(cè)與分析,旨在促進(jìn)無人機(jī)遙感技術(shù)在溫排水遙感監(jiān)測(cè)領(lǐng)域的應(yīng)用。
基于2019年12月11日Landsat8數(shù)據(jù)的B6,B5,B2波段和全色波段,將影像融合為15 m空間分辨率,疊加地理要素后成圖,可見福清核電廠位于福建省福清市三山鎮(zhèn)南端,東、西、南三面環(huán)海,東北側(cè)與陸地連接,周邊道路發(fā)達(dá),交通較為便利(圖1)。廠址共規(guī)劃6臺(tái)百萬千瓦級(jí)核電機(jī)組,1—4號(hào)機(jī)組為二代改進(jìn)型壓水堆機(jī)組,5和6號(hào)機(jī)組為我國自主研發(fā)的“華龍一號(hào)”核電機(jī)組。2014年11月以來,福清核電站1—4號(hào)機(jī)組相繼投產(chǎn),2021年1月30日全球首臺(tái)“華龍一號(hào)”機(jī)組投入商業(yè)運(yùn)行。

圖1 福清核電地理位置示意圖(2019年12月11日)Fig.1 Location of the Fuqing nuclear power plant
福清核電廠址區(qū)域?qū)俚湫蛠啛釒ШQ笮约撅L(fēng)氣候,海域多年平均氣溫為22.3 ℃; 廠址溫排水海域潮流為規(guī)則半日潮,每個(gè)潮汐日有2次高潮和2次低潮,運(yùn)動(dòng)形式為往復(fù)流,漲潮流向偏西北,落潮流向偏東南,落潮流歷時(shí)稍長(zhǎng)于漲潮流歷時(shí)。
本文數(shù)據(jù)源為L(zhǎng)andsat8和無人機(jī)熱紅外數(shù)據(jù)。Landsat8數(shù)據(jù)具體參數(shù)和技術(shù)指標(biāo)參見文獻(xiàn)[17],數(shù)據(jù)時(shí)相分別為2013年8月4日、2014年1月27日、2018年12月1日、2019年7月29日、2019年12月11日、2020年8月7日,無人機(jī)數(shù)據(jù)時(shí)相為2019年7月15日。為驗(yàn)證反演結(jié)果,獲取了2019年7月15日、2019年7月29日近同步海面測(cè)量數(shù)據(jù)。熱紅外波段參數(shù)、過境時(shí)間、機(jī)組運(yùn)行工況、所處季節(jié)和潮態(tài)見表1。根據(jù)Landsat8波段設(shè)置情況,對(duì)數(shù)據(jù)進(jìn)行處理,包括幾何精糾正、工作區(qū)裁剪、水陸分離等。幾何精糾正基于地形數(shù)據(jù),誤差控制在1個(gè)像元(15 m)之內(nèi); 海陸分離利用Landsat8 OLI綠光波段和短波紅外波段對(duì)陸地和水體反射率差異開展波段比值獲取; 各影像按照研究區(qū)矢量進(jìn)行裁剪。

表1 不同傳感器熱紅外數(shù)據(jù)參數(shù)信息對(duì)比Tab.1 Parameter comparison for thermal band of different sensors
學(xué)者針對(duì)Landsat8的2個(gè)熱紅外波段開展了劈窗算法研究[18-19],但因第11波段存在嚴(yán)重的條帶噪聲,使得該方法存在較大的誤差。本文采用輻射傳輸方程算法,對(duì)其第10波段開展溫度反演。
在無云情況下,不考慮大氣對(duì)電磁波的散射,水平大氣各種組分混合均勻,設(shè)溫度為TS,則衛(wèi)星傳感器獲取的大氣頂層輻射可表示為:
Lλ=ελτλLλ(TS)+(1-ελ)τλLλatm↓+Lλatm↑,
(1)
式中:Lλ為傳感器接收到的大氣頂層輻射,可由傳感器輻射定標(biāo)獲取;TS為地物實(shí)際溫度;Lλ(TS)為地物溫度為TS時(shí)的黑體輻射;ελ為地物的比輻射率;Lλatm↓和Lλatm↑分別是大氣下行輻射和大氣上行輻射;τλ為波長(zhǎng)為λ時(shí)地表和傳感器之間的大氣透射率。
將式(1)變形,可獲取地物黑體輻射Lλ(TS)時(shí)與大氣頂層輻射Lλ及各項(xiàng)大氣作用參數(shù)的關(guān)系,將各參數(shù)逐一解算,即可求得Lλ(TS),公式為:
(2)
1)大氣頂層輻射Lλ的計(jì)算。對(duì)Landsat8第10波段灰度值輻射定標(biāo)獲取。定標(biāo)公式為:
Lλ=a·DN+b,
(3)
式中:DN為像元灰度值;a和b為L(zhǎng)andsat8第10波段的定標(biāo)系數(shù),可以直接從元數(shù)據(jù)中獲取。
2)海表比輻射率ελ。與其他地物不同,本次溫度反演主要針對(duì)海水表層進(jìn)行,因其比輻射率變化不大,接近黑體[20],本次研究取定值0.995。
3)大氣參數(shù)的獲取。Lλatm↓,Lλatm↑和τλ均與大氣成分含量、溫度、濕度、壓力密切有關(guān)。根據(jù)MODTRAN輻射傳輸模型,輸入衛(wèi)星過境時(shí)刻核電氣象站獲取的溫度、壓力、相對(duì)濕度等氣象數(shù)據(jù)、平均高程數(shù)據(jù)和NCEP大氣廓線數(shù)據(jù),開展大氣校正[21-22],即可獲取上述參數(shù)。
4)海表溫度TS的計(jì)算。在以上計(jì)算完成后,獲取溫度TS時(shí)的黑體輻射Lλ(TS),根據(jù)普朗克法則,求得TS,公式為:
(4)
式中K1和K2分別為校正系數(shù),均可從Landsat8數(shù)據(jù)元數(shù)據(jù)獲取,通過波段運(yùn)算,獲得福清核電附近海域海面溫度場(chǎng)分布數(shù)據(jù)(圖2、表2)。



圖2 熱紅外溫度場(chǎng)圖

表2 溫度場(chǎng)監(jiān)測(cè)結(jié)果概況Tab.2 Monitoring results of the distribution of thermal infrared temperature (℃)
無人機(jī)平臺(tái)為國產(chǎn)F300固定翼無人機(jī),熱紅外成像儀為美國產(chǎn)FLIR TAU2,無人機(jī)數(shù)據(jù)單張圖像大小為640像素×512像素,溫度分辨率可達(dá)0.05 ℃。無人機(jī)平臺(tái)系統(tǒng)配備50 Hz高精度GNSS板卡,支持PPK,RTK以及PPK和RTK融合作業(yè)模式。
根據(jù)氣象情況,2019年7月15日開展了無人機(jī)航飛,航高為380 m,航速為60 km/h,地面分辨率為0.5 m,航向重疊80%,旁向重疊55%,共飛行12架次,獲取有效影像數(shù)量8 252張,航程約350 km。
無人機(jī)數(shù)據(jù)包括陸地和海域2個(gè)部分,陸地紋理信息豐富,通過基于特征匹配的方法進(jìn)行拼接[23]。海域紋理較弱,首先利用無人機(jī)照片POS數(shù)據(jù)點(diǎn)位信息、確定單個(gè)無人機(jī)影像的4個(gè)角點(diǎn)的地面坐標(biāo),利用坐標(biāo)拼接圖像[24],其次匹配架次之間、影像之間存在的色彩差異,采用勻色算法進(jìn)行影像色彩調(diào)整,實(shí)現(xiàn)色彩一致。對(duì)于圖像中錯(cuò)誤的區(qū)域,進(jìn)行影像的位置調(diào)整。經(jīng)處理后,可以得到色彩相對(duì)一致,位置準(zhǔn)確的圖像(圖3)。圖4為無人機(jī)航測(cè)數(shù)據(jù)與實(shí)際數(shù)據(jù)的擬合,圖5為獲得的水體表面溫度。由于FLIR TAU2已經(jīng)進(jìn)行了輻射校正,可獲取熱紅外波段輻射亮度值和溫度值。因海域輻射受多方面影響,測(cè)溫型機(jī)芯獲取的溫度結(jié)果不可避免會(huì)產(chǎn)生偏差。因此,在無人機(jī)監(jiān)測(cè)的同時(shí),開展了近

圖3 無人機(jī)熱紅外影像拼接圖Fig.3 Mosaic of the UAV thermal infrared image

圖4 無人機(jī)航拍影像的溫度標(biāo)定Fig.4 Temperature calibration for UAV image

圖5 2019年7月15日海面溫度場(chǎng)分布Fig.5 Distribution of thermal infraredtemperature in Jul. 15, 2019
同步的海溫測(cè)量,用于糾正偏差,獲取真實(shí)海面溫度。對(duì)同步獲取的測(cè)量數(shù)據(jù)及航拍數(shù)據(jù)進(jìn)行點(diǎn)對(duì)點(diǎn)的時(shí)空匹配,時(shí)間尺度設(shè)定為±30 s,空間匹配尺度設(shè)定為1像元,共獲取數(shù)據(jù)134組,將實(shí)測(cè)數(shù)據(jù)和無人機(jī)航測(cè)溫度數(shù)據(jù)進(jìn)行擬合(圖4),獲取到水體真實(shí)的表面溫度(圖5)。結(jié)果表明2019年7月15日福清核電排次口周邊溫度為30.1~34.1 ℃,取水口周邊溫度為29.3~30.1 ℃。
為了保證反演結(jié)果的準(zhǔn)確性,在2019年7月15日無人機(jī)航測(cè)和2019年7月29日Landsat8衛(wèi)星過境時(shí)刻,開展了近同步海面溫度測(cè)量。測(cè)量設(shè)備為JENCO 6010M水質(zhì)測(cè)量?jī)x,標(biāo)定后儀器溫度測(cè)量精度為0.1 ℃。測(cè)量時(shí)每50~100 m間距進(jìn)行一次測(cè)溫,以排水口為中心,測(cè)點(diǎn)布置近密遠(yuǎn)疏,控制溫度變化。
對(duì)熱紅外遙感數(shù)據(jù)反演值和實(shí)測(cè)值使用最小二乘法進(jìn)行線性回歸擬合,以檢查其相關(guān)程度。為衡量實(shí)際值與理論預(yù)測(cè)值的偏離程度引進(jìn)了相關(guān)系數(shù)r、標(biāo)準(zhǔn)誤差Ω和殘差R,并對(duì)其進(jìn)行標(biāo)準(zhǔn)化處理。數(shù)據(jù)處理時(shí)定義實(shí)測(cè)值為x,溫度反演值為y,則回歸直線方程為y=a+bx,主要參數(shù)的計(jì)算公式為:
(5)
(6)
Ry=yi-y=yi-a-bxi。
(7)

無人機(jī)數(shù)據(jù)溫度反演值與實(shí)測(cè)值的擬合后結(jié)果為y= 0.891x+ 3.353 8,擬合后回歸系數(shù)的平方值R2為 0.914,標(biāo)準(zhǔn)誤差為0.540 7(圖6),殘差最大為4.05 ℃,位于排水口附近,80%的數(shù)據(jù)集中在(-0.5,0.5) ℃,殘差絕對(duì)值的平均值為0.28 ℃。Landsat8溫度反演結(jié)果與實(shí)測(cè)值擬合后結(jié)果為y= 0.814 1x+ 6.125 7,擬合后回歸系數(shù)的平方值R2為0.897 1,標(biāo)準(zhǔn)誤差為0.679 2(圖7),殘差最大為1.37 ℃,同樣位于排水口附近,70%的數(shù)據(jù)集中在(-0.5,0.5) ℃,殘差絕對(duì)值的平均值為0.44 ℃。綜上,除排水口附近因水動(dòng)力條件較強(qiáng),溫度變化略大外,其余處測(cè)點(diǎn)和實(shí)測(cè)關(guān)系具有較強(qiáng)的線性回歸關(guān)系,從線性關(guān)系程度、誤差、殘差上可以反映出熱紅外遙感反演結(jié)果準(zhǔn)確可靠。


2013年8月4日和2014年1月6日福清核電廠無機(jī)組運(yùn)行,排水口、取水口附近海水溫度與鄰近區(qū)域海水差異較小,不同季節(jié)溫度差別較大(圖2(a)—(b)),夏季和冬季溫度分別在23.0~24.0 ℃和9.0~10.0 ℃之間,從排水口到鄰近海域溫度緩慢過渡。
2018年后4景Landsat8數(shù)據(jù)獲取時(shí)4臺(tái)機(jī)組均為滿功率運(yùn)行,無人機(jī)航飛時(shí)2臺(tái)機(jī)組滿功率,2臺(tái)機(jī)組降功率運(yùn)行。溫度反演結(jié)果顯示(圖2(c)—(f)),福清核電運(yùn)行后,排水口周邊海域溫度場(chǎng)具有明顯的水溫分異與溫升現(xiàn)象,相似潮態(tài)下溫度梯度及空間分布特征相似: 核電排水口附近海域受溫排水影響,距離排水口越近溫度越高,遠(yuǎn)離排水口后,溫度逐漸降低,到達(dá)一定距離后,變化趨緩并逐漸穩(wěn)定。
與運(yùn)行前數(shù)據(jù)相比,夏季(圖2(c)和(d))核電排水口附近海域溫度范圍為30.0~35.8 ℃,取水口附近海域溫度范圍為29.3~31.8 ℃; 冬季(圖2(e)和(f))核電排水口附近海域溫度范圍為16.5~23.8 ℃,取水口附近海域溫度范圍為16.5~19.7 ℃。
核電周邊海域整體溫度受季節(jié)影響,夏季比冬季溫度場(chǎng)高10~15 ℃。溫度場(chǎng)分布特征則明顯受到了海域漲、落潮的影響,落潮時(shí)(圖2(c),(e),圖5),高溫冷卻水出排水口后,主要沿潮水向南、東南方向展布,漲潮則相反(圖2(f)),主要向西北方向和東北方向展布(圖2(d)和(f)),隨著遠(yuǎn)離排水口逐漸降溫,最終與環(huán)境溫度變化一致。
為獲取溫排水對(duì)周邊水域造成的熱影響,首先要剔除海域環(huán)境本底溫度。對(duì)于福清核電附近海域,由于遠(yuǎn)近海區(qū)域存在一定的溫度梯度,無法選取遠(yuǎn)海的溫度作為背景溫度。經(jīng)對(duì)比研究核電運(yùn)行前溫度場(chǎng)數(shù)據(jù)(圖2(a)和(b)),核電排水口和取水口附近溫度基本一致,且不同季節(jié)同步變化。因此,本次本底溫度取值選取取水口北側(cè)不受溫排水影響海域平均溫度作為本底溫度,取值海域范圍如圖8(a)中紅框所示。將海域溫度場(chǎng)數(shù)據(jù)扣除本底溫度,得到核電溫排水形成的溫度場(chǎng)熱影響分布范圍。根據(jù)實(shí)際情況,采用高于本底水溫溫度的0.1 ℃,0.5 ℃,1.0 ℃,2.0 ℃,3.0 ℃和4.0℃溫度,劃分出6個(gè)等級(jí),并分別進(jìn)行編碼(圖8)。根據(jù)各級(jí)水溫像元分布情況和數(shù)量,計(jì)算不同溫升級(jí)別面積(圖9)。



圖8 溫升編碼圖

圖9 不同時(shí)相溫升面積對(duì)比Fig.9 Contrast chart of temperature rising area at different time
溫升編碼圖(圖8)顯示,漲、落潮溫升分布特征與溫度場(chǎng)分布特征一致: 落潮時(shí),溫排水整體向東南方向展布,夏季最大溫升影響范圍為5.95 km,冬季最大溫升影響范圍為5.99 km; 漲潮時(shí)主要向西北方向展布,隨著潮流方向差異,向東北方向偏轉(zhuǎn),夏季最大影響范圍為5.80 km,冬季最大影響范圍為2.38 km。
溫升面積統(tǒng)計(jì)結(jié)果顯示(圖9),相同季節(jié),相同工況條件下,落潮時(shí)刻各級(jí)別的溫升范圍均比漲潮時(shí)刻要大; 相似潮態(tài)、相同工況條件下,冬季溫升分布范圍總體小于夏季。
綜合分析認(rèn)為,漲潮時(shí),外海水體隨漲潮流流入溫排水影響海域,與溫排水快速混合,導(dǎo)致溫排水溫度快速降低,溫排水隨潮水繼續(xù)向西北方向流動(dòng),海水進(jìn)一步混合,溫度繼續(xù)降低; 落潮時(shí)刻,水體流速相對(duì)較慢,靠近排水口處,高溫?zé)崴磁c低溫海水充分混合,仍保持向南流動(dòng)特征。遠(yuǎn)離排水口,溫排水總體隨著落潮主流向東南方向泄出,因潮流的拖曳作用,排水口附近的高溫水體更容易向外擴(kuò)散,因?yàn)榛旌纤肯鄬?duì)漲潮較少,流速相對(duì)較慢,造成溫升分布面積相對(duì)較大,也更易出現(xiàn)4 ℃以上溫升。
因核電所在的興化灣海域,冬夏季潮差、流速均相似[25],但夏季溫度明顯比冬季高,不利于溫排水溫度擴(kuò)散,造成夏季溫升范圍總體比冬季要大。
受排水量、具體潮汐狀態(tài)、氣象條件等因素影響,核電廠溫排水對(duì)海溫的影響在不斷變化。2019年7月15日無人機(jī)航飛數(shù)據(jù)因2臺(tái)機(jī)組降功率運(yùn)行,溫升分布比2019年7月29日衛(wèi)星數(shù)據(jù)溫升范圍要小。無人機(jī)溫升監(jiān)測(cè)結(jié)果顯示,北側(cè)存在殘留的小面積紊動(dòng)熱團(tuán),東南側(cè)低值溫升呈現(xiàn)出一定的分岔現(xiàn)象和不連續(xù)帶狀分布,監(jiān)測(cè)到相對(duì)明顯的4 ℃以上溫升(圖8(e))。
Landsat8衛(wèi)星遙感數(shù)據(jù)具有固定時(shí)間分辨率(16 d)和空間分辨率(100 m),可瞬時(shí)獲取研究海域數(shù)據(jù)(30 s以內(nèi))。與其相比,無人機(jī)航測(cè)更加機(jī)動(dòng)靈活,可根據(jù)氣象、機(jī)組工況、空域情況,調(diào)整數(shù)據(jù)獲取時(shí)間,且空間分辨率能達(dá)到亞米級(jí),但獲取時(shí)間相對(duì)較長(zhǎng)(1 h以上)。由此可以看出,無人機(jī)可在特定時(shí)段對(duì)典型潮態(tài)下溫排水分布特征開展監(jiān)測(cè),且空間分辨率0.5 m的熱紅外數(shù)據(jù),相比于空間分辨率100 m的Landsat8衛(wèi)星數(shù)據(jù),可更清晰反映溫度場(chǎng)細(xì)節(jié)問題。但因單個(gè)照片覆蓋寬度有限,數(shù)據(jù)獲取無法瞬時(shí)完成,溫度場(chǎng)整體性較差,溫排水影響區(qū)外數(shù)據(jù)海域溫度存在不一致的情況,表現(xiàn)為溫度條帶現(xiàn)象。
因無人機(jī)航測(cè)數(shù)據(jù)和Landsat8衛(wèi)星數(shù)據(jù),獲取時(shí)氣象、潮汐狀態(tài)存在不一致的情況,兩者之間對(duì)比不可避免受到干擾。為更直觀反映兩者之間關(guān)系,客觀開展對(duì)比分析,下一步工作應(yīng)在衛(wèi)星過境前后,同時(shí)開展無人機(jī)航測(cè)和近同步的海上測(cè)量工作。
1)溫度反演過程中,根據(jù)水體近似黑體的特性,海水比輻射率取定值0.995,但海水的比輻射率會(huì)隨泥沙含量、海浪狀況等條件發(fā)生改變。相應(yīng)的研究表明,平靜海面的比輻射率可能會(huì)有所下降[8]。因此,比輻射率的取值在一定程度上會(huì)影響到反演的溫度精度。
2)反演值與測(cè)量值擬合產(chǎn)生的誤差除與測(cè)量設(shè)備、測(cè)量操作有關(guān)外,還與實(shí)測(cè)時(shí)間和遙感數(shù)據(jù)獲取時(shí)間不完全一致存在關(guān)系。實(shí)測(cè)數(shù)據(jù)在衛(wèi)星過境前后30 min內(nèi)或無人機(jī)測(cè)量期間獲取,而衛(wèi)星數(shù)據(jù)在不足1 min內(nèi)獲取,無人機(jī)數(shù)據(jù)亦無法與實(shí)測(cè)數(shù)據(jù)完全同步,造成該誤差的存在,尤其在排水口附近因水動(dòng)力條件復(fù)雜,誤差相應(yīng)也會(huì)更大。
1)通過對(duì)反演結(jié)果與海上測(cè)量數(shù)據(jù)的回歸分析,證明Landsat8數(shù)據(jù)利用輻射傳輸方程算法進(jìn)行溫度反演和測(cè)溫型機(jī)芯獲取的溫度結(jié)果可靠。由于測(cè)量不完全同步,2種數(shù)據(jù)之間存在著一定誤差。
2)不同時(shí)相Landsat8遙感數(shù)據(jù)顯示,核電運(yùn)行前,附近海域溫度較為均勻,除自然增溫外,無明顯溫度分異現(xiàn)象; 核電運(yùn)行后,受核電溫排水影響,附近海域具有明顯溫升現(xiàn)象,且不同季節(jié)和不同潮汐條件下,溫升分布特征不同。
3)福清核電廠溫排水溫升在各個(gè)方向的分布受潮態(tài)控制。落潮時(shí),溫升主要集中在排水口東南側(cè)。漲潮時(shí),溫排水總體向西北方向分布。落潮時(shí)刻的溫升范圍要大于漲潮,夏季溫升分布大于冬季同潮態(tài)溫升分布。
4)因受限于飛行效率,無人機(jī)技術(shù)在核電溫排水監(jiān)測(cè)領(lǐng)域的應(yīng)用受到一定制約,但其可以實(shí)現(xiàn)對(duì)溫排水的實(shí)時(shí)監(jiān)測(cè),識(shí)別溫排水分布的細(xì)節(jié)。