水电站水能特性预测方法及终端设备与流程

1.本发明属于水电站水能计算技术领域，具体涉及水电站水能特性预测方法及终端设备。

背景技术：

2.我国早期投产水电机组受当时设计水平、材料加工、制造工艺、安装质量的限制，水轮机普遍存在转轮型号不匹配、效率低、运行区域狭窄(调节能力不足)、机组整体稳定性差等缺陷，严重影响了水电站的发电效益和运行安全性，因此，近年开始对这类早期投产的机组进行改造。
3.近年来，随着计算流体动力学(computational fluid dynamics，简称cfd)、有限元计算、模型测试及真机测试等技术的发展，我国的水轮机设计、材料、加工工艺等水平均有了长足进步，相比20年以前，在相同条件下设计的水轮机的能量、效率及稳定性指标均有了极大的提升。为了指导老水电机组的提质改造，有必要根据现有的水轮机指标如出力、效率等，对老的水电机组改造后预期的水力及能量特性进行预测，分析机组改造的可行性及经济效益，指导机组改造工作的开展，提升水电机组的综合效益。
4.目前对于水电站的水能计算采用的是电站综合出力系数k的方法，将水轮机的效率、流道的水力损失对水轮机出力的影响放入出力系数k值当中，根据一定的经验曲线确定不同水头、流量下的出力系数的取值，从而计算得到水电站的水力和出力特性。该方法具有计算简单的优点，但是出力系数的取值一般取定值或是根据经验曲线取值，造成计算结果不精确。

技术实现要素：

