郭華， 王明， 姚雨暘， 常暢,4

1 中國(guó)自然資源航空物探遙感中心， 北京 100083 2 吉林大學(xué)， 長(zhǎng)春 130026 3 中國(guó)地質(zhì)大學(xué)(北京)， 北京 100084 4 中國(guó)地質(zhì)科學(xué)院地質(zhì)研究所， 北京 100037

0 引言

作為新興產(chǎn)業(yè)的無(wú)人機(jī)航空物探在基礎(chǔ)地質(zhì)調(diào)查、能源資源調(diào)查、礦產(chǎn)勘查等領(lǐng)域具有廣闊的應(yīng)用前景，受到了許多國(guó)家普遍關(guān)注并相繼開(kāi)始了基于固定翼或直升機(jī)的無(wú)人機(jī)平臺(tái)的航空地球物理勘探相關(guān)的裝備研發(fā)(Lum et al., 2005; Kurvinen et al., 2005; Partner, 2006; 張洪瑞和范正國(guó), 2007; 熊盛青, 2009; P?ll?nen et al., 2009; Barnard, 2008, 2010; Koyama, 2013; Laliberte, 2013; Hitoshi et al., 2013; Stoll, 2013; 于顯利等, 2012; 李文杰等, 2014; 胥值禮等, 2016; 郭華和謝宏，2018；李飛等, 2018；張富明等，2019；鞠星等，2020a，b；謝小國(guó)等，2021).無(wú)人機(jī)航空物探技術(shù)在國(guó)外的發(fā)展以及相關(guān)應(yīng)用相對(duì)國(guó)內(nèi)成熟度較高，已經(jīng)在許多領(lǐng)域開(kāi)展了試驗(yàn)或應(yīng)用(李文杰等, 2014)；近十年國(guó)內(nèi)中國(guó)自然資源航空物探遙感中心、中國(guó)地質(zhì)科學(xué)院及中國(guó)科學(xué)院遙感與數(shù)字地球研究所、中船重工第七一五研究所等多家單位先后開(kāi)展的基于固定翼或直升機(jī)的無(wú)人機(jī)航磁測(cè)量技術(shù)的研究工作尚處于研究階段.2012—2015年，為滿足國(guó)內(nèi)地質(zhì)調(diào)查和礦產(chǎn)資源勘探的迫切需求，中國(guó)地質(zhì)調(diào)查局組織中國(guó)地質(zhì)科學(xué)院地球物理地球化學(xué)勘查研究所、核工業(yè)航測(cè)遙感中心和航天空氣動(dòng)力技術(shù)研究院等多家科研單位共同研制完成了世界上第一套基于“彩虹-3”無(wú)人機(jī)航空地球物理勘探的綜合測(cè)量系統(tǒng)，自此無(wú)人機(jī)開(kāi)始了在地質(zhì)資源調(diào)查領(lǐng)域的民用征程.基于“彩虹-3”無(wú)人機(jī)航空地球物理勘探測(cè)量系統(tǒng)相對(duì)于以往的航空地球物理勘探測(cè)量系統(tǒng)具備較強(qiáng)的環(huán)境適應(yīng)能力、較高的安全保障、優(yōu)秀的夜航巡航能力，目前為止在我國(guó)及贊比亞等地“彩虹-3”無(wú)人機(jī)航空地球物理綜合測(cè)量系統(tǒng)累計(jì)完成16.5萬(wàn)測(cè)線公里的調(diào)查任務(wù)，填補(bǔ)了國(guó)內(nèi)外無(wú)人航空地球物理測(cè)量領(lǐng)域裝備的空白(中國(guó)地質(zhì)科學(xué)院, 2014; 郭華和謝宏, 2018).

