孫 瀟，陳建國(guó)，2，房曉龍，馬秋石

(1．中國(guó)地質(zhì)大學(xué)(武漢)資源學(xué)院，湖北 武漢 430074;2．中國(guó)地質(zhì)大學(xué)地質(zhì)過(guò)程與礦產(chǎn)資源國(guó)家重點(diǎn)實(shí)驗(yàn)室，湖北 武漢 430074)

云南羊拉銅礦三維地質(zhì)建模及儲(chǔ)量計(jì)算

孫 瀟1，陳建國(guó)1，2，房曉龍1，馬秋石1

(1．中國(guó)地質(zhì)大學(xué)(武漢)資源學(xué)院，湖北 武漢 430074;2．中國(guó)地質(zhì)大學(xué)地質(zhì)過(guò)程與礦產(chǎn)資源國(guó)家重點(diǎn)實(shí)驗(yàn)室，湖北 武漢 430074)

以云南羊拉銅礦的斷層、地層、巖體、礦體為研究對(duì)象，在分析研究建模技術(shù)方案的基礎(chǔ)上，依據(jù)收集到的鉆孔柱狀圖、勘探線剖面圖、中段地質(zhì)圖和地形地質(zhì)圖，結(jié)合野外踏勘與認(rèn)識(shí)，利用Micromine軟件建立了研究區(qū)的地形、地層、斷層以及礦體的三維實(shí)體模型，直觀展示了區(qū)內(nèi)地表、斷層、地層、巖體以及礦體的空間展布以及相互的空間位置關(guān)系。運(yùn)用克里格法對(duì)里農(nóng)礦段的銅礦體進(jìn)行品位估值，獲得了礦體的儲(chǔ)量。

三維地質(zhì)建模;資源量計(jì)算;Micromine軟件;羊拉銅礦;云南德欽

0 引言

三維地質(zhì)建模這一概念是由加拿大Houlding提出的。目前，隨著其理論逐漸成熟和完善以及計(jì)算機(jī)硬件軟件的快速發(fā)展，它已經(jīng)成為研究礦床的又一新途徑。一些礦業(yè)軟件公司推出了各自的產(chǎn)品，國(guó)外較成熟的軟件有：Micromine、Surpac、Datamine、Vulcan 等(孫豁然等，2007;朱超等，2010)．國(guó)內(nèi)對(duì)該方面的研究起步較晚，但也取得了一定的成果，如中地?cái)?shù)碼集團(tuán)推出的基于MapGIS平臺(tái)的資源量估算與礦體三維建模系統(tǒng)、中國(guó)地質(zhì)大學(xué)(武漢)開(kāi)發(fā)的GeoView軟件、中國(guó)地質(zhì)科學(xué)院開(kāi)發(fā)的探礦者軟件Minexplorer、北京東澳達(dá)科技有限公司自主研發(fā)的3DMine礦業(yè)工程軟件等。

羊拉銅礦近些年產(chǎn)量不斷增長(zhǎng)，為確保后續(xù)資源，加強(qiáng)羊拉地區(qū)找礦研究，快速找到資源儲(chǔ)備區(qū)尤為重要。為使羊拉銅礦節(jié)約成本，充分利用資源，對(duì)已有資料進(jìn)行再次開(kāi)發(fā)利用，運(yùn)用新技術(shù)成為一條擴(kuò)大找礦成果的重要途徑?；趪?guó)內(nèi)外三維地質(zhì)建模的研究，筆者運(yùn)用Micromine軟件對(duì)羊拉銅礦構(gòu)建了礦山三維實(shí)體模型和屬性數(shù)據(jù)庫(kù)，將礦山的空間數(shù)據(jù)實(shí)現(xiàn)數(shù)字化存儲(chǔ)和表達(dá)，并對(duì)已知礦體進(jìn)行了資源量估算，為礦山的進(jìn)一步開(kāi)采提供參考。

1 研究區(qū)概況

云南德欽羊拉銅礦地處滇、川、藏3省交界，雪域青藏高原南緣，橫斷山脈北段，屬高山深切割地區(qū)，地形為向東傾斜的單向坡，坡度范圍為25°～50°，局部地區(qū)達(dá)70°，高差約2 300多 m。

