一种基于深度学习的风电调频逐步惯性控制方法

一种基于深度学习的风电调频逐步惯性控制方法
1.技术领域
2.本发明涉及风电调频技术领域，尤其涉及一种基于深度学习的风电调频逐步惯性控制方法。


背景技术：

3.为了应对气候、环境、能源等多方面问题，可再生能源在电力系统中的比例不断增大已经成为了一种趋势。风电作为一种重要的可再生能源，得到了稳定、持续的发展。风电在提供清洁能源的同时，也使得电力系统惯量下降，削弱了系统应对有功功率不平衡的能力，给系统的频率控制质量与频率稳定带来严峻挑战。风力涡轮机一般采用最大功率点跟踪控制，通过电力电子变流器接入电网，与系统频率解耦，无法通过释放或吸收能量响应功率偏差，不具备惯量响应特性，有功功率扰动下不能够为电网主动提供惯量支撑。为了提高系统频率控制的质量与效率，有必要让风力发电机参与这种控制，并通过这种方式补偿整个系统惯性的降低。近期，新的电网规范也要求风电场必须有助于电力系统频率控制。为了满足要求，当电网中出现有功功率不平衡时，可以使用逐步惯性控制参与电网频率控制，它能够迅速地提供短暂的功率支撑。
4.在使用逐步惯性控制进行调频过程中，风机通过释放转子动能向电网提供一段时间长度的增发功率，转子转速也随之下降。如果风机转子转速下降至其限定下限，风机将与电力系统断开连接。因此，风机只能通过释放有限的转子动能额外提供给电网一个短暂的功率支撑。逐步惯性控制也必须保证风机转子在达到最低转速前终止功率超发。然而，突然地终止功率超发会导致电网频率出现新的下降，也称为二次频率跌落（secondary frequency drop，sfd）。在逐步惯性控制下，调频效果好坏与二次频率跌落的大小有着复杂的联系，在保证良好调频效果的前提下，为了减轻二次频率跌落，可以采用不同的方法。
5.目前，有人提出将风电调频策略设计为分段函数，并使用粒子群算法获得不同风速区间下的分段函数值，可以有效改善sfd现象，但分段函数的使用也导致电网频率出现多次跌落现象，产生更多的频率波动，不利于频率快速恢复。还有人提出使用高压直流链路补偿风力涡轮机在终端的功率降低，由于高压直流链路的快速响应特性，功率衰减补偿技术在缓解二次频率跌落方面表现出良好的性能。然而，这种方法需要大量的投资。目前的逐步惯性控制下的改进方法，在面对不同负荷扰动事件时，难以快速、高效、经济地提供相应频率控制方案，发挥出逐步惯性控制的全部实力；难以满足风电调频时对于优秀逐步惯性控制的要求。


技术实现要素：