“彩虹-4”無(wú)人機(jī)是中國(guó)航天科技集團(tuán)公司下屬的航天彩虹無(wú)人機(jī)股份有限公司在“彩虹-3”無(wú)人機(jī)基礎(chǔ)上研發(fā)的一種中程察打/偵察探測(cè)(民用版)一體無(wú)人機(jī)；其性能上，較“彩虹-3”無(wú)人機(jī)有了很大的提升(http:∥baike.baidu.com/view/9615585.html)，其翼展18.0 m，最大起飛重量1330 kg，最大使用半徑2000 km，輕載100 kg時(shí)航時(shí)約30 h，最大裝載量可達(dá)345 kg，航時(shí)約12 h，是中空長(zhǎng)航時(shí)無(wú)人機(jī)(鞠星等，2020b；崔志強(qiáng)等，2016)，但其作為無(wú)人機(jī)航空物探搭載平臺(tái)在國(guó)內(nèi)外行業(yè)領(lǐng)域內(nèi)屬于空白.在航空磁測(cè)領(lǐng)域，與傳統(tǒng)的航空磁總場(chǎng)測(cè)量相比，矢量測(cè)量所獲得的地磁場(chǎng)模量大小和方向信息能夠有效減少反演的多解性，提高地下磁性異常體的探測(cè)精度(林君等，2017).對(duì)于航空磁矢量測(cè)量，國(guó)內(nèi)外相繼開(kāi)展了相關(guān)的研究(Christensen and Dransfield, 2002；孫昂等，2017；Xie et al., 2020)，而基于“彩虹-4”固定翼無(wú)人機(jī)航磁三分量測(cè)量系統(tǒng)的研發(fā)在國(guó)內(nèi)外也屬于首次.

中國(guó)自然資源航空物探遙感中心在國(guó)家重點(diǎn)研發(fā)項(xiàng)目“航空磁場(chǎng)測(cè)量技術(shù)系統(tǒng)研制”支持下，與吉林大學(xué)、中國(guó)航天空氣動(dòng)力技術(shù)研究院以及中國(guó)地質(zhì)科學(xué)院物化探研究所等多家科研院所合作首次開(kāi)展了“彩虹-4”無(wú)人機(jī)航磁三分量測(cè)量系統(tǒng)研發(fā)工作.通過(guò)開(kāi)展“彩虹-4”無(wú)人機(jī)平臺(tái)改裝與系統(tǒng)集成、無(wú)人機(jī)航磁三分量磁補(bǔ)償?shù)群诵暮完P(guān)鍵技術(shù)研究，以及飛行測(cè)試試驗(yàn)(飛行平臺(tái)靜態(tài)干擾測(cè)試、飛機(jī)磁本底測(cè)試、高空八邊形磁干擾補(bǔ)償飛行、低空重復(fù)線飛行以及實(shí)際測(cè)區(qū)飛行)等工作，成功地獲取了高質(zhì)量的“彩虹-4”無(wú)人機(jī)磁本底數(shù)據(jù)和無(wú)人機(jī)三分量磁測(cè)數(shù)據(jù).2019年在國(guó)內(nèi)西北某地圓滿完成試驗(yàn)飛行，“彩虹-4”無(wú)人機(jī)航磁三分量測(cè)量系統(tǒng)各項(xiàng)指標(biāo)滿足《航空磁測(cè)技術(shù)規(guī)范》(DZ/T 0142-2010)(中華人民共和國(guó)國(guó)土資源部, 2010)的要求，測(cè)量精度達(dá)到25 nT的長(zhǎng)航時(shí)、大航程、超視距、高精度飛控設(shè)計(jì)指標(biāo)，填補(bǔ)了國(guó)內(nèi)外“彩虹-4”無(wú)人機(jī)航磁三分量測(cè)量裝備的空白，形成了一套高安全、高效率、高質(zhì)量、全天時(shí)、長(zhǎng)航時(shí)、超低空的智能化“彩虹-4”無(wú)人機(jī)航磁三分量測(cè)量方法技術(shù)，為我國(guó)地質(zhì)調(diào)查和礦產(chǎn)資源勘探增添了先進(jìn)、安全、可靠的技術(shù)和裝備，促進(jìn)了我國(guó)基于無(wú)人機(jī)的航空物探勘探產(chǎn)業(yè)的進(jìn)一步發(fā)展，研究成果具有很好的應(yīng)用前景.

1 測(cè)量系統(tǒng)組成

1.1 航磁三分量測(cè)量原理

圖1所示，地磁場(chǎng)7要素T(地磁場(chǎng)總強(qiáng)度)、X(地磁場(chǎng)北向分量)、Y(地磁場(chǎng)東向分量)、Z(地磁場(chǎng)的垂直分量)、H(地磁場(chǎng)水平分量)、I(地磁傾角)、D(地磁偏角)空間關(guān)系為：

圖1 地磁要素示意圖Fig.1 The diagram of seven elements of geomagnetic field

