祁連山和首都圈衛星熱紅外背景場變化特征

溫少妍，屈春燕，單新建，閆麗莉，宋冬梅

(1.新疆維吾爾自治區地震局，烏魯木齊 830011;2.中國地震局地質研究所地震動力學國家重點實驗室，北京 100029;3.中國石油大學(華東)，青島 266555)

0 引言

衛星熱紅外遙感技術具有動態性強、覆蓋面積廣、連續穩定及可獲得大面積的連續場信息等優勢。很多學者就利用該技術對地震監測進行研究，并取得了較好的結果[1－3]。研究表明，震前熱紅外增溫出現的時間長短不一，從震前幾天到震前3個月都有存在，一般增溫幅度可達幾度[4]。震前這種增溫異常出現的時間和變化幅度特點，為利用多年背景場數據判別當日觀測圖像中的異常區域提供了基礎。然而，這些研究大多是圍繞震例展開的，專門針對衛星熱紅外亮溫背景場的研究較少。

衛星熱紅外圖像受地表各種環境因素的影響，正常動態亮溫背景場的建立是有效提取震前熱紅外亮溫異常的關鍵。國內學者已認識到背景場的建立直接影響到熱紅外異常的識別與提取，目前研究重心正在從異常的識別轉向背景場的建立[5－8]。陳順云等[5]借用分離窗算法和小波分析方法提取了中國地表亮溫的年變基準場，為進一步提取與分析異常場奠定基礎;屈春燕等[6]在分析地震活動熱紅外異常影響因素中指出建立紅外亮溫變化的基準模型及各類非震因素的扣除模型，將有可能較好地排除非震因素的影響，可有效檢測出與地震活動有關的異常;郭衛英等[7]定量分析了在地震活動相對平靜的年份內，活動斷裂造成了地形地貌差異與紅外輻射的動態變化，并認為地形地貌因素對斷裂帶的紅外亮溫也產生一定作用;溫少妍等[8]給出了區域背景場建立的方法，得出地形地貌及構造活動具有很好的相關性。許多地震的熱紅外圖像具有復雜性，有些是有增溫異常而無地震，有些則是有地震而無增溫異常。到目前為止，還不能完全解釋這種現象，在機理研究方面還需深入探索。

本文采用NOAA/AVHRR極軌氣象衛星數據，選定祁連山和首都圈2個地震重點監測區作為研究區，建立不同時間尺度的紅外亮溫背景場，并分析正常情況下背景場的時空演化規律及其與地形地貌的關系，為有效識別與地震有關的熱紅外異常信息提供方法依據。

1 數據源與研究方法

1.1 數據源及其預處理

為了降低太陽輻射對地表亮溫的影響以便更真實地反映地表熱輻射狀態的變化，本文選取NOAA衛星第5通道夜間2:00—5:00時的熱紅外遙感圖像。在亮溫背景場計算之前，對亮溫圖像進行幾何精糾正(誤差在1個像元以內)，并去除云和噪聲，從而確保計算出的亮溫背景場能夠真實反映無震時的亮溫正常動態變化[8]。

1.2 背景場建立方法

正常動態亮溫背景場以構造活動區多年相同時相觀測的均值和標準差作為表征量，計算時對多年累積的數據按地理坐標進行時間序列上的統計分析，即

式中:N 為可用數據年限;Ti(x，y，t)為位置(x，y)處第 i年 t時刻的亮溫觀測值;T(x，y，t)和 σ(x，y，t)分別為位置(x，y)處t時刻N年內，與研究區地震活動不相關的自然噪聲源影響下的亮溫平均值和標準差。

2 祁連山地區紅外亮溫背景場

2.1 研究區概況

