用户名: 密码: 验证码:
MEMS-IMU误差分析补偿与实验研究
详细信息    本馆镜像全文|  推荐本文 |  |   获取CNKI官网全文
摘要
为了首次在国内将MEMS-IMU应用在某战术级项目的姿态测量上,本文研究了在较为复杂的实用环境下,MEMS-IMU及其中的MEMS惯性传感器的输出误差的机理及相应的补偿方法。由于组成MEMS-IMU的MEMS惯性传感器的工作方式与其它类型的传感器大不相同,因此开展的MEMS惯性传感器在实验条件下的误差机理研究,是补偿方法的依据。
     首先,本文介绍了MEMS-IMU的工作原理和研制的战术级精度MEMS-IMU样机,分析并测量了样机的各项误差,包括MEMS惯性传感器的误差与集成产生的误差,评估了在实验条件下各项误差导致的导航解算误差的大小,得到了本文的研究重点是传感器温度漂移误差与加速度灵敏度误差的抑制。
     然后,本文提出了用以分析微结构误差导致的输出误差的统计参数法,建立了微传感器输出误差与制造过程产生的尺寸误差的联系,可以用以预测给定制造工艺条件下传感器所能达到的误差指标。接着在考虑加工工艺步骤的基础上,分析了微加速度计温度漂移的各种可能的原因,重点分析了硅玻璃热膨胀系数失配导致的相对位移与微加速度计加工误差导致的U形梁刚度不对称综合作用的贡献,并对温度漂移的进行解析计算;建立了微加速度计结构-静电-温度多物理场耦合有限元模型并进行了温度漂移的分析仿真,还进行了温度实验,并对分析,仿真和实验的结果进行了对比验证。
     接着,本文分析了产生双极解耦微陀螺加速度灵敏度与温度漂移可能的原因,重点分析了在驱动力和惯性力作用下加工随机误差导致微陀螺U形梁刚度不一致使得中心质量块发生了扭转,使检测电容差不为零的现象。建立了微陀螺结构有限元模型,得到了微陀螺温度漂移与加速度灵敏度的仿真结果,再通过温度实验与离心实验对误差项的分析,仿真和实验结果进行对比验证。
     最后,本文针对于使用中的具体问题展开了实验研究。提出了一种MEMS-IMU整机标定补偿方法。提出了使用磁场传感器测量超过微陀螺量程的角速率的方法。介绍了进行的飞行实验,并结合飞行实验的结果验证了本文中提出的整机误差标定补偿方法,磁场测量大角速率的方法,温度漂移误差和加速度灵敏度误差的补偿方法。实验取得圆满成功,达到要求的姿态测量精度。
In order to implement the attitude measurement of the tactical ballistic trajectoryweapon, this dissertation focused on the research of the Micro-Electro-MechanicalSystems Inertial Measurement Unit(MEMS-IMU) error mechanism and correspondingcompensation method under the actual combination condition, including large angle rate,acceleration and temperature variation. Since the principal of the MEMS inertial sensorwhich comprised the MEMS-IMU are different from the other kinds of inertial sensors,the traditional compensation models are not optimized. The research on theMEMS-IMU error mechanism is the basis of the compensation method.
     This dissertation firstly introduced the principle of MEMS-IMU and the fabricatedtactical grade accuracy MEMS-IMU. The errors of the MEMS-IMU were estimatedunder the certain condition to provide the importance of certain errors includingtemperature drift and g sensitivity. Part of this dissertation would focus on theinterpretation of these errors.
     Then, this dissertation proposed a statistical parameter method to analyze thesensor output error caused by the manufactured dimension variation. It connected thevery beginning of the sensor fabrication and output of the sensor. It also can predict therange of the output error with the certain fabrication process. The possibleinterpretations of the accelerometer thermal drift based on the width parameters and thefabrication process were analyzed. The contribution of the stiffness asymmetry of theU-springs of the structure and relative displacement caused by the mismatch in thermalexpansion coefficients between the Pyrex glass substrate and heavily boron-dopedsilicon structure was investigated, modeled and simulated. The analytical model andmultiphysics simulation model were established to calculate and simulate the thermaldrift. The temperature experiments were carried out to compare and verify the analysisresult.
     Next, this dissertation analyzed the possible error mechanism of the thermal driftand g-sensitivity of doubly decoupled MEMS gyroscope. The phenomenon of thenon-zero detection capacity difference caused by the twist of the center mass, whichwas caused by the acceleration force and driven force and the U-spring width variation,was mainly investigated. The finite element model of the gyroscope’s structure was established and simulated. The simulation result of zero rate output, thermal drift andg-sensitivity were given and compared with the analysis value and experiment result.
     Furthermore, this dissertation proposed a magnetometer aid over-range rotationrate combination measurement method which can compensate the large angle rate errorcaused by the gyroscope. The overall MEMS-IMU calibration method was proposedbased on the analysis of the error principle above. A compensation method for thecalibration of the MEMS-IMU was also designed to improve the performance ofMEMS-IMU.
     Finally, this dissertation introduced the successful flight experiment which verifiedthe entire work of this dissertation. The overall compensation method, themagnetometer aid method and the analysis of the error of MEMS-IMU were verified bythis experiment.
