首页 / 专利库 / 水能 / 水轮机 / 基于图论的水电站引水发电系统小波动稳定性分析方法

基于图论的电站引水发电系统小波动稳定性分析方法

阅读:567发布:2023-03-09

专利汇可以提供基于图论的电站引水发电系统小波动稳定性分析方法专利检索,专利查询,专利分析的服务。并且本 发明 公开了一种基于图论的 水 电 站引水发电系统小 波动 稳定性 分析方法,该方法首先 指定 系统的状态变量顺序,然后基于图论列出状态变量方程,将状态变量方程组合成微分方程组,通过矩阵变换得到系统 状态方程 ,计算其中矩阵的特征值,进而判断系统的稳定性。本发明基于水电站引水发电系统小波动稳定分析理论,寻求状态方程的规律性特征,建立水电站引水发电系统小波动稳定性分析图论原理,并通过指定状态变量顺序以及矩阵变换的方法,给出状态变量系数矩阵通用的求解方法,该方法可以准确快速得到状态变量的系数矩阵,分析系统的稳定性特征及优化系统的参数。,下面是基于图论的电站引水发电系统小波动稳定性分析方法专利的具体信息内容。

1.一种基于图论的电站引水发电系统小波动稳定性分析方法,其特征在于包括以下步骤:
步骤一,指定系统的状态变量顺序;
步骤二,基于图论列出状态变量方程;
步骤三,将状态变量方程组合成矩阵方程;
步骤四,矩阵方程通过矩阵变换得到系统状态方程
步骤五,计算系统状态方程中矩阵的特征值;
步骤六,由特征值判断系统的稳定性。
2.根据权利要求1所述的基于图论的水电站引水发电系统小波动稳定性分析方法,其特征在于:所述步骤一中状态变量的顺序为 μ、q、z,其中 为机组转速变化相对值, 为机组开度变化相对值, 为流量变化相对值, 为调压室水位
或其他变量变化相对值。
3.根据权利要求1所述的基于图论的水电站引水发电系统小波动稳定性分析方法,其特征在于:所述步骤二中系统状态方程包括列机组运动方程、调速器方程、管道段水流动方程、调压室连续方程或其他变量方程。
4.根据权利要求3所述的基于图论的水电站引水发电系统小波动稳定性分析方法,其特征在于:所述机组运动方程为:
其中 为第i台机组相对转速,μi为第i台机组开度变化相对值,qi为第i台机组流量变化相对值,xi为系统负荷相对波动值,Spi为第i台机组负荷自调节系数,Tmi为第i台机组加速时间常数,S8i、S9i、S10i为已知水轮机特性参数,t为时间。
5.根据权利要求3所述的基于图论的水电站引水发电系统小波动稳定性分析方法,其特征在于:所述调速器方程为:
其中bpi、bti、Tdi、Tni分别为第i(i=1~N)台机组对应调速器的永态调差率、暂态转差系数、缓冲时间常数、微分时间常数。
6.根据权利要求3所述的基于图论的水电站引水发电系统小波动稳定性分析方法,其特征在于:利用图论原理直接给出小波动稳定分析的水道系统水流动力方程:
其中L为管段长度,A为管段面积,α为对应管段水头损失系数,Q为对应管段内流量,Hu,Hd分别是该段管道的进口水头和出口水头。
7.根据权利要求1所述的基于图论的水电站引水发电系统小波动稳定性分析方法,其特征在于:所述矩阵方程为
其中:E、F、R为系统相关的系数矩阵,F=EA,R=EBx。
8.根据权利要求2所述的基于图论的水电站引水发电系统小波动稳定性分析方法,其特征在于:所述系统状态方程为
其中: 为微分算子,Y为状态变量,A和B为系数矩阵,x为相对扰动量。

说明书全文

基于图论的电站引水发电系统小波动稳定性分析方法

技术领域

[0001] 本发明属于水利水电工程,特别是一种基于图论的水电站引水发电系统小波动稳定性分析方法。

背景技术

