陳婉佳，高孝杰

（1. 成都理工大學 地球科學學院，四川 成都610059；

2 . 國土資源部 國土資源信息技術與應用重點實驗室,四川 成都 610059）

基于RS和GIS的龍門山南段活動構造信息提取

陳婉佳1,2，高孝杰1

（1. 成都理工大學 地球科學學院，四川 成都610059；

2 . 國土資源部 國土資源信息技術與應用重點實驗室,四川 成都 610059）

以龍門山南段為研究區(qū)，基于ETM和DEM數據，應用主成分分析顯示斷裂相關地層、水系信息差異,定向濾波分析斷裂位置及活動性。通過山體陰影分析、坡度分析，發(fā)現研究區(qū)域地勢及地表坡度變化強烈的地區(qū)是斷裂構造較發(fā)育或活動的區(qū)域，為活動構造的解譯提供依據。初步建立研究區(qū)活動構造解譯標志，對區(qū)內的線性構造進行解譯，分析研究區(qū)主要活動構造分布特征和地質地貌特征。

活動構造；ETM；DEM；龍門山南段

研究區(qū)域位于龍門山斷裂帶南段,地理坐標為北緯30°02'43"～30°38'47"，東經102°41'53"～103°35'23"。該區(qū)是龍門山推覆構造帶、康滇南北帶和北西向甘孜地槽褶皺帶交接地區(qū)，構造活動明顯。

1 數據處理與活動斷裂解譯標志

1.1 數據源及圖像預處理

數據源為Landsat-7/ETM和ASTER GDEM數據。本文應用ASTGTM_N29E102T、ASTGTM_N29E103Y兩幅DEM數據和L71130039_03920020405景遙感影像，覆蓋工作區(qū)。影像紋理清晰、色調柔和、層次感好、無掃描壞線、云量少、整體效果較好。首先通過幾何校正、741波段假彩色合成、圖像融合、圖像裁剪等處理得到研究區(qū)的遙感影像。

1.2 遙感圖像處理方法

斷裂活動在地質地貌上有明顯表現，通過對元數據進行增強處理，提取相關信息，從而達到識別斷裂的目的。

1.2.1 主成分分析

在遙感圖像中，各波段之間常常高度相關。經過主成分分析，除去波段之間的冗余信息，將多波段的圖像信息壓縮到少數幾個互不相關的波段，對斷裂相關信息進行增強。主成分分析后1、2、3波段合成圖上能夠清晰地看到受活動斷裂控制、淺色調的第4系地層呈串珠狀分布。水系彎曲的方向和第4系地層的長軸方向，都指示著斷裂的運動方向。

1.2.2 纓帽變換

纓帽變換是根據多光譜遙感中土壤、植被等信息在多光譜空間中的信息分布結構，對圖像所作的經驗性線性正交變換。由于不同巖性的含水性不同以及植被不一樣，經纓帽變化的1、2、3成分組合圖像上，突出了水系特征[1-3]。圖1為中里斷層纓帽變化后的圖像，可以明顯看到上盤受斷層控制而成梳狀的水系，以及下盤不規(guī)則的紋理特征。且在中里附近，斷層切割丘陵，丘陵兩側被錯切。

圖1 中里斷層纓帽變換后123波段組合圖

1.2.3 定向濾波

定向濾波是有選擇性地增強特定方向的圖像特征。通過設置，使得特定方向的線性影像和細微紋理得到突出，其他灰度變化信息被抑制，從而有利于線性構造的解譯[4]。

本文通過對ETM全色波段進3×3模板、 0°～90°之間每隔30°的定向濾波，得到如圖2的0°（正北）、N 30°E、N60°E定向濾波結果?？梢钥吹剑?N60°W 向，沿著斷裂的山脊錯斷以及沿線的溝谷、河流被明顯增強，明暗交替的線性紋理的扭動方向指示了斷裂的運動方向，而對 0°、30°突出效果相對不明顯，表明與斷裂展布方向一致的定向濾波可以突出相關的地形和解譯標志信息。

圖2 3個方向的定向濾波圖

1.3 DEM增強方法

活動構造與現代地貌密切相關，其發(fā)展在一定程度上塑造了目前的地形地貌格局。因此，利用DEM進行區(qū)域地形地貌分析，能夠突出活動構造方面的特征，提高識別能力，同時驗證了其他增強方法的正確性，提高了可信度[5]。

1.3.1 山體陰影分析

研究區(qū)北西部分屬于山區(qū)，在斷裂發(fā)育處往往有許多微地貌呈線狀排列的負地形，如深切溝谷、陡坎等，直接觀察遙感圖像則難以分辨。利用DEM數據提取不同方向的山體陰影圖，可增加山體的立體感，對于該方向垂直的負地形也有突出顯示作用，有助于相應斷裂的解譯[3]。

以北北東向斷裂為例，如圖3所示。較之圖b，a即當太陽方位角為50°（約垂直于斷裂的展布方向）時，斷裂的線狀影紋被尖銳突顯。在圖上可以十分清晰地看到線性排列的陡坎，表明太陽方位角與斷裂延伸方向垂直的山體陰影分析，能突出相關斷裂解譯標志信息。

圖3 不同太陽方位角山體陰影圖

1.3.2 坡度分析

基于ASTER GDEM數據生成坡度圖對研究區(qū)進行坡度分析。研究區(qū)坡度值為0°～85.300 9°，從紅到藍坡度值逐漸增大。結合遙感圖像，并對呈線狀分布的坡度突變帶及坡度較大地區(qū)的分析研究，再來提取斷裂信息。以五龍斷層為例，構成了明顯的坡度分界線。圖4中，北西側主體為藍色，為切割破碎帶的陡坡；南西側主要為黃色和紅色，為相對平坦的夷平面。

五龍斷裂構造活動強烈而且形成較早，斷裂活動控制的兩盤區(qū)域的地形差異反差大，表現在坡度上為一盤坡度較大，另一盤坡度相對較小，沿斷裂呈線狀延伸，2盤之間坡度變化率較高。通過受控斷裂地形坡度的突然變化以及線性延伸可以有效識別活動斷裂的位置走向，并對其他增強方法顯示的斷裂進行驗證。

圖4 五龍斷裂的坡度分析圖

1.4 斷裂解譯標志

活動斷裂的解譯標志是基于遙感圖像提取活動構造信息的依據。受斷裂控制的地質體、水體等地物，其反射來自太陽電磁波輻射的能力具有差異性，在遙感影像上表現為不同色調和不同形態(tài)特征[6]?；顒訑嗔呀庾g標志主要有斷層三角面、斷層陡坎、地層標志、水系展布特征變化、山脊溝谷的錯動。

2 研究區(qū)活動斷裂遙感信息特征

龍門山斷裂帶南段由北西向南東依次為五龍斷裂帶（F1）；小關子-下井溪斷裂帶（F2）；中林-雙石斷裂帶（F3）；蒙頂山斷層（F4）；白云斷層（F5）。這5條主干斷裂現在仍具有一定的活動性，有著地震活動、線性地貌、地形陡坎等表現，如圖5所示。

2.1 五龍斷裂帶

五龍斷裂的遙感影像圖上可以看到溝谷錯斷和山脊錯斷的現象。同時，其山體陰影圖上，白色為水系，紅色則是斷裂，整體上可以看到斷裂旁的水系發(fā)生扭曲，圖上明暗交替的紋理指示了斷裂的活動性。因此，初判斷五龍斷裂屬于右旋走滑，這條斷裂向北東延伸，與龍門山中段的茂汶斷裂或映秀斷裂相接，向南西終止在鮮水河斷裂帶上[7]。

2.2 小關子-下井溪斷裂帶

在主成分分析后，1、2、3波段合成圖中，可以清晰地看到水系的扭轉。水系由于受到斷裂的控制顯示出一致性彎曲，這指示出斷裂的活動。在斷裂的南東側，也可以看到線性異常現象，也顯示出了斷裂的活動性。

2.3 雙石-西嶺斷裂帶

該斷裂由3部分組成，并相互錯疊在平面上呈兩端散開、中間收斂的束狀，反映出以逆沖推覆為主兼具右行走滑的性質。它的遙感信息特征分為3段來描述:第1段主要位于斷裂西南段，從天全縣的白果樹經沙坪村向東北延伸至蘆山縣羅家壩與第2段匯合。該段的遙感影像圖上，可以清晰地看到線性陡坎，將斷裂帶痕跡突顯出來，且斷裂上、下盤地層特征有明顯差異。斷裂下盤紋理規(guī)則，而上盤顯示為不規(guī)則運動，初步判斷該段主要為活動性逆斷層。第2段主要位于斷裂中段，從蘆山縣的羅家壩向東北經雙石鎮(zhèn)一直延伸到圍塔地區(qū)。從山體陰影圖上，雙石鎮(zhèn)和寶勝鄉(xiāng)中間位置可以明顯看到斷裂錯移。在雙石-大川鎮(zhèn)部分，斷裂帶經過丘陵，表現出山脊錯斷現象。同時在該段上部可見明顯的陡坎，且兩側巖層存在顯著差異，斷層三角面和山脊錯斷等顯示出該段的走滑特征。第3段則位于斷裂東北段。由于構造作用，衍生出樹枝狀的次生斷層，發(fā)育于須家河組。在纓帽變化后1、2、3波段組合圖上，可以看到呈線狀排列的一系列陡坎。亮綠色為河流，受斷裂控制，河流呈肘狀變化。而殘存地層多表現為規(guī)模大小不一的短片或構造透鏡體，在西嶺鎮(zhèn)左側部分，大的斷片在平面上可長達5 km[8]。

