用户名: 密码: 验证码:
基于脉搏波分析的麻醉手术中无创生理参数监测系统
详细信息    本馆镜像全文|  推荐本文 |  |   获取CNKI官网全文
摘要
麻醉手术中的低血压及伤害感受的监测及预测是非常重要的,它可以帮助麻醉医生及时了解病人生理状态的变化。目前市场上已有的麻醉监测设备大多只能提供简单的测量参数显示和报警功能,这些设备中的参数虽多,但却没有建立起这些参数与生理状态之间的联系,也就难以发挥麻醉医生们的主观能动性。随着近几十年来对光电容积脉搏波(PPG)所做的大量临床研究,PPG已被广泛的用于麻醉生理参数监测,该方法安全无创、可靠性高,所以本论文根据PPG研究的生理学基础,设计了基于PPG的麻醉手术中无创生理参数监测系统。
     本论文的主要研究工作如下:
     1.设计了基于光电容积脉搏波的麻醉手术中无创生理参数监测系统,可以完成麻醉过程相关的心率(HR)、血氧饱和度(SpO2)、收缩压、舒张压、平均压、脉搏波传导时间(PTT)、衰减时间常数(RC)、光电容积脉搏波幅度(PPGA)、下降支面积比(ARD)等参数的同步实时监测。其中,RC和ARD均为本系统所独有的生理参数。
     2.设计实现了以nRF24L01模块为核心的低功耗无线传输方案。该方案以nRF系列芯片内嵌的Enhanced ShockBurst协议为基础,设计开发了短距离无线传输网络协议(NWSN),实现了数据传输的低功耗。该方案与载波监听多路访问协议(CSMA)相比,极大地降低了数据收发的电流,总功耗约为后者的1/3。同时该方案利用了无线天然无漏电流、无需隔离的特点,节省了硬件成本,降低了设计复杂度。
     3.提出了麻醉监测过程中脉搏波传导时间变异性(PTTV)的分析算法,并首次将其应用于麻醉手术中无创生理参数监测系统。本论文提出了一个新的光电容积脉搏波(PPG)信号特征点检测算法。此算法主要采用窗口化的局部查找方法,再与一阶导数和逻辑判断相结合来检测特征点。当接受区间为24.0ms时,该算法的检出率高于95%。由心电图(ECG)信号的R波和PPG信号的P波位置可以计算得到PTT,进而将其作为PTTV的原始输入,采用自回归(AR)模型建立PTTV分析算法。
     4.首次将PTTV分析方法应用于麻醉低血压研究。在浙江大学医学院附属妇产科医院,经过24例临床实验,结果表明,甚低频与低频功率的比值(VLF/LF)以及低频功率与总功率的比值(LF/ALL)这两个参数的P值分别为:PVLf/LF=0.028、PLF/ALL=0.031,说明这两个参数在低血压组与正常组中有明显差异,具有预测麻醉低血压发生的可行性。
     该系统被应用于麻醉低血压及伤害感受的临床研究,初步实现了对低血压及伤害感受的有效评估。本论文所设计之系统严格按照国家相关标准要求开发,具有很好的临床应用前景。
