基于SWAT的流域下垫面变化的水文响应研究
详细信息    本馆镜像全文|  推荐本文 |  |   获取CNKI官网全文
摘要
水资源与水环境问题是全球环境问题的研究热点。流域下垫面是影响水文循环的重要因素之一,而人类不合理的土地利用方式等活动使下垫面的结构发生改变,引起了水文过程的变化。研究下垫面变化的水文响应对水资源保护和生态规划具有重要的理论和现实意义。
     随着计算机、遥感和地理信息系统技术的发展,分布式水文模型已成为研究水文过程的重要手段。由美国农业部农业研究中心开发的SWAT (Soil and Water Assessment Tool)模型,以其强大的物理结构和多方面的模拟功能,在水文问题研究中发挥了重要作用。
     本文选择海河流域的漳河上游地区作为研究区域,应用SWAT模型对该区进行径流模拟,并在此基础上,通过设置各种下垫面变化情景,模拟下垫面变化的水文过程,研究水文要素对下垫面变化的响应程度和变化规律。主要研究成果如下:
     (1)搜集了流域的SWAT模型基础数据资料,通过图件资料处理、基础参数计算和文件格式转换,构建模型数据库。
     (2)比较和分析流域的子流域划分数量与径流量模拟结果关系,发现子流域划分数量对水系生成和径流模拟存在一定影响。结果显示,在漳河上游流域以300km2作为划分子流域的面积阂值,共划分37个子流域为子流域划分水平的最佳状态,能较真实反映流域的径流情况。
     (3)确定了CN2、SOL_AWC等8个敏感型参数,使率定过程合理有效。比较分析从上游到下游逐次校准和直接校准流域总出口两种校准方法的差异性,发现从上游到下游逐次校准更能保证上游子流域的模拟效率,具有更好的模拟结果。
     (4)模型在2001-2004年校准期R2和NS均大于0.71,结果较为理想,而在2006-2009年验证期评价系数不能完全达标,主要原因为由于人类活动导致2006-2009年的实测资料难以反映流域自然状态下真实的水文情况和水资源利用数据资料的短缺。
     (5)基于流域原土地利用类型,应用土地极端法和空间配置法,设置不同的下垫面变化情景,在2001-2004年时间段进行水文模拟,定量评估流域下垫面变化对流域陆地和河道两个阶段中多个水文要素(包括蒸散发、入渗、地表径流、产水量和河道径流量、蒸发损失、传输损失等)的影响。结果显示,城市用地的变化对各水文要素影响剧烈,城市化和耕地化在不同程度上减少入渗,增大地表径流、产水量、河道径流、河道传输损失。林地有很好的涵养水分作用,为减少水土流失,在流域内进行适当的退耕还林是很有必要的。
Water resources and water environmental problems are research hotspots of global environmental issues. The watershed underlying surface, which is one of important factors that influence the hydrological cycle, has been changed due to improper human activities, and it results in the changing of hydrological process. Therefore, studying on hydrological response to basin underlying surface change is significant for protecting water resources and formulating ecology plans.
     With the development of computer, remote sensing and GIS technologies, distributed hydrological model has become the important means of hydrological process studies. The physically-based SWAT(Soil and Water Assessment Tool) model, which is developed by USDA, has played an important role in hydrological research.
     In this paper, the upper Zhang river watershed of Hai river basin is selected. The SWAT model is applied to the streamflow simulation of this area, and then the hydrological response to basin underlying surface change in different scenes is simulated and analyzed. The main conclusions of the studv are as follows:
     (1) The basic data of upper Zhang river watershed are collected and the database of model is constructed by processing maps, calculating parameters and changing file format.
     (2) The relationship between the different partition numbers of subbasin and the results of streamflow simulation is analyzed and evaluated. The results indicate that the partition level of subbasin causes influence in river network generation and streamflow simulation. In order to ensure the natural conditions of hydrological process, the upper Zhang river watershed is delineated37subbasins with area threshold value of300km2
     (3) The values of8sensitive parameters are determined, which make the calibration process more reasonable and effective. Through comparison of the results of two calibration methods, the paper finds that calibrating from upstream to downstream has better simulation result than calibrating in total outlet directly.
     (4) The relative coefficient and Nash-Suttcliffe coefficient of streamflow in calibration period from2001to2004are satisfactory, while in verification period from2006to2009are unsatisfactory. The main reason is that the observational data from2006to2009are difficult to reflect the hydrological process in natural conditions and the data of water resources utilization is lacked.
     (5) Different underlying surface scenes are set up using extreme and spatial configuration methods. Through the simulation of hydrological process in each scene, the paper manages to analyze qualitatively the effects of underlying surface change on many hydrological elements, including of surface runoff, water yield, evapotranspiration and percolation in land phase and streamflow, evaporation losses, transmission losses in routing phase. The results indicate that the change of urban land causes dramatic influence in each hydrological element. The increased urban and farmland lead to reduced percolation, and increased surface runoff, water yield and streamflow. Forest has the capacity for water retention. In order to reduce water
