基于ASTER和WorldView-3數據在新疆坡北地區蝕變礦物信息提取對比研究

王平平　王婷　趙慧　姚虎　尹藝霖

摘? ?要：利用ASTER和WorldView-3數據進行預處理，通過借鑒和分析ASTER數據礦物指數、礦物指數法理論、蝕變礦物光譜特征規律，總結提出WorldView-3數據蝕變礦物指數。通過觀察研究區蝕變礦物信息分布情況，對比分析ASTER和WorldView-3數據蝕變礦物信息提取的適用性。ASTER數據更適用于赤鐵礦、褐鐵礦，黃鉀鐵礬礦物信息提取，WorldView-3數據更適用于白云母、綠泥石、綠簾石、菱鐵礦、方解石礦物的信息提取。同時，結合研究區地質特征，預測了成礦遠景區?；贏STER和WorldView-3數據蝕變礦物指數，實現兩種數據的優勢互補，完善了WorldView-3數據蝕變礦物信息提取。預測成礦遠景區對進一步找礦提供了參考。

關鍵詞：ASTER；WorldView-3；礦物光譜指數；蝕變礦物信息；成礦遠景

遙感影像儲存了大量地物信息，與找礦相關的礦化作用所產生的圍巖蝕變存儲在影像上的信息，可通過數字圖像處理和技術增強提取遙感影像上的成礦信息，從而快速、有效圈定成礦遠景區，為遙感地質找礦指明方向[1]。目前，較常使用的蝕變礦物信息提取方法有比值法、主成分分析法、光譜角法、混合匹配濾波法、礦物光譜指數法（比值法）等。

在遙感蝕變信息提取方面，ASTER數據應用較多，并取得了較好的效果。利用ASTER數據可見光近紅外和短波紅外波段提出了羥基、高嶺石、明礬石、方解石等4個礦物光譜指數[2]。使用ASTER數據采用光譜分析法提取與金礦化有關的蝕變礦物[3]。利用ASTER與Sentinel-2A融合數據提取云南普朗銅礦化蝕變信息[4]。使用ASTER數據采用光譜角分類法提取大平掌銅礦區外圍遙感礦化蝕變信息[5]。使用ASTER數據采用主成分分析法提取新疆西昆侖塔什庫爾干地區老并鐵礦區遙感蝕變異常信息[6]。利用ASTER可見光-近紅外波段和短波紅外波段，對校正后遙感數據采用波段比值分析、波段組合分析和主成分分析來提取粘土礦物、綠泥石和方解石化蝕變。結合野外勘查結果發現，主成分分析、波段組合分析能初步劃分礦區綠泥石化、碳酸鹽化和粘土化蝕變帶及硅化帶[7]。本文使用ASTER數據，對穆龍套金礦床及外圍地區開展遙感蝕變異常提?。豪貌ǘ伪戎捣?、礦物指數法提取硅化蝕變，利用主成分分析法提取MgOH、鐵染和碳酸鹽化蝕變信息?；贏rcGIS平臺，將提取出的遙感異常信息、搜集到的已知礦點及重砂異常疊加分析。結果表明，遙感異常密集區與已知礦點吻合度較高，這對研究及尋找同類型礦床有一定參考價值[8]。以雅魯藏布江縫合帶念扎金礦為研究區，通過對與Fe3+，Fe2+，ALOH，MgOH，硅酸鹽、碳酸鹽等成分相關特征的礦物光譜特征分析，運用“波段比值法＋假彩色合成法”進行細致的巖性單元分類與地質填圖，研究結果表明，基于ASTER光譜特征的巖性填圖和蝕變信息提取方法進行找礦預測，在西藏高海拔地區具天然的適用性[9]?；贠LI數據和ASTER數據，利用主成分分析法（PCA）對都蘭縣通突爾鉛鋅多金屬礦區進行羥基蝕變信息提取取得了良好效果，已知礦點與提取的羥基蝕變異常吻合程度較高[10]。利用主成分分析法對Landsat8和ASTER數據分別提取黃鐵礦、褐鐵礦等鐵染蝕變礦物和絹云母等羥基蝕變礦物，通過對Aster和Landsat8兩種數據提取的蝕變信息進行空間疊加分析，定量分析對比兩種數據的蝕變信息提取效果，具提高提取蝕變信息準確性的優勢[11]。閆柏琨和王潤生編制了ASTER遙感數據處理與礦物填圖技術指南，建立了鐵氧化物、硅酸鹽鐵、黃鉀鐵礬、高嶺石、白云母/蒙脫石、明礬石、碳酸鹽/綠泥石/綠簾石、SiO2含量、碳酸鹽、硫酸鹽礦物信息識別光譜指數[12]，為ASTER數據提取礦物信息提供指引。