Hypotension and nociceptive monitoring in anesthesia is very important, which can help the anesthesiologist to keep abreast of changes in the patient's physiological state. Anesthesia monitoring equipments on the market are mostly only able to provide some simple parameters and alarm functions, which did not establish the link between parameters and physiological states. So it is difficult to play the subjective initiative of anesthesia doctors. In recent decades, a large number of clinical studies about photoplethysmography (PPG) were done, PPG has been extensive used in anesthesia physiological parameters monitoring. Based on these studies, a non-invasive physiological parameters monitoring system is designed.
     The main research works are as follows:
     A. A non-invasive physiological parameters monitoring system during anesthesia based on photoplethysmography is designed, which can monitor heart rate, oxygen saturation, systolic pressure, diastolic pressure, mean pressure, pulse wave conduction time, decay time constant. PPG amplitude, area ratio of descending limb and other parameters in the same time. Especially, decay time constant and area ratio of descending limb are the unique physiological parameters in this system.
     B. A low power wireless transmission scheme based on nRF24L01module is designed. The scheme uses a wireless network protocol base on Enhanced ShockBurst. which is an on-chip Agreement of nRF series, to achieve low-power data transmission. The scheme greatly reduceed the current of the data sent and received. The total power consumption is about1/3of the carrier sense multiple access protocol (CSMA). The same time, because there is no leak current through wireless, so the scheme saved the hardware costs and reduced the design complexity.
     C. A algorithm about monitoring pulse transmit time variability (PTTV) during anesthesia is proposed, which is applied to non-invasive physiological parameters monitoring system during anesthesia at the first time. Meanwhile, a new PPG signal feature detection algorithm is also proposed. This algorithm is the combination of the local search method in main window and first order derivative with logical judgment. When the acceptance interval is24.0ms. the detection rate of this algorithm is higher than95%. Pulse transmit time(PTT) can be calculated by the P wave position of PPG signal and the R wave position of electrocardiogram (ECG) signal, which is the original input of PTTV analysis algorithm. Finally, an auto regressive model is used in this algorithm. The algorithm has a reliable performance after the experiments.
     D. The PTTV analysis method is the first time used in hypotension research during anesthesia.24cases of clinical trials have been done in Women's Hospital School of Medicine Zhejiang University. The results showed that the P value of the ratio of very low frequency power and low frequency power (VLF/LF) and the ratio of low-frequency power and total power (LF ALL) were:PVLF/LF=0.028. PLF/ALL=0.031. which indicates that these parameters have feasibility of anesthesia hypotension prediction.
     The system was used in clinical research about hypotension and nociceptive during anesthesia, which realized the effective assessment of hypotension and nociceptive preliminarily. The system is strictly in accordance with the relevant national standards, and it also has high reliability and good prospects for clinical application.
