用户名: 密码: 验证码:
CO_2驱油非达西渗流理论模型研究及应用
详细信息    本馆镜像全文|  推荐本文 |  |   获取CNKI官网全文
摘要
我国低渗透油田储量资源丰富,但是受目前技术条件的限制,其采收率一般只能达到20%左右。为提高低渗透油田采收率和产量,迫切需要研究新的驱油剂和驱油方式。国内外大量研究结果表明,CO2驱可以较大幅度地提高油气采收率,而且没有污染,有利于环境保护,因此采用CO2驱油开发低渗透油田具有广阔的前景。
     为准确预测不同CO2注入方案下动态指标的变化规律,并确定最优注入方案,需要进行数值模拟研究。但现有的商业化CO2驱数值模拟软件都是基于达西定律开发的,没有考虑启动压力梯度的影响,不能够完整地反映低渗透油藏CO2驱油的全过程,有必要建立考虑非达西渗流规律的CO2驱数值模拟模型。
     本文首先开展了CO2驱油室内实验,测定了CO2驱油的最小混相压力和启动压力梯度,然后通过研究低渗透油藏CO2驱油机理以及流体的非达西渗流规律,在现有组分模型的基础上,考虑油、气、水三相启动压力梯度,建立了低渗透油藏CO2驱油三维三相多组分非达西渗流模型,采用有限差分方法和IMPES方法进行了数值求解,同时开发出一套适用于低渗透油藏CO2驱油的数值模拟软件,并通过室内驱油实验验证了该模型以及自主开发软件的正确性和适用性。
     根据树101区块和树16区块的地质特征及开发现状,应用Petrel地质建模软件分别建立了这两个区块的相控地质模型,然后利用自主开发的数值模拟软件对树101区块进行了CO2驱油数值模拟研究,对树16区块进行了CO2吞吐数值模拟研究。
     树101区块CO2驱油数值模拟研究结果表明,该区块以线性注气方式为最优井网;结合经济效益评价结果,连续注入、单段塞注入和水气交替注入三种CO2注入方式中水气交替注入为最优注入方式,其最优注入参数为:气水比为1:1,段塞大小为0.03PV,日注入量为200m3(即液态CO2年注入量为0.050PV),CO2总注入量为0.3PV,后续水驱年注入量为0.12PV。该注入方式20年末采出程度为25.15%,比基础方案提高6.15%。
     树16区块的CO2吞吐数值模拟研究结果表明,以换油率和增油量作为评价指标,优选11Y67-611井最佳吞吐参数为:日注入量80m3,周期注入量100m3,焖井时间30天,吞吐2个周期,增油量和换油率分别为101.31m3和1.01m3.m-3;优选1Y69-56井最佳吞吐参数为:日注入量60m3,周期注入量135m3,焖井时间30天,吞吐3个周期,增油量和换油率分别为154.38m3和1.14m3.m-3。
     上述研究结果表明,基于C02驱油非达西渗流理论模型的数值模拟技术能够准确地模拟榆树林油田CO2驱油过程,适用于榆树林油田注气开发方案设计,对其他同类型油田应用CO2驱提高原油采收率具有重要的指导意义。
