Landsat ETM+遙感數(shù)據(jù)在尋找臨沂地?zé)嶂械膽?yīng)用

劉效才，江泳，劉孝陽(yáng)，劉琨

(1.山東省第七地質(zhì)礦產(chǎn)勘察院，山東 臨沂　276000；2.臨沂市國(guó)土資源局，山東 臨沂　276000；3.中國(guó)地質(zhì)大學(xué)(北京)土地科學(xué)技術(shù)學(xué)院，北京　100083)

Landsat ETM+遙感數(shù)據(jù)在尋找臨沂地?zé)嶂械膽?yīng)用

劉效才1，江泳2，劉孝陽(yáng)3，劉琨1

(1.山東省第七地質(zhì)礦產(chǎn)勘察院，山東 臨沂276000；2.臨沂市國(guó)土資源局，山東 臨沂276000；3.中國(guó)地質(zhì)大學(xué)(北京)土地科學(xué)技術(shù)學(xué)院，北京100083)

沂沭斷裂帶內(nèi)蘊(yùn)藏著豐富的地?zé)豳Y源，以landsat ETM+為主要遙感數(shù)據(jù)，利用輻射傳輸方程反演出整個(gè)臨沂地區(qū)地表溫度，對(duì)灰色溫度圖像進(jìn)行交互式拉伸突顯出33℃及以上溫度異常區(qū)，結(jié)合ETM+真彩色圖像以及Google影像，排除城鎮(zhèn)建筑、農(nóng)村居民點(diǎn)、河岸灘、裸露巖體等異常區(qū)域，將所得異常結(jié)果與地?zé)岬刭|(zhì)背景圖進(jìn)行疊加比對(duì)，并從區(qū)域地質(zhì)角度對(duì)其進(jìn)行分析，以期找到地?zé)釢摲鼌^(qū)。

Landsat ETM+；溫度反演；地?zé)?；臨沂

引文格式：劉效才，江泳，劉孝陽(yáng)，等.Landsat ETM+遙感數(shù)據(jù)在尋找臨沂地?zé)嶂械膽?yīng)用[J].山東國(guó)土資源，2015，31(5)：76-81.LIU Xiaocai，LIU Xiaoyang, JIANG Yong，etc. Application of Landsat ETM+ Remote Sensing Data in Searching Geothermal Field in Linyi City[J].Shandong Land and Resources,2015，31(5)：76-81.

0　引言

地?zé)崾翘N(yùn)藏在地球內(nèi)部的一種自然熱能，其熱量主要來(lái)自地球內(nèi)部放射性元素的蛻變。從能源的可持續(xù)發(fā)展戰(zhàn)略看，地?zé)豳Y源是一種發(fā)展前景非常好的綠色能源，它分布廣泛而且性質(zhì)較穩(wěn)定，與煤、石油、天然氣等傳統(tǒng)能源相比，具有廉價(jià)、潔凈和可持續(xù)利用等優(yōu)點(diǎn)。據(jù)估算，地球蘊(yùn)藏的地?zé)崮芗s14.5×10J，折合4948萬(wàn)億t標(biāo)準(zhǔn)煤，而全球煤的埋藏量?jī)H為1萬(wàn)億t[1]。加上近年來(lái)全球氣候變暖，霧霾天氣頻發(fā)等環(huán)境污染問(wèn)題日益突出，地?zé)豳Y源開發(fā)利用的巨大前景備受各國(guó)重視?？茖W(xué)開發(fā)利用地?zé)崮苣軌驅(qū)崿F(xiàn)經(jīng)濟(jì)效益社會(huì)效益和環(huán)境效益的統(tǒng)一，有著各方面的重要意義[2]。

比較成熟的地?zé)豳Y源的探測(cè)方法有地質(zhì)學(xué)方法、地球物理勘探方法以及地?zé)岙惓，F(xiàn)象調(diào)查法、地球化學(xué)法等[3]。但這些傳統(tǒng)方法周期長(zhǎng)，資金投入大，如果對(duì)探測(cè)區(qū)地?zé)岷暧^環(huán)境條件不明，就會(huì)出現(xiàn)投資風(fēng)險(xiǎn)大或盲目開采的情況，對(duì)地?zé)豳Y源的可持續(xù)發(fā)展是一種很大的危害[4]。隨著科學(xué)技術(shù)的不斷發(fā)展，遙感技術(shù)的不斷成熟，其在地?zé)豳Y源探測(cè)與調(diào)查中的優(yōu)越性越來(lái)越明顯：可快速、經(jīng)濟(jì)、大區(qū)域的定量研究熱量資源和熱場(chǎng)分布變化。郝艷梅等于2001年探索了遙感技術(shù)在小區(qū)域應(yīng)用的方法和途徑,揭示了遙感技術(shù)的廣闊應(yīng)用前景[5]；喬玉良等于2002年利用TM6作為數(shù)據(jù)源對(duì)快速尋找地下熱水進(jìn)行了探究[6]；楊波等于2003年利用遙感技術(shù)發(fā)現(xiàn)了騰沖西南地區(qū)地?zé)岣患囊?guī)律性，提出了地?zé)崽锏倪b感影像模式[7]；章方等于2012年利用Landsat 7 ETM+遙感影像對(duì)長(zhǎng)春地區(qū)地表溫度異常進(jìn)行了遙感分析[8]。

1　研究區(qū)域與數(shù)據(jù)資料

1.1研究區(qū)域概況

臨沂市近黃海，東連日照，西接棗莊、濟(jì)寧、泰安，北靠淄博、濰坊，南鄰江蘇。南北最大長(zhǎng)距228km，東西最大寬度161km，總面積17191.2km2，地勢(shì)西北高東南低，是山東省面積最大的地級(jí)市。臨沂市地質(zhì)構(gòu)造復(fù)雜，地層發(fā)育比較齊全，從太古宇至新生界，除上奧陶統(tǒng)、志留系、泥盆系、下石炭統(tǒng)、三疊系及中、下侏羅統(tǒng)、古新統(tǒng)、新近系缺失以外，其他各期地層都有發(fā)育。主要構(gòu)造以郯廬斷裂帶(境內(nèi)稱沂沭斷裂帶)為主。沂沭斷裂帶經(jīng)郯城北向延伸，縱貫臨沂市，以斷裂為界，臨西為魯西臺(tái)背斜的一部分，屬華北地臺(tái)范疇，臨東為膠南隆起的一部分，屬揚(yáng)子大陸塊范疇。

沂沭斷裂帶由4條主干斷裂組成，自西向東依次為鄌郚-葛溝斷裂、沂水-湯頭斷裂、安丘-莒縣斷裂及昌邑-大店斷裂[9](圖1)。該斷裂帶不但控制了地層展布，巖漿活動(dòng)也沿?cái)嗔亚秩?，同時(shí)把深部熱液沿?cái)嗔严驕\部運(yùn)移，為地?zé)嵝纬商峁┝送ǖ篮蛢?chǔ)存空間[10]。復(fù)雜的地質(zhì)構(gòu)造,為區(qū)內(nèi)地?zé)豳Y源形成提供了地質(zhì)條件,是地?zé)豳Y源較豐富地區(qū)[11]，臨沂也因此被國(guó)土資源部授予“中國(guó)溫泉之城”的美稱。目前已發(fā)現(xiàn)多處溫泉，如湯頭地?zé)崽?、銅井地?zé)崽?、松山地?zé)崽锏?。由于近年?lái)勘查工作一直未能取得突破性進(jìn)展，加上溫泉長(zhǎng)期過(guò)量開采，多處溫泉出現(xiàn)不同程度的枯竭，因此方向正確的溫泉勘查工作迫在眉睫。

1—第四系；2—新近系；3—古近系；4—白堊系；5—奧陶系；6—寒武系；7—沂水巖群；8—花崗閃長(zhǎng)巖；9—閃長(zhǎng)玢巖；10—二長(zhǎng)花崗巖；11—韌性剪切帶；F1—甘霖?cái)嗔?；F2—燕甘斷裂；F3—蒙山斷裂；F4—新泰-垛莊斷裂；F5—銅冶店-孫祖斷裂；F6—相邸斷裂圖1　臨沂市地?zé)岬刭|(zhì)背景圖鞏貴仁等修改)

1.2數(shù)據(jù)資料