6.本发明提供了一种基于深度学习的风电调频逐步惯性控制方法，其特征在于，包括以下三个步骤：
步骤1、在时域仿真中使用金鹰优化算法获得最优风电调频逐步惯性控制参数并生成数据集；步骤2、基于堆叠式降噪自动编码器对数据集进行特征提取；步骤3、基于深度神经网络学习数据特征生成最优风电调频逐步惯性控制方案。
7.进一步地，所述步骤1的最优风电调频逐步惯性控制参数是在风电调频中使用逐步惯性控制策略时所使用的。
8.逐步惯性控制（stepwise inertial control，sic）是一种风电参与电力系统频率调整的控制策略，包括两个主要阶段：短时超发阶段和转速恢复阶段。图2、图3分别表示固定风电占比、恒风速的条件下使用逐步惯性控制策略时，风机转子转速与输出功率的关系以及逐步惯性控制过程中输出功率的变化。具体参与调频过程如下：1）风机正常工作在mppt曲线上的a点，输出的电磁功率为，转子转速为，输出功率。t0时刻，系统突然出现功率不平衡，导致频率下降，风机监测到频率下降超过死区限值，随即切换到逐步惯性控制模式，立即增加输出功率，输出功率，输出功率从a点到提升到b点，进入短时超发阶段。
9.2）在短时超发阶段，风机的输出功率在一段时间内保持不变，即输出功率一直为，从b点持续到c点。此阶段风机输出的电磁功率大于机械功率，转子随即减速，此时风机的转子运动方程为：
ꢀꢀꢀꢀꢀꢀꢀꢀꢀ
（1）式中：指的是风机惯性常数。
10.值得注意的是，短时超发阶段必须在转子过度减速之前终止，以避免风机发生转子失速事件。因此，短时超发阶段的持续时间的取值需保证转子转速在时刻未达到最小转速，一般为0.7p.u.。
11.3）时刻，短时超发阶段结束，风机输出功率突然下降，此时，转子转速为，输出功率由c点突降到d点，风机进入到恢复阶段。
12.4）风机进入转速恢复阶段，此时风机在最大功率点跟踪模式下运行，风机输出功率为，产生的系统功率缺额为：，该阶段风机的机械功率大于电磁功率，转子开始加速，经过的恢复时间，在时刻转子转速由恢复到，风机输出的电磁功率也将随着转子转速的上升恢复至，输出功率沿mppt轨迹从d点缓慢回到a点。
13.进一步地，所述步骤1中获得最优风电调频逐步惯性控制参数是在时域仿真中使用金鹰优化算法获得的。
14.金鹰优化算法（golden eagle optimizer，geo）的灵感来自于金鹰的智慧，具有寻优能力强，收敛速度快等优点。金鹰狩猎时在不同的盘旋轨迹阶段调整速度。在狩猎的初期，其表现出更多的巡游、寻找猎物的倾向，在后期，其表现出更多的攻击倾向。金鹰通过调整这两个部分使得其在最短的时间、可行的区域内捕捉到最佳的猎物。将此类行为编写为数学模型，并以突出探索和全局优化的方法开发。
15.盘旋动作：金鹰优化算法基于金鹰的盘旋运动，每只金鹰都记得它目前所到过的最佳位置。金鹰想要攻击猎物和巡航寻找更好食物。在每次迭代中，每只金鹰随机选择另一只金鹰的猎物，并在其目前访问过的最佳位置盘旋。金鹰也可以选择圈出自己的记忆。
16.猎物选择：在每次迭代中，每只金雕必须选择一个猎物来执行巡航和攻击操作。在金鹰优化算法中，猎物被建模为金雕群迄今找到的最佳解决方案。每只金鹰都能记住目前为止找到的最佳解决方案。在每次迭代中，每个搜索代理从整个群体的记忆中选择目标猎物。然后计算每只金鹰选定相应猎物的攻击和巡航向量。如果通过攻击和巡航向量计算的新位置优于内存中的前一个位置，则更新内存。猎物选择策略在金鹰优化算法中起着重要作用。选择可以用一种基本的方式进行，每只金鹰只在自己的记忆中选择猎物。为了使金鹰更好地探索这片土地，提出了随机的一对一映射方案，每只金鹰从其他群体成员的记忆中随机选择当前迭代中的猎物。值得注意的是，被选中的猎物不一定是最近或最远的猎物。在这个方案中，记忆中的每一个猎物，都被分配或映射给唯一的金鹰。然后每只金鹰对选定的猎物进行攻击和巡航。
17.攻击行为：攻击行为可以通过一个矢量来模拟，从金鹰当前的位置开始，到金鹰记忆中猎物的位置结束。通过式（2）可以计算出金鹰的攻击矢量：
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
（2）式中：为第i只金鹰的攻击向量；为金鹰f目前所到达的最佳地点(猎物)；为第i只金鹰的当前位置。
18.巡航（探索）行为：巡航向量是在攻击向量的基础上计算出的。巡航向量是圆的切向量，垂直于攻击向量。巡航也可以被认为是金鹰相对于猎物的线性速度。n维巡航矢量位于与圆相切的超平面内，所以，为了计算巡航矢量，必先计算相切超平面的方程。n维超平面的方程可以由超平面上的任意一点与一个垂直与该超平面的向量确定，这个向量称为超平面的法向量。超平面方程在n维空间中的标量形式如式（3）所示：
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
（3）式中：为法向量；为变量向量；为超平面上的任意点；。把（攻击向量）看作超平面的法线，那么（在迭代中的金鹰巡航向量）就可以根据式（4）表示出它所属的超平面。
19.ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
（4）式中：为攻击向量；为决策/设计变量向量；为被选中的猎物的位置。
20.为了在巡航超平面上找到一个随机矢量，我们必须先在这个超平面上找到一个随机目的点c点，而不是我们已有的点(金鹰的当前位置)。注意其巡航矢量的起点是金鹰的当前位置。由于超平面比它们的环境空间小一维，不能简单地生成一个随机的1
×
n点。n维空间中的简单随机点不能保证位于巡航超平面上。位于n维巡航超平面上的一个新点具有n-1个自由度，意味着n-1个维度可以自由选择，但超平面方程决定了最后的1维度，如式（3）所示。必须选择最后一个维数使其满足超平面方程，因此，有n-1个自由变量和1个固定变量。
使用以下步骤找到一个位于金鹰的巡航超平面上的随机n维目的地点c，具体步骤如下：步骤1）从n个变量中随机选择一个变量作为固定变量；步骤2）将随机值赋给除第k个变量以外的所有变量，因为第k个变量是固定的。
21.步骤3）用方程（5）求出固定变量的值。
22.ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
（5）式中：为目标点c的第k个元素；为攻击向量的第j个元素；k为固定变量的编号。在巡航超平面上找到随机目标点，巡航超平面上目的点的一般表示如式（6）所示：
ꢀꢀ
（6）转移到新位置：金鹰的位移由攻击和向量组成。我们将金鹰迭代的步长向量定义为方程（7）和（8）。
23.ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
（7）
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
（8）式中：为攻击系数；为巡航系数；、为[0,1]内的随机向量。因此，通过将迭代t中的步进向量与迭代t中的位置相加，即可计算出迭代t 1中金雕的位置：
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
（9）式中：为金鹰的第t 1次的位置；为金鹰的第t次的位置；为金鹰的移动的步长大小。
[0024]
如果金鹰新位置的适应性比它记忆中的位置更好，这只鹰的记忆会随着新位置而更新。否则，记忆保持完整，但是鹰会在新的位置上。在新的迭代中，每只金鹰从种群中随机选择一只金鹰，围绕其最佳访问位置旋转，计算攻击向量、巡航向量，最后计算下一次迭代的步长向量和新位置。这个循环将一直执行，直到满足任何终止条件。式（7）中有两个系数，即攻击系数和巡航系数，用来控制攻击向量、巡航向量对步长向量的影响，具体见式（10）：
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
（10）式中：t表示当前迭代；t表示最大迭代；和为攻击倾向的初始值和终止值；和为巡航倾向的初始值和终止值：进一步地，所述步骤1中在时域仿真中使用金鹰优化算法获得最优风电调频逐步惯性控制参数并生成数据集的基本步骤如下：步骤1）初始化金鹰数量；
步骤2）计算适应度函数，初始化群体记忆；步骤3）初始化、；步骤4）计算式（10），更新、；步骤5）计算式（2），从种群的记忆计算攻击向量中随机选择猎物；步骤6）计算巡航向量（式（3）-（6）)，计算步长向量(式（7）-（9）)更新位置（式（9）)评估新位置的适应度函数；步骤7）更新最优解及最优位置；步骤8）判断是否满足最大迭代次数，若满足，则输出最优金鹰位置和全局最优解，否则，返回步骤4重新迭代计算；进一步地，所述步骤2中的基于堆叠式降噪自动编码器对数据集进行特征提取是由堆叠式降噪自动编码器进行的。
[0025]
自动编码器（autoencoder，ae）如图4所示，自动编码器由编码器和解码器组成，编码器、解码器具有隐藏层:通过编码器将输入转换为隐含层中的潜在表征，然后通过解码器将内部表征转换为输出，输出相当于是对于输入的重构，要尽可能接近输入。
[0026]
选择输入样本集合，由n组样本组成，其中，n为样本组数，n为每组样本中样本个数。设隐含层特征向量集合为，由n组特征向量由n组特征向量组成，其中，m为每组特征向量中向量个数，则与的编码关系为
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
（11）式中：为输入层与隐含层的权值矩阵；为输入层与隐含层偏置矩阵；为编码器的神经元激活函数，通常采用sigmoid函数，其具有良好的特征辨识度。
[0027]
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
（12）式中：为输入向量。
[0028]
解码器是编码器的逆运算，以隐含层的特征向量作为输入向量，设为输出向量集合，共有n组，维数为n,则解码器的表达式为
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
（13）式中：为隐含层与输出层的权值矩阵；为隐含层与输出层的偏置矩阵；为解码器的神经元激活函数。
[0029]
自动编码器通过最小化输出向量与输入向量之间的重构误差来达到特征学习的目的，利用梯度下降算法不断调整网络权值和偏置，降低重构误差。但是自动编码器的学习可能仅仅简单地保留原始输入数据信息，不能确保获得一种有效的特征信息表达。
[0030]
降噪自动编码器（denoising autoencoder，dae）是对自动编码器的改进，其架构如图5所示。在自动编码器的基础上，通过将一部分输入值中增加噪声，训练它以恢复原始的无噪声输入。选用降噪自动编码器后隐含层可学习出更加有效的数据表征，挖掘更深的隐藏信息。
[0031]
堆叠式降噪自动编码器（stacked denoising autoencoder，sdae）由多个dae组成，如图6所示。wn为第n-1层隐含层与第n层隐含层的权重矩阵，bn为第n-1层隐含层与第n
层隐含层的偏置矩阵。通过无监督贪婪逐层预训练和有监督细化训练两个步骤完成训练，实现特征学习。第一层dae的输入为原始数据，隐含层输出数据作为上层dae输入数据，且每一层的训练前都要对输入数据添加噪声。堆叠式降噪自动编码器通过训练可以高效地提取数据的高阶特征，更优地拟合复杂函数，提升获得网络参数的速度，提高神经网络的深层特征学习能力。
[0032]
进一步地，所述步骤2中基于堆叠式降噪自动编码器对数据集进行特征提取的基本步骤如下：步骤1）对在时域仿真中使用金鹰优化算法获得的数据集进行归一化处理；步骤2）划分训练集与验证集，将80%的原始数据作为训练集，剩下的20%作为测试数据，生成训练和测试数据；步骤3）构建一个堆叠式降噪自动编码器模型，选择隐含层层数及神经元个数；步骤4）逐层贪婪无监督预训练，以添加噪声后的数据作为输入数据，并采用dropout技术，随机选取部分神经元暂时不工作；步骤5）有监督细化训练，对网络的权重和偏置进行微调，直到迭代次数达到设定值。
[0033]
进一步地，所述步骤3中的基于深度神经网络学习数据特征生成最优风电调频逐步惯性控制方案是在深度神经网络的基础上进行的。
[0034]
深度神经网络（deep neural network，dnn）是具有多个隐藏层的神经网络，又被称为深度前馈网络（feedforward neural network，dfn）。按不同层的位置划分，dnn内部的神经网络层可以分为：输入层，隐藏层和输出层，一般第一层是输入层，最后一层是输出层，而中间的层数都是隐藏层。层与层之间是全连接的，即第i层的任意一个神经元一定与第i 1层的任意一个神经元相连，基本结构如图7所示。复杂的dnn从局部模型来看，同感知机一样即一个线性关系，再加上一个激活函数。
[0035]
进一步地，所述步骤3中基于深度神经网络学习数据特征生成最优风电调频逐步惯性控制方案的基本步骤如下：步骤1）将堆叠式降噪自动编码器提取的特征当作深度神经网络的输入；步骤2）划分训练集与验证集，将80%的原始数据作为训练集，剩下的20%作为测试数据，生成训练和测试数据；步骤3）构建一个深度神经网络模型，选择隐含层层数和每个隐含层的神经元个数；步骤4）使用adam优化网络的输入权重矩阵、隐含层特征向量与输出权重矩阵，并对隐含层选择性增加正则化；步骤5）将mse作为评价指标，判断其是否不再降低，若是满足则完成训练，否则返回步骤3；步骤6）保存训练好的模型，当下次遇到负荷扰动事件时，直接使用模型生成最优风电调频逐步惯性控制方案。
[0036]
本发明采用以上技术方案与现有技术相比，具有以下优点：1、在进行风电调频逐步惯性控制时，与考虑因素单一的风电调频策略相比较，本专利所用方法综合考虑不同风速、风电占比、负荷扰动量对系统频率特性的影响，参考因素
更加全面，能够高效地获得多场景下的最优逐步惯性控制参数，满足了风电调频在电力系统中不同场景下的使用需求，有着良好的泛化能力。
[0037]
2、本专利所用方法通过堆叠式降噪自动编码器的特征提取以及深度神经网络学习得到的最优逐步惯性参数具有极高的准确性，满足了风电调频实际操作时快速实现频率控制的需求。与时域仿真方法对比，所用方法基于深度学习计算速度快的优点，极大地节省了进行在线决策的频率控制时间，提高了在线频率控制的效率。
[0038]
3、与其他逐步惯性控制方法相比，本专利所用方法通过上述技术方案，可获得极高精确度的最优逐步惯性控制参数，得到更佳的调频效果。能够高效地参与电力系统的频率控制，及时抑制频率跌落，有利于补偿系统惯性，提高系统频率控制的质量与效率，对电力系统运行的安全性、稳定性有着重大意义。
附图说明
[0039]
图1为摘要附图；图2为恒风速、固定风电占比的条件下逐步惯性控制策略风机转子转速与输出功率的关系图；图3为逐步惯性控制过程中输出功率随时间的变化图；图4为自动编码器结构图；图5为降噪自动编码器结构图；图6为堆叠式降噪自动编码器结构图；图7为深度神经网络结构图；图8为ieee 9节点测试系统；图9为基于深度学习的风电调频逐步惯性控制方法预测结果示意图；图10为不同逐步惯性控制结果示意图。
[0040]
具体实施方式
[0041]
下面结合附图对本发明的技术方案做进一步的详细说明：本发明可以以许多不同的形式实现，而不应当认为限于这里所述的实施例。相反，提供这些实施例以便使本公开透彻且完整，并且将向本领域技术人员充分表达本发明的范围。
[0042]
本发明专利提供了一种基于深度学习的风电调频逐步惯性控制方法，其特征在于，包括以下三个步骤：步骤1、在时域仿真中使用金鹰优化算法获得最优风电调频逐步惯性控制参数并生成数据集；步骤2、基于堆叠式降噪自动编码器对数据集进行特征提取；步骤3、基于深度神经网络学习数据特征生成最优风电调频逐步惯性控制方案。
[0043]
步骤1的最优风电调频逐步惯性控制参数是在风电调频中使用逐步惯性控制策略时所使用的。
[0044]
逐步惯性控制（stepwise inertial control，sic）是一种风电参与电力系统频率调整的控制策略，包括两个主要阶段：短时超发阶段和转速恢复阶段。图2、图3分别表示固定风电占比、恒风速的条件下使用逐步惯性控制策略时，风机转子转速与输出功率的关系
以及逐步惯性控制过程中输出功率的变化。具体参与调频过程如下：1）风机正常工作在mppt曲线上的a点，输出的电磁功率为，转子转速为，输出功率。t0时刻，系统突然出现功率不平衡，导致频率下降，风机监测到频率下降超过死区限值，随即切换到逐步惯性控制模式，立即增加输出功率，输出功率，输出功率从a点到提升到b点，进入短时超发阶段。
[0045]
2）在短时超发阶段，风机的输出功率在一段时间内保持不变，即输出功率一直为，从b点持续到c点。此阶段风机输出的电磁功率大于机械功率，转子随即减速，此时风机的转子运动方程为：
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
（1）式中：指的是风机惯性常数。
[0046]
值得注意的是，短时超发阶段必须在转子过度减速之前终止，以避免风机发生转子失速事件。因此，短时超发阶段的持续时间的取值需保证转子转速在时刻未达到最小转速，一般为0.7p.u.。
[0047]
3）时刻，短时超发阶段结束，风机输出功率突然下降，此时，转子转速为，输出功率由c点突降到d点，风机进入到恢复阶段。
[0048]
4）风机进入转速恢复阶段，此时风机在最大功率点跟踪模式下运行，风机输出功率为，产生的系统功率缺额为：，该阶段风机的机械功率大于电磁功率，转子开始加速，经过的恢复时间，在时刻转子转速由恢复到，风机输出的电磁功率也将随着转子转速的上升恢复至，输出功率沿mppt轨迹从d点缓慢回到a点。
[0049]
步骤1中获得最优风电调频逐步惯性控制参数是在时域仿真中使用金鹰优化算法获得的。
[0050]
金鹰优化算法（golden eagle optimizer，geo）的灵感来自于金鹰的智慧，具有寻优能力强，收敛速度快等优点。金鹰狩猎时在不同的盘旋轨迹阶段调整速度。在狩猎的初期，其表现出更多的巡游、寻找猎物的倾向，在后期，其表现出更多的攻击倾向。金鹰通过调整这两个部分使得其在最短的时间、可行的区域内捕捉到最佳的猎物。将此类行为编写为数学模型，并以突出探索和全局优化的方法开发。
[0051]
盘旋动作：金鹰优化算法基于金鹰的盘旋运动，每只金鹰都记得它目前所到过的最佳位置。金鹰想要攻击猎物和巡航寻找更好食物。在每次迭代中，每只金鹰随机选择另一只金鹰的猎物，并在其目前访问过的最佳位置盘旋。金鹰也可以选择圈出自己的记忆。
[0052]
猎物选择：在每次迭代中，每只金雕必须选择一个猎物来执行巡航和攻击操作。在金鹰优化算法中，猎物被建模为金雕群迄今找到的最佳解决方案。每只金鹰都能记住目前为止找到的最佳解决方案。在每次迭代中，每个搜索代理从整个群体的记忆中选择目标猎物。然后计算每只金鹰选定相应猎物的攻击和巡航向量。如果通过攻击和巡航向量计算的新位置优于内存中的前一个位置，则更新内存。猎物选择策略在金鹰优化算法中起着重要作用。选择可以用一种基本的方式进行，每只金鹰只在自己的记忆中选择猎物。为了使金鹰
更好地探索这片土地，提出了随机的一对一映射方案，每只金鹰从其他群体成员的记忆中随机选择当前迭代中的猎物。值得注意的是，被选中的猎物不一定是最近或最远的猎物。在这个方案中，记忆中的每一个猎物，都被分配或映射给唯一的金鹰。然后每只金鹰对选定的猎物进行攻击和巡航。
[0053]
攻击行为：攻击行为可以通过一个矢量来模拟，从金鹰当前的位置开始，到金鹰记忆中猎物的位置结束。通过式（2）可以计算出金鹰的攻击矢量：
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
（2）式中：为第i只金鹰的攻击向量；为金鹰f目前所到达的最佳地点(猎物)；为第i只金鹰的当前位置。
[0054]
巡航（探索）行为：巡航向量是在攻击向量的基础上计算出的。巡航向量是圆的切向量，垂直于攻击向量。巡航也可以被认为是金鹰相对于猎物的线性速度。n维巡航矢量位于与圆相切的超平面内，所以，为了计算巡航矢量，必先计算相切超平面的方程。n维超平面的方程可以由超平面上的任意一点与一个垂直与该超平面的向量确定，这个向量称为超平面的法向量。超平面方程在n维空间中的标量形式如式（3）所示：
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
（3）式中：为法向量；为变量向量；为超平面上的任意点；。把（攻击向量）看作超平面的法线，那么（在迭代中的金鹰巡航向量）就可以根据式（4）表示出它所属的超平面。
[0055]
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
（4）式中：为攻击向量；为决策/设计变量向量；为被选中的猎物的位置。
[0056]
为了在巡航超平面上找到一个随机矢量，我们必须先在这个超平面上找到一个随机目的点c点，而不是我们已有的点(金鹰的当前位置)。注意其巡航矢量的起点是金鹰的当前位置。由于超平面比它们的环境空间小一维，不能简单地生成一个随机的1
×
n点。n维空间中的简单随机点不能保证位于巡航超平面上。位于n维巡航超平面上的一个新点具有n-1个自由度，意味着n-1个维度可以自由选择，但超平面方程决定了最后的1维度，如式（3）所示。必须选择最后一个维数使其满足超平面方程，因此，有n-1个自由变量和1个固定变量。使用以下步骤找到一个位于金鹰的巡航超平面上的随机n维目的地点c，具体步骤如下：步骤1：从n个变量中随机选择一个变量作为固定变量；步骤2：将随机值赋给除第k个变量以外的所有变量，因为第k个变量是固定的。
[0057]
步骤3：用方程（5）求出固定变量的值。
[0058]
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
（5）式中：为目标点c的第k个元素；为攻击向量的第j个元素；k为固定变量的编号。在巡航超平面上找到随机目标点，巡航超平面上目的点的一般表示如式（6）所示：
ꢀꢀ
（6）转移到新位置：金鹰的位移由攻击和向量组成。我们将金鹰迭代的步长向量定义为方程（7）和（8）。
[0059]
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
（7）
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
（8）式中：为攻击系数；为巡航系数；、为[0,1]内的随机向量。因此，通过将迭代t中的步进向量与迭代t中的位置相加，即可计算出迭代t 1中金雕的位置：
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
（9）式中：为金鹰的第t 1次的位置；为金鹰的第t次的位置；为金鹰的移动的步长大小。
[0060]
如果金鹰新位置的适应性比它记忆中的位置更好，这只鹰的记忆会随着新位置而更新。否则，记忆保持完整，但是鹰会在新的位置上。在新的迭代中，每只金鹰从种群中随机选择一只金鹰，围绕其最佳访问位置旋转，计算攻击向量、巡航向量，最后计算下一次迭代的步长向量和新位置。这个循环将一直执行，直到满足任何终止条件。式（7）中有两个系数，即攻击系数和巡航系数，用来控制攻击向量、巡航向量对步长向量的影响，具体见式（10）：
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
（10）式中：t表示当前迭代；t表示最大迭代；和为攻击倾向的初始值和终止值；和为巡航倾向的初始值和终止值：进一步地，所述步骤1中在时域仿真中使用金鹰优化算法获得最优风电调频逐步惯性控制参数并生成数据集的基本步骤如下：步骤1）初始化金鹰数量；步骤2）计算适应度函数，初始化群体记忆；步骤3）初始化、；步骤4）计算式（10），更新、；步骤5）计算式（2），从种群的记忆计算攻击向量中随机选择猎物；步骤6）计算巡航向量（式（3）-（6）)，计算步长向量(式（7）-（9）)更新位置（式（9）)评估新位置的适应度函数；步骤7）更新最优解及最优位置；步骤8）判断是否满足最大迭代次数，若满足，则输出最优金鹰位置和全局最优解，否则，返回步骤4重新迭代计算；
进一步地，所述步骤2中的基于堆叠式降噪自动编码器对数据集进行特征提取是由堆叠式降噪自动编码器进行的。
[0061]
自动编码器（autoencoder，ae）如图4所示，自动编码器由编码器和解码器组成，编码器、解码器具有隐藏层:通过编码器将输入转换为隐含层中的潜在表征，然后通过解码器将内部表征转换为输出，输出相当于是对于输入的重构，要尽可能接近输入。
[0062]
选择输入样本集合，由n组样本组成，其中，n为样本组数，n为每组样本中样本个数。设隐含层特征向量集合为，由n组特征向量由n组特征向量组成，其中，m为每组特征向量中向量个数，则与的编码关系为
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
（11）式中：为输入层与隐含层的权值矩阵；为输入层与隐含层偏置矩阵；为编码器的神经元激活函数，通常采用sigmoid函数，其具有良好的特征辨识度。
[0063]
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
（12）式中：为输入向量。
[0064]
解码器是编码器的逆运算，以隐含层的特征向量作为输入向量，设为输出向量集合，共有n组，维数为n,则解码器的表达式为
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
（13）式中：为隐含层与输出层的权值矩阵；为隐含层与输出层的偏置矩阵；为解码器的神经元激活函数。
[0065]
自动编码器通过最小化输出向量与输入向量之间的重构误差来达到特征学习的目的，利用梯度下降算法不断调整网络权值和偏置，降低重构误差。但是自动编码器的学习可能仅仅简单地保留原始输入数据信息，不能确保获得一种有效的特征信息表达。
[0066]
降噪自动编码器（denoising autoencoder，dae）是对自动编码器的改进，其架构如图5所示。在自动编码器的基础上，通过将一部分输入值中增加噪声，训练它以恢复原始的无噪声输入。选用降噪自动编码器后隐含层可学习出更加有效的数据表征，挖掘更深的隐藏信息。
[0067]
堆叠式降噪自动编码器（stacked denoising autoencoder，sdae）由多个dae组成，如图6所示。wn为第n-1层隐含层与第n层隐含层的权重矩阵，bn为第n-1层隐含层与第n层隐含层的偏置矩阵。通过无监督贪婪逐层预训练和有监督细化训练两个步骤完成训练，实现特征学习。第一层dae的输入为原始数据，隐含层输出数据作为上层dae输入数据，且每一层的训练前都要对输入数据添加噪声。堆叠式降噪自动编码器通过训练可以高效地提取数据的高阶特征，更优地拟合复杂函数，提升获得网络参数的速度，提高神经网络的深层特征学习能力。
[0068]
步骤2中基于堆叠式降噪自动编码器对数据集进行特征提取的基本步骤如下：步骤1）对在时域仿真中使用金鹰优化算法获得的数据集进行归一化处理；步骤2）划分训练集与验证集，将80%的原始数据作为训练集，剩下的20%作为测试数据，生成训练和测试数据；
步骤3）构建一个堆叠式降噪自动编码器模型，选择隐含层层数及神经元个数；步骤4）逐层贪婪无监督预训练，以添加噪声后的数据作为输入数据，并采用dropout技术，随机选取部分神经元暂时不工作；步骤5）有监督细化训练，对网络的权重和偏置进行微调，直到迭代次数达到设定值。
[0069]
步骤3中的基于深度神经网络学习数据特征生成最优风电调频逐步惯性控制方案是在深度神经网络的基础上进行的。
[0070]
深度神经网络（deep neural network，dnn）是具有多个隐藏层的神经网络，又被称为深度前馈网络（feedforward neural network，dfn）。按不同层的位置划分，dnn内部的神经网络层可以分为：输入层，隐藏层和输出层，一般第一层是输入层，最后一层是输出层，而中间的层数都是隐藏层。层与层之间是全连接的，即第i层的任意一个神经元一定与第i 1层的任意一个神经元相连，基本结构如图7所示。复杂的dnn从局部模型来看，同感知机一样即一个线性关系，再加上一个激活函数。
[0071]
步骤3中基于深度神经网络学习数据特征生成最优风电调频逐步惯性控制方案的基本步骤如下：步骤1）将堆叠式降噪自动编码器提取的特征当作深度神经网络的输入；步骤2）划分训练集与验证集，将80%的原始数据作为训练集，剩下的20%作为测试数据，生成训练和测试数据；步骤3）构建一个深度神经网络模型，选择隐含层层数和每个隐含层的神经元个数；步骤4）使用adam优化网络的输入权重矩阵、隐含层特征向量与输出权重矩阵，并对隐含层选择性增加正则化；步骤5）将mse作为评价指标，判断其是否不再降低，若是满足则完成训练，否则返回步骤3；步骤6）保存训练好的模型，当下次遇到负荷扰动事件时，直接使用模型生成最优风电调频逐步惯性控制方案。
[0072]
将ieee 9节点系统作为算例测试系统，将风力涡轮机模型连接到系统中的线路l3上，如图8所示。所涉及到的研究都是在风力涡轮机以低于额定功率仿真运行中得到的。为了满足在不同场景下使用此方法的需求，考虑风速、风电占比、负荷扰动量对系统频率特性的影响。分别从这三个方面对sic进行参数设置。此算例中，设置风速4m/s至11m/s，每种情况递增1m/s。设置风电占比从5%至60%，每种情况递增5%。负荷扰动量设置从1.005至1.25，每种情况递增0.005。共计4800种情况。最优sic参数组合预测结果示意图如图9所示。并在典型场景中添加与其他调频控制方法的效果对比。系统中的场景为风速5m/s、风电占比25%、负荷扰动量0.035p.u.。在此场景中，与其他sic控制以及无风电参与的频率控制进行对比，频率控制效果见图10所示。
[0073]
以上实施方式仅为说明本发明的技术思想，并不用于限制本发明的保护范围，凡是按照本发明提出的技术思想，在本发明技术方案基础上所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。