祁連山研究區位于 E 98°～108°，N 32°～40°之間(圖1)，地處青藏高原東北部邊緣，西起祁連山北緣斷裂，東至六盤山斷裂，北接阿拉善地塊，南達東昆侖斷裂一線，包括祁連山脈的主體及河西走廊新生代陸盆地群。該區受到印度板塊和歐亞板塊的強烈碰撞和擠壓，形成祁連山北緣左旋逆走滑沖斷裂帶，主要斷裂帶有祁連山北緣斷裂、東昆侖斷裂、海原斷裂和西秦嶺北緣斷裂，多為NWW走向。區內河西走廊地處甘肅西北部地區，位于祁連山和龍首山之間，是祁連山北緣的山前拗陷，地震活動具有高強度和高頻度的特點[9]。

圖1 祁連山研究區范圍及地質構造背景圖Fig.1 Region and tectonic background of Qilian Mountains

2.2 背景場時空特征

利用自主研發的衛星紅外亮溫地震信息處理軟件對祁連山研究區2003—2011年的NOAA圖像進行不同時間尺度的紅外亮溫背景場計算，分析其時空演化特征。圖2可反映出紅外亮溫的冬春夏秋季節性變化特征。

圖2 祁連山2010年4個月份的亮溫背景場Fig.2 Brightness temperature background field of Qilian Mountains in four months of 2010

從圖2看出:研究區4個月份的亮溫背景場的變化主要受季節和地形的控制;其次，受活動斷裂帶的影響也較明顯。對2010年背景場亮溫均值進行統計，發現月份亮溫值呈現夏高冬低的年變特征，1—7月的亮溫均值逐漸升高，7月達到最高值11℃，8月亮溫開始逐漸下降，12月的亮溫值約－17℃。由于平均海拔較高，呈現整體平均亮溫偏低。在空間上，整個研究區以祁連山北緣斷裂為界分成2個不同的溫度區域，斷裂以北的河西走廊地帶，為相對高溫區，斷裂以南的青藏高原為相對低溫區。

圖3 祁連山DEM剖面與2010年7月的亮溫剖面對比Fig.3 Contrast between DEM profile and brightness temperature profile in July 2010

從圖3可以看出，亮溫變化與高程變化呈顯著負相關關系。在紅外亮溫剖面線處，地面高程最大值4 825 m時對應的亮溫值約2℃，地面高程最小值1 294 m時對應的亮溫值約21℃。說明該區地面高程每增加100 m，其亮溫降低幅度約0.54℃。這與我國氣溫直減率0.63℃/100 m基本一致。

活動構造斷裂帶在紅外圖像上有明顯反映，表現為亮溫高于周圍背景的線性條帶及不同溫度區域的分界帶。斷裂帶處的亮溫變化也符合季節規律的變化，總體表現出夏高冬低的特征。此外，河流和湖泊(如青海湖)在紅外圖像上清晰可見，往往是紅外圖像上較周邊地物為相對高亮溫的標志性地物，而且由于水體的熱容量和熱慣量比較大，其溫度變化不如陸表的劇烈。因此，通常也可以將這些水體的溫度作為斷裂構造活動變化的參照。

圖4所示為多年月的亮溫背景場圖。多年月的亮溫背景場的變化隨季節的變化表現出夏高冬低的特征，同時亮溫變化也受到地形地貌和斷裂構造的控制。對比單年月的背景場不難發現，在多年月的背景場圖像上斷裂構造的線性影像特征和高亮溫條帶不如單年月的背景場圖像上的構造特征清晰，且亮溫值明顯變小。原因在于參與多年月的背景場統計的圖像越多，亮溫背景場圖像上的一些突變信息就會被平滑。參與計算的年份越多，所得到的亮溫背景值則越可靠，就越能真實反映區域性亮溫變化特征。這也表明，可以認定多年月背景場得出的亮溫變化基準值可以為局部目標亮溫變化檢測提供參考，為斷裂活動及地震所引起的增溫異常識別提供穩定的動態亮溫背景場。