該文所用遙感數(shù)據(jù)為2003年5月11日的Landsat7 ETM+遙感影像，衛(wèi)星過(guò)境時(shí)間為上午10點(diǎn)31分，圖像質(zhì)量較好，清晰無(wú)云。所用氣象統(tǒng)計(jì)資料來(lái)源于中國(guó)氣象局。

2　地表溫度反演

2.1輻射定標(biāo)

ETM+數(shù)據(jù)是以灰度值來(lái)表示的，灰度值并非實(shí)際的反射或者輻射強(qiáng)度值，但它們之間存在著一定的換算關(guān)系，在進(jìn)行溫度反演之前，必須要將灰度值轉(zhuǎn)化為傳感器接受到的輻射強(qiáng)度值，即所謂的輻射定標(biāo)。對(duì)于ETM+數(shù)據(jù)，輻射定標(biāo)的轉(zhuǎn)換公式可以簡(jiǎn)化為

L(λ)=gain×QDN+offset

式中：QDN為像元灰度值；gain和offset分別為圖像的增益和偏置，是數(shù)據(jù)頭文件中提供的定標(biāo)系數(shù)。

2.2地表比輻射率計(jì)算

物體的比輻射率是物體向外輻射電磁波的能力表征，這一基本參數(shù)主要取決于地表物質(zhì)結(jié)構(gòu)和觀測(cè)波段，與觀測(cè)角度等因素也存在關(guān)系。ETM+熱紅外波段分辨率為60m，大都為混合像元，在大尺度上對(duì)比輻射率精確測(cè)量的難度很大，目前一般采用間接的方法來(lái)獲取。地球表面不同區(qū)域的地表結(jié)構(gòu)雖然很復(fù)雜，但從衛(wèi)星像元的尺度來(lái)看,可以大體視作由3種類型構(gòu)成：水面、城鎮(zhèn)和自然表面[12]。根據(jù)VandeGriend和Owe的實(shí)測(cè)數(shù)據(jù)研究，ε與植被覆蓋度存在密切關(guān)系。該文將水體像元的比輻射率賦值為0.995，自然表面和城鎮(zhèn)像元的比輻射率估算則分別根據(jù)下式進(jìn)行計(jì)算：

式中：εsurface和εbuilding分別代表自然表面像元和城鎮(zhèn)像元的比輻射率；Fv代表植被覆蓋度，指植被冠層的垂直投影面積與土壤總面積之比。采用的是混合像元分解法，將整景影像的地類大致分為水體、植被和建筑，具體的計(jì)算公式為：

FV=(NDVI- NDVIS)/(NDVIV-NDVIS)

式中：NDVI為歸一化差異植被指數(shù)，NDVIv為100%植被覆蓋區(qū)域的NDVI值，對(duì)應(yīng)研究區(qū)域NDVI最大值；NDVIs為100%裸土覆蓋區(qū)域的NDVI值，對(duì)應(yīng)研究區(qū)域NDVI最小值。然而影像不可避免地會(huì)受到噪聲的影響，可能對(duì)NDVI值產(chǎn)生一定的影響，為了保證估算精度，該文利用不同波段彩色合成分別識(shí)別出植被與土壤，經(jīng)過(guò)統(tǒng)計(jì)分析，得到NDVIv=0.5，NDVIs=0。且有，當(dāng)某個(gè)像元的NDVI大于0.5時(shí)，F(xiàn)V取值為1；當(dāng)NDVI小于0時(shí)，F(xiàn)V取值為0。

2.3計(jì)算相同溫度下黑體的輻射亮度值

在熱紅外窗口，大氣并非是完全透明的。傳感器接收的熱紅外輻射除地表信息以外，還受大氣狀況的影響。由于熱紅外波段大氣散射一般可以忽略不計(jì)，熱紅外輻射的大氣傳輸方程應(yīng)考慮大氣的吸收和發(fā)射[13]。Chandrasekhar[14]把傳感器接收的熱紅外輻射表達(dá)為3項(xiàng)之和，即大氣上行輻射項(xiàng)，經(jīng)大氣吸收衰減后到達(dá)傳感器的輻射，大氣下行輻射項(xiàng)經(jīng)地表反射后通過(guò)大氣傳輸?shù)絺鞲衅鞯妮椛?。?jiǎn)化Chandrasekhar的計(jì)算表達(dá)式得到輻射傳輸方程：