WorldView-3數據國內外研究較少。2014年，數字地球公司發射了WorldView-3衛星，為目前市場上空間分辨率最高的多光譜商業遙感衛星，具較高的光譜分辨率。據WorldView-3數據參數設置，從AVIRIS數據中模擬出WorldView-3短波紅外數據，進行了蝕變礦物信息提取試驗[13]。使用混合調諧匹配濾波（MTMF）方法，利用在軌、模擬的WorldView-3短波紅外數據和AVIRIS短波紅外數據進行礦物填圖，發現在軌WorldView-3數據與AVIRIS數據提取結果相當，優于WorldView-3模擬數據[14]。在新疆坡北地區使用WorldView-3數據進行蝕變信息提取方法研究，發現WorldView-3數據豐富的短波紅外波段，非常適合識別含羥基蝕變，對于區分含ALOH蝕變和含MgOH蝕變有較好的效果，并提出了[15]?；鵚orldView-3數據，建立了羥基指數、Al-OH指數、Mg-OH指數、Fe-OH指數、鐵染指數、碳酸鹽指數、高嶺石指數、方解石指數、菱鐵礦指數、黃鉀鐵礬指數等10個礦物指數模型[16]。在內蒙古衛境地區鈾礦勘查中使用WorldView-3數據，采用主成分分析法提取蝕變礦物信息，對蝕變礦物信息和基礎地質資料進行綜合分析，圈出一處鈾成礦有利地段，并對部分區域開展野外查證，取得了較好的找礦效果，為該區后續鈾礦勘查提供了參考[17]?；赪orldView-3（16）波段遙感數據，利用主成分分析方法提取了鐵染、ALOH、MgOH和碳酸鹽異常信息。進一步查明巴音戈壁盆地本巴圖地區蝕變類型及分布，為該地區鈾礦勘查提供線索[18]?？偟膩碚f，WorldView-3數據在蝕變礦物信息提取方面研究較少，建立WorldView-3數據礦物指數的方法有限，因此加強對WorldView-3數據這方面的研究顯得至關重要。

1? 研究區地質特征

研究區位于新疆維吾爾自治區巴音郭楞蒙古自治州若羌縣，面積約162 km2。范圍為40°25′53″N～40°32′50″N、91°31′44″E～91°46′37″E。地表徑流及植被不發育，巖石出露較好，是開展遙感地質調查工作的有利地區[19]。礦產資源主要有金、銅、錳等。大地構造位于塔里木盆地北緣，北山裂谷西南緣及天山成礦帶南緣之間（圖1）。地質構造復雜，斷裂、褶皺構造較為發育。區內斷裂以NEE向為主，次為NW、NE向兩組斷裂，斷裂在地表無明顯活動痕跡，多為隱伏小斷層。區內發育的NE向斷裂對礦體起到破壞作用，且控制區內金礦分布。區內地層出露單一，主要為第四系、中石炭統茅頭山組上亞組和下石炭統紅柳園組。據紅柳園組的巖石類型，紅柳園組可以分為2個亞組，分別是紅柳園組下亞組和上亞組。其中，紅柳園組下亞組主要由碳質綠泥石片巖、玄武巖、泥硅質和鈣質板巖、石英片巖、銅礦體、金礦體組成。紅柳園組上亞組主要由碳質板巖、長石巖屑砂巖和蝕變凝灰巖組成[20-21]。據前人研究可知，紅柳園組下亞組為主要的含礦地層，大青山金礦產于該組地層中[22]。巖漿巖主要為花崗巖、閃長巖及酸性巖脈（圖3-A）。