引文
[1] Stephen D. Senturia, Microsystem Design[M], Kluwer Academic Publishers,2003
    [2] www.analog.com
    [3]陈哲,捷联惯导系统原理[M].北京:宇航出版社,1986
    [4]邓正隆.惯性导航原理[M].哈尔滨:哈尔滨工业大学出版社,1994
    [5] K.R. Inertial Navigation System Analysis[M]. New York: Wiles-Interscience,1971
    [6]张树侠,孙静.捷联式惯性导航系统[M].北京:国防工业出版社,1992
    [7]刘俊,微惯性技术[M],电子工业出版社,2006
    [8] William C.Tang, MEMS Program at DARPA[C], MEMS Components and Application forIndustry, Automobiles, Aerospace and Communication, Proceeding of SPIE Vol.4559(2001)
    [9] David W.Panhorst,Vicki Lefevre and Lynne K.Rider, Micro Electro-MechanicalSystems(MEMS) Inertial Measurement Unit(IMU)Common Guidance Program[J],Ferroelectrics,2006,342:p205~211
    [10]李荣冰,刘建业,曾庆化,华冰,基于MEMS技术的微型惯性导航系统的发展现状[J],中国惯性技术学报,2004,12(6),pp88-94
    [11] Gun-hard, low-cost MEMS IMU for the most demanding projectile guidance requirements-HG1930Datasheet, http://www51.honeywell.com/aero/common/documents/myaerospace-catalog-documents/Missiles-Munitions/HG1930_Datasheet.pdf
    [12] MEMS IMU AIST-330/AIST-330-G Datasheet,http://www.isense.ru/documents/aist_330_datasheet_(eng).pdf
    [13] DMU02-Six Degrees of Freedom Dynamics Measurement Unit Datasheet,http://www.sssj.co.jp/en/products/gyro/pdf/DMU02_SPEC.pdf
    [14] Six Degrees of Freedom Inertial Sensor:ADIS16367Datasheet,http://www.analog.com/static/imported-files/data_sheets/ADIS16367.pdf
    [15]牛小骥,微机械姿态测量单元及其用于卫星电视天线稳定的研究[D],清华大学博士论文,2002
    [16]李荣冰,刘建业,孙永荣,MEMS-IMU构型设计及惯性器件安装误差标定方法[J],中国惯性技术学报,2007,15(5):526-529
    [17] Chang-Sun Yoo, Iee-Ki Ahn, Low cost GPS/INS sensor fusion system for UAVnavigation[C], Digital Avionics Systems Conference,2003
    [18] Meiling Wang, Data Fusion of MEMS IMU/GPS Integrated System for Autonomous LandVehicle[C],2006IEEE International Conference on Information Acquisition
    [19]尚捷,MEMS-IMU及其与GPS组合系统设计与实验研究[D],清华大学工学博士论文,2005
    [20]华冰,刘建业,李荣冰,赵伟,刘瑞华,余度MEMS-IMU/GPS组合导航系统[J],南京航空航天大学学报,2007,39(5):570-575
    [21]张琳,初海彬,张乃通.卫星导航用户机空时滤波抗干扰技术[J].北京航空航天大学学报,2007,30(1):76-80.
    [22]汪捷,赵学军,导航对抗背景GPS信号宽带压制干扰分析[J].海军工程大学学报,2010,22(6):30-34
    [23]向培胜, GPS抗干扰技术综述[J],电讯技术,2009,49(10):93-96
    [24]李跃,邱致和.导航与定位[M].第2版.北京:国防工业出版社,2008
    [25] Jun Liu, Yunbo Shi, Wendong Zhang, Micro Inertial Measurement Unit based integratedvelocity strapdown testing system[J], Sensors and Actuators A:Physical,2004,A112(1):44-48
    [26] Rong Zhu, Zhaoying Zhou, Sha Li, Xuefeng Sun, A novel miniature azimuth-level detectorbased on MEMS[C],2001Microelectromechanical Systems Conference Berkeley, CA, USA,24-26August2001IEEE Circuits&Syst. Soc
    [27] Su, K.;Ren, D.H.;You, Z.;Zhou, Q, Application of MEMS-IMU/magnetometer integratedsystem on the attitude determination of micro satellite[C], Intelligent Mechatronics andAutomation,2004
    [28]李荣冰,刘建业,段方,牛新元,MEMS-INS微型飞行器姿态确定系统的实现研究[J],应用科学学报,2006,24(6):618-622
    [29] A. Kourepenis, J. Borenstein, J. Connelly, R. Elliott, P. Ward, M. Weinberg, Performance OfMems Inertial Sensors[C], IEEE PLANS,1998, California
    [30] Pascal Zwahlen, Anne-Marie Nguyen, Yufeng Dong et al. Navigation grade MEMSaccelerometer[C], IEEE MEMS,2010,631-634..
    [31] Byoung-doo Choi,Sangjun Park,Hyoungho Ko, The first sub-deg/HR bias stability,silicon-microfabricated gyroscope[C], The13th International Conference on Solid-StateSensors, Actuators and Microsystems,2005
    [32]王忆锋, MEMS和导航级MEMS惯性传感器[J],飞航导弹,2008,8:53-56
    [33] K I Lee, H Takao,K Sawada, et al., Analysis and experimental verification of thermal drift ina constant temperature control type three-axis accelerometer for high temperatures with anovel composition of Wheatstone bridge[C], Micro Electro Mechanical Systems, MEMS2004:241-244.
    [34] M I Ferguson, D Keymeulen, C Peay, et al., Effect of Temperature on MEMS VibratoryRate Gyroscopeence[C],1-6.
    [35]陈怀,张嵘,周斌,陈志勇,微机械陀螺仪温度特性及补偿算法研究[J],传感器技术,2004,10:25-27
    [36] FANG JC, LI JL, Integrated Model and Compensation of Thermal Errors of SiliconMicroelectromechanical Gyroscope [J], IEEE Transactions On Instrumentation AndMeasurement, VOL.58, NO.9,2009:2923-2930
    [37]周斌,微机械陀螺数字化技术研究[D],清华大学博士论文,2003.
    [38] CRG20Digital Angular Rate Sensor Datasheet, CRG20-00-0100-110Rev8, DCR617538
    [39]任大海,顾启泰,毛刚,尤政,刘学斌,微型惯性测量组合标定技术[J],清华大学学报(自然科学版),2001.8
    [40] He Xiufeng, Liu Jianye, Analysis Of Lever Arm Effects In GPS/IMU Integration System[J],Transactions of Nanjing University of Aeronautics&Astronau,2002,1
    [41] Shin Eun-Hwan. Estimation Techniques for Low-Cost Inertial Navigation[D]. Calgary,Alberta: Doctor Dissertation, MAY,2005
    [42]孙宏伟,房建成,盛蔚,一种基于MEMS的微惯性测量单元标定补偿方法[J],北京航空航天大学学报,第34卷第4期,2008:439-442.
    [43]尹文,MEMS-IMU微惯性测量单元误差建模与补偿技术[D],国防科学技术大学硕士论文,2007
    [44]李杰,洪惠惠,张文栋, MEMS微惯性测量组合标定技术研究[J],传感技术学报,21卷7期,2008:85-89.
    [45]罗兵,吴美平,尹文,曹聚亮,微机械陀螺温度系数的快速标定方法[J],传感技术学报,2010,10
    [46] Shin E-H and El-Sheimy N2002A new calibration method forstrapdown inertial navigationsystems[J]. Z. Vermess.1271-10
    [47] Z F Syed, P Aggarwal, C Goodall, et al., A new multi-position calibration method forMEMS inertial navigation systems[J], Meas. Sci. Technol.18(2007)1897-1907.
    [48] Tsuchiya Toshiyuki, Funabashi Hirofumi, A z-axis differential capacitive SOI accelerometerwith vertical comb electrodes[J], Sensors and Actuators A116(2004)378-383
    [49] Burrer Chr, Esteve J, A novel resonant silicon accelerometer in bulk-micromachiningtechnology[J], Sensors and Actuators A46-47(1995)185-189
    [50] Hierold C, N her U, Scheiter T, Mensching B, Steger M, Hildebrandt A, Tielert, R, A pureCMOS surface-micromachined integrated accelerometer[J], Sensors and Actuators A57(1996)111-116
    [51] Chih-ming Sun, ChuanWei Wang, Weileun Fang, On the sensitivity improvement of CMOScapacitive accelerometer[J], Sensors and Actuators A141(2008)347-352
    [52] Chau, K. H.-L., et al.. An integrated force-balanced capacitive accelerometer for low-gapplications[J]. Sensors and Actuators A,54,(1996)472-476
    [53] T. Meydan. Recent trends in linear and angular accelerometers [J]. Sensors and actuator A59,1997:43-50.
    [54]胡雪梅.微机械电容式加速度计的系统分析[D].西安电子科技大学,2006.3.
    [55] Bernstein J. An overview of MEMS inertial sensing technology[J]. Sensors,2003,(2)
    [56]李疆,高钟毓,董景新.微加速度计摆片质量的测量方法研究[J],机械工程学报,2004,40(3):115-118.
    [57]何晓平,张德,文贵印.静电伺服微加速度计的量程设计[J].中国惯性技术学报.2003,11(6):90-93.
    [58] B. S. Davis, Using low-cost mems accelerometers and gyroscopes as strapdown IMUs onrolling projectiles[C], PLANS, Apr.1998, pp.594-601.
    [59]吴天准,董景新,刘云峰.电容式微机械加速度计闭环系统的零偏[J].清华大学学报(自然科学版),2005,45(2):201-204.
    [60]吴天准,董景新,赵长德,冯宇翔.梳齿式微机械加速度计闭环系统的线性度[J].清华大学学报(自然科学版),2006,45(2):201-204.
    [61] Gang Dai, Mei Li, Wei Su and Beibei Shao. Systematic Error Analysis and Modeling ofMEMS Close-loop Capacitive Accelerometer [C]. Proceedings of the5th IEEE InternationalConference on Nano/Micro Egineering and Molecular Systems, Xiamen, P. R. China,NEMS2010, January20-23,2010:92-96.
    [62]刘云峰,董景新.静电力反馈微加速度计的吸合失效[J].机械工程学报,2006,42(12):128-131.
    [63] Wei Tech Ang,Pradeep K. Khosla, and Cameron N. Riviere, Nonlinear Regression Model ofa Low-g MEMS Accelerometer[J], IEEE Sensors Journal,vol.7,no.1,Jan2007,pp81-86
    [64]翁彦雯,董景新,刘云峰,微机械加速度计的温度特性实验研究[J],航空精密制造技术,2004,4:40-43
    [65]符强,马炳和,常洪龙,姿态参照系统中硅微机械加速度计的温度补偿[J],仪表技术与传感器,2009,1:15-17
    [66] F. Mohd-Yasin, N. Zaiyadi, D.J. Nagel, D.S. Ong, C.E. Korman, A.R. Faidz, Noise andreliability measurement of a three-axis micro-accelerometer[J], Microelectronic Engineering,86(2009):991-995
    [67] Single-Axis, High-g, iMEMS Accelerometers ASXL78Datasheet,http://www.analog.com/static/imported-files/data_sheets/ADXL78.pdf
    [68] www.colibrys.com
    [69] www.silicondesigns.com
    [70] www.summitinstruments.com
    [71] Lee, Kyung Il, Takao Hidekuni, Sawada Kazuaki, Ishida Makoto, Low temperaturedependence three-axis accelerometer for high temperature environments with temperaturecontrol of SOI piezoresistors[J], Sensors and Actuators A,104(2003),53-60
    [72] Christian Falconi, Marco Fratini, CMOS microsystems temperature control[J], Sensors andActuators B129(2008)59-66
    [73] H. Lakdawala and G.K. Fedder, Temperature stabilization of CMOS capacitiveaccelerometers [J], J. Micromech. Microeng.14(2004):559-566
    [74] Charles Kitchin, Make Wide Temperature Range, Ultralow Drift Accelerometers Using LowCost Crystal Ovens, AN-385, Application note, Analog Device.
    [75] Hyoungho Ko and Dong-il Dan Cho, Highly programmable temperature compensatedreadout circuit for capacitive microaccelerometer[J], Sensors and Actuators A,158,(2010):72-83
    [76] Lin He, Yong-Ping Xu, Anping Qiu, Folded silicon resonant accelerometer with temperaturecompensation[C], Sensors,2004. Proceedings of IEEE,512-515
    [77]王成伟,闫桂珍,朱泳.ICP硅深槽刻蚀中的线宽控制问题研究[J].微纳电子技术,2003,7-8:104-107.
    [78] Xuesong Liu, Chengwei Wang, Yong Zhu, Guizhen Yan. Vertical profiles and CD losscontrol in deep RIE technology[C].7th International Conference on Solid-State andIntegrated Circuits Technology Proceedings, ICSICT2004.2004:1848-1851.
    [79]张鉴,黄庆安,李伟华, MEMS工艺中反应离子深刻蚀硅片的数值模型研究[J],传感技术学报,2006,5:130-133
    [80]李伟剑.微机电系统的多域耦合分析与多学科设计优化[D].西北工业大学博士论文.2004:10-11.
    [81] Michael A. Gray,Multiphysics Technologies,Computing in Science andEngineering[J],March2008, vol.10(2):82-84
    [82] S. Rajendran a, K.M. Liew, Design and simulation of an angular-rate vibratingmicrogyroscope[J], Sensors and Actuators A116(2004)241-256
    [83] Mrigank Sharma, Design And Modeling Of Advanced Gyroscopes[D], Master Degree thesisof the University of British Columbia,2006
    [84] Said Emre Alper, Kivanc Azgin, Tayfun Akin, A high-performance silicon-on-insulatorMEMS gyroscope operating at atmospheric pressure[J], Sensors and Actuators A135(2007)34-42
    [85]唐海林,刘显学,周浩等,双框架解耦微陀螺机理样机研制[J],信息与电子工程,2010,8(6):716-719
    [86] Huikai Xie,Fedder, G.K., Fabrication, characterization, and analysis of a DRIECMOS-MEMS gyroscope[J], IEEE Sensors Journal,2003,3(5):622-631
    [87] Ajit Sharma, M. Faisal Zaman, Farrokh Ayazi, A0.2/hr Micro-Gyroscope with AutomaticCMOS Mode Matching[C],2007IEEE International Solid-State Circuits Conference,2007
    [88] W.Geiger, W. U. Butt, A. Gaiber et al. Decoupled microgyros and the design principleDAVED[J], Sensors and actuators A95(2002)239-249
    [89] Said Emre Alper, Tayfun Akin, A symmetric surface micromachined gyroscope withdecoupled oscillation modes[J], Sensors and Actuators A97-98(2002)347-358
    [90] Precision Navigation and Pointing Gyroscope CRM200Technical Datasheet,http://www.pinpoint-gyro.com/wp-content/uploads/2011/03/PinPoint-CRM200-00-0100-132-Rev-4.pdf
    [91] High Stability, Low Noise Vibration Rejecting Yaw Rate Gyro ADXRS646Datasheet,http://www.analog.com/static/imported-files/data_sheets/ADXRS646.pdf
    [92] STIM202Multi-Axis Gyro Module Datasheet,http://www.sensonor.com/media/70534/datasheet_stim202_ts1439_r8.pdf
    [93] David H. Titterton, John L.Weston, Strapdown Inertial Navigation Technology2ndEdition[M], Published by IEE,2004
    [94]刘梅,周百令,硅微陀螺机械热噪声研究[J],仪器仪表学报,2006, S2:145-146
    [95] Wang Hao, Weifeng Tian, Modeling the Random Drift of Micro-Machined Gyroscope withNeural Network[J], Neural Processing Letters(2005)22:235-247
    [96]孙畅,富立,刘文丽,成竹,微陀螺零偏在线补偿方法及仿真[J],弹箭与制导学报,2008,3:58-60
    [97]李建利,房建成,改进的MEMS陀螺静态误差模型及标定方法[J],宇航学报,2007,28(6):1614-1618
    [98]郑翔,杨福俊,裘安萍,何小元,双线性振动微陀螺特征频率光学测试和有限元分析[J],工程力学,2009,26(3):195-201
    [99] Alfredo Cigada, Elisabetta Leo, Marcello Vanali, Electrical method to measure the dynamicbehaviour and the quadrature error of a MEMS gyroscope sensor[J], Sensors and ActuatorsA134(2007)88-97
    [100] Sitaraman V.Iyer, Modeling and simulation of non-idealities in a z-axis cmos-memsgyroscope[D], PhD thesis of Caernegie Mellon University,2003
    [101]施芹,提高硅微陀螺仪性能若干关键技术研究[D],东南大学博士论文,2005
    [102] Y.Y. Bao, C.L. Yung, Modelling and compensation of quadrature error for silicon MEMSmicrogyroscope[C], Proceedings of Power Electronics and Drive Systems,2001. Vol.2:871-876
    [103] Nan-Chyuan Tsai, Chung-Yang Sue, Chih-Che Lin, Design and dynamics of an innovativemicro gyroscope against coupling effects[J], Microsyst Technol(2008)14:295-306
    [104]施芹,硅微陀螺仪的机械耦合误差分析[J],光学精密工程,2008.10(5):894-898
    [105] Tao Jiang, Anlin Wang, Jiwei Jiao, Guangjun Liu, Detection capacitance analysis methodfor tuning fork micromachined gyroscope based on elastic body model[J], Sensors andActuators A128(2006)52-59
    [106]杨军,高钟毓,张嵘,陈志勇,周斌,线振动硅微机械陀螺结构误差参数分离和辨识[J],中国惯性技术学报,2007,15(3):327-333
    [107] Shcheglov, K., Evans, C., Gutierrez, R., Tang, T.K., Temperature dependent characteristicsof the JPL silicon MEMS gyroscope[C],2000IEEE Aerospace Conference Proceedings, Vol.1,403-411
    [108] A.Kulygin, C. Kirsch, P. Schwarz, U. Schmid, H. Seidel, Study on the lower resolution limitand the temperature-dependent performance of a surface micromachined gyroscope[C],MEMS2008, Tucson, AZ, USA, January13-17,2008
    [109] Guangjun Liu, Anlin Wang, Tao Jiang, Jiwei Jiao,Jong-Bok Jang, Effects of environmentaltemperature on the performance of a micromachined gyroscope[J], Microsyst Technol(2008)14:199-204
    [110] C C Painter and A M Shkel, Structural and thermal modeling of a z-axis rate integratinggyroscope[J], J. Micromech. Microeng.13(2003)229-237
    [111]赵晓辉,伊国兴,王常虹,硅微陀螺温度漂移补偿研究[J],传感器与微系统,2008,27(10):48-50
    [112]罗兵,尹文,江明明,吴美平,微机械陀螺零偏温度梯度建模[J],电光与控制,2008,15(10):29-31
    [113]高世桥,刘海鹏.微机电系统力学[M].国防工业出版社,2008
    [114]方靖,顾启泰,丁天怀,捷联式微惯性导航系统的标定方法[J],弹箭与制导学报,2007,27(5):29-32
    [115] Wasilik M, Chen N, Deep reactive ion etch conditioning recipe[C], Micromachining andMicrofabrication Process Technology IX, Proceeding of SPIE, Vol5342,(2004):103-110
    [116] Comsol Material Database. http://www.comsol.com
    [117]李童杰,刘云峰,董景新,刘继奎,微加速度计温度特性及敏感元件自恒温方案[J],清华大学学报(自然科学版),2010,50(7):1013-1017
    [118] A. Wolfenden and A.C. Anthony, Measurement Of Young’S Modulus And Damping ForHexoloy Sg Silicon Carbide[J], Young’s Modulus And Damping,1999,41(6):611-615
    [119] Sitaraman Iyer and Tamal Mukherjee. Simulation of Manufacturing Variations in a Z-axisCMOS-MEMS Gyroscope[C], Technical Proceedings of the2002International Conferenceon Modeling and Simulation of Microsystems,2002, Vol.1:186-189

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

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

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