式中：ε為地表輻射率；TS為地表真實(shí)溫度；B(TS)為普朗克定律推到得到的黑體在TS的熱輻射亮度；τ為大氣在熱紅外波段的透過(guò)率。則溫度為T的黑體在熱紅外波段的輻射亮度B(TS)為：

根據(jù)遙感圖像時(shí)間以及中央經(jīng)緯度等信息，在NASA官網(wǎng)模型中查詢式中所需參數(shù)值。

2.4反演地表溫度

在獲取溫度為TS的黑體在熱紅外波段的輻射亮度后，根據(jù)普朗克公式的反函數(shù)，求得地表真實(shí)溫度TS=K2/ln(K1/B(TS)+1)。對(duì)于ETM+數(shù)據(jù)，K1=666.09W/(m2·sr·μm)，K2=1282.71K。溫度反演結(jié)果見圖2。

圖2　臨沂市地表溫度圖

2.5反演地表溫度驗(yàn)證

該文所用輻射傳輸方程法反演出臨沂地區(qū)最高溫度38.17℃，最低溫度16.58℃，大部分地區(qū)溫度介于20～30℃之間。由于沒有衛(wèi)星過(guò)境時(shí)的同步氣象觀測(cè)數(shù)據(jù)，所以不能對(duì)地表溫度反演的結(jié)果進(jìn)行客觀的評(píng)價(jià)[15]。該文只是采用中國(guó)氣象局統(tǒng)計(jì)的往年5月份氣溫資料進(jìn)行簡(jiǎn)單比對(duì)。據(jù)統(tǒng)計(jì)，臨沂地區(qū)5月份溫度極端最低溫度為3℃(1961年)，極端最高溫度為38℃(1982年)。從圖3可以看出，平均最高溫度約為25℃，平均最低溫度約為14℃。由于陸地與大氣熱交換的不平衡性，一般情況下認(rèn)為白天陸地溫度會(huì)略高于大氣溫度，水域溫度略低于大氣溫度，與所反演結(jié)果大體吻合。從溫度空間分布來(lái)看，河流、水庫(kù)區(qū)域溫度較低，北部山區(qū)和城鎮(zhèn)地區(qū)溫度較高，溫度分布符合邏輯，所以認(rèn)為反演結(jié)果準(zhǔn)確。

圖3　臨沂市5月份溫度走勢(shì)圖

3　溫度異常區(qū)分析

3.1溫度異常區(qū)提取

為了便于提取溫度異常信息，該文對(duì)灰度值溫度圖像進(jìn)行交互式拉伸，突顯33℃以上異常區(qū)，并用ENVI軟件與ETM+真彩色圖像進(jìn)行地理鏈接。由于ETM影響分辨率為30m，不能很細(xì)致地觀察地表成分，只能對(duì)溫度異常區(qū)進(jìn)行粗篩選，排除一些顯而易見的人為或地形地貌因素干擾。為了更好地對(duì)異常區(qū)進(jìn)行分析，該文將交互式拉伸處理后的圖像加載到Global mapper中，對(duì)溫度異常區(qū)進(jìn)行數(shù)字化，圈出共計(jì)58個(gè)區(qū)域，并導(dǎo)出為KMZ矢量格式，利用Google Earth高分辨率和歷史圖像等優(yōu)勢(shì)，對(duì)其進(jìn)一步比較分析。

經(jīng)認(rèn)真比對(duì)分析，該文將溫度異常區(qū)劃分為兩類：A類溫度異常區(qū)和B類溫度異常區(qū)。A類溫度異常區(qū)指由于城鎮(zhèn)建筑、村莊、河岸灘、裸露山體等因素影響下呈現(xiàn)的溫度異常區(qū)，這種溫度異常并非地?zé)嵋蛩貙?dǎo)致。B類溫度異常區(qū)排除人為或地形地貌等因素對(duì)溫度的影響，溫度異常原因尚不明確，待考證。兩者共計(jì)58處，其中A類溫度異常區(qū)39處，B類溫度異常區(qū)19處，對(duì)B類溫度異常區(qū)進(jìn)行編號(hào)B1-B19。分析結(jié)果如圖4所示。

圖4　溫度異常區(qū)在Google Earth上的顯示

3.2B類溫度異常區(qū)與臨沂市地?zé)豳Y源開發(fā)利用規(guī)劃圖的疊加分析

