东江流域径流演变规律研究.docxVIP

东江流域径流演变规律研究 在过去的20年中，由于世界气候变化的影响，特别是人类不合理的水和土资源的开发利用，该流域的水和土壤资源发生了变化，下游区的水量被带到上级，这是下游供水量不足的主要原因。与此同时, 伴随植被的退化与减少, 流域水土流失、洪涝灾害事件频繁发生, 给流域的生态安全造成极大威胁。因此, 本研究以东江流域为例, 对其进行径流分析, 了解其演变规律, 以期为下游及整个东江流域的水资源合理配置奠定坚实的基础。 1 江苏河流域分布 东江地处湿润热带地区, 位于珠江三角洲的东部, 是广东省重要的三大水系之一。地理坐标为东经113°52′~115°52′, 北纬22°33′~25°14′。东江的上游有两条支流, 分别为东源浔乌水和西源九曲河, 都发源于赣南。源头河水向南流进入广东省境内, 然后汇入枫树坝水库;在河源县有支流新丰江注入。东江的主要支流有新丰江、西枝江、增江等, 其干流流经的城市有河源、惠州、东莞等, 在东莞石龙经东江河网汇入狮子洋, 东江全长523km, 河道的平均坡降为0.39‰, 流域面积为3.32万km2, 占全珠江流域面积7.3%。东江年径流总量为280亿m3, 占珠江总径流量8.3%, 在博罗站测得多年年径流量只有229亿m3, 夏季最丰, 干流径流量占全年径流量的51.7%;冬季最少, 只占7.6%。流域内多年平均降雨量为1500~2400mm, 平均1795mm。东江流域蒸发分布趋势:西南多, 东北少, 流域年蒸发能力在980~1150mm, 平均1100mm, 其中博罗站年均蒸发能力为1150mm, 为流域最高值, 龙川站年均蒸发能力为980mm, 为流域最低值。东江流域径流年际变化不大, 年径流变差系数最高出现在博罗地区, 仅为0.256, 这对水资源的稳定利用极为有利。但年内变化较大, 使东江在年内受汛期和枯期两方面的影响。 2 数据与流分析 2.1 年际变化特征分析 本研究通过选取东江流域代表站博罗站的年平均径流量进行流域的年径流演变规律分析, 进而揭示东江流域径流的年际变化规律。分析内容包括:径流丰枯预测、变异性分析及趋势性分析。径流丰枯预测采用马尔科夫链;突变性分析采用累积距平曲线联合T检验法;趋势性分析采用Kendall秩次检验法模型法。 2.2 年平均流量为5年平均流量为5. 水文数据采用东江 (博罗站) 1963~2005年逐月流量观测资料, 通过模比系数Kp判别, 计算得出东江 (博罗站) 各年平均流量的丰枯程度 (表1) 。 3 结果与分析 3.1 步状态转移概率矩阵 本研究应用马尔科夫链的基本原理和方法来研究径流序列的变化规律, 从而预报其未来变化趋势。如果序列变量的状态是可数的, 假如有N个, 那么从状态i经过一步状态转移到状态j, 都有不同状态发生的可能, 若用Pij表示这种可能性, 则称Pij为一步状态概率。将这些概率依序排列成一个矩阵———称为1步状态转移概率矩阵, 用P表示为: 该矩阵具有两个特性: (1) 矩阵每个元素均非负; (2) 矩阵每行元素之和等于1。若序列在时刻t0处于状态i, 经过n步转移, 在时刻tn处于状态j, 这种转移的可能性的数量指标称为n步转移概率, 记为: 则n步转移概率矩阵记为: 一步状态转移概率矩阵可用来预报次年的定性状态, 只要统计出一步状态转移概率矩阵元素值, 就可由关系式P (n) =Pn计算出n步状态转移概率矩阵元素值, 由此可预报n年后的定性状态。根据式 (3) 计算, 结果见表2。 由表2可知, 东江流域2004年、2005年丰枯状态分别为偏枯、正常的可能性最大。 3.2 特征点的初步识别—径流变异点分析 目前, 对序列的变异性分析大多都运用了大量的统计公式, 方法较为复杂, 本研究运用累计距平曲线-滑动T检验联合识别法。用累计距平曲线法对变点做初步诊断, 找出变异点, 再利用滑动t检验法进行精确识别。 累积距平是一种常用的判断变化趋势的方法。对于水文样本序列X1, X2, ……Xn, 在某一时刻t的累积距平表示为: 点绘出St~t曲线, 可进行初步变点分析。滑动t检验法:对已知的年径流样本序列, 选定某一年份, 分别取其前和后相邻的连续n1和n2年的年径流量值计算统计量T值。 3.2.1 年平均流量的变异点 由图1可知, 博罗站年平均径流量累积距平曲线最小值出现在1974年, 最大值出现在2001年以后, 因资料长度限制无法确定最大值出现的年份, 故初步选取1974年为东江博罗站年平均径流量的变异点。 3.2.2 计算事件卷型 博罗站在1974年前取10年 (1965~1973年) 和后取10年 (1975~1984年) , 分别计算出x1=217亿m3, x2=268亿m3, S1=36.05, S2=70.3,