陸建寶,季國松,張鐘華,劉 浩
(1.貴州省地質(zhì)礦產(chǎn)勘查開發(fā)局117地質(zhì)大隊,貴州 貴陽 550018;2.貴州省地質(zhì)礦產(chǎn)勘查開發(fā)局一○二地質(zhì)大隊,貴州 遵義 563000)
近年來,隨著礦產(chǎn)勘查工作的深入推進,找礦工作逐漸向深部進軍,找礦難度越來越大,找礦風(fēng)險越來高,因此科學(xué)開展成礦預(yù)測非常必要,那就必須采用多手段、多方法相結(jié)合[1],盡可能提取多元化成礦信息。通過目前的勘查成果,結(jié)合區(qū)內(nèi)開展的礦產(chǎn)調(diào)查工作。本文通過成礦構(gòu)造及成礦結(jié)構(gòu)面分析,結(jié)合地球化學(xué)信息提取、地球物理信息提取、遙感蝕變信息提取,科學(xué)分析各項手段方法的應(yīng)用效果及適用性,為區(qū)內(nèi)下步找礦研究優(yōu)選找礦方法的組合,為區(qū)內(nèi)礦產(chǎn)調(diào)查及深部勘查提供依據(jù)。
研究區(qū)位于貴州西南部,大地構(gòu)造上位于揚子陸塊區(qū)上揚子陸塊之南盤江—右江裂陷盆地北部[2-3]。區(qū)域上地層以沉積巖為主,區(qū)域巖漿巖不甚發(fā)育。
研究區(qū)內(nèi)中二疊世為開闊臺地碳酸鹽巖沉積為主,晚二疊世為潮坪相含煤細碎屑巖沉積;早三疊世印度階以淺海陸棚相為主,早三疊奧倫階以開闊臺地—局限臺地相碳酸鹽巖沉積為主,中三疊世以潮坪相碳酸鹽巖沉積為主;晚三疊世至早期以陸棚邊緣盆地相碎屑巖沉積為主,晚三疊世卡尼階中期為淺海陸棚相濁流沉積。白堊系分布局限,為內(nèi)陸河湖相砂礫巖沉積[4-6]。
研究區(qū)構(gòu)造以NE向、NW向、近EW向、近SN向四組,區(qū)域上礦產(chǎn)以低溫?zé)嵋盒?微細粒浸染型)金礦為主,金礦多產(chǎn)于大型背斜軸部有利成礦結(jié)構(gòu)面及成礦構(gòu)造中,主要包括斷裂裂隙、界面。研究區(qū)內(nèi)已發(fā)現(xiàn)了戈塘大型金礦、萬人洞、海子等小型金礦點(圖1),金礦主要受上二疊統(tǒng)龍?zhí)督M與中二疊統(tǒng)茅口組不整合界面控制[7]。賦礦巖石以細沙巖、粉砂巖、灰?guī)r及其形成的角礫巖為主。圍巖蝕變以硅化、黃鐵礦化、碳酸鹽化為主[8]。研究區(qū)成礦結(jié)構(gòu)面主要為Si/Ca面,其物性明顯差異,可通過物探方法提取有利成礦部位。區(qū)域上金礦主成礦元素Au與As、Sb、Hg等元素組合為特征,可根據(jù)元素性質(zhì)采用地球化學(xué)提取找礦信息[9-11]。
圖1 戈塘地區(qū)地質(zhì)圖Fig.1 The geological map of Getang area
研究區(qū)主要以布格重力剩余負值異常為主,在區(qū)內(nèi)北東側(cè)及西南側(cè),布格重力剩余異常軸向以NW向為其主軸方向;在中部,布格重力剩余異常較為凌亂,呈團塊狀展布,局部呈近NE向展布,從布格重力剩余異常幅值看,布格重力剩余異常極小值位于區(qū)內(nèi)西南側(cè)萬屯—魯屯一帶,向NE方向異常幅值略增大,但總體仍呈現(xiàn)出布格重力剩余負值異常,推斷區(qū)域內(nèi)引起布格重力剩余負值異常的主要因素是區(qū)域內(nèi)密度較低、規(guī)模較大、埋藏深度較淺的龍?zhí)督M引起。
基于地質(zhì)體的地球物理特性,當(dāng)斷裂產(chǎn)生后,會使地質(zhì)體在三度空間發(fā)生位移和錯斷時斷裂兩側(cè)的地層的物性產(chǎn)生變化;斷裂反映的是地層物性界面的陡變帶,斷裂的規(guī)模越大,兩種物性界面的陡變帶規(guī)模也越大、物性差異也越大,地球物理異常梯級帶形態(tài)也越明顯。
根據(jù)研究區(qū)的重力場,結(jié)合地質(zhì)以及其他資料,推斷研究區(qū)的斷裂構(gòu)造。斷裂構(gòu)造推斷與圈定的主要依據(jù)包括布格重力異常、剩余重力異常、水平方向?qū)?shù)、垂向二階導(dǎo)數(shù)、水平總梯度模等。
1)在布格重力異常平面圖及布格重力剩余異常平面圖上的標(biāo)志:走向明顯的重力梯級帶;狹長的帶狀異常;異常帶的水平錯位或平移;異常等值線的同向扭曲;重力場發(fā)生明顯變化的分界線;重力高與重力低之間的線性過渡帶。
2)導(dǎo)數(shù)水平方向?qū)?shù)是強化斷裂構(gòu)造信息的一種行之有效的方法,不同方向的水平導(dǎo)數(shù)可以突出與之垂直的異常展布特征,從而較準(zhǔn)確的推測斷裂構(gòu)造在地面的投影位置。在布格重力水平方向?qū)?shù)異常平面圖上的標(biāo)志:線性排列的同號異常的極值連線;走向明顯的線性異常軸連線;線性異常的平移或錯位。
3)在重力垂向二階導(dǎo)數(shù)異常平面圖上的標(biāo)志:線性異常邊界的平直零值線。
4)水平總梯度重力異常計算是將重力梯級帶轉(zhuǎn)換為梯度值,其極值帶更好地對應(yīng)了斷裂位置,從而提高了對斷裂的平面分辨能力。重力水平總梯度異常極大值位置標(biāo)示著斷裂的位置,其幅值大小反映了斷裂的規(guī)模,極大值走向突變和錯斷代表斷裂被切割和錯。在重力水平總梯度模異常平面圖上的標(biāo)志:線性排列的極大值連線;極大值走向的突變和錯斷。
5)斷裂產(chǎn)生,使地質(zhì)體在三度空間發(fā)生位移和錯段,地層之間的密度產(chǎn)生變化,斷裂反應(yīng)的是地層密度界面的陡變帶,表現(xiàn)為重力異常梯級帶形態(tài)。通過異常梯級帶濾波增強技術(shù),可以突出重力梯級帶變化的效果,然后進行水平總梯度處理,能更準(zhǔn)確有效的提取斷裂信息。
根據(jù)上述原則并結(jié)合地質(zhì)以及其他資料,推斷研究區(qū)的斷裂構(gòu)造8條(圖2、表1),布格重力異常主要表現(xiàn)為密集重力異常梯級帶,剩余重力異常主要表現(xiàn)為正負異常的分界線;重力水平一階導(dǎo)數(shù)45°方向上表現(xiàn)為串珠狀的負異常帶;重力垂向二階導(dǎo)數(shù)異常上表現(xiàn)為零值線的連線;重力水平總梯度模呈串珠狀重力高值異常帶等特征,推斷解譯標(biāo)志明顯。其中,推斷的F2、F3深大斷裂分別控制了上二疊統(tǒng)厚度及沉積相的變化,可能為同沉積斷裂。F4斷裂與地表出露海馬谷斷裂大致一致,為一顆區(qū)域性大斷裂,區(qū)內(nèi)內(nèi)已發(fā)現(xiàn)的金礦床(點)分布于該斷裂2 km影響范圍內(nèi),可能為區(qū)內(nèi)的導(dǎo)礦構(gòu)造。
表1 重力推斷主要斷裂統(tǒng)計Table 1 List of the major faults inferred by gravity anomaly
圖2 研究區(qū)布格重力異常及推斷斷裂圖Fig.2 Bouguer gravity anomaly and inferred fault diagram in the study area1—布格重力異常等值線(m/s2) 2—重力推斷斷裂及編號
本次對異常及蝕變帶采用可控源音頻大地電磁測深測量(CSAMT)。以WK4線為例,根據(jù)數(shù)據(jù)分析顯示,有效探測深度達到了1200 m。結(jié)合地質(zhì)圖與物性表可以推斷1000/WK4~2480/WK4點之間高程1100 m處為茅口組與龍?zhí)督M不整合界面。2600/WK4號點附近的高低阻梯級帶傾向與地質(zhì)上的上河壩斷裂類似,推測該異常是南東傾向的上河壩斷裂的反應(yīng)。推斷2600/WK4~3560/WK4號點之間的高程1250 m處為茅口組與龍?zhí)督M不整合界面。3560/WK4號點附近的臺階狀異常其傾向與魯溝斷裂的傾向一致,推測此異常為魯溝斷裂的反應(yīng)。3560/WK4號點高程1200 m龍?zhí)督M與茅口組不整合界面(圖3)。
圖3 WK4線電阻率斷面與地質(zhì)剖面對照圖Fig.3 The comparison map of resistivity section and geological section in WK4 line1—下三疊統(tǒng)飛仙關(guān)組 2—上二疊統(tǒng)龍?zhí)督M 3—中二疊統(tǒng)茅口組 4—粉砂巖5—黏土巖 6—粉砂質(zhì)黏土巖 7—灰?guī)r 8—正斷層 9—逆斷層
綜上,剖面由北至南橫穿上河壩斷層、魯溝斷層,物探成果總體上能反映該區(qū)域的地質(zhì)、構(gòu)造特征。上河壩斷裂表現(xiàn)為高低阻梯級帶,傾向SE,傾角70°~80°,延伸400余米。魯溝斷裂表現(xiàn)與上河壩斷裂類似,延伸500 m,斷裂表現(xiàn)與地質(zhì)較吻合。同時綜合推斷出深部有利成礦構(gòu)造及成礦有利界面[9],即龍?zhí)督M與茅口組不整合界面,對尋找隱伏的主要巖性突變界面具有較好指示作用。該方法對提取0~500 m的成礦構(gòu)造及成礦結(jié)構(gòu)面有較好的指示作用。
對研究區(qū)成礦有利地段開展了寬頻大地電磁測深。以WK9線為例,WK09線反演結(jié)果所示(圖4),從電性結(jié)構(gòu)上,以靠近小海子附近的122號測點為界,在高程500 m以上的淺部地層的左右兩側(cè)可分電性差異明顯的塊體,右側(cè)表現(xiàn)為連片的高阻塊體,推斷為二疊系茅口組碳酸巖地層;左側(cè)可再分為上下2個差異明顯的塊體,上部整體表現(xiàn)為中高電阻率塊體,推斷為下三疊統(tǒng)飛仙關(guān)組(T1f)灰?guī)r,下部整體表現(xiàn)為低電阻率塊體,推斷為下三疊統(tǒng)飛仙關(guān)組—上二疊統(tǒng)龍?zhí)督M(T1f~P3l)的碎屑巖。
圖4 WK9線寬頻電阻率斷面與推斷解釋圖Fig.4 The inferred interpretation diagram and wideband resistivity section in WK9 line1—下三疊統(tǒng)水寧鎮(zhèn)組 2—下三疊統(tǒng)飛仙關(guān)組 3—上二疊統(tǒng)龍?zhí)督M 4—中二疊統(tǒng)茅口組5—實推測地質(zhì)界線 6—實推測斷層 7—實測點位置 8—施工驗證鉆孔
根據(jù)WK9線的地質(zhì)推斷解譯結(jié)果所示,剖面128測點深部(高程1200~1400 m)位置存在的一個相對低阻帶應(yīng)為多期次活動海馬谷主斷層(F1)顯示,根據(jù)電阻率梯度帶的變化特征,推斷海馬谷斷層的總體傾向為SE方向,傾角非常陡立,大致為70°~80°,斷層的性質(zhì)為逆斷層。124號測點左側(cè)高程400~500 m處的低、高電阻率梯度帶推斷為海馬谷斷層下盤P3/P2有利賦礦層位的顯示,該P3/P2界面向北存在逐漸抬升的趨勢,其對應(yīng)的高低電阻率梯度帶顯示不連續(xù),存在局部的“凹凸”地形的彎曲,在剖面110~112號測深位置“隱約”存在一個相對“隆起點”,其埋深在高程約500 m處,距離地表約1000 m。由于F4逆斷層的作用,對應(yīng)100~104號測深點整體位于逆斷層的下盤,相對于上盤,其深部P3/P2界面存在相對向下的距離位移。
通過寬頻大地電磁測深,對于埋藏在地表1 km以下的構(gòu)造和界面均反映比較清楚,綜合推斷出深部有利成礦構(gòu)造及成礦有利界面,包括海馬谷斷裂及龍?zhí)督M與茅口組不整合界面,對尋找隱伏的主要巖性突變界面具有較好指示作用。
在研究區(qū)內(nèi)橫穿海馬谷斷裂及萬人洞斷層做了復(fù)電阻率法剖面測量(CR),對數(shù)據(jù)做2.5D反演,形成的電阻率(ρ)斷面圖和極化率(IP)斷面圖(圖5)。WF1線在1400/WF1~1500/WF1點間和3100/WF1~3200/WF1點間存在兩處低阻異常帶。這兩處異常分別對應(yīng)海馬谷斷裂帶和萬人洞斷裂帶。所以,推斷此兩處異常為海馬谷斷裂和萬人洞斷裂引起。兩低阻帶都有較高的極化率異常對應(yīng)。說明在兩條斷裂帶中含有高極化物體。就調(diào)查區(qū)而言,高極化物有兩種:一種為金屬礦物,一種為含煤(或含炭質(zhì))地層。對比地質(zhì)剖面圖,可以發(fā)現(xiàn)在斷裂帶中含有明顯礦化物,在海馬谷斷裂南東側(cè)附近出露有龍?zhí)督M,而萬人洞附近主要為茅口組和棲霞組。結(jié)合物性資料,認(rèn)為1400/WF1~1800/WF1點間的高極化異常為斷裂帶中礦化體和龍?zhí)督M含煤地層共同形成,3000/WF1點以南的高極化異常為礦化異常的可能性較高。通過復(fù)電阻率法剖面測量測深,對于斷裂構(gòu)造比較清楚的反映,同時對于含金屬礦物的礦種具有直接指導(dǎo)意義。
圖5 WF1線 CR法反演電阻率斷面、極化率斷面與地質(zhì)剖面對照圖
搜集研究區(qū)1∶5萬土壤地球化學(xué)資料數(shù)據(jù)分析,圈定Au單元素異常18個(圖6)。其中Au1單元素異常規(guī)模大,強度高,峰值高;Au5、Au6、Au7、Au8、Au9、Au10、Au11、Au15等7個Au單元素異常具有一定規(guī)模和強度(均達三級);Au2、Au3、Au12、Au13、Au18等5個Au單元素異常具有一定規(guī)模,但異常強度較低(以一級為主,少部分達二級)。其余Au14、Au16、Au17等3個Au單元素異常強度低,規(guī)模較小。
圖6 研究區(qū)土壤地球化學(xué)圖
引起區(qū)內(nèi)Au元素異常的主要因素:① 產(chǎn)于茅口組與龍?zhí)督M界面的蝕變體,使其形成規(guī)模大、強度高、峰值高的異常;② 產(chǎn)于茅口組頂部巖溶洼地、漏斗控制的第四系風(fēng)化殘積型金礦,使其在茅口組中形成高背景、規(guī)模大、強度高、峰值高的異常;③ 成巖作用期或成巖期后,構(gòu)造及熱液活動引起地層局部Au元素富集形成高背含量或弱異常;④ 表生風(fēng)化淋漓作用使Au元素在土壤及半風(fēng)化巖石中富集形成弱異常。
通過土壤地球化學(xué)異常分析,結(jié)合工程揭露情況,土壤地球化學(xué)信息提取方法對于尋找氧化礦具有直接指示意義,對于尋找埋藏較淺的原生礦有較好的指示信息。
為了進一步查證引起土壤化探異常的原因,確定異常源,針對海馬谷斷裂、萬人洞斷層、干寨斷層等NE向蝕變帶試驗性開展了構(gòu)造地球化學(xué)采樣測試,針對性的沿斷裂蝕變帶采集蝕變礦物、方解石細脈、斷層泥等樣品作測試分析。以HG01為例,對海馬谷斷裂北西盤圈出的低緩異常(Au元素異常峰值為15×10-9)開展了構(gòu)造地球化學(xué)采樣測試分析,其Au元素異常峰值為146.2×10-9,與As、Sb、Hg、W、Mo元素套合較好(圖7)。進一步確定異常源及靶位,同時進一步證實了異常是由含礦熱液引起。目前通過鉆孔驗證,揭露出含金蝕變帶,該方法采樣介質(zhì)為與成礦有關(guān)的物質(zhì)或蝕變礦物,相對于土壤具有直接指導(dǎo)作用,同時強化了弱異常信息的提取,對于基巖出露較好的區(qū)域?qū)ふ覠嵋盒徒鸬V具有較強的指導(dǎo)作用。不但強化了異常或找礦信息,而且大大減少了采樣量和分析費用,提高了工作效率和找礦效果。
圖7 HG01構(gòu)造地球化學(xué)剖面與土壤地球化學(xué)異常疊合圖Fig.7 Superposition diagram of geochemical structure profile and geochemical soil anomaly in HG01 line1—下三疊統(tǒng)水寧鎮(zhèn)組第三段 2—水寧鎮(zhèn)組第二段 3—水寧鎮(zhèn)組第一段 4—下三疊統(tǒng)飛仙關(guān)組第四段5—地質(zhì)界線 6—斷層 7—構(gòu)造地球化學(xué)采樣點 8—土壤化探異常(Au/10-9)
地電化學(xué)測量是以地下巖石中的離子動態(tài)平衡狀態(tài)為基礎(chǔ)的地球化學(xué)方法。礦石的電化學(xué)溶解作用及溶解了的離子,在電流驅(qū)動下向上遷移并在測點處形成積累,通過測定積累隨時間的變化規(guī)律,即可指示隱伏礦體賦存位置。在人工電流場的作用下,使深部處在動態(tài)平衡的離子產(chǎn)生兩極分化,陰陽離子分別向無窮遠極和接收極遷移。離子的遷移并非指從深部直接遷移至地表,而是相鄰單元體積內(nèi)的離子沿電流場的方向移出和移入,并最終形成接收器中的陽離子積累。接收器中陽離子積累的濃度與供電時間、介質(zhì)中離子的濃度及離子的遷移速度成正消長關(guān)系。
地氣的概念是1982年由瑞典蘭德大學(xué)物理系Kristiansson博士和波立登礦業(yè)公司勘探部Malmqvist博士共同提出的,我國童純菡(1989)首次開始了地氣法隱伏金礦的試驗,隨著王學(xué)求等人對深穿透地球化學(xué)的研究表明,地球內(nèi)部存在著垂直運移的上升氣流,當(dāng)它流經(jīng)礦體或巖層時,將其中元素的納米微粒攜帶并遷移至地表,從而在礦體上方形成了成礦元素、伴生元素的地氣異常。
橫穿干寨斷裂開展了地電化學(xué)及地氣測量,以HG03為例,提取了各元素不同程度的異常,且異常明顯,異常重合性較好。尤其是19~21號采樣點上的Au元素呈單峰異常,異常寬度為60 m,異常峰值為178×10-9,與As、Hg元素組合異常近似完全套合(圖8),其與地氣剖面和快速測試分析儀器圈定的異常套合較好,地電化學(xué)異常強度較大。通過對異常進行查證,在干寨斷層附近發(fā)現(xiàn)一條礦化蝕變帶,經(jīng)取樣分析具金礦化顯示。該方法對于深覆蓋區(qū),尋找隱伏礦床具有較好的實用性[12]。
圖8 HD3線地電、地氣化學(xué)剖面與地質(zhì)剖面對照圖Fig.8 The comparison diagram of geoelectricity and gas geochemical profile and geological profile in HD3 line1—下三疊統(tǒng)水寧鎮(zhèn)組第三段 2—水寧鎮(zhèn)組第二段 3—水寧鎮(zhèn)組第一段 4—下三疊統(tǒng)飛仙關(guān)組第四段5—砂巖 6—黏土巖 7—粉砂質(zhì)黏土巖 8—灰?guī)r 9—正斷層 10—含金礦化蝕變帶
遙感礦化蝕變信息提取,是利用遙感圖像處理技術(shù),直接從遙感圖像數(shù)據(jù)中提取反映礦化圍巖蝕變的波譜信息,獲取區(qū)域羥基及鐵染等與礦化密切相關(guān)的資料信息[13],為進一步確定找礦靶區(qū)提供依據(jù)。
對確定的蝕變異常分量,參考化探異常分級的辦法,以σ(標(biāo)準(zhǔn)離差)做為尺度,用數(shù)倍σ值做為閾值,限定異常水平。將主分量分析結(jié)果的均值(X)理解為代表區(qū)域背景,利用(X+kσ)劃分異常強度等級。本次k值分別取4、3.5、3時,羥基異常從強到弱分為1~3級;當(dāng)k值分別取5、4.5、4時,鐵染異常從強到弱分為一級、二級、三級。
研究區(qū)內(nèi)的金礦礦化蝕變類型主要有硅化、黃鐵礦化、碳酸鹽化、石膏化、黏土化等。其中黃鐵礦化主要通過鐵染(Fe3+)異常提取來反應(yīng),而碳酸鹽化礦物含有羥基(—OH),主要通過羥基(—OH)異常提取來反應(yīng)碳酸鹽化。采用TM數(shù)據(jù)的(1,4,5,7)波段和TM(1,3,4,5)波段進行主成分分析,分別獲取羥基異常及鐵染異常信息,根據(jù)其中特征向量特點,分別選取指示性吻合的主成分分量PC4為羥基蝕變異常及鐵染蝕變異常,并通過設(shè)定閥值,限定異常水平,對進行異常分級圈定,羥基異常與鐵染異常從強到弱均分為一級、二級、三級[15]。
遙感蝕變異常主要分布在戈塘背斜南東翼,位于海馬谷斷層和魯溝斷層之間(圖9),與解譯的線性構(gòu)造和環(huán)形構(gòu)造的相關(guān)度較高,也與已知金礦點吻合度高。羥基與鐵染蝕變異常主要分布在二疊系茅口組、龍?zhí)督M中,與斷層及線性構(gòu)造呈一定的相關(guān)性。
圖9 研究區(qū)遙感蝕變異常分布圖Fig.9 Distribution map of remote sensing alteration anomalies in the study area1—中三疊統(tǒng) 2—下三疊統(tǒng) 3—上二疊統(tǒng) 4—中二疊統(tǒng) 5—地層界線 6—正斷層 7—逆斷層 8—斷層9—一級鐵染異常 10—二級鐵染異常 11—三級鐵染異常 12—一級羥基異常 13—二級羥基異常 14—三級羥基異常
對研究區(qū)開展方法性實驗,研究區(qū)內(nèi)地質(zhì)、地球物理、地球化學(xué)、礦化蝕變等多元地學(xué)數(shù)據(jù)的找礦信息提取和綜合處理技術(shù)研究,分析其方法應(yīng)用效果及其適用性,形成了一套多元信息綜合找礦的技術(shù)。得出的結(jié)論如下:
1)通過物探方法試驗,重力能解決深大斷裂或隱伏巖體,即可以推斷成礦地質(zhì)體或流體運移通道位置;對于淺部構(gòu)造可控源音頻大地電磁測深測量能較好的反映斷裂構(gòu)造規(guī)模及延伸、巖性界面埋深及起伏情況;對于深部構(gòu)造采用寬頻大地電磁測深測量能較好的反映斷裂構(gòu)造規(guī)模及延伸、巖性界面埋深及起伏情況。
2)通過構(gòu)造地球化學(xué)的新方法試驗工作,相對于土壤地球化學(xué)能直接提取與成礦相關(guān)的找礦信息,能較好的提取弱異常信息,直接指導(dǎo)找礦靶位,對尋找隱伏礦體有一定指導(dǎo)作用。通過地電、地氣地球化學(xué)的新方法試驗工作,發(fā)現(xiàn)地電、地氣化學(xué)剖面測量能提取與成礦有關(guān)的指示信息,強化了弱異常,對強覆蓋區(qū)尋找隱伏礦體有一定指導(dǎo)作用。
3)研究區(qū)內(nèi)的金礦礦化蝕變類型主要有硅化、黃鐵礦化、碳酸鹽化、石膏化、黏土化等。其中黃鐵礦化主要通過鐵染(Fe3+)異常提取來反應(yīng),而碳酸鹽化礦物含有羥基(—OH),主要通過羥基(—OH)異常提取來反應(yīng)碳酸鹽化。采用TM數(shù)據(jù)進行主成分分析,可以獲取羥基異常及鐵染異常信息,可以較好的指示礦化蝕變情況及分布范圍。