圖5為臨沂市政府所編制下發(fā)《中國(guó)溫泉之城——山東臨沂發(fā)展建設(shè)規(guī)劃(2013—2020年)》的地?zé)衢_發(fā)利用規(guī)劃略圖，包含2013年前已經(jīng)開發(fā)的溫泉示范區(qū)，2013—2015年溫泉規(guī)劃開發(fā)區(qū)以及2016—2020年溫泉規(guī)劃開發(fā)區(qū)。其中已開發(fā)示范區(qū)已完成勘察驗(yàn)證，規(guī)劃開發(fā)區(qū)為地?zé)岙惓^(qū)。

圖5　臨沂市“一軸兩翼”地?zé)豳Y源開發(fā)利用規(guī)劃略圖

將該文所得B類溫度異常區(qū)與之疊加分析，共發(fā)現(xiàn)6處B類溫度異常區(qū)完全處于已開發(fā)示范區(qū)或規(guī)劃開發(fā)區(qū)范圍內(nèi)，編號(hào)分別為B1，B5，B8，B14，B15，B19，它們顯示出的溫度異常面積大都較小，呈現(xiàn)為點(diǎn)狀分布，應(yīng)為地?zé)崧额^或者地?zé)嶂行膮^(qū)。B4，B10，B12三處溫度異常區(qū)與溫泉已開發(fā)示范區(qū)或規(guī)劃開發(fā)區(qū)有大部分重疊，它們的面積較大，呈現(xiàn)長(zhǎng)條帶狀。其中B10受到北部費(fèi)縣處溫泉和區(qū)內(nèi)山體的共同影響，表現(xiàn)出溫度異常形狀狹長(zhǎng)且面積較大；B12處于城市建筑密度較大區(qū)域，受其與地?zé)岬墓餐绊?，表現(xiàn)出的面積也較大。部分地?zé)釁^(qū)在熱紅外遙感溫度圖像上得到一定程度的體現(xiàn)，表現(xiàn)為熱異常，但由于受到熱源埋深、溫度和地表覆蓋情況等因素的影響，并非所有地?zé)釁^(qū)都會(huì)在地表呈現(xiàn)出熱異常，有些地?zé)釁^(qū)未能監(jiān)測(cè)的到。

3.3B類溫度異常區(qū)與地?zé)岬刭|(zhì)背景圖的疊加分析

排除A類溫度異常區(qū)和上文中所列與溫泉已開發(fā)或規(guī)劃開發(fā)區(qū)重疊部分的異常區(qū)，該文利用圖像疊加層功能，將所剩B類溫度異常區(qū)與臨沂市地?zé)岬刭|(zhì)背景圖在Google Earth中進(jìn)行疊加，從區(qū)域地質(zhì)角度對(duì)其進(jìn)行分析判斷地?zé)岽嬖诘目赡苄?。?jīng)分析，共發(fā)現(xiàn)3處B類溫度異常區(qū)具有存在地?zé)岬目赡苄?，且均處于韌性剪切帶或者斷裂帶上，編號(hào)分別為B2，B3，B9(圖6)。

圖6　溫度異常區(qū)與地質(zhì)圖在Google Earth上的疊加

3.3.1B2號(hào)異常區(qū)

該區(qū)位于蒙陰縣東部興旺莊村—魏家?guī)r路村一帶，呈NW—SE向條帶狀分布，處于銅冶店-孫祖斷裂上。

銅冶店-孫祖斷裂為NW向旋轉(zhuǎn)構(gòu)造帶的臨北帚狀構(gòu)造，全長(zhǎng)100km，走向NW，傾向NE；斷裂帶內(nèi)構(gòu)造片巖、糜棱巖、角礫巖發(fā)育，斷面上斜沖擦痕清晰，斷盤上地層擠壓牽引明顯，不利于地?zé)崴魍?，也是阻水?dāng)嗔选盤(弧外)下降由泰山巖群構(gòu)成，S盤(弧內(nèi))上升由古生界、中生界組成；控制著中、新界凹陷盆地沉積，其特點(diǎn)為北斷南超，地層?xùn)|老西新；由一組多期活動(dòng)的左行逆斷層組成；活動(dòng)時(shí)間為侏羅紀(jì)至新近紀(jì)，白堊紀(jì)活動(dòng)最強(qiáng)烈。該斷裂位于水文地質(zhì)條件良好，富水性強(qiáng)的沂沭斷裂帶以西中低山丘陵水文地質(zhì)區(qū)，該區(qū)域分布大面積的碳酸鹽類裂隙巖溶水，含水層主要由灰?guī)r組成，厚度大，分布廣，裂隙巖溶發(fā)育，彼此連通，為地下水提供了充裕的賦存空間。銅冶店-孫祖斷裂與NNE向斷裂交會(huì)處有非常好的成熱條件，同時(shí)由于鄌郚-葛溝斷裂為由北向南的疏水?dāng)嗔?，故NNE向斷裂與其交會(huì)部位為地下水匯集地區(qū)，地?zé)崴簇S沛。