根據(jù)坐標(biāo)系不同，可以將7個(gè)分量劃分為：直角坐標(biāo)系中有X、Y、Z，球坐標(biāo)系中有H、I、D，柱坐標(biāo)系中有Z、H、D；測(cè)得其中一組就可以求出其他分量(王一，2015).

1.2 系統(tǒng)組成

“彩虹-4”無(wú)人機(jī)航磁三分量測(cè)量系統(tǒng)主要由“彩虹-4”無(wú)人機(jī)平臺(tái)、三分量磁力儀、光泵磁力儀、高精度慣導(dǎo)系統(tǒng)、高度計(jì)、數(shù)據(jù)采集和收錄系統(tǒng)組成，如圖2所示.

圖2 “彩虹-4”無(wú)人機(jī)航磁三分量測(cè)量系統(tǒng)組成示意圖Fig.2 The diagram of aeromagnetic three-component detction system for CH-4 UAV

1.3 三分量磁力儀

三分量磁測(cè)系統(tǒng)所用三軸型號(hào)是英國(guó)Bartington公司生產(chǎn)的Mag-03MS100，性能指標(biāo)詳見(jiàn)表1；慣導(dǎo)系統(tǒng)采用的是加拿大NovAtel公司生產(chǎn)的SPAN-IGM-S1，性能見(jiàn)表2.

表1 三軸磁力儀性能指標(biāo)Table 1 Performance index of triaxial magnetometer

表2 慣導(dǎo)系統(tǒng)性能指標(biāo)Table 2 Performance index of inertial navigation system

2 系統(tǒng)試驗(yàn)測(cè)試

2.1 飛機(jī)磁本底測(cè)試

無(wú)人機(jī)的機(jī)體是由鐵磁性材料和非鐵磁性材料構(gòu)成的，鐵磁性材料是機(jī)體上的主要磁干擾源.首先，為了最大限度消除機(jī)體上的磁干擾，我們采用鋁合金、碳纖維復(fù)合材料等無(wú)磁性材料對(duì)無(wú)人機(jī)進(jìn)行改裝；其次，測(cè)試無(wú)人機(jī)的磁本底數(shù)據(jù)，以便確定三分量磁測(cè)系統(tǒng)安裝位置.根據(jù)機(jī)體的長(zhǎng)度和寬度設(shè)計(jì)的測(cè)量區(qū)域?yàn)?8 m×12 m，網(wǎng)格化間隔為1 m，共216個(gè)測(cè)點(diǎn)(高全明, 2020).

無(wú)人機(jī)機(jī)體磁本底測(cè)量共分三步：第一步在沒(méi)有無(wú)人機(jī)停放的情況下，利用光泵磁力儀測(cè)量216個(gè)測(cè)點(diǎn)的地磁場(chǎng)強(qiáng)度和對(duì)應(yīng)點(diǎn)的磁日變數(shù)據(jù)；第二步在無(wú)人機(jī)停放的情況下，重復(fù)第一步的操作；第三步利用第二步測(cè)量的磁場(chǎng)數(shù)據(jù)減去第一步測(cè)量的磁場(chǎng)數(shù)據(jù)和磁日變數(shù)據(jù)，即為無(wú)人機(jī)機(jī)體磁本底數(shù)據(jù)(圖3).

圖3 “彩虹-4”無(wú)人機(jī)機(jī)體左半部分磁本底測(cè)試結(jié)果Fig.3 Magnetic background test results of left half for CH-4 UAV

2.2 飛機(jī)發(fā)動(dòng)干擾測(cè)試

測(cè)試無(wú)人機(jī)在飛行過(guò)程，發(fā)動(dòng)機(jī)對(duì)三軸磁力儀產(chǎn)生的磁干擾影響.主要測(cè)試和分析無(wú)人機(jī)的發(fā)動(dòng)機(jī)處于“停止”、“啟動(dòng)”和“臨飛”三種狀態(tài)下產(chǎn)生的磁干擾情況 (高全明, 2020)，如圖4所示.

圖4 發(fā)動(dòng)機(jī)在不同工作狀態(tài)下的磁干擾曲線Fig.4 The curve of magnetic interference for CH-4 engine under different working conditions

