08。黄河河源区枯季径流模型的研究及应用.docVIP

黄河河源区枯季径流模型的研究及应用刘根生孙胜利苟万德（黄委会上游水文水资源局，甘肃 兰州 730030）摘要：通过对黄河河源区枯季径流演变规律进行分析研究，在此基础上开发研制适合于河源区的长期径流预报模型，利用该模型对黄河上游地区未来径流演变趋势进行预测预报，给水量调度和水资源合理配置提供依据。 关键词：枯季；径流模型；研究；应用 黄河河源区（唐乃亥站以上流域）是黄河的主要来水区，径流量约占兰州以上来水量的65%。但九十年代以来，黄河河源区主要产流区降水、来水持续偏少，湖泊、沼泽面积减少，流域内生态发生了较大的变化，水资源短缺问题显得更加突出。针对这种情况，有必要对黄河河源区枯季径流进行分析、研究，开发适合流域特点的径流预报模型，为有关部门提供及时、准确的长期径流预测预报服务，作为水量调度和水资源合理配置依据。 流域概况 黄河河源区唐乃亥以上流域，介于东经95°00′∽103°30′、北纬32°19′∽36°08′之间, 面积12.20万km2。流域内海拔大都在3000m以上，巴颜喀拉山、阿尼玛卿山、岷山等山脉分布其间，众多河川、盆地、丘陵相间。海拔4000m以上的高山顶颠，大多岩石裸露，而山麓却绿草成茵。境内海拔最高处为阿尼玛卿山主峰，高6282m，主要由冰川覆盖，冰川面积191.95 km2,冰川消融量约为唐乃亥以上年径流量的1.0%。 枯季径流数学模型 流域内径流的变化既受确定性因素的影响，又受到随机性因素的影响，但这种过程是错综复杂的。通常采用确定性的、线性模型来逼近系统的输出。 在流域条件基本不变的条件下，流域径流的响应与其在在流域上发生的变化有其确定的因素，通常采用常微分方程表示： 当初始条件、在t=0时的值已知时，方程有唯一解。 如果用等距采样的离散序列、，表示一个离散系统，则可以通过对上式的离散化处理，得出离散系统模型。 其定解初始条件、为已知值。此式就是离散线性系统的差分方程模型。 微分、差分模型的全部特性反映在模型的阶数n和系数、 上。若b0=0，则k时刻输出变量仅为k时刻以前各输入、输出变量的函数，与k时刻输入变量无关，这称为显式差分格式；否则，则为隐式差分格式。 考虑到差分模型对系统模拟的不精确，存在一个模型噪声，依照线性叠加原理，方程将改为 模型参数的确定 模型识别 对于任何一种系统模型，如果转化为线性系统数学模型，符合以下通用的形式，就可以实现参数最小二乘法则进行识别。既 式中：为系统输入矩阵，维；为系统输出向量，维；为系统参数向量，维；为模型噪声向量，维。是模型参数的个数，是系统的输入、输出变量的观测数据组数。 影响因子的择取 黄河河源区枯水期，降水主要以雪或雨加雪形式，径流补给主要以流域蓄水和地下水组成，但影响径流的因子很多，要把众多因子都包含在内，是不可能的也是没有必要的，在满足精度的前提下，希望在整个系统中选取尽量少的因子。在实际应用中，采取数学方法提取方差贡献大的因子，剔除作用小的因子，从而简化模型参数。 模型参数估计 模型参数采用最小二乘法估计，如果系统参差满足一个零均值的白噪声，则用识别的参数还原计算系统的输出估计值与已知值的误差平方和达到最小。相应的参数就称为最小二乘法估计值。其最小二乘法准则的目标函数为 根据最小二乘法准则，推导线性系统数学模型的参数计算公式为： 依据极值理论，达到极小值的位置在处，既 则 上式称为线性系统数学模型的正规方程组，容易求出它的解向量为 上式不随时间变化，因此，模型是时不变的。但实际上，由于流域响应是时变系统，各次径流发展的不同阶段，其运动规律是各不相同的，故模型的系数是变动的，它是时间的函数。为达到随时间变化，必须构造一个反馈器，以便能从径流的发展变化过程中不断地获取新的信息，修正值。 成果检验 模型模拟情况 根据建立的枯季径流模型，对1956-2000年黄河河源区枯水期径流进行了模拟检验，从模拟的45年枯水期径流情况来看，根据长期径流预报规范，以相对误差以不超过30%为依据，对45年枯水期径流进行模拟，总合格率达到97%。1956-2000年枯水期径流模拟情况 100% 100% 100% 87% 4.2 模型预报情况 利用模型对2001年11月至次年4月径流进行试预报，从预报结果看，11月至次年3月模型预报结果完全满足精度要求，4月份径流预报值偏大，详细结果见表。2001-2002年枯水期径流预报情况2001～2002年364 375 3 12 194 189 -2 1 122 155 27 2 139 147 6 3 171 203 19 4 243 359 48 5 结语 通过对黄河河源区长期径流未来演变趋势的分析和模拟，从模拟和预报的结果可以看出，11月到次年3月模型模拟和