水电站径流调节的数值解法.docx

2010年第4期东北水利水电水利科研

·49·

[文章编号]1002-0624(2010)04-0049-04

水电站径流调节的数值解法

张晓波，郑雄伟，陈昌军

(浙江省水利水电勘测设计院，浙江杭州310002)

[摘要]文章推导了水电站径流调节的多种数值计算方法，并从数学原理和工程实用方面对各数值解法进行了比较，推荐简便、高效的不动点迭代法应用于径流调节计算。

[关键词]水电站；径流调节；不动点迭代法[中图分类号]TV74

随着计算机技术的广泛引用，水电站径流调节计算方法已由传统的列表试算法、图解法转向数值计算方法。大连理工大学于1981年首次提出采用改进欧拉法进行数值求解，受此启发，笔者将水电站径流调节计算的数值解法进行拓展，分别采用二分法、改进欧拉法、不动点迭代法、牛顿迭代法对水电站径流调节进行数值求解推导，并从迭代效率、计算精度、实用性等方面进行比较，推荐实用、简便的不动点迭代法应用于径流调节计算。

1计算原理

水电站径流调节基于以下基本方程1:

(1)(2)N(t)=koq(Z上(t)-Z

(1)(2)

式(1)为水量平衡方程，Q(t)为入库流量；q(t)为下泄流量；V(t)为水库库容。为便于公式表示，采用三次多项式拟合库水位和库容关系[1],如式(3):

V(t)=V(Z上(t))

=a?+a?Z上(t)+a?z2上(t)+aZ上(t)(3)

需要指出，地形复杂的水库其水位库容曲线很难用式(3)完整地表达，很多时候需采用分段函数拟合，程序计算应采用插值法。

式(2)为水电站出力计算公式，N(t)为出力；Z上为库水位；Z下为电站尾水位；ko为出力系数，ko=9.81η,η为机组出力效率系数，可取一定值；DH(t)为水头损失，为便于公式表示，表达为q(t)的二次函数，如式(4):

[文献标识码]A

DH(t)=k?g2(t)(4)

需要指出，实际发电中，机组出力效率是变化的，水头损失也受出力不均匀的影响，式(4)不能反映这些影响，程序计算应考虑这些具体问题。

电站尾水位也可拟合为q(t)的n次函数，如式(5):

Z下(t)=qg(t),n=0,1,2,…(5)

联解式(1),(2)及相关公式，显然难以直接求解。对该非线性方程，列举以下方法求解之。

2数值解法

2.1改进欧拉法

(1)数学原理。欧拉法迭代原理如图1,对于方程y=y(x),可以从初始点P(xo,yo)出发，沿该点切线方向推进到x=x?上一点P,然后再从P依该点切线方向推进到x=x?上一点P?依次类推囚，记y(x)=f(xx,y(x?),有求解关系：

y(x+1)=y(xA)+△xf(x,y.)(6)欧拉法相当粗略，多次推进后，欧拉法计算的

Pn显著偏离实际曲线，因此提出了改进欧拉法2:y(x+)=y(xa)+△xf(v?ya)(7)

(2)应用求解。联立式(1),(2),并代入辅助方程式(3),(4),(5),使得未知数q(t)和Z上(t)显现(以下Z上(t)简化为Z(t)),式(1)可变形为：

(9)

·50·

水利科研东北水利水电2010年第4期

图1欧拉迭代法

为水库面积，为库水位的函数。

式(2)可变形为：

(10)

将式(10)代入式(9),可得：

(11)

式(8)为非线性的微分方程。大连理工大学陈守煜教授将式(11)看作常微分方程的初值问题求解13,采用改进欧拉法求解，陈守煜教授对式(11)的改进欧拉求解公式为：

2(t+1)=Z(t)+△f(t,z(t))

(12)

应用式(12),(13)可以对水电站径流调节直接求解。q=f1(Z)无法直接求解，需要迭代法求解。

(3)解法评价。改进欧拉法使得所求的末时段的Z(t)显式表示，无需迭代，方法简便。但是，正是由于改进欧拉法是非线性方程的二阶近似，只有二阶精度。对于在步长△t较大时(如逐月径流调节),误差较大，陈守煜教授在其文献[3]中的示例，已经说明了该法精度不及传统的迭代法，因为传统的迭代法可控制误差精度，欧拉法正是以牺牲精度来实现其简单的计算方法。而且，发电