非一致性年径流序列频率分析计算
|
摘要
由于受到气候变化和人类活动的频繁影响,流域(或区域)内的下垫面条件发生了较大变化,破坏了天然水文序列的一致性,使其不能直接采用传统的水文频率计算方法进行计算(谢平等2005a)。本文以陕北及关中地区的非一致性年径流资料为例,采用三种不同的方法进行了频率计算研究,以期为该地区水利工程的规划、设计、施工和运行管理提供水文依据。主要研究内容和结论如下:
(1)水文变异综合诊断。采用水文变异综合诊断方法对研究区15个水文测站的年径流序列进行诊断。初步诊断包括过程线法、累积曲线斜率差异幅度分析法和Hurst系数法。详细诊断包括趋势详细诊断和跳跃详细诊断:对于趋势变异,采用线性趋势相关系数检验法、斯皮尔曼(Spearman)秩次相关检验法和肯德尔(Kendall)秩次相关检验法进行判断;对于跳跃变异,采用Lee-Heghinian法、有序聚类法、最优信息二分割法、R/S法、Brown-Forsythe法、滑动F法,滑动T法、滑动秩和法、滑动游程法和Mann-Kendall法进行判断。综合诊断包括趋势综合和跳跃综合,并根据效率系数确定序列变异形式。结合实际水文调查分析,对变异形式和结论进行确认,进而得到最可能的变异诊断结果。经诊断,研究区15个水文测站中,有7个测站的的年径流序列发生变异,其中属于趋势变异的测站为赵石窑站;属于跳跃变异的测站和变异年份分别为:神木站1979年、枣园站1996年、黄陵站1985年、林家村站1985年、张家山站1985年和状头站1985年。
(2)基于随机方法的非一致性年径流频率计算。根据水文变异综合诊断的结果,对赵石窑站采用基于趋势分析的非一致性年径流序列频率计算方法进行计算;对神木站、枣园站、黄陵站、林家村站、张家山站和状头站采用基于跳跃分析的非一致性年径流序列频率计算方法进行计算,得到了过去、现在和未来条件下满足P-Ⅲ型分布的统计参数。通过对7个变异测站的统计参数以及频率曲线的分析可知, 7个测站的年径流量在未来条件下存在不同程度的减少,也就是说,其年径流序列的总体变化趋势是逐年减小的。
(3)基于条件概率的非一致性年径流序列频率计算。对基于条件概率的非一致性洪水频率计算方法进行转换和改进,得到研究区基于条件概率的非一致性年径流频率计算方法:首先,根据水文变异综合诊断结果,对于非一致性年径流序列,按其变异点的位置将序列进行划分。其次,对概率权重和条件概率进行估计。本文对条件概率的估计假设了两种情况,即条件概率服从对数正态分布或服从P-Ⅲ分布。最后,采用ABS、OLS和WLS准则对频率计算结果进行拟合优度检验。从研究区7个变异测站的频率曲线拟合优度情况来看,基于条件概率的非一致性年径流频率计算方法能满足参数估计精度的要求。
(4)基于混合分布的非一致性年径流序列频率计算。对基于混合分布的非一致性洪水频率计算方法进行转换和改进,得到研究区基于混合分布的非一致性年径流频率计算方法,并将其应用到研究区6个跳跃变异测站的年径流资料上。首先,假设非一致性年径流水文序列的变异点前后分别服从于两个分布,全序列服从于由这两个分布组成的混合分布,本文采用的两对数正态混合分布和两P-Ⅲ混合分布。其次,采用模拟退火算法进行参数估计。最后,采用ABS、OLS和WLS准则对频率计算结果进行拟合优度检验。
基于混合分布的非一致性年径流频率分析计算过程,包括统计分析和成因分析,使其较之传统的频率计算方法更具合理性;其次,混合分布的子分布可以根据拟合优度的要求来选取,这奠定了其方法上的灵活性;再次,采用模拟退火算法估计混合分布参数,可以满足各种分布参数求解的要求,不受参数个数的限制;最后,从拟合优度来看,混合分布模型优于条件概率模型。
Due to the frequent human activities and the climate change, partial climate and underlying surface conditions were changed considerablly. All this change destroyed the consistence of natural hydrological data, which made it cannot use traditional frequency analysis methods. This research used three methods to calculate the frequency of inconsistent annual runoff serise in Shanbei and Guanzhong region, in hopes to provide hydrological design values for these regions. The main contents and the conclusions are as follows:
(1) Comprehensive hydrology change point analysis. Adopt the comprehensive hydrology change point analysis method to determine the change-point of inconsistent annual runoff in the regions. Preliminary diagnosis includes the process line method, the differences of cumulative slope method and Hurst coefficient method. Detailed diagnosis includes trend detailed diagnosis and jump point for diagnosis of detailed. Three methods for trend detailed diagnosis were enployed: correlation coefficient method, Spearman correlation test and Kendall rank correlation test. Using 11 methods for jump detailed diagnosis: Lee-Heghinian method, orderly clustering method, two segmentation method of optimal information, R/S method, Brown-Forsythe method, sliding F method, sliding T method, sliding rank method, sliding run method and Mann-Kendall method. Comprehensive diagnosis includes trend integrated diagnosis and jump integrated diagnosis, and the change-point can be determined according to efficiency coefficients and practical survey analysis. Through analysis, there are 7 stations’annual runoff series hydrological changed: Zhaoshiyao belong to trend type, and the rest of the annual runoff series belong to jump type, and the change point is Shenmu 1979, Zaoyuan 1996, Huangling 1985, Linjiacun 1985, Zhangjiashan 1985 and Zhuangtou 1985.
(2) For this 7 variation stations, using the inconsistent annual runoff frequency calculation method which based on random analysis, we can gain the statistical parameters of P-Ⅲdistribution in the past, present and future. And the results indicated that the annual runoff series reduced year by year.
(3) Frequency calculation of inconsistent annual runoff series based on conditional probability. Firstly, according to change-point, classified inconsistent annual runoff series as two parts. Secondly, eatimate the probability weight and conditional probability(in this paper probability include LN(2)and P-Ⅲ). Finally, by use ABS, OLS and WLS standard test the goodness-of-fit, we can find the frequency curve fitting well.
(4) Frequency calculation of inconsistent annual runoff series based on mixed distribution. Firstly, according to change-point, clasiified inconsistent annual runoff series as two parts. Assumpt this two series follow differently distribute, and the whole series follows mixed distribution(mixed-two lognormal distribution and mixed-two P-Ⅲdistribution). Secondly, use the simulated annealing algorithm estimate the parameters of mixed distribution. Finally, ABS, OLS and WLS standard were tested the goodness-of-fit, we can find the frequency curve fitting well.
The mixed distribution contains statistic analysisy and casue-effect analysis, is more rational and flexibal than traditonal methods. Acorrding to calculation results, the mixed distributiom model is better than to conditional probability model.
(1)水文变异综合诊断。采用水文变异综合诊断方法对研究区15个水文测站的年径流序列进行诊断。初步诊断包括过程线法、累积曲线斜率差异幅度分析法和Hurst系数法。详细诊断包括趋势详细诊断和跳跃详细诊断:对于趋势变异,采用线性趋势相关系数检验法、斯皮尔曼(Spearman)秩次相关检验法和肯德尔(Kendall)秩次相关检验法进行判断;对于跳跃变异,采用Lee-Heghinian法、有序聚类法、最优信息二分割法、R/S法、Brown-Forsythe法、滑动F法,滑动T法、滑动秩和法、滑动游程法和Mann-Kendall法进行判断。综合诊断包括趋势综合和跳跃综合,并根据效率系数确定序列变异形式。结合实际水文调查分析,对变异形式和结论进行确认,进而得到最可能的变异诊断结果。经诊断,研究区15个水文测站中,有7个测站的的年径流序列发生变异,其中属于趋势变异的测站为赵石窑站;属于跳跃变异的测站和变异年份分别为:神木站1979年、枣园站1996年、黄陵站1985年、林家村站1985年、张家山站1985年和状头站1985年。
(2)基于随机方法的非一致性年径流频率计算。根据水文变异综合诊断的结果,对赵石窑站采用基于趋势分析的非一致性年径流序列频率计算方法进行计算;对神木站、枣园站、黄陵站、林家村站、张家山站和状头站采用基于跳跃分析的非一致性年径流序列频率计算方法进行计算,得到了过去、现在和未来条件下满足P-Ⅲ型分布的统计参数。通过对7个变异测站的统计参数以及频率曲线的分析可知, 7个测站的年径流量在未来条件下存在不同程度的减少,也就是说,其年径流序列的总体变化趋势是逐年减小的。
(3)基于条件概率的非一致性年径流序列频率计算。对基于条件概率的非一致性洪水频率计算方法进行转换和改进,得到研究区基于条件概率的非一致性年径流频率计算方法:首先,根据水文变异综合诊断结果,对于非一致性年径流序列,按其变异点的位置将序列进行划分。其次,对概率权重和条件概率进行估计。本文对条件概率的估计假设了两种情况,即条件概率服从对数正态分布或服从P-Ⅲ分布。最后,采用ABS、OLS和WLS准则对频率计算结果进行拟合优度检验。从研究区7个变异测站的频率曲线拟合优度情况来看,基于条件概率的非一致性年径流频率计算方法能满足参数估计精度的要求。
(4)基于混合分布的非一致性年径流序列频率计算。对基于混合分布的非一致性洪水频率计算方法进行转换和改进,得到研究区基于混合分布的非一致性年径流频率计算方法,并将其应用到研究区6个跳跃变异测站的年径流资料上。首先,假设非一致性年径流水文序列的变异点前后分别服从于两个分布,全序列服从于由这两个分布组成的混合分布,本文采用的两对数正态混合分布和两P-Ⅲ混合分布。其次,采用模拟退火算法进行参数估计。最后,采用ABS、OLS和WLS准则对频率计算结果进行拟合优度检验。
基于混合分布的非一致性年径流频率分析计算过程,包括统计分析和成因分析,使其较之传统的频率计算方法更具合理性;其次,混合分布的子分布可以根据拟合优度的要求来选取,这奠定了其方法上的灵活性;再次,采用模拟退火算法估计混合分布参数,可以满足各种分布参数求解的要求,不受参数个数的限制;最后,从拟合优度来看,混合分布模型优于条件概率模型。
Due to the frequent human activities and the climate change, partial climate and underlying surface conditions were changed considerablly. All this change destroyed the consistence of natural hydrological data, which made it cannot use traditional frequency analysis methods. This research used three methods to calculate the frequency of inconsistent annual runoff serise in Shanbei and Guanzhong region, in hopes to provide hydrological design values for these regions. The main contents and the conclusions are as follows:
(1) Comprehensive hydrology change point analysis. Adopt the comprehensive hydrology change point analysis method to determine the change-point of inconsistent annual runoff in the regions. Preliminary diagnosis includes the process line method, the differences of cumulative slope method and Hurst coefficient method. Detailed diagnosis includes trend detailed diagnosis and jump point for diagnosis of detailed. Three methods for trend detailed diagnosis were enployed: correlation coefficient method, Spearman correlation test and Kendall rank correlation test. Using 11 methods for jump detailed diagnosis: Lee-Heghinian method, orderly clustering method, two segmentation method of optimal information, R/S method, Brown-Forsythe method, sliding F method, sliding T method, sliding rank method, sliding run method and Mann-Kendall method. Comprehensive diagnosis includes trend integrated diagnosis and jump integrated diagnosis, and the change-point can be determined according to efficiency coefficients and practical survey analysis. Through analysis, there are 7 stations’annual runoff series hydrological changed: Zhaoshiyao belong to trend type, and the rest of the annual runoff series belong to jump type, and the change point is Shenmu 1979, Zaoyuan 1996, Huangling 1985, Linjiacun 1985, Zhangjiashan 1985 and Zhuangtou 1985.
(2) For this 7 variation stations, using the inconsistent annual runoff frequency calculation method which based on random analysis, we can gain the statistical parameters of P-Ⅲdistribution in the past, present and future. And the results indicated that the annual runoff series reduced year by year.
(3) Frequency calculation of inconsistent annual runoff series based on conditional probability. Firstly, according to change-point, classified inconsistent annual runoff series as two parts. Secondly, eatimate the probability weight and conditional probability(in this paper probability include LN(2)and P-Ⅲ). Finally, by use ABS, OLS and WLS standard test the goodness-of-fit, we can find the frequency curve fitting well.
(4) Frequency calculation of inconsistent annual runoff series based on mixed distribution. Firstly, according to change-point, clasiified inconsistent annual runoff series as two parts. Assumpt this two series follow differently distribute, and the whole series follows mixed distribution(mixed-two lognormal distribution and mixed-two P-Ⅲdistribution). Secondly, use the simulated annealing algorithm estimate the parameters of mixed distribution. Finally, ABS, OLS and WLS standard were tested the goodness-of-fit, we can find the frequency curve fitting well.
The mixed distribution contains statistic analysisy and casue-effect analysis, is more rational and flexibal than traditonal methods. Acorrding to calculation results, the mixed distributiom model is better than to conditional probability model.
引文
陈广才,谢平. 2006.水文变异的滑动F识别与检验方法.水文, 26(2):57～60
陈望春,赵立锋.2002.权函数法在P-Ⅲ型分布中的应用.浙江水利科技,4:46~47
陈晓林,汪四水.2007.一类混合正态分布参数估计的EM算法和数据扩张.苏州大学学报(自然科学版), 23(3):13~18
崔维成,徐向东,邱强.1998.一种快速计算随机变量函数均值与标准差的新方法.船舶力学,(6):50~60
丁晶,邓育仁.1988.随机水文学.成都:成都科技大学出版社:64~110
董亚非,李兆元. 1995.关中地区气候资源特征分析.陕西气象,6:34~36
关中平原水资源考察小组.1990.关中平原水资源考察报告.杨凌:西北农业大学
韩小雷.2008.粒子群—模拟退火融合算法及其在函数优化中的应用.武汉:武汉理工大学出版社
胡四一.2006.水文学及水资源的学科前沿和研究要求.水文, 26(3): 1~4
金光炎.1993.水文水资源随机分析.北京:中国水利水电出版社
金光炎.1999.水文频率分析述评.水科学进展,10(3): 319~327
金光炎.2004.洪水频率分布模型的类别、比较和选用.安徽水利科技, (2): 36~38
雷红富,谢平,陈广才,李晶.2007.水文序列变异点检验方法的性能比较分析.水电能源科学, 25(4):36~40
李宏伟,宋松柏.2009.蚁群算法在水文频率曲线参数计算中的应用.人民黄河, 31(4):38-40
李立民,李玉波.2007.西安至宝鸡铁路客运专线可行性研究报告.西安:铁道第一勘察设计院
李述山,闫保英.1999.两正态混合分布的一种参数统计方法.山东轻工业学院学报, 13(2):73~77
李艳,陈晓宏,张鹏飞.2007.基于相关积分法的水文时间序列变异性度量.水文, 27(3):6~9
李志,曹明明.2003.关中地区水资源问题及对策研究.水土保持学报, 17(3):159~161
梁国文,佟伟.2008.水文频率曲线与水文随机过程.黑龙江水利科技, 5(36):162
梁艳琴.1994.用现状还原年径流模拟径流系列的探讨.海河水利, (3) :39～41
刘昌民,成立.2000.黄河干流下游断流的径流序列分析.地理学报, 55(3):257～265
刘光祖.2000.概率论与数理统计.北京:高等教育出版社
刘俊民,孙勇.2006.关中地区水资源特性分析及可持续利用对策.人民黄河, 28(1):39~41
陆中央.2000.关于年径流量系列的还原计算问题.水文, 20(6):9~12
马志明,刘瑞元,习丽.2007.多个子总体混合分布的参数估计.西北民族大学学报(自然科学版),28(65):11~15
曲小红,石生新,荣丰涛.1999.降雨径流系列的一致性分析.山西水利科技, (1):34~37
任树梅,李靖.2005.工程水文与水利计算.北京:中国农业出版社
桑燕芳,王栋.2008.水文时间序列周期识别的新思路与两种新方法.水科学进展,19(3):412~417
陕西省地方志编撰委员会.1999.陕西省志.第十三卷.水利志.西安:陕西人民出版社:65~69
沈宏.2003.天然径流还原计算方法初步探讨.水利规划与设计, (3):15~18
史鉴,陈兆丰,邢大韦.2002.关中地区水资源合理开发利用与生态环境保护.郑州:黄河水利出版社
史卫东.2001.水文资料系列的一致性分析.甘肃水利水电技术, 37(1):22~25
宋承新.1999.径流还原计算的综合修正法.水文,(2):46～48
宋松柏,蔡焕杰,粟晓玲.2005.专门水文学.杨凌:西北农林科技大学出版社
宋松柏,康艳,荆萍.2008.水文频率曲线参数优化估计研究.西北农林科技大学(自然科学版), 36(4):193~198
宋松柏,康艳.2008.3种智能优化算法在设计洪水频率曲线适线法中的应用.西北农林科技大学学报:自然科学版, 36(2):205-209
宋晓波,王文圣,欧源.2007.非一致性年径流分析计算初探.黑龙江水专学报,34(2):4~61
王巧平.2003.天然年径流量系列一致性修正方法的改进.水利规划与设计, (2):38~40
王双银,向友珍,朱晓群,黄燕荣.2006.基于EXCEL的水文频率计算软件开发.西北农林科技大学学报(自然科学版), 34(4): 113~116
王双银,向友珍.2004.水文统计学.西北农林科技大学(校内教材)
王双银,张鑫.2005.水资源评价.西北农林科技大学(校内教材):17~28
王文圣,金菊良,李跃清.2007.水文随机模拟进展.水科学进展, 18(5): 768~771
王孝礼,胡宝清,夏军.2002.水文时序趋势与变异点的R/S分析法.武汉大学学报(工学版), 35(2):10～12
夏军,穆宏强,邱训平.2001.水文序列的时间变异性分析.长江职工大学学报, 18(3):1～4
夏军.2004.水问题的复杂性与不确定性研究与发展.北京:中国水利水电出版社
肖宜,夏军,申明亮,胡向阳.2001.差异信息理论在水文时间序列变异点诊断中的应用.中国农村水利水电, (11):28~30
谢平,陈广才,韩淑敏,武方圆.2006.从潮白河年径流频率分布变化看北京市水资源安全问题.长江流域资源与环境, 15(6):713～717.
谢平,郑泽全.2003.水文频率计算有约束加权适线法.武汉水利电力大学学报, 33(1):49~52
谢平,陈广才,夏军.2005a.变化环境下非一致性年径流序列的水文频率计算原理.武汉大学学报(工学版), 38(6):6~15
谢平,陈广才,李德,朱勇.2005b.水文变异综合诊断方法及其应用研究.水电能源科学,23(2):11~14
谢平,陈广才,夏军.2005c.基于趋势分析的非一致性年径流序列频率计算方法.中国水利学会第二届青年科技论坛论文集
谢平,陈广才,夏军.2005d.基于降雨径流关系分析的非一致性年径流序列频率计算方法.中国水利学会第二届青年科技论坛论文集.
谢平,陈广才,雷红富,李晶.2007a.论变化环境下的地表水资源评价方法.水资源研究,28(3):1~3
谢平,朱勇,陈广才,雷红富,李晶.2007b.考虑土地利用/覆被变化的集总式流域水文模型及应用研究.山地学报,25(3):257~264
谢平,陈广才,雷红富,韩淑敏,夏军,朱勇.2009.变化环境下地表水资源评价方法.北京:科学出版社,135~149.
谢云.1999.模拟退火算法的原理及实现.高等学校计算数学学报, (03): 212-218.
熊立华,周芬,肖义,郭生练.2005.水文时间序列变异点分析的贝叶斯方法.水电能源科学, 23(2):11～14
秀峰,阮本清.2001.对权函数法求概率分布参数的讨论.水利学报, (11): 34~40
薛联青,崔广柏,陈凯麒.2002.非平稳时间序列的动态水位神经网络预报模型.湖泊科学,14(1):19~24
杨珂玲,韩慧芳.2006.两混合正态分布的参数估计方法.黄冈师范学院学报,26(3):16~20.
杨士荣,戴少辉,王英凯.1997.山西省大中型水库年径流变化趋势分析.山西水利科技, 3(17):37~42
叶守泽,詹道江.2000.工程水文学.北京:中国水利水电出版社:173～200
袁超.2008.渭河流域主要河流水文干旱特性研究[硕士学位论文].杨凌:西北农林科技大学
张士峰.2004.混合正态分布参数极大似然估计的EM算法.飞行器测控学报,23(4):47~52
张学成,王玲,高贵成,郭喜有.2001.黄河流域降雨径流关系动态分析.水利水电技术,32(8):1~5
张一驰,周成虎,李宝林.2005.基于Brown~Forsythe检验的水文序列变异点识别.地理研究, 24(5):741~748
赵碧波,吕书君.2000.论陕西关中地区的水资源开发利用.地下水, 22(3):93~95
郑泽全,谢平,蔡伟.2001.小波变换在非平稳水文时间序列分析中的初步应用.水电能源科学,19(3): 49～51
周芬.2005.Kendall检验在水文序列趋势分析中的比较研究.人民珠江,(2):35~37
周芬,郭生练,肖义,徐长江.2003.P-Ⅲ型分布参数估计方法的比较研究.水电能源科学,21(3):10~13.
周光辉,王铁牛.2006.洪水频率计算新方法探讨.湖南水利水电, (2): 58~59
周至革.2001.一种计算随机变量函数均值和标准差的方法.机械强度,23(1):107~110
周忠学.2007.陕北黄土高原土地利用变化与社会经济发展关系及效益评价[硕士学位论文].西安:陕西师范大学
WMO&UNESCO. 2001.水资源评价.国家能力评估手册.李世明等翻译.郑州:黄河水利出版社.
Ramachandra Rao, Khaled H. Hamed. 2000. Flood Frequency Analysis. United States of American Alila Y, Mtiraoui A. 2002.Implications of heterogeneous flood-frequency distributions on traditional stream-discharge prediction techniques. Hydrological Processes, (16): 1065~1084
Arie B Z, Benjamin A. 1998. Distribution of Occurrence Days of Runoff Peaks in Israel and its Shifts for High~Discharge Events. Water Resources Management, (12):313~324
Aziz A, Burn D H. 2006. Trends and variability in the hydrological regime of the Mackenzie River Basin. Journal of Hydrology, 319(1): 282~294
Burn H D, Hag Elnur M A. 2002. Detection of hydrologic trends and variability. Journal of Hydrology, (225):107~122
Cohen A C. 1967. Estimation in mixtures of two normal distribution. Technometrics, 9(1):15~28
Fiorentino M, Arora K, Singh V P. 1987. The tow-component extreme value distribution for flood frequency anaylsis: derivation of a new estimation method. Stochasitic Hydrology and Hydraulics, (1): 199~208
Hamed K H. 2008. Trend detection in hydrologic data: The Mann–Kendall trend test under the scaling hypothesis. Journal of Hydrology: (13):350~363
Jiang S. and Kececioglu D. 1992. Maximum likelihood estimates from censored data for mixed Weibull distribution. IEEE Transactions on Reliability, 41(2):248-255
Kundzewicz Z W. 2004. Change detection in hydrological records-a review of the methodology. Hydrological Sciences, 49(1):7~19
Leytham K M. 1984. Maximum likelihood estimates for the parameters of mixed distributions. Water Resources Research, 20(7):896~902
Mann H B. 1945. Non-parametric test against trend. Econometric, (13):245~259
Moran PAP. 1959. The theory of storage. Methuen: London Perreault L, Bernier J,Bobee B,Parent E. 2000. Bayesian change-point analysis in hydrometeorological time series. Part1. Comparison of change-point models and forecasting, Journal of Hydrology, (235):221~241
Perreault L, Bernier J, Bobee B, Parent E. 2000. Bayesian change-point analysis inhydrometeorological time series. Part1. The normal model revisited, Journal of Hydrology, (235): 221~241
Peter J M. 2001. Evaluation of Mixed-Population Flood-Frequency Analysis. Journal of Hydrological Engineering: 6(1): 62~70
Rossi J, Fiorentino M, Versace P. 1984. Tow-component extreme value distribution for flood frequency analysis. Water Resources Research, 20(7):847~856
Singh K P, Sinclair R A. 1972. Tow~distribution method for flood frequency analysis. Journal of Hydraulics Division, ASCE98, No.HY1:29~44
Singh K P. 1987. A versatile flood frequency methodology. Water International, 12(3):139~145
Singh K P, Nakashima M. 1981. A new methodology for flood frequency analysis with objective detection and modification of outliers/inliers. Illinois State Water Survey.Report No.272:145.
Singh V P, Wang S X, Zhang L. 2005. Frequency analysis of nonidentically distributed hydrologic flood data. Journal of Hydrology, (307): 175~195
Waylen P, Woo M K.1982. Prediction of annual floods generated by mixed processes. Water Resources Research, 18(4):1283~1286
Woo M K, Thorne R. 2002. Comment on‘Detection of hydrologic trends and variability’by Burn, Journal of Hydrology, (255):107~122
陈望春,赵立锋.2002.权函数法在P-Ⅲ型分布中的应用.浙江水利科技,4:46~47
陈晓林,汪四水.2007.一类混合正态分布参数估计的EM算法和数据扩张.苏州大学学报(自然科学版), 23(3):13~18
崔维成,徐向东,邱强.1998.一种快速计算随机变量函数均值与标准差的新方法.船舶力学,(6):50~60
丁晶,邓育仁.1988.随机水文学.成都:成都科技大学出版社:64~110
董亚非,李兆元. 1995.关中地区气候资源特征分析.陕西气象,6:34~36
关中平原水资源考察小组.1990.关中平原水资源考察报告.杨凌:西北农业大学
韩小雷.2008.粒子群—模拟退火融合算法及其在函数优化中的应用.武汉:武汉理工大学出版社
胡四一.2006.水文学及水资源的学科前沿和研究要求.水文, 26(3): 1~4
金光炎.1993.水文水资源随机分析.北京:中国水利水电出版社
金光炎.1999.水文频率分析述评.水科学进展,10(3): 319~327
金光炎.2004.洪水频率分布模型的类别、比较和选用.安徽水利科技, (2): 36~38
雷红富,谢平,陈广才,李晶.2007.水文序列变异点检验方法的性能比较分析.水电能源科学, 25(4):36~40
李宏伟,宋松柏.2009.蚁群算法在水文频率曲线参数计算中的应用.人民黄河, 31(4):38-40
李立民,李玉波.2007.西安至宝鸡铁路客运专线可行性研究报告.西安:铁道第一勘察设计院
李述山,闫保英.1999.两正态混合分布的一种参数统计方法.山东轻工业学院学报, 13(2):73~77
李艳,陈晓宏,张鹏飞.2007.基于相关积分法的水文时间序列变异性度量.水文, 27(3):6~9
李志,曹明明.2003.关中地区水资源问题及对策研究.水土保持学报, 17(3):159~161
梁国文,佟伟.2008.水文频率曲线与水文随机过程.黑龙江水利科技, 5(36):162
梁艳琴.1994.用现状还原年径流模拟径流系列的探讨.海河水利, (3) :39～41
刘昌民,成立.2000.黄河干流下游断流的径流序列分析.地理学报, 55(3):257～265
刘光祖.2000.概率论与数理统计.北京:高等教育出版社
刘俊民,孙勇.2006.关中地区水资源特性分析及可持续利用对策.人民黄河, 28(1):39~41
陆中央.2000.关于年径流量系列的还原计算问题.水文, 20(6):9~12
马志明,刘瑞元,习丽.2007.多个子总体混合分布的参数估计.西北民族大学学报(自然科学版),28(65):11~15
曲小红,石生新,荣丰涛.1999.降雨径流系列的一致性分析.山西水利科技, (1):34~37
任树梅,李靖.2005.工程水文与水利计算.北京:中国农业出版社
桑燕芳,王栋.2008.水文时间序列周期识别的新思路与两种新方法.水科学进展,19(3):412~417
陕西省地方志编撰委员会.1999.陕西省志.第十三卷.水利志.西安:陕西人民出版社:65~69
沈宏.2003.天然径流还原计算方法初步探讨.水利规划与设计, (3):15~18
史鉴,陈兆丰,邢大韦.2002.关中地区水资源合理开发利用与生态环境保护.郑州:黄河水利出版社
史卫东.2001.水文资料系列的一致性分析.甘肃水利水电技术, 37(1):22~25
宋承新.1999.径流还原计算的综合修正法.水文,(2):46～48
宋松柏,蔡焕杰,粟晓玲.2005.专门水文学.杨凌:西北农林科技大学出版社
宋松柏,康艳,荆萍.2008.水文频率曲线参数优化估计研究.西北农林科技大学(自然科学版), 36(4):193~198
宋松柏,康艳.2008.3种智能优化算法在设计洪水频率曲线适线法中的应用.西北农林科技大学学报:自然科学版, 36(2):205-209
宋晓波,王文圣,欧源.2007.非一致性年径流分析计算初探.黑龙江水专学报,34(2):4~61
王巧平.2003.天然年径流量系列一致性修正方法的改进.水利规划与设计, (2):38~40
王双银,向友珍,朱晓群,黄燕荣.2006.基于EXCEL的水文频率计算软件开发.西北农林科技大学学报(自然科学版), 34(4): 113~116
王双银,向友珍.2004.水文统计学.西北农林科技大学(校内教材)
王双银,张鑫.2005.水资源评价.西北农林科技大学(校内教材):17~28
王文圣,金菊良,李跃清.2007.水文随机模拟进展.水科学进展, 18(5): 768~771
王孝礼,胡宝清,夏军.2002.水文时序趋势与变异点的R/S分析法.武汉大学学报(工学版), 35(2):10～12
夏军,穆宏强,邱训平.2001.水文序列的时间变异性分析.长江职工大学学报, 18(3):1～4
夏军.2004.水问题的复杂性与不确定性研究与发展.北京:中国水利水电出版社
肖宜,夏军,申明亮,胡向阳.2001.差异信息理论在水文时间序列变异点诊断中的应用.中国农村水利水电, (11):28~30
谢平,陈广才,韩淑敏,武方圆.2006.从潮白河年径流频率分布变化看北京市水资源安全问题.长江流域资源与环境, 15(6):713～717.
谢平,郑泽全.2003.水文频率计算有约束加权适线法.武汉水利电力大学学报, 33(1):49~52
谢平,陈广才,夏军.2005a.变化环境下非一致性年径流序列的水文频率计算原理.武汉大学学报(工学版), 38(6):6~15
谢平,陈广才,李德,朱勇.2005b.水文变异综合诊断方法及其应用研究.水电能源科学,23(2):11~14
谢平,陈广才,夏军.2005c.基于趋势分析的非一致性年径流序列频率计算方法.中国水利学会第二届青年科技论坛论文集
谢平,陈广才,夏军.2005d.基于降雨径流关系分析的非一致性年径流序列频率计算方法.中国水利学会第二届青年科技论坛论文集.
谢平,陈广才,雷红富,李晶.2007a.论变化环境下的地表水资源评价方法.水资源研究,28(3):1~3
谢平,朱勇,陈广才,雷红富,李晶.2007b.考虑土地利用/覆被变化的集总式流域水文模型及应用研究.山地学报,25(3):257~264
谢平,陈广才,雷红富,韩淑敏,夏军,朱勇.2009.变化环境下地表水资源评价方法.北京:科学出版社,135~149.
谢云.1999.模拟退火算法的原理及实现.高等学校计算数学学报, (03): 212-218.
熊立华,周芬,肖义,郭生练.2005.水文时间序列变异点分析的贝叶斯方法.水电能源科学, 23(2):11～14
秀峰,阮本清.2001.对权函数法求概率分布参数的讨论.水利学报, (11): 34~40
薛联青,崔广柏,陈凯麒.2002.非平稳时间序列的动态水位神经网络预报模型.湖泊科学,14(1):19~24
杨珂玲,韩慧芳.2006.两混合正态分布的参数估计方法.黄冈师范学院学报,26(3):16~20.
杨士荣,戴少辉,王英凯.1997.山西省大中型水库年径流变化趋势分析.山西水利科技, 3(17):37~42
叶守泽,詹道江.2000.工程水文学.北京:中国水利水电出版社:173～200
袁超.2008.渭河流域主要河流水文干旱特性研究[硕士学位论文].杨凌:西北农林科技大学
张士峰.2004.混合正态分布参数极大似然估计的EM算法.飞行器测控学报,23(4):47~52
张学成,王玲,高贵成,郭喜有.2001.黄河流域降雨径流关系动态分析.水利水电技术,32(8):1~5
张一驰,周成虎,李宝林.2005.基于Brown~Forsythe检验的水文序列变异点识别.地理研究, 24(5):741~748
赵碧波,吕书君.2000.论陕西关中地区的水资源开发利用.地下水, 22(3):93~95
郑泽全,谢平,蔡伟.2001.小波变换在非平稳水文时间序列分析中的初步应用.水电能源科学,19(3): 49～51
周芬.2005.Kendall检验在水文序列趋势分析中的比较研究.人民珠江,(2):35~37
周芬,郭生练,肖义,徐长江.2003.P-Ⅲ型分布参数估计方法的比较研究.水电能源科学,21(3):10~13.
周光辉,王铁牛.2006.洪水频率计算新方法探讨.湖南水利水电, (2): 58~59
周至革.2001.一种计算随机变量函数均值和标准差的方法.机械强度,23(1):107~110
周忠学.2007.陕北黄土高原土地利用变化与社会经济发展关系及效益评价[硕士学位论文].西安:陕西师范大学
WMO&UNESCO. 2001.水资源评价.国家能力评估手册.李世明等翻译.郑州:黄河水利出版社.
Ramachandra Rao, Khaled H. Hamed. 2000. Flood Frequency Analysis. United States of American Alila Y, Mtiraoui A. 2002.Implications of heterogeneous flood-frequency distributions on traditional stream-discharge prediction techniques. Hydrological Processes, (16): 1065~1084
Arie B Z, Benjamin A. 1998. Distribution of Occurrence Days of Runoff Peaks in Israel and its Shifts for High~Discharge Events. Water Resources Management, (12):313~324
Aziz A, Burn D H. 2006. Trends and variability in the hydrological regime of the Mackenzie River Basin. Journal of Hydrology, 319(1): 282~294
Burn H D, Hag Elnur M A. 2002. Detection of hydrologic trends and variability. Journal of Hydrology, (225):107~122
Cohen A C. 1967. Estimation in mixtures of two normal distribution. Technometrics, 9(1):15~28
Fiorentino M, Arora K, Singh V P. 1987. The tow-component extreme value distribution for flood frequency anaylsis: derivation of a new estimation method. Stochasitic Hydrology and Hydraulics, (1): 199~208
Hamed K H. 2008. Trend detection in hydrologic data: The Mann–Kendall trend test under the scaling hypothesis. Journal of Hydrology: (13):350~363
Jiang S. and Kececioglu D. 1992. Maximum likelihood estimates from censored data for mixed Weibull distribution. IEEE Transactions on Reliability, 41(2):248-255
Kundzewicz Z W. 2004. Change detection in hydrological records-a review of the methodology. Hydrological Sciences, 49(1):7~19
Leytham K M. 1984. Maximum likelihood estimates for the parameters of mixed distributions. Water Resources Research, 20(7):896~902
Mann H B. 1945. Non-parametric test against trend. Econometric, (13):245~259
Moran PAP. 1959. The theory of storage. Methuen: London Perreault L, Bernier J,Bobee B,Parent E. 2000. Bayesian change-point analysis in hydrometeorological time series. Part1. Comparison of change-point models and forecasting, Journal of Hydrology, (235):221~241
Perreault L, Bernier J, Bobee B, Parent E. 2000. Bayesian change-point analysis inhydrometeorological time series. Part1. The normal model revisited, Journal of Hydrology, (235): 221~241
Peter J M. 2001. Evaluation of Mixed-Population Flood-Frequency Analysis. Journal of Hydrological Engineering: 6(1): 62~70
Rossi J, Fiorentino M, Versace P. 1984. Tow-component extreme value distribution for flood frequency analysis. Water Resources Research, 20(7):847~856
Singh K P, Sinclair R A. 1972. Tow~distribution method for flood frequency analysis. Journal of Hydraulics Division, ASCE98, No.HY1:29~44
Singh K P. 1987. A versatile flood frequency methodology. Water International, 12(3):139~145
Singh K P, Nakashima M. 1981. A new methodology for flood frequency analysis with objective detection and modification of outliers/inliers. Illinois State Water Survey.Report No.272:145.
Singh V P, Wang S X, Zhang L. 2005. Frequency analysis of nonidentically distributed hydrologic flood data. Journal of Hydrology, (307): 175~195
Waylen P, Woo M K.1982. Prediction of annual floods generated by mixed processes. Water Resources Research, 18(4):1283~1286
Woo M K, Thorne R. 2002. Comment on‘Detection of hydrologic trends and variability’by Burn, Journal of Hydrology, (255):107~122