羊拉礦區(qū)地處金沙江接合帶中，區(qū)內(nèi)主要出露泥盆系(D)、下石炭統(tǒng)貝吾組(C1b)、上三疊統(tǒng)(T3)、第三系(E)和第四系(Q)地層。主要含礦地層為泥盆系中—上統(tǒng)里農(nóng)組一段(D2+3l1)，其次為泥盆系下統(tǒng)江邊組三段(D1j3)、二段(D1j2)。該區(qū)的礦體產(chǎn)出于(斑)巖體及其外接觸帶的矽卡巖，角巖化變質(zhì)碎屑巖中。含礦巖性主要為透輝矽卡巖、石榴石矽卡巖，角巖化變質(zhì)石英砂巖、花崗閃長(zhǎng)巖次之。

本區(qū)域內(nèi)的斷層較為發(fā)育，規(guī)模最大的F4斷層斜穿礦區(qū)中部，呈北東向展布。除F4外還有近北東向的F6、F8、F10斷層及層間破碎帶。

里農(nóng)礦段共圈出層狀、似層狀及透鏡狀的礦體12 個(gè)，其中工業(yè)礦體有 KT2(上)、KT2(下)、KT2-7、KT3、KT4、KT4-1、KT5 7個(gè)，礦體總體走向近南北，向西緩傾，一般傾角180°～35°，主要產(chǎn)出于里農(nóng)巖體西側(cè)的外接觸帶中。路農(nóng)礦段KT5-1礦體產(chǎn)于上部矽卡巖礦化層之上部，賦存于地層D2+3l1與D2+3l2之層間破碎帶。KT1礦體產(chǎn)于下部矽卡巖礦化層的頂部，加仁巖體東側(cè)內(nèi)外接觸帶。KT4礦體賦存于D2+3l1與D1j3之層間破碎帶，往北被F4斷裂錯(cuò)失(曾禮傳，2010;胡光龍等，2008)(圖1)。

圖1 研究區(qū)地質(zhì)簡(jiǎn)圖

2 三維建模技術(shù)方案

三維地質(zhì)建模涉及地表建模、勘探數(shù)據(jù)庫(kù)、斷層建模、地層建模、礦體建模等(圖2)。地表模型是利用地表高程數(shù)據(jù)建立不規(guī)則三角網(wǎng)(TIN)形成面模型?？碧綌?shù)據(jù)庫(kù)是利用軟件將整理好的井口信息表、測(cè)斜數(shù)據(jù)表、巖層數(shù)據(jù)表以及取樣分析表相互關(guān)聯(lián)建立的。斷層模型是通過(guò)將數(shù)據(jù)中的斷層線依次用不規(guī)則三角網(wǎng)連成的面模型。巖體和地層是通過(guò)將線連結(jié)成面，再運(yùn)用面切割體的方法建立的。礦體模型的建立采用剖面圈礦，然后將各剖面上的礦體輪廓對(duì)應(yīng)用線框連結(jié)起來(lái)的方法建立。綜合地質(zhì)模型的建模步驟是：先將收集到的數(shù)據(jù)分析、整理并建立數(shù)據(jù)庫(kù)對(duì)其進(jìn)行存儲(chǔ)、管理及操作。然后，根據(jù)研究區(qū)的Xmin、Xmax、Ymax、Ymin、Zmin、Zmax構(gòu)建立方體線框，用 DTM 切割該體建立三維空間研究區(qū)范圍線框。在此基礎(chǔ)上，對(duì)各地質(zhì)體運(yùn)用相應(yīng)的方法建模，為使相鄰地質(zhì)體的模型能夠達(dá)到無(wú)縫貼合，對(duì)各地質(zhì)體建模應(yīng)遵循從建立斷層模型→建立巖體模型→建立地層模型→建立礦體模型這樣一個(gè)順序。

圖2 地質(zhì)體三維建模流程圖

3 羊拉三維地質(zhì)模型建立

3．1 Micromine軟件簡(jiǎn)介

Micromine軟件是由澳大利亞Micromine礦業(yè)軟件公司開(kāi)發(fā)的一套專門處理勘探和采礦數(shù)據(jù)的商業(yè)三維礦業(yè)軟件。它具有以下特點(diǎn)：(1)具有良好的用戶操作界面和與其他相關(guān)軟件的數(shù)據(jù)接口;(2)能夠有效地管理、分析勘探資料及地質(zhì)編錄數(shù)據(jù)，并能將其進(jìn)行三維可視化表達(dá);(3)能高效便捷地在三維空間中圈定礦體，并提供多種方法進(jìn)行資源量計(jì)算及動(dòng)態(tài)管理分析。

