徐龍生,路 瀚,黃天統(tǒng),苑東洋
(1. 核工業(yè)二三〇研究所,湖南 長沙 410000; 2. 中國地質(zhì)大學(武漢)地球物理與空間信息學院,湖北 武漢 430074; 3. 中國冶金地質(zhì)總局西北地質(zhì)勘查院,陜西 西安 710119)
多元統(tǒng)計分析是從經(jīng)典統(tǒng)計學中發(fā)展起來的一個分支,是一種綜合分析方法,它能夠在多個對象和多個指標互相關(guān)聯(lián)的情況下分析統(tǒng)計規(guī)律,很適合地球科學研究的特點。在地質(zhì)科學中應(yīng)用比較多的有相關(guān)分析、回歸分析、聚類分析、多元方差分析、主成分分析、因子分析、判別分析與聚類分析等。
統(tǒng)計分析是系統(tǒng)論述統(tǒng)計理論和方法的科學,因其具有較好的實用性,可應(yīng)用到地球物理、地球化學研究中。許多學者[1-6]將多元統(tǒng)計分析應(yīng)用到地球物理、地球化學、水文地質(zhì)、煤炭地質(zhì)、石油地質(zhì)等地球科學領(lǐng)域中,并取得了非常好的效果。通過統(tǒng)計方法在地球物理與地球勘查中的應(yīng)用,展現(xiàn)出統(tǒng)計方法在相關(guān)地球物理與地球化學研究中取得的成果,同時驗證多元統(tǒng)計方法的優(yōu)越性和應(yīng)用廣泛性。
目前的多元統(tǒng)計分析作為一類應(yīng)用廣泛、新興的數(shù)據(jù)處理方法,在各個學科中均得到較好應(yīng)用,其并沒有展現(xiàn)出復雜系統(tǒng)中結(jié)構(gòu)形成和演化的動力學特性,亦沒有探索出復雜系統(tǒng)中的地球物理學與地球化學原理[7]。
應(yīng)用多元統(tǒng)計分析方法(回歸分析、相關(guān)分析、聚類分析、因子分析、判別分析)基本思想是將研究區(qū)所有局部異常(樣品)進行分類,根據(jù)分類后某些異常對應(yīng)于某些已知礦(已知樣品),則同類中其他異常(未知樣品)也是該已知礦母體的反映,從而作出成礦預(yù)測[8]。
(1)
式(1)中:n為樣本量,xi、yi、x、y分別為兩個變量的觀測值和均值。若b→0,表明兩變量是正相關(guān);若b<0,表明兩變量是負相關(guān)。另外,b的大小還表明若變量x發(fā)生一個單位變化時,變量y的變化量對于總體回歸系數(shù)為0的原假設(shè),也可以用t統(tǒng)計量進行檢驗。當t檢驗顯著時,拒絕原假設(shè),即總體回歸系數(shù)不為0,兩個變量呈線性關(guān)系;當t檢驗不顯著時,不能拒絕原假設(shè),兩個變量不是線性相關(guān)。
在研究重力異常和高程的關(guān)系時,常常用回歸分析來研究高程對重力布格異常的影響。在高海拔地區(qū),尤其是青藏高原地區(qū)這種影響特別明顯,故而一元線性回歸分析就有很重要的意義。
相關(guān)分析是研究兩個或兩個以上處于同等地位的隨機變量間相關(guān)關(guān)系統(tǒng)計分析的方法。相關(guān)分析與回歸分析之間的區(qū)別:回歸分析側(cè)重于研究隨機變量間的依賴關(guān)系,以便用一個變量去預(yù)測另一個變量;相關(guān)分析側(cè)重于發(fā)現(xiàn)隨機變量間的相關(guān)特性[9]。相關(guān)分析在地質(zhì)礦產(chǎn)勘查中,尤其在地球物理、地球化學等方面都有廣泛應(yīng)用。
相關(guān)系數(shù)是按積差方法計算,是以兩個變量與各自平均值的離差為基礎(chǔ),通過兩個離差相乘來反映兩個變量之間的相關(guān)程度。重點研究線性的單相關(guān)系數(shù),確定相關(guān)關(guān)系的數(shù)學表達式,確定因變量估計值誤差的程度[10]。
相關(guān)分析現(xiàn)象之間的相互聯(lián)系,經(jīng)常表現(xiàn)為一定關(guān)系,其中一個或若干個起著影響作用的變量用x表示,與之相關(guān)的另一現(xiàn)象變量用y表示。通常用相關(guān)系數(shù)r表征這兩個量的相互關(guān)聯(lián)程度,計算公式[11]:
(2)
相關(guān)系數(shù)r的取值范圍:-1≤r≤1。r的絕對值愈大,變量間的等級相關(guān)程度愈大[1];r>0為正相關(guān),r<0為負相關(guān);|r|=0表示不存在相關(guān)關(guān)系;|r|=1表示完全線性相關(guān);0<|r|<1表示存在不同程度線性相關(guān)。
在地球科學研究應(yīng)用中,|r|≤0.3為不存在線性相關(guān);0.3<|r|≤0.5為低度線性相關(guān);0.5<|r|≤0.8為顯著線性相關(guān);|r|→0.8為高度線性相關(guān)。
統(tǒng)計分析方法就是運用數(shù)學思維方法,構(gòu)建一定的數(shù)學模型,分析通過各種途徑獲取的不同數(shù)據(jù)及資料進行數(shù)理統(tǒng)計研究,形成定量解釋成果。統(tǒng)計分析方法是廣泛使用的現(xiàn)代科學方法,是一種較科學、精確和客觀的測評方法。
成秋明院士提出多重分形與地質(zhì)統(tǒng)計學方法用于勘查地球化學異??臻g結(jié)構(gòu)和奇異性分析??辈榈厍蚧瘜W和地球物理場的局部空間結(jié)構(gòu)變化性包括空間自相關(guān)性以及奇異性,前者可通過地質(zhì)統(tǒng)計學中常用的變異函數(shù)來實現(xiàn);后者可用多重分形模型進行刻畫,具有自相似性或統(tǒng)計自相似性的多重分形分布的奇異性可以反映地球化學元素在巖石等介質(zhì)中的局部富集和貧化規(guī)律[2]。而多重分形插值和估計方法可以同時度量以上兩種局部結(jié)構(gòu)性質(zhì)(空間自相關(guān)性以及奇異性),不僅能夠進行空間數(shù)據(jù)插值,同時還能保持和增強數(shù)據(jù)的局部結(jié)構(gòu)信息,對于地球化學和地球物理異常分析和識別具有巨大的引用價值[12]。
因此,地統(tǒng)計學近年來成為空間分析的重要工具,其應(yīng)用領(lǐng)域廣泛分布在自然科學的眾多領(lǐng)域[13]。
物探工作中地形對重力異常的影響不容忽視。地形的隆起或是凹陷,將導致重力異常圖在形態(tài)上的直接變化,影響對異常的客觀解釋。
將布格重力異常和地形兩個物理量進行相關(guān)分析,可用來反映研究區(qū)內(nèi)地形起伏對重力異常的影響情況,地形改正是否完全以及改算方法是否完善[14]。這里所需要的重力異常數(shù)據(jù)文件和高程數(shù)據(jù)文件必須是以相同的網(wǎng)格間距網(wǎng)格化的、具有相同橫縱坐標的“*grd”網(wǎng)格化文件。
山西平順縣龍溪鎮(zhèn)新天嶺礦業(yè)鐵礦勘查區(qū)地處太行山斷隆南部,太行山復式背斜西翼,最高海拔1 782 m,最低海拔1 393 m,在研究區(qū)約8 km2的范圍,地形高程變化范圍約389 m,地形較為復雜(見圖1a),地形改正是研究區(qū)較為重要的工作。本次地形改正使用收集的1∶10 000的DEM高程數(shù)據(jù)進行中區(qū)地改,遠區(qū)地改使用1∶50 000的DEM高程數(shù)據(jù),最大限度地減小地形變化對重力場的影響,增強地改效果。研究區(qū)南北部是重力異常富集區(qū),因而在研究區(qū)南北部是各截取一個剖面(見圖1b)分別進行了5 km和20 km的地形改正,用于探討本次地改效果(見圖2)。
圖1 研究區(qū)海拔高程地形(a)、研究區(qū)地改剖面位置(b)
圖2 1-1′剖面地改綜合異常(a)、2-2′剖面地改綜合異常(b)
經(jīng)過地形相關(guān)分析可知,1號和2號剖面地形線(黑色線)均顯示單邊上行的形態(tài),不同半徑(5 km和20 km)的地形改正曲線均反映出地形改正的削高填低的作用,而且小半徑地改的這種作用更明顯,因此20 km地改半徑后的布格重力異常包含了更多區(qū)域重力場,其剩余異常特別是重力高的范圍相對較大,反映了更多的深部地質(zhì)體的作用。由前述地質(zhì)成礦的描述可知,研究區(qū)礦體和巖體范圍和埋深有限,礦體小、多為豆莢狀或透鏡體,結(jié)合兩個剖面的地改計算對比結(jié)果來看,選取5 km的地改半徑比較合適。從研究區(qū)來看,研究區(qū)地形圖(見圖1a)出現(xiàn)東西兩側(cè)地勢低、中部明顯隆起的形態(tài),經(jīng)過地形改正及其他各項改正后的布格異常圖顯示了自西向東逐漸升高的重力場特點,地形和布格重力場沒有出現(xiàn)明顯的相關(guān)性,本次地改較為有效地消除了地形對重力場的影響,而且經(jīng)過半徑為5 km地形改正后的布格重力異常,更多地顯示出淺部地質(zhì)體的重力場,壓制了深部因素引起的重力場。
根據(jù)布格異常與地形的相關(guān)分析,遠區(qū)地改做5 km的改正半徑較為合理。背景重力場整體是東高西低,平均場值為-1.8×10-5m/s2,變化范圍為-2.8×10-5~0.8×10-5m/s2。研究區(qū)布格重力場等值線見圖3。
圖3 布格重力異常
相關(guān)分析得到重力場具有顯著特點,中部重力場以梯級帶形式將研究區(qū)重力場一分為二,西低東高,單向遞增,重力異常梯級帶走向基本是北北東走向。由前述可知研究區(qū)發(fā)育軸向為北北東走向的背斜,背斜軸部大致位于重力梯級帶的東部。這種重力場特點顯然反映出了研究區(qū)的背斜地質(zhì)構(gòu)造形態(tài)。西部重力低是傾向為西的背斜西翼的反映,體現(xiàn)出夾持在第一和第二巖漿巖帶之間的向斜構(gòu)造形態(tài);東部重力高則反映出傾向東的背斜東翼,但東部整個重力高值區(qū)很明顯具有南北兩種形態(tài)的重力場,南部異常形態(tài)較為規(guī)整,賦值也高,可以看作一個較為完整的異常;北部異常較為凌亂賦值也相對低,基本屬于低緩異常,突出了背斜東翼南北的地質(zhì)構(gòu)造包括礦產(chǎn)地質(zhì)等差別較大。重力場最低處位于研究區(qū)西北角,最低值為-2.8×10-5m/s2,反映出向斜構(gòu)造向西延伸,其核部位于研究區(qū)外圍的西側(cè);位于研究區(qū)東南部的重力場最高值為0.8×10-5m/s2,重力值最高處周圍的重力場為環(huán)形異常,多個局部異常等值線圈閉范圍較大,是研究區(qū)異常值變化最為劇烈的部位,位于鐵礦05、06井口處及其以北地區(qū),很明顯是研究區(qū)已知礦體的體現(xiàn)。
重力勘探是基于密度差異而進行的,掌握研究區(qū)內(nèi)各地層及巖石的密度特征是重力解釋與推斷工作的基礎(chǔ)。研究區(qū)巖石種類較多,主要以矽卡巖為主,其次為大理巖、二長花崗巖等。工作期間未對巖性開展密度樣采集工作,主要對礦化的矽卡巖進行了系統(tǒng)的密度樣采集工作,對不同巖礦石的不同地球物理參數(shù)進行統(tǒng)計分析。
應(yīng)用多元統(tǒng)計分析方法,主要用統(tǒng)計分析、正態(tài)分布等方法,綜合統(tǒng)計表中的統(tǒng)計數(shù)據(jù)分析,繪制了物性參數(shù)統(tǒng)計曲線圖(見圖4),鐵多金屬礦矽卡巖和磁鐵礦矽卡巖磁性強高于其他巖石,矽卡巖性礦體存在高磁特性。研究區(qū)基本以古生界地層為主要目標地層,最老的志留奧陶系灘間山群(OT3)至最新的二疊系打柴溝組(P1-2d),老地層的花崗巖基本具有無磁性、低極化率、高電阻率、低密度的特性,只有電性參數(shù)與矽卡巖礦石頭一致,其他均不一樣。而重點地層上石炭統(tǒng)締敖蘇組(C2d)是矽卡巖礦化的重點層位。整個底層從老至新,磁性從無到強,電阻率從高到低,密度由低至高,綜合來看,整體體現(xiàn)了青藏高原巖石圈板塊具有輕殼重幔,軟殼硬幔,熱殼冷幔的地球物理特征。
圖4 研究區(qū)電阻率、密度、極化率的平均值
綜上所述,牛苦頭地區(qū)的矽卡巖型礦床具有“三高一低”的物性特征,即強磁特征、高密度、高電阻率、低極化率,為找礦提供了巖礦石物性基礎(chǔ)。
3.3.1 重力異常與測點高程的相關(guān)性分析
采用一元線性相關(guān)分析法,對研究區(qū)布格重力異常與測點高程進行了相關(guān)分析。重力異常與高程之間的關(guān)系一般用線性回歸方程表示為:
G=k×H+b
(3)
式(3)中:G為重力異常值;H為海拔高程,m;k為回歸方程的斜率。它反映重力值(G)對高程(H)的敏感程度。k為正,說明高度改正值大于布格改正值,或密度改正值小于中間層實際的平均密度改正值2.67 g/cm3;k為負,則高度改正值小于布格改正值,或密度改正值大于中間層實際的平均密度改正值2.67 g/cm3。b為方程的截矩,反映整個地殼內(nèi)物質(zhì)平均密度的相對(相對于其他塊體)大小?;貧w方程的可靠性用相關(guān)系數(shù)R表示。
3.3.2 重力場與高程的相關(guān)性分析及地質(zhì)意義
布格重力異常與高程的統(tǒng)計分析(見圖5),從一個方面反映了地殼內(nèi)的結(jié)構(gòu)特征,為推斷解釋提供統(tǒng)計信息。圖5為中間層密度2.67 g/cm3,高度起算面不同布格重力異常值與高程相關(guān)性散點圖,圖5a為改正到大地水準面的布格重力異常值與高程回歸分析散點圖,二者為正相關(guān),斜率0.20,相關(guān)系數(shù)為0.64,截距為-506.513;圖5b以羅布莎重力基點為改算點(4 219 m)的布格重力異常值與高程回歸分析散點圖,二者為負相關(guān),斜率-0.003,相關(guān)系數(shù)為-0.14,截距為-505.630;圖5c是以研究區(qū)內(nèi)乃東縣重力基點為改算點(3 522 m)的布格重力異常值與高程相關(guān)分析散點圖,斜率-0.001,相關(guān)系數(shù)為-0.03,截距為-558.055,說明高程與重力場值不相關(guān),布格重力異常值不受地形起伏的影響,說明重力場值與近地表地層巖石密度不均勻性關(guān)系密切。
通過以上線性回歸分析及相關(guān)分析,結(jié)合研究區(qū)內(nèi)地層、巖體的出露情況及研究區(qū)內(nèi)采測的地層巖(礦)石、鉆孔巖芯密度統(tǒng)計結(jié)果,結(jié)合地形地貌特征,形成以下結(jié)論:①布格重力異常與高程相關(guān)性分析,其相關(guān)系數(shù)、斜率、截距的大小,與進行統(tǒng)計現(xiàn)象變量的數(shù)量的多少及現(xiàn)象變量選取區(qū)段有關(guān)。只改變現(xiàn)象變量的數(shù)量和現(xiàn)象變量中H的高度段,得到的結(jié)果相差很大,布格重力異常與高程為呈顯著線性相關(guān),進行比較,可以看出異常范圍、形態(tài)大體一致,說明布格重力異常是近地表地層巖石密度的不均勻性引起的[16-17]。②從宏觀考慮,依據(jù)重力場特征,重力場分區(qū)、斷裂構(gòu)造劃分、巖體圈定等使用的布格重力異常數(shù)據(jù),為《區(qū)域重力調(diào)查規(guī)范》(DZ/T 0082—2021)中的計算公式計算的布格重力異常值,中間層改正密度采用2.67 g/cm3,繪制重力基礎(chǔ)圖件,及重力位場轉(zhuǎn)換圖件;③針對異常區(qū)重力剖面查證和小范圍內(nèi)的高精度重力調(diào)查,結(jié)合研究區(qū)內(nèi)地層、巖礦石密度資料,采用《大比例尺重力勘查規(guī)范》(DZ/T 0171—1997)中的計算公式計算布格重力異常值,中間層改正密度采用2.80 g/cm3,以羅布莎重力基點(4 219 m)為改算點(見圖6),繪制相應(yīng)的基礎(chǔ)圖件和位場轉(zhuǎn)換圖件。
多元統(tǒng)計分析在高海拔地區(qū)、大區(qū)域、大比例尺地球物理異常解釋中具有必不可少的關(guān)鍵作用。
布格重力異常在高原山區(qū)和平原區(qū)通常為負值,海洋區(qū)為正值,即布格重力異常和地形高程呈負相關(guān)。本研究區(qū)南北橫跨祁連山區(qū),地形起伏大,因此以構(gòu)造特征為主,結(jié)合平均海拔高程及地形起伏程度,將研究區(qū)分為四個區(qū)域分別進行布格重力異常與高程關(guān)系的統(tǒng)計分析(見圖7)。
圖7 布格重力異常與高程相關(guān)分析
通過相關(guān)性分析統(tǒng)計研究,得出以下3個結(jié)論:①布格重力異常與高程為負相關(guān),由北向南,布格重力異常與高程負相關(guān)系數(shù)分別為-0.655、-0.690、-0.534、-0.481,其絕對值不斷減小,說明布格重力異常與高程相關(guān)程度在研究區(qū)北部較南部緊密,南部布格重力異??赡芘c應(yīng)力或巖石圈異常密度分布之間的關(guān)系更加密切;②自空重力異常與高程為正相關(guān),各區(qū)域相關(guān)系數(shù)都大于0.5,為緊密相關(guān);③各區(qū)域曲線斜率不同,反映各區(qū)塊地殼內(nèi)部結(jié)構(gòu)不同,物質(zhì)密度分布不均勻[18]。
地球化學找礦現(xiàn)已發(fā)展成為尋找殘坡積層下隱伏礦床的一種有效的勘查地球化學方法。在地球化學測量實際工作中,因測試元素較多造成測試數(shù)據(jù)繁雜龐大,對土壤地球化學異常特征總結(jié)及解釋評價帶來諸多不便,從而影響找礦效果。因此,利用數(shù)學地質(zhì)中的數(shù)理統(tǒng)計基本原理,對土壤樣品測試的數(shù)據(jù)進行多元統(tǒng)計分析,對相關(guān)元素組合分類,探索影響地球化學異常特征的主要因素,可有效排除地球化學異常中的干擾因素,得到成礦指示元素組合特征,再結(jié)合研究區(qū)地質(zhì)特征進行篩選,最終提取最佳元素組合異常[19]。
元素組合特征是地球化學特征中特別重要的一個方面。將勘查區(qū)土壤地球化學測量因子分析結(jié)果經(jīng)過方差極大旋轉(zhuǎn)獲得正交旋轉(zhuǎn)因子矩陣,正交旋轉(zhuǎn)因子載荷矩陣可以使各元素在某因子上有較高載荷,而在其他因子上只有較小的載荷,使每個因子的元素組合更加明顯,地質(zhì)意義更加清楚。同時進行R型聚類分析,反映元素的共生組合關(guān)系。本次報告使用SPSS軟件對土壤樣品各元素進行因子分析與聚類分析。
利用SPSS軟件對地球化學元素數(shù)據(jù)進行統(tǒng)計分析,采用主成分提取,使用有Kaiser標準化的正交旋轉(zhuǎn)法,旋轉(zhuǎn)在4次迭代后收斂,采用3個因子正交旋轉(zhuǎn)因子矩陣,根據(jù)研究區(qū)域的元素分布特征與地質(zhì)綜合因素,確定元素組合因子得分低限值絕對值>0.5。元素相關(guān)矩陣見表1,正交旋轉(zhuǎn)因子矩陣見表2。
表1 元素相關(guān)矩陣值向量間的相關(guān)性
表2 正交旋轉(zhuǎn)因子矩陣
由表2可以看出,3個因子基本反映了研究區(qū)元素組合特征及成礦作用特點。因子1:Cu-Ni高溫成礦元素組合和Cu-Zn中溫成礦元素組合,推斷與勘查區(qū)南部出露的玄武巖體有關(guān)。因子2:Ag-Pb中溫成礦元素組合,可能與花崗斑巖脈體侵入有關(guān)。因子3:Au單獨元素組合,應(yīng)與淺成低溫熱液礦化作用有關(guān),且成礦期單一,可能只受到一次主要成礦熱液活動影響。
利用SPSS軟件進行R型聚類分析(系統(tǒng)聚類分析),方法采用組間連接,度量標準采用距離的平方。
從勘查區(qū)元素的R型聚類分析結(jié)果(見圖8)可以看出,當X軸(距離平方)選取16時,可分為3簇。第1簇:Pb-Mo-Ag組合,這3種元素共生可能與脈巖侵入有關(guān);第2簇:Cu-Zn-Ni為一組中高溫組合,其富集主要與玄武巖有關(guān),在有利區(qū)塊具有富集成礦的可能;第3簇:Au與其他元素組相關(guān)性表現(xiàn)較弱。
圖8 R型聚類分析
綜合研究區(qū)土壤測量因子分析和聚類分析結(jié)果,將研究區(qū)7個元素劃分為3個元素組合。①中高溫成礦元素組合:Cu-Zn-Ni;②中溫元素組合:Pb-Mo-Ag;③Au與其他元素親和性較差,呈獨立因子組合。對研究區(qū)所有地球化學元素的含量進行統(tǒng)計分析,計算出平均值、變化系數(shù)、標準離差及富集系數(shù)等各項地球化學指標,了解區(qū)內(nèi)各元素的分布與變化規(guī)律。平均值是一個集中性參數(shù),反應(yīng)勘查區(qū)各元素含量的總體水平;標準離差與變化系數(shù)均為離散型參數(shù),反應(yīng)樣品含量的離散程度,但變化系數(shù)為相對參數(shù),可用于不同元素之間對比;富集系數(shù)為平均值與地殼元素豐度的比值,能反應(yīng)元素在區(qū)域上的分散、富集規(guī)律。
多元統(tǒng)計分析方法有很多,在各個學科都得以廣泛應(yīng)用,尤其在地球科學中也有著重要的應(yīng)用價值。通過前面的研究主要形成以下幾個方面的結(jié)論。
1)地球物理勘探工作中研究區(qū)及外圍地形對布格重力異常的影響不容忽視,地形的隆起或者凹陷,都將引起布格重力異常曲線圖在形態(tài)上、走向上、異常圈閉直接變化,影響到研究者對異常從客觀上的解釋,線性相關(guān)分析方法有著不可取代的積極作用。
2)巖石地球物理數(shù)據(jù)統(tǒng)計分析過程中,統(tǒng)計分析、相關(guān)分析與正態(tài)分布有很不錯的效果,若結(jié)合專業(yè)分析處理軟件,將能夠?qū)r礦石地球物參數(shù)在橫向上、縱向上都有很明顯效果,分析結(jié)果對地球物理異常解釋至關(guān)重要。
3)重力資料數(shù)據(jù)處理中,相關(guān)分析、回歸分析在布格重力異常與高程的關(guān)系分析、中間層改正的密度選取都是一個必經(jīng)過程,分析結(jié)果對重力異常最終數(shù)據(jù)改正值的選取,異常結(jié)果是否合理,都具很重要的作用,尤其在青藏高原的大比例尺,其余地區(qū)的小比例尺重力工作中,應(yīng)用相關(guān)分析、回歸分析選取中間層密度,以及高程改正起點的選取都將決定解釋的精度和準確性。
4)在地球化學勘查中運用因子分析、聚類分析、統(tǒng)計分析、相關(guān)元素分析等統(tǒng)計方法,輔助SPSS統(tǒng)計軟件,結(jié)合勘查異常區(qū)地質(zhì)及構(gòu)造特征,可以獲取礦化特征的元素共生組合因子,依據(jù)異常特征及組合因子研究成礦特征和預(yù)測成礦規(guī)律。
5)實際工作中依據(jù)研究區(qū)地球化學數(shù)據(jù)相關(guān)分析、聚類分析和因子分析成果,揭示研究區(qū)地質(zhì)異常中元素共生組合規(guī)律,這將是判斷礦體剝蝕程度及成礦重要指示元素。利用這些組合異常特征,進而預(yù)測成礦特征和規(guī)律。