内蒙古北山地区成矿信息提取技术及成矿预测研究
详细信息    本馆镜像全文|  推荐本文 |  |   获取CNKI官网全文
摘要
内蒙古北山地区位于哈萨克斯坦板块与华北-塔里木板块的交接部位,区内古生代地层出露广泛,总体构造以近东西、北西向为主,岩浆作用强烈,具备大规模成矿的条件和潜力。由于区内自然条件艰苦,地质工作程度很低,加上近十几年找矿投入严重不足,致使大量找矿线索得不到进一步评价。因此,为实现内蒙古北山地区金属矿勘查的重大突破,有效缓解资源短缺,支持地区经济发展,在该区加强快速找矿评价方法研究具有重要的现实意义。
     本论文充分收集了内蒙古北山地区黑鹰山、红柳大泉、六驮山、路井、石板井和五道明六幅1:20万的地质矿产标准图幅,1:5万(部分1:10万)的航磁数据和1:100万的重力数据,1:20万化探数据扫面中测定的8种化学元素以及遥感数据,同时在MapGis平台上进行数字化。对于物探数据,使用了先进的Geosoft软件进行了处理,除了重力布格异常、航磁化极异常、方向导数异常、相应的空间延拓异常,还应用了视磁化率填图、解析信号分析、梯级带滤波增强技术等数据处理方法,提取出与基性岩体、中酸性岩体有关的信息,圈定了多处磁性块体,对本区的深部构造、莫氏界面、基底起伏进行了分析研究。对Au、Ag、Cu、Cr、Mo、Ni、Pb、Zn等8种化学元素的化探数据进行了背景场和异常场的分离,确定了化探异常的空间分布。使用ENVI软件处理遥感图像,开展了遥感多光谱蚀变信息提取找矿研究,并划分了区内环形和线性构造。在区域地质成矿理论的指导下,使用MRAS矿产资源预测软件,利用“证据权”法进行了成矿预测,为内蒙古北山地区指出了金属矿产的找矿方向。论文取得的主要成果如下:
     1.使用Geosoft软件处理区域重磁数据,首次根据重磁资料的地质解释,将本区划分为哈萨克斯坦板块东缘构造区和塔里木板块北缘两大Ⅰ级构造区以及6大Ⅱ级构造小区。推断划分北西向(近东西向)深大断裂构造带5条和北东向断裂构造7条。初步判断、圈定多处磁性块体,并区分出基性岩体、中酸性岩体的分布。
     2.在MRAS矿产资源评价软件平台上确定地球化学8种主要元素的异常下限、均值、标准差,并进行了相关性分析和聚类分析。研究了主要元素的组合特征和分布趋势,划分出5个地球化学异常区和数处有价值的高值异常区(找矿靶区)。
     3.使用ENVI软件处理多光谱卫星Aster数据,通过图像大气校正、遥感数据最大噪声分离、混合像元分解、矿物端元选取、比值法、特征主组分法、光谱角法、光谱特征拟合等诸多方法,系统提取了研究区铁氧化物、铝羟基类蚀变、黑云母、绿泥绿帘石类蚀变,针对重点研究区又提取了硅化、萤石化、云英岩化等蚀变信息。
     4.在MapGis平台上建立了区域地质、化探、航磁和遥感蚀变信息数据库,利用这些多元信息,在MRAS矿产资源评价软件平台上展开了矿产资源综合信息成矿预测。在“证据权”法成矿预测的基础上,结合区域成矿规律,共圈出13个找矿靶区,为内蒙北山地区金属矿勘查指出了方向。
     5.在对各类主要矿种成矿模式研究的基础上,提出找矿前景较大的矿种主要是铁、钼(铜)、钨(锡)、锑、金和铬。对铁矿床的寻找主攻方向是海底火山喷流沉积型,其次是矽卡岩型。钼(铜)矿床主攻方向是斑岩型。钨(锡)、锑矿床主攻方向是石英脉型和蚀变岩型。金矿主攻方向是石英脉型和蚀变岩型。铬铁矿主攻方向是镁质超基性岩浆晚期熔离型。