3.3.2B3號(hào)異常區(qū)

該區(qū)位于沂南縣蘭香埠村—司馬村一帶，呈NNE向條帶狀分布，處于沂水-湯頭斷裂帶上。

沂水-湯頭斷裂走向15°，構(gòu)成蘇村-馬站地塹的東界和汞丹山凸起的西界，斷裂南段出露不好，北段沂水境內(nèi)出露好。在新生代早期構(gòu)造處于開放系統(tǒng)，火山活動(dòng)頻繁，在該斷裂附近發(fā)育一系列派生斷裂，在斷裂交會(huì)處對(duì)地?zé)嵝纬蓸O為有利，湯頭溫泉就在該斷裂帶上。

3.3.3B9號(hào)異常區(qū)

該區(qū)位于費(fèi)縣盛家?guī)X—安樂窩村一帶，呈NE向條帶狀，處在甘霖?cái)嗔烟幐浇?；區(qū)內(nèi)地層多以寒武紀(jì)海相沉積巖為主，在中生代印支期被中酸性花崗巖沿裂隙或順層侵入。推斷區(qū)內(nèi)中生代侵入巖根部及甘霖?cái)嗔烟峁┝说責(zé)豳Y源的熱源及通道，存在較好的地?zé)峥辈榍熬啊?/p>

4　結(jié)論

該文從熱紅外遙感和區(qū)域地質(zhì)角度對(duì)臨沂市地?zé)岱植紶顩r進(jìn)行了簡(jiǎn)要分析，共有部分溫泉已開發(fā)區(qū)或者規(guī)劃開發(fā)區(qū)顯示出了溫度異常，在遙感溫度圖像得到了很好的驗(yàn)證。另外發(fā)現(xiàn)3處地?zé)釢摲鼌^(qū)，對(duì)以后臨沂地區(qū)地?zé)峥睖y(cè)具有一定的導(dǎo)向性作用。在實(shí)際應(yīng)用中，此法與同類技術(shù)相比具有以下優(yōu)勢(shì)：

(1)利用熱紅外波段進(jìn)行地表溫度反演，確切的量化出整個(gè)臨沂地區(qū)地表溫度，而以往遙感地?zé)嵯嚓P(guān)研究往往是通過(guò)波段差值運(yùn)算和增強(qiáng)處理等步驟定性的區(qū)分地表溫度高低，對(duì)于溫度異常區(qū)缺乏定量分析。相比之下，通過(guò)溫度反演這一過(guò)程使得研究更加精確。

(2)對(duì)于灰色溫度圖像進(jìn)行交互式拉伸，凸顯出33℃以上區(qū)域，而使得其他溫度區(qū)作為黑色背景值不予顯示，使得圖像更具有可讀性。

(3)該文將33℃以上溫度異常區(qū)進(jìn)行了矢量化，并利用Google Earth具有高分辨率和歷史圖像等優(yōu)勢(shì)排除人為因素或地形地貌等引起的溫度異常區(qū)，把所得到的結(jié)果分別與地?zé)嵋?guī)劃圖和地?zé)岬刭|(zhì)背景圖進(jìn)行疊加分析，通過(guò)2次篩選增加了結(jié)果的準(zhǔn)確性和可靠性。

但由于缺乏物探、鉆探以及磁異常等實(shí)際工作，對(duì)得出的3處地?zé)釢摲鼌^(qū)也只屬于初探，尚需等待驗(yàn)證。遙感雖快捷方便，可節(jié)省大量人力物力，但畢竟還不能完全代替?zhèn)鹘y(tǒng)的地?zé)峥睖y(cè)技術(shù)。遙感技術(shù)必須與其他技術(shù)方法相配合,與專業(yè)知識(shí)及專家思想相結(jié)合才能取得比較好的應(yīng)用效果。

