Pg电子平台：基于能量函数设计抑制器的方法与流程
技术领域本发明属于电力系统控制领域，具体涉及一种基于能量函数设计抑制器的方法。
现代电力系统是一个典型的高维数、非线性复杂动态系统，其正常运行时不断遭受各种各样的扰动(如故障、切机、切负荷等)，机电扰动在系统中的传播可能会影响到电力系统的安全稳定运行。因此，如何有效地保证电力系统的安全稳定运行一直受到人们的广泛关注。电力系统的各种控制装置是保证系统安全稳定运行的重要手段，现在电力系统中抑制器的设计方法主要是将非线性系统在运行点附近进行线性化。但电力系统的强非线性，使得局部线性化的方法设计的抑制器不能适应电力系统受到大扰动后对暂态性能的要求。近年来，非线性控制已经广泛用于电力系统中，并取得了丰富的成果，能量函数是非线性理论的重要组成部分，基于能量函数的Lyapunov理论曾受到了广泛关注并用于抑制器的设计之中。但是，对于如何构造Lyapunov函数一直没有系统理论，这也限制了Lyapunov函数的推广与应用。Hamilton系统广泛存在于数理科学、生命科学和工程科学的各个领域，是当今非线性科学研究的重要方向。近年来，人们已经将Hamilton能量控制理论引入电力系统的非线性控制领域。但目前许多文献在构造Hamilton能量函数时多数采用经验法完成，并没有深入研究其构造方法，且现有的扰动抑制器系数确定采用经验法完成，无法提供详细的理论依据。
针对现有技术中的上述不足，本发明提供了一种基于能量函数设计抑制器的方法，采用该方法设计的抑制器能够产生附加阻尼，有效抑制机电扰动，保证系统的安全稳定运行。为了达到上述发明目的，本发明采用的技术方案为：提供一种基于能量函数设计抑制器的方法，其包括以下步骤：建立带励磁和汽门开度控制的电力系统状态空间模型：x=f(x)+g1(x)u+g2(x)w]]
其中，x＝(x1，x2，…，xn)T∈Rn，其为n维状态变量；u＝(u1，u2，…，um)T∈Rm，其为m维控制变量；w∈Rr，其为r维扰动变量；f(x)∈Rn，其为n维状态变量；g1(x)∈Rn×Rm，g2(x)∈Rn×Rr，g1(x)和g2(x)均为微分流形；对建立的电力系统状态空间模型进行泰勒展开，并基于常值实现方法建立能量函数：H(x)=i=1nxi(0)xiri(x1,x2,...,xn)dxi]]
其中，i＝1…n，f(x)∈Rn，f(x)为n维状态变量，(·)T表示转置，x＝(x1，x2，…，xn)T∈Rn，其为n维状态变量；(x1(0)，x2(0)，…，xn(0))T为任一初始点；基于能量函数表述电力系统状态空间模型：x=J(x)-R(x)H+g1(x)u+g2(x)w]]
其中，J(x)为n×n的反对称阵；R(x)为n×n的半正定矩阵；u＝(u1，u2，…，um)T∈Rm，其为m维控制变量；w∈Rr，其为r维扰动变量；f(x)∈Rn，其为n维状态变量；g1(x)∈Rn×Rm，g2(x)∈Rn×Rr，g1(x)和g2(x)均为微分流形；对采用能量函数表述后的电力系统状态空间模型设计抑制器：u=-1/2hT(x)h(x)+1/(22)Img1T(x)H]]
其中，h(x)为抑制器系数；(·)T表示转置；g1(x)∈Rn×Rm，其为微分流形；γ为设定的扰动抑制系数。本发明的有益效果为：采用本方法设计抑制器时，考虑了汽门开度的影响和常值实现法，避开了Hamilton能量函数的一般充分条件，扩大了常值实现的应用范围。在确定抑制器系数时，通过特征根进行求取，为励磁和汽门相协调的扰动抑制器系数确定提供了明确方法，使得抑制器系数确定不再仅仅依靠经验法。采用本方法设计的抑制器在电力系统发生扰动后，抑制器能够产生附加阻尼，有效抑制机电扰动，保证系统的安全稳定运行；扰动消失后，励磁和汽门相协调的扰动抑制器使得系统恢复到初始平衡点。另外，将本方法设计的抑制器用于控制单机无穷大系统时，其优于单独的励磁控制，加入协调控制后，系统在20s内恢复到初始平衡点，比单独励磁控制抑制机电扰动所需时间减少3s，暂态变量Pe变化幅度减37.5％，系统稳定性增强；抑制器作用于3机9节点系统时，其控制效果优于调速器控制，比有功功率振荡幅度减小了40％。附图说明图1为基于能量函数设计抑制器的方法的流程图。图2为单机无穷大系统的母线为本发明特征根随抑制器系数的变化规律图。图4为发电机采用本方法设计的抑制器后的功角δ仿真曲线为发电机采用本方法设计的抑制器后的有功功率Pe仿真曲线为发电机采用本方法设计的抑制器后，其q轴暂态电势E′q仿真曲线节点系统的母线采用本方法设计的抑制器后的功角δ仿真曲线采用本方法设计的抑制器后的有功功率Pe仿真曲线采用本方法设计的抑制器后的q轴暂态电势E′q仿真曲线采用本方法设计的抑制器后的机械功率Pm仿真曲线图。具体实施方式下面对本发明的具体实施方式进行描述，以便于本技术领域的技术人员理解本发明，但应该清楚，本发明不限于具体实施方式的范围，对本技术领域的普通技术人员来讲，只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内，这些变化是显而易见的，一切利用本发明构思的发明创造均在保护之列。该基于能量函数设计抑制器的方法是针对单机无穷大系统设计的抑制器，不过其也可以用于与其他系统中，单机无穷大系统的母线所示，该方法包括以下步骤：建立带励磁和汽门开度控制的电力系统(单机无穷大系统)状态空间模型：x=f(x)+g1(x)u+g2(x)w]]
其中，x＝(x1，x2，…，xn)T∈Rn，其为n维状态变量；u＝(u1，u2，…，um)T∈Rm，其为m维控制变量；w∈Rr，其为r维扰动变量；f(x)∈Rn，其为n维状态变量；g1(x)∈Rn×Rm，g2(x)∈Rn×Rr，g1(x)和g2(x)均为微分流形；为x的导数。在本发明的一个实施例中，当所述电力系统为单机无穷大系统时，其状态空间模型为：=-0=-DM(-0)+0MPm-0EqVsMxdsin+w1Eq=-1TdEq+1Td0xd-xdxdVscos+1Td0uf+w2Pm=1Ts(-Pm+Pm0)+1Tsup+w3]]
其中，w1、w2和w3是外部扰动，Pm是机械功率，Pm0是机械功率标定值，up是汽门开度控制输入，Ts是系统常数，δ是发电机的功角，ω是转子角速度，ω0＝2πf0，Eq′是q轴内部暂态电压，xd是d轴绕组自感抗，x′d是d轴绕组暂态阻抗，uf是发电机励磁电路电压，M是发电机转子转动惯量，D是阻尼常数，Td0是励磁电路时间常数，Td′是定子闭路时励磁绕组时间常数，Pm是机械功率，Vs是无穷大机组电压，x′dΣ＝x′d+xT+0.5xL，xT是变压器电抗，xL是输电线路电抗。对建立的电力系统状态空间模型进行泰勒展开，并基于常值实现方法建立能量函数：H(x)=i=1nxi(0)xiri(x1,x2,...,xn)dxi]]
其中，i＝1…n，f(x)∈Rn，f(x)为n维状态变量，(·)T表示转置，x＝(x1，x2，…，xn)T∈Rn，其为n维状态变量；(x1(0)，x2(0)，…，xn(0))T为任一初始点。在本发明的一个实施例中，上述对建立的电力系统状态空间模型进行泰勒展开，并基于常值实现方法建立能量函数步骤进一步包括：将建立的电力系统状态空间模型采用泰勒方式展开，选取其常数项和一次项；使电力系统状态空间模型对应状态矩阵为常数，状态矩阵元素满足的方程为：其中，f(x)∈Rn，其为n维状态变量，(·)T表示转置，Xn(x)为n维向量场。采用状态矩阵与Hamilton算子得到能量函数，Hamilton算子为：H=Nf(x)=X1T,...,XnTTf(x)=r1(x),...,rn(x)T]]
其中，i＝1…n，f(x)∈Rn，f(x)为n维状态变量，(·)T表示转置，x＝(x1，x2，…，xn)T∈Rn，其为n维状态变量；(x1(0)，x2(0)，…，xn(0))T为任一初始点。接着，基于能量函数表述电力系统状态空间模型：x=J(x)-R(x)H+g1(x)u+g2(x)w]]
其中，J(x)为n×n的反对称阵；R(x)为n×n的半正定矩阵；u＝(u1，u2，…，um)T∈Rm，其为m维控制变量；w∈Rr，其为r维扰动变量；f(x)∈Rn，其为n维状态变量；g1(x)∈Rn×Rm，g2(x)∈Rn×Rr，g1(x)和g2(x)均为微分流形；为x的导数。对采用能量函数表述后的电力系统状态空间模型设计抑制器：u=-1/2hT(x)h(x)+1/(22)Img1T(x)H]]
其中，h(x)为抑制器系数；(·)T表示转置；g1(x)∈Rn×Rm，其为微分流形；γ为设定的扰动抑制系数。当电力系统为单机无穷大系统时，针对其设计的抑制器也会根据单机无穷大系统的具体参数进行细化，此时对采用能量函数表述后的单机无穷大系统设计抑制器为：uf=Td02(h12+12)(0MxdVscos-Td00TdM(xd-xd)Eq+w0TdoM(xd-xd)V0)+Td0V0up=Ts02M(h22+12)(-Pm+Pm0)-Ts(-0)]]
其中，Td0是励磁电路时间常数，T′d是定子闭路时励磁绕组时间常数，h1和h2均是h(x)的对角线元素，γ为设定的扰动抑制系数，ω是转子角速度，ω0＝2πf0，xd是d轴绕组自感抗，x′d是d轴绕组暂态阻抗，M是发电机转子转动惯量，D是阻尼常数，E′q是q轴内部暂态电压，Pm是机械功率，Pm0是机械功率标定值，Td0V0为常值励磁。在本发明的一个实施例中，本方法采用特征根求取抑制器系数，加入抑制器后特征根状态矩阵为：A=0100-j2M-DM-j1M0M-j40-j300-Ts0-j5]]
其中，j6＝Ts，Ts是系统常数，δ是发电机的功角，ω是转子角速度，ω0＝2πf0，Eq′是q轴内部暂态电压，Eq′0为Eq′的初始值，xd是d轴绕组自感抗，x′d是d轴绕组暂态阻抗，M是发电机转子转动惯量，D是阻尼常数，Td0是励磁电路时间常数，T′d是定子闭路时励磁绕组时间常数，Vs是无穷大机组电压，x′dΣ＝x′d+xT+0.5xL，xT是变压器电抗，xL是输电线路电抗，γ为设定的扰动抑制系数，h1和h2均是h(x)的对角线元素。其中的设定的扰动抑制系数为：γ≥γ*，其中，M是发电机转子转动惯量，D是阻尼常数。下面结合图1至图11对采用本方法设计的抑制器的应用原理作详细的描述。将采用本方法设计的抑制器应用于单机无穷大系统：单机无穷大系统的母线所示，选取一组参数进行仿线，D，Vs，xd，x′d，x′dΣ，Td0，Td′，Ts
如您需求助技术专家，请点此查看客服电线.聚合物绝缘材料老化 2.电力系统可靠性分析
基于结合图割(Graph Cut)的Vibe改进算法的运动目标实时提取方法