There are rich resources of low permeability oilfields in China, whose recovery can only reach 20% commonly limited by the current technical conditions. In order to improve oil recovery and production of the low permeability oilfield, new displacement agents and methods should be urgently investigated. Domestic and foreign research results indicate that CO2 flooding can greatly improve oil recovery, and has no pollution favorable to environment protection, so using CO2 flooding to develop the low permeability oilfield has wide prospect.
     In order to accurately predict the variations in dynamic indexes of different CO2 injection schemes and determine the optimal scheme, numerical simulation research should be conducted. But the existing commercial CO2 flooding numerical simulation software was developed based on the Darcy flow without considering the influence of actuating pressure gradient, which could not fully reflect the whole CO2 displacement process of low permeability reservoirs, and therefore it is necessary to establish a CO2 flooding numerical simulation model considering non-Darcy flow pattern.
     This paper firstly conducted CO2 flooding lab experiments to measure minimal miscibility pressure and actuating pressure gradient, and then on the basis of the studies on CO2 flooding mechanisms, the non-Darcy flow pattern of low permeability reservoirs and the existing component model, a three-dimensional three-phase multi-composition percolation mathematical model suitable for CO2 flooding of the low permeability reservoirs has been established, which considered the actuating pressure gradient of oil, gas and water. After that, numerical solution was brought out by finite-difference method and IMPES method, and a set of CO2 flooding numerical simulation software suitable for low permeability reservoirs was then developed. Finally correctness and applicability of the model and developed software both were verified by core experiment.
     According to geological features and present development situation of Shu101 block and Shu16 blcok, facies-controlled geological models of the two blcoks were established by Petrel geological modeling software. And then using the self-developed numarical simulation software, this paper conducted numerical simulation research on CO2 flooding of Shu101 block and CO2 huff and puff of Shu 16 block.
     The results of numerical simulation research on CO2 flooding of Shu101 block showed that linear gas injection mode between the two designed well pattern adjustment schemes. Combained with the economic benefit evaluation results, water alternating gas injection mode among the three CO2 injection modes of continuous gas injection, single slug injection and water alternating gas injection has been confirmed to be the best, whose optimum parameters were that 0.3PV of CO2 injection quantity,0.03PV of slug size,1-1 of gas-water ratio,200m3 of liquid CO2 injection rate (i.e.0.050PV of liquid CO2 annual injection rate), and 0.12PV of subsequent water annual injection rate. The recovery degree of this scheme could reach 25.15 percent in 20 years, which could increase 6.15 percent than the foundation scheme.
     The results of numerical simulation research on CO2 huff and puff of Shu 16 block showed that taking oil-draining rate and oil increment as evaluation indexes, the optimum huff and puff parameters of well 11Y67-611 were that80 m3/d of injection rate,100 m3 of cycle injection volume,30 days of soak time and two cycle of puff and huff, while the oil increment and the oil-draining rate were separately 101.31 m3 and 1.01 m3-m-3. The optimum huff and puff parameters of well 11Y69-56 were that 60 m3/d of injection rate,135 m3 of cycle injection volume,30 days of soak time and three cycle of puff and huff, while the oil increment and the oil-draining rate were separately 154.38 m3 and 1.14 m3·m-3.
     The above reaserch results showed that the numerical simulaton technique based on the non-Darcy flow theorical model of CO2 flooding could simulate CO2 displacement process in Yushulin oilfield accurately and was suitable for designing the development scheme of Yushulin oilfield, which has important guiding significance for other same type oilfields to improve oil recovery by CO2 flooding.