該區域圍巖蝕變極發育，礦體主要產于兩種礦化蝕變帶中：①賦存于石英脈中的金礦（化）體。礦區內發育的石英脈與區域構造線一致，當石英脈周圍巖石黃鐵礦化、黃鉀鐵礬化、碳酸鹽化及絹云母蝕變較強時，金元素富集，形成金礦（化）體；②賦存于NE向斷裂構造碎裂巖破碎帶中的銅金礦體。發育較強的絹云母化、褐鐵礦化、黃鉀鐵釩化、綠泥石化等蝕變，其中絹云母、褐鐵礦化與銅金關系密切[21]。黃鐵礦在地表易被氧化為黃鉀鐵礬及褐鐵礦，結合前人研究可知，區內主要蝕變礦物為褐鐵礦、赤鐵礦、黃鉀鐵礬、云母、綠泥石、綠簾石、菱鐵礦、方解石等[23]。

2? 數據源與預處理

2.1? 數據源特征

本次使用ASTER數據包含3個VNIR波段，6個SWIR波段和5個TIR波段，波長范圍分別為0.52～0.86 μm，1.6～2.43 μm和8.125～11.65 μm，空間分辨率分別為15 m，30 m和90 m[24-26]。WorldView-3數據包括8個VNIR和8個SWIR波段，波長范圍分別為0.4～1.04 μm和1.21～2.365 μm，空間分辨率分別為1.24 m和3.7 m[27]。由于許可限制，WorldView-3數據可見光近紅外波段和短波紅外波段的空間分辨率為2 m和7.5 m（表1）。

2.2? 數據預處理

本次使用的ASTER遙感數據級別為L1T級，在數據生產時，已進行幾何校正，消除了由傳感器影響而造成的輻射畸變，因此，輻射校正只需消除大氣對輻射畸變的影響。ASTER數據輻射利用ENVI5.2中的輻射定標工具進行定標。利用ENVI5.2中的FLASSH模塊進行大氣校正，該模塊要求同時對可見光近紅外波段與短波紅外波段進行處理。因此，將ASTER數據的SWIR波段空間分辨率重采樣到15 m，與VNIR波段一致，并與VNIR波段按順序進行打包，進行大氣校正[25]。

本次使用的WorldView-3遙感數據級別為L2A級，在數據生產時，已進行幾何校正，消除了由于傳感器的影響造成的輻射畸變，因此，輻射校正只需消除大氣對輻射畸變的影響。輻射校正過程主要分為兩步，分別為輻射定標與大氣校正。數據的輻射定標，利用ENVI5.2中的輻射定標工具對WorldView-3數據的VNIR和SWIR波段進行定標，將遙感影像中的像元亮度值（DN值）轉化為輻射亮度值。在大氣校正方面，利用6S大氣校正模型對VNIR和SWIR進行大氣校正，6S模型能夠模擬遙感數據在太陽-目標-傳感器的傳輸過程中受到不同情況下的大氣影響[28-30]。

3? 蝕變礦物光譜參考

3.1? 礦物指數法理論基礎

礦物光譜指數法用于礦物識別和填圖。該方法可去除影像陰影，有效增強礦物信息，是礦物信息提取的常用方法。通過分析蝕變礦物的波譜特征，找出斜率變化幅度最大的波譜范圍，明確相對意義的反射峰和吸收谷，利用反射峰除以吸收谷（當反射峰為2個時，反射峰加和乘以0.5，再除以吸收谷）形成礦物光譜指數，運算后的圖像可很好的突出礦物信息[12]。

3.2? ASTER數據光譜參考

閆柏琨和王潤生在ASTER遙感數據處理與礦物填圖技術指南中提出赤鐵礦、褐鐵礦在band1至band2之間斜率變化幅度最大。因此，赤鐵礦、褐鐵礦的礦物指數為band2/band1；黃鉀鐵礬在band6至band8之間有一個明顯的吸收谷band7，且斜率變化幅度大。因此，黃鉀鐵礬的礦物指數為（band6+band8）×0.5/band7；白云母在band5至band7之間有一個明顯的吸收谷band6，且斜率變化幅度大，因此，白云母礦物指數為（band5+band7）×0.5/band6；綠泥石、綠簾石、菱鐵礦、方解石在band7至band9之間有一個明顯的吸收谷band8。因此，綠泥石、綠簾石、菱鐵礦、方解石的礦物指數為（band7+band9）×0.5/band8（表2）。為更加直觀的展示礦物吸收反射特征，將光譜重采樣至ASTER波段（圖2-A）。

3.3? WorldView-3數據光譜參考