從圖4中可以看到，當(dāng)無(wú)人機(jī)的發(fā)動(dòng)機(jī)啟動(dòng)過(guò)程中，磁場(chǎng)曲線上會(huì)出現(xiàn)一個(gè)明顯地“階躍響應(yīng)”變化，這個(gè)變化隨著無(wú)人機(jī)的發(fā)動(dòng)機(jī)轉(zhuǎn)速的提升而不斷增大，直到發(fā)動(dòng)機(jī)轉(zhuǎn)速達(dá)到“臨飛”狀態(tài)后開(kāi)始趨于平穩(wěn)，整個(gè)過(guò)程測(cè)量的數(shù)值，變化幅度約為60 nT.因此，當(dāng)無(wú)人機(jī)的發(fā)動(dòng)機(jī)轉(zhuǎn)速保持不變或勻速飛行狀態(tài)時(shí)，可視為其產(chǎn)生的磁干擾源是穩(wěn)定.

2.3 重復(fù)性飛行測(cè)試

重復(fù)線飛行測(cè)試的目的是為了檢驗(yàn)“彩虹-4”無(wú)人機(jī)航磁三分量磁測(cè)系統(tǒng)的動(dòng)態(tài)性能 (高全明, 2020)，進(jìn)而保證磁測(cè)數(shù)據(jù)的質(zhì)量；重復(fù)線飛行測(cè)試是以測(cè)線和重復(fù)線測(cè)量的磁測(cè)數(shù)據(jù)曲線的重合度為量化評(píng)估指標(biāo).

圖5所示為兩條重復(fù)測(cè)線飛行軌跡，長(zhǎng)度約30 km，偏移距離不超過(guò)20 m，平均飛行高度分別為612.03 m和605.18 m，高度偏差小于7 m.

圖5 重復(fù)線試驗(yàn)飛行航跡圖Fig.5 The diagram of repeated line test flight

圖6所示為重復(fù)線試驗(yàn)獲取的航磁三分量數(shù)據(jù)對(duì)比結(jié)果.從圖中可以直觀地看出兩條重復(fù)線上對(duì)應(yīng)點(diǎn)的x、y和z分量的磁測(cè)數(shù)據(jù)曲線變化趨勢(shì)有較好的重合度，內(nèi)符合精度ε分別為5.43 nT、11.61 nT和11.19 nT；試驗(yàn)結(jié)果驗(yàn)證了基于“彩虹-4”無(wú)人機(jī)的航磁三分量測(cè)量系統(tǒng)具有很好的一致性和穩(wěn)定性.

2.4 磁補(bǔ)償飛行試驗(yàn)

磁補(bǔ)償飛行試驗(yàn)是將未進(jìn)行磁補(bǔ)償處理的航磁x、y、z三分量數(shù)據(jù)由機(jī)體坐標(biāo)系轉(zhuǎn)換到空間地理坐標(biāo)系下 (高全明, 2020)，如圖7所示.圖中顯示出，坐標(biāo)轉(zhuǎn)換后三分量數(shù)據(jù)變化幅度很明顯，主要是由機(jī)體產(chǎn)生的磁干擾所導(dǎo)致的，特別是z分量數(shù)據(jù)變化幅度尤為明顯，這說(shuō)明了進(jìn)行磁補(bǔ)償是十分必要的.

根據(jù)《航空磁測(cè)技術(shù)規(guī)范》(DZ/T 0142-2010)(中華人民共和國(guó)國(guó)土資源部, 2010)，磁補(bǔ)償飛行試驗(yàn)在敦煌市機(jī)場(chǎng)附近磁場(chǎng)變化穩(wěn)定的區(qū)域進(jìn)行，磁場(chǎng)梯度值不超過(guò)60 pT/m，磁環(huán)境滿足磁補(bǔ)償飛行試驗(yàn)要求.在磁補(bǔ)償飛行試驗(yàn)中“彩虹-4”無(wú)人機(jī)按照順時(shí)針?lè)较蜓刂?°→90°→180°→270°和45°→135°→225°→315°的矩形閉合框飛行(圖8)，邊框的航線的長(zhǎng)度約為16 km，平均飛行高度為3974.46 m，高度變化小于15 m；在每一邊上依次完成航向、橫滾和俯仰三個(gè)補(bǔ)償飛行動(dòng)作，運(yùn)動(dòng)幅度變化為±5°(圖9)，用時(shí)約40 min.如圖10所示，磁補(bǔ)償飛行時(shí)慣導(dǎo)系統(tǒng)記錄姿態(tài)角度的變化情況.