5.本发明提供了一种水电站水能特性预测方法及终端设备，解决了现有技术中出力系数的取值一般取定值或是根据经验曲线取值，导致计算所得水电站的水能结果不精确的技术问题。
6.本发明内容的第一方面公开了一种水电站水能特性预测方法，包括：
7.步骤1、获取水电站中每台机组对应的输水系统的水头损失系数；
8.步骤2、根据每台机组的水轮机的运转特性曲线，确定每个所述水轮机的水头与其最大流量之间的关系；
9.步骤3、根据每个所述水轮机的水头与其最大流量之间的关系、所述水头损失系数和研究时间内所述水电站的水情数据，迭代计算得到研究时间内每个时间段下每个水轮机的净水头和流量；
10.步骤4、根据每个水轮机的净水头和流量，确定所述水电站在所述研究时间内的预期电量。
11.优选地，所述步骤3具体包括：
12.步骤31、获取研究时间内所述水电站的逐小时时间段的水情数据，根据所述逐小
时时间段的水情数据，确定每个时间段内每个所述水轮机的毛水头，将所述毛水头设置为水头迭代的初始值；
13.步骤32、根据每个所述水轮机的水头与其最大流量之间的关系，确定插值得到的第一水头下每个所述水轮机的最大流量；
14.步骤33、根据所述最大流量和所述水头损失系数，计算每个所述水轮机的净水头；判断所述净水头与所述第一水头是否满足预设条件，如是将所述净水头记为所述第一水头，重复步骤32和步骤33，如否，则所述净水头为所求所述水轮机的净水头，所述净水头对应的流量为所求所述水轮机的流量。
15.优选地，所述步骤33具体包括：
16.步骤331、根据所述最大流量和水库的出库流量，确定每个所述水轮机的流量；
17.步骤332、根据每个所述水轮机的流量和所述水头损失系数，确定每个所述水轮机的净水头；
18.步骤333、判断所述净水头与所述第一水头是否满足预设条件，如是将所述净水头记为所述第一水头，重复步骤32和步骤33，如否，则所述净水头为所求所述水轮机的净水头，所述净水头对应的流量为所求所述水轮机的流量。
19.优选地，所述步骤332具体包括：
20.根据第一公式确定每个所述水轮机的净水头，所述第一公式为：
21.h＝h
g-α
t
q222.式中，h为净水头，hg为水头迭代的初始值，α
t
为局部水头损失系数，q为水轮机的流量。
23.优选地，所述预设条件如第二公式，所述第二公式为：
[0024][0025]
式中，h为净水头，h
(k)
为第一水头，ε为迭代误差。
[0026]
优选地，所述步骤4具体包括：
[0027]
步骤41、根据每个水轮机的净水头和流量，确定每个水轮机的效率；
[0028]
步骤42、根据所述净水头、流量和效率，确定每个所述时间内的水轮机的出力；
[0029]
步骤43、根据所述水轮机的出力，确定每个所述时间段内所述水电站的出力；
[0030]
步骤44、根据每个所述时间段内所述水电站的出力，确定所述水电站在所述研究时间内的预期电量。
[0031]
优选地，所述步骤41具体包括：
[0032]
步骤411、确定所述水轮机的净水头对应的流量数列和效率数列；
[0033]
步骤412、根据所述水轮机的流量在所述流量序列中的位置和所述效率数列，确定每个水轮机的效率。
[0034]
优选地，所述步骤411具体包括：
[0035]
根据第三公式确定所述水轮机的净水头对应的流量数列，所述第三公式为：
[0036]
[0037]
式中，q
h,i
为所述水轮机的净水头在第i条导叶等开度线上的流量，q
i,k
为第i条导叶等开度线上第k个水头对应的流量，q
i,k+1
为第i条导叶等开度线上第k+1个水头对应的流量，h为所述水轮机的净水头，hk和h
k+1
为根据所述水轮机的运转特性曲线确定的两个相邻的水头，所述h位于所述hk和所述h
k+1
之间；
[0038]
根据第四公式确定所述水轮机的净水头对应的效率数列，所述第四公式为：
[0039][0040]
式中，η
h,i
为所述水轮机的净水头在第i条导叶等开度线上的效率，η
i,k
为第i条导叶等开度线上第k个水头对应的效率，η
i,k+1
为第i条导叶等开度线上第k+1个水头对应的效率，h为所述水轮机的净水头，hk和h
k+1
为根据所述水轮机的运转特性曲线确定的两个相邻的水头，所述h位于所述hk和所述h
k+1
之间。
[0041]
优选地，所述步骤412具体为：
[0042]
根据第五公式确定每个水轮机的效率，所述第五公式为：
[0043][0044]
式中，η为所述水轮机的效率，q为所述水轮机的流量，q
h,i
和q
h,i+1
为所述流量序列中的相邻的两个流量，所述q位于所述q
h,i
和所述q
h,i+1
之间，η
h,i
为所述水轮机的净水头在第i条导叶等开度线上的效率，η
h,i+1
为所述水轮机的净水头在第i+1条导叶等开度线上的效率。
[0045]
本发明内容的第二方面公开了一种终端设备，包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述方法的步骤。
[0046]
本发明相对于现有技术，具有如下有益效果：
[0047]
本发明在水电站水能预测计算中根据水电站输水系统的设计，准确计算各台机组输水系统的水力损失系数。根据水情信息中上下游水位及水库出库流量，利用水轮机运转特性曲线、输水系统的水力损失系数采用迭代求解的方法准确计算各时间段的水轮机工作水头、流量、水轮机效率及水轮机出力，从而可以准确预测整个水电站机组改造后的水能特性。
附图说明
[0048]
图1为本发明实施例的水电站水能特性预测方法的流程图；
[0049]
图2为本发明实施例的水轮机运转特性曲线图；
[0050]
图3为本发明实施例的步骤3的详细流程图。
具体实施方式
[0051]
下文将结合附图以及具体实施案例对本发明的技术方案做更进一步的详细说明。应当了解，下列实施例仅为示例性地说明和解释本发明，而不应被解释为对本发明保护范围的限制。凡基于本发明上述内容所实现的技术均涵盖在本发明旨在保护的范围内。
[0052]
水电机组改造的可行性及经济性分析需要对机组改造后水电站的水力及出力特性进行一个较准确的计算，包括水轮机工作水头、水轮机出力，进而对改造后水电站的经济效益进行评价。根据水轮机出力公式n＝9.81qhη，式中n为机组出力，mw；η为水轮机效率，％；q为机组发电流量，m3/s；h为水轮机工作水头，m。对于已建成的水电站，其水情信息主要包括水库出库流量，水库上下游水位，要预测改造后水电机组水力及出力情况，需要根据水情信息计算出机组的工作水头和机组出力，从而计算整个电站的发电量。目前主要采用的方法为水能计算的简易模型，将水轮机出力公式简化为n＝kqh,其中定义k为电站综合出力系数k＝9.81η，一般的计算中取k为一固定值或是根据经验曲线取值，计算方法较粗略，且没有考虑不同工况下水轮机效率、流道内水力损失的差异，导致计算所得水电站的水能结果不精确。
[0053]
本发明提出了一种较为精确的水电站能量特性的预测方法，根据水情信息数据、输水系统水力损失系数计算、水轮机运转特性曲线插值计算得到准确的水轮机工作水头和出力，从而预测整个水电站机组改造后的水能特性。
[0054]
具体地，本发明内容的第一方面公开了一种水电站水能特性预测方法，如图1所示，包括：
[0055]
步骤1、获取水电站中每台机组对应的输水系统的水头损失系数。
[0056]
水电站输水系统的水头损失系数主要包括沿程水头损失和局部水头损失。其中，沿程水头损失系数αf可按公式(1)计算：
[0057][0058]
式中n为输水系统糙率，混凝土衬砌取0.013，钢衬取0.012；l为输水系统的长度，可根据水电站的设计图纸获得；r为输水系统断面的水力半径，根据水电站的设计图纸获得；a为输水系统断面的面积，根据水电站的设计图纸获得。
[0059]
水电站的局部水头损失主要为进出水口及分叉管的水头损失，根据《水电站调压室设计规范》(nb/t 35021-2014)可获得不同进、出水口体型及分叉管型式下的系数ξ，按公式(2)计算得到局部水头损失系数αm：
[0060][0061]
式中g为重力加速度取9.81，a为输水系统断面的面积，根据水电站的设计图纸获得。
[0062]
将各台机组对应的输水系统沿程损失系统和局部损失系数相加计算得到各台机组输水系统的总水头损失系数α
t
＝αf+αm。
[0063]
步骤2、根据每台机组的水轮机的运转特性曲线，确定每个水轮机的水头与其最大流量之间的关系。
[0064]
获取水轮机运转特性曲线如图2所示，沿不同导叶等开度线，获取水轮机水头、水轮机流量，水轮机效率，保证各导叶等开度线上水轮机水头取值点一致，则构建水轮机流量矩阵、水轮机效率矩阵。
[0065]
假设水轮机运转特性曲线上有n条导叶等开度线，各导叶等开度线上导叶开度为
{y1,y2,
…
,yn}，获取m个水轮机水头h，得到{h1,h2,
…
,hm}，则构建的水轮机流量矩阵q如公式(3)、水轮机效率η参数矩阵如公式(4)：
[0066][0067][0068]
式(3)中的q
ij
表示第i个导叶等开度线、第j个水头处水轮机的流量，式(4)中的η
ij
表示第i个导叶等开度线、第j个水头处水轮机的效率。
[0069]
将流量矩阵、效率矩阵输入预设的水轮机有功功率转换模型，即可获取水轮机有功功率矩阵。其中，水轮机有功功率转换模型如公式(5)：
[0070]nij
＝9.81hiq
ij
η
ij
ꢀꢀꢀꢀꢀꢀꢀꢀ
(5)
[0071]
则得到的水轮机的有功功率矩阵n
ij
可以表示为公式(6)的形式：
[0072][0073]
式(6)中，n
ij
为第i个导叶等开度线、第j个水头处水轮机的有功功率。
[0074]
步骤3、根据每个水轮机的水头与其最大流量之间的关系、水头损失系数和研究时间内水电站的水情数据，迭代计算得到研究时间内每个时间段下每个水轮机的净水头和流量，其流程图如图3所示，具体包括：
[0075]
步骤31、获取研究时间内水电站的逐小时时间段的水情数据，根据逐小时时间段的水情数据，确定每个时间段内每个水轮机的毛水头，将毛水头设置为水头迭代的初始值。
[0076]
本发明中的水情数据包括时间段t，上游水位z
up
，下游水位z
down
和水库的出库流量q
out
。
[0077]
则时间段t内水轮机的毛水头hg为：hg＝z
up-z
down
；
[0078]
设水头迭代的初始值h
(0)
为hg，即h
(0)
＝hg。
[0079]
步骤32、根据每个水轮机的水头与其最大流量之间的关系，确定插值得到的第一水头下每个水轮机的最大流量。
[0080]
由公式(3)中水头{h1,h2,
…
,hm}及最大流量{q
n,1
,q
n,2
,
…
,q
n,m
}的关系，插值得到第一水头下h
(k)
的水轮机最大流量
[0081]
步骤33、根据最大流量和水头损失系数，计算每个水轮机的净水头；判断净水头与
第一水头是否满足预设条件，如是将净水头记为第一水头，重复步骤32和步骤33，如否，则净水头为所求水轮机的净水头，净水头对应的流量为所求水轮机的流量，具体包括：
[0082]
步骤331、根据最大流量和水库的出库流量，确定每个水轮机的流量。
[0083]
计算并向上取整得到台数初步计算值n1,如果n1大于等于该电站总的装机台数，则取n1等于总装机台数。
[0084]
根据n1和水库的出库流量q
out
计算水轮机单机流量q为：
[0085]
步骤332、根据每个水轮机的流量和水头损失系数，确定每个水轮机的净水头，具体包括：
[0086]
根据公式(7)确定每个水轮机的净水头：
[0087]
h＝h
g-α
t
q2ꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(7)
[0088]
式中，h为净水头，hg为水头迭代的初始值，α
t
为水头损失系数，q为水轮机的流量。
[0089]
步骤333、判断净水头与第一水头是否满足预设条件，如是将净水头记为第一水头(设h
(k+1)
＝h)，重复步骤32和步骤33，反复迭代试算直至符合精度要求，如否，则净水头为所求水轮机的净水头，净水头对应的流量为所求水轮机的流量。
[0090]
其中的预设条件如公式(8)：
[0091][0092]
式中，h为净水头，h
(k)
为第一水头，ε为迭代误差，可设为0.01。
[0093]
步骤4、根据每个水轮机的净水头和流量，确定水电站在研究时间内的预期电量，具体包括：
[0094]
步骤41、根据每个水轮机的净水头和流量，利用公式(3)和公式(4)，采用线性插值确定每个水轮机的效率，具体包括：
[0095]
步骤411、确定水轮机的净水头对应的流量数列和效率数列，具体包括：
[0096]
判断水轮机净水头h在{h1,h2,
…
,hm}中的位置，如果在hk，h
k+1
之间，根据公式(9)计算出当前水头h对应的流量数列{q
h,1
,q
h,2
,
…
,q
h,n
}和根据公式(10)计算出当前水头h对应的水轮机效率数列{η
h,1
,η
h,2
,
…
,η
h,n
}：
[0097][0098]
式中，q
h,i
为水轮机的净水头在第i条导叶等开度线上的流量，q
i,k
为第i条导叶等开度线上第k个水头对应的流量，q
i,k+1
为第i条导叶等开度线上第k+1个水头对应的流量，h为水轮机的净水头，hk和h
k+1
为根据水轮机的运转特性曲线确定的两个相邻的水头，h位于hk和h
k+1
之间。
[0099]
水轮机的净水头对应的效率数列如公式(10)：
[0100]
[0101]
式中，η
h,i
为水轮机的净水头在第i条导叶等开度线上的效率，η
i,k
为第i条导叶等开度线上第k个水头对应的效率，η
i,k+1
为第i条导叶等开度线上第k+1个水头对应的效率，h为水轮机的净水头，hk和h
k+1
为根据水轮机的运转特性曲线确定的两个相邻的水头，h位于hk和h
k+1
之间。
[0102]
步骤412、根据效率数列和水轮机的流量在流量序列中的位置，确定每个水轮机的效率，具体包括：
[0103]
判断水轮机流量q在{q
h,1
,q
h,2
,
…
,q
h,n
}中的位置，如果在q
h,i
，q
h,i+1
之间，根据公式(11)计算出当前流量q下对应的水轮机效率η：
[0104][0105]
式中，η为水轮机的效率，q为水轮机的流量，q
h,i
和q
h,i+1
为流量序列中的相邻的两个流量，q位于q
h,i
和q
h,i+1
之间，η
h,i
为水轮机的净水头在第i条导叶等开度线上的效率，η
h,i+1
为水轮机的净水头在第i+1条导叶等开度线上的效率。
[0106]
步骤42、根据净水头、流量和效率，确定每个时间内的水轮机的出力。
[0107]
计算出各时间段的各台机组水轮机工作水头h
t
，水轮机流量q
t
，水轮机效率η
t
，计算各时段水轮机出力：n
t
＝9.81h
tqt
η
t
。
[0108]
步骤43、根据水轮机的出力，确定每个时间段内水电站的出力。
[0109]
根据各时段的机组台数n1可以求得水电站总的出力：
[0110]
步骤44、根据每个时间段内水电站的出力，确定水电站在研究时间内的预期电量e
t
：
[0111][0112]
本发明要求解决的技术问题是通过水电站的水情数据及水轮机的运转特性曲线对水轮机改造后的水电站的水能特性进行准确的预测。步骤包括：1)根据水电站输水系统的设计准确计算各台机组输水系统的水力损失系数；2)根据水轮机运转特性曲线构建水轮机的流量矩阵、效率矩阵和出力矩阵；3)基于水电站各时段的水情数据及水轮机运转特性曲线，迭代计算得到各时段下的水轮机的工作水头；4)根据水轮机运转特性曲线插值计算得到各时段的准确的水轮机引用流量和水轮机效率；5)计算各时段水轮机出力及预期电量。
[0113]
本发明在水电站水能预测计算中根据水电站输水系统的设计，准确计算各台机组输水系统的水力损失系数。根据水情信息上下游水位及水库出库流量，利用水轮机运转特性曲线、输水系统的水力损失系数采用迭代求解的方法准确计算各时间段的水轮机工作水头、流量、水轮机效率及水轮机出力，从而可以准确预测整个水电站机组改造后的水能特性。
[0114]
本发明内容的第二方面公开了一种终端设备，包括存储器、处理器以及存储在存储器中并可在处理器上运行的计算机程序，处理器执行计算机程序时实现上述方法的步骤。
[0115]
上述只是本发明的较佳实施例，并非对本发明作任何形式上的限制。虽然本发明已以较佳实施例揭露如上，然而并非用以限定本发明。任何熟悉本领域的技术人员，在不脱离本发明技术方案范围的情况下，都可利用上述揭示的技术内容对本发明技术方案做出许多可能的变动和修饰，或修改为等同变化的等效实施例。因此，凡是未脱离本发明技术方案的内容，依据本发明技术实质对以上实施例所做的任何简单修改、等同变化及修饰，均应落在本发明技术方案保护的范围内。