引文
[1]于布为.2009年麻醉科的回顾与展望[J].上海医学.2010,33(1):24-25.
    [2]龚忠兵,王保华,刘红等.麻醉监护[J].中国医疗器械信息.2003,9(4):9-13.
    [3]孙家骧.现代麻醉监测技术和发展[J].引进国外医药技术与设备.1998,6:26-31.
    [4]丁蔚萍.医用监护仪的发展近况[J].科技信息(吉林).1995,8:10.
    [5]陈利佳.麻醉信息系统的实现及意义[J].重庆医学.2008,21:1269-1271.
    [6]Mindray. iPM9800[Z].
    [7]Comen. STAR8000C[Z].
    [8]Goldway. UT4000F Pro[Z].
    [9]陈宁,李邦翅.临床麻醉仪器设备新进展[J].医疗卫生装备.2012.33(2):78-80.
    [10]冯成才,闵红星.综合措施防治剖宫产术中低血压的临床观察[J].宁夏医科大学学报.2009,31(5):633-635.
    [11]王德智.仰卧位低血压综合征的诊断与防治[J].中国实用妇科与产科杂志.1995.11(1):23.
    [12]袁桂莲,许我先.妇产科综合征[G].北京:人民卫生出版社,1989:229.
    [13]曹缵孙,陈晓燕.妇产科综合征[G].北京:人民卫生出版社,2003:249-251.
    [14]邓晓勤,陈江辉.腰-硬联合麻醉下剖宫产手术低血压原因分析[J].医学信息.2010,5:1106-1107.
    [15]邓晓勤,陈江辉.腰-硬联合麻醉下剖宫产手术低血压原因分析[J].医学信息.2010(5):1106-1107.
    [16]Mccrae A F. Wildsmith J A. Prevention and treatment of hypotension during central neural block[J]. Br J Anaesth.1993,70:672.
    [17]Brooker R F, Butterworth J F T, Kitzman D W. et al. Treatment of hypotension after hyperbaric tetracaine spinal anesthesia, a randomized, double-blind, cross-over comparison of phenylephrine and epinephrine[J]. Anesthesiology,1997.86:797.
    [18]胡祥芹,王春梅.有创血压监测的研究进展[J].护理研究.2008.22(1):193-195.
    [19]张政波:吴太虎.无创血压测量技术与进展[J].中国医疗器械杂志.2003.27(3):196-199.
    [20]Willam K. Etal. Apparatus and method for noninvasive blood pressure measurement[J], United States Patent.1996,5533511.
    [21]向海燕,俞梦孙,刘娟.用有创血压验证示波法测量血压的准确性[J].北京生物医学工程.2005.24(3):195-198.
    [22]汤池,杨国胜,席涛.无创连续血压测量方法的研究进展[J].医疗卫生装备.2004.25(10):26-27.
    [23]Heard S O, Lisbon A. Toth I. et al. An evaluation of a new continuous blood pressure monitoring system in critically ill patients[J]. Journal of Clinical Anesthesia.2000.12(7): 509-518.
    [24]Maguire M. Ward T. Markham C. et al. A comparative study in the use of brachial photoplethysmography and the QRS complex as timing references in determination of pulse transit time[J]. Engineering in Medicine and Biology Society.2001. Proceedings of the 23rd Annual International Conference of the IEEE.2001.1:215-218.
    [25]Jago J R, Murray A. Repeatability of peripheral pulse measurements on ears, fingers and toes using photoelectric plethysmography[J]. Clin Phys Physiol Meas.1988.9(4):319-330.
    [26]杜娟.生理信号变异性分析及其在高血压发展中的应用研究[D].浙江大学.2006.
    [27]Chamchad D. Arkoosh V A. Horrow J C. et al. Using heart rate variability to stratify risk of obstetric patients undergoing spinal anesthesia[J]. Anesth Analg.2004.99:1818-1821.
    [28]Hanss R. Bein B. Ledowski T. et al. Heart rate variability predicts severe hypotension after spinal anesthesia for elective cesarean delivery [J]. Anesthesiology.2005,102:1086-1093.
    [29]Hanss R, Bein B. Weseloh H. et al. Heart rate variability predicts severe hypotension after spinal anesthesia[J]. Anesthesiology.2006.104:537-545.
    [30]Huang C J, Kuok C H. Kuo T B. et al. Pre-operative measurement of heart rate variability predicts hypotension during general anesthesia[J]. Acta Anaesthesiol Scand.2006.50:542-548.
    [31]Fujiwara Y. Kurokawa S. Asakura Y. et al. Correlation between heart rate variability and haemodynamic fluctuation during induction of general anaesthesia:comparison between linear and non-linear analysis[J]. Anaesthesia.2007.62:117-121.
    [32]Fujiwara Y. Ito H, Asakura Y, et al. Preoperative ultra short-term entropy predicts arterial blood pressure fluctuation during the induction of anesthesia[J]. Anesth Analg.2007,104: 853-856.
    [33]Guignard B. Monitoring analgesia[J]. Best Practice & Research Clinical Anaesthesiology. 2006,20(1):161-180.
    [34]Hoymork S C. Antinociceptive monitors:tools or fools? [J]. Acta Anaesthesiologica Scandinavica.2008,52:1035-1037.
    [35]庄心良.曾因明,陈伯銮.现代麻醉学[M].北京:人民卫生出版社,2003.
    [36]Rantanen M. Yli-Hankala A. Van G M. et al. Novel multiparameter approach for measurement of nociception at skin incision during general anaesthesia[J]. British Journal of Anaesthesia.2006.96:367-376.
    [37]Rundshagen I. Schrtoder T. Prichep L S. et al. Changes in cortical activity during inducti-on of anaesthesia with thiopental/fentanyl and tracheal intubation:a quantitative electroencep-halographic analysis[J]. British Journal of Anaesthesia.2004.92:33-38.
    [38]Kawaquchi M, Takamatsu I. Masui K. et al. Effect of landiolol on bispectral index and spectral entropy responses to tracheal intubation during propofol anaesthesia[J]. British Journal of Anaesthesia.2008.101:273-278.
    [39]Leslie K, Sessler D I. Smith W D. et al. Prediction of movement during propofol/nitrous oxide anesthesia:performance of concentration, electroencephalographic, papillary, and hemodynamic indicators[J]. Anesthesiology.1996.84:52-63.
    [40]Dincklage F V, Send K. Hackbarth M. et al. Comparison of the nociceptive flexion reflex threshold and the bispectral index as monitors of movement responses to noxious stimuli under propofol mono-anaesthesia[J]. British Journal of Anaesthesia.2009,102:244-250.
    [41]Dincklage F V, Velten H. Rehberg B. et al. Monitoring of the responsiveness to noxious stimuli during sevoflurane mono-anaesthesia by using RIII reflex threshold and bispectral index[J], British Journal of Anaesthesia.2010,104:740-745.
    [42]Storm H. Changes in skin conductance as a tool to monitor nociceptive stimulation and pain[J]. Current Opinion in Anaesthesiology.2008.21:796-804.
    [43]Ledowski T, Ang B. Schmarbeck T. et al. Monitoring of sympathetic tone to assess postoperative pain:skin conductance vs surgical stress index[J]. Anaesthesia,2009.64: 727-731.
    [44]Loggia M L. Juneau M, Bushnell M C. Autonomic responses to heat pain:Heart rate, skin conductance, and their relation to verbal ratings and stimulus intensity[J]. Pain.2011,152: 592-598.
    [45]Ledowski T, Pascoe E, Ang B, et al. Monitoring of intra-operative nociception:skin conductance and surgical stress index versus stress hormone plasma levels[J]. Anaesthesia.2010, 65:1001-1006.
    [46]Kushki A. Fairley J. Merja S, et al. Comparison of blood volume pulse and skin conductance responses to mental and affective stimuli at different anatomical sites[J]. Physiological Measurement.2011.32:1529.
    [47]Nijboer J A. Dorlas J C. Comparison of plethysmograms taken from finger and pinna during anaesthesia[J]. British Journal of Anaesthesia.1985.57:531-534.
    [48]Luginbuhl M, Reichlin F. Sigurdsson G H. et al. Prediction of the haemodynamic response to tracheal intubation:comparison of laser-Doppler skin vasomotor reflex and pulse wave reflex[J]. British Journal of Anaesthesia.2002,89:389-397.
    [49]骆璇,常昌远.脉搏波的解析模型及初步临床应用研究[J].数理医药学杂志.1999.4:293-295.
    [50]杭燕南.沈坚,孙大金,无创性血氧饱和度检测及其临床意义[J].国际麻醉学与复苏杂志.1987.8(2):107-111.
    [51]张伯英.指端光电容积脉搏波评价麻醉患者循环功能稳定性的初步研究[J].生物医学工程与临床.2010.14(6):142-145.
    [52]张伯英,罗晓民.基于光电容积血流脉搏波特征参数监测术术中应激反应临床研究[J].生物医学工程与临床.2011.15(6):533-538.
    [53]王燕,张松.杨益民.et a1.基于容积血流脉搏波的心血管血流参数监护模块的研制[J].北京生物医学工程.2006,2:148-150.
    [54]罗志昌,张松,杨益民.脉搏波的工程分析与临床应用[M].北京:科学出版社.2006.
    [55]Geddes L A, Voelz M H. Babbs C F, et al. Pulse Transit-Time as an Indicator of Arterial Blood-Pressure[J]. Psychophysiology.1981,18(1):71-74.
    [56]Chen W, Kobayashi T. Ichikawa S, et al. Continuous estimation of systolic blood pressure using the pulse arrival time and intermittent calibration[J]. Medical & Biological Engineering & Computing.2000,38(5):569-574.
    [57]Lin Y H, Ko PCI. Wang H Y, et al. Estimation of beat-to-beat systolic blood pressure using pulse arrive time and pulse width derived from the photoplethysmogram[J]. Annual International Conference of the IEEE Engineering in Medicine and Biology Society IEEE Engineering in Medicine and Biology Society Conference.2004.5:3456-3458.
    [58]向海燕,俞梦孙.用脉搏波传导时间实现血压连续测量[J].医疗卫生装备.2006,27(2):19-21.
    [59]Chen W, Kobayashi T, Iehikawa S. et al. Continuous estimation of systolie blood pressure using the Pulse arrival time and intermittent calibration [J]. Med Boil Eng Comput.2000,38(5):569-574.
    [60]Jago J R. Murray A. Repeatability of peripheral pulse measurements on ears, fingers and toes using photoelectric plethysmography[J]. Clinical Physics and Physiological Measurement. 1988.9:319.
    [61]Li L, Wang Z Z. Study on interval variability of arterial pulse [J]. Proc. of the First Joint BMES/EMBS Conference:Atlanta.1999:223.
    [62]Patterson R. Belalcazar A. RR and R wave to finger pulse interval changes as a function of respiration rate and body position using VCRS [J]. Proceedings of the 19th Annual International Conference of the IEEE Chicago.1997:2190-2192.
    [63]全晓莉,廖磊.秦浪.外周阻力的测量与临床研究[J].中国医学物理学杂志.2008.25(2):598-600.
    [64]高磊,李卫华,动脉顺应性的研究进展[J].心血管病学进展.2007.28(1):74-76.
    [65]孙磊.血管系统仿真建模与脉搏波分析研究[D].杭州:浙江大学.2008.
    [66]周毅,苏镇培.生理系统建模仿真的研究与展望[J].现代医学仪器与应用.2002.14(1):12-15.
    [67]柳兆荣,李惜惜.弹性腔理论及其在心血管系统分析中的应用[M].北京:科学出版社,1987.
    [68]De W R B P, Schreuder J J. Den Van B P C M, et al. An evaluation of cardiac output by five arterial pulse contour techniques during cardiac surgery [J].62.2007.8(760-768).
    [69]European Parliament And Of Council. COUNCIL DIRECTIVE 93/42/EEC [S].2007.
    [70]张坤,曹鸣.一种基于小波变换的心电去噪算法[J]现代生物医学进展.2009,19:44-47.
    [71]刘宝华,杜连涛,王国永.脉搏波中呼吸干扰的去除[J].燕山大学学报.2008.32(1):53-57.
    [72]张翔,周严.基于自适应滤波的心电图中呼吸干扰的抑制[J].生物医学工程学杂志2011.28(1):33-36.
    [73]郭艳卫.谢建华.刘海波.单电源心电图机前置放大器的设计[J].吉林大学学报:信息科学版.2011.2:116-120.
    [74]管青华.浅谈医疗设备电气安全[J].医疗设备信息.2006.21(3):76-77.
    [75]崔志刚,刘坚,高虹.医疗设备电气安全的检测与探讨[J]医疗卫生装备.2010.31(5):110-113.
    [76]The International Electrotechnical Commission. IEC 60601-1:1988 Medical electrical equipment-Part 1:General requirements for safety[S].1988.
    [77]The International Electrotechnical Commission. IEC 60601-2-34:2000 Medical electrical equipment Part 2-34:Particular requirements for the safety, including essential performance, of invasive blood pressure monitoring equipment[S].2000.
    [78]中国国家标准化管理委员会.GB 9706.1-2007医用电气设备第1部分:安全通用要求[S].2007.
    [79]中国国家标准化管理委员会.YY 0668-2008医用电气设备第2-49部分:多参数患者监护设备安全专用要求[S].2008.
    [80]中国国家标准化管理委员会.GB 9706.25-2005医用巴气设备第2-27部分:心电监护设备安全专用要求[S].2005.
    [81]李雨明,张宜川,潘全亮.电介质强度及其测试[J].中国医疗器械杂志.2007,31(4):67-69.
    [82]徐美英.麻醉与主要器官的循环[J].国外医学麻醉学与复苏分册.1994,15(4):257-260.
    [83]谈维洁,王文清,罗惠兰,et a1.从指尖脉搏图的变化看心脏介入治疗对微循环的影响[J].17.2001,2(82-86).
    [84]闻大翔,周颖,王珊娟,et a1.体动和低灌注影响脉搏血氧饱和度准确性的临床研究[J].临床麻醉学杂志.2004,20(6):333-335.
    [85]Webster J G. Design of Pulse Oximeters[M]. Bristol:Institute of Physics Publishing,1997.
    [86]Allen J. Photoplethysmography and its application in clinical physiological measurement [J]. Physiol. Meas.2007,28:1-39.
    [87]Allen J, Murray A. Variability of photoplethysmography peripheral pulse measurements at the ears, thumbs, and toes[J]. IEEE. Proc. Sci. Meas. Technol.2000.147:403-407.
    [88]张兰保.触电电流对人体的伤害[J].电工技术.1992,7:63-64.
    [89]百度百科.串行接口http://baike.baidu.com/view/69108.htm[Z].
    [90]夏建春.标准并行接口在数据采集中的应用[J].科技创新导报.2007.32:29.
    [91]谢力.物联网中的几种短距离无线传输技术浅析[J].四川省通信学会2010年学术年会论文集.2010.
    [92]Adi.采用iCoupler技术、符合IEC 60601-1标准的数字隔离器[J].电子技术应用.2008.6:10.
    [93]兰颖,李刚.iCoupler数字隔离技术及其应用[J].电子技术应用.2006,11:89-93.
    [94]Bilstrup K. A Preliminary Study of Wireless Body Area Networks [R]. School of Information Science, Computer and Electrical Engineering, Halmstad University,2008.
    [95]Alexandros P. Nikolaos G B. A Survey on Wearable Sensor-Based Systems for Health Monitoring and Prognosis[J]. IEEE TRANSACTIONS ON SYSTEMS, MAN, AND CYBERNETICS.2010.40(1):1-12.
    [96]Marinkovic S J, Popovici E M. Spagnol C, et al. Energy-Efficient Low Duty Cycle MAC Protocol for Wireless Body Area Networks[J]. IEEE Trans. Inf. Technol. Biomed.2009,13(6): 915-925.
    [97]王芳,吴效明.基于ZigBee的无线传输脉搏检测系统设计[J].医疗卫生装备.2009.30(6):21-22.
    [98]Ying Z, Hannan X. Bluetooth-Based Sensor Networks for Remotely Monitoring the Physiological Signals of a Patient[J]. IEEE Trans. Inf. Technol. Biomed.2009,13(6): 1040-1048.
    [99]Xiuquan F, Weihong C. Shuming Y, et al. A Wireless Implantable Sensor Network System for In Vivo Monitoring of Physiological Signals[J]. IEEE Trans. Inf. Technol. Biomed.2011. 15(4):577-584.
    [100]卫红春.信息系统分析与设计[M].北京:清华大学出版社.2006.
    [101]Aboy M. Mcnames J. Thong T. et al. An automatic beat detection algorithm for pressure signals[J]. IEEE transactions on biomedical engineering.2005,52:1662-1670.
    [102]Ramsey M W. Stewart W R. Jones C J H. Real-time measurement of pulse wave velocity from arterial pressure waveforms[J]. Medical and Biological Engineering and Computing.1995. 29(4):631-645.
    [103]Glindecrantz K, Lilja H. New software QRS detector algorithm suitable for real time application with low signal-to-noise ratios[J]. J.Biomed.Eng.1988.10(3):280-284.
    [104]李向军,陈裕泉.QRS波群时频检测方法的新进展[J].国外医学:生物医学工程分册.2005.28(5):281-286.
    [105]Kinias P. Fozzard H A. Norusis M J. A real-time pressure algorithm[J]. Computers in Biology and Medicine.1981.11:211-220.
    [106]Karamanoglu M. A system for analysis of arterial blood pressure waveforms in humans [J]. Computers and Biomedical Research.1997.30:244-255.
    [107]Zong W. Heldt T, Moody G B. et al. An open-source algorithm to detect onset of arterial blood pressure pulses[J]. Computers in Cardiology.2003.30:259-262.
    [108]Li B N. Dong M C. Vai M I. On an automatic delineator for arterial blood pressure waveforms [J]. Biomedical signal processing and control.2009,5:76-81.
    [109]Jago J R. Murray A. Repeatability of peripheral pulse measurements on ears, fingers and toes using photoelectric plethysmography[J]. Clinical Physics and Physiological Measurement. 1988,9:139.
    [110]Charles D S,Fox M D. Pulse wave velocity for cardiovascular characterization[J]. Proceedings of the IEEE 1997 23rd Northeast:Durham.1997:77-78.
    [111]冷建华.傅里叶变换[M].北京:清华大学出版社,2005.
    [112]Kitney R, Fulton T, Mcdonald A, et al. Transient interactions between blood pressure, respiration and heart rate in man[J]. J Biomed Eng.1985,7:217-224.
    [113]Marple L. A new autoregressive spectrum analysis algorithm[J]. IEEE Transactions on Signal Processing.1980,28:441-454.
    [114]Hou L X. Wei M, Wang X, et al. A novel non-iterative shape method for estimating the decay time constant of the finger photoplethysmographic pulse [J]. Journal of Zhejiang University-SCIENCE A (Applied Physics & Engineering).2011,12:438-445.
    [115]Schnider T W, Minto C F, Gambus P L, et al. The influence of method of administration and covariates on the pharmacokinetics of propofol in adult volunteers [J]. Anesthesiology.1998. 88:1170-1182.
    [116]Minto C F, Schnider T W, Egan T D, et al. Influence of age and gender on the pharmacokinetics and pharmacodynamics of remifentanil[J]. Anesthesiology.1997,86:10-23.

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

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

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