圖6 航磁三分量數(shù)據(jù)重復(fù)性對(duì)比圖Fig.6 The comparison diagram of repeated lines for aeromagnetic three-component

圖7 坐標(biāo)轉(zhuǎn)換前后航磁三分量磁測(cè)數(shù)據(jù)對(duì)比圖Fig.7 The comparison diagram of before and after coordinate transformation for aeromagnetic three-component

同時(shí)，為了評(píng)價(jià)航磁三分量磁補(bǔ)償方法的可靠性，文中采用了布谷鳥(niǎo)搜索算法進(jìn)行計(jì)算，并分別與基于無(wú)跡卡爾曼濾波算法(UKF)和遺傳算法(GA)的磁補(bǔ)償方法進(jìn)行對(duì)比分析 (高全明, 2020)，結(jié)果如圖11所示.按照《航空磁測(cè)技術(shù)規(guī)范》(DZ/T 0142-2010)(中華人民共和國(guó)國(guó)土資源部, 2010)的要求，應(yīng)用CS補(bǔ)償算法的航磁x、y、z三分量的動(dòng)態(tài)四階差分噪聲水平均小于25 nT，并且從異常曲線可以看出，補(bǔ)償后的曲線更平滑，補(bǔ)償效果更加明顯.

圖8 磁補(bǔ)償飛行軌跡示意圖Fig.8 The diagram of magnetic compensation flight

圖9 磁補(bǔ)償飛行時(shí)飛機(jī)動(dòng)作示意圖Fig.9 The diagram of aircraft action during magnetic compensation flight

3 實(shí)際飛行試驗(yàn)

在敦煌機(jī)場(chǎng)附近開(kāi)展區(qū)域磁測(cè)飛行試驗(yàn)，測(cè)線軌跡如圖12所示，共飛行了6條測(cè)線，每條測(cè)線長(zhǎng)度約為10.7 km，線距為520 m，飛行高度約450 m，高度偏差不超過(guò)20 m，滿足《航空磁測(cè)技術(shù)規(guī)范》(DZ/T 0142-2010)(中華人民共和國(guó)國(guó)土資源部, 2010)規(guī)范要求.

圖10 “彩虹-4”無(wú)人機(jī)磁補(bǔ)償飛行姿態(tài)角度變化Fig.10 The diagram of magnetic compensation flight attitude angle change for CH-4 UAV

圖11 磁干擾補(bǔ)償結(jié)果Fig.11 The diagram of magnetic interference compensation results

圖12 實(shí)際飛行試驗(yàn)測(cè)線軌跡示意圖Fig.12 The diagram of line distribution in actual flight test

圖13 磁補(bǔ)償前航磁三分量合成總場(chǎng)剖面圖Fig.13 The profile of total field calculated by aeromagnetic three-component before magnetic compensation

圖13為未進(jìn)行磁補(bǔ)償處理的航磁三分量磁測(cè)數(shù)據(jù)合成總場(chǎng)剖面圖.可以看出，由每條測(cè)線的航磁三分量數(shù)據(jù)轉(zhuǎn)換的航磁總場(chǎng)曲線變化幅度較大、規(guī)律性和相關(guān)性差.

為了評(píng)估航磁三分量數(shù)據(jù)的精度，我們將由磁補(bǔ)償處理后的航磁三分量數(shù)據(jù)轉(zhuǎn)換的總場(chǎng)數(shù)據(jù)與同機(jī)CS-3銫光泵磁力儀測(cè)得的航磁總場(chǎng)數(shù)據(jù)進(jìn)行對(duì)比分析.圖14顯示了分別經(jīng)過(guò)磁補(bǔ)償處理與日變改正的實(shí)測(cè)總場(chǎng)數(shù)據(jù)與三分量合成總場(chǎng)數(shù)據(jù)，航磁三分量數(shù)據(jù)換算的磁總場(chǎng)變化趨勢(shì)與實(shí)測(cè)磁總場(chǎng)變化整體趨勢(shì)基本一致，異常幅值、異常形態(tài)與實(shí)測(cè)總場(chǎng)吻合；但前者在梯度過(guò)渡帶內(nèi)的磁梯度變化更強(qiáng)烈、主要是由磁干擾補(bǔ)償處理與各項(xiàng)誤差校正不徹底導(dǎo)致的.同時(shí)，數(shù)據(jù)的顯示的梯級(jí)帶由測(cè)線飛行高度差異所引起，因此，需要做進(jìn)一步的調(diào)平等處理.