圖4 2003—2011年多年月的祁連山亮溫背景場Fig.4 Multi－ year month brightness temperature background field of Qilian Mountains from 2003 to 2011

3 首都圈地區紅外亮溫背景場

3.1 研究區概況

首都圈研究區地理范圍為 E 112°～118°，N 38°～42°(圖5)，位于華北拗陷盆地的北緣，包括山西斷陷盆地及華北平原的大部分。

圖5 首都圈研究區范圍及地質構造背景圖Fig.5 Region and tectonic background of Capital Zone

山西斷陷盆地斷裂系自北向南有延懷盆地、大同盆地、忻定盆地和太原盆地等主要斷陷盆地，呈右行斜列，位于鄂爾多斯塊體東側，NNE走向。控制盆地邊界的斷裂以NNE，NE和NEE走向為主，它們在第四紀的活動不僅控制了盆地的發育，還控制了盆地內次級構造單元的形成。區內一系列受伸展斷層控制的掀斜塊體在區域應力場的作用下，形成了典型的盆－嶺式構造[10]，一系列NNE走向的盆地和斷裂分布與強震活動密切相關，形成了以燕山山脈隆起為背景、面對華北平原廣大沉降帶強烈活動的格局[11]。

3.2 背景場時空特征

選取1999年和2003—2011年共10 a的NOAA/AVHRR衛星熱紅外夜間圖像數據，計算首都圈地區不同時間尺度的紅外亮溫背景場(圖6)，可以看出首都圈亮溫背景場的變化主要受季節和地形因素的控制，其次受斷裂構造的影響。在時間上，亮溫變化符合夏高冬低的季節變化特征;在空間上，整體表現為東南部華北平原的亮溫高于西北部山區，而且西北部的亮溫隨地面高程的增加而逐漸降低，說明亮溫與地面高程呈顯著負相關關系。

圖6 1999年和2003—2011年多年月的首都圈亮溫背景場Fig.6 Multi－year mouth brightness temperature background field of Capital Zone of 1999 and from 2003 to 2011

紅外亮溫剖面(2010年7月)與DEM高程剖面對比曲線(圖7)反映出亮溫與地形有明顯的負相關關系。圖7橢圓a所處位置為北京市區，橢圓b所處位置為天津市區，在這2處的亮溫形成2個峰值，這是由城市熱島效應所致。在紅外亮溫剖面線上高程最大值1 950 m所對應的亮溫值約9℃;在不考慮城市熱島效應的情況下，高程最小值0 m所對應的亮溫值約20℃。說明該區地面高程每增加100 m，其亮溫降低約0.56℃。東南部華北平原夏季亮溫最高溫約23℃，冬季亮溫最低溫約－13℃;西北部山區夏季最高溫約13℃，冬季最低溫約－20℃;平原地區的亮溫季節差大于西北部山區的，分析表明在一定區域和季節內，地面高程的變化對氣溫有較大影響。

圖7 首都圈DEM剖面與2010年7月的亮溫剖面對比Fig.7 Contrast between DEM profile and brightness temperature profile in July 2010

活動構造斷裂帶在紅外圖像上有明顯表現，大都表現為亮溫高于周圍背景，其走向呈與斷裂帶一致的線性條帶。斷裂帶內的亮溫變化也符合季節變化規律，總體表現為夏高冬低的特征。渤海灣海域亮溫隨季節變化的幅度小于東南部華北平原，而西北部山區亮溫隨季節變化的幅度相對較大，這主要與地表地物的熱輻射特性及地理緯度有關，而且水體具有較大的比熱容和熱慣量，受太陽輻射的影響相對較小。1999年和2003—2011年多年月4個月份的亮溫背景場反映了一定時期內亮溫穩定變化的狀態，可以為局部目標亮溫變化檢測提供參照背景，同時可作為震前紅外亮溫異常識別和檢測的基準。

4 對比分析

祁連山和首都圈研究區所處的緯度大體一致，但地形地貌及地理環境相差很大。前者位于青藏高原東北緣，地勢復雜，高程起伏大，植被稀少;后者位于我國華北地塊，地勢相對平坦，植被覆蓋好。

通過以上分析看出，2個研究區的紅外亮溫時空演化規律基本一致，但也有差異性。相同的時空特征主要表現為:空間上，亮溫背景場的變化主要受地形地貌和斷裂帶的控制;時間上，主要受季節變化規律的影響，亮溫最高值出現在7月，亮溫最低值出現在1月;夏季溫差總體趨勢大于冬季溫差。

由于地理環境的不同，2個研究區的亮溫值存在一定差別(圖8)。

圖8 2009年和2010年月份亮溫背景場變化Fig.8 Change of month brightness temperature background field in 2009 and 2010

從圖8可以看出:月份背景場總體都顯示出首都圈研究區的亮溫值高于祁連山地區;2個地區的溫差變化不穩定，在(－2℃，12℃)范圍內，每年具體的溫差各不相同，規律性不強，原因在于氣象復雜多變，反映到月份尺度上變化更是反復無常(如在2009年11月祁連山的亮溫值反而比首都圈地區高2℃，這可能是由于氣候變化，引起了亮溫變化)。因此，用單年份資料做背景來判斷異常，可能難以排出氣候等因素的影響;而以多年份的平均亮溫場做背景才能更加穩定地反映氣候特征影響亮溫變化的規律。

2個研究區亮溫背景場的變化符合季節規律變化，亮溫變化相對單年份的穩定(圖9)。首都圈的亮溫值明顯高于祁連山地區，二者間夏半年溫差(最大溫差9℃)大于冬半年溫差(最大溫差3℃)，這是由于祁連山地區的平均高程(2 536 m)大于首都圈地區的平均高程(902 m)。結果表明，亮溫與地面高程呈負相關關系，高程每增加100 m，亮溫就降低約0.21～0.63℃。這與方精云[12]所指出的我國氣溫直減率在0.25～0.63℃基本一致。另外，首都圈亮溫的變化較祁連山穩定，祁連山冬半年出現小幅度浮動，這可能是由于祁連山地區地形復雜引起的。

圖9 2003—2011年多年月平均亮溫背景場對比Fig.9 Annual variation contrast of multi－year month brightness temperature background field

5 結論

衛星熱紅外遙感技術為研究地震及斷裂活動提供了一種有效的觀測手段，在基于紅外資料識別提取地震異常的研究中，多年尺度紅外背景場的建立是異常判別的前提。本文基于NOAA衛星圖像數據建立了2個研究區的紅外亮溫背景場，同時分析亮溫變化與地形地貌、活動構造及季節變化的關系，為后期地震紅外異常判別提供依據，具有較為重要的研究價值。通過對2個不同區域的紅外亮溫背景場時空演化特征的分析，得到以下結論:

1)紅外亮溫背景場受季節變化影響最大，年變化規律明顯。夏半年亮溫逐漸升高，亮溫最高值出現在7月份;冬半年亮溫逐漸降低，亮溫在1月達到最低。這為排除氣候變化因素的影響識別地震異常提供了參考依據。

2)紅外亮溫受地形因素較大，與地面高程呈顯著的負相關關系。地面高程每增加100 m，亮溫降低約0.21～0.63℃，這與我國氣溫直減率在0.25～0.63℃基本一致。對于地形地貌復雜的地區，亮溫變化曲線呈現不穩定特征。

3)活動斷裂帶在亮溫圖像上具有明顯界線特征，呈高亮溫線性條帶，或亮溫分界帶。河流和湖泊等水系在紅外亮溫圖像上清晰可見，表現較周圍地物明顯的高亮溫標志。

4)多年平均亮溫背景場具有較好的穩定性。多年平均背景場平滑了氣候等突變信息，呈現出穩定性較強的規律性變化特征。多年平均亮溫背景場可作為亮溫變化的基準場，為斷裂活動和地震所引起的增溫異常檢測提供穩定的動態參考。

[1]Dimitar O，Nevin B，Thomas L，et al.Satellite thermal IR phenomena associated with some of the major earthquakes in 1999－2003[J].Physics and Chemistry of the Earth，2006，31(4/9):154－163.

[2]Saraf A K，Rawat V，Banerjee P，et al.Satellite detection of earthquake thermal infrared precursors in Iran[J].Nat Hazards，2008，47(1):119－135.

[3]馬 瑾，單新建.利用遙感技術研究斷層現今活動的探索——以瑪尼地震前后斷層相互作用為例[J].地震地質，2000，22(3):210－215.Ma J，Shan X J.An attempt to study fault activity using remote sensing technology:A case of the Mani earthquake[J].Seismology and Geology，2000，22(3):210－215.

[4]荊 鳳，申旭輝，張鐵寶，等.與地震有關的活動斷裂帶紅外輻射變化特征[J].國土資源遙感，2013，25(1):56－60.Jing F，Shen X H，Zhang T B，et al.Variation characteristics in infrared radiation of active fault zone related to earthquakes[J].Remote Sensing for Land and Resources，2013，25(1):56－60.

[5]陳順云，馬 瑾，劉培洵，等.中國地表亮度溫度年變基準場[J].地震地質，2004，26(3):528－538.Chen S Y，Ma J，Liu P X，et al.Normal annual variation field of land surface brightness temperature in China[J].Seismology and Geology，2004，26(3):528－538.

[6]屈春燕，單新建，馬 瑾.地震活動熱紅外異常的影響因素分析[J].地震研究，2007，30(2):113－119.Qu C Y，Shan X J，Ma J.Influence factors to seismic infrared anomaly[J].Journal of Seismological Research，2007，30(2):113－119.

[7]郭衛英，馬 瑾，單新建.活動斷裂帶的地形地貌差異與紅外亮溫年變特征的研究[J].地球物理學進展，2008，23(5):1437－1443.Guo W Y，Ma J，Shan X J.Researches on the yearly characteristic of IR brightness temperature and the different hypsography on the active faults[J].Progress in Geophysics，2008，23(5):1437－1443.

[8]溫少妍，屈春燕，閆麗莉，等.區域紅外亮溫背景場建立方法及初步應用[J].地震，2011，31(2):59－67.Wen S Y，Qu C Y，Yan L L，et al.Method for constructing regional brightness temperature background field and its preliminary application[J].Earthquake，2011，31(2):59－67.

[9]董治平，馮建剛，王 先，等.祁連山構造帶地震遷移活動特征分析[J].國際地震動態，2007，346(10):8－14.Dong Z P，Feng J G，Wang X，et al.Characteristics of earthquakes migration along the Qilianshan structural zone[J].Recent Developments in World Seismology，2007，346(10):8－14.

[10]韓月萍，楊國華，董運宏，等.首都圈地區現今地殼垂直形變特征[J].華北地震科學，2009，27(1):13－16.Han Y P，Yang G H，Dong Y H，et al.The current characteristics of vertical crustal deformation of the Capital Circal[J].North China Earthquake Sciences，2009，27(1):13－16.

[11]王若柏，顧國華，徐 杰，等.張家口－渤海地震構造帶的地殼形變研究[J].地震地質，2004，12(4):586－595.Wang R B，Gu G H，Xu J，et al.Discussion on characteristics of crustal deformation along the Zhangjiakou－Bohai seismotectonic zone[J].Seismology and Geology，2004，12(4):586－ 595.

[12]方精云.我國氣溫直減率分布規律的研究[J].科學通報，1992(9):817－820.Fang J Y.Study on distribution rules of air temperature lapse－rate in China[J].Chinese Science Bulletin，1992(9):817－ 820.