引文
[1]高诞源,叶寿征,张君友,等.水文下垫面分析与分类初探[J].水文,1999,4:13-18.[2]Roland E,Schulze.Modeling Hydrological Responses to Land Use and Climate Change:A Southern African Perspective[J].Ambio,2000,29(l):13-17.[3]陈军锋,李秀彬.森林植被变化对流域水文影响的争论[J].自然资源学报,2001,16(5):476-480.[4]张志强,王礼先等.森林植被影响径流形成机制研究进展[J].自然资源学报,2001,16(1):79-84.[5]陈军锋,李秀彬.土地覆被变化的水文响应模拟研究.应用生态学报,2004,15(5):833-836.[6]马雪华.四川米亚罗地区高山冷杉林水文作用的研究[J].林业科学,1987,23(3):253-265.[7]Paulo R.L,Wolfram K.F,Nilson A.V.N,et al.Real evaporanspiration and transpiration throughat ropical rain forest in central Amazonia as estimated by the water balancemethod[J].Forestecology and management,1995,73(1-3):185-195.[8]康文星,等.杉木人工林水量平衡和蒸散的研究[J].植物生态学报与地植物学学报,1992,16(2):187-196.[9]Marcos H.C,Aurelie B,Jeffrey.A.C.Effects of large-scale changes in land cover on the discharge of the Tocantins River,Southeastern Amazonia[J] Journal of Hydrology,2003,28(3)206-217.[10]夏军,乔云峰,宋献方等.岔巴沟流域不同下垫面对降雨径流关系影响规律分析[J].资源科学,2007,29(1):70-76.[11]王根绪,张钰,刘桂民,等.马营河流域1967-2000年土地利用变化对河流径流的影响[J].中国科学D辑地球科学,2005,35(7):671-681.[12]杨君,杨敬坡.区域土地利用变化对水资源的影响研究[J].河北省科学院学报.2009,26(3):31-35.[13]徐振辞,郭永辰.城市不同下垫面条件的降雨径流模拟试验研究[J].南水北调与水利科技,2007,5(1):64-66.[14]Nassire E L,Jean R.A flood Damage Model for Flood Plain Studies [J].Water Resources Bulletin.1987,23(3):253-266.[15]陈莹,许有鹏,尹义星.土地利用/覆被变化下的暴雨径流过程模拟分析[J].地理科学,2009,29(1):117-123.[16]冉茂玉.论城市化的水文效应[J].四川师范大学学报(自然科学版),2000,23(4):436-439.[17]Paul M J,Meyer J L.Streams in the urban landscape [J]. Annual Review of Ecology and Systematics,2001,32:333-365.[18]杨舒媛,严登华,李杨,等.金沙江中上游下垫面变化的水文响应过程[J].人民长江,2008,39(15):39-41.[19]谢平,朱勇,陈广才,等.考虑土地利用/覆被变化的集总式流域水文模型及应用[J].山地学报,2007,25(3):257-264.[20]胡梦珺,刘文兆,赵姚阳.黄土高原农、林、草地水量平衡异同比较分析[J].干旱地区农业研??究.2003,21(4):113-116.[21]张北赢,徐学选,白晓华.黄土丘陵区不同土地利用方式下土壤水分分析[J].干旱地区农业研究2006,24(2):96-99.[22]王建,吴发启,孟秦倩等.不同利用类型土壤水分下渗特征试验研究[J].干旱地区农业研究,2006,24(6):159-162.[23]李丽娟,姜德娟,李九一等.土地利用/覆被变化的水文效应研究进展[J].自然资源学报,2007,22(2):211-224.[24]张蕾娜,李秀彬.用水文特征参数变化表征人类活动的水文效应初探—以云州水库流域为例[J].资源科学,2004,26(2):62-67.[25]Onstad C A,Jamieson D GModelling the effects of land use modifications on runoff[J].Water Resources Research,1970,6(5):1287-1295.[26]Arnold J QSrinivasan R,Muttiah R S,Williams,J R.Large area hydrologic modeling and assessment part I:Model developmen[J] Journal of the American Water Resources Association,1998,34(l):73-89.[27]http://swatmodel.tamu.edu/conferences.[28]Neitsch S L,Arnold J QKiniry J R,et al.Soil and water assessment tool theoretical documentation version2005[M].Grassland,Soil and Water Research Laboratory,Agricultural Research Service,2005.[29]Gassman P W,Reyes M R,Green C H,et al.The soil and water assessment tool:historical development,applications,and future research directions [J].Transactions of the AS ABE,2007,50(4):1211-1250.[30]Sophocleous M A,Koelliker J K,Govindaraju R S,et al,Ramireddygarie S R,Perkinsa S P.Integrated numerical modeling for basinwide water management:The case of the Rattlesnake Creek basin in south-central Kansans[J] Journal of Hydrology,1999,214:179-196.[31]Eckhardt K,Haverkamp S,Fohrer N,et al.SWAT-G,a version of SWAT99.2modified for application to low mountain range catchments [J].Physics and Chemistry of the Earth,2002,27(9-10):641-644.[32]Griensven A V,Bauwens W.Application and evaluation of ESWAT on the Dender basin and the Wister Lake basin[J].Hydrological Processes,2005,(19):827-838.[33]Krysanova V,Hatterman F,Wechsung F.Development of the ecohydrological model SWIM for regional impact studies and vulnerability assessment[J].Hydrological Processes,2005,19(3):763-783.[34]Du B,Arnold J G,Saleh A,et al.Development and application of SWAT to landscapes with tiles and potholes[J].Transactions of the ASABE,2005,48(3):1121-1133.[35]Du B,Saleh AJaynes D B,et al.Evaluation of SWAT in simulating nitrate nitrogen and atrazing fates in a watershed with tiles and potholes [J].Transactions of the ASABE,2006,49(4):949-959.[36]Pohlert T,Huisman J A,Breuer L,et al.Integration of a detailed biogeochemical model into SWAT for improved nitrogen predictions—Model development,sensitivity,and GLUE analysis [J].Ecological Model,2007,203:215-228.[37]Watson B M,McKeown R A,Gordon P,et al.Modification of SWAT for modelling streamflow from forested watersheds on the Canadian Boreal Plain [J] Journal of Environmental Engineering and Science,2008,7(15):145-159.[38]Allen P M,Arnold J QSkipwith W.Prediction of channel degradation rates in urbanizing watersheds [J].Hydrological Sciences Journal,2008,53(5):1013-1029.[39]Easton Z M,Fuka D R,Walter M T,et al.Re-conceptualizing the Soil and Water Assessment Tool(SWAT) model to predict runoff from variable source areas [J] Journal of Hydrology,2008,348:279-291.[40]Kim N W,Chung M,Won Y S,et al.Development and application of the integrated SWAT-MODFLOW model[J] Journal of Hydrology,2008,356:1-16.[41]Ouessar M,Bruggeman A,Abdelli F,et al.Modelling water-harvesting systems in the arid south of Tunisia using SWAT[J].Hydrology and Earth System Scinces,2009,13:2003-2021.[42]罗毅.分布式生态水文学模型研究取得重大进展:SWATMOD.2K4[J]中国西部环境和生态科学简报,2004,1(6):2-8.[43]张银辉,罗毅.基于分布式水文学模型的内蒙古河套灌区水循环特征研究[J].资源科学,2009,31(5):763-771.[44]李丹,薛联青,郝振纯.基于SWAT模型的流域面源污染模拟影响分析[J].环境污染与防治,2008,30(3):4-7.[45]张永勇,王中根,于磊,夏军,陈向东SWAT水质模块的扩展及其在海河流域典型区的应用[J].资源科学,2009,31(1):94-100.[46]李硕,康杰伟,王志华.基于输入文件定制的SWAT模型集成应用方法研究[J].地理与地理信息科学,2010,26(4):16-19.[47]张雪刚,毛媛媛,董家瑞,李致家SWAT模型与MODFLOW模型的耦合计算及应用[J].水资源保护,2010,26(3):49-52.[48]Milewski A,Sultan M,Yan E,et al.A remote sensing solution for estimating runoff and recharge in arid environments [J] Journal of Hydrology,2009,373:1-14.[49]于磊,邱殿明.基于SWAT模型的漳卫南流域水量模拟[J].吉林大学学报(地球科学版),2007,37(5):949-954.[50]Falkenmark M.Coping with Water Scarcity under Rapid Population Growth[R].Conference of SADC Minister,Pretoria,1995:23-24.[51]Falkenmark M,Rockstrom J.The new blue and green water paradigm:Breaking new ground for water resources planning and management[J] Journal of Water Resources Planning and Management,2006,132:129-132.[52]Ringersma JJBatjes N,Dent D.Green Water:Definitions and Data for Assessment (ISRIC Report2003/02)[R].Wageningen:World Soil Information Centre,2003.[53]Schuol J,Abbaspour K C,H Yang,et al.Modeling blue and green water availability in Africa[J].Water Resources Research,2008,44:1-18.[54]Faramarzi M, Abbaspour K C,Schulin R,et al.Modelling blue and green water resources availability in Iran[J].Hydrological Process,2008,23(3):486-501.[55]吴洪涛,武春友,郝芳华,金英学.绿水的多角度评估及其在碧流河上游地区的应用[J].资源科学,2009,31(3):420-428.[56]甄婷婷,徐宗学,程磊,王洁.蓝水绿水资源量估算方法及时空分布规律研究—以卢氏流域为例[J].??资源科学,2010,32(6):1177-1183.[57]杨菁荟,张万昌SWAT模型及其在水环境非点源污染中的应用研究进展[J].水土保持研究,2009,16(5):261-266.[58]Lam Q D,Schmalz B,Fohrer N.Modelling point and diffuse source pollution of nitrate in a rural lowland catchment using the SWAT model[J].Agricultural Water Management,2010,97:317-325.[59]Chiang L,Chaubey I.Differentiating impacts of land use changes from pasture management in a CEAP watershed using SWAT model[J].Transactions of the ASABE,2010,53(5):1569-1584.[60]唐莉华,林文婧等.基于SWAT的温榆河流域非点源污染模拟与分析[J].水力发电学报,2010,29(4):6-13.[61]Coffey R,Cummins E,Bhreathnach N,et al. Development of a pathogen transport model for Irish catchments using SWAT[J].Agricultural Water Management,2010,97(1):101-111.[62]Coffey R,Cummins E,Flaherty V O,et al.Analysis of the soil and water assessment tool(SWAT) to model Cryptosporidium in surface water sources[J].Biosystems Engineering,2010,106:303-314.[63]Ghaffari G,Keesstra S,Ghodousi J,et al.SWAT-simulated hydrological impact of land-use change in the Zanjanrood basin,Northwest Iran[J].Hydrological Processes,2010,24:892-903.[64]Jha M K,Schilling K E,Gassman P W,et al.Targeting land-use change for nitrate-nitrogen load reductions in an agricultural watershed[J] Journal of Soil snd Water conservation,2010,65(6):342-352.[65]秦耀民,胥彦玲,李怀恩.基于SWAT模型的黑河流域不同土地利用情景的非点源污染研究[J].环境科学学报,2009,29(2):440-448.[66]王艳君,吕宏军,施雅风,等.城市化流域的土地利用变化对水文过程的影响—以秦淮河流域为例[J].自然资源学报,2009,24(1):30-36.[67]Ficklin D L,Luo Y Z,Luedeling E,et al.Climate change sensitivity assessment of a highly agricultural watershed using SWAT[J] Journal of Hydrology,2009,374:16-29.[68]赵芳芳,徐宗学.黄河源区未来气候变化的水文响应[J].资源科学,2009,31(5):722-730.[69]Franczyk J,Chang H J.The effects of climate change and urbanization on the runoff of the Rock Creek basin in the Portland metropolitan area,Oregon,USA[J].Hydrological Processes,2009,23(6):805-815.[70]张圣微,雷玉平等土地覆被和气候变化对拉萨河流域径流量的影响[J].水资源保护,2010,26(2):39-44.[71]Chaubey I,Chiang L,Gitau M W,et al.Effectiveness of best management practices in improving water quality in a pasture-dominated watershed[J] Journal of Soil snd Water conservation,2010,65(6):424-437.[72]Tuppad P,Kannan N,Srinivasan R, et al.Simulation of agricultural management alternatives for watershed protection [J].Water Resour Manage,2010,24:3115-3144.[73]苏保林,王建平,贾海峰,等.密云水库流域非点源污染识别[J].清华大学学报(自然科学版),2006,46(3):360-365.[74]Yang Q,Meng F R,ZhaoZ Y, et al.Assessing the impacts of flow diversion terraces on stream water and sediment yields at a watershed level using SWAT model [J]. Agriculture, Ecosystems and Environment,2009,132:23-31.[75]Sahua M,Gu R R.Modeling the effects of riparian buffer zone and contour strips on stream water quality[J].Ecological Engineering,2009,35:1167-1177.[76]王良民,王彦辉.植被过滤带的研究和应用进展[J].应用生态学报,2008,19(9):2074-2080.[77]Immerzeel W W,Gaur A,et al.Integrating remote sensing and a process-based hydrological model to evaluate water use and productivity in a south Indian catchment[J]. Agricultural Water Management,2008,95:11-24.[78]Faramarzi M,Yang H,Schulin R,et al.Modeling wheat yield and crop water productivity in Iran:Implications of agricultural water management for wheat production[J].Agricultural Water Management,2010,97:1861-1875.[79]代俊峰,崔远来.基于SWAT的灌区分布式水文模型—Ⅰ.模型构建的原理与方法[J].水利学报,2009,40(2):145-152.[80]代俊峰,崔远来.基于SWAT的灌区分布式水文模型—Ⅱ.模型应用[J].水利学报,2009,40(3):311-318.[81]谢先红,崔远来.灌溉水利用效率随尺度变化规律分布式模拟[J].水科学进展,2010,21(5):681-689.[82]Srinivasan R,Zhang X,Arnold J.SWAT ungauged:Hydrological budget and crop yield predictions in the upper Mississippi River basin[J].Transactions of the ASABE,2010,53(5):1533-1546.[83]Baskaran L,Jager H I,Schweizer P E, et al.Progress toward evaluating the sustainability of switch-grass as a bioenergy crop using the SWAT model[J].Transactions of the AS ABE,2010,53(5):1547-1556.[84]Narasimhan B,Srinivasan R,Bednarz S, et al comprehensive modeling approach for reservoir water quality assessment and management due to point and nonpoint source pollution[J].Transactions of the AS ABE,2010,53(5):1605-1617.[85]Gassman P W,Osei E,Saleh A,et al.Alternative practices for sediment and nutrient loss control on livestock farms in northeast Iowa[J].Agriculture,Ecosystems&Environment,2006,117(2-3):135-144.[86]Cools J,Broekx S,Vandenberghe et al.Coupling a hydrological water quality model and an economic optimization model to set up a cost-effective emission reduction scenario for nitrogen[J].Environmental Modelling&Software,2010,26:44-51.[87]Sidney A,James H,James LV.Model Evaluation Process.USDA-EPA. http://www.epa.gov/scipoly/sap/1998/july/lpart5.pdf.[88]于涛,孟伟,Edwin Ongley,等.我国非点源负荷研究中的问题探讨[J].环境科学学报,2008,28(3):401-407.[89]Douglas-Mankin K R,Srinivasan R,Arnold J G.Soil and water assessment tool(SWAT)model:current development and applications [J].Transactions of the ASABE,2010,53(5):1423-1431.[90]Di Luzio M,Arnold J GFormulation of a hybrid calibration approach for a physically based dist ributed model with NEXRAD data input [J] Journal of Hydrology,2004,298:136-154.[91]Lee K S,Chung E S,Kim Y O.Integrated watershed management for mitigating streamflow depletion in an urbanized watershed in Korea[J].Physics and Chemistry of the Earth,2008,33:382-394.[92]Yang W H,Wang X X,Liu Y B,et al.Simulated environmental effects of wetland restoration scenarios??in a typical Canadian prairie watershed[J].Wetlands Ecol Manage,2010,18:269-279.[93]Lee M,Park G A,Park M J,et al.Evaluation of non-point source pollution reduction by applying Best Management Practices using a SWAT model and QuickBird high resolution satellite imagery[J] Journal of Environmental Sciences,2010,22(6):826-833.[94]徐爱兰,王鹏.基于SWAT模型的圩区农业非点源污染模拟[J].环境监控与预警,2010,2(1):38-43.[95]刘昌明,李道峰等.基于DEM的分布式水文模型在大尺度流域应用研究[J].地理科学进展,2003,22(5):437-445.[96]Sexton A,Sadeghi A,etal Using NEXRAD and rain gauge precipitation data for hydrologic calibration of SWAT in a northeastern watershed[J].Transactions of the ASABE,2010,53(5):1501-1510.[97]仕玉治,张弛,周惠成,等SWAT模型在水稻灌区的改进及应用研究[J].水电能源科学,2010,28(7):18-22.[98]苏保林,王建平,贾海峰,等.密云水库流域非点源模型系统[J].清华大学学报(自然科学版),2006,46(3):355-359.[99]桑学锋,周祖昊,秦大庸,等.改进的SWAT模型在强人类活动地区的应用[J].水利学报,2008,39(12):1377-1383.[100]宋艳华,马金辉SWAT模型在陇西黄土高原地区的适用性研究[J].干旱区地理,2007,30(6):933-938.[101]张思聪,刘铭环.竹竿河流域面源污染模拟计算和分析[J].水力发电学报,2006,25(5):51-56.[102]张利平,秦琳琳,胡志芳,等.南水北调中线工程水源区水文循环过程对气候变化的响应[J].水利学报,2010,41(11):1261-1271.[103]赖格英,于革,桂峰.太湖流域营养物质输移模拟评估的初步研究[J].中国科学D辑地球科学,2005,35(增刊Ⅱ):121-130.[104]荣琨,陈兴伟.晋江西溪流域土地利用变化对非点源污染影响的SWAT模拟[J].农业环境科学学报,2009,28(7):1488-1493.[105]郑璟,方伟华,史培军等.快速城市化地区土地利用变化对流域水文过程影响的模拟研究—以深圳市布吉河流域为例[J].自然资源学报,2009,24(9):1560-1572.[106]王中根,刘昌明,黄友波SWAT模型的原理、结构及应用研究[J].地理科学进展,2003,22(1):79-86.[107]唐晓晖,唐曙暇.海河流域漳河上游水环境的治理措施[J].江淮水利科技,2007,3:40-41.[108]牛巧荣,申志文,漳河上游水环境恢复初探[J].河北水利,2002.6,62-66.[109]郝芳华,程红光,杨胜天,非点源污染模型:理论方法与应用[M].中国环境科学出版社,2006.[110]郝芳华,张雪松,程红光等.分布式水文模型亚流域合理划分水平刍议[J].水土保持学报,2003,17(4):75-78[111]Manguerra H B, Engel B A.Hydrologic Parameterization of Watersheds for Runoff Prediction Using SWAT[J] Journal of the American Water Resources Association,1998,34(5):1149-1162.[112]Morris M D.Factorial sampling plans for preliminary computational experiments [J]. Technometrics,1991,33(2):161-174.[113]Duan Q Y,Gupta V K,Sorooshian S.Shuffled complex evolution approach for effective and??efficient global minimization[J] Journal of Optimization Theory and Application,1993,76(3):501-521.[114]Motovilov Y QGottschalk L,England K,et al.Validation of a distributed hydrological model against spatial observations [J]. Agricultural and Forest Meteorology,1999,(98-99):257-277.[115]Braemort K S,Arabi M,Frankenberger J R,et al. Modeling long-term water quality impact of structural BMPs[J].Transactions of the ASABE,2006,49(2):367-384.[116]Lenhart T, Eckhardt K, Fohrer N, et al. Comparison of two different approaches of sensitivity analysis [J].Physics and Chemistry of the Earth,2002,27:645-654.[117]朱新军,王中根SWAT模型在漳卫河流域应用研究[J].地理科学进展.2006,25(5):105-110.

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

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

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