包含时变时滞的电动汽车与发电机组协调频率控制方法与流程

本发明涉及电力系统频率调节领域，尤其涉及一种包含时变时滞的电动汽车集群与发电机组协调频率控制方法。

背景技术：

在2015年巴黎气候大会上，来自世界各地的近200个国家和地区一致同意通过《巴黎协议》，为2020年后全球应对气候变化行动做出安排。在电力行业，为了达到节能减排的目的，大量的可再生能源正在或已经接入到电力系统中。在2020年英国风力发电的比例将达到总发电量的30％，而同期我国在“十三五”规划中也提出风电累积装机总量将达到210gw。由于可再生能源的间歇性导致负荷和发电功率不平衡，进而引起的调频问题是本发明的研究重点。

以电动汽车集群为代表的负荷需求侧响应具有绿色环保、响应速度快、能够向电网反供电vehicle-to-grid(v2g)、价格便宜等优点。但是，由于通讯延迟和电动汽车在地理上分布广泛，在采集电动汽车状态或对电动汽车单体发送控制信号时存在明显的延迟，这种延迟会对系统稳定性产生不利影响。因此，降低时滞对电力系统调频的危害并确定系统的稳定裕度是利用电动汽车集群进行调频的重要问题。此外，一些可再生能源如风力发电机和太阳能发电，当系统频率发生问题时不能向系统提供惯性支持，由此带来的系统惯性不确定性在本发明中也有所研究。

目前，计算系统稳定裕度方法比较流行的是自由权矩阵方法。它具有易于理解和使用简单、保守性小的特点，但由于其在计算时引入了大量需要求解的未知矩阵，因此具有计算效率慢的缺点。所以需要提出一种计算效率高且能同时处理时变时滞和参数不确定性的稳定判据。

技术实现要素：

针对现有技术，本发明提出一种包含时变时滞的电动汽车集群与发电机组协调频率控制方法。该方法是一种充分利用电动汽车集群参与负荷频率控制，使其与发电机组进行协调控制的电力系统调频新方法。在系统频率超过范围后，以电动汽车集群具有更高的参与度为宗旨进行调频。基于不含自由权矩阵的方法，本发明提出一种能同时处理时变时滞和系统惯性不确定性的系统稳定性新判据，分析了时变时滞和系统惯性的不确定性对系统频率的影响，并且利用电动汽车集群中的pi控制器来提高系统的稳定裕度。

为了解决上述技术问题，本发明提出的一种包含时变时滞的电动汽车与发电机组协调频率控制方法，电力系统包括多个发电机组和多个电动汽车集群，并包括以下步骤：

步骤一、构建电力系统负荷频率控制的数学模型，包括构建发电机负荷频率控制的动态数学模型及电动汽车负荷频率控制的动态数学模型，具体内容如下

步骤1-1)构建电力系统中发电机负荷频率控制的动态数学模型，首先，将电力系统中所有发电机组等效为单台发电机，则等效发电机的动态数学模型即为电力系统中发电机负荷频率控制的动态数学模型，表达式如下：

式(1)中：△f是电力系统的频率偏差；△y是等效发电机调速器阀门增量；△pr是等效发电机蒸汽轮机的输出；△pm是等效发电机的功率输出；是电力系统需要的等效发电机输出的功率值；req是等效发电机调速器的下垂系数；tg是等效发电机调速器的执行时间常数；t1和t2表示等效发电机调速器的暂态时间常数；tt表示等效发电机涡轮机的时间常数；d是电力系统中的阻尼系数；heq是电力系统的标称惯性常数；

步骤1-2)构建电动汽车负荷频率控制的动态数学模型，

第i个电动汽车集群负荷频率控制的动态数学模型为一阶模型，该一阶模型包括一个pi控制器；

式(2)中，i＝0,1,…,n，表示电动汽车集群的个数；rev是电动汽车集群的功率调节的下垂系数；tai是电动汽车集群的充放电时间常数；ni是电动汽车集群中的电动汽车数目；是电动汽车集群的输出功率；是电力系统需要电动汽车集群输出的功率值；kp是pi控制器的比例系数，ki是pi控制器的积分系数；kai是电动汽车集群的充放电系数；

式(3)至式(5)中：是电动汽车集群的充电系数；是电动汽车集群的放电系数；soci,j是电动汽车实时的电池状态；φsoci,j是soci,j的概率密度函数；

步骤1-3)在电动汽车集群负荷频率控制的动态数学模型中引入时变时滞后，电力系统频率偏差的导数为

式(6)中：di(t)为电力系统中的时变时滞，且同时满足：

0＝d0(t)≤di(t)≤τi,i＝1,2,…,n(7)

0≤τ1≤τ2≤…≤τn＝τmax(9)

式(7)、式(8)和式(9)中，d0(t)是i为0时的时变时滞，其值为0；τi是第i个时变时滞的上限，τmax是时变时滞上限的最大值；μi是第i个时变时滞随时间的变化率；

步骤1-4)考虑电力系统的惯性不确定性，得出系统惯性鲁棒值

式(10)中：heq是电力系统的标称惯性常数，keq∈[-1,0]是系统惯性鲁棒系数；

步骤1-5)以电力系统状态方程表达电力系统负荷频率控制的数学模型：

式(11)中：t表示时间变量；x＝[△y,△pr,△pm,△f,△pa1,△pa2,…,△pai,…,△pan-1,△pan,∫△f]^t为状态变量，属于n+5维实数向量为状态变量对时间的导数；ai，为系统系数矩阵，i＝0,1,…,n，其中，n表示时滞环节数目；di(t)为系统的时滞常数；x(t-τi)为时滞状态变量；h(t，ξ)为状态变量的历史轨迹；上述代数变量均属于实数域△ai为电力系统状态方程的不确定部分，且满足

[△a0,△a1,...,△an]＝df[e0,e1,...,en](12)

式(12)中，是常数矩阵；f是一个未知的实矩阵且满足

f^tf≤i(13)

式(13)中，符号t表示矩阵的转置，i表示单位矩阵；

综上，可得非时滞系数矩阵a0，时滞系数矩阵ai，i＝1,2,…,n的具体数值如下：

当i＝0时，将a0用分块矩阵的形式表示，即其中

当i＝1,2,…,n时，将ai用分块矩阵的形式表示，其中，

其中当i＝0,1,2,…,n时属于4×4维实数矩阵当i＝0,1,2,…,n时属于4×(n+1)维实数矩阵当i＝0,1,2,…,n时属于(n+1)×4维实数矩阵当i＝0,1,2,…,n时属于(n+1)×(n+1)维实数矩阵

常数矩阵d的具体数值如下：

将d用分块矩阵的形式表示，其中，

常数矩阵ei，i＝0,1,2,…,n的具体数值如下：

当i＝0时，将e0用分块矩阵的形式表示，即其中

当i＝1,2,…,n时，将ei用分块矩阵的形式表示，其中，

步骤二、由电力系统中的负荷需求和总发电量计算得到系统总不平衡功率△pt；通过收集到的电动汽车状态估计电动汽车个体的功率可控容量，继而到电动汽车能提供的功率响应容量其中是电动汽车能向电力系统提供的放电功率响应容量，是电动汽车能向电力系统提供的充电功率响应容量；通过发电机机组的备用信息估计发电机机组能提供的功率响应容量；

步骤三：对电动汽车和发电机组进行协调控制，将电力系统的频率偏差控制在允许的范围内；包括：

步骤3-1)判断电力系统的频率偏差值δf是否越限；若越限，则顺序执行步骤3-2)；否则执行步骤五；

步骤3-2)判断电力系统的频率偏差值δf是否越上限△f^u；若越上限，则顺序执行步骤3-3)；否则执行步骤3-4)；

步骤3-3)根据步骤二获得的系统总不平衡功率△pt判断电动汽车集群能提供的功率响应容量是否充足，若是，则令且否则，

执行步骤四；

步骤3-4)根据步骤二获得的系统总不平衡功率△pt判断电动汽车集群能提供的功率响应容量是否充足，若是，则令且否则，

步骤四：利用下述稳定性判据1确定电力系统负荷频率控制的数学模型在步骤3-3)或步骤3-4)情景下的稳定裕度；

判据1：式(11)所示的电力系统负荷频率控制的数学模型中，对于时变时滞的上限τi，如果存在一个标量ε>0,对称正定矩阵半正定矩阵i＝1,2,…,n,且该电力系统负荷频率控制的数学模型满足如下条件，则该电力系统是渐进稳定的：

ξ2n+2,2n+2＝-h

ξ2n+3,2n+3＝-εi

ξ1,2n+3＝pd

ξ2n+2,2n+3＝hd

步骤五、判断该电力系统当前时滞是否在稳定裕度之内，若是，返回步骤二，对下一时刻进行分析；

否则，利用pso算法求解式(19)所示的优化问题，得到电动汽车负荷频率控制动态数学模型中pi控制器的比例系数kp和积分系数ki，从而使电力系统当前时滞在稳定裕度之内，并返回步骤四；

进一步讲，步骤五中，利用pso算法求解式(19)所示的优化问题，得到电动汽车负荷频率控制动态数学模型中pi控制器的比例系数kp和积分系数ki的具体过程如下：

步骤5-1)利用pso算法随机产生pi控制器的比例系数kp和积分系数ki，并利用公式(20)和(21)来更新粒子的位置ki,i+1、kp,i+1；

步骤5-2)利用二分法和判据1计算ki,i+1、kp,i+1对应的时滞稳定裕度τi；计算ki,i+1、kp,i+1对应的δpev是否满足条件，如果满足，则执行步骤5-3，如果不满足返回步骤5-1；

步骤5-3)利用公式(22)和(23)更新粒子的速度v_ki,i+1、v_kp,i+1；

步骤5-4)比较步骤5-2)获得的时滞稳定裕度τi和当前粒子位置pbest对应的时滞稳定裕度以更新当前粒子位置pbest；比较全局粒子位置gbest和当前粒子位置pbest对应的时滞稳定裕度以更新全局粒子位置gbest；

步骤5-5)判断是否到达收敛条件或者最大迭代次数，如果是，则令pi控制器的积分系数最优值和pi控制器的比例系数最优值如果不是，返回步骤5-1以保证求得满足条件的时滞稳定裕度。

与现有技术相比，本发明的有益效果是：

本发明协调控制策略是在电动汽车参与系统调频服务的基础上，电动汽车和发电机组协调控制为核心的电力系统调频方法。为了充分利用电动汽车在频率响应中的优点，电动汽车具有较高的响应优先级。在协调控制策略中，本发明考虑了存在于电动汽车中的通信延迟和系统惯性的不确定性，提出了能够同时处理两者的系统渐进稳定性新判据。并且在电动汽车控制中心中采用了一组pi控制器来增加系统的稳定性。

附图说明

图1是本发明中实现电动汽车和发电机组协调控制的系统框图；

图2是本发明实现协调控制策略的原理图；

图3是本发明中发电机组控制器模型；

图4是本发明中电动汽车控制器模型。

具体实施方式

下面结合附图对本发明进行详细描述。

本发明的设计思路是，为实现电动汽车和发电机组协调控制，由系统负荷频率控制中心根据系统频率偏移量δf向电动汽车和发电机组下发控制信号和如图1所示，电动汽车控制器和发电机组控制器实际的响应功率△pev和△pm，与总不平衡功率△pt共同作用，经过系统惯性与负荷阻尼作用后，来达到调节系统频率偏差的作用。和由系统负荷平率控制中心通过本发明提出的协调控制策略进行计算，△pm和△pev由图3和图4所示的电动汽车和发电机组控制器模型计算所得。

本发明提出的一种包含时变时滞的电动汽车与发电机组协调频率控制方法，其中，电力系统包括多个发电机组和多个电动汽车集群。

本发明提出的协调控制策略由电动汽车负荷频率控制和发电机负荷频率控制实现。在系统频率偏差δf发生变化时，电动汽车和发电机组的输出功率响应按照本发明提出的协调控制策略进行控制。具体步骤如下：

步骤一、构建电力系统负荷频率控制的数学模型，包括构建发电机负荷频率控制的动态数学模型及电动汽车负荷频率控制的动态数学模型，具体内容如下

步骤1-1)构建电力系统中发电机负荷频率控制的动态数学模型，首先，将电力系统中所有发电机组等效为单台发电机，如图3所示，等效发电机包含调速器、涡轮机、惯性及阻尼，则等效发电机的动态数学模型即为电力系统中发电机负荷频率控制的动态数学模型，表达式如下：

式(1)中：△f是电力系统的频率偏差；△y是等效发电机调速器阀门增量；△pr是等效发电机蒸汽轮机的输出；△pm是等效发电机的功率输出；是电力系统需要的等效发电机输出的功率值；req是等效发电机调速器的下垂系数；tg是等效发电机调速器的执行时间常数；t1和t2表示等效发电机调速器的暂态时间常数；tt表示等效发电机涡轮机的时间常数；d是电力系统中的阻尼系数；heq是电力系统的标称惯性常数；

如图3所示，由协调控制策略得出系统需要发电机组输出的功率值协调控制策略在发电机组上的实现是通过将作为发电机组调速器的参考值来调节发电机组的输出功率，发电机组实际的输出功率记为△pm。由于发电机组的功率调节是通过增加或减小机组出力来实现的，具有响应速度较慢且成本较高的特点，所以系统调度中心对发电机组的控制优先级低于电动汽车。

步骤1-2)构建电动汽车负荷频率控制的动态数学模型，

如图4所示，第i个电动汽车集群负荷频率控制的动态数学模型为一阶模型，该一阶模型包括一个pi控制器；

式(2)中，i＝0,1,…,n，表示电动汽车集群的个数；rev是电动汽车集群的功率调节的下垂系数；tai是电动汽车集群的充放电时间常数；ni是电动汽车集群中的电动汽车数目；是电动汽车集群的输出功率；是电力系统需要电动汽车集群输出的功率值；kp是pi控制器的比例系数，ki是pi控制器的积分系数；kai是电动汽车集群的充放电系数；

式(3)至式(5)中：是电动汽车集群的充电系数；是电动汽车集群的放电系数；soci，j是电动汽车实时的电池状态；φsoci，j是soci，j的概率密度函数；

步骤1-3)在电动汽车集群负荷频率控制的动态数学模型中引入时变时滞后，电力系统频率偏差的导数为

式(6)中：di(t)为电力系统中的时变时滞，且同时满足：

0＝d0(t)≤di(t)≤τi,i＝1,2,…,n(7)

0≤τ1≤τ2≤…≤τn＝τmax(9)

式(7)、式(8)和式(9)中，d0(t)是i为0时的时变时滞，其值为0；τi是第i个时变时滞的上限，τmax是时变时滞上限的最大值；μi是第i个时变时滞随时间的变化率；

步骤1-4)考虑电力系统的惯性不确定性，得出系统惯性鲁棒值

式(10)中：heq是电力系统的标称惯性常数，keq∈[-1,0]是系统惯性鲁棒系数；

步骤1-5)以电力系统状态方程表达电力系统负荷频率控制的数学模型：

式(11)中：t表示时间变量；x＝[△y,△pr,△pm,△f,△pa1,△pa2,…,△pai,…,△pan-1,△pan,∫△f]^t为状态变量，属于n+5维实数向量为状态变量对时间的导数；ai，为系统系数矩阵，i＝0,1,…,n，其中，n表示时滞环节数目；di(t)为系统的时滞常数；x(t-τi)为时滞状态变量；h(t，ξ)为状态变量的历史轨迹；上述代数变量均属于实数域△ai为电力系统状态方程的不确定部分，且满足

[△a0,△a1,...,△an]＝df[e0,e1,...,en](12)

式(12)中，是常数矩阵；f是一个未知的实矩阵且满足

f^tf≤i(13)

式(13)中，符号t表示矩阵的转置，i表示单位矩阵；

综上，可得非时滞系数矩阵a0，时滞系数矩阵ai，i＝1,2,…,n的具体数值如下：

当i＝0时，将a0用分块矩阵的形式表示，即其中

当i＝1,2,…,n时，将ai用分块矩阵的形式表示，其中，

其中当i＝0,1,2,…,n时属于4×4维实数矩阵当i＝0,1,2,…,n时属于4×(n+1)维实数矩阵当i＝0,1,2,…,n时属于(n+1)×4维实数矩阵当i＝0,1,2,…,n时属于(n+1)×(n+1)维实数矩阵

常数矩阵d的具体数值如下：

将d用分块矩阵的形式表示，其中，

常数矩阵ei，i＝0,1,2,…,n的具体数值如下：

当i＝0时，将e0用分块矩阵的形式表示，即其中

当i＝1,2,…,n时，将ei用分块矩阵的形式表示，其中，

步骤二、由电力系统中的负荷需求和总发电量计算得到系统总不平衡功率△pt；通过收集到的电动汽车状态估计电动汽车个体的功率可控容量，继而到电动汽车能提供的功率响应容量其中是电动汽车能向电力系统提供的放电功率响应容量，是电动汽车能向电力系统提供的充电功率响应容量；通过发电机机组的备用信息估计发电机机组能提供的功率响应容量；

如图2所示，电动汽车的功率响应方式如下：

每个电动汽车集群控制一个电动汽车组。将电动汽车集群集合记作ωa＝{a1,a2,…,ai,…,an}，其中n是电动汽车集群个数；ai是第i个电动汽车集群，其中包含ni台电动汽车，记为ai＝{ei,1,ei,2,…,ei,j,…,ei,ni},i＝1,2,…,n.ei,j∈ai是ai中的第j台电动汽车。将电动汽车单体和集群的状态信息记为si,j和si＝[si,1,si,2,…,si,j,…,si,ni],i＝1,2,…,n，，其中表示电池容量、初始电池状态(soc)、实时soc、充电模式等信息。得到ei,j的可控功率记为其中表示ei,j作为电力系统的发电端可以向电网注入的有功功率；表示ei,j作为电力系统的负荷可以增加的有功功率负荷。进一步，计算出第i个电动汽车集群的可控功率

该式中，表示ai作为电力系统的发电端可以向电网注入的有功功率；表示ai作为电力系统的负荷可以增加的有功功率负荷。

进而，电动汽车的可控功率记为

该式中，表示电动汽车作为电力系统的发电端可以向电网注入的有功功率；表示电动汽车作为电力系统的负荷可以增加的有功功率负荷。

步骤三：对电动汽车和发电机组进行协调控制，将电力系统的频率偏差控制在允许的范围内；包括：

步骤3-1)判断电力系统的频率偏差值δf是否越限；若越限，则顺序执行步骤3-2)；否则执行步骤五；

步骤3-2)判断电力系统的频率偏差值δf是否越上限△f^u；若越上限，则顺序执行步骤3-3)；否则执行步骤3-4)；

步骤3-3)根据步骤二获得的系统总不平衡功率△pt判断电动汽车集群能提供的功率响应容量是否充足，若是，则令且否则，

执行步骤四；

步骤3-4)根据步骤二获得的系统总不平衡功率△pt判断电动汽车集群能提供的功率响应容量是否充足，若是，则令且否则，

如图2所示，电动汽车收到后，控制信号的分解到各个单台电动汽车的方式如下：

按以下方式分解为是电力系统需要电动汽车集群ai输出的功率值，进一步，可得电力系统需要每台电动汽车ei,j输出的功率值

电动汽车集群根据按公式(2)产生电动汽车集群的输出功率进而得到电动汽车的输出功率

在电动汽车参与频率响应的过程中，其通信存在时变时滞，且系统惯性存在不确定性。两种不确定性都会给系统的安全稳定运行造成威胁。为了对系统稳定性给予量化，将电力系统的时滞稳定裕度作为衡量电力系统稳定性的重要指标。稳定裕度越大，则系统的稳定区域越大，说明系统的稳定性越强。

步骤四：利用下述稳定性判据1确定电力系统负荷频率控制的数学模型在步骤3-3)或步骤3-4)情景下的稳定裕度；

判据1：式(11)所示的电力系统负荷频率控制的数学模型中，对于时变时滞的上限τi，如果存在一个标量ε>0,对称正定矩阵半正定矩阵i＝1,2,…,n,且该电力系统负荷频率控制的数学模型满足如下条件，则该电力系统是渐进稳定的：

ξ2n+2,2n+2＝-h

ξ2n+3,2n+3＝-εi

ξ1,2n+3＝pd

ξ2n+2,2n+3＝hd

步骤五、判断该电力系统当前时滞是否在稳定裕度之内，若是，返回步骤二，对下一时刻进行分析；

否则，利用pso算法求解式(19)所示的优化问题，得到电动汽车负荷频率控制动态数学模型中pi控制器的比例系数kp和积分系数ki，从而使电力系统当前时滞在稳定裕度之内，并返回步骤四；

步骤五中，利用pso算法求解式(19)所示的优化问题，得到电动汽车负荷频率控制动态数学模型中pi控制器的比例系数kp和积分系数ki的具体过程如下：

步骤5-1)利用pso算法随机产生pi控制器的比例系数kp和积分系数ki，并利用公式(20)和(21)来更新粒子的位置ki,i+1、kp,i+1；

步骤5-2)利用二分法和判据1计算ki,i+1、kp,i+1对应的时滞稳定裕度τi；计算ki,i+1、kp,i+1对应的δpev是否满足条件，如果满足，则执行步骤5-3，如果不满足返回步骤5-1；

步骤5-3)利用公式(22)和(23)更新粒子的速度v_ki,i+1、v_kp,i+1；

步骤5-4)比较步骤5-2)获得的时滞稳定裕度τi和当前粒子位置pbest对应的时滞稳定裕度以更新当前粒子位置pbest；比较全局粒子位置gbest和当前粒子位置pbest对应的时滞稳定裕度以更新全局粒子位置gbest；

步骤5-5)判断是否到达收敛条件或者最大迭代次数，如果是，则令pi控制器的积分系数最优值和pi控制器的比例系数最优值如果不是，返回步骤5-1以保证求得满足条件的时滞稳定裕度。

综上，本发明实现协调控制策略主要由两大控制器组成：电动汽车负荷频率控制器，简称为ev-lfc；发电机负荷频率控制器，简称为pp-lfc。如图2所示，在协调控制中，ev-lfc有较高的参与度并提供功率支撑，进一步考虑其时变时滞和系统惯性不确定性后得到系统动态模型，进而利用新判据对系统稳定裕度进行计算，并提出利用pi控制器提高系统时滞稳定性。