CN104503228A - 一种水轮机调速器功率模式下一次调频稳定域的确定方法 - Google Patents
一种水轮机调速器功率模式下一次调频稳定域的确定方法 Download PDFInfo
- Publication number
- CN104503228A CN104503228A CN201410662711.5A CN201410662711A CN104503228A CN 104503228 A CN104503228 A CN 104503228A CN 201410662711 A CN201410662711 A CN 201410662711A CN 104503228 A CN104503228 A CN 104503228A
- Authority
- CN
- China
- Prior art keywords
- hydrogovernor
- frequency modulation
- formula
- model
- power mode
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Landscapes
- Control Of Water Turbines (AREA)
Abstract
本发明公开一种水轮机调速器功率模式下一次调频稳定域的确定方法,步骤包括:1)建立功率模式下一次调频的非线性水轮机调速器模型,非线性水轮机调速器模型包括PID控制环节以及非线性环节,根据非线性水轮机调速器模型构建水轮机调节系统模型;2)将PID控制环节的PID参数作为分岔参数,并根据水轮机调节系统模型求解发生hopf分岔时PID参数满足的分岔条件;由分岔条件得到水轮机调速器在功率模式下一次调频稳定域的边界,根据稳定域的边界确定稳定域的范围。本发明能够有效确定水轮机调速器在功率模式下一次调频的稳定域,具有实现方法简单、精度高的优点。
Description
技术领域
本发明涉及水电机组控制技术领域,尤其涉及一种水轮机调速器功率模式下一次调频稳定域的确定方法。
背景技术
电网频率是反映电力系统安全稳定运行和电能质量的重要指标,如果电网频率的变动超过了允许值,对用户和发电厂的运行都将产生不利影响,同时带来极大的经济损失。电力系统频率的稳定主要取决于电力系统内发出功率和消耗功率的平衡,目前电网频率的控制手段主要有一次调频(PFR)、二次调频(SFR)、高频切机、自动低频减负荷等,其中一次调频由于响应速度快,已成为维持电网频率稳定的重要手段。相对于火电厂和核电厂,水电厂具有运行方式灵活、调节性能优越的特点,在电力系统中承担了主要的调频任务,同时由于水电机组的一次调频是由调速器自身来完成,响应速度极快,因此水电机组的一次调频性能对于电网频率的稳定具有至关重要的影响作用。
水轮机调速器在并网状态下一般运行于开度模式或功率模式下,相比于开度模式,功率模式由于在运行时使得调速器形成功率闭环控制,能够有效避免一次调频与AGC的矛盾,并且具有负荷调节速度快、品质优的优点,因此目前并网运行的大型水电机组均优先在功率模式下运行。水轮机调速器在功率模式下运行时,由于一次调频由调速器自身来完成,通过调速器完成功率闭环控制,而整个水轮机调节系统为一非最小相位系统,如果参数整定不合理,就会影响整个系统的稳定性,导致失稳或产生周期振荡。因此为合理整定参数,保证一次调频过程的稳定,需要确定以调速器参数为影响参数的整个系统的稳定域。
水轮机调速器功率模式下一次调频系统具有明显的非线性特征,但目前水轮机调节系统稳定域的分析方法大多基于线性理论,即将非线性系统近似为线性系统进行分析,该类方法使得整个系统产生了本质变化,因而不适用于准确分析非线性的一次调频系统。另外,目前水轮机调速器在功率模式下进行一次调频的技术研究尚不成熟,对于水轮机调速器在功率模式一次调频稳定域的研究较少,目前的研究成果中也均未能给出明确的稳定域求解方法,因而提供一种有效的水轮机调速器功率模式下一次调频稳定域的确定方法是亟待解决的问题。
发明内容
本发明要解决的技术问题是克服现有技术的不足,提供一种水轮机调速器功率模式下一次调频稳定域的确定方法,具有实现方法简单且精度高的优点。
为解决上述技术问题,本发明提出的技术方案为:
一种水轮机调速器功率模式下一次调频稳定域的确定方法,步骤包括:
1)建立功率模式下一次调频的非线性水轮机调速器模型,所述非线性水轮机调速器模型包括PID控制环节以及非线性环节,根据所述非线性水轮机调速器模型构建水轮机调节系统模型;
2)将所述PID控制环节的PID参数作为分岔参数,并根据所述水轮机调节系统模型求解发生hopf分岔时所述PID参数满足的分岔条件;由所述分岔条件得到水轮机调速器在功率模式下一次调频稳定域的边界,根据所述稳定域的边界确定稳定域的范围。
作为本发明的进一步改进:所述非线性水轮机调速器模型的输入参数包括机组频率、机组功率以及给定的频率、功率,所述非线性水轮机调速器模型将输入参数依次经过人工死区非线性环节、PID控制环节、调节限幅非线性环节以及电液随动系统后输出导叶开度。
作为本发明的进一步改进:所述人工死区非线性环节包括人工频率死区以及人工功率死区,所述人工频率死区以及人工功率死区分别采用多项式函数表示;所述调节限幅非线性环节采用双曲函数表示;所述电液随动系统采用如式(1)所示的一阶线性模型表示;
式(1)中,Y1、Y分别为电液随动系统的输入、输出导叶开度,Ty为接力器反应时间常数。
作为本发明的进一步改进:所述步骤1)还包括非线性水轮机调速器模型参数核定流程,具体实施方法为:通过对水轮机调速器进行静态试验核定所述机组频率、所述PID参数、所述人工死区非线性环节的绝对值以及调节限幅非线性环节的饱和非线性系数,辨识所述电液随动系统中接力器反应时间常数Ty。
作为本发明的进一步改进:所述水轮机调节系统模型是结合水轮机模型、引水系统模型以及所述非线性水轮机调速器模型构建得到。
作为本发明的进一步改进:所述水轮机模型采用如式(2)所示的线性模型;所述引水系统模型采用如式(3)所示的刚性水击模型;
q=eqxx+eqyy+eqhh
(2)
mt=exx+eyy+ehh
式(2)中,q为水轮机流量的标幺值,mt为水轮机动力力矩的标幺值,x为机组转速的标幺值,y为导叶开度的标幺值,h为水轮机工作水头的标幺值;eqx、eqy、eqh和ex、ey和eh分别为传递系数;
式(3)中,Tw为水流加速时间常数。
作为本发明的进一步改进,所述步骤1)还包括水轮机模型与引水系统模型参数计算流程,具体实施方法为:通过水轮机的综合特性曲线获取所述水轮机模型中传递系数;通过模型辨识方法辨识出所述引水系统模型中水流加速时间常数Tw。
作为本发明的进一步改进,所述步骤2)中求解发生hopf分岔时所述PID参数满足的分岔条件的具体实施步骤为:
2.1)根据所述水轮机调节系统模型构建如式(4)所示的状态方程;
式(4)中,x为状态变量,为状态响应,μ为分岔参数且即为PID参数,Rn为n维的实数集;
2.2)构建如式(5)所示的平衡点方程,得到状态方程在平衡点的Jacobian矩阵;由所述Jacobian矩阵的特征方程的系数构造如式(6)所示的Hurwitz行列式;
A(μ)=Dx(0,μ) (5)
式(5)中,Dx(0,μ)表示所述状态方程在平衡点f(x,μ)=0处含分岔参数μ的方程组的系数矩阵,A(μ)表示所述状态方程在平衡点的Jacobian矩阵;
式(6)中,a1~am为Jacobian矩阵A(μ)的特征方程的系数,且m=1,2,…,n,n为维数,Δm为m维的Hurwitz行列式,其中如果i>n则ai=0;
2.3)由所述Hurwitz行列式得到发生hopf分岔时PID参数满足的分岔条件,所述分岔条件的表达式如式(7)所示;
式(7)中,μc表示发生hopf分岔时的分岔参数。
作为本发明的进一步改进:所述步骤2)中采用PID参数中缓冲时间常数Td和暂态转差系数bt作为分岔参数。
与现有技术相比,本发明的优点在于:
1)本发明通过建立包含非线性环节的非线性水轮机调速器模型,引入hopf分岔理论分析非线性的水轮机调节系统的稳定性,由发生hopf分岔时PID参数满足的分岔条件得到水轮机调速器在功率模式下一次调频稳定域的边界,从而确定稳定域的范围,有效解决了一次调频稳定域的确定问题,实现方法简单且精度高。
2)本发明考虑人工死区以及调节限幅的非线性特征建立非线性水轮机调速器模型,能够准确表述功率模式下的一次调频系统的非线性特性,从而获得准确的一次调频稳定域范围。
3)本发明采用连续的数学函数表示非线性环节,人工死区非线性环节采用多项式函数近似逼近、调节限幅非线性环节采用表示双曲函数近似逼近,有效的表述非线性环节的特征同时便于对非线性系统进行频域分析。
附图说明
图1是本实施例水轮机调速器功率模式下一次调频稳定域的确定方法实现流程示意图。
图2是本实施例中水轮机调速器在功率模式下一次调频的控制原理示意图。
图3是本实施例中水轮机调速器模型的结构示意图。
图4是本实施例中水轮机模型的结构示意图。
图5是本实施例中引水系统模型的结构示意图。
图6是本发明具体实施例中水轮机调速器功率模式下一次调频稳定域的结果示意图。
具体实施方式
以下结合说明书附图和具体优选的实施例对本发明作进一步描述,但并不因此而限制本发明的保护范围。
本发明水轮机调速器功率模式下一次调频稳定域的确定方法,其核心思想是通过考虑非线性特征建立功率模式下一次调频的非线性水轮机调速器模型,再引入hopf分岔理论对非线性的水轮机调节系统进行分岔分析,进而求解得到水轮机调速器功率模式下一次调频稳定域的范围。
如图1所示,本实施水轮机调速器功率模式下一次调频稳定域的确定方法,步骤包括:
1)建模:建立功率模式下一次调频的非线性水轮机调速器模型,非线性水轮机调速器模型包括PID控制环节以及非线性环节,根据非线性水轮机调速器模型构建水轮机调节系统模型;
2)hopf分岔分析:将PID控制环节的PID参数作为分岔参数,并根据水轮机调节系统模型求解发生hopf分岔时PID参数满足的分岔条件;由分岔条件得到水轮机调速器在功率模式下一次调频稳定域的边界,根据稳定域的边界确定稳定域的范围。
由于水轮机调速器功率模式下一次调频系统具有明显的非线性特征,而混沌理论揭示了非线性系统不同运动状态之间的联系和转化,因而混沌理论能够用于一次调频非线性系统的稳定性分析。hopf分岔理论是混沌分岔理论中的一种,也是其中一种最基本的分岔问题,能够分析系统中稳定或不稳定的振荡及其对应的相空间稳定或不稳定极限环。本实施例通过引入hopf分岔理论能够有效分析得到水轮机调速器在功率模式下一次调频系统的稳定性,解决目前未能确定一次调频稳定域的问题,以合理的整定一次调频参数,避免水轮机调节系统产生周期振荡,从而确保电网的安全稳定运行。
本实施例中水轮机调节系统控制原理如图2所示,在功率模式下一次调频过程中,水轮机调速器的输入为机组频率及水轮机有功功率,调速器根据机组频率和有功功率反馈调节水轮机的导叶开度,水轮机导叶开度变化引起了流道内流量、水头的变化及水轮机动力矩(有功功率)的变化,最终机组将稳定在新的平衡状态。
由于一次调频人工死区和调节限幅的存在,水轮机调速器具有非线性特性,本实施例中考虑一次调频人工死区、调节限幅的非线性特征建立功率模型下一次调频的非线性水轮机调节器模型,能够更为准确的表述功率模式下的一次调频系统的非线性特性,从而获得准确的一次调频稳定域。如图3所示,非线性水轮机调速器模型的输入参数包括机组频率、机组功率以及给定的频率、功率,非线性水轮机调速器模型将输入参数依次经过人工死区非线性环节、PID控制环节、调节限幅非线性环节以及电液随动系统后输出导叶开度,其中人工死区非线性环节包括人工频率死区以及人工功率死区。非线性水轮机调速器模型中各环节可采用如下表达式表示,其中人工频率死区表达式为:
Δf1=fc-fg (8)
式(8)、(9)中,Δf1、Δf2分别为人工频率死区的输入、输出频率,fc为频率给定值,fg为机组频率,Ef为人工频率死区的绝对值。
人工功率死区表达式为:
Δp1=pc-pg (10)
式(10)、(11)中,Δp1、Δp2分别为人工功率死区的输入、输出功率pc为功率给定值,pg为机组功率;Ep为人工功率死区的绝对值。
PID控制环节的表达式为:
YIN=Δf2+Δp2ep (12)
式(12)、(13)中,YIN、YPID分别为PID控制环节的输入、输出值,ep为功率模式永态转差系数,Kp、KI、KD分别为比例、积分、微分增益,其中Td为缓冲时间常数,bt为暂态转差系数,Tn为微分时间常数。
调节限幅非线性环节的表达式为:
式(14)中,Y1表示调节限幅非线性环节的输出值,Ym,Y0为饱和非线性系数。
本实施例中,电液随动系统采用如(1)所示的一阶模型表示:
式(1)中,Y1、Y分别为电液随动系统的输入、输出导叶开度,Ty为接力器反应时间常数。
由于式(9)、(11)、(14)所表示的人工死区线性环节与调节限幅线性环节是采用非连续的分段函数表示,为进行频域分析,则需要将人工死区非线性环节和调节限幅值非线性环节采用连续函数进行近似描述。本实施例中,人工频率死区以及人工功率死区分别采用多项式函数表示,调节限幅非线性环节采用双曲函数表示,其中人工频率死区以及人工功率死区分别采用如式(15)、式(16)所示的多项式进行近似逼近,调节限幅非线性环节采用如式(17)所示的双曲函数进行近似逼近;
Δf2=anΔf1 n+an-1Δf1 n-1+......+a1Δf1+a0 (15)
Δp2=bnΔp1 n+bn-1Δp1 n-1+......+b1Δp1+b0 (16)
式(15)中a0~an为人工频率死区逼近系数;式(16)中b0~bn为人工功率死区逼近系数。
式(17)中c、d为调节限幅非线性环节逼近系数。
本实施例中,人工频率死区逼近系数a0~an、人工功率死区逼近系数b0~bn和双曲函数的系数(c、d)分别通过多组静态实验并采用最小二乘法预先计算得到,具体实施步骤为:
①对水轮机调速器进行静态试验,检测机组频率变化值在[-0.2Hz,0.2Hz]范围内的一组数据,即人工频率死区的输入、输出频率Δf1,i,Δf2,i(i=0,1,2,...);检测功率变化相对值在[-10%,10%]范围内的一组数据,即人工功率死区的输入、输出功率Δp1,i,Δp2,i(i=0,1,2,...);检测导叶开度变化相对值在[-10%,10%]范围内的一组数据,即电液随动系统输入、输出值Y1,i,YPID,i(i=0,1,2,...);
②根据测得的Δf1,i,Δf2,i及Δp1,i,Δp2,i采用三次多项式(n=3)表示并以最小二乘法进行拟合,计算得到人工频率死区逼近系数a0,a1,a2,a3及人工功率死区逼近系数b0,b1,b2,b3;根据测得的Y1,i,YPID,i,采用双曲函数表示并以最小二乘法进行拟合,计算得到非线性环节逼近系数c、d。
本实施例中,步骤1)还包括非线性水轮机调速器模型参数核定流程,具体实施方法为:通过对水轮机调速器进行静态试验核定机组频率、PID参数、人工死区非线性环节中人工频率死区的绝对值Ef,人工功率死区的绝对值Ep以及调节限幅非线性环节中的饱和非线性系数Ym、Y0;通过试验辨识电液随动系统中接力器反应时间常数Ty。通过核定参数保证与实际水轮机调速器的一致性,具体实现方法为在仿真模型中采用实际调速器设定的参数进行仿真,根据仿真和实测曲线对参数进行校核。
本实施例中,水轮机调节系统模型结合水轮机模型、引水系统型、发电机模型以及非线性水轮机调速器模型构建得到,其中水轮机模型采用线性模型,引水系统模型采用刚性水击模型。
由于一次调频过渡过程为小扰动过程,因此水轮机模型可采用六个传递系数表示的简化模型表示,如图4所示,本实施例中水轮机模型采用如式(2)所示的线性模型:
q=eqxx+eqyy+eqhh
(2)
mt=exx+eyy+ehh
式(2)中,q为水轮机流量的标幺值,mt为水轮机动力力矩的标幺值,x为机组转速的标幺值,y为导叶开度的标幺值,h为水轮机工作水头的标幺值。
如图5所示,本实施例中引水系统模型采用如式(3)所示的刚性水击模型:
式(3)中,Tw为水流加速时间常数。
本实施例中,步骤1)还包括水轮机模型与引水系统模型参数计算流程,具体实施方法为:通过水轮机的综合特性曲线获取所述水轮机模型中传递系数(eqx,eqy,eqh,ex,ey,eh);通过模型辨识方法辨识出引水系统模型中水流加速时间常数Tw。在试验中退出人工频率/功率死区及导叶开度调节限幅非线性环节,由所得到的试验结果即可辨识出引水系统模型中的水流加速时间常数Tw的值。
本实施例中,建立水轮机调节系统模型后,引入hopf分岔理论对非线性的水轮机调节系统进行稳定性分析,通过hopf分岔的代数判据来判断系统hopf分岔的存在性及分岔值,从而求出发生hopf分岔时作为分岔参数的PID参数应满足的关系,得到以PID参数为影响因素的水轮机调速器功率模式下一次调频稳定域。
本实施例中,步骤2)中求解发生hopf分岔时PID参数满足的分岔条件的具体实施步骤为:
2.1)根据水轮机调节系统模型构建如式(4)所示的状态方程;
式(4)中,x为状态变量,为状态响应,μ为分岔参数且为PID参数,Rn为n维的实数集;
2.2)构建如式(5)所示的平衡点方程,得到状态方程在平衡点的Jacobian矩阵;由Jacobian矩阵的特征方程的系数构造如式(6)所示的Hurwitz行列式;
A(μ)=Dx(0,μ) (5)
式(5)中,Dx(0,μ)表示状态方程在平衡点f(x,μ)=0处含分岔参数μ的方程组的系数矩阵,A(μ)表示状态方程在平衡点的Jacobian矩阵;
式(6)中,a1~am为Jacobian矩阵A(μ)的特征方程的系数,且m=1,2,…,n,n为维数,Δm为m维的Hurwitz行列式,其中如果i>n则ai=0;
2.3)由Hurwitz行列式得到发生hopf分岔时PID参数满足的分岔条件,分岔条件的表达式如式(7)所示;
式(7)中,μc表示发生hopf分岔时的分岔参数。
由于Jacobian矩阵A(μ)的特征方程展开为λn+a1(μ)λn-1+…+an-1(μ)λ+an(μ)=0,
其中λ为Jacobian矩阵A(μ)的特征值,a1(μ)…an(μ)为含分岔参数μ的系数,n为维数,由a1(μ)…an(μ)即可构造得到如(6)所示的Hurwitz行列式。
本实施例中,步骤2)中采用PID参数中缓冲时间常数Td和暂态转差系数bt作为分岔参数,当水轮机调节系统的状态方程发生分岔时暂态转差系数bt和缓冲时间常数Td所满足的条件即为水轮机调速器在功率模式下一次调频稳定域的边界,即满足Δn-3(bt、Td)>0、Δn-1(bt、Td)=0的取值点。由稳定域的边界即可确定得到暂态转差系数bt和缓冲时间常数Td的稳定域,即暂态转差系数bt和缓冲时间常数Td的取值范围,在稳定域内整定暂态转差系数bt和缓冲时间常数Td就能够保证系统的稳定。
本实施例中,由得到的分岔条件具体采用两点法迭代求解稳定域的边界,以下以求解式Δn-3(bt、Td)>0为例说明求解的具体步骤:
①确定参数一(Td、bt其中任意一个参数)并设为x1;
②确定参数二(设为y)的变化范围在[y1,y2]之间,且保证Δn-3(x1、y1)Δn-3(x1、y2)<0;
③令y10=y1,y20=y2,对于k=0,1,....,M循环执行以下a)~c)步骤:
a)计算若Δn-3(x1、yk)<ε,所得到的yk则为满足式(1)的解,停止计算并退出,否则转至执行b),其中ε为预先给定的精度水平;
b)判断是否满足Δn-3(x1、y1,k)Δn-3(x1、yk)<0,如果是则令y1,k+1=y1,k,y2,k+1=yk,如果Δn-3(x1、y1,k)Δn-3(x1、yk)>0,则令y1,k+1=yk,y2,k+1=y2,k,转入执行c);
c)判断是否满足k=M,如果是则输出M次迭代不成功的信息,否则返回执行步骤a)。
④根据预先设定的计算步长重新确定参数,返回执行步骤①,直到计算得到x=x1,.....xn时满足Δn-3(bt、Td)>0的每一个x对应的yk,即得到满足Δn-3(bt、Td)>0的bt和Td。
本实施例具体可通过计算机编程按上述方法求解得到所有满足Δn-3(bt、Td)>0、Δn-1(bt、Td)=0的bt、Td,即可得到由bt和Td的取值点构成的稳定域边界。
如图6所示为本发明应用于具体实施例中得到的稳定域,附图中曲线表示稳定域边界,即为满足分岔条件的bt、Td的取值点,在稳定域边界内bt、Td的取值范围即为稳定域范围。
上述只是本发明的较佳实施例,并非对本发明作任何形式上的限制。虽然本发明已以较佳实施例揭露如上,然而并非用以限定本发明。任何熟悉本领域的技术人员,在不脱离本发明技术方案范围的情况下,都可利用上述揭示的技术内容对本发明技术方案做出许多可能的变动和修饰,或修改为等同变化的等效实施例。因此,凡是未脱离本发明技术方案的内容,依据本发明技术实质对以上实施例所做的任何简单修改、等同变化及修饰,均应落在本发明技术方案保护的范围内。
Claims (9)
1.一种水轮机调速器功率模式下一次调频稳定域的确定方法,其特征在于,步骤包括:
1)建立功率模式下一次调频的非线性水轮机调速器模型,所述非线性水轮机调速器模型包括PID控制环节以及非线性环节,根据所述非线性水轮机调速器模型构建水轮机调节系统模型;
2)将所述PID控制环节的PID参数作为分岔参数,并根据所述水轮机调节系统模型求解发生hopf分岔时所述PID参数满足的分岔条件;由所述分岔条件得到水轮机调速器在功率模式下一次调频稳定域的边界,根据所述稳定域的边界确定稳定域的范围。
2.根据权利要求1所述的水轮机调速器功率模式下一次调频稳定域的确定方法,其特征在于:所述非线性水轮机调速器模型的输入参数包括机组频率、机组功率以及给定的频率、功率,所述非线性水轮机调速器模型将输入参数依次经过人工死区非线性环节、PID控制环节、调节限幅非线性环节以及电液随动系统后输出导叶开度。
3.根据权利要求2中所述的水轮机调速器功率模式下一次调频稳定域的确定方法,其特征在于:所述人工死区非线性环节包括人工频率死区以及人工功率死区,所述人工频率死区以及人工功率死区分别采用多项式函数表示;所述调节限幅非线性环节采用双曲函数表示;所述电液随动系统采用如式(1)所示的一阶线性模型表示;
式(1)中,Y1、Y分别为电液随动系统的输入、输出导叶开度,Ty为接力器反应时间常数。
4.根据权利要求3所述的水轮机调速器功率模式下一次调频稳定域的确定方法,其特征在于:所述步骤1)还包括非线性水轮机调速器模型参数核定流程,具体实施方法为:通过对水轮机调速器进行静态试验核定所述机组频率、所述PID参数、所述人工死区非线性环节的绝对值以及调节限幅非线性环节的饱和非线性系数,辨识所述电液随动系统中接力器反应时间常数Ty。
5.根据权利要求4所述的水轮机调速器功率模式下一次调频稳定域的确定方法,其特征在于:所述水轮机调节系统模型是结合水轮机模型、引水系统模型以及所述非线性水轮机调速器模型构建得到。
6.根据权利要求5所述的水轮机调速器功率模式下一次调频稳定域的确定方法,其特征在于:所述水轮机模型采用如式(2)所示的线性模型;所述引水系统模型采用如式(3)所示的刚性水击模型;
q=eqxx+eqyy+eqhh
(2)
mt=exx+eyy+ehh
式(2)中,q为水轮机流量的标幺值,mt为水轮机动力力矩的标幺值,x为机组转速的标幺值,y为导叶开度的标幺值,h为水轮机工作水头的标幺值;eqx、eqy、eqh和ex、ey和eh分别为传递系数;
式(3)中,Tw为水流加速时间常数。
7.根据权利要求6所述的水轮机调速器功率模式下一次调频稳定域的确定方法,其特征在于,所述步骤1)还包括水轮机模型与引水系统模型参数计算流程,具体实施方法为:通过水轮机的综合特性曲线获取所述水轮机模型中传递系数;通过模型辨识方法辨识出所述引水系统模型中水流加速时间常数Tw。
8.根据权利要求1~7中任意一项所述的水轮机调速器功率模式下一次调频稳定域的确定方法,其特征在于,所述步骤2)中求解发生hopf分岔时所述PID参数满足的分岔条件的具体实施步骤为:
2.1)根据所述水轮机调节系统模型构建如式(4)所示的状态方程;
式(4)中,x为状态变量,为状态响应,μ为分岔参数且即为PID参数,Rn为n维的实数集;
2.2)构建如式(5)所示的平衡点方程,得到状态方程在平衡点的Jacobian矩阵;由所述Jacobian矩阵的特征方程的系数构造如式(6)所示的Hurwitz行列式;
A(μ)=Dx(0,μ) (5)
式(5)中,Dx(0,μ)表示所述状态方程在平衡点f(x,μ)=0处含分岔参数μ的方程组的系数矩阵,A(μ)表示所述状态方程在平衡点的Jacobian矩阵;
式(6)中,a1~am为Jacobian矩阵A(μ)的特征方程的系数,且m=1,2,…,n,n为维数,△m为m维的Hurwitz行列式,其中如果i>n则ai=0;
2.3)由所述Hurwitz行列式得到发生hopf分岔时PID参数满足的分岔条件,所述分岔条件的表达式如式(7)所示;
式(7)中,μc表示发生hopf分岔时的分岔参数。
9.根据权利要求8所述的水轮机调速器功率模式下一次调频稳定域的确定方法,其特征在于:所述步骤2)中采用PID参数中缓冲时间常数Td和暂态转差系数bt作为分岔参数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410662711.5A CN104503228B (zh) | 2014-11-19 | 2014-11-19 | 一种水轮机调速器功率模式下一次调频稳定域的确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410662711.5A CN104503228B (zh) | 2014-11-19 | 2014-11-19 | 一种水轮机调速器功率模式下一次调频稳定域的确定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104503228A true CN104503228A (zh) | 2015-04-08 |
CN104503228B CN104503228B (zh) | 2017-02-08 |
Family
ID=52944632
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410662711.5A Active CN104503228B (zh) | 2014-11-19 | 2014-11-19 | 一种水轮机调速器功率模式下一次调频稳定域的确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104503228B (zh) |
Cited By (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104806438A (zh) * | 2015-05-04 | 2015-07-29 | 河海大学 | 非最小相位控制的水轮机组转速调节系统及设计方法 |
CN105065191A (zh) * | 2015-07-16 | 2015-11-18 | 国家电网公司 | 一种加快高水头水电机组减负荷后系统稳定的方法 |
CN105634009A (zh) * | 2016-01-22 | 2016-06-01 | 重庆大唐国际武隆水电开发有限公司 | 一种水电厂调频调节方法 |
CN105717793A (zh) * | 2016-02-03 | 2016-06-29 | 国网山东省电力公司电力科学研究院 | 一种确保一次调频动作方向的优化控制方法及系统 |
CN105868497A (zh) * | 2016-04-19 | 2016-08-17 | 国家电网公司 | 一种带变顶高尾水隧洞水轮机调节系统仿真建模方法及模型 |
CN106842913A (zh) * | 2016-12-02 | 2017-06-13 | 上海电机学院 | 一种基于随机概率分布控制的水轮机调节系统 |
CN107102554A (zh) * | 2017-06-22 | 2017-08-29 | 北京石油化工学院 | 一种磁悬浮球形飞轮不平衡振动抑制方法 |
CN107800146A (zh) * | 2017-11-16 | 2018-03-13 | 国网四川省电力公司电力科学研究院 | 兼顾一次调频和超低频振荡抑制的调速器参数优化方法 |
CN107942664A (zh) * | 2017-11-23 | 2018-04-20 | 中国南方电网有限责任公司 | 一种基于灵敏度分析的水轮机调速器参数整定方法及系统 |
CN108512233A (zh) * | 2018-05-17 | 2018-09-07 | 贵州电网有限责任公司 | 一种水电机组基于实际水头一次调频计算方法 |
CN110492483A (zh) * | 2019-09-17 | 2019-11-22 | 华北电力大学 | 一种一次调频反馈通道非线性环节参数配置方法和系统 |
CN110970911A (zh) * | 2019-12-13 | 2020-04-07 | 四川省电力工业调整试验所 | 一种开度模式下agc和一次调频相互叠加的控制方法 |
CN112003317A (zh) * | 2020-06-05 | 2020-11-27 | 国电南瑞科技股份有限公司 | 一种水电机组功率调节输入输出优化方法及装置 |
CN112270088A (zh) * | 2020-10-26 | 2021-01-26 | 中国长江电力股份有限公司 | 集群水轮发电机组仿真模型及其建模方法 |
CN112436554A (zh) * | 2020-11-29 | 2021-03-02 | 国家电网有限公司 | 一种水力发电并网仿真方法及系统 |
CN113224772A (zh) * | 2021-04-20 | 2021-08-06 | 华北电力大学 | 一种考虑切换环节的水电系统超低频频率振荡抑制方法 |
CN114204609A (zh) * | 2021-11-05 | 2022-03-18 | 华能澜沧江水电股份有限公司 | 一种基于有功功率调节死区的一次调频转嫁调节方法 |
CN114329835A (zh) * | 2021-12-30 | 2022-04-12 | 华中科技大学 | 一种多机共管道水轮机调节系统稳定域确定方法及系统 |
CN114352469A (zh) * | 2021-12-08 | 2022-04-15 | 国网江西省电力有限公司电力科学研究院 | 一种高阶水轮机调节系统稳定域计算方法及系统 |
CN115085287A (zh) * | 2022-07-20 | 2022-09-20 | 华北电力大学 | 一种考虑调速器限幅环节频率稳定约束的建模方法 |
CN116047915A (zh) * | 2023-02-21 | 2023-05-02 | 陕西铁路工程职业技术学院 | 水轮机全负荷工况自适应控制方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS6463307A (en) * | 1987-09-03 | 1989-03-09 | Kubota Ltd | Coupling device for tractor and working machine |
-
2014
- 2014-11-19 CN CN201410662711.5A patent/CN104503228B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS6463307A (en) * | 1987-09-03 | 1989-03-09 | Kubota Ltd | Coupling device for tractor and working machine |
Non-Patent Citations (5)
Title |
---|
丁聪: "混流式水轮机调节系统建模及非线性动力分析", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
叶波 等: "水轮机调速器在功率闭环调节模式下的一次调频试验", 《水电站机电技术》 * |
田伦: "水轮机调速器功率模式在索风营电厂的运行实践", 《华电技术》 * |
袁璞 等: "考虑限幅环节的水轮机调节系统建模与分析", 《人民长江》 * |
鲍海艳 等: "基于水电站运行稳定性的调压室设置条件探讨", 《水力发电学报》 * |
Cited By (34)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104806438A (zh) * | 2015-05-04 | 2015-07-29 | 河海大学 | 非最小相位控制的水轮机组转速调节系统及设计方法 |
CN105065191B (zh) * | 2015-07-16 | 2017-07-28 | 国家电网公司 | 一种加快高水头水电机组减负荷后系统稳定的方法 |
CN105065191A (zh) * | 2015-07-16 | 2015-11-18 | 国家电网公司 | 一种加快高水头水电机组减负荷后系统稳定的方法 |
CN105634009A (zh) * | 2016-01-22 | 2016-06-01 | 重庆大唐国际武隆水电开发有限公司 | 一种水电厂调频调节方法 |
CN105717793A (zh) * | 2016-02-03 | 2016-06-29 | 国网山东省电力公司电力科学研究院 | 一种确保一次调频动作方向的优化控制方法及系统 |
CN105717793B (zh) * | 2016-02-03 | 2018-04-06 | 国网山东省电力公司电力科学研究院 | 一种确保一次调频动作方向的优化控制方法及系统 |
CN105868497B (zh) * | 2016-04-19 | 2018-11-23 | 国家电网公司 | 一种带变顶高尾水隧洞水轮机调节系统仿真建模方法及模型 |
CN105868497A (zh) * | 2016-04-19 | 2016-08-17 | 国家电网公司 | 一种带变顶高尾水隧洞水轮机调节系统仿真建模方法及模型 |
CN106842913A (zh) * | 2016-12-02 | 2017-06-13 | 上海电机学院 | 一种基于随机概率分布控制的水轮机调节系统 |
CN106842913B (zh) * | 2016-12-02 | 2020-01-17 | 上海电机学院 | 一种基于随机概率分布控制的水轮机调节系统 |
CN107102554A (zh) * | 2017-06-22 | 2017-08-29 | 北京石油化工学院 | 一种磁悬浮球形飞轮不平衡振动抑制方法 |
CN107800146A (zh) * | 2017-11-16 | 2018-03-13 | 国网四川省电力公司电力科学研究院 | 兼顾一次调频和超低频振荡抑制的调速器参数优化方法 |
CN107800146B (zh) * | 2017-11-16 | 2020-02-21 | 国网四川省电力公司电力科学研究院 | 兼顾一次调频和超低频振荡抑制的调速器参数优化方法 |
CN107942664A (zh) * | 2017-11-23 | 2018-04-20 | 中国南方电网有限责任公司 | 一种基于灵敏度分析的水轮机调速器参数整定方法及系统 |
CN107942664B (zh) * | 2017-11-23 | 2020-11-06 | 中国南方电网有限责任公司 | 一种基于灵敏度分析的水轮机调速器参数整定方法及系统 |
CN108512233A (zh) * | 2018-05-17 | 2018-09-07 | 贵州电网有限责任公司 | 一种水电机组基于实际水头一次调频计算方法 |
CN108512233B (zh) * | 2018-05-17 | 2021-02-26 | 贵州电网有限责任公司 | 一种水电机组基于实际水头一次调频计算方法 |
CN110492483A (zh) * | 2019-09-17 | 2019-11-22 | 华北电力大学 | 一种一次调频反馈通道非线性环节参数配置方法和系统 |
CN110970911A (zh) * | 2019-12-13 | 2020-04-07 | 四川省电力工业调整试验所 | 一种开度模式下agc和一次调频相互叠加的控制方法 |
CN112003317A (zh) * | 2020-06-05 | 2020-11-27 | 国电南瑞科技股份有限公司 | 一种水电机组功率调节输入输出优化方法及装置 |
CN112270088B (zh) * | 2020-10-26 | 2023-07-18 | 中国长江电力股份有限公司 | 集群水轮发电机组仿真模型及其建模方法 |
CN112270088A (zh) * | 2020-10-26 | 2021-01-26 | 中国长江电力股份有限公司 | 集群水轮发电机组仿真模型及其建模方法 |
CN112436554A (zh) * | 2020-11-29 | 2021-03-02 | 国家电网有限公司 | 一种水力发电并网仿真方法及系统 |
CN112436554B (zh) * | 2020-11-29 | 2024-01-19 | 国家电网有限公司 | 一种水力发电并网仿真方法及系统 |
CN113224772A (zh) * | 2021-04-20 | 2021-08-06 | 华北电力大学 | 一种考虑切换环节的水电系统超低频频率振荡抑制方法 |
CN114204609A (zh) * | 2021-11-05 | 2022-03-18 | 华能澜沧江水电股份有限公司 | 一种基于有功功率调节死区的一次调频转嫁调节方法 |
CN114204609B (zh) * | 2021-11-05 | 2024-01-19 | 华能澜沧江水电股份有限公司 | 一种基于有功功率调节死区的一次调频转嫁调节方法 |
CN114352469A (zh) * | 2021-12-08 | 2022-04-15 | 国网江西省电力有限公司电力科学研究院 | 一种高阶水轮机调节系统稳定域计算方法及系统 |
CN114352469B (zh) * | 2021-12-08 | 2024-02-09 | 国网江西省电力有限公司电力科学研究院 | 一种高阶水轮机调节系统稳定域计算方法及系统 |
CN114329835A (zh) * | 2021-12-30 | 2022-04-12 | 华中科技大学 | 一种多机共管道水轮机调节系统稳定域确定方法及系统 |
CN114329835B (zh) * | 2021-12-30 | 2024-05-28 | 华中科技大学 | 一种多机共管道水轮机调节系统稳定域确定方法及系统 |
CN115085287A (zh) * | 2022-07-20 | 2022-09-20 | 华北电力大学 | 一种考虑调速器限幅环节频率稳定约束的建模方法 |
CN116047915A (zh) * | 2023-02-21 | 2023-05-02 | 陕西铁路工程职业技术学院 | 水轮机全负荷工况自适应控制方法 |
CN116047915B (zh) * | 2023-02-21 | 2023-08-15 | 陕西铁路工程职业技术学院 | 水轮机全负荷工况自适应控制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104503228B (zh) | 2017-02-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104503228B (zh) | 一种水轮机调速器功率模式下一次调频稳定域的确定方法 | |
CN103268728B (zh) | 一种电力系统动态仿真压水反应堆控制系统模型构建方法 | |
CN101446807B (zh) | 一种电力系统仿真中火电厂调速系统模型的实现方法 | |
CN109190145B (zh) | 一种包含功率控制系统的核电机组仿真模型 | |
CN105068424B (zh) | 一种适用于电力系统分析的转桨式水轮机调节系统动态模型 | |
CN105005654B (zh) | 一种包含引水系统的调速系统数模混合实时仿真方法 | |
Zhang et al. | A novel surface-cluster approach towards transient modeling of hydro-turbine governing systems in the start-up process | |
CN105863948B (zh) | 一种带变顶高尾水隧洞水轮机调速器变参数控制方法 | |
CN106532744A (zh) | 一种基于自抗扰控制的光伏电站抑制电网低频振荡的方法 | |
CN104389733A (zh) | 基于不确定性模型的水轮机pid调速器控制参数整定方法 | |
Gil-González et al. | Passivity-based control and stability analysis for hydro-turbine governing systems | |
Zhang et al. | Detailed modelling and parameters optimisation analysis on governing system of hydro‐turbine generator unit | |
CN106777944A (zh) | 一种水电机组经直流送出系统的调速器参数整定方法 | |
Beus et al. | Application of an adaptive model predictive control algorithm on the Pelton turbine governor control | |
CN107065556A (zh) | 一种堆芯机组变功率运行策略优化方案的自动搜索方法 | |
CN104500150B (zh) | 一种汽轮机调速器执行机构分段线性模型及参数获取方法 | |
CN105449698A (zh) | 一种新型水电机组负荷与频率控制器 | |
CN106327057A (zh) | 一种水电机组一次调频理论动作电量的计算系统及方法 | |
CN101718248A (zh) | 水轮发电机组主动稳定控制方法 | |
Beus et al. | Application of model predictive control algorithm on a hydro turbine governor control | |
Raach et al. | ℋ∞ controller design for closed-loop wake redirection | |
Beus et al. | Practical implementation of a hydro power unit active power regulation based on an MPC algorithm | |
CN105065191A (zh) | 一种加快高水头水电机组减负荷后系统稳定的方法 | |
CN105022286B (zh) | 一种水轮发电机组主接力器反应时间常数测试方法 | |
CN106528978A (zh) | 一种先进压水堆核电机组调速系统的动态仿真方法和系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |