摘要
水资源是人类生存与发展最宝贵的基础性自然资源,同时又是国民经济与社会发展的战略性自然资源。水资源短缺和生态环境恶化已经成为喀斯特流域所面临的严重问题。喀斯特流域可溶性的双重含水介质、独特的水文地貌结构及其产生的功能效应,使得其水系发育、水文动态上表现出与非喀斯特流域的巨大差异。水文过程是受地理、气候、环境、人类活动等等因素综合作用的结果,是复杂的动态过程,水文时间序列表现出高度非线性和多时间尺度特性,包含确定成分和随机成分。尺度问题,是水文水资源领域一个热点问题。研究喀斯特地区水文时间序列的多时间尺度问题,可以加强对水文尺度变化规律的认识,指导水资源的合理开发利用,协调人口、资源、环境和经济的持续发展。
本文采用统计学、小波理论、频谱分析等方法,对喀斯特流域水文时间序列的多时间尺度问题进行了研究。主要研究成果和结论如下:
(1)本研究以典型的喀斯特地区湘江流域为例,采用Morlet小波,对年径流量时间序列进行了多时间尺度分析,获得了周期性规律。年径流时间序列存在26a时间尺度的主周期变化,同时还存在7a时间尺度的次周期变化,在这两种时间尺度下,径流演变呈现出较明显的丰、枯交替振荡。
(2)利用滑动平均法和频谱分析方法对水文序列进行多时间尺度研究,各种分析方法的结果基本一致,也进一步说明小波理论所分析的湘江流域年径流量周期是可靠的。
(3)从气候、人类活动方面重点分析了喀斯特流域年径流序列演变的原因,并探讨其未来演变趋势。发现年径流量变化主要受年平均降水量的影响;人类活动的各种影响,总体上会使径流量略有增加。依据《气候变化国家评估报告》,未来贵州省喀斯特流域年径流量很可能呈增加趋势。
The water resources is the most precious foundational natural resources, is the strategic natural resource which the national economy and the society develops at the same time. Shortage of water resources and eco-environment degradation were the serious situation in karst drainage basins. There are doublely soluble aquifer media, unique hydro-geomorphological structure and the functional effects in karst drainage basins, whitch lead to a huge difference from non-karst drainage basins about the water system development and the hydrological dynamic. Hydrological process is a complex dynamic process, because it is affected by geography, climate, environment, human activities and so on. Hydrological time series demonstrated highly nonlinear and multiple time scales features, including identify and random compositions. Scale is a hot issue in the field of hydrology and water resources. Using the research of multiple time scales of hydrological time series in karst areas, we can enhanced the understanding of the laws on hydrological scale changes, guide the rational development and utilization of water resources, coordinate sustainable development including population, resources, environment and economy.
In this article, using statistics, wavelet theory, spectrum analysis method, the multiple time scales of hydrological time series in karst drainage basins are studied. Main research results and conclusions are as follows:
(1)In this study, taking the typical karst areas Xiangjiang River Basin as an example, using Morlet wavelet, the multiple time scales of the annual runoff time series analysed, the cyclical regularity is abtained. The annual flow time series has a mainly time scale cycle, which is 26 years. There is also a secondary time scale cycle of 7 years. At the two scales, runoff evolution shows obviously oscillation of rain and draught.
(2)The multiple time scales of hydrological series are studied through the sliding average method and spectral analysis, different analysis methods have basically the same results, whitch shows that runoff cycle of the Xiangjiang River basin analysed with the wavelet theory is reliable.
(3)The reasons of the evolution of annual runoff series are importantly analysed about the climate and human activities in karst drainage Basins, and the trend of their evolution in the future is discussed. The results show that the runoff changes are mainly affected by the average annual precipitation, the overall runoff will increase slightly on the various impacts of the human activities. Based on the "National Climate Change Assessment Report", the annual runoff will be likely to increase in karst drainage basins in Guizhou Province.
引文
[1] A·拉尔斯登,H·S·维尔夫.数字计算机上用的数学方法[M].上海:上海科学技术出版社:1983.
[2]蔡运龙.中国西南喀斯特山区的生态重建与农林牧业发展:研究现状与趋势[J].资源科学.1999(5), 37-40.
[3]常福宣,丁晶.嘉陵江流域洪水区域分析[J].长江流域资源与环境,2001,10(5):473-479.
[4]常福宣,丁晶.年最大洪峰区域变化的标度性[J].四川大学学报,(工程科学版)2001,33(1):5-8.
[5]常福宣,丁晶.降雨随历时变化标度性质的探讨[J].长江流域资源与环境,2002,11(1):79-83.
[6]陈传友,王春元,窦以松.水资源与可持续发展[M].北京:中国科学技术出版社,1999.86-110.
[7]陈喜.水文尺度和水文过程模拟研究[D].河海大学博士学位论文,2000.
[8]陈喜,陈永勤.日雨量随机解集研究[J].水利学报,2001(4):47-52.
[9]陈晓宏等.珠江三角洲河网区水位特征空间变异性研究[J].水利学报,2004,10:36-42.
[10]程琳,朱晓原,何惠等.国际水文合作与教育培训[J].水文,2006,26(3):81-84.
[11]丁晶,王文生,金菊良.论水文学中的尺度分析[J].四川大学学报(工程科学版),2003,35(3):9-13.
[12]丁晶,王文圣,金菊良.水文尺度分析[J].四川大学学报(工程科学版). 2003,35(2):1-4.
[13]丁晶.随机水文学[M].成都:成都科技大学出版社,1998.
[14]董长虹,高志,余啸海.Matlab 小波分析工具箱原理与应用[M].北京:国防工业出版社,2004.
[15]飞思科技产品研发中心.小波分析理论与 MATLAB7 实现[M].北京:电子工业出版社,2006.
[85]冯尚友.水资源持续利用与管理导论[M].北京:科学出版社,2000.133-139.
[16]傅国斌,李丽娟,刘昌明.遥感水文应用中的尺度问题[J].地球科学进展,2001, 16(6):755-760.
[17]高贵龙,邓自民,熊康宁,等.喀斯特的呼唤与希望[M].贵阳:贵州科技出版社,2003:1-165.
[18]郭彬锋,尹伟平.利用频谱分析法预测地下水动态[J].勘察科学技术,1993,(3):41-44.
[19]韩至钧,金占省.贵州省水文地质志[M].北京:地震出版社,1996.
[20]郝振纯,任立良,刘新仁.陆气耦合水文模型研究,淮河流域能量与水分循环研究(一)[M].北京:气象出版社,1999.
[21]胡昌华,张军波,夏军,等.基于 MATLAB 的系统分析与设计——小波分析[M].西安:西安电子科技大学出版社,2000.
[22]胡学龙,吴镇扬.数字信号处理教学指导[M]. 北京:高等教育出版社, 2007.
[23]黄友波,夏军,郑冬燕,等.应用频谱法分析水文时间序列的代表性[J].吉林水利,2002,20(10):1-4.
[24]黄忠恕.波谱分析方法及其在水文气象中的应用[M].北京:气象出版社,1983.
[25]蒋国荣,陈奕德,张韧,等.中国区域月降水量的多时间尺度特征[J].解放军理工大学学报(自然科学版),2002(4):77-81.
[26]孔兰,梁虹,姜维东.喀斯特流域枯水研究—以贵州省为例[J].水土保持研究,2007,14(5):288-292.
[27]孔兰,梁虹,贺向辉,等.喀斯特流域水资源问题及可持续利用对策[J].中国农村水利水电,2008,(3):17-19.
[28]孔兰,梁虹,姜维东.南明河流域水文序列预测方法研究[A].环境变化与水安全—第五届中国水论坛论文集[C].2008:778-782.
[29]孔兰,梁虹,黄法苏.喀斯特流域径流量时序演变特征分析——以贵州省为例.[J].中国岩溶.2007,26(4):341-346.
[30]孔兰,梁虹,黄法苏,等.喀斯特流域径流量多时间尺度的小波分析.[J].人民长江2008,39(5):17-18.
[31]冷疏影,宋长青,赵楚年,等.关于地理学科“十五”重点项目的思考[J].地球科学进展,2000,15(6):745-746.
[32]李长兴.论流域水文尺度化和相似性[J].水利学报,1995(1):277-283.
[33]李后强,程光钺.分形与分维[M].成都:四川教育出版社,1990.
[34]李继兵.Excel 数据处理与分析[M].北京:中国青年出版社,2006.
[35]李眉眉,丁晶,王文圣.基于混沌理论的径流降尺度分析[J].四川大学学报(工程科学版),2004,36(3):14-17.
[36]李平,卢文喜,王福林.频谱分析在挠力河流域年降水量中的应用[J].水文,2007,(27):25-27.
[37]李双成,蔡运龙.地理尺度转换若干问题的初步探讨[J].地理研究,2005,24(1):11-18.
[38]李贤彬,丁晶,李后强.水文时间序列的子波分析法[J].水科学进展,1999,10(2):144-149.
[39]李艳.水文系统变异的尺度研究[D].中山大学博士学位论文,2007.
[40]李玉柱,许炳南.贵州短期气候预测技术[M].北京:气象出版社,2001.149-154.
[41]梁虹.喀斯特流域尺度与枯水流量初步研究[J].贵州师范大学学报,1997,15(3):1-5.
[42]梁虹.喀斯特流域空间尺度对洪、枯水水文特征值影响初探?以贵州省河流为例[J].中国岩溶,1997,16(2):121-129.
[42]梁虹,王剑.喀斯特地区流域岩性差异与洪、枯水水文特征值相关分析?以贵州省河流为例[J].中国岩溶,1998,17(1):67-73.
[44]林振山,邓子旺.子波气候诊断技术的研究[M].北京:气象出版社,1991,22-25.
[45]刘建梅,裴铁璠.水文尺度转换研究进展[J].应用生态学报,2003,14(12):2305-2310.
[46]刘建梅,王安志,裴铁璠,等.杂谷脑河径流趋势及周期变化特征的小波分析[J].北京林业大学学报,2005,27(4):49-55.
[47]刘俊萍,田峰巍,黄强,等.基于小波分析的黄河河川径流变化规律研究[J].自然科学进展,2003,13(4):383-387.
[48]刘贤赵.论水文尺度问题[J].西北林学院学报,2004,19(3):84-88.
[49]刘新仁.大尺度水文模拟若干问题的探讨[A].第三届全球能量与水循环国际科学大会“淮河流域能量与水文循环研究”[C].北京:气象出版社,1999.
[50]龙兴明,周静.连续小波变换的一维信号检测[J].重庆邮电学院学报,2004,16(3):77-80.
[51]卢岱.SPSS for Windows 统计方法[M].北京:电子工业出版社.2006.
[56]卢晓宁,邓伟,张树清,等.霍林河中游径流量序列的多时间尺度特征及其效应分析[J].自然资源学报,2006,21(5):819-826.
[57]孟斌,王劲峰.地理数据尺度转换方法研究进展[J].地理学报,2005,2(60):277-288.
[58]彭贤伟,熊康宁.喀斯特地区可持续发展指标体系构建与可持续发展能力评价[C].中国西南(贵州)喀斯特生态环境治理与持续发展咨询会议论文集,2003:198-199.
[59]齐善忠,王涛,罗芳等.黑河流域水资源可持续利用对策[J].中国人口.资源与环境, 2004. 14(6):58-60.
[60]苏凤阁,郝振纯.水文模型中雨量资料的解集分析及应用[J].气候与环境研究,2001,6(2):261-266.
[61]苏金明,张莲花,刘波,等编.MATLAB 工具箱应用[M].北京:电子工业出版社,2004.
[62]谭浩强.FORTRAN 程序设计教程[M].北京:清华大学出版社,2000.
[63]王红瑞,叶乐天,刘昌明,等.水文序列小波周期分析中存在的问题及改进方式[J].自然科学进展,2006,16(8):1002-1008.
[64]王文圣,丁晶,向红莲.小波分析在水文学中的应用研究及展望[J].水科学进展,2002,13(4):515-517.
[65]王文圣,丁晶,向红莲.水文时间序列多时间尺度分析的小波变换法.四川大学学报(工程科学版),2002,34(6):14-17.
[66]王文圣,丁晶,李跃清.水文小波分析[M].北京:化学工业出版社,2005.
[67]王在高,梁虹,杨明德.喀斯特流域地貌类型对枯水径流的影响?以贵州省为例[J].地理研究,2002,21(4):441-448.
[68]魏林宏.时空尺度对洪水模拟的影响研究[D].河海大学博士学位论文,2004.
[69]乌建国.景观生态学一格局、过程、尺度与等级[M].北京:高等教育出版社,2000.
[70]吴险峰,刘昌明.流域水文模型研究的若干进展[J].地理科学进展,2002,21(4):341-348.
[71]夏军.《水文非线性系统理论与方法》[M].武昌:武汉大学出版社,2002.
[72]夏军.水文尺度问题[J].水利学报,1993,(5):32-37.
[73]邢纪元.谐波分析在短期气侯预测中的应用[J].山东气象,2005,25(2):14-15.
[74]熊立华,郭生练.分布式流域水文模型[M].北京:中国水利水电出版社,2004.
[75]徐建华.现代地理学中的数学方法[M].北京:高等教育出版社.2006.
[76]徐娜.警惕气候变暖——气候变化国家评估报告简介[J].中国减灾,2007(2):18-19.
[77]徐英,陈亚新,史海滨,等. 土壤水盐空间变异尺度效应的研究[J].农业工程学报,2004,20(2):1-5.
[78]严登华,王浩,王建华,等.国际水文计划发展与中国水资源研究体系构建[J].地理学报,2004,59(2):249-259.
[79]杨福生.小波变换的工程分析与应用[M].北京:科学出版社,2001,1-6.
[80]杨明德,谭明,梁虹.喀斯特流域水文地貌系统[M].北京:地质出版社,1998.1-19.
[81]杨明德.喀斯特研究——杨明德论文选集[M].贵阳:贵州民族出版社,2003.
[82]叶守泽,夏军.水文科学研究的世纪回眸与展望[J].水科学进展,2002,13(1):93-104.
[83]叶泽刚,汪娜娟,周文波.降水特性空间变异性初步研究[J].水利学报.1994,(12):7-13.
[84]詹道江,叶守泽.工程水文学[M].北京:中国水利水电出版社,2000.
[85]张超,杨秉赓.计量地理学基础[M].北京:高等教育出版社,1991.
[86]张少文,丁晶,廖杰,等.基于小波的黄河上游天然年径流变化特性分析[J].四川大学学报(工程科学版),2004,36(3):32-37.
[87]张少文.黄河流域天然年径流变化特性分析及其预测[D].四川大学博士学位论文,2005.
[88]张学真,刘燕.灞河出山径流序列变化的小波分析[J].水资源保护.2006,22(3): 12-15.
[89]赵文武,傅伯杰,陈利顶.尺度推绎研究中的几点基本问题[J].地球科学进展. 2002,17(6):905-911.
[90]钟晔,金昌杰,裴铁璠.水文尺度转换探讨[J].应用生态学报,2005,16(8):1537-1540.
[91]周忠发,安裕伦.贵州省水土流失遥感现状调查及空间变化分析[J].水土保持通报,2000,20(6):23-25.
[92]左其亭,王中根.《现代水文学》[M].郑州:黄河水利出版社,2006,
[93]左其亭,王中根.中国西部流域水循环重大科学问题及研究展望[J].西北水资源与水工程,2001,12(4):1-5.
[1] http://www.chinaihp.org/main.htm.
[2]http://www.icsu.org/index.php.
[3]http://www.igbp.kva.se/.
[4]http://www.unesco.org/water/ihp/.
[5]A.Martinez-Cob. Multivariate geostatistical analysis of evapotranspiration and precipitation in mountainous terrain[J]. Journal of Hydrology. 1996,174:19-35.
[6]B Sivakumar.A chaotic approach to rainfall disaggregation[J]. Water Resources Research,2001,37(1):61-72.
[7]Becker A,Braun P.Disaggregation,aggregation and spatial scaling in hydrological modelling[J].Joumal of Hydrology,1999,217(3-4):239—252.
[8]Bloschl G,Sivapalan M (Eds.). Special issue on scale issues in hydrological modeling [J].Hydrology Processes,1995,9:251-290.
[9]Bloschl G.,Sivapalan M.. Scale issues in hydrological modeling:a review[J].Hydrological Processes,1995,9:251-290.
[10]Cao,C.Y.,and Lam, N.S-N. Understanding the scale and resolution effects in remote sensing and GIS,Scaling in remote sensing[M].ed by Quanttrochi,D. A.,and Goodchild,M. F.,CRC Press,1997,57-73.
[11]CiarapicaL,Todini E. TOPKAPI:A model for the representation of the rainfall-runoff process at different scales[J].Hydrol Process,2002,(16):207-229.
[12]Connolly RD,Silburn DM. Distributed parameter hydrology model(ANSWERS) applied to a range of catchment scales using rainfall simulator dataⅡ.Application to spatially uniform catchments[J].J Hydrol,1995(172):105-125.
[13]Cushman JH. Unifying the concepts of scale,instrumentation,and stochastics in the development of multiphase transport theory[J].Water Resour Res. 1984,20:1668-1676.
[14]Dsvid Labat.Recent advances in wavelet analyses:Part 1.A review of concepts[J].Journal of Hydrology, 2005,314,275-288.
[15]Entekhabi D,Eagleson P.. Land surface hydrology parameterization for atmospheric general circulation models including subgridscale spatial variability[J].J Clim,1989,2(8):816-831.
[16]Gupta V K,Waymine E. Multiscaling properties of special rainfall and river flow distribution[J].J Geophysical Res.,1990,95(3):1999-2009.
[17]Gupta VK,Mesa O J,Dawdy D R. Multiscaling theory of flood peaks:Regional quantile analysis[J]. Water Resource Res.,1994,30(12):3405-3421.
[18]Hall, David;Fitzpatrick, Caroline S.B.Spectral analysis of pressure variations duringcombined air and water backwash of rapid gravity filters[J].Water research ,1999,33(17):3666-3672.
[19]James C.I. Dooge,Scale Problems in Hydrology,in:Reflections on Hydrology,American Geophysical Union,1997.
[20]Kirkby MJ,Baird AJ,Diamond SM,etal.A physically basically based process model for hydrology,ecology and land degradation interactions[A]. In:Brandt CJ,Thornes JB,eds. Mediterranean Desertification and LandUse[C].Chichester:Wiley,1996.303-354.
[21]Klemes,V..Conceptualisation and scale in hydrology[J].Journal of Hydrology,1983,65,1-23.
[22]Kumar P,Foufoula-Georgiou E. A multicomponent decomposition of spatial rainfall fields 1. segregation of Large-and Small-Scale features using Wavelet transforms[J].Water Resources Research.1993,29(8):2515-2532.
[23]Kumar P, Foufoula-Georgiou E. A multicomponent decomposition of spatial rainfall fields 1. segregation of Large-and Small-Scale features using wavelet transforms[J].Water Resources Research.1993,29(8):2515-2532.
[24]Lam,N.,Quattrochi,D. A..On the issues of scale,resolution,and fractal analysis in the mapping sciences[J]. Prof. Geogr.,1992,(44):88-98.
[25]Margriet Nakken. Wavelet analysis of rainfall-runoff variability isolating climatic from anthropogenic patterns[J]. Environmental Modelling & Software,1999,14:283-295.
[26]Rigon R,Rodriguez-Iturbe I,Maritan A,et al. On Hack.s law[J].WatResourRes,1996.11 (32):3367-3374.
[27]Roland Schulze. Transcending scales of space and time in impact studies of climate and climate change on agrohydrological responses[J]. Agriculture,Ecosystems and Environment,2001,82:185-212.
[28]Sposito G. Scale dependence and scale invariability in hydrology[M]. Cambridge: Cambridge University Press,1998.6-31.
[29]Sten Bergstrom,Phil Graham L. On the scale problem in hydrological modeling[J]. Journal of Hydrology,1998,211:253-265.
[30]Turgay Partal, Murat K. Long-term trend analysis using discrete wavelet components of annual precipitations measurements in Marmara region (Turkey)[J]. Physics and Chemistryof the Earth. 2006,31: 1189–1200.
[31]Venckp V , Foufoula-GeorgiouE. Energy decomposition of rainfall in the time-frequency-scale domain using wavelet packets[J]. Journal of Hydrology,1996,27(3):263-271.