3．2 勘探數(shù)據(jù)庫(kù)的建立及鉆孔的三維可視化

數(shù)據(jù)庫(kù)的建立是地質(zhì)體三維建模的基礎(chǔ)，經(jīng)過(guò)對(duì)收集的羊拉銅礦的勘探工程資料(包括鉆孔柱狀銅礦圖和平硐素描圖)進(jìn)行匯總、分析，運(yùn)用Micromine軟件建立羊拉勘探數(shù)據(jù)庫(kù)，該數(shù)據(jù)庫(kù)中包括井口信息表、測(cè)斜數(shù)據(jù)表、巖層數(shù)據(jù)表以及取樣分析表共4個(gè)表文件。在數(shù)據(jù)錄入并更新該數(shù)據(jù)庫(kù)后，通過(guò)Micromine的鉆孔軌跡可視化功能將這錄入的133個(gè)鉆孔的空間形態(tài)展示在三維空間中，根據(jù)鉆孔中不同巖性或不同品位值填充不同的顏色和花紋。

3．3 地表模型的建立

地表模型的建立有利于三維空間研究區(qū)范圍的確定，并對(duì)近地表的地質(zhì)體模型的建立進(jìn)行空間范圍限定。研究區(qū)1∶1萬(wàn)地形地質(zhì)圖和里農(nóng)礦段(等高距40 m)、路農(nóng)礦段1∶2 000地質(zhì)地形圖中的等高線數(shù)據(jù)(等高距10 m)為研究區(qū)DTM模型的建立提供了數(shù)據(jù)源。該區(qū)高程數(shù)據(jù)的特點(diǎn)是高精度數(shù)據(jù)僅覆蓋局部區(qū)域，為使模型達(dá)到較好的模擬效果，將兩種精度的數(shù)據(jù)相結(jié)合：(1)在MapGIS中將等高線數(shù)據(jù)從各地質(zhì)地形圖文件中提取出來(lái)導(dǎo)入到同一個(gè)Micromine的線文件中，將重疊區(qū)域中40 m間距的等高線刪除;(2)用Micromine中DTM模塊下的創(chuàng)建表面工具對(duì)等高線加密，使等高線的間距統(tǒng)一為10 m;(3)運(yùn)用創(chuàng)建DTM表面工具依據(jù)加密后的等高線數(shù)據(jù)創(chuàng)建DTM表面。

3．4 斷層模型的建立

斷層的出現(xiàn)會(huì)使地質(zhì)體之間的空間關(guān)系變得復(fù)雜，它對(duì)其他地質(zhì)體模型的建立具有重要的參考作用，因而先于其他地質(zhì)體建模。斷層數(shù)據(jù)需要從包含有斷層信息的地質(zhì)地形圖、區(qū)內(nèi)構(gòu)造綱要圖、勘探線剖面圖、中段平面圖中提取。在將這些圖件轉(zhuǎn)換為DXF格式或MapGIS明碼格式導(dǎo)入Micromine后，運(yùn)用“測(cè)量→坐標(biāo)系統(tǒng)變換→平面”工具，利用兩個(gè)已知點(diǎn)對(duì)圖件進(jìn)行坐標(biāo)系統(tǒng)變換，最后將各文件及鉆孔數(shù)據(jù)中描述斷層的數(shù)據(jù)提取出來(lái)，按斷層編號(hào)分別存儲(chǔ)。用線框?qū)⑼晃募械木€兩兩連接構(gòu)成面，再依據(jù)勘探地質(zhì)報(bào)告中該斷層的描述(包括斷層的性質(zhì)與規(guī)模、兩盤(pán)的位移、走向及傾角)對(duì)所建立的斷層面進(jìn)行修編，圖3中展示了區(qū)內(nèi)已建成的斷層模型。

3．5 巖體及地層模型的建立

巖體是一種“有根”的地質(zhì)體，在勘探過(guò)程中，工程往往只揭露一定深度巖體的信息，對(duì)它的建模一般采用將鉆孔數(shù)據(jù)和勘探線剖面圖、中段平面圖以及地形地質(zhì)圖中描述巖體上表面空間位置的線數(shù)據(jù)用不規(guī)則三角網(wǎng)(TIN)連成面，然后用該面切割研究區(qū)三維線框模型，所得面下的線框模型為巖體模型。由于巖體常會(huì)侵入地層中，因此，地層建模晚于巖體建模。地層模型的建立同樣采用面切體的方法進(jìn)行，用各地層的接觸面依次“切割”研究區(qū)三維線框模型，得到各地層模型(圖4)。