技术特征：

1.一种水电站水能特性预测方法，其特征是，包括：步骤1、获取水电站中每台机组对应的输水系统的水头损失系数；步骤2、根据每台机组的水轮机的运转特性曲线，确定每个所述水轮机的水头与其最大流量之间的关系；步骤3、根据每个所述水轮机的水头与其最大流量之间的关系、所述水头损失系数和研究时间内所述水电站的水情数据，迭代计算得到研究时间内每个时间段下每个水轮机的净水头和流量；步骤4、根据每个水轮机的净水头和流量，确定所述水电站在所述研究时间内的预期电量。2.如权利要求1所述的方法，其特征是，所述步骤3具体包括：步骤31、获取研究时间内所述水电站的逐小时时间段的水情数据，根据所述逐小时时间段的水情数据，确定每个时间段内每个所述水轮机的毛水头，将所述毛水头设置为水头迭代的初始值；步骤32、根据每个所述水轮机的水头与其最大流量之间的关系，确定插值得到的第一水头下每个所述水轮机的最大流量；步骤33、根据所述最大流量和所述水头损失系数，计算每个所述水轮机的净水头；判断所述净水头与所述第一水头是否满足预设条件，如是将所述净水头记为所述第一水头，重复步骤32和步骤33，如否，则所述净水头为所求所述水轮机的净水头，所述净水头对应的流量为所求所述水轮机的流量。3.如权利要求2所述的方法，其特征是，所述步骤33具体包括：步骤331、根据所述最大流量和水库的出库流量，确定每个所述水轮机的流量；步骤332、根据每个所述水轮机的流量和所述水头损失系数，确定每个所述水轮机的净水头；步骤333、判断所述净水头与所述第一水头是否满足预设条件，如是将所述净水头记为所述第一水头，重复步骤32和步骤33，如否，则所述净水头为所求所述水轮机的净水头，所述净水头对应的流量为所求所述水轮机的流量。4.如权利要求3所述的方法，其特征是，所述步骤332具体包括：根据第一公式确定每个所述水轮机的净水头，所述第一公式为：h＝h
g-α
t
q2式中，h为净水头，h
g
为水头迭代的初始值，α
t
为水头损失系数，q为水轮机的流量。5.如权利要求3所述的方法，其特征是，所述预设条件如第二公式，所述第二公式为：式中，h为净水头，h
(k)
为第一水头，ε为迭代误差。6.如权利要求1-5任一项所述的方法，其特征是，所述步骤4具体包括：步骤41、根据每个水轮机的净水头和流量，确定每个水轮机的效率；步骤42、根据所述净水头、流量和效率，确定每个所述时间内的水轮机的出力；步骤43、根据所述水轮机的出力，确定每个所述时间段内所述水电站的出力；
步骤44、根据每个所述时间段内所述水电站的出力，确定所述水电站在所述研究时间内的预期电量。7.如权利要求6所述的方法，其特征是，所述步骤41具体包括：步骤411、确定所述水轮机的净水头对应的流量数列和效率数列；步骤412、根据所述效率数列和所述水轮机的流量在所述流量序列中的位置，确定每个水轮机的效率。8.如权利要求7所述的方法，其特征是，所述步骤411具体包括：根据第三公式确定所述水轮机的净水头对应的流量数列，所述第三公式为：式中，q
h,i
为所述水轮机的净水头在第i条导叶等开度线上的流量，q
i,k
为第i条导叶等开度线上第k个水头对应的流量，q
i,k+1
为第i条导叶等开度线上第k+1个水头对应的流量，h为所述水轮机的净水头，h
k
和h
k+1
为根据所述水轮机的运转特性曲线确定的两个相邻的水头，所述h位于所述h
k
和所述h
k+1
之间；根据第四公式确定所述水轮机的净水头对应的效率数列，所述第四公式为：式中，η
h,i
为所述水轮机的净水头在第i条导叶等开度线上的效率，η
i,k
为第i条导叶等开度线上第k个水头对应的效率，η
i,k+1
为第i条导叶等开度线上第k+1个水头对应的效率，h为所述水轮机的净水头，h
k
和h
k+1
为根据所述水轮机的运转特性曲线确定的两个相邻的水头，所述h位于所述h
k
和所述h
k+1
之间。9.如权利要求8所述的方法，其特征是，所述步骤412具体为：根据第五公式确定每个水轮机的效率，所述第五公式为：式中，η为所述水轮机的效率，q为所述水轮机的流量，q
h,i
和q
h,i+1
为所述流量序列中的相邻的两个流量，所述q位于所述q
h,i
和所述q
h,i+1
之间，η
h,i
为所述水轮机的净水头在第i条导叶等开度线上的效率，η
h,i+1
为所述水轮机的净水头在第i+1条导叶等开度线上的效率。10.一种终端设备，包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序，其特征是，所述处理器执行所述计算机程序时实现如权利要求1至9任一项所述方法的步骤。

技术总结

本发明公开了一种水电站水能特性预测方法及终端设备，其中方法包括获取水电站中每台机组对应的输水系统的水头损失系数；根据每台机组的水轮机的运转特性曲线，确定每个水轮机的水头与其最大流量之间的关系；根据每个水轮机的水头与其最大流量之间的关系、水头损失系数和研究时间内水电站的水情数据，迭代计算得到研究时间内每个时间段下每个水轮机的净水头和流量；根据每个水轮机的净水头和流量，确定水电站在研究时间内的预期电量。本发明可根据水情信息，利用水轮机运转特性曲线、输水系统的水力损失系数采用迭代求解的方法准确计算各时间段的水轮机工作水头、流量、效率及水轮机出力，从而准确预测整个水电站机组改造后的水能特性。的水能特性。的水能特性。