王永立

（1.中國(guó)地質(zhì)大學(xué)（北京），北京 100083；2.天津市地質(zhì)調(diào)查研究院，天津 300191）

0 引言

淺層地溫能是指在太陽(yáng)能和地?zé)崃骶C合作用下，存在于地球表層恒溫帶至200m埋深中土壤、巖石和地下水的低溫地?zé)崮埽ㄒ话銣囟鹊陀?5℃）。20世紀(jì)70年代，地?zé)崮芤咽翘旖蚴械难a(bǔ)充能源。淺層地溫能分布于地球表層，是可再生綠色低溫能源。

在基巖裸露地區(qū)，可通過(guò)地質(zhì)調(diào)查推斷淺層地溫能的賦存情況；在第四紀(jì)地層覆蓋地區(qū)，通過(guò)物探、鉆探等方法進(jìn)行解釋推斷。熱紅外遙感檢測(cè)可在短時(shí)間內(nèi)獲取大面積地面溫度場(chǎng)信息[3]，受地面條件限制小，可以有效應(yīng)用于淺層地溫能的調(diào)查。

1 研究區(qū)概況

研究區(qū)位于華北平原東部，北緯38°40′～39°00′，東經(jīng)117°20′～118°00′，陸域面積積約2270km2。熱分布面積800 km2，熱水總儲(chǔ)量200×108t。2010年濱海新區(qū)淺層地溫能利用面積為112.09×104m2。

本區(qū)共圈定了9個(gè)地?zé)岙惓^(qū)。在200m深度的溫度21～28°C，非地?zé)岙惓^(qū)為17～19°C。

第四系地下水淺層地溫能可劃分為4個(gè)含水組。第Ⅰ含水組為全新統(tǒng)和上更新統(tǒng)，埋深一般70m；第Ⅱ含水組為中更新統(tǒng)，底界埋深180～220m；第Ⅲ含水組為下更新統(tǒng)上部，底界埋深290～310m；第Ⅳ含水組為下更新統(tǒng)，底界埋深370～430m（圖1）。

圖1 天津?yàn)I海新區(qū)地表溫度與地?zé)岙惓７植紙D

2 地表溫度反演

目前針對(duì)Landsat 熱紅外波段數(shù)據(jù)進(jìn)行地表溫度反演的算法較多，覃志豪等人根據(jù)地表熱輻射傳導(dǎo)方程，推導(dǎo)出一個(gè)精度較高的演算方法，把大氣和地表的影響直接包括在演算公式中。由于該算法用于僅有一個(gè)熱紅外通道的遙感數(shù)據(jù)，故稱為單窗算法。該算法的優(yōu)點(diǎn)僅需要用地表輻射率、大氣透射率和大氣平均溫度3個(gè)參數(shù)進(jìn)行地表溫度的演算[4][5]。

式中，Ts是地表實(shí)際溫度；Tb是行星亮度溫度；Ta是大氣等效溫度；a 和b 是參考系數(shù)，分別為-67.355351,0.458606；ε為地表輻射率；τ大氣透過(guò)率。

本次遙感圖像采用兩景美國(guó)陸地衛(wèi)星Landsat7 ETM+數(shù)據(jù)，條帶號(hào)122，行編號(hào)33。獲取時(shí)間為2001年5月12日和2002年2月8日，包括與兩期影像相應(yīng)的大氣探空數(shù)據(jù)、地表氣溫、相對(duì)濕度等資料。

2.1 輻射校正

由于地表熱輻射在其傳導(dǎo)過(guò)程中受到大氣和輻射面的多重影響，Landsat傳感器所觀測(cè)到的熱輻射強(qiáng)度已經(jīng)轉(zhuǎn)化為相應(yīng)的灰度值，使直接利用Landsat的原數(shù)據(jù)進(jìn)行區(qū)域分析所得到的結(jié)論存在很大的偏差[5]。通常Landsat ETM+6數(shù)據(jù)是以灰度值(DN值)來(lái)表示，DN值越大，亮度越大，表示地表熱輻射強(qiáng)度越大，溫度越高。ETM+6數(shù)據(jù)中求算亮度溫度的過(guò)程包括把DN值轉(zhuǎn)化為相應(yīng)的熱輻射強(qiáng)度值，然后推算所對(duì)應(yīng)的地表輻射溫度。陸地衛(wèi)星傳感器在設(shè)計(jì)時(shí)就考慮到如何把所接收到的輻射強(qiáng)度轉(zhuǎn)化為DN值。對(duì)于ETM+數(shù)據(jù)，所接收到的輻射強(qiáng)度與其DN值存在如下關(guān)系：