2.4 蒙頂山斷層

在研究區(qū)內，斷層南起雅安隴西鄉(xiāng)的騎山廟，往北東延至名山縣官田壩，至觀音場附近，再往北東斷層被第4系掩蓋，在地表已不明顯，表現為陡立帶。在研究區(qū)內，斷層長約10 km，在地貌上表現為小陡坎或鞍狀，如圖6。

圖5 龍門山南段主要活動斷裂分布圖

圖6 蒙頂山斷層的330°山體陰影圖

由于構造運動性質的差異，在名山縣的蒙頂山地區(qū)造成河流襲奪及深切河曲現象，留下了清楚的地貌現象。如圖7所示。由于蒙頂山背斜及蒙頂山斷層受第四紀強烈褶皺作用和蒙頂山斷層的逆沖推覆作用，使蒙頂山地區(qū)快速隆升，造成蒙頂山西側的青衣江支流中里河襲奪了名山河[9]。

圖7 名山縣蒙頂山地區(qū)ETM影像圖(實線為現今河流；虛線為被襲奪河流)

2.5 白云斷層

斷層南起名山縣中峰寺，往東延伸至浦江縣白云鄉(xiāng)。斷層地表出露較好，在遙感圖像上呈現清楚的線性構造。在地貌上該斷層帶表現十分明顯，沿斷層的上升盤形成數km的斷層陡坎。斷層活動性表現為切錯朝陽鄉(xiāng)東側砂子硬山體，并在馬鞍石-余山一代形成線狀斷層陡崖，沿斷層多有崩塌、滑坡等地質災害活動。遙感影像顯示為向西斜沖，水平方向向右行錯動。第四紀活動性較強，曾有中小地震發(fā)生，通過研究該斷裂帶的地震地質條件，有助于闡明該斷裂帶發(fā)生中小地震的原因及強度。

3 結 語

遙感技術以其客觀性、真實性等特點，在龍門山南段活動構造信息提取過程中通過多種圖像處理方法，增強了活動斷裂特征，建立了解譯標志，是一種快速、經濟、實用的手段。地理信息系統(tǒng)技術以其強大的空間分析功能，實現了基于DEM數據的活動構造信息提取，如山體陰影分析、坡度分析等，成為斷裂活動性研究的依據，減少了野外勘察的工作量。

[1] Oguchi T, Aoki T, Matsuta N. Identification of an Active Fault in the Japanese Alps from DEM-based Hillshading[J].Computers & Geosciences， 2003(29): 885-891

[2] 付碧宏,張松林,謝小平,等. 阿爾金斷裂系西段——康西瓦斷裂的晚第四紀構造地貌特征研究[J]. 第四紀研究, 2006 (2): 228-235

[3] 陳文凱,張景發(fā),姜文亮,等. 基于TM和DEM的茅山地區(qū)斷裂構造解譯[J].地殼構造與地殼應力文集,2007(19):67-75

[4] 劉小利,李雪,李井岡,等. ETM和DEM數據在斷裂解譯中的應用[J].大地測量與地球動力學,2012,32(6): 50-53

[5] 袁小祥. 多源遙感數據在活動構造信息提取中的應用研究[D].北京:中國地震局地震預測研究所,2011

[6] 單新建,葉紅,陳國光,等. 數字遙感圖像處理的多源數據融合方法在地質學中的一些應用[J].地震地質,1999, 21(4): 465-472

[7] 張岳橋,李海龍. 龍門山斷裂帶西南段晚第四紀活動性調查與分析[J].第四紀研究, 2010(2): 699-710

[8] 林茂炳,陳運則. 龍門山南段雙石斷裂的特征及地質意義[J].成都:成都理工學院學報, 1996,23(2) : 64-68

[9] 李漢武. 四川名山——蒲江地區(qū)新構造特征及活動性分析[D].成都:成都理工大學,2011

P208

B

1672-4623（2015）01-0119-04

10.3969/j.issn.1672-4623.2015.01.039

陳婉佳,碩士，主要從事測繪、遙感地質與地理信息系統(tǒng)研究。

2014-01-03。

項目來源：教育部博士點基金資助項目（20095122110003）。