張國(guó)利,趙更新,滕菲,王得啟,匡海陽(yáng),黃忠峰,陳楓
(1.中國(guó)地質(zhì)調(diào)查局天津地質(zhì)調(diào)查中心,天津,300170;2.河北地球物理勘查院,廊坊,065000)
高精度重力測(cè)量在冀東鐵礦整裝勘查區(qū)查找隱伏鐵礦中的應(yīng)用
張國(guó)利1,趙更新1,滕菲1,王得啟2,匡海陽(yáng)1,黃忠峰1,陳楓1
(1.中國(guó)地質(zhì)調(diào)查局天津地質(zhì)調(diào)查中心,天津,300170;2.河北地球物理勘查院,廊坊,065000)
在研究冀東鐵礦外圍覆蓋區(qū)的地質(zhì)概況和地球物理特征,野外高精度重力調(diào)查的基礎(chǔ)上,結(jié)合重磁剖面,利用二度半重磁聯(lián)合反演,在重磁異常經(jīng)鉆孔驗(yàn)證見(jiàn)鐵礦厚度20 m左右,認(rèn)為使用剩余重力異常是尋找深部隱伏鐵礦的有效方法,冀東鐵礦外圍具有很大找礦潛力。
冀東鐵礦;剩余重力異常;重磁剖面;重磁異常;
鐵礦勘查以航磁為基礎(chǔ),通過(guò)地面檢查和鉆探驗(yàn)證,然后轉(zhuǎn)入普詳查和勘探,此方法發(fā)現(xiàn)了我國(guó)80%鐵礦[1]。因此,大部分人認(rèn)為磁法可解決鐵礦找礦問(wèn)題,并提出了航(地)磁測(cè)量的鐵礦異常標(biāo)志是高強(qiáng)度、大梯度、尖峰異常。但實(shí)際上也有一些鐵礦找礦問(wèn)題沒(méi)有得到有效解決。實(shí)踐表明,一些埋藏較深、產(chǎn)狀復(fù)雜的規(guī)模鐵礦,單靠磁法找礦難度加大,出現(xiàn)了一些難以解決的新問(wèn)題。這是因?yàn)榇女惓榇排紭O子異常,異常本來(lái)就復(fù)雜,空間位置難確定。同時(shí)會(huì)受剩磁影響?;瘶O后也定位不準(zhǔn),不可靠;南北成帶分布時(shí),獨(dú)立磁性體的磁異常可能發(fā)生正負(fù)異常疊加,造成異常改變、弱化或消失,對(duì)研究成礦規(guī)律不太有利。這是制約隱伏鐵礦找礦的重要技術(shù)問(wèn)題。
地球重力場(chǎng)是地球系統(tǒng)物質(zhì)屬性產(chǎn)生的一個(gè)最基本的物理場(chǎng),反映由地球各圈層相互作用和動(dòng)力過(guò)程決定的物質(zhì)空間分布、運(yùn)動(dòng)和變化[2],它是重力勘探的理論基礎(chǔ)[3]。重力勘探是測(cè)量與圍巖有密度差異的地質(zhì)體在其周?chē)鸬闹亓Ξ惓R源_定這些地質(zhì)體存在的空間位置大小和形狀從而對(duì)工作地區(qū)的地質(zhì)構(gòu)造做出判斷的一種地球物理勘探方法。重力勘探應(yīng)用范圍廣泛,它既適用于如大地構(gòu)造分區(qū)、地殼與上地幔深部構(gòu)造、區(qū)域地質(zhì)調(diào)查等大尺度勘探,也適合探測(cè)古墓[4]、礦山采空區(qū)[5]、溶洞等小尺度目標(biāo)體[6,7]。自從20世紀(jì)30年代重力被用于石油勘探以來(lái),經(jīng)過(guò)幾十年的發(fā)展,已經(jīng)取得巨大成就。我國(guó)1936年首次利用重力進(jìn)行礦產(chǎn)勘探,之后的很長(zhǎng)一段時(shí)間,重力勘探的主要應(yīng)用領(lǐng)域是油氣和固體礦產(chǎn)勘查及區(qū)域重力調(diào)查[8~10]?,F(xiàn)代數(shù)理理論和計(jì)算機(jī)技術(shù)的發(fā)展,完善了各項(xiàng)改正方法。20世紀(jì)90年代GPS測(cè)地技術(shù)的應(yīng)用使重力勘探取得了突破性的進(jìn)展,利用微重力探測(cè)小尺度的目標(biāo)體成為可能,拓寬了重力勘探的應(yīng)用領(lǐng)域。重力測(cè)量方法已成為解決地質(zhì)問(wèn)題最好的方法之一。
因此,針對(duì)傳統(tǒng)的方法(磁法)的局限性和問(wèn)題,結(jié)合重力測(cè)量方法的優(yōu)越性,采用重磁結(jié)合,充分利用重力,結(jié)合電法可能是鐵礦深部勘查研究的必由之路。
航磁測(cè)量時(shí)發(fā)現(xiàn)冀東鐵礦外圍灤縣、灤南為低緩磁異常區(qū),工區(qū)以往重力工作程度低、工作比例尺都小于1:20萬(wàn),且重力資料年代久遠(yuǎn)、精度都較低不利于利用。針對(duì)這些問(wèn)題,2010~2012年在第四系覆蓋區(qū)的低緩航磁異常布置了高精度1:5萬(wàn)航磁調(diào)查,面積約1800 km2。通過(guò)布格重力異常特征分析了區(qū)域地質(zhì)構(gòu)造基本格架,對(duì)重力資料的系統(tǒng)研究,推斷昌黎斷裂以北為太古宙隆起區(qū),主要由賦存含鐵建造
的太古宙沉積變質(zhì)基底(遷西群?)引起,其內(nèi)有大量的南北向展布的隱伏含鐵建造。其他地勘單位在異常區(qū)鉆探驗(yàn)證,在700 m深處見(jiàn)20 m厚度的鐵礦[2]。
調(diào)查區(qū)主要涉及的行政縣是遵化、遷安、灤縣、灤南等。該區(qū)基底為太古宙變質(zhì)巖系,其上為華北地臺(tái)型的沉積蓋層,包括中上元古界、古生界、中生界、新生界地層。與鐵礦有關(guān)的地層主要是太古宇,包括3個(gè)群。1)遷西群:主要分布于遷安水廠-平林鎮(zhèn)-婁子山一帶,包括上川、三屯營(yíng)組,為一套麻粒巖相含鐵礦層的表殼巖系。巖石組合主要為暗色中?;月榱r和含紫蘇輝石黑云斜長(zhǎng)變粒巖,夾透輝斜長(zhǎng)角閃巖及含輝石磁鐵石英巖,變質(zhì)程度為麻粒巖相。遷安的水廠、孟家溝等大型鐵礦均與該層位有關(guān)。2)遵化群:出露于遵化-寬城-青龍一帶,包括拉馬溝組和跑馬廠組,巖石組合主要為斜長(zhǎng)角閃巖、黑云角閃斜長(zhǎng)片麻巖和斜長(zhǎng)變粒巖,夾磁鐵石英巖,變質(zhì)程度達(dá)高角閃巖相。寬城豆子溝鐵礦等產(chǎn)于此層位中。3)灤縣群:主要分布于盧龍及灤縣一帶,指原單塔子群白廟組、八道河群三門(mén)店組、灤縣群及雙山子群中的部分表殼巖。主要為黑云斜長(zhǎng)變粒巖,夾斜長(zhǎng)角閃巖和磁鐵石英巖;中上部夾有片巖,局部夾有大理巖透鏡體;上部魯杖子巖組主要為薄層狀斜長(zhǎng)角閃巖夾各種片巖及(黑云)角閃斜長(zhǎng)變粒巖等。該巖群變質(zhì)程度為低角閃巖相-高綠片巖相。司家營(yíng)、馬城等大中型鐵礦皆產(chǎn)于此層位內(nèi)。
2.1 區(qū)域布格重力異常特征
依據(jù)1:20萬(wàn)布格重力異常圖(圖1),冀東地區(qū)主要處在北東東向重力高值帶上,是我國(guó)東部地區(qū)少數(shù)幾個(gè)重力正異常區(qū)之一,也是地殼減薄區(qū),局部的區(qū)域重力異常主要反映了基底隆起的分布。
調(diào)查區(qū)重力異常呈北東向分布,在重力高和重力低異??缍葞?,重力值由北向南呈降低趨勢(shì),在調(diào)查區(qū)東南部又逐漸升高。
2.2 區(qū)域磁場(chǎng)特征
1:20萬(wàn)航磁△T等值線平面圖(圖2)上顯示,調(diào)查區(qū)航磁異常走向?yàn)榻鼥|西向,北部處于正負(fù)磁異常的伴生負(fù)磁場(chǎng)中,航磁△T最低值約為-200 nT左右。冀東鐵礦則處在正、負(fù)磁異常的過(guò)渡帶上或正磁異常的西側(cè)邊部,推斷正磁異常是司家營(yíng)式沉積變質(zhì)型鐵礦的反映,負(fù)磁場(chǎng)是礦致異常的伴生異常。
調(diào)查區(qū)中部處于正負(fù)磁異常的伴生正磁場(chǎng)中,為東西走向,航磁△T最高值為140 nT左右,大部分為80 nT左右的低緩磁異常,面積約800 km2。調(diào)查區(qū)南部處于正負(fù)磁異常的伴生正磁場(chǎng)中,為東西走向,最高值為100 nT左右,大部分為40 nT左右的低緩磁異常,面積較大。
2.3 巖(礦)石物性特征
本區(qū)鐵礦石與圍巖之間有明顯的磁性和密度差異。圍巖只具微弱磁性、低密度,礦石為強(qiáng)磁性、高密度。圍巖黑云變粒巖和花崗質(zhì)混合巖的磁化率與剩磁基本為零,磁化率相差30 000—90 000單位,剩磁相差5 000—10 000單位。圍巖密度也較低,花崗質(zhì)混合巖密度平均值為2.66×103kg/m3,片麻狀混合巖密度平均值為2.69×103kg/m3,黑云變粒巖密度平均值為2.72×103kg/m3,與礦石之間密度差約為0.8× 103kg/m3左右,差異明顯。當(dāng)全鐵品位在25%~40%之間時(shí),磁性與含鐵品位呈正變。但含鐵品位大于40%時(shí),其磁性則略有下降,且無(wú)規(guī)律。另外,鐵礦石的磁性還與氧化程度密切相關(guān),氧化程度越深,磁性越弱。主要巖礦石的磁性和密度參數(shù)列于表1。
3.1 重力儀器準(zhǔn)備及性能試驗(yàn)
野外工作開(kāi)始之前對(duì)儀器進(jìn)行了充分調(diào)試和性能試驗(yàn),重力儀器包括格值標(biāo)定、靜態(tài)試驗(yàn)、動(dòng)態(tài)試驗(yàn)及一致性試驗(yàn),測(cè)地儀器天寶雙頻5800型GPS儀器進(jìn)行一次全面檢查和性X能鑒定,儀器合格的方可投入生產(chǎn)。
3.2 重力基點(diǎn)網(wǎng)、GPS測(cè)量控制網(wǎng)的聯(lián)建
根據(jù)設(shè)計(jì)方案,結(jié)合地形、交通條件和高程等情況進(jìn)行了布點(diǎn)。支、基點(diǎn)選擇在地基穩(wěn)固、聯(lián)測(cè)方便、重力水平梯度小、近期內(nèi)不被占用的地方,避免設(shè)在較大的陡崖、河堤、大樹(shù)底下以及有震源之處。
3.3 野外測(cè)點(diǎn)重力值觀測(cè)及坐標(biāo)觀測(cè)
測(cè)點(diǎn)點(diǎn)位布設(shè)有關(guān)技術(shù)要求執(zhí)行。野外作業(yè)中采用地形圖定點(diǎn)與手持導(dǎo)航儀相結(jié)合的方法布設(shè)測(cè)點(diǎn)。定點(diǎn)是以1:5萬(wàn)地形圖為工作底圖布設(shè)測(cè)點(diǎn)。為減小地形影響,測(cè)點(diǎn)盡量布設(shè)在25 m范圍內(nèi)平坦的地帶。
測(cè)點(diǎn)重力觀測(cè)采用單程觀測(cè)法,起閉于基點(diǎn)上。首先在基點(diǎn)上進(jìn)行基—輔—基觀測(cè),基點(diǎn)三次讀數(shù),輔基二次讀數(shù),最大與最小讀數(shù)之差小于0.005格;然后進(jìn)行測(cè)點(diǎn)觀測(cè),每個(gè)測(cè)點(diǎn)讀二次數(shù),兩次讀數(shù)之間差值小于0.005格。平均數(shù)采用四舍五入法記錄,最后閉合于重力基點(diǎn)。要求重力測(cè)點(diǎn)觀測(cè)的每個(gè)閉合段的零點(diǎn)位移值不大于重力測(cè)點(diǎn)觀測(cè)精度的2-3倍。
3.4 質(zhì)量檢查
為了衡量野外觀測(cè)結(jié)果質(zhì)量,認(rèn)真完成了檢查觀測(cè)工作。質(zhì)量檢查點(diǎn)隨外業(yè)工作同步進(jìn)行,力爭(zhēng)時(shí)間和空間上均勻分布;近區(qū)地改也要實(shí)行加密方位復(fù)測(cè)的方法觀測(cè)和評(píng)價(jià)近區(qū)地改的精度。質(zhì)量檢查工作采用“一同三不同”方法,即同一點(diǎn)位、不同時(shí)間、不同儀器、不同操作員的檢查;近區(qū)改正采用一同三不同的原則進(jìn)行野外實(shí)測(cè)進(jìn)行。中區(qū)地形改正采用不同人進(jìn)行重新讀圖和重新計(jì)算檢查。遠(yuǎn)區(qū)地改精度采用區(qū)域數(shù)據(jù)庫(kù)的精度評(píng)價(jià)質(zhì)量,檢查工作量占總工作量的3%,質(zhì)量檢查點(diǎn)布設(shè)在每個(gè)工作單元并均勻分布。
表1 冀東地區(qū)巖礦石物性參數(shù)表Table 1 Physical parameters of rocks and minerals table Jidong area
3.5 資料整理
數(shù)據(jù)整理嚴(yán)格按照重力勘查規(guī)范[11]的要求,統(tǒng)一采用2000國(guó)家重力基本網(wǎng)系統(tǒng),統(tǒng)一采用1954年北京坐標(biāo)系和1985國(guó)家高程基準(zhǔn),統(tǒng)一采用國(guó)際大地測(cè)量協(xié)會(huì)(IAG)推薦的1980年公式計(jì)算正常重力值,統(tǒng)一采用規(guī)范規(guī)定的公式進(jìn)行布格改正和中間層改正,密度統(tǒng)一采用2.67 g/cm3,統(tǒng)一采用166.7 km的半徑進(jìn)行地形改正的“五統(tǒng)一”要求執(zhí)行。最后統(tǒng)計(jì)重力觀測(cè)精度26微伽,GPS平面誤差0.22 m,高程誤差0.103 m,布格總精度55微伽。
根據(jù)觀測(cè)重力值得到的重力異?;虿几裰亓Ξ惓#藦牡乇淼缴畈克忻芏炔痪鶆蛞鸬闹亓π?yīng),重力異常是所有這些重力效應(yīng)的總和或疊加[3]。數(shù)據(jù)處理的目的是通過(guò)不同的數(shù)據(jù)處理手段,有效地克服或壓制不同干擾異常,從疊加場(chǎng)中分離或突出某些目標(biāo)物的場(chǎng),根據(jù)地球物理模型用數(shù)學(xué)變換等手段將目標(biāo)場(chǎng)提取出來(lái),達(dá)到突出區(qū)域重力場(chǎng)信息、突出與強(qiáng)化斷裂帶異常信息、突出局部重力異常信息,順利達(dá)到完成區(qū)域重力場(chǎng)特征分析,斷裂構(gòu)造劃分與分析,圈閉巖體,分析礦產(chǎn)有利區(qū)域等地質(zhì)任務(wù)。
根據(jù)本區(qū)內(nèi)斷裂構(gòu)造發(fā)育,巖漿巖分布廣泛,我們選用如下數(shù)據(jù)處理方法:解析延拓(向上延拓)以壓制淺部干擾反映深部信息,并定性推斷斷裂及巖體深部的賦存狀況,在上延的同時(shí)進(jìn)行不同方向的水平導(dǎo)數(shù)計(jì)算和垂向二次導(dǎo)數(shù)的計(jì)算,以反映不同深處斷裂、巖體的賦存狀況。
對(duì)圈定的異常進(jìn)行統(tǒng)一編號(hào)、分類(lèi)、描述,并建立異??ㄆ?,內(nèi)容包括異常編號(hào)、異常位置、異常特征、地質(zhì)背景概況、初步推斷及分類(lèi)。
4.1 數(shù)據(jù)預(yù)處理
為了避免或減少某些數(shù)據(jù)處理方法帶來(lái)的邊緣損失或邊界效應(yīng)、及正確勾繪圖幅邊部的重力異常等值線,在數(shù)據(jù)處理前利用相鄰圖幅1:20區(qū)域重力調(diào)查成果資料,對(duì)各圖幅四周均進(jìn)行一定范圍的擴(kuò)充。
4.2 數(shù)據(jù)網(wǎng)格化
為適應(yīng)計(jì)算機(jī)數(shù)據(jù)處理及繪圖要求,對(duì)擴(kuò)充后的重力數(shù)據(jù)進(jìn)行了網(wǎng)格化,具體工作時(shí)使用Sufer軟件采用“Kring泛克立格法”方法進(jìn)行網(wǎng)格化,初步采用網(wǎng)格距為500 m×250 m2,搜索半徑初步定為1000 m。
4.3 區(qū)域背景場(chǎng)及剩余異常的分離
一般來(lái)說(shuō),重力場(chǎng)三級(jí)異常主要體現(xiàn)地塊內(nèi)部隆起與斷陷的密度分布關(guān)系,它表現(xiàn)的是布格重力異常的局部特征,這些局部特征可以在異常分離后的剩余場(chǎng)中得到突出。為求得重力局部異常,需要通過(guò)數(shù)據(jù)處理求取區(qū)域場(chǎng)和局部剩余異常。1:5萬(wàn)重力異常采用窗口滑動(dòng)平均法求取區(qū)域場(chǎng)。窗口大小通過(guò)實(shí)驗(yàn)來(lái)具體確定。網(wǎng)格化數(shù)據(jù)直接減去區(qū)域場(chǎng)數(shù)據(jù)求取重力剩余異常。
4.4 求導(dǎo)
為了突出不同走向的斷裂、脈巖位置和寬大地質(zhì)體的邊界線,對(duì)布格重力異常數(shù)據(jù)作水平導(dǎo)數(shù)(方向?qū)?shù))處理。求導(dǎo)方向依據(jù)研究區(qū)內(nèi)實(shí)際的和推斷的構(gòu)造方向來(lái)確定,初步按0、45、90、135度四個(gè)方向的求取導(dǎo)數(shù)。
4.5 重力異常向上延拓
重力異常的向上延拓是指由觀測(cè)平面上已知重力異常換算得到較高的另一些平面上的重力異常。向上延拓將使深源或者大尺度密度體和淺源或小尺度密度體引起的異常都衰減,但是淺源或小尺度密度體引起的異常衰減速度比深源或者大尺度密度體引起的異常(區(qū)域異常)快,因此向上延拓一定高度后的異常能夠較為明顯地反映區(qū)域異常的特征。
5.1 布格重力異常特征
在重力布格異常圖(圖3)上明顯看到,測(cè)區(qū)北部基底抬升形成隆起,南部地區(qū)基底下沉形成坳陷。測(cè)區(qū)東北有明顯的重力高,是由昌黎斷裂引起。
5.2 剩余重力異常特征
在高精度剩余重力異常平面等值線圖上(圖4),剩余重力異常在測(cè)區(qū)中部呈北東走向,與已知昌黎斷裂基本一致,測(cè)區(qū)北部剩余重力異常呈近東西向展布,剩余重力異常圖也反映了部分的次級(jí)斷裂,走向基本也是北東向。
根據(jù)收集到的部分地質(zhì)鉆孔資料(河北物理勘察院),進(jìn)行了初步反演,反演結(jié)果推斷重磁異常為朝西緩傾的磁性體引起(圖5),頂部埋深約0~150 m
左右。
利用其它方法給出的條件,如測(cè)井資料,進(jìn)行約束反演,可以收到較好的效果[12]。因此有待進(jìn)一步搜集相關(guān)地質(zhì)、物性資料,及鉆孔資料進(jìn)行驗(yàn)證,然后再進(jìn)行反演,需要不斷改進(jìn)。
從地質(zhì)背景、地球物理特征及控礦地層、巖石特征及構(gòu)造環(huán)境等方面的分析,研究區(qū)有良好的成礦地質(zhì)條件,具有尋找大中型沉積變質(zhì)型鐵礦的巨大潛力。
圖3 布格重力異常影像圖Fig.3 Bouguer gravity anomaly image
圖4 構(gòu)造斷裂初步推斷圖Fig.4 Preliminary structural faults inferred map
圖5 1號(hào)線重磁綜合剖面圖Fig.5 Line 1 sectional view of the gravity and magnetic integrated
結(jié)合地質(zhì)研究,利用剩余重力異常及重磁綜合剖面,進(jìn)行二度半重磁聯(lián)合反演,可以大致確定礦體的位置,使得在冀東鐵礦外圍的找礦工作有所突破,高精度重力勘探配合重磁剖面是尋找隱伏磁鐵礦較為有效的找礦手段。
[1]孫文柯.中國(guó)礦床發(fā)現(xiàn)史.物化探卷[M].北京:地質(zhì)出版社,2002,10.
[2]寧津生,王正濤.地球重力場(chǎng)研究現(xiàn)狀與進(jìn)展[J].高端論壇, 2013,38(1):1-7.
[3]曾華霖.重力場(chǎng)與重力勘探[M],北京:地質(zhì)出版社,2005:1
[4]袁炳強(qiáng),楊明生,劉士毅等.高精度重力測(cè)量探測(cè)秦始皇帝陵地宮[J].地球科學(xué)-中國(guó)地質(zhì)大學(xué)學(xué)報(bào),2005,30(4):498-502.
[5]張善法,孟令順,杜曉娟,等.高精度重力測(cè)量在金礦采空區(qū)探測(cè)中的應(yīng)用研究[J].地球物理學(xué)進(jìn)展,2009,24(2): 590-595
[6]王懋基.中國(guó)重力勘探的新進(jìn)展[J].地球物理學(xué)報(bào),1994(S1):353-360
[7]滕龍,倪四道,張寶松,等.CORS系統(tǒng)在重力調(diào)查中的應(yīng)用[J].物探與化探,2013,37(6):1018-1022.
[8]王懋基,蔡鑫,涂承林.中國(guó)重力勘探的發(fā)展與展望[J].地球物理學(xué)報(bào),1997(S1):292-298.
[9]張赤軍,邊少鋒.關(guān)于油氣勘探開(kāi)發(fā)中的高程及其確定方法[J].物探與化探,2005(04):342-346.
[10]周?chē)?guó)藩.油氣勘查中重力方法技術(shù)的新進(jìn)展[J].物探與化探,1992,16(1):14-16.
[11]孫文珂,丁鵬飛.DZ/T 0171—1997大比例尺重力勘查規(guī)范[S].地質(zhì)礦產(chǎn)部,1997.
[12]劉光鼎,郝天珧.應(yīng)用地球物理方法尋找隱伏礦床[J].地球物理學(xué)報(bào),1995,36(6):850-854.
Application of the High Precision Gravity Survey in Prospeting Concealed Deposit in Jidong Iron Ore Peripheral Coverage Area
ZHANG Guo-li1,ZHAO Geng-xin1,TENG Fei1,WANG De-qi2,KUANG Hai-yang1,HUANG Zhong-feng1,CHEN Feng1
(1.Tianjin Center of China Geological Survey,Tianjin,300170,China; 2.Hebei Institute of Geophysical Exploration,Langfang Hebei,065000,China)
This paper introduced the process of the high precision gravity survey in Jidong iron ore peripheral coverage area.Combined with the gravity and magnetic profile,using second degree semi-gravity and magnetic data joint inversion,the gravity and magnetic anomaly by drilling verification finds iron thickness of 20 meters.This overview summarized that the residual gravity anomaly is effective method in search for deep concealed ore,and Jidong iron ore periphery had great prospecting potential.
Jidong iron ore;residual gravity anomaly;gravity and magnetic profile;gravity and magnetic anomaly;
P631.1
A
1672-4135(2013)04-0046-06
2013-10-25
國(guó)家地質(zhì)大調(diào)查項(xiàng)目:冀東鐵礦外圍1:5萬(wàn)重力調(diào)查(1212010050230)
張國(guó)利(1980-),男,高級(jí)工程師,中國(guó)地質(zhì)調(diào)查局天津地質(zhì)調(diào)查中心,主要從事應(yīng)用地球物理研究,Email: zgl003523@163.com。