圖3 斷層模型圖

圖4 地質(zhì)體模型圖

3．6 礦體模型的建立

礦體模型(圖5)以鉆孔數(shù)據(jù)為支撐，同時(shí)參考現(xiàn)有勘探線剖面圖、中段平面圖等成果資料。礦體模型的建立需要經(jīng)過(guò)以下幾個(gè)關(guān)鍵環(huán)節(jié)：(1)需要根據(jù)鉆孔中的樣品分析數(shù)據(jù)，按照羊拉銅礦指定的礦體圈定原則以及對(duì)礦體空間分布特征的分析在三維空間中圈定各礦體在剖面上的輪廓線。(2)利用Micromine中“創(chuàng)建線框”工具將相鄰的兩個(gè)剖面上同一礦體的輪廓線依次連結(jié)。(3)相鄰兩個(gè)工程，一個(gè)見(jiàn)礦，另一個(gè)不見(jiàn)礦，礦體推至兩工程間距1/2處尖滅。(4)建立近地表的礦體模型時(shí)需要用DTM進(jìn)行礦體邊界的限定，對(duì)于高出地表的部分需要通過(guò)布爾運(yùn)算對(duì)礦體模型進(jìn)行修正。當(dāng)?shù)V體出現(xiàn)分支閉合現(xiàn)象時(shí)，可以通過(guò)建立輔助輪廓線和連結(jié)線來(lái)幫助礦體線框的連結(jié)。當(dāng)出現(xiàn)礦體被斷層切割的情況時(shí)，需要借助布爾運(yùn)算實(shí)現(xiàn)斷層模型對(duì)礦體模型的切割，以確保二者之間具有正確的空間拓?fù)潢P(guān)系。建模過(guò)程中還可以通過(guò)建立“連結(jié)線”、抽稀或加密輪廓線上的點(diǎn)來(lái)避免交叉三角形的出現(xiàn)和線框模型失真或扭曲。

圖5 里農(nóng)礦段礦體模型圖

4 羊拉里農(nóng)礦段銅礦體資源量估算

4．1 礦體空塊模型的建立

對(duì)于實(shí)體的屬性信息線框模型無(wú)法完整地表達(dá)，這就需要建立塊體模型。綜合考慮研究區(qū)地質(zhì)復(fù)雜程度、礦體連續(xù)程度、礦體開(kāi)采方式、勘探工程間距以及計(jì)算機(jī)硬件條件，確定了本次研究中塊體模型體元的尺寸：北×東×高=6 m×10 m×2 m。

4．2 數(shù)據(jù)預(yù)處理

數(shù)據(jù)的預(yù)處理包括數(shù)據(jù)轉(zhuǎn)換、特高值處理以及組合樣長(zhǎng)。運(yùn)用克里格插值的前提是樣本數(shù)據(jù)服從正態(tài)分布，經(jīng)分析，原始數(shù)據(jù)的自然對(duì)數(shù)服從正態(tài)分布，故對(duì)數(shù)據(jù)進(jìn)行自然對(duì)數(shù)變換。在估值過(guò)程中，特高值對(duì)估值結(jié)果影響明顯，為此，筆者采用頻率曲線法界定特高品位，將曲線拐點(diǎn)處的品位值取為特高值的下限，對(duì)那些獨(dú)立出現(xiàn)的高值用該下限值替代。資源量估算前需要對(duì)樣長(zhǎng)進(jìn)行組合，組合時(shí)應(yīng)最少拆分原始樣長(zhǎng)，故而選取原始樣長(zhǎng)的眾數(shù)1．5 m作為組合樣長(zhǎng)，并依據(jù)經(jīng)驗(yàn)將0．5 m作為最小組合樣長(zhǎng)。

4．3 變差函數(shù)分析