[0002] 大型水电站水轮发电机组的运行稳定性对机组本身和整个电系统的安全稳定都有较大的影响,确保水电站水力机械系统的运行稳定性及具备良好的调节品质是水电站
引水发电系统设计中一项重要的课题。随着特高压直流输电技术的应用以及出现特殊工
况,水电站运行中则可能会出现孤岛运行方式,特别是西部一些特大型的引水式电站,其水道系统较长,并设置了大型调压室,小波动问题较为突出。水利水电规划设计总院在《水电站输水发电系统调节保证设计专题报告编制暂行规定》中明确提出对于小波动计算原则上
应按照孤网运行条件进行。
[0003] 水轮机的调节系统是水力、机械及电气相互耦合的控制系统,属于非线性调节系统,涉及到的领域及理论较多。对于该系统的研究主要采用理论分析和数值模拟两种方法。
采用数值模拟方法虽然能够尽量多的考虑到系统的非线性因素,但数值模拟存在数学模型
复杂,编程困难以及数值误差影响,特别是在小扰动的情况下,会对计算结果产生较大影
响,另外,数值模拟基于时域,因此不方便对系统的参数及调速器的参数进行敏感性分析。
对于非线性系统,采用现有的自动控制理论分析系统的稳定性十分困难,因此,基于状态空间法的水轮机调节系统线性化模型及刚性水体模型仍然是目前分析水电站引水发电系统
小波动稳定性的一种重要方法。前人在该方面的研究成果较为丰富,对水电站的引水发电
系统的设计起到了重要的指导作用。但是随着水电站布置形式的日益复杂,在进行系统小
波动稳定分析时,模型阶数较高,解耦过程较为困难,且推导过程中容易出错,往往细微的差别就可能导致结论发生改变。因此如何准确快速得到状态方程的系数矩阵,对于分析复
杂系统的小波动稳定分析具有重要的意义。

发明内容

