
Citation: | YANG Qing-hua, WU Xiao-juan, XIAO Zheng-hao, et al., 2015. TECHNOLOGY RESEARCH ON MINERAL EXTRACTION BY USING CASI/SASI AIRBORNE HYPERSPECTRAL DATA. Journal of Geomechanics, 21 (2): 241-251. |
2008年3月21日发生在新疆于田(东经81°32',北纬35°28',深度23 km) [1]的7.3级地震是继2001年11月昆仑山口西8.1级地震后中国内陆发生的最大一次7级以上地震,打破中国大陆7级以上地震平静多年的现象。地震是构造活动的集中反映,会伴随着较大的地表形变及断裂的明显活动,能够放大构造的正常活动方式,同震变形场是这种放大作用在地表的直观反映,对于认识发震断裂运动性质,研究邻近构造活动性具有重要意义。同震形变监测和发震构造活动性的研究方法主要有现场测量、构造形迹分析、震源机制解和数值模拟等。然而在青藏高原西北缘恶劣的自然环境中,常规方法无法对于田地震开展及时有效的同震形变测量,凸显出差分干涉雷达(InSAR,Differential Interferometry Synthetic Aperture Radar)的优势[2~4]。本文介绍了InSAR观测原理,并重点分析了通过干涉雷达对于田地震同震形变场的观测结果,据此分析发震机制及邻近构造带的活动性。
于田7.3级地震宏观震中位于青藏高原西北缘的阿什库勒地区(见图 1)。青藏高原北缘的2条大断裂自第四纪以来的运动方式与滑移速率对认识青藏高原西北部变形作用和机制,以及对地震的预测具有重要意义。研究区断裂活动的复杂性源于印度板块冲击欧亚板块,青藏高原开始对大陆岩石圈碰撞挤压,在塔里木盆地古老岩石的阻挡下,使青藏高原地壳缩短和高原隆升。由于青藏高原西部的块体向北推进,阿尔金断裂及其南盘山体隆起强烈,山前、山间拗陷剧烈,断错地貌明显,震区及附近的构造变形活动异常剧烈,尤以垂直运动最为显著,运动幅度达5000~7000 m[5]。地壳活动块体间的相互作用,使得该地区的新构造活动状态更为复杂,垂直、水平运动和内部变形均很明显,导致了高原边缘地区的构造以挤压、拉张为主,兼有走滑运动[6]。大的区域构造应力方向为南北向,所以形成的几个大的断裂带都是东西向的,如:阿尔金断裂带、康西瓦断裂带、大红柳滩—郭扎错断裂、普鲁断裂等。同时于田地震就发生在康西瓦断裂和阿尔金断裂的交汇处,所以通过震源机制解和其他的一些方法都认为于田地震所处的断裂走向应该是东西向,是阿尔金深断裂左旋扭错,致使断裂带西南端的破裂构造继续活动的结果。
干涉雷达作为一项新的地球物理技术,是以合成孔径雷达复数图像的相位信息获取地表变形信息的技术,测量精度可以达到厘米级甚至更高。InSAR技术的应用为地震地质研究开拓了一条新的途径。InSAR技术可以使我们直接获取大范围的、连续空间覆盖的断层位移和运动速率等定量化的基础数据,为深入了解断裂活动方式、区域地壳稳定性、地震重复间隔、潜在震源区等问题做数据支撑。InSAR技术不受大空间尺度和复杂自然因素影响,可以获得极为经济的垂直形变信息,结合GPS的高精度水平形变信息,在时间轴上进行矢量叠加,可以得到区域地表的三维动态位移场、速度场和应变场,对于监测震前的微量形变、震时的同震位错和震后的变形回弹等地形形变研究有着极其重要的意义,可为地震预测提供丰富可靠的定量信息。
地形变InSAR监测至少需要根据同一区域两次重复的雷达数据获取,雷达传感器的回波信号携带了地物后向散射体的相位和强度信息,计算同一区域不同时间获取的两景单视复数(SLC,Single Look Complex)雷达影像(φm及φs)的相位差生成干涉图Δφ,该干涉图中既包含了两次成像期间地表相对运动的相位信息(φdef),也含有成像区域的地形信息(φtopo)、观测向斜距信息(φflat)、传感器轨道误差(φorbit)、大气效应(φatmos)和其他去相干因素的相位(φnoise)值,公式表示为:
|
(1) |
差分干涉的基本任务就是从干涉图中提取有用的φdef信息。式(1)中的地形相位φtopo可以采用数字高程模型或多轨观测方法去除,观测向斜距φflat属于系统观测常量,通过卫星姿态参数校正,其他相位信息引起的测量误差在厘米级,是阻碍InSAR测量精度提高的主要因素,目前尚未有完善的解决手段。处理后的地形变干涉相位信息(φdef)与沿传感器视线向(LOS,Line of Sight)地表变形ΔR的关系为:
|
(2) |
式中λ为雷达波波长,cm。
可见对于相同的地表变形ΔR,雷达波长与干涉条纹的密度成反比,对于同一变形,波长分别为X波段(3.0 cm,TerraSAR-X & COSMO-SkyMed)、C波段(5.6 cm,EnviSAT ASAR & ERS1/2)和L波段(23.5 cm,ALOS卫星的PALSAR传感器)的干涉条纹如图 2所示。干涉图像中每个干涉条纹代表波长一半的变形量,即X波段图像上的一个条纹为1.5 cm变形,C波段为2.8 cm变形,L波段为11.8 cm变形[3~4],也就是波长越短相对测量精度越高。但短波雷达数据受环境影响显著,不易形成干涉,因此在应用中应结合自然环境和研究精度需求确定采用的波段。
本文中采用“3.21”于田地震前后的PALSAR雷达数据,结合数字高程模型(DEM)数据,对阿什库勒盆地及其邻近构造带进行“二通道+ DEM”的干涉监测(见图 3)。于田地震影响区主要为高原山区与现代冰川,雷达回波中噪声较大,ALOS卫星的PALSAR雷达数据拥有较长的L波段(23.6 cm),在2次成像中雷达后向散射保持很好的相关性,可以获得较高质量的全区干涉图(每个条纹周期代表 11.8 cm的变形)。
InSAR监测地震同震形变场,必须要有震前和震后观测的雷达数据。按ALOS卫星的观测策略,震前一年内库存了至少1次该地区的PALSAR观测数据。在“3.21”震后日本宇航局(JAXA)启动应急补充观测,构成了完整的阿什库勒地区地震前后的干涉像对,本文中应用了震前震后的3对雷达数据。数据的轨道号、获取时间、垂直基线长度等信息见表 1。
![]() |
据日本航空航天局(JAXA)的全球观测策略,目前在于田地区PALSAR雷达数据大部分为升轨时拍摄,垂向入射角均为34.3°,视线向为北东81.84°,相应的雷达视线向变形量(dLOS)与垂直(dup)、北(dn)、东(de)变形量关系为[3] :
|
(3) |
可以看出,垂直形变与东西向形变的敏感度相差不大,但两者对视线向的贡献为相反方向。如果野外调查获得了地表水平和垂直位错量,通过公式(3),再结合断层的走向和雷达成像几何关系,可以计算出实际的视线向位错量,从而实现野外观测与D-InSAR的互相检验和比较。图 4将PALSAR干涉雷达形变图像附在地形图上,是缠绕在一个相位周期(-π~π弧度)内的D-InSAR观测结果,相当于视线向变形的等值线,PALSAR数据波长为23.6 cm,每个条纹代表 11.8 cm的视线向形变。这其中还包括InSAR的地形误差、轨道误差、大气延迟误差和电离层干扰误差等,但相对于本次地震形变量居于次要地位。
通过图 4可见发震断裂走向总体为南北向,此地震破裂主要以正断层破裂方式为主,兼有左旋走滑运动的倾滑破裂特征:西向为上盘发生沉降,东向为下盘发生隆起。同时,根据图 4可以粗略地算出该发震断裂的长度为25.6 km,明显的地震形变影响范围为2500 km2,西侧上盘的视线向形变为200.6 cm (11.8 cm/条纹× 17条纹),按入射角可分解为上盘下沉量166.5 cm,上盘水平位移量112.3 cm。虽然位移量较大,但此次于田地震所发生的位置人员稀少,所造成的损失也比较小,没有人员伤亡。通过InSAR干涉测量方法,能准确地确定断裂走向和地震形变影响范围,为分析青藏高原西北缘的构造变形和动力学机制提供准确、客观的依据。
Shan等(2009) [7]用EnviSAT卫星的C波段(5.6 cm)雷达数据对于田地震也进行了干涉观测(见图 5)。通过比较可见C波段干涉图中,发震断裂的总体走向也是南北向,其东盘的干涉条纹比较清晰,而西盘的条纹模糊,这与本文中所做出的干涉条纹情况刚好相反。在研究地表变形的时候,可以综合图 4、图 5这2个干涉图进行分析。
干涉图产生差异的原因: ①由于卫星的类型不同,本文使用的是日本ALOS卫星的PALSAR L波雷达数据,而干涉图 5用的是欧空局EnviSAT卫星的ASAR-C波雷达数据; ②卫星的拍摄角度不同,本文采用的ALOS卫星PALSAR数据是升轨数据,雷达视线向为近东向(见图 4); EnviSAT卫星的ASAR是降轨观测,视线近西向(见图 5); ③数据拍摄的时间不同,这对形成干涉的效果也会产生一定的影响。
结合目前地震前后各方面的观测结果,可以从InSAR变形场中解译出重要的构造信息:由于青藏高原西南部的块体向北推进,阿尔金断裂及其南盘山体隆起强烈,山前、山间拗陷剧烈,断错地貌明显; 同时地壳活动块体间的相互作用,使得这个地区的新构造活动状态更为复杂,垂直、水平运动和内部变形均很明显。于田地震所处的断裂走向不是东西向,而是和大的作用力方向一致的南北向,这并不矛盾。在大的南北向作用力作用下,大断裂的走向基本上都是东西向的,特别是板块碰撞边界线更是如此; 但于田地震发生处与板块碰撞边界线还有一些距离,在东西向阿尔金断裂和康西瓦断裂的共同作用下,阿尔金断裂和康西瓦断裂的中间地带是有可能形成南北向破裂的,同时与发震断裂近平行发育的玉龙喀什河的走向也佐证了这条断裂的走向。通过震源机制解分析这次地震就是由于阿尔金断裂的左旋扭错,致使阿尔金断裂西南端发生破裂,只是破裂的方向是与阿尔金断裂垂直的南北向。这条南北向断裂也正反映了该地区构造活动的复杂性,说明断裂之前也存在很强的相互影响,本次于田7.4级地震之后的一个多月就发生了5月12日的四川汶川8级地震。
通过ALOS卫星的PLSAR雷达数据InSAR形变观测,本文获得了于田地震的同震地表形变图像,对于理解该地震的发震机理及该区域的构造变形具有重要意义。
2008年3月21日于田Ms 7.3级地震的宏观震中位于康西瓦断裂东南端的南北2分支(大红柳滩断裂和墓士山南麓断裂)和阿尔金断裂西南端帚状的3个分支交汇的三角地带。于田地震引发了阿尔金断裂的一个近南北向分支断裂的同震地表破裂,破裂长度为25.6 km。根据干涉雷达变形图像推测于田地震震源机制解为正断层,断层倾向西,上盘的最大运动幅度在200 cm以上。于田地震所在的两大断裂交汇处的构造应力场以近南北向挤压为主。
虽然距离得出精确的断层破裂模型还有较长的路要走,但是从原始的干涉形变图像上已经可以得出一些重要的地质和地球物理信息,这些信息对于深入研究地震构造都具有重要意义。下一步需要研究宽幅InSAR同震形变的观测,获取青藏高原北缘更大范围内对本次地震的变形反应,以期待全面分析同震和震后应力应变场演化等更为重要的科学问题。
[1] |
Clark R N, Swayze G A, King T V V. Imaging spectroscopy: Earth and planetary remote sensing with the USGS tetracorder and expert systems[J]. Journal of Geophysical Research, 2003, 1008(E12): 5131. http://cat.inist.fr/?aModele=afficheN&cpsidt=15473253
|
[2] |
唐攀科, 李永丽, 李国斌, 等.成像光谱遥感技术及其在地质中的应用[J].矿产与地质, 2006, (4):160~164. http://www.cnki.com.cn/Article/CJFDTOTAL-KCYD200602012.htm
TANG Pan-ke, LI Yong-li, LI Guo-bin, et al. Imaging spectrometry remote sensing technology and its applications in geology[J]. Mineral Resources and Geology, 2006, (4): 160~164. http://www.cnki.com.cn/Article/CJFDTOTAL-KCYD200602012.htm
|
[3] |
冯孝良, 管仕平, 牟传龙, 等.西藏甲马铜多金属矿床的岩浆热液交代成因:地质与地球化学证据[J].地质地球化学, 2001, 29(4):40~48. http://www.cnki.com.cn/Article/CJFDTOTAL-DZDQ200104006.htm
FENG Xiao-liang, GUAN Shi-ping, MU Chuan-long, et al. Geological characteristics and genesis of the Jiama copper polymetallic deposit in Tibet[J]. Geology Geochemistry, 2001, 29(4): 40~48. http://www.cnki.com.cn/Article/CJFDTOTAL-DZDQ200104006.htm
|
[4] |
崔廷伟, 马毅, 张杰.航空高光谱遥感的发展与应用[J].遥感技术与应用, 2003, 18(2):118~122. doi: 10.11873/j.issn.1004-0323.2003.2.118
CUI Ting-wei, MA Yi, ZHANG Jie. The development and applications of the airborne hyperspectral remote sensing[J]. Remote Sensing Technology and Application, 2003, 18(2): 118~122. doi: 10.11873/j.issn.1004-0323.2003.2.118
|
[5] |
甘甫平, 王润生, 马蔼乃.基于特征谱带的高光谱遥感矿物谱系识别[J].地学前缘, 2003, 10(2):445~454. http://www.cnki.com.cn/Article/CJFDTOTAL-DXQY200302033.htm
GAN Fu-ping, WANG Run-sheng, MA Ai-nai. Spectral Identification Tree (SIT) for mineral extraction based on spectral characteristics of minerals[J]. Earth Science Frontiers, 2003, 10(2): 445~454. http://www.cnki.com.cn/Article/CJFDTOTAL-DXQY200302033.htm
|
[6] |
燕守勋, 张兵, 赵永超, 等.高光谱遥感岩矿识别填图的技术流程与主要技术方法综述[J].遥感技术与应用, 2004, 19(1):52~63. http://www.cnki.com.cn/Article/CJFDTOTAL-YGJS200401012.htm
YAN Shou-xun, ZHANG Bing, ZHAO Yong-chao, et al. Summarizing the technical flow and main approaches for discrimination and mapping of rocks and minerals using hyperspectral remote sensing[J]. Remote Sensing Technology and Application, 2004, 19(1): 52~63. http://www.cnki.com.cn/Article/CJFDTOTAL-YGJS200401012.htm
|
[7] |
王润生.高光谱遥感的物质组分和物质成分反演的应用分析[J].地球信息科学, 2009, 11(3):261~267. http://www.cnki.com.cn/Article/CJFDTOTAL-DQXX200903002.htm
WANG Run-sheng. Spectral identification and inversion of composition and component of objects with hyperspectral remote sensing[J]. Geo-information Science, 2009, 11(3): 261~267. http://www.cnki.com.cn/Article/CJFDTOTAL-DQXX200903002.htm
|
[8] |
叶发旺, 刘德长, 赵英俊.CASI/SASI航空高光谱遥感测量系统及其在铀矿勘查中的初步应用[J].世界核地质科学, 2011, 28(4):231~236. http://www.cnki.com.cn/Article/CJFDTOTAL-GWYD201104010.htm
YE Fa-wang, LIU De-chang, ZHAO Ying-jun. Airborne hyper-spectral survey system CASI/SASI and its preliminary application in uranium exploration[J]. World Nuclear Geoscience, 2011, 28(4): 231~236. http://www.cnki.com.cn/Article/CJFDTOTAL-GWYD201104010.htm
|
[9] |
赵泽辉, 郭召杰, 王毅.甘肃北山柳园地区花岗岩类的年代学、地球化学特征及构造意义[J].岩石学报, 2007, 23(8):1847~1860. http://www.cnki.com.cn/Article/CJFDTOTAL-YSXB200708006.htm
ZHAO Ze-hui, GUO Zhao-jie, WANG Yi. Geochronology, geochemical characteristics and tectonic implications of the granitoids from Liuyuan area, Beishan, Gansu Province, northwest China[J]. Acta Petrologica Sinica, 2007, 23(8): 1847~1860. http://www.cnki.com.cn/Article/CJFDTOTAL-YSXB200708006.htm
|
[10] |
江思宏, 聂凤军.北山地区花岗岩类成因的Nd同位素制约[J].地质学报, 2006, 80(6):8261~842. http://www.cnki.com.cn/Article/CJFDTOTAL-DZXE200606005.htm
JIANG Si-hong, NIE Feng-jun. Nd-isotope constraints on origin of granitoids in Beishan Mountain area[J]. Acta Geologica Sinica, 2006, 80(6): 8261~842. http://www.cnki.com.cn/Article/CJFDTOTAL-DZXE200606005.htm
|
[11] |
左国朝, 刘义科, 刘春燕.甘新蒙北山地区构造格局及演化[J].甘肃地质学报, 2003, 1(12):1~15. http://www.cnki.com.cn/Article/CJFDTOTAL-GSDZ200301000.htm
ZUO Guo-chao, LIU Yi-ke, LIU Chun-yan. Framework and evolution of the tectonic structure in Beishan area across Gansu Province, Xinjiang Autonomous Region and Inner Mongolia Autonomous Region[J]. Acta Geologica Gansu, 2003, 1(12): 1~15. http://www.cnki.com.cn/Article/CJFDTOTAL-GSDZ200301000.htm
|
[12] |
赵泽辉, 郭召杰, 韩宝福, 等.新疆东部—甘肃北山地区二叠纪玄武岩对比研究及其构造意义[J].岩石学报, 2006, 22(5):1279~1293. http://www.cnki.com.cn/Article/CJFDTOTAL-YSXB200605019.htm
ZHAO Ze-hui, GUO Zhao-jie, HAN Bao-fu, et al. Comparative study on Permian basalts from eastern Xinjiang-Beishan area of Gansu Province and its tectonic implications[J]. Acta Petrologica Sinica, 2006, 22(5): 1279~1293. http://www.cnki.com.cn/Article/CJFDTOTAL-YSXB200605019.htm
|
[13] |
梅华林, 李惠民, 陆松年, 等.甘肃柳园地区花岗质岩石时代及成因[J].岩石矿物学杂志, 1999, 18(1):14~17. http://www.cnki.com.cn/Article/CJFDTOTAL-YSKW901.002.htm
MEI Hua-lin, LI Hui-min, LU Song-nian, et al. The age and origin of the Liuyuan granitoid, northwestern Gansu[J]. Acta Petrrologica Et Mineralogica, 1999, 18(1): 14~17. http://www.cnki.com.cn/Article/CJFDTOTAL-YSKW901.002.htm
|
[14] |
王润生, 甘甫平, 闫柏坤, 等.高光谱矿物填图技术及应用研究[J].国土资源遥感, 2010, (1):1~13. doi: 10.6046/gtzyyg.2010.01.01
WANG Run-sheng, GAN Fu-ping, YAN Bai-kun, et al. Hyperspectral mineral mapping and its application[J]. Remote Sensing for Land & Resources, 2010, (1): 1~13. doi: 10.6046/gtzyyg.2010.01.01
|
[15] |
Chaudhry F, Wu C C, Liu W M, et al. Pixel purity index-based algorithms for end member extraction from hyperspectral imagery[J]. Transworld Research Network37/661 (2), FortP.O., Trivandrum-695 023, Kerala, India.2006.
|
[16] |
Bateson A, Curtiss B. A method for manual endmember selection and spectral unmixing[J]. Remote Sensing of Environment, 1996, 55: 229~243. doi: 10.1016/S0034-4257(95)00177-8
|
[17] |
Craig M D. Minimum-volume transforms for remotely sensed data[J]. IEEE Transaction on Geoscience and Remote Sensing, 1994, 32(3): 542~552. doi: 10.1109/36.297973
|
[18] |
Nascimento J M P, Dias J M B. Vertex component analysis: A fast algorithm to unmix hyperspectral data[J]. IEEE Transactions on Geoscience and Remote Sensing, 2005, 43: 898~910. doi: 10.1109/TGRS.2005.844293
|
[19] |
Winter M E. Comparison of approaches for determining end-members in hyperspectral data[C]//IEEE Aerospace Conference Proceedings. IEEE, 2000: 305~313.
|
[20] |
Bowles J H, Palmadesso P J, Antoniades J A, et al. Use of filter vectors in hyperspectral data analysis[C]//Strojnik M, Andresen B F. Infrared Spaceborne Remote Sensing: Proceedings of SPIE. San Diego. SPIE, 1995: 148~157.
|
[21] |
Berman M, Kiiveri H, Lagerstrom R, et al. ICE: A statistical approach to identifying endmembers in hyperspectral images[J]. IEEE Transactions on Geoscience and Remote Sensing, 2004, 42: 2085~2095. doi: 10.1109/TGRS.2004.835299
|
[22] |
吴波, 汪小钦, 张良培.端元光谱自动提取的总体最小二乘迭代分解[J].武汉大学学报:信息科学版, 2008, 33(5):457~460. http://www.cnki.com.cn/Article/CJFDTOTAL-WHCH200805006.htm
WU Bo, WANG Xiao-qin, ZHANG Liang-pei. Iterative abstraction of endmember based on Total Least Square for Mixture Pixel Decomposition[J]. Geomatics and Information Science of Wuhan University, 2008, 33(5): 457~460. http://www.cnki.com.cn/Article/CJFDTOTAL-WHCH200805006.htm
|
[23] |
Plaza A, Martinez P, Perez R, et al. Spatial/spectral endmember extraction by multidimensional morphological operations[J]. IEEE Transactions on Geoscience and Remote Sensing, 2002, 40: 2025~2041. doi: 10.1109/TGRS.2002.802494
|
[24] |
Roggea D M, Rivard B, Zhang J, et al. Integration of spatial-spectral information for the improved extraction of endmembers[J]. Remote Sensing of Environment, 2007, 110: 287~303. doi: 10.1016/j.rse.2007.02.019
|