克里格估值法是通過(guò)變差函數(shù)對(duì)區(qū)域變量進(jìn)行結(jié)構(gòu)分析和樣品品位估值的(沈步明等，1994)。計(jì)算礦體的走向、傾向以及礦體的厚度方向上的實(shí)驗(yàn)變差函數(shù)，并進(jìn)行實(shí)驗(yàn)變差函數(shù)與理論模型的擬合。球型模型是資源量估值中最常用的模型，選擇球狀模型進(jìn)行擬合，以確定各方向上變異函數(shù)的塊金值、基臺(tái)值和變程。根據(jù)表1中的各變差函數(shù)理論模型參數(shù)，得到各方向?qū)嶒?yàn)變差函數(shù)和理論模型的擬合結(jié)果，如圖6、圖7、圖8所示。

表1 變差函數(shù)理論模型參數(shù)

圖6 走向方向的變差函數(shù)曲線

圖7 傾向方向的變差函數(shù)曲線

圖8 厚度向方向的變差函數(shù)曲線

4．4 搜索參數(shù)的設(shè)置

估算過(guò)程中，需要確定搜索參數(shù)，它對(duì)克里格估值結(jié)果有很重要的影響。通過(guò)劃分扇區(qū)限定每個(gè)扇區(qū)內(nèi)最少、最多數(shù)據(jù)點(diǎn)來(lái)排除由勘探工程的分布不均勻?qū)е碌臉悠反丶挠绊?孫玉建，2008)。設(shè)定每個(gè)扇區(qū)最多有6個(gè)點(diǎn)參與估值，這6個(gè)點(diǎn)為與待估點(diǎn)距離最近的6個(gè)點(diǎn)，最少參與估值的點(diǎn)數(shù)為2，少于2個(gè)點(diǎn)則待估點(diǎn)不被估算。筆者設(shè)置的搜索半徑分別為50、100、160 m，搜索半徑為160 m時(shí)所有的待估塊段都被插值完畢。

4．5 資源量估算

運(yùn)用普通克里格法對(duì)里農(nóng)礦段礦體進(jìn)行Cu品位估值，將品位邊界值分別設(shè)定為0．2%、0．3%、0.5%，礦石體積質(zhì)量分?jǐn)?shù)如表2所示，經(jīng)計(jì)算，得出里農(nóng)礦段的銅礦石總量為3 830．3萬(wàn)t，Cu金屬總量為18.3萬(wàn)t。

表2 礦石體積質(zhì)量分?jǐn)?shù)

用另外一種估算方法——距離冪次反比法再次對(duì)里農(nóng)礦體進(jìn)行估值和資源儲(chǔ)量計(jì)算，計(jì)算結(jié)果如表3所示。采用與克立格法相一致的搜索參數(shù)進(jìn)行品位估值，將其所得結(jié)果與所采用普通克里格方法估算所得的資源量結(jié)果進(jìn)行比較，以此來(lái)衡量運(yùn)用克里格法所估算結(jié)果的可靠程度。對(duì)比結(jié)果表明，兩種計(jì)算結(jié)果的相對(duì)偏差僅為2．31%，證明其可靠程度較好。圖9為運(yùn)用克里格法對(duì)Cu品位進(jìn)行差值后里農(nóng)礦段礦體Cu品位的分布圖。

表3 兩種方法資源量估算比較

圖9 里農(nóng)礦段礦體Cu品位分布圖

5 結(jié)論

以云南羊拉銅礦里農(nóng)礦段、路農(nóng)礦段為研究對(duì)象，在分析研究建模技術(shù)方案的基礎(chǔ)上，依據(jù)所收集到的鉆孔的柱狀圖、勘探線剖面圖、中段地質(zhì)圖和地形地質(zhì)圖，結(jié)合野外踏勘與認(rèn)識(shí)，利用Micromine軟件建立了研究區(qū)的地形、地層、斷層以及礦體的三維實(shí)體模型，直觀展示了區(qū)內(nèi)地表、斷層、地層、巖體以及礦體的空間展布以及相互的空間位置關(guān)系，對(duì)分析各地質(zhì)體之間的相互關(guān)系以及礦體的空間展布特征提供了幫助。并運(yùn)用克里格法對(duì)里農(nóng)礦段的銅礦體進(jìn)行品位估值，獲得了礦體的儲(chǔ)量，對(duì)礦山今后的生產(chǎn)具有一定的指導(dǎo)意義。

戴碧波，王李管，賈明濤，等．2007．三維數(shù)字建模技術(shù)在某銅礦山中的應(yīng)用［J］．地質(zhì)與勘探，43(3)：97 －101．