據礦物指數法和ASTER數據礦物指數，提出WorldView-3數據礦物指數，彌補WorldView-3數據在礦物信息提取方面的不足，為開展礦物信息提取提供參考。赤鐵礦+褐鐵礦在band3至band4之間斜率變化幅度最大。因此，赤鐵礦、褐鐵礦的礦物指數為band4/band3；黃鉀鐵礬在band14至band16之間有一個明顯的吸收谷band15，且斜率變化幅度大，因此，黃鉀鐵礬的礦物指數為（band14+band16）×0.5/band15。白云母在band13至band15之間有一個明顯的吸收谷band14，且斜率變化幅度大。因此，白云母的礦物指數為（band13+band15）×0.5/band14。與ASTER數據相比，WorldView-3數據缺少與ASTER對應的2.4 μm（ASTERband9），無法計算綠泥石、綠簾石、菱鐵礦、方解石礦物指數，因此參考Kruse F A提出的Worldview-3數據碳酸鹽指數band14∕（band15+band16）[14]（表2）。為更加直觀的展示礦物吸收反射特征，將光譜重采樣至Worldview-3波段（圖2-B-C）。

4? 蝕變礦物提取結果與討論

本文利用ENVI軟件中的Band Math工具，按上述礦物指數，計算得到ASTER和Worldview-3數據赤鐵礦、褐鐵礦；黃鉀鐵礬；白云母；綠泥石、綠簾石、菱鐵礦、方解石礦物圖像，并使用Raster Color Slices工具將圖像數據進行顏色分割，數值越大表示礦物含量越高，因此，保留數值相對高的顏色，將其轉換為shp格式，加入ArcGIS結合遙感地質解譯圖對蝕變礦物信息進行優選，最終得到ASTER和World-view-3數據蝕變礦物信息。

四類蝕變礦物信息提取結果在ASTER和World-view-3數據中分布規律大體相似，主要呈團塊狀分布在研究區西南部，呈條帶狀分布在研究區東部區域的頂部和底部。分述如下：赤鐵礦、褐鐵礦在ASTER數據中提取效果優于Worldview-3數據（礦物分布較少，下同）；黃鉀鐵礬在ASTER數據中提取效果優于Worldview-3數據；白云母在Worldview-3數據中的提取效果優于ASTER數據；綠泥石、綠簾石、菱鐵礦、方解石在Worldview-3數據中的提取效果優于ASTER數據。因此，在提取蝕變礦物信息時可結合使用ASTER和Worldview-3數據，赤鐵礦、褐鐵礦，黃鉀鐵礬可應用ASTER數據礦物指數提取蝕變礦物信息，白云母，綠泥石、綠簾石、菱鐵礦、方解石可應用Worldview-3數據礦物指數提取蝕變礦物信息（圖3-A-C）。

同時，Worldview-3數據在研究區中部偏東區域提取了呈斑狀分布、NE向展布的4類蝕變礦物組合信息，疊加效果較好（圖3）。該處蝕變礦物信息與遙感地質解譯圖中的含礦地層、石英脈、NE向斷裂構造疊加較好，該處可能有成礦遠景，由于不具備實地調查條件，無法佐證其準確性，僅做參考。

5? 結論

（1） 本文基于ASTER和Worldview-3數據，利用前人提出的礦物指數法理論基礎和ASTER數據礦物指數，分析蝕變礦物波譜特征，總結歸納提出了Worldview-3數據蝕變礦物指數，赤鐵礦、褐鐵礦ASTER和Worldview-3數據礦物指數分別為band2/band1和band4/band3；黃鉀鐵礬ASTER和Worldview-3數據礦物指數分別為（band6+band8）×0.5/band7和（band14+band16）×0.5/band15；白云母ASTER和Worldview-3數據礦物指數分別為（band5+band7）×0.5/band6和（band13+band15） ×0.5/band14；綠泥石、綠簾石、菱鐵礦、方解石ASTER和Worldview-3數據礦物指數分別為（band7+band9）×0.5/band8和band14∕（band15+band16）。通過觀察研究區蝕變礦物信息分布情況，對比分析了ASTER和WorldView-3數據蝕變礦物信息提取的適用性，ASTER數據更適用于赤鐵礦、褐鐵礦，黃鉀鐵礬礦物信息提取，WorldView-3數據更適用于白云母，綠泥石、綠簾石、菱鐵礦、方解石礦物信息提取?；贏STER和WorldView-3數據蝕變礦物指數，實現了兩種數據的優勢互補，較好地彌補了Worldview-3數據在蝕變礦物信息提取方法上的不足，為進一步開展Worldview-3數據蝕變礦物信息提取奠定基礎。