引文
[1]王多伟,郝瑞芬,田殿龙.低渗透油田注气提高采收率研究现状[J].中国高新技术企业,2000,9:46-47.
    [2]李士伦,周守信,杜建芬,等.国内外注气提高石油采收率技术回顾与展望[J].油气地质与采收率,2002,9(2):236-251.
    [3]王建业,陈道学,杨汉林.特低渗透油藏高效开发技术[J].内蒙古石油化工,2003,(29):119-120.
    [4]李士伦,郭平,孙雷.改变单一的注水开发模式,因地制宜发展注气提高采收率技术[C]//中国石油天然气集团公司油气田开发工作会议论文集(1999上册).北京:石油工业出版社,1999.
    [5]李中锋,何顺利.低渗透储层非达西渗流机理探讨[J].特种油气藏,2005,12(2):35-38.
    [6]张光明,廖锐全,徐永高,等.特低渗透油藏流体渗流特征试验研究[J].江汉石油学院学报,2004,26(1):88-89.
    [7]汪伟英,张磊,杨文新,等.低渗透油藏非线性渗流规律研究[J].新疆石油天然气,2008,4(3):53-55.
    [8]孙志刚.低渗透砂岩储层启动压力及水驱油效率影响因素实验研究[D].西北大学,2005.
    [9]兰玉波,肖鲁川,赵玉珍.描述特低渗透油层两相渗流特征的改进方法[J].大庆石油地质与开发,2002,21(2):65-67.
    [10]李涛.特低渗油藏水驱油规律研究[D].中国石油大学,2008.
    [11]Zhang Ji-Cheng, Song Kao-Ping, Liu Li, et al. Investigation on Mechanisms of Polymer Enhanced Oil Recovery by Nuclear Magnetic Resonance and Microscopic Theoretical Analysis [J]. Chinese Physics Letters, May,2008,25(5):1750-1752.
    [12]Zhang Ji-Cheng, Song Kao-Ping. Considerable potential remains after Daqing polymer flood [J]. Oil & Gas Journal, Sept.15,2008,106(35):61-66.
    [13]Zhang Ji-Cheng, Shu Fang-Chun. Gas injection enhances recovery from polymer flood pilot[C]. Oil & Gas Journal, May 21,2007,105(19):45-47.
    [14]Zhang Ji-Cheng, Liu Li, Song Kao-Ping. Neural Approach for Calculating Permeability of Porous Medium [C]. Chinese Physics Letters,2006.4,23(4):1009-1011.
    [15]Zhang Ji-Cheng, Song Kao-ping, Du Ruo-Xia. Dongguo casing failures derive from many causes. Oil & Gas Journal, Jan.21,2008,106(3):51-54.
    [16]吴景春,袁满,张继成,等.大庆东部低渗透油藏单相流体低速非达西渗流特征[J].大庆石油学院学报,1999,23(2):83-84.
    [17]陈永敏.低速非达西渗流现象的实验论证[J].重庆大学学报:自然科学版,200023(增刊):59-60.
    [18]阎庆来.单相均质液体低速渗流机理及流动规律[C].第二届全国流体力学学术会议论文集.北京:科学出版社,1983.
    [19]阎庆来,何秋轩,尉立岗,等.低渗透油层中单相液体渗流特征的实验研究[J].西安石油学院院报,1990,5(2):1-6.
    [20]姚约东,葛家理.低渗透油层非达西渗流规律的研究[J].新疆石油地质,2000,21(3):213-216.
    [21]杨正明,尚根华,迟砾,等.计算启动压力梯度的经验公式[J].试采技术,1999,20(3):30-31.
    [22]依呷,唐海,吕栋梁.低渗气藏启动压力梯度研究与分析[J].海洋石油,2006,26(3):51-54.
    [23]郝斐,程林松,李春兰,等.考虑启动压力梯度的低渗透油藏不稳定渗流模型[J].石油钻采工艺,2006,28(5):58-60.
    [24]郝斐,程林松,李春兰.特低渗透油藏启动压力梯度研究[J].西南石油学院学报,2006,28(6):29-33.
    [25]李爱芬,张磊,李春芹,等.特低渗透油藏气驱启动压力实验研究[J].断块油气田,2008,15(5):75-78.
    [26]尚根华.低渗透非线性渗流规律研究[D].中国科学院渗流流体力学研究所,2004.
    [27]黄延章.低渗透油层非线性渗流特征[J].特种油气藏,1997,4(1):9-14.
    [28]邓英尔,刘慈群.低渗油藏非线性渗流规律数学模型及其应用[J].石油学报,2001,22(4):72-76.
    [29]Young L C. A generalized Compositional Approach for Reservoir Simulation [J]. SPE10516,1983.
    [30]Coats K H. An Equation of State Compositional Model [J]. SPE8284,1980.
    [31]施文,桓冠仁,李福垲.一维三相全组份混相驱模型及其应用[J].石油学报,1995,16(4):75-83.
    [32]M.Noh, L.W.Lake, S.L.Bryant, and A Aragua-Martinez Implications of Coupling Fractional Flow and Geochemistry for CO2 Injection in Aquifers [J].SPE 2004:89341.
    [33]D.J.Wood, L.W.Lake, SPE, and R.T.Johns. A Screening Model for CO2 Flooding and Storage in Gulf Coast Reservoirs Based on Dimensionless Groups [J].SPE 2006:100021.
    [34]L.K.Folger and S.N.Guillot. A Case Study of the Development of the Sundown Slaughter Unit CO2 Flood [J]. SPE35189.
    [35]G.D.Fullbright, G.P.Hild, T.A.Korf, etc. Evolution of Conformance Improvement Efforts in a Major CO2 WAG Injection Project [J]. SPE35361.
    [36]Yih-BorChang, CoatsBK, NolenJS. A compositional model for CO2 flood including CO2 Solubility in water [J]. SPEREE,1998:155-160.
    [37]程林松,李春兰,陈月明.CO2在水中溶解的数值模拟方法[J].石油大学学报(自然科学版),1998,2(22):38-41.
    [38]张烈辉,张红梅,鲁友常,等.混相驱动态的拟四组分模型[J].西南石油学院学报,1999,21(4):64-66.
    [39]刘昌贵,孙雷,李士伦,等.多相渗流的几种数学模型及相互关系[J].西南石油学院学报.2002,24(1):64-66.
    [40]侯健.一种基于流线方法的CO2混相驱数学模型[J].应用数学和力学,2004,25(6):635-641.
    [41]计有权.C02混相驱多组分多相非等温数学模拟[D].吉林大学,2005.
    [42]盛金昌.C02驱油过程中多场全耦合数学模型[J].岩石力学与工程学报,2006,25(9):1893-1897.
    [43]程杰成,朱维耀,姜洪福.特低渗透油藏C02驱油多相渗流理论模型研究及应用[J].石油学报,2008,29(2):246-251.
    [44]沈平平,黄磊.二氧化碳-原油多相多组分渗流机理研究[J].石油学报,2009,30(2):247-251.
    [45]高海涛.低渗透油藏CO2驱油渗流规律宏观描述研究[D].中国石油大学,2009.
    [46]韩洪宝,程林松,张明禄.特低渗油藏考虑启动压力梯度的物理模拟及数值模拟方法[J].石油大学学报(自然科学版),2004,6(28):49-53.
    [47]赵国忠.变启动压力梯度三维三相渗流数值模拟方法[J].石油学报,2006,(27):119-123.
    [48]M. R. Islam, A. Chakma. Storage and utilization of CO2 in petroleum reservoirs-A simulation study [J].Energy Conversion and Management,1993,34(5):1205-1212.
    [49]John D R, Reid B G·A literature analysis of the WAG injectivity abnormalities in the CO2 process[C]·SPE 73830,2001:375-386.
    [50]Betty J. Felber, U.S. Department of Energy. Selected U.S. Department of Energy EOR Technology Applications[C]. SPE 89452,2004.
    [51]Gates G L, Caraway W H. Solvent Stimulation of Viscous Crude Oil Production [J]. SPE 3680,1972:3-4.
    [52]Patton J T, K H Coat, Ken Spence. Carbon Dioxide well Stimulation:Patel-A Parametric Study [J]. Journal of Petroleum Technology,1982:1-3.
    [53]M.A. Meadows, Fugro Seismic Imaging. Time-Lapse Seismic Modeling and Inversion of CO2 Saturation for Storage and Enhanced Oil Recovery[C]. Offshore Technology Conference,5-8 May 2008, Houston, Texas, USA.
    [54]Fixm N H, Beaty W G. Evaluation of CO2 flood performance, Spring "A" Sand, Ne Purdy Unit, Gavin Country, SPE 12665,1984.
    [55]曹学良,郭平,杨学峰,等.低渗透油藏注气提高采收率前景分析[J].天然气工业,26(3):100-102.
    [56]谈士海,张文正.非混相CO2驱油在油田增产中的应用[J].石油钻探技术,2001,29(2):58-60.
    [57]Monger T G, Coma J M. A laboratory and Field Evaluation of the CO2 Huff-Puff Process for Light Oil Recovery [J]. SPE 15501,1988:2-6.
    [58]Palmer F S, Landry R W. Design and Implementation of Immiscible Carbon Dioxide Displacement Projects (CO2 Huff-Puff) in South Louisiana [C], SPE15497,1986:5-8.
    [59]高正龙,杨铁梅,张俊法,等.文184断块二氧化碳混相驱及油藏数值模拟研究[J].江汉石油学院学报,2000,22(4):61-62.
    [60]雷友忠.低渗透油藏注CO2提高采收率技术与应用[D].西南石油大学,2006.
    [61]关振良,谢丛姣,齐冉,等.二氧化碳驱提高石油采收率数值模拟研究[J].天然气工业,2007,27(4):142-144.
    [62]史华.正理庄高89-1区块CO2混相驱参数优化数值模拟研究[J].海洋石油,2008,28(1):68-73.
    [63]宋道万.二氧化碳混相驱数值模拟结果的主要影响因素[J].油气地质与采收率,2008,15(4):72-74.
    [64]杨光璐.高246块二氧化碳非混相驱可行性分析及数值模拟[J].油气地质与采收率,2008,15(6):56-58.
    [65]张兴勇,王新海,江山,等.低渗透油藏中高含水期CO2驱数值模拟研究[J].长江大学学报(自然科学版),2010,7(2):189-191.
    [66]廖海婴.腰英台DB34井区CO2驱替油藏数值模拟研究[J].西安石油大学学报(自然科学版),2010,25(5):50-53.
    [67]付美龙,叶成,熊帆,等.茨34块CO2驱参数优化与方案设计[J].钻采工艺,2011,34(1):56-58.
    [68]高慧,张丽华,郎兆新.CO2吞吐开采稠油油藏的数值模拟研究[J].油气采收率技术,2000,7(4):13-15.
    [69]罗瑞兰,程林松.稠油油藏多轮次蒸汽吞吐后注CO2的可行性[J].新疆石油地质2004,25(2):182-184.
    [70]熊钰,曾明,武毅.兴北S3块CO2吞吐可行性及注采参数优化研究[J].西南石油学院学报,2005,27(2):34-36.
    [71]何应付,梅士盛,杨正明,等.苏丹Palogue油田稠油CO2吞吐开发影响因素数值模拟分析[J].特种油气藏,2006,13(1):64-67.
    [72]庞进,孙雷,孙良田,等.WC54井区C02单井吞吐强化采油室内实验及数值模拟研究[J].天然气勘探与开发,2006,29(2):44-48.
    [73]孙国强,孙雷,孙良田,等.CO2吞吐工艺操作参数的整体优化设计[J].钻采工艺,2006,29(4):47-50.
    [74]张国强,孙雷,孙良田,等.小断块单元油藏单井多周期CO2吞吐强化采油研究[J].大庆石油地质与开发,2006,25(5):85-89.
    [75]曾贤辉,熊玉滨,吴应川,等.文留油田CO2吞吐方案优化研究[J].内蒙古石油化 工2006,3:106-107.
    [76]欧瑾,刘学伟,何应付,等.苏丹Samaa稠油油藏CO2吞吐开采影响因素分析[J].大庆石油学院学报,2007,31(2):34-36.
    [77]陈民锋,姜汉桥,吴应川,等.复杂断块油藏单井CO2吞吐注采参数优化研究[J].油气田地面工程,2008,27(11):1-3.
    [78]F.M.orr. and M.K.Silva. Effect of oil composition on minimum Miscibility Pressure-Part2 correlation [C]. SPE14150,1987.
    [79][美]M.A.克林斯.二氧化碳驱油机理及工程设计[M].北京:石油工业出版社,1989:100-102.
    [80]李士伦,郭平,王仲林,等著.中低渗透油藏注气提高采收率理论及应用[M].北京市:石油工业出版社,2007:59-60.
    [81]秦同洛.关于低渗透油田的开发问题[J].断块油气田,1994,1(3):21-23.
    [82]Holm L W, Josendal V A. Mechanisms of oil displacement by carbon dioxide [J]. J Petrol Technology,1974,12:1427-1436.
    [83]Holm L W, Josendal V A. Effect of oil composition on miscible type displacement by carbon dioxide [J]. Soc Petron Engrs J,1982,2:87-98.
    [84]杨胜来.注C02气开采提高原油采收率的试验研究[D].中国石油大学博士后出站报告,1998.
    [85]Preditis, J., Wolle, C.R., Notz, P.K. A laboratory and field injectivity study CO2 WAG in the san Andres formation of west Texas[C]. SPE22653,1993.
    [86]谷丽冰,李治平,欧瑾.利用二氧化碳提高原油采收率研究进展[J].中国矿业,2007,16(10):66-67.
    [87]Christensen J R, Stenby E H. Skauge A. Review of WAG Field Experience [C]. SPE39883,1998.
    [88]路向伟,路佩丽.利用CO2非混相驱提高采收率的机理及应用现状[J].石油地质与工程,2007,21(2):58-61.
    [89]Grigg RB, Siagian UWR. Understanding and Exploiting Four-Phase Flow in Low Temperature CO2 Flood[C]. SPE39790,1998.
    [90]祝春生,程林松.低渗透油藏CO2驱提高原油采收率评价研究[J].钻采工艺,2007,30(6):55-61.
    [91]赵明国.特低渗透油藏CO2驱室内实验及数值模拟研究[D].大庆石油学院,2008.
    [92]刘一江,刘积松,黄忠桥,等.聚合物和二氧化碳驱油技术[M].北京:中国石化出版社,2001,98.
    [93]王进安,袁广均,张军,等.长岩心注二氧化碳驱油物理模拟实验研究[J].特种油气藏,2001,8(2):75-78.
    [94]郭平,李苗.低渗透砂岩油藏注CO2混相条件研究[J].石油与天然气地质,2007,28(5):689-690.
    [95]Christman P G, Gorell S B.A Comparison of Laboratory and Field-Observed CO2 Tertiary Injectivity [J]. JPT (Feb 1990).
    [96]Luo, R.Cheng L.S. and Peng, J.C. Feasibility Study of CO2 Injection for Heavy oil Reservoir after Cyclic Stream Stimulation:Liaohe Oil Field Test[C]. SPE97462.
    [97]郭万奎,廖广志.注气提高采收率技术[M].北京:石油工业出版社,2003:53-57.
    [98]Desch JB, Larsen WK. Enhanced Oil Recovery by CO2 Miscible displacement in the Little Knife Field, Billings County, North Dakota[C]. SPE10696,1984.
    [99]杨承志,岳清山,沈平平.混相驱提高石油采收率[M].北京:石油工业出版社,1989:52-63.
    [100]梅士盛.苏丹Palogue油田注CO2采油实验和数值模拟研究[D].中国科学院渗速流体力学研究所,2006.
    [101]莫增敏,石家雄.非混相CO2驱提高重油采收率[J].国外油田工程,2001,1:7-11.
    [102]林波.油藏条件下CO2的生成与CO2驱油的实验研究[D].西南石油学院,2001,40-41.
    [103]Claridge, E.L. CO2 flooding strategy in a communicating layered reservoir [J]. JPT DEC.1982.
    [104]Langston M V, Hoadlej S F, Young D N. Definitive CO2 Flooding Response in the SACROC Unit[C]. SPE71321.
    [105]Dai, K.K. and Orr, R.M. Jr. Prediction of CO2 Flood Performance:Interaction of Behavior with Microscopic Pore Structure Heterogeneity [J]. SPERE (Nov 1997).
    [106]张怀文,张翠林,多力坤.CO2吞吐采油工艺技术研究[J].新疆石油科技,2006,4(16):19-21.
    [107]Roper M KJ r, Pope G A. Analysis of tertiary injectivity of carbon dioxide [J]. SPE23974,1992.
    [108]刘一江,刘积松,黄忠桥,等.聚合物和二氧化碳驱油技术[M].北京:中国石化出版社,2001,98.
    [109]李淑霞,谷建伟.油藏数值模拟基础[M].北京:中国石油大学出版社,2009:44-46.
    [110]郑浩,马春华.基于正交试验法的低渗透油藏超前注水影响因素分析[J].2007,35(5):90-92.
    [111]茆诗松,周纪芗.概率论与数理统计(第二版)[M].北京:中国统计出版社,2000:232-245.
    [112]陈家鼎,孙山泽,李东风,刘力平.数理统计学讲义(第二版)[M].北京:高等教育出版社,2003:272-286.
    [113]贺文华,文竹青,娄涛.试验设计与数据处理[M].长沙市:国防科技大学出版社,2002:62-101.
    [114]董如何,肖必华,方永水.正交试验设计的理论分析方法及应用[J].安徽建筑工 业学院学报:自然科学版,2004,12(6):103-06.
    [115]韩兴刚,徐文,刘海锋.正交试验法在油气田开发方案优化设计中的应用[J].天然气工业,2005,25(4):116-18.
    [116]马成良,张海军,李素平.现代试验设计优化方法及应用[M].郑州市:郑州大学出版社,2007:121-68.
    [117]周德民,贺翠敏.社会统计学导论[M].长沙市:中南大学出版社,2004:3-11.
    [118]范勇.杏十三区聚合物驱试验数值模拟研究[D].大庆石油学院,2009:10-12.
    [119]井洋.欧力坨子油田开发效果评价及综合调整研究[D].大庆石油学院,2009,33-34.
    [120]战菲,宋考平,尚文涛,等.低渗透油藏单井C02吞吐参数优选研究[J].特种油气藏,2010,17(5):70-72.

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

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

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