圖14 航磁總場(chǎng)數(shù)據(jù)對(duì)比(a) 實(shí)測(cè)的航磁總場(chǎng)數(shù)據(jù)； (b) 航磁三分量數(shù)據(jù)換算的總場(chǎng)數(shù)據(jù).Fig.14 The comparison diagram of aeromagnetic total field obtained by different methods(a) Total field measured by optical-pump magnetometer; (b) Total field calculated by aeromagnetic three-component.

圖15 航磁三分量磁測(cè)數(shù)據(jù)等值線圖(a) x分量等值線圖； (b) y分量等值線圖； (c) z分量等值線圖.Fig.15 The contour map of aeromagnetic three-component magnetic(a) The contour map of x component； (b) The contour map of y component； (c) The contour map of z component.

圖15所示為航磁x、y、z三個(gè)分量的磁場(chǎng)等值線圖.由圖15a和圖15b的磁場(chǎng)特征可以看出，沿著北東向存在一個(gè)明顯的航磁梯度變化帶，與圖14中航磁總場(chǎng)反映的梯度變化帶保持著很好對(duì)應(yīng)關(guān)系；圖15c中z分量的磁場(chǎng)特征與圖14中航磁總場(chǎng)特征相比，二者的變化趨勢(shì)一致且高度吻合.因此，航磁三分量數(shù)據(jù)不但可以很好的表征航磁總場(chǎng)特征，還能提供豐富的x、y、z三個(gè)方向的磁場(chǎng)分量信息.基于彩虹-4固定翼飛行平臺(tái)的航磁三分量測(cè)量系統(tǒng)的成功試驗(yàn)飛行，獲取了高精度實(shí)測(cè)三分量數(shù)據(jù)集，為后續(xù)的反演與解釋工作提供豐富的數(shù)據(jù)基礎(chǔ)，有效降低反演的多解性.同時(shí)，依托彩虹-4固定翼無(wú)人機(jī)長(zhǎng)航時(shí)、適應(yīng)夜航、大載重、超低空、適應(yīng)無(wú)磁化改造等特性，此航磁測(cè)量系統(tǒng)的研發(fā)對(duì)于我國(guó)大面積、實(shí)用化的航磁多參量測(cè)量飛行具有十分重要的意義.

4 結(jié)論

通過(guò)技術(shù)攻關(guān)，我國(guó)首次成功地將“彩虹-4”無(wú)人機(jī)改裝集成應(yīng)用到航空物探領(lǐng)域，首次成功地研發(fā)出基于“彩虹-4”無(wú)人機(jī)平臺(tái)的航磁三分量測(cè)量系統(tǒng)，填補(bǔ)了國(guó)內(nèi)外航空地球物理勘探領(lǐng)域內(nèi)裝備的空白.實(shí)際試驗(yàn)飛行驗(yàn)證了“彩虹-4”無(wú)人機(jī)航磁三分量測(cè)量系統(tǒng)各項(xiàng)指標(biāo)滿足《航空磁測(cè)技術(shù)規(guī)范》(DZ/T 0142-2010)(中華人民共和國(guó)國(guó)土資源部, 2010)要求，測(cè)量精度達(dá)到25 nT的長(zhǎng)航時(shí)、大航程、超視距、高精度飛控設(shè)計(jì)指標(biāo)，形成了一套高安全、高效率、高質(zhì)量、全天時(shí)、長(zhǎng)航時(shí)、超低空的智能化“彩虹-4”無(wú)人機(jī)航磁三分量測(cè)量方法技術(shù)，為我國(guó)地質(zhì)調(diào)查和礦產(chǎn)資源勘探增添了先進(jìn)、安全、可靠的技術(shù)和裝備，促進(jìn)了我國(guó)基于無(wú)人機(jī)的航空物探勘探產(chǎn)業(yè)的進(jìn)一步發(fā)展，研究成果具有很好的應(yīng)用前景.

致謝本文研究工作得到了中國(guó)地質(zhì)科學(xué)院物化探研究所劉俊杰高級(jí)工程師、吉林大學(xué)趙靜教授、高全明博士的大力支持，在此表示衷心的感謝.