[1]中國(guó)國(guó)土資源報(bào).國(guó)際視野下的地?zé)豳Y源開發(fā)利用[J].國(guó)土資源，2014，(1):56-57.

[2]吳新雄.科學(xué)、清潔、高效、可持續(xù)地發(fā)展地?zé)崮躘J].能源，2014，(3):38-39.

[3]曹凱,江南.基于TM6的地?zé)豳Y源的熱紅外遙感探查模型研究[J].遙感信息,2006，(2):18-21.

[4]喬玉良,呂芳,王棟,等.衛(wèi)星遙感技術(shù)在尋找地下熱水中的應(yīng)用[J].航天返回與遙感,2007,23(4):58-62.

[5]郝艷梅,高新法,張義文.遙感技術(shù)在下莊地區(qū)地?zé)豳Y源開發(fā)中的應(yīng)用研究[J].國(guó)土資源遙感,2001，(1):19-24.

[6]喬玉良.應(yīng)用紅外遙感尋找地下熱水初探[J].航天返回與遙感,2002,23(1):33-35.

[7]楊波,吳德文,賴健清,等.遙感技術(shù)在騰沖西南地區(qū)地?zé)豳Y源研究預(yù)測(cè)中的應(yīng)用[J].國(guó)土資源遙感,2003，(2):23-26.

[8]章方,王艷江,楊曉紅,等.長(zhǎng)春地區(qū)地表溫度異常遙感分析[J].測(cè)繪與空間地理信息,2012,35(9):97-99.

[9]王志才,賈榮光,孫昭民,等.沂沭斷裂帶安丘-莒縣斷裂安丘-朱里段幾何結(jié)構(gòu)與活動(dòng)特征[J].地震地質(zhì),2005,27(2):212-220.

[10]李小鳳.郯廬斷裂帶[M].北京:地質(zhì)出版社,2000.

[11]鞏貴仁,唐孟武,逯光明,等.遙感技術(shù)在臨沂地區(qū)地?zé)豳Y源勘查中的應(yīng)用[J].山東國(guó)土資源,2008,24(10):17-20.

[12]覃志豪,李文娟,徐斌,等.利用Landsat TM6反演地表溫度所需地表輻射率參數(shù)的估計(jì)方法[J].海洋科學(xué)進(jìn)展,2004,22(B10):129-137.

[13]田國(guó)良,柳欽火.熱紅外遙感(第2版)[M].北京：電子工業(yè)出版社，2014.

[14]Chandrasekhar S. Radiative Transfer. New York: Dover Publ Inc,1960.

[15]宮阿都,江樟焰,李京,等.基于Landsat TM圖像的北京城市地表溫度遙感反演研究[J].遙感信息,2005，(3):18-20.

Application of Landsat ETM+ Remote Sensing Data in Searching Geothermal Field in Linyi City

LIU Xiaocai1，LIU Xiaoyang1, JIANG Yong2, LIU Kun1

(1.No.7 Exploration Institute of Geology and Mineral Resources, Shandong Linyi 276000, China; 2. Linyi Bureau of Land and Resources, Shandong Linyi 276000, China; 3. Land Science and Technology College of China University of Geosciences (Beijing), Beijing 100083, China)

Yishu fault zone is rich in geothermal resources. By using Landsat ETM+ as the main remote sensing data, and inversing the whole surface temperature in Linyi city by using radiative transfer equation. Gray temperature images have been interactively stretched in order to show temperature abnormal areas in 33℃ and the above temperature. Combining with ETM+ true color images and Google images, excluding town, rural residential buildings, river beach, bare rock and other abnormal areas, the already gained abnormal results have been compared with geothermal geologic background images, and analyze it from the angle of regional geology in order to find potential geothermal regions.

Landsat ETM+; temperature inversion; geothermal

2014-07-13；

2014-09-24；編輯：曹麗麗

劉孝才(1983—)，男，山東臨沂人，工程師，主要從事地質(zhì)礦產(chǎn)工作；E-mail:285410775@qq.com

P208

B