式中Lλ為TM傳感器接收到的輻射亮度,單位為w?m-2?Sr-1?μm-1，Lmax和Lmin為TM遙感器所接收到的最大和最小輻射強(qiáng)度。當(dāng)ETM+6低增益時(shí)，Lmax=17.04，Lmin=0.0；當(dāng)高增益時(shí)，Lmax=12.65，Lmin=3.2。然后根據(jù)輻射亮度推算衛(wèi)星高度傳感器上所探測(cè)到的像元亮度溫度Tb，其公式為:

2.2 大氣等效溫度及大氣透射率

大氣透射率主要取決于大氣水分含量，而大氣等效溫度取決于大氣剖面的氣溫分布和大氣水分狀況。影響大氣透射率的因素較多，氣壓、氣溫、氣溶膠含量、大氣水分含量、CO2等對(duì)熱輻射傳導(dǎo)均有不同程度的作用。因此求算確定大氣透射率比較復(fù)雜，需要較詳細(xì)的大氣剖面數(shù)據(jù)。一般來(lái)說(shuō)，準(zhǔn)確氣透射率需要進(jìn)行大氣模擬，目前較普遍使用的大氣模擬程序有LOWTRAN、MODTRAN和6S等。這種大氣模擬需要很詳細(xì)的大氣剖面數(shù)據(jù)；而大多數(shù)情況下不具備這樣的條件，從而使大氣模擬法難以實(shí)施。

大氣等效溫度主要取決于大氣剖面氣溫分布和大氣狀態(tài)。由于衛(wèi)星飛空的時(shí)間很短，一般直接利用覃志豪總結(jié)出的大氣等效溫度Ta的估計(jì)方程[6]:

中緯度夏季平均大氣作用溫度:

中緯度冬季平均大氣作用溫度:

其中，To為近地面(一般為2m處)的氣溫，To和Ta的單位均為K。這些，在標(biāo)準(zhǔn)大氣狀態(tài)下(天空晴朗、沒(méi)有渦旋作用)，大氣平均溫度是地面線性函數(shù)。

大氣透射率的變化主要取決于大氣水分含量的動(dòng)態(tài)變化，其動(dòng)態(tài)變化不大對(duì)大氣透射率的變化沒(méi)有顯著影響[6]。因此，大氣水分含量就成為透射率估計(jì)的主要因素。在高溫和低溫不同的溫度季節(jié)，其大氣透射率估計(jì)方程不同（表1）。

根據(jù)楊景梅[7]等研究，大氣水汽含量可以通過(guò)與地面水汽壓之間的關(guān)系確定。

其中，e是絕對(duì)水汽壓，可以表示為：

表1 大氣透射率估計(jì)方程

式中，RH是相對(duì)濕度，T0是氣溫，a、b是常數(shù)分別為17.67K、29.65K。

2.3 地表比輻射率確定

地表比輻射率是單窗算法的一個(gè)重要參數(shù)[8]。通過(guò)歸一化植被指數(shù)(NDVI)獲得地表比輻射率是在沒(méi)有實(shí)時(shí)的參考數(shù)據(jù)下一種有效方法。

Enric-Valor等[9]提出了一種使用NDVI計(jì)算植被覆蓋率，然后推導(dǎo)計(jì)算地表比輻射率的方法：

式中，P v表示植被覆蓋度，取N DV Iv=0.70，NDVIs=0.05。當(dāng)某個(gè)像元的NDVI值大于0.70時(shí)，Pv取值為1；當(dāng)NDVI小于0.05時(shí)，Pv取值為0。ρNIR、ρR分別表示傳感器的近紅外波段與紅波段的地面反射率值。

3 結(jié)果分析

通過(guò)上述計(jì)算得出地表亮度值、大氣等效溫度、大氣透射率和比輻射率，進(jìn)而求得地表實(shí)際溫度。將研究區(qū)的斷裂構(gòu)造、地?zé)岙惓Ｔ贕IS軟件中進(jìn)行疊加得到地表溫度與地?zé)岙惓７植紙D（圖1）。

通過(guò)對(duì)比分析可以看出，這些已探明的橋沽、山嶺子、萬(wàn)家碼頭和沙井子等地?zé)岙惓^(qū)在熱紅外遙感影像上都顯示出較高的溫度值，而且大都沿著活動(dòng)斷裂的控制方向，特別是在活動(dòng)斷裂的交叉點(diǎn)，溫度異常尤為明顯。而斷裂為地下熱能向上傳輸提供了便利通道，說(shuō)明地表溫度異常可能是由地下熱流通道斷裂上升到淺部所引起的。

