降雨资料一致性检验实例计算分析
2015-10-29王琼
王 琼
(佳木斯水文局,黑龙江 佳木斯154002)
0 前 言
水文资料系列的一致性是指该数据系列在进行统计特征值分析时其短系列长度样本的代表性与总体样本的代表性相近,也就是说资料各长度系列的统计特征相近,可以作为一个整体来使用,否则就应剔除或修正某一短系列资料。
实际工作中经常遇到雨量站迁移或流域气候环境演变而出现的降雨资料是否具有一致性的问题,该问题是资料能否按一个统一总体来进行统计特征计算分析的前提。
有关文献上介绍的常规分析资料系列一致性的方法主要是假定系列已不存在一致性,通过点绘各短系列资料的降雨~径流相关点据,最后将前期降雨点据修正到近期降雨点群上来。该方法实际上就是对前期降雨数据直接进行了一致性修正,而没有对其是否偏离了允许限度进行检验,方法结果虽合理,但工作量较大,对一些总体系列符合检验标准的数据做了重复劳动,降低了整体工作效率[1-3]。
本文利用假设检验中t 检验、F 检验方法,对不同短系列降雨数据的总体性进行检验。如果通过检验,说明资料系列可视为同一总体,即总体资料具有一致性,并与采用一致性修正法的资料总体的统计特征值进行对比分析,来判断该方法的有效性。
1 假设检验介绍
首先提出原假设,选择统计量,从总体中抽取样本找出统计量。统计量的分布常用正态分布、t 分布和F 分布。确定接受或舍弃区间。对原假设做出接受或舍弃的判断,利用算出的统计量与选定临界值比较。若统计量不属于接受区间,则舍弃;反之则接受[4-6]。
F 检验法通过比较两组数据的方差,判断精密度是否存在显著性差异。F 检验法计算公式为:
式中:x、y、m、n 分别代表系列样本数据和样本个数,且F 值服从第一自由度m -1 和第二自由度n -1 的F 分布。
根据F 分布表查得临界值,如果F <F表,表明两组数据没有显著差异。t 检验用t 分布理论判断系列间存在差异的概率。
即先假定两个总体平均数之间没有显著差异;然后计算统计量t 值,本文有两组样本,其统计量T值的计算公式为:
式中:x、y、m、n 分别代表两系列样本数据和样本个数,其自由度为m + n -2,查表得出t 临界值并与计算t 值进行比较,若计算︱t︱值小于临界值则表明二系列数据差异不显著,通过检验。
2 实例分析
红旗雨量站设立于1963年,1993年站址变动迁移15 km。降雨资料系列自1963—2011年共49 a,测站迁移前有30 a连续资料,迁移后有19 a连续资料,先对其迁站前后资料进行一致性分析。
利用Excel 进行系列数据的f、t 检验,Excel 的“数据分析”功能α =0.05,查表可知,F年降雨量=1.9<F临界值=2.11;︱t︱年降雨量= 0.2 <t临界值=2.01,数据系列各项统计值均小于临界值,通过检验,说明该站迁移前后降雨资料可以视为同一总体连续应用,也说明具有一致性。F - 检验成果表见表1,t 检验成果表见表2。
表1 F-检验成果表
表2 t 检验成果表
3 结 语
数据一致性分析是选择和利用水文资料进一步分析计算的基础,选择f 检验和t 检验进行资料一致性分析的优势是利用Excel 计算方便快捷,对优选参数的取舍可以多次试算,这在手工计算时要耗费很大的精力。
文中实例数据经一致性修正法进行修正后的统计量计算结果与原系列计算的统计量值基本一致(过程略)。说明经本次检验已经可以按一致性资料进行处理,不必再进行一致性修正[7-9]。
本文所述方法仅局限于以上两种假设检验均通过时各系列资料才可视为具有一致性的总体性,且资料系列不宜过长,资料应来源于同一总体,否则会出现通过伪检验的错误结论。文中所述方法仅供参考,实际工作时建议采用多种方法进行对比分析。
[1]杨城芳. 地表水资源与水文分析[M]. 北京:水利水电出版社,1990:26 -30.
[2]缪铨生. 概率与数理统计[M]. 上海:华东师范大学出版社,2005:15 -19.
[3]祁萌,谢俊莹,潘建波. 无资料地区陆地蒸发量的推求[J]. 华北水利水电学院学报,2007(04):36 -38.
[4]石鹏,李冬梅,张宇发. 利用降雨资料推求无资料地区的径流量[J]. 贵州气象,2007(04):58 -60.
[5]冯伟,张万军. 陆面蒸发估算的CRAE 模型探讨[J]. 中国农学通报,2005(08):98 -100.
[6]顾本文,戴丛蕊. 云南省水、陆面蒸发量计算方法的初步研究[J]. 中国农业气象,2004(02):45 -48.
[7]高文义,全永哲. 无资料地区水资源量计算方法的探讨[J]. 吉林水利,2003(10):78 -80.
[8]胡顺军,田长彦,周宏飞,李慧. 渭干河灌区陆面蒸发量估算[J]. 干旱区地理,2000(01):42 -45.
[9]晋华,赵喜云. 利用气象资料确定蒸发量的新方法[J].太原理工大学学报,1999(05):36 -40.