胡光龍，姜華，蔣靖，等．2008．云南德欽羊拉銅礦區(qū)礦床控礦因素及找礦遠(yuǎn)景［J］．金屬礦山，(12)：87－89．

黎楓佶．2010．西藏墨竹工卡縣甲瑪銅多金屬礦三維模型構(gòu)建和資源量估算——基于Micromine的應(yīng)用［D］．成都：成都理工大學(xué)．

盧大超，付友山．2010．三維礦產(chǎn)資源評(píng)價(jià)軟件Micromine在金屬礦山資源儲(chǔ)量估算中的應(yīng)用——以吉林舒蘭季德鉬礦為例［J］．世界地質(zhì)，29(3)：450 －457．

沈步明，沈遠(yuǎn)超．1994．金屬礦山地質(zhì)數(shù)據(jù)庫(kù)與地質(zhì)統(tǒng)計(jì)學(xué)［M］．北京：科學(xué)出版社．

孫豁然，徐帥．2007．論數(shù)字礦山［J］．金屬礦山，(2)：1 －5．

孫玉建．2008．地質(zhì)統(tǒng)計(jì)學(xué)在固體礦產(chǎn)資源評(píng)價(jià)中的若干問(wèn)題研究［D］．北京：中國(guó)地質(zhì)大學(xué)(北京)．

張新宇，肖克炎，劉光勝，等．2006．阿舍勒銅礦可視化儲(chǔ)量計(jì)算的指示克里格法應(yīng)用研究［J］．吉林大學(xué)學(xué)報(bào)：地球科學(xué)版，36(2)：305 －308．

鄒艷紅，戴塔根，毛先成．2008．廣西大廠礦田銅坑深部隱伏礦體立體定量預(yù)測(cè)建模與可視化研究［J］．地質(zhì)與勘探，44(2)：62 －66．

張思科，倪晉宇，高萬(wàn)里，等．2009．三維地質(zhì)建模技術(shù)方法研究——以東昆侖造山帶為例［J］．地質(zhì)力學(xué)學(xué)報(bào)，15(2)：201 －208．

朱超，吳仲雄，張?jiān)妴ⅲ?010．?dāng)?shù)字礦山的研究現(xiàn)狀和發(fā)展趨勢(shì)［J］．現(xiàn)代礦業(yè)，(2)：25 －27．

曾禮傳．2010．云南德欽羊拉銅礦區(qū)地層研究［J］．云南地質(zhì)，29(1)：84 －89．

ANDREA ZANCHI， SALVIFRANCESCAC， ZANCHETTA STEFANO，et al．2009．3D reconstruction of complex geological bodies：Examples from the Alps［J］．Computers ＆Geosciences，35(1)：49 －69．

CHEN JIANPING，LV PENG，WU WEN，et al．2007．A 3D prediction method for blind orebody based on 3D visualization model and its application［J］．Earth Science Frontiers，14(5)：54 －62．

3D geological modeling and reserves estimation in Yangla Copper Deposit in Yunnan

SUN Xiao1，CHEN Jian-guo1，2，F(xiàn)ANG Xiao-long1，MA Qiu-shi1

(1．Faculty of Earth Resources，China University of Geosciences(Wuhan)，Wuhan 430074，China;2．State Key Laboratory of Geological Processes and Mineral Resources of China University of Geosciences，Wuhan 430074，China)

Taking the fault，strata，rock mass and ore body in Yangla Copper Deposit as study targets，based on the analysis and study modeling technology plan，the available of borehole histogram，prospecting line profile，middle geological map and topographic geologic maps and in combination with field exploration and knowledge，the authors established a 3D solid model on landform，strata，fault and ore body in the study area by the use of Micromine software，directly exhibited the spatial distribution and mutual spatial locations for the surface，faults，strata，rock mass and ore body in the ore district．The authors estimated the grade of Linong ore section with Kriging method and acquired the reserve of the ore body．

3D geological modeling;Resource estimation;Micromine software;Yangla Copper Deposit;Deqin，Yunnan

TP319;P618．41

A

1674－3636(2012)03－0326－07

10．3969/j．issn．1674-3636．2012．03．326

2012－05－20;編輯：蔣艷

云銅橫向合作項(xiàng)目(2010026410)資助成果

孫瀟(1988— )，女，碩士研究生，地質(zhì)工程專業(yè)，E-mail：sunxiao-76116@163．com