（2） 本文利用ASTER和Worldview-3數據提取的蝕變礦物信息，結合研究區內遙感地質解譯圖，地質背景、含礦地層、構造特征和石英脈分布情況，預測研究區內東北部偏中區域具有一定成礦遠景，由于條件限制，無法進行野外驗證，成礦遠景僅做推斷和參考。

參考文獻

[1] 李寧，趙龍飛.淺議現代遙感技術在地質找礦中的應用[J].科技創新與應用，2013，1（29）：25-25.

[2] Ninomiya Y.A stabilized vegetation index and several mineralogic indices defined for ASTER VNIR and SWIR data[C].Geoscience and Remote Sensing Symposium，2003，3：1552-1554.

[3] Amer R，El Mezayen A，Hasanein M.ASTER spectral analysis for alteration minerals associated with gold mineralization[J].Ore Geology Reviews，2016，1（75）： 239-251.

[4] 陳琪，趙志芳，姜琦剛，等.基于ASTER與Sentinel-2A融合數據的云南普朗銅礦化蝕變信息提取[J].地質與勘探，2021，57（4）：728-738.

[5] 胡官兵，劉舫，黨偉，等.ASTER遙感數據在云南大平掌銅礦外圍勘查中的應用[J].礦產勘查，2018，9（3）：464-469.

[6] 武慧智，陳俊魁，白朝軍，等.基于ASTER與Landsat8的遙感蝕變信息提取——以新疆西昆侖塔什庫爾干地區老并鐵礦區為例[J].礦產勘查，2019，10（10）：2667-2673.

[7] 李霖龍，李士輝，宋立， 等.ASTER數據在遙感蝕變填圖中的應用研究——以內蒙古花敖包特銀鉛鋅多金屬礦床為例[J].礦床地質，2021，40（02）：384-397.

[8] 茹菲娜·阿力木江，陳川，等.基于ASTER數據的遙感蝕變信息提取——以西天山穆龍套地區為例[J].地球學報，2022，43（02）：235-245.

[9] 宋伊圩，王鵬，連琛芹，等.基于ASTER光譜特征的巖性填圖和蝕變信息提?。耗钤鸬V例析[J].西北地質，2021，54（02）：126-136.

[10] 魯立輝，許榮科，鄭有業，等.基于OLI和ASTER數據的蝕變信息提取分析與對比——以青海都蘭通突爾鉛鋅多金屬礦區為例[J].地質與勘探，2019，55（02）：600-607.

[11] 王曦，程三友，林海星，等.基于Aster和Landsat8數據在青海賽什騰地區蝕變信息提取研究[J].地質論評，2022，68（01）：262-280.

[12] 閆柏琨，王潤生.ASTER遙感數據處理與礦物填圖技術指南[R].北京：中國國土資源航空物探遙感中心，2010.

[13] Kruse F A，Perry S L.Mineral mapping using simulated Worldview-3 short-wave-infrared imagery[J].Remote Sensing，2013，5（6）：2688-2703.

[14] Kruse F A，Baugh W M，Perry S L.Validation of DigitalGlobe WorldView-3 Earth imaging satellite shortwave infrared bands for mineral mapping[J].Journal of Applied Remote Sensing，2015，9（1）：307-310.

[15] Yaqin Sun，Shufang Tian，Baogang Di.Extracting mineral alteration information using WorldView-3 data[J].Geoscience Frontiers，2017，8（5）：1051-1062.

[16] 孫婭琴.WorldView-3數據處理與蝕變信息提取方法研究[D].中國地質大學（北京），2018.

[17] 張元濤，余長發，潘蔚，等.WorldView-3影像與ASTER熱紅外影像在內蒙古衛境地區鈾礦勘查中的應用[J].鈾礦地質，2020，36（05）：408-417.

[18] 童勤龍，葉發旺，秦明寬，等.WorldView-3遙感數據在巴音戈壁盆地本巴圖地區鈾礦勘查中的應用[J].地質論評，2022，68（06）：2349-2364.

