长江宜昌站径流量、输沙量趋势变化及突变分析.pdf

第四篇 应用实践与行动措施 长江宜昌站径流量、输沙量趋势变化及突变分析’ 郭瑾李琼芳马颖夏自强邹振华 (河海大学水资源环境学院南京210098) 摘要本文以长江流域重要拉{}4站宜昌站为研究站。应用1882—2005年日流量和1950—2060年日奢沙量、 2601—2004年月平均音沙量、2(105年日平均悬移质j葑沙串等青科，借骂}Mann—Kendall非参数检验方法对年径 洗、汛期及非汛期径流、年输沙量、汛期厦非汛期输沙量的趋势变化和宪变进行了分析。趋势分析表明宜昌站年 径流量、汛期与非汛期径流减小照势非常明显，且年径流量自20世纪80年代波动幅度械小。年输沙t，汛期与 非汛期辅涉量减小趋势非常明显。宜昌站年径巍量在1926年发生突变，减小趋势要加明显。年输沙量减小趋势 在1993年发生突变，下降程度更加显著。 关键词长江 宜昌站趋势分析蹙变分析Mann—Kendall 1引言 长江是我国第一太河，有充沛的水量和丰富的物种资源，为长江流域经济社会的发展提供了必要的保 障。但随着长江流域水资源开发利用程度的提高和气候变化，长江的水沙特性发生了变化，从而对长江生态 系统的稳定性与健康产生了影响。本文以长江宜昌水文站为研究站，研究与分析宜昌站水沙特性的变化规 律。该内容是研究变化环境下的水文循环和泥沙循环的主要内容之一．能为长江生态系统的健康保护提供 依据。 2研究区域 宜昌站位于长江上、中游的交界处，是长江上游干流出1：3控制站，控制上游流域面积100．6万时，宜 昌断面上游6．0 km为葛洲坝水利枢纽，上游约44 km有三峡水利枢纽工程，上游水利工程的兴建，改变河 流的径流量和其原有的季节分配，与径流相关的若干环境因子也将随之改变，诸如输沙量、水力学特征，水 质、温度、水体的自净能力等。鉴于宜昌站的特殊地理位置，有必要对其水沙变化规律进行研究。 3研究方法 趋势是指在一段时间参数变化的方向(上升、下降或横向延伸)。突变是自然界一种正常的现象，即描 述某自然现象的时间系列从一种稳定态(或稳定持续的变化趋势)跳跃式地转变到另一种稳定态(或稳定 持续的变化趋势)的现象，它表现为该系列在从一个统计特征到另一个统计特征的急剧变化”J。 本文采用Mann—Kendall(M—K)非参数统计检验方法，对宜昌站的径流量及输沙量资料序列进行趋 势分析及突变检验。 假设有一时间序列如下：置．五，玛，?，J。其趋势检验统汁量公式是： s=∑∑sign(X．一葺) 其中，Bi印()为符号函数。当X。一盖，小于、等于或大于零时，sign(置一置)分别为一1、0或1；M—K 统计量公式5大于、等于、小于零时分别为： r(s一1)／／n(Ⅱ一1)(2n+5)／18 s0 扯{o一——s_o1(s+1)／Vn(n一1)(2n+5)／18 s0 Z为正值表示增加趋势，负值表示减少趋势。Z的绝对值在大于等于1．28、1．64、2．32时表示分别通过 ·基金项目：国家自然科学基金资助项目。 第一作者简介：郭瑾(1981一)，男，山西长治人。硕士．水文学及永资源专业。E—mail：gjl9810268@163．oom 第十丸章长江流域 了信度90％、95％、99％显著性检验。 当用M—K法来检测时间序列的突变时，假设有一时问序列如下：X。，五，蜀，?。以，构造一秩序列 I mi．m．表示置Xj(1句≤i)的样本累积数。定义dk：dt=∑m．(2≤‘《Ⅳ) 也均值以及方差定义如F： E池】：掣 4 Var[以】：业』磬!生盟(2≤k≤．v) 在时间序列随机独立假定下。定义统计量：眠：锱吣1,2“3．．∽ ~，VarL吼J 这里峨呈标准正态分布，给定一显著性水平‰。查JE态分布表得到临界值to，当l UFkII％I， 表明序列存在一个明显的增长或减少趋势，所有u‘将组成一条曲线c。．通过信度检验可知其是否具有明 显的变化趋势。 把以上方法应用到反序列中．重复上述计算过程，并使计算值乘以一l，得出叩。。Ⅷ。在图中表示为 岛，若cI和G的交点位于信度线之间，则此点可能就屉突变点的开始”7。 本文选用宜昌站1882～2005年的日流量、1950—2000年的日古沙量、2001—2004年月平均含沙量、 2005年日平均息移质输沙率等资料，分别计算得到宜昌站年径流、年输沙量、汛期与非汛期径流量和输沙 量时闸序列，并借助Mann—Kendnll(^I—K)非参数统计检验方法对这些时间序划进行了趋势和突变分析。 4径流M—K法成果及讨论 4．1年径流量 宜昌站1882—2005年的年径流序列资料的分析成果见图1。由图1可以看出宜昌年径流量呈下降趋势。 宜昌站M—K趋势检验统计值为一2．98，超过置信度99％时的I临界值下限一2，32，表明径流量减少趋势明