分布在北部的看才莊地?zé)岙惓^(qū)，屬于于黃驊凹陷的寧河凸起構(gòu)造單元，地溫梯度變化區(qū)間3.5～6.5℃/100m，NE向的楊家泊斷裂構(gòu)成北部控?zé)徇吔?，?nèi)部為近SN向的楊家泊東斷裂縱穿，而到南部則被晚更新世活動(dòng)的漢沽斷裂所橫切。在2002年冬季遙感圖像上，異常區(qū)北部呈現(xiàn)地表溫度高值，普遍高于背景值。

西北部的橋沽異常區(qū)，位于黃驊坳陷的寧河凸起和北塘坳陷的接觸部位，被近EW向漢沽斷裂和NE向的楊家泊斷裂、大田莊斷裂控制。在兩期遙感圖像上，除了薊運(yùn)河地表溫度呈現(xiàn)低值外，其它區(qū)域都高于背景值。

在西部的山嶺子地?zé)岙惓^(qū)，位于滄縣隆起的潘莊凸起構(gòu)造單元，地溫梯度變化區(qū)間3.5～7.0°C/100m，由滄東斷裂、大寺斷裂、海河斷裂和張貴莊斷裂幾組斷裂所控制。在兩期圖像上地表溫度都呈現(xiàn)高值。

西南部的萬(wàn)家碼頭地?zé)岙惓^(qū)，位于滄縣隆起的小韓莊凸起構(gòu)造單元，地溫梯度變化區(qū)間3.5～7.0℃/100m，在區(qū)內(nèi)受到滄東斷裂、小韓莊斷裂和上古林?jǐn)嗔训目刂?。在兩期遙感圖像上，都不同程度的出現(xiàn)高值區(qū)，在板橋農(nóng)場(chǎng)、中塘鎮(zhèn)和大臺(tái)子出現(xiàn)極值；受獨(dú)流減河和北大港水庫(kù)的影響，地表溫度在異常區(qū)東部呈現(xiàn)低值。

南部的沙井子地?zé)岙惓^(qū)，處于黃驊坳陷的港西凸起構(gòu)造單元，受港東斷裂帶的控制。在遙感圖像上，地表溫度普遍不高，僅在太平村呈現(xiàn)高值，推測(cè)受到北大港水庫(kù)及周邊河流水系的影響。

4 結(jié)語(yǔ)

地表溫度與地下地?zé)岙惓＞哂幸欢ǖ南嚓P(guān)性。建議在地表溫度高值區(qū)采用物探和鉆探等方法開(kāi)展地?zé)嵴{(diào)查。

天津?yàn)I海新區(qū)道路網(wǎng)密集，地表土層主要由近代海侵海退過(guò)程中之海水浸漬而成，具有較高的鹽分，在遙感圖像上顯示較高的反射率，反演出地表溫度較高。因此，必須選用其它手段剔除非淺層地溫能引發(fā)地表溫度升高等干擾因素，有效發(fā)揮遙感信息的作用。

[1]朱家玲等.地?zé)崮荛_(kāi)發(fā)與應(yīng)用技術(shù)[M].北京:化學(xué)工業(yè)出版社,2006.

[2]汪 集,馬偉斌,龔宇烈.地?zé)崂眉夹g(shù)—可再生能源叢書,北京:中國(guó)化學(xué)工業(yè)出版社,2006.

[3]李小文,汪駿發(fā),王錦地等，多角度與熱紅外對(duì)地遙感[M].北京:科學(xué)出版社,2001.

[4]覃志豪,Zhang Minghua, Arnon, Arnon Karnieli, Pedro Berliner等.用陸地衛(wèi)星TM6數(shù)據(jù)演算地表溫度的單窗算法[J].地理學(xué)報(bào).2001,56(4):456～466.

[5]Qin Z, Karnieli A, Berliner P.A mono-window algorithm for retrieving land surface temperature from Landsat TM data and its application to the Israel-Egypt border region[J].International Journal of Remote Sensing,2001,22(18):3719～3746.

[6]覃志豪,Li Wenjuan, Zhang Minghua etc.單窗算法的大氣參數(shù)估計(jì)方法[J].國(guó)土資源遙感.2003,56(2):37～43.

[7]楊景梅,邱金恒.我國(guó)可降水量同地面水汽壓關(guān)系的經(jīng)驗(yàn)表達(dá)式[J].大氣科學(xué).1996,20(5):620～626.

[8]張兆明,何國(guó)金,肖榮波等.基于RS與GIS的北京市熱島研究[J].地球科學(xué)與環(huán)境學(xué)報(bào).2007,29(1):107～109.

[9]Enric-Valor,Vicente-Caselles.從NDVI估算地表發(fā)射率-以歐洲、非洲和南美洲為例[J].環(huán)境遙感.1996(57):167～184.

[10]Enric-Valor,Vicente-Caselles.Mapping land surface emissivity from NDVI: application to European, African and south American areas[J].Remote sensing of environm ent,1996(57):167～184.