一种保证电力系统稳定性的采样周期优化方法与流程

本发明涉及电力系统控制技术领域，更具体地，涉及一种保证电力系统稳定性的采样周期优化方法。

背景技术：

现代电力系统趋于向规模化和复杂化的方向发展，电网规模不断扩大，输变电容量和电压等级不断提高。大区电网的互联可实现水火互济补偿和跨区域补偿调节，有利于电网错峰、机组互为备用、调节余缺、紧急事故支援等。虽然大电网互联带来了诸多好处，但是电网规模庞大、运行特性复杂、现有控制措施局限等原因使互联电网在发生故障时可能波及很广的范围并造成极大的危害。世界范围内已经发生了多起严重的互联电网大停电事故，如2003年美国东北部和加拿大联合电网大面积停电事故,2005年莫斯科大面积停电事故,2006年西欧大陆互联电网大面积停电事故，造成了巨大的经济损失和极大的社会影响。因此，大电网的安全和稳定己经成为电力系统领域研究的焦点。

为了实现大系统的安全稳定运行，降低大规模停电风险，有必要从全局范围的角度对电力系统运行进行监测、分析和控制。近年来，相量测量单元(pmu)技术和广域测量系统(wams)受到广泛的关注，为从全局角度分析和控制电力系统提供了可能。pmu可安装于电网各处且采集的数据可进行同步处理，加上飞速发展的通信网络和控制技术，以及现代电网稳定性的需要，基于同步相量技术的广域稳定控制己成为电力系统的重要发展方向。在保证电力系统大范围稳定的前提下，对系统各处进行数据采样分析是非常必要的。广域信号在由不同通信介质(如光纤、电话线、数字微波、卫星等)组成的wams通信网络中传输和处理，通信网络的信道容量是有限的。因此，在保证系统稳定性的前提下，尽可能的增大采样周期，降低采样频率，减少网络传输数据，从而达到降低本的目的。

技术实现要素：

本发明要解决的技术问题是针对现有技术不足和缺陷，提供一种保证电力系统稳定性的采样周期计算方法，估算出系统稳定条件下的最大采样周期。

为实现上述目的，本发明采用以下技术方案：

一种保证电力系统稳定的采样周期优化方法，其特征在于，包括以下步骤：

步骤1，建立n机电力系统的状态模型，具体表示为：其中x(t)∈rⁿ为系统状态，x(tk)表示采样tk时刻的状态向量；a、a1表示合适维数的常数矩阵；φ(t)表示初始状态向量函数；h＝tk+1-tk为采样周期；

步骤2，根据步骤1的电力系统状态模型构造如下形式的李亚谱诺夫泛函v(xt)，v(xt)具体表示为

步骤3，求解李亚谱诺夫泛函导数根据李亚谱诺夫稳定性理论，得到保证电力系统稳定的充分条件：如果存在对称矩阵p>0，r1>0，r2>0，z，以及合适维数矩阵x，q，mi，ni，yj(i＝1,2,3，j＝1,2,3,4)，使得

则系统渐近稳定，其中，

π4＝[γ^t00]^t，π5＝[γ^t0]^t，π6＝[0γ^t]^t，π7＝e1-e2，π8＝e5-e1，π9＝e1+e2-2e3，π10＝e1-e2-6e3+6e4，π11＝e5+e1-2e6，π12＝e5-e1+6e6-6e7，π13＝ae3+a1e2，π14＝ae6+a1e2，π15＝ae4+a1e2，π16＝ae7+a1e2，π17＝2e3-2e2，π18＝2e5-2e6，γ＝ae1+a1e2，ei＝[0n×(i-1)nin0n×(7-i)n](i＝1,2,…,7)；

步骤4，求解步骤3中的稳定性条件并最大化h，即可以得到保证系统稳定的最大采样周期hmax。

本发明的有益效果为在保证系统稳定的前提下，在电力系统向量测量单元pmu端进行信号采样记录时，计算并获取最大采样周期，使得pmu端信号发送频率减少，以减轻网络传输压力降低数据采样次数，减少数据传输量，降低运营成本,具有一定的前瞻性、经济性。

附图说明

图1为本发明实施例的优化方法流程图。

具体实施方式

下面结合具体实施方式对本发明作进一步的说明。

如图1所示，本发明实施提供一种保证电力系统稳定性的最大采样周期优化方法，包括以下步骤：

s101：建立多机电力系统模型。

对于单机无穷大系统，选取状态向量x1＝[δδδωδe'qδefd]^t，δ表示功角；ω表示角速度，e'q表示q轴电抗后暂态电势；efd表示励磁系统的输出电势。

假设在tk时刻采样获取状态采样数据为

在电力系统平衡点处，对系统方程线性化后的状态空间表达式为：

其中，ta为avr控制增益，ka为avr时间常数，m＝2h表示发电机的惯性时间常数矩阵，d表示阻尼系数矩阵，vs为无穷大母线端电压，xe表示线路阻抗；xd和x'd分别为发电机稳态、暂态电抗，ωb＝2πf0表示额定角速度，f0表示基准频率。

那么对于n机电力系统，状态变量表示为x＝[x1x2…xn]^t，其中xi＝[δδiδωiδe'qiδefdi]^t。经过推导可得出采样周期h下的n机电力系统的状态表达式为其中x(t)表示系统状态，x(tk)表示采样tk时刻的状态；a、a1表示合适维数的常数矩阵；φ(t)表示初始向量函数。h＝tk+1-tk为采样周期。

s102：基于上述n机电力系统状态模型，构造如下形式的李亚谱诺夫泛函：

s103：对上述的李亚谱诺夫泛函求导，获得n机电力系统稳定判据

定义ξ(t)＝[x(t)x(tk)v1(t)v2(t)x(tk+1)v3(t)v4(t)]^t，π4＝[γ^t00]^t，π5＝[γ^t0]^t，π6＝[0γ^t]^t，π7＝e1-e2，π8＝e5-e1，π9＝e1+e2-2e3，π10＝e1-e2-6e3+6e4，π11＝e5+e1-2e6，π12＝e5-e1+6e6-6e7，π13＝ae3+a1e2，π14＝ae6+a1e2，π15＝ae4+a1e2，π16＝ae7+a1e2，π17＝2e3-2e2，π18＝2e5-2e6，γ＝ae1+a1e2，ei＝[0n×(i-1)nin0n×(7-i)n](i＝1,2,…,7)；求李亚谱诺夫泛函导数得到

用以下两个积分不等式替换中的积分项：

并将以下四个零等式的右边部分加入到泛函导数中：

0＝2ξ^t(t)y1[(t-tk)∏3-∏7]ξ(t)，

0＝2ξ^t(t)y2[(tk+1-t)∏14-π8]ξ(t)，

0＝2ξ^t(t)y3[(t-tk)π15-π17]ξ(t)，

0＝2ξ^t(t)y4[(tk+1-t)π16-∏18]ξ(t)，

整理后得到其中,

若和则由李亚谱诺夫稳定性理论可知，系统渐近稳定。

由schur补引理可知，和分别等价于

和

综上，得到系统的稳定性判据。

s104：运用线性矩阵不等式工具箱，求解上述稳定性条件并最大化采样周期h，即可以得到保证系统稳定的最大采样周期hmax。