Beishan area, Inner-Mongolia is located in the joint between Kazakstan plate and NorthChina plate. Paleozoic layer distribute abroad and magma active violence so it have goodcondition and great potential for large metallogenesis. Many clues for prospecting in Beishan areacould not be used to actualize more appreciation due to hard natural surrounding and lower gradeof geology survey and shortage of prospecting financing. Hence, it have important realisticsignification to strengthen fast prospecting approach research because it can mitigate resourcesshortage effectively and support regional economy develop fleetly and it also can realize greatbreach for Inner-Mongolia's metallic mine prospecting.
     Many multi-source geological data was collected in this research, these data include: 1:20000 standard mineral map sheet such as Heiyingshan,Hongliudaquan,Liutuoshan,Lujin,Shibanjinand Wudaoming,1:50 000(partial 1:100 000) aero magnetic data, 1:10 000 000 gravity data,1:20 0000 geochemical data about 8 kinds of element and remote sensing data. Whole these data wasdigitalized in MapGis software. For geophysical data, all data was processed in advancedgeophysical software-Geosoft. The process method include susceptibility mapping, analytic signaland grade flitter enhance method suppressed noise effectively and stand out useful information.By these methods, some information about basic rock and acidic rock was extracted and manymagnetic plates were also confirmed. The deep structure, Moho discontinuity and undulation ofbase were also be researched in this paper. Background field and anomaly field forAu,Ag,Cu,Cr,Mo,Ni,Pb and Zn were separated and the space distribution of geochemicalanomalies were confirmed. The prospecting alternation information extract from multi-spectrumremote sensing were carry out based on image process in ENVI software and circle wisedstructure and lineal structure in research area were plot out. Under the guide of regional geologicalmetallogenic theory, the "evidence weight method" was used to metallogenic prognosis in mineralresources forecast software-MRAS.The forecast result point out the direction for metallic mineprospecting in Beishan area,Inner-Mongolia. The main production of this paper as follows:
     1. Regional gravity and magnetic dada were processed in Geosoft software, the research areawas divided into 2 bigⅠstructure area and 6Ⅱstructure area among east edgestructure area of Kazakstan plate and north edge of Talimu plate based on the geologicalinterpret ion from gravity and magnetic data. 5 deep west north (east west) faults and 7 eastnorth faults were inferred. Many magnetic plates were estimated primarily and distribution ofbasic rock and acidic rock were distinguished.
     2. On the mineral resource forecast software MRAS, anomaly lower limit, average, STDEV of8 kinds of geochemical element were confirmed and their correlation and clusteringrelationship were also be analyzed. Combination character and distributional trend of mainelement were researched and 5 high concentration anomalies and several valuable highanomalies were distinguished.
     3. ENVI software was utilized to process multi-spectrum satellite Aster data. The processmethod include Spectral feature fitting, Spectral fitting filtering, Multi-Range SFF, radiomethod, Maximum Noise Fraction and so on. By these method, some kinds of alternationinformation include iron oxide and hydroxide and biotite and epidote were extracted from multi-spectrum information and contraposed key research area some other alternationinformation include silification and fluorite and greisenization were also be extracted.
     4. A multi-source geological information spatial database and attribute database was built basedon MAPGIS software and it realized synthetical management effectively for complexgeological data of Beishan. On the mineral resource forecast software -MARS, syntheticalinformation forecast was completed by utilized this multi-sources information. On the basisof "evidence weight method",integrated with regional metallogenic rules, 13 target forprospecting were confrmed and this result point out the direction for metallic prospectingof Beishan,Inner-Mogolia.
     5. One the basis of metallogenic pattern research for all main kinds of metallic mine, themain mine type that have good foreground for prospecting in Beishan area are iron ore,molybdenum(copper), tungsten(stannum), antimony, golden and chrome. The mainprospecting direction for iron mine are submarine plume deposit type and skam type,for molybdenum(copper) is porphyry type,for ungsten(stannum) and antimony are quartzesvein type and alternation rock type, for golden is also quartzes vein type and alternationrock type, for chrome iron is ultrabasic magmatic fusion separated type.
引文
1. TURCOTTED L. A fractal approach to the relationship between ore grade and tonnage [J]. Economic Geology, 1986, (81): 1528-1532.
    2. ZHAO Pengda, CHEN Jjianping, CHEN Jianguo, et al. 2001. Introduction to "Three—step"Digital Mineral Prospecting and assessment. Proceedings of the International Symposium, 1-5.
    3. CHEN Qiu-ming, AGTERBERQ F P, BALLANTYNE S B. The separation of geochemical anomalies from background by fractal methods[J]. Journal of Geochemical Exploration, 1994, 51: 109-130.
    4. Y Zheng, Q Zhang, G Zuo, Y Wang, et al. 1996. Great Jurassic thrust sheets in Beishan (North Mountains)-Gobiareas of China and southern Mongolia [J]. Journal of Structural Geology [J]. Vol. 18. No. 9: 1111-1126.
    5. Agterberg F P, Bonham-Carter G F, Wright D F. 1990.Statistical Pattern Integration for Mineral Exploration, In: Gaal G, Merriam D F eds, Computer Applications for Mineral Exploration in Resource Exploration, Oxford: Pergamon Press, 1-21.
    6. Bonham-Carter G F, Agterberg F P. 1990. Weights of evidence: a new approach to mapping mineral potential, statistical applications in the earth sciences. Canada: Geological Survey of Canada.
    7. Mandelbrot B B. 1982.The fractal geometry ofNature. W. H. Freeman and Co., New York.
    8. Boardman, J. W., 1998, Leveraging the High Dimensionality of AVIRIS Data for Improved Sub-Pixel Target Unmixing and Rejection of False Positives: Mixture Tuned Matched filtering, Summaries of the Seventh Annual JPL Airborne Earth Science Workshop, Pasadena, CA.
    9. Boardman, J. W., Kruse, F. A., and Green, R. O., 1995, Mapping target signatures via partial unmixing of AVIRIS data: in Summaries, Fifth JPL Airborne Earth Science Workshop, JPL Publication 95-1, v. 1, pp. 23-26.
    10. Lee J. B., Woodyatt S., and Berman M., 1990, Enhancement of high spectral resolution remote-sensing data by a noise-Adjust principal components transform, IEEE Trans. Geosci. Remote Sens., 28(3): 295-304.
    11. Green A. A., Berman M., Switzer P., et al, 1988, A transformation for ordering multispectral data in terms of image quality with implications for noise removal, IEEE Trans. Geosci. Remote Sens., 26(1): 65-74.
    12. Hsu, S.-K., 2002, Imaging magnetic sources using Euler's equation: Geophysical Prospecting, 50, 15-25.
    13. Janet E. Tilden, David A. Ponce, Jonathan M. G. Glen, and Kathleen D. Gans Gravity and magnetic data along a seismic refraction-reflection line in northwest Nevada and northeast California Open-File Report-U. S. Geological Survey (2005), OF 2005-1446: 12 pp.
    14.G.F.Bonham-Carter,F.P.Agterberg,and D.F.Wright.Nova Scotia地区金矿勘探地质资料的综合利用[M].地理信息系统及其在地学领域中的应用译文集,地矿部科技司计算机办公室,1993.
    15. G. P. Watson, A. N. Rencz and G. F. Bonham-Carte, Geographic information system are being applied to mineral resource assessment in Northern New Brunswick [J]. GEOS 1989, 18 (1).
    16. Lesley Wyborn et al., Using GIS formineral potential evalution in Areaswith fewknowmineral occurrences [J]. The second forum on GIS in the geosciences, AGSO, p. 199-211, 1995.
    17. Lyle A. Burgess, Recent applications and research into mineral prospective mapping using GIS [J]. Proceeding of third national forum on GIS in the geosciences, AGSO, p. 121-129, 1997.
    18. C. M. Knox-Robinson, L. A. I. Wyborn, Towards a holistic exploration strategy: Using geographic information systems as a tool to enhance exploration [J]. Australian journal of earth sciences, 1997, 44: 453-463.
    19. Roest, W. R., Verhoef, J., and Pilkington, M., 1992, Magnetic interpretation using the 3-D analytic signal, Geophysics, 57~116.
    20. Baranov W. A new method for interpretation of aeromagnetic maps: pseudo—gravimetric anomalies. Geophysics, 1957, 22(2): 359~383.
    21. L. E. Reed Gravity and magnetic three-dimensional (3D) modelling; Discover Abitibi Initiative Open File Report-Ontario Geological Survey (2005), 6163: 40 pp.
    22. Yaoguo Li 3D inversion of potential-field data and its applications in exploration problems The Bulletin of the Houston Geological Society (March 2005), 47(7): 23.
    23. V. I. Starostenko, V. N. Zavoisky, and O. V. Legostaeva The forward problem of magnetic survey for 3-D bodies of complex configuration with anisotropic magnetic susceptibility Izvestiya-Russian Academy of Sciences. Physics of the Solid Earth (July 2005), 41(7): 585-593.
    24. Philip J. Heath, Stewart Greenhalgh, and Nicholas G. Direen Modelling gravity and magnetic gradient tensor responses for exploration within the regolith Exploration Geophysics (Melbourne) (December 2005), 36(4): 357-364.
    25. G. V. Prostolupov, V. M. Novoselitskiy, V. N. Koneshov, and G P. Shcherbinina Gravity and magnetic field interpretation based on transformations of horizontal gradients in the VECTOR system Izvestiya-Russian Academy of Sciences. Physics of the Solid Earth (June 2006), 42(6): 530-535.
    26. Rui Z Y, Liu Y L, Wang L S, et al. 2002. The porphyry copper belt and its tectonic framework of eastern Tianshan, Xinjiang[J]. Acta Geologica Sinica, 6(1): 83~94 (in Chinese with English abstract).
    27. Xu X W, Wen C Q, Mao Y S, et al. 1998. Ore types and genes significance of Xiaorequanzi copper deposit, Xinjiang[J]. Minera Deposits, 17 (Supp.): 679~682 (in Chinese).
    28. Yan W Y. 1985. Characteristics of early Carboniferous are volcanic and minerals product in eastern rianshan Mts[J]. Xinjiang Geol., 3(2)49~58(in Chinese with English abstract).
    29. Yang J G, Ji J S and Lu D Y. 1998. Geological characteristics and minerogenetic condition of Xitan gold deposit, east Tianshan[J]. Mineral Deposit, 17 (Supp): 349~352 (in Chinese with English abstract).
    30. Zhang L C, Zeng Z R and Hart Z X. 1994. Study on geochemistry of Kangguer gold deposit in Xinjiang[J]. Contrib. Geol. & Explr., 9 (4): 12~21 (in Chinese with English abstract).
    31. Zhang L C, Ji J S, Li H Q, et al. 2000. Geochemical characteristics and source of two_type ore_forming fluids in Kangguertage gold ore belt, east Tianshan[J]. Acta Petrologica Sinica, 16(4): 535~541(in Chinese with English abstract).
    32.地质部航空物探大队九0五队,内蒙古额济纳旗地区航磁普查报告.1980年.
    33.地质部航空物探大队几0五队.甘肃北山地区航空物探结果报告.1966年.
    34.甘肃省地质局物探队.1内蒙黑鹰山地区铁矿普详查物探工作结果报告.1966年.
    35.内甘肃省二0一物探队.蒙黑鹰山铁矿及其外围—带物(化)探成果报告.1959年.
    36.美国地质训查局,美国本土矿产资源评价计划系列报告之二[M].地理信息系统及其在地学领域中的应用译文集,地矿部科技司计算机办公室,1993.
    37.美国地质调查局,美国本土矿产资源评价计划系列报告之一[M].地理信息系统及其在地学领域中的应用译文集,地矿部科技司计算机办公室,1993.
    38.内蒙古自治区地质训查院.内蒙古1:20万黑鹰山幅区域重力调查技术说明书.2002年.
    39.王世称,陈永良,夏立显.综合信息矿产预测理论与方法.北京:科学出版社,2000.
    40.赵鹏大,陈建平,张寿庭.“三联式”成矿预测新进展.地学前缘,2003,10(2):454~463.
    41.赵鹏大,陈永清,等.地质异常成矿预测理论与实践.武汉:中国地质大学出版社,1999,2~10.
    42.陈祖伊,张学权,张昭明等.1993北山区域地质发展史和金矿区域成矿模式 铀矿地质,9(2):65~75.
    43.崔惠文,陈祖伊.1996.甘肃北山地区金矿地质.北京:地质出版社.
    44.黄典豪,董群英,甘志贤.1997.额济纳旗老硐沟氧化—淋滤型金矿床成矿特征.贵金属地质,6(2):93~100.
    45.何世平,任秉琛,姚文光等.甘肃内蒙古北山地区构造单元划分.西北地质,35(4):30~40.
    46.刘雪亚,王荃.1995.中国西部北山造山带的大地构造及其演化.地学研究,第28号,37~48.
    47.穆治国,刘驰,黄宝玲等.1994甘肃北山古生代造山带地壳演化的同位素和稀土元素同位素特征.北京大学学报,30(2):202~214.
    48.聂凤军,江思宏,赵省民等.2000北山地区金矿床类型和生成演化.长春科技大学学报,第30卷,金矿专辑,21~26.
    49.聂凤军,江思宏,刘妍等.2002阿拉善东七一山萤石矿床Sm-Nd同位素年龄及其地质意义.矿床地质21(1):10~15.
    50.宋学信.1989.黑鹰山铁矿床的岩石学、地球化学特征及成因.中国地质科学院矿床地质研究所所刊,1:114~205.
    51.王伏泉.1996白山堂铜矿床两期有关岩体的Rb-Sr等时线年龄及其稀土配分特征.矿物岩石地球化学通报,15(3):187~190.
    52.王伏泉.1996北山地区海相火山岩中铜—多金属成矿作用的构造背景.大地构造与成矿学,20(4):319~318.
    53.王玉往等.1997北山地区各时代火山岩组合特征及分布.中国区域地质,16(3):297~304.
    54.杨合群,李英,杨建国等.2006.北山造山代的基本成矿特征.西北地质,39(2):78~95.
    55.杨敏之.2003.甘肃北山富铁矿床类型、成矿系列及成矿预测.地质找矿论丛,18(1):6~12.
    56.左国朝,李茂松,等.1996.甘蒙北山地区早古生代岩石圈形成与演化.兰州:甘肃科学技术出版社.
    57.左国朝等.1987.北山地区白云山—洗肠井一带早古生代蛇绿岩带.
    58.推覆构造及蛇绿混杂岩体的发现,《中国北方板块构造论文集》,第二集,北京:地质出版社,51~57.
    59.左国朝等.1987.北山洗肠井地区寒武系—中奥陶统浊积岩及其火山-沉积组合特征.沉积学报,5(2):63~72.
    60.左国朝,刘春燕,冯永忠等.1992.北山地区早山一带澄江花岗岩的发现及其地球化学特征.现代地质学研究文集(上),北京:地质出版社68~73.
    61.周国庆,赵建新,李献华.2000.内蒙古月牙山蛇绿岩特征及形成的构造背景:地球化学和Sr-Nd同位素制约.地球化学,29(2):108~119.
    62.赵安生,李景春.1993.甘肃北山北带金矿区域成矿特征。贵金属地质,2(2):104~109.
    63.朱裕生.1984.矿产资源评价方法学导论.北京:地质出版社.
    64.卢作祥,范永香,刘辅臣.1981.成矿规律和成矿预测学.武汉地质学院勘探教研室.
    65.朱裕生,肖克炎等.1997.成矿预测方法.北京地质出版社.
    66.孙殿卿,高庆华.1987.隐伏矿床预测.北京:地质出版社.
    67.聂风军,江思宏,等.北山地区金属矿床成矿规律及找矿方向.2002,北京:地质出版社.
    68.翟裕生,邓军,彭润民.中国区域成矿若干问题探讨.矿床地质,1999,18(4):323~332.
    69.赵鹏大.地质异常与成矿预测.当代矿产资源勘查平均的理论与方法.北京:地震出版社,1999,98~106.
    70.刘承祚.1999.地质专家系统与成矿预测[M].当代矿产资源勘查评价的理论与方法,北京:地震出版社.
    71.张义等.中蒙边境欧玉陶勒盖大型铜—金矿床的发现及对找矿勘查工作的启示.地质通报,2003,22(9):708~712.
    72.刘石年.1993.成矿预测学.长沙:中南工业大学出版社.
    73.肖克炎,张晓华,等.2000.矿产资源GIS评价系统.北京:地质出版社,1~142.
    74.袁峰,周涛发,岳书仓.2003.GIS矿产资源预测中的证据权重法.黄金地质,9(3):75~77.
    75.范永香,曾键年,刘伟.2004.论成矿预测的理论体系.湖北地矿,18(1):4~10.
    76.欧阳玉飞,黄满湘.2006.利用多元信息进行成矿预测.西部探矿工程,7:139~142.
    77.黄海峰,姚书振,丁振举.2003.基于GIS的证据权重法在成矿预测中的应用—以甘肃省岷县-礼县地区的金矿预测为例.地质科技情报,22(3):77~82.
    78.张寿庭,赵鹏大,等.2003.多目标矿产预测评价及其研究意义.成都理工大学学报(自然科学版),30(5):441~446.
    79.郝社峰,刘汉湖,杨武年.2005.多元信息成矿预测模型分析及实例应用.江苏地质,29(4):234~238.
    80.申维.2002.分形混沌与矿产预测.北京地质出版社.
    81.李宝强,孙泽坤.2004.区域地球化学异常信息提取方法研讨.西北地质,37(1),102~107.
    82.李随民,姚书振.2005.基于MAPGIS的分形方法确定化探异常下限.地球学报,26(2):187~190.
    83.左国朝等.北山地区早古生代板块构造特征.地质科学,1990,4:305~314.
    84.聂风军,江思宏.中蒙边境塔林大型金矿化带发现对我们的启示.内蒙古地质,2000,95(2):20~23.
    85.聂风车等.中蒙边境及邻区斑岩型铜矿床地质特征及成因.矿床地质,2004,23(2):176~189.
    86.左国朝,刘义科,刘春燕.2003.甘新蒙北山地区构造格局及演化.甘肃地质学报,12(1),1~15.
    87.李士勤,苏君伯,詹定省.1971.1:20万区域地质调查报告(白云鄂博幅)[R].内蒙古呼和浩特:内蒙古区调队.
    88.王金龙,李福亮.1995.内蒙古达茂旗赛乌素金矿床及其外围找矿预测[R].内蒙古呼和浩特:武警黄金第二支队.
    89.尹忠.2002.内蒙古金矿床分布规律及找矿方向[J].内蒙古地质,(4):8214.
    90.彭巨贵,张发荣,赵福昌.2006.甘新蒙北山地区成矿带划分和基本特征对比.地球科学与环境学报,28(3),11~18.
    91.陈毓川.1999.中国主要成矿区带矿产资源远景评价[M].北京:地质出版社.
    92.姜常义等.2006.新疆北山地区中坡山北镁铁质岩体岩石地球化学与岩石成因.岩石学报,22(1),115~139.
    93.王惠等.2007.内蒙古西部北山地区甜水井—青山一带双堡塘组时代讨论.地质通报,26(2),174~182.
    94.鲍庆中,等.2006.内蒙古东南部西乌珠穆沁旗地区石炭纪—二叠纪岩石地层和层序地层.地质通报,25(2):572~579.
    95.赵月明.2002.内蒙古北山重点成矿远景区主要地质特征与矿产资源潜力评价方向.内蒙古地质,1,1~5.
    96.王福同等.2001.物探在新疆土屋地区铜矿找矿中的应用—兼谈斑岩铜矿“三位一体”的找矿模式[J].中国地质,28(3),3.
    97.陈伯林等.2000.北山地区发现韧性剪切带型金矿床[J].中国区域地质,19(3),8.
    98.朱创业.2000.成矿系统研究现状及发展趋势.成都理工学院学报,27(1),50~53.
    99.翟裕生.1999.关于矿床学研究前景的探讨[J].矿床地质,18(2):146~152.
    100.李人澍.1996.成矿系统分析的理论与实践[M].北京:地质出版社,1~240.
    101.胡云中,侯增谦.1998.当代主要金属矿产资源勘查与研究的发展态势(上)[J].地质科技管理,(1):1~7.
    102.翟裕生.1999.论成矿系统[J].地学前缘.6(1):13~26.
    103.内蒙古自治区地质调查院.内蒙古石板井地区1:20万区域重力调查技术说明书.2004年.
    104.严良俊,胡文宝,姚长利.2006.重磁资料面积处理中的滤波增强技术与应用.勘探地球物理进展,29(4):102-104.
    105.王芳.2004.重磁勘探方法新技术.地质与资源,13(3),184~186.
    106.管志宁,郝天姚,姚长利.2002.21世纪重力与磁法勘探的展望.地球物理学进展,17(2),237~244.
    107.邵积东.1998.内蒙古大地构造分区及其特征.内蒙古地质,2:1~23.
    108.王世称,陈永清.成矿系列预测的基本原则及特点.地质找矿论从,1994,9(4):79~85.

© 2004-2018 中国地质图书馆版权所有 京ICP备05064691号 京公网安备11010802017129号

地址:北京市海淀区学院路29号 邮编:100083

电话:办公室:(+86 10)66554848;文献借阅、咨询服务、科技查新:66554700