[19] 段磊，秦友平，滕海清，等.若羌縣昌運三峰山金礦地質特征[J].新疆有色金屬，2014，37（1）：44-47.

[20] 王啟宴，李山，范家盟，等.新疆若羌縣三峰山銅礦區地質特征及礦床成因初探[J].山東國土資源，2013，29（7）：20-24.

[21] 趙月平.新疆若羌縣三峰山銅金礦床地質特征及成因分析[J].新疆有色金屬，2014，37（B05）：70-73.

[22] 徐衛東，王興保，李生全，等.新疆若羌縣大青山金礦構造控礦特征[J].地質與勘探，2006（1）：34-37.

[23] Ester Shalimba. Geology and Geochemistry of the Sanqingshan Cu-Au Deposit in Ruoqiang County， Xinjiang， China[D].China University of Geosciences （Beijing），2014.

[24] Salati S，van Ruitenbeek F，van der Meer F，et al.Detection of alteration induced by onshore gas seeps from ASTER and WorldView-2 Data[J].Remote Sensing， 2014，6（4）：3188-3209.

[25] Bedini E.Mineral mapping in the Kap Simpson complex，central East Greenland，using HyMap and ASTER remote sensing data[J].Advances in Space Research，2011，47（1）：60-73.

[26] Shahriari H，Honarmand M，Ranjbar H.Comparison of multi-temporal ASTER images for hydrothermal alteration mapping using a fractal-aided SAM method[J].International Journal of Remote Sensing，2015，36（5）：1271-1289.

[27] Wang T，Zhang H，Lin H，et al.Textural-Spectral Feature-Based Species Classification of Mangroves in Mai Po Nature Reserve from Worldview-3 Imagery[J].Remote Sensing，2015，8（1）：24-24.

[28] 秦其明，朱黎江.基于6S模型的可見光，近紅外遙感數據的大氣校正[J].北京大學學報：自然科學版，2004，40（4）：611-618.

[29] 張勇.基于6S輻射傳輸模型的大氣校正研究與應用[D].中南大學，2014.

[30] 周婷婷.遙感影像輻射校正研究與應用[D].福建師范大學，2010.

Comparative Study on Extracting Altered Mineral Information in the Pobei

Area of Xinjiang Based on ASTER and WorldView-3 Data

Wang Pingping1， Kong Lingfeng1， Zhao Hui2， Yao Hu3， Yin Yilin4

（1.Henan Provincial Natural Resources Monitoring and Land Regulation Institute，Zhengzhou，Henan，450016，China;

2.Kunming University of Science and Technology， Kunming，Yunnan，650093，China;3.Yunnan Open University，

Kunming，Yunnan，650000，China;4.The First Geological Team of Henan Province Non-ferrous Metals

Geological Mineral Resources Bureau，Zhengzhou，Henan，450016，China）

Abstract： This article used ASTER and WorldView-3 data to preprocess the data. By referencing and analyzing the mineral index of ASTER data， the theory of mineral index method， and the spectral characteristics of altered minerals， the WorldView-3 data altered mineral index was summarized and proposed. By observing the distribution of altered mineral information in the study area， the applicability of extracting altered mineral information from ASTER and WorldView-3 data was compared and analyzed. ASTER data is more suitable for extracting mineral information from hematite， limonite， and jarosite， while WorldView-3 data is more suitable for extracting mineral information from muscovite， chlorite， epidote， siderite， and calcite. At the same time， based on the geological characteristics of the study area， a prospective mineralization area was predicted.Based on ASTER and WorldView-3 data on altered mineral indices， the complementary advantages of the two types of data have been achieved， improving the extraction of altered mineral information from WorldView-3 data. The predicted mineralization prospect area provides a reference for further exploration.

Key words： ASTER; WorldView-3; Mineral spectral index; Altered mineral information; Ore-forming prospect

項目資助：云南省教育廳科學研究基金項目（2023J0800）資助

收稿日期：2023-08-01；修訂日期：2023-11-06

第一作者簡介：王平平（1990-），女，河南魯山人，碩士研究生，助理工程師，2017年畢業于中國地質大學（北京）資源與環境遙感專業，主要研究方向為遙感地質應用；E-mail： 770697872@qq.com

通訊作者：姚虎（1990-），男，博士研究生，講師，研究方向為環境地質；E-mail： dhsyaohu@163.com