[0004] 针对现有技术存在的问题,本发明提出了一种基于水电站水力-机械系统小波动稳定分析理论,寻求状态方程的规律性特征,建立水电站引水发电系统小波动稳定性分析
图论原理,并通过指定状态变量顺序以及矩阵变换的方法,给出状态变量系数矩阵通用的
求解方法,该方法可以准确快速得到状态变量的系数矩阵,从而可以判定出系统的稳定性
特征。
[0005] 为了解决上述技术问题,本发明采用如下的技术方案:
[0006] 步骤一,指定系统的状态变量顺序;
[0007] 步骤二,基于图论列出状态变量方程;
[0008] 步骤三,将状态变量方程组合成矩阵方程;
[0009] 步骤四,矩阵方程通过矩阵变换得到系统状态方程;
[0010] 步骤五,计算系统状态方程中矩阵的特征值;
[0011] 步骤六,由特征值判断系统的稳定性。
[0012] 具体地,所述步骤一中状态变量的顺序为 μ、q、z,其中 为机组转速变化相对值, 为机组开度变化相对值, 为流量变化相对值, 为调压室水
位或其他变量变化相对值。
[0013] 具体地,所述步骤二中系统状态方程包括列机组运动方程、调速器方程、管道段水流动力方程、调压室连续方程或其他变量方程。
[0014] 具体地,所述机组运动方程为:
[0015]
[0016] 其中 为第i台机组相对转速,μi为第i台机组开度变化相对值,qi为第i台机组流量变化相对值,xi为系统负荷相对波动值,Spi为第i台机组负荷自调节系数,Tmi为第i台机组加速时间常数,S8i、S9i、S10i为已知水轮机特性参数,t为时间。
[0017] 具体地,所述调速器方程为:
[0018]
[0019] 其中bpi、bti、Tdi、Tni分别为第i(i=1~N)台机组对应调速器的永态调差率、暂态转差系数、缓冲时间常数、微分时间常数。
[0020] 具体地,利用图论原理直接给出小波动稳定分析的水道系统水流动力方程:
[0021]
[0022] 其中L为管段长度,A为管段面积,α为对应管段水头损失系数,Q为对应管段内流量,Hu,Hd分别是该段管道的进口水头和出口水头。
[0023] 具体地,所述矩阵方程为
[0024] 其中:E、F、R为系统相关的系数矩阵,F=E A,R=E B x
[0025] 具体地,所述系统状态方程为
[0026] 其中: 为微分算子,Y为状态变量,A和B为系数矩阵,x为相对扰动量
[0027] 上述稳定性分析方法的推算过程如下:
[0028] (一)系统稳定性特征及状态矩阵:
[0029] 系统状态方程表达式如下:
[0030]
[0031] 其中: 为微分算子,Y为状态变量,A和B为系数矩阵,x为相对扰动量。
[0032] 水电站水力一机械系统的小波动稳定性取决于上式中矩阵A的特征值λn=σn+iωn(n为状态变量个数),当A的所有特征值的实部均为负值,即σn<0,系统才是稳
定的,否则系统不稳定。当系统状态变量较少时,系数矩阵A可直接通过数学推导得出。当系统的状态变量较多时,状态方程解耦较为困难,则很难通过数学推导直接得出系数矩阵
A。
[0033] (二)状态矩阵变换:
[0034] 将系统状态矩阵表达式两边同时左乘一系数矩阵E,进行变换得如下表达式:
[0035]
[0036] 其中:E、F、R为系统相关的系数矩阵,F=E A,R=E B x。上式即为关于状态变量的常系数微分方程组。相对于系数矩阵A,矩阵E、F、R更容易直接得出,且具有一定的规律性,便于编写通用程序。当系数矩阵E、F、R得出后,只需在上两侧同时左乘系数矩阵E的逆矩阵,则可直接得出系统状态矩阵表达式。
[0037] 为了寻求系数矩阵E、F、R的规律性特征,便于编写通用程序,指定状态变量的顺序 μ、q、z。各变量意义如下: 为机组转速变化相对值; 为机组开度变化相对值; 为流量变化相对值, 为调压室水位(或其他变量)变化相对值。N
为机组总数;NQ为流量变量总数、W为调压室水位(或其他变量)的总数。因此状态变量的
总数M为2N+NQ+W,须联立M个微分方程。
[0038] (三)根据水电站机组台数列出机组运动方程,表达式如下:
[0039]
[0040] 其中:为第i台机组相对转速,μi为第i台机组开度变化相对值,qi为第i台机组流量变化相对值,xi为系统负荷相对波动值,Spi为第i台机组负荷自调节系数,Tmi为第i台机组加速时间常数,S8i、S9i、S10i为已知水轮机特性参数,t为时间。
[0041] (四)PID型调速器方程,表达式如下:
[0042]
[0043] 其中bpi、bti、Tdi、Tni分别为第i(i=1~N)台机组对应调速器的永态调差率、暂态转差系数、缓冲时间常数、微分时间常数。
[0044] (五)将(三)中的表达式两边对时间求导数,带入(四)中的表达式并消去对每台机组列方程并按照指定的顺序排列得以下调速器表达式:
[0045]
[0046] (六)管道段的水流动力方程,表达式如下:
[0047]
[0048] 其中L为管段长度,A为管段面积,α为对应管段水头损失系数,Q为对应管段内流量,Hu,Hd分别是该段管道的进口水头和出口水头。根据水流动力方程,可列出关于流量变量的微分方程。由于不同水电站的水道布置差别较大,很难像机组及调速器方程一样给
出统一的表达式,本发明借助有向线性图,描述水电站引水发电系统小波动稳定性分析图
论原理。具体见实施例
[0049] (七)其他状态变量的微分方程:
[0050] 小波动稳定分析的状态变量除了机组转速、开度、流量以外,有些还会涉及到调压室水位以及其他的变量等等,如图2中的调压室水位,则需要调压室的连续方程的表达式,具体表达式见实施例。
[0051] (八)状态方程系数矩阵构建:
[0052] 将所有方程写成(二)中的矩阵方程的形式:
[0053]
[0054] 式 中 ei,i= 1、 ei+N,i=Tni+Tdi+Tni×Tdi×fi,i、ei+N,i+N=(bti+bpi)Tdi+Tni×Tdi×fi,i+N、Ei+N,i+2N=Tni×Tdi×fi,i+2N、fi+N,i=-1、fi+N,i+N=-bpi,i=1~N(N为机组数)。该部分为机组及调速器方程对应的系数表达式,对于调速器采用PID及PI(Tn=0)调节模式的水电站,其表达式不变,即系数矩阵上
部分系数表达式不变。e2N+i,j=BXCY、f2N+i,i=-S6,i、f2N+i,N+i=-S7,i、f2N+i,j=DXCY、f2N+i,2N+i=f′2N+i,2N+i-S5,i(i=1~N、j=2N+1~NQ),该部分为基于图论的管道水流动力方程对应的系数表达式,其中B,C和D的表达式以及如何得到在实施例中均作了详细说明,须指出的是,由于水流动力方程中用状态变量表示机组水头相对变化值ζ,因此引入了系数f2N+i,i=-S6,i、f2N+i,N+i=-S7,i以及f2N+i,2N+i=f′ 2N+i,2N+i-S5,i,f′2N+i,2N+i在式f2N+i,j=DXCY已赋值。系数矩阵中最后面一部分为除机组变量 以外的状态变量对应的系数,对于
其他流量部分,仍然可以利用图论原理直接列出水流动力方程,其他的变量则根据相关的
物理意义给出相应的微分方程,再将各项按照指定的顺序排列,分别赋给系数矩阵中对应
的系数,矩阵中未赋值的系数均为0。则系统的系数矩阵E、F、R已经得出,只需在上式两侧同时左乘系数矩阵E的逆矩阵,则可直接得出判别系统稳定的系数矩阵A。
附图说明:
[0055] 图1为本发明中引水发电系统小波动稳定性分析原理图;
[0056] 图2为本发明中具体实施方式所用水电站引水发电系统布置简图;
[0057] 图3为本发明中算例分析得出的系统稳定图;
[0058] 图4为本发明中算例分析得出系统发散图。具体实施方式:
[0059] 图1为本发明中引水发电系统小波动稳定性分析原理图,步骤为(1)指定状态变量顺序、(2)列出所有机组运动方程、(3)列出所有调速器方程、(4)基于图论原理列出所
管道系统方程、(5)列出其他状态变量方程、(6)将各方程组合成 形式、(7)
两侧同时左乘系数矩阵E的逆矩阵得到 (8)求得系数矩阵A的特征值λn=
σn+iωn,(9)判断系统的稳定性特征。
[0060] 图2为某水电站引水发电系统布置简图,该电站上游引水系统采用“单管单机”,尾水系统采用“二机一室一洞”的的布置形式,尾水调压室形式采用阻抗式,该系统的状态变量为 μ1,μ2,q1,q2,q3,zsu,共8个独立变量。根据小波动稳定性分析的图论原理,给每段管道编号,再根据连续性原理,将流量变量标示到对应的每段管道上,如图2所示。以1#机组单元为例,从上游水库至下游水库列出有压管道水流动力方程。结合水流连
续方程、恒定流方程对水流动力方程进行化简并在稳定工况点处线性化,由于篇幅限制,略去化简过程,得出:
[0061]
[0062] 令: i=1~K(K为管道编号总数); i=1~N,j=1-NQ(N为机组总数,NQ为流量变量总数);Di=-2αiQ′i0,i=1~K(K为管道编号总数),下表含‘0’的为稳定工况下的已知值。上式可写为:
[0063]
[0064] 对照上式及图2可以看出,方程中状态变量的系数与水道布置存在着一一对应关系。如图2所示,1#机组对应的引水管道中,状态变量q1对应的管段为1、2;q3对应的管段为5;因此2#机组对应引水管道水流动力方程根据图2可直接写为:
[0065]
[0066] 由于调压室的存在,除了机组流量对应的状态变量q1、q2,还有另外一个关于流量的状态变量q3。同样,仍然可以利用图论原理直接列出关于这个状态变量的水流动力方程。以调压室为起点,下游水库为终点列水流动力方程。
[0067]
[0068] 式中, 为调压室水位变化相对值。 i=1~NQ。上式中含有调压室水位变量,因而还需联列其他微分方程。如图2中的调压室水位,利用调压室的连续方程得:
[0069]
[0070] 式中Asu为调压室的截面积。
[0071] 由三个水流动力方程、一个调压室连续方程、两个机组运动方程以及两个调速器方程即可构成图2所示系统的8个关于状态变量的一阶微分方程。相应的状态变量系数矩
阵构建如下:
[0072]
[0073] 只需在上式两侧同时左乘系数矩阵E的逆矩阵,则可直接得出判别系统稳定的系数矩阵A。
[0074] 选取特征计算工况(上游水位1903.3m,下游水位1820.0m,机组额定出力运行,每台机组减少2%负荷。)及两组调速器参数A(bt=0.3,Td=8,Tn=1.5)和B(bt=0.5,Td=12,Tn=1.5)。将各已知参数代入(八)中系数表达式,直接得到系统的系数矩阵E、
F、R,两侧同时左乘系数矩阵E的逆矩阵,则可得出判别系统稳定的系数矩阵A。利用QR方
法得出系统的特征值如下:
[0075] 表1系统特征值
[0076]
[0077] 由表1说明,该工况下当调速器参数选择A参数时,系统稳定;当选择B参数时,系统发散。为说明结果的正确性,利用较为成熟的水电站过渡过程计算程序对该工况进行模拟,给出调压室水位变化过程线,结果如图3~4所示。基于图论的水电站水力机械系统小
波动稳定性分析方法得出的结论与水电站过渡过程计算程序模拟结论一致,说明该方法是
正确的。
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