本发明实施方式提供一种用于对光伏发电站进行无功优化的方法、系统及存储介质，属于光伏发电站调控技术领域。所述方法包括：接收从发无功/电压指令中获取的无功目标值；采集光伏发电站的每个节点电压，其中，节点电压包括光伏单元节点电压和并网点电压；对每个节点电压进行分析；根据分析结果来控制包括无功补偿装置和光伏逆变器的无功电源的动作先后顺序。该方法、系统及存储介质能够解决现有技术中光伏发电站输出无功功率时节点电压偏差较大的问题。
1.一种用于对光伏发电站进行无功优化的方法，其特征在于，所述方法包括：接收从发无功/电压指令中获取的无功目标值；
采集所述光伏发电站的每个节点电压，其中，所述节点电压包括光伏单元节点电压和并网点电压；
在判断每个所述节点电压均位于预设的稳定裕度区间内的情况下，采用光伏逆变器输出无功功率，在所述光伏逆变器输出的无功功率达到限值的情况下，采用无功补偿装置输出无功功率；
在判断部分节点电压达到所述稳定裕度区间的上限Us+或下限Us‑且所述光伏发电站的无功功率未达到预设的容性/感性无功目标值的情况下，采用所述光伏逆变器输出无功功率，在所述光伏逆变器输出的无功功率达到限值情况下，再采用无功补偿装置输出无功功率；
在判断达到所述稳定裕度区间的上限Us+或下限Us‑的部分节点电压为所述光伏单元节点电压且所述光伏发电站的无功功率未达到预设的感性/容性无功目标值的情况下，同时采用所述无功补偿装置和所述光伏逆变器输出无功功率；
在判断达到所述稳定裕度区间的上限Us+或下限Us‑的部分节点电压为所述并网点电压且所述光伏发电站的无功功率未达到所述感性/容性无功目标值的情况下，采用所述无功补偿装置输出无功功率，在所述无功补偿装置输出的无功功率达到限值情况下，再采用所述光伏逆变器输出无功功率；
在判断所述节点电压超出所述稳定裕度区间且未超出预设的极值的情况下，采用所述无功补偿装置输出无功功率，在所述无功补偿装置输出的无功功率达到限值情况下，再采用所述光伏逆变器输出无功功率；
判断输出的所述光伏发电站输出的无功功率是否达到所述无功目标值且每个所述节点电压是否满足所述发无功/电压指令的要求；
在判断所述光伏发电站输出的无功功率未达到所述无功目标值和/或至少一个节点电压未满足所述发无功/电压指令的要求的情况下，再次执行采集每个所述节点电压并进行分析和输出无功功率的操作。
2.根据权利要求1所述的方法，其特征在于，所述方法进一步包括：计算用于控制所述光伏逆变器和/或所述无功补偿装置的无功优化模型；
根据所述无功优化模型控制所述光伏逆变器和/或无功补偿装置以输出无功功率。
3.根据权利要求2所述的方法，其特征在于，所述计算用于控制所述光伏逆变器和/或所述无功补偿装置的无功优化模型进一步包括：确定目标函数；
4.根据权利要求3所述的方法，其特征在于，所述确定目标函数包括：根据公式(1)计算用于指示所述光伏发电站的并网点实际输出的无功功率与从所述无功/电压指令中获取的与该并网点对应的目标无功功率的差值的第一目标子函数，f1＝Q1‑QM， (1)其中，Q1为所述并网点实际输出的无功功率，QM为从所述无功/电压指令中获取的与该并网点对应的目标无功功率；
根据公式(2)计算用于指示所述光伏单元节点电压偏差的第二目标子函数，其中，Ui为一个光伏单元i的节点电压的幅值，Uav为所述光伏单元节点电压的平均值，SPV为所有所述光伏单元的集合；
根据公式(3)计算用于指示所述光伏发电站有功功率损耗的第三目标子函数，其中，Ui、θi分别为一个节点电压i的电压幅值和相角,Uj为另一个节点电压j的电压幅值；Gij为一个所述光伏发电站的一个节点电压i和另一个节点电压j导纳的实部；θij是所述一个节点电压i和所述另一个节点电压j的相角差值；NB为所述光伏发电站的所述节点电压的数量的总和；SNB为所述光伏发电站的所有所述节点电压的集合；
根据公式(4)计算用于指示光伏发电站的无功电源动作的优先级的第四目标子函数，其中，a1、a2、a3、a4和a5分别为预设的调节因子，Us+为所述稳定裕度区间的上限，Qsvd为所述无功补偿装置输出的无功功率，Qsvgmin、Qsvgmax分别为所述无功补偿装置输出无功功率的下限值和上限值，Q1为所述并网点实际输出的无功功率，Qm为从所述无功/电压指令中获取的与该并网点对应的目标无功功率，Uim为光伏发电站的各个节点电压的最大值，Uav为所述平均值，Qsvg为所述无功补偿装置输出的无功功率，U1为所述并网点电压，Us‑为所述稳定裕度区间的下限；
根据公式(5)计算所述目标函数，minf(x)＝w1f1+w2f2+w3f3+w4f4， (5)其中，f1为所述第一目标子函数，f2为所述第二目标子函数，f3为所述第三目标子函数,f4为第四目标子函数，w1、w2、w3和w4分别为所述第一目标子函数、所述第二目标子函数、所述第三目标子函数和所述第四目标子函数的权重系数。
5.根据权利要求4所述的方法，其特征在于，所述确定约束条件包括：根据不等式(1)计算所述光伏逆变器的无功调节能力的第一约束，其中，QPV,i为一个光伏逆变器i输出的无功功率， 为一个光伏逆变器i输出的无功功率的下限值， 为所述光伏逆变器i输出的无功功率的上限值；
根据不等式(2)计算所述节点电压的第二约束，其中，Ui为一个节点电压i， 为一个节点电压i的下限值， 为一个节点电压i的上限值，SNB为所有所述节点电压的集合；
根据不等式(3)计算所述无功补偿装置的调节能力的第三约束，其中，Bc,i为一个无功补偿装置i输出的无功功率， 为一个无功补偿装置i输出的无功功率的下限值， 为一个无功补偿装置i输出的无功功率的上限值，SC为所有所述无功补偿装置的集合。
6.根据权利要求5所述的方法，其特征在于，所述确定约束条件包括：根据等式(1)确定所述约束条件的等式约束，其中：PGi、QGi分别为一个节点电压i所在节点的有功、无功注入；PDi、QDi分别为一个节点电压i所在节点的有功、无功负荷,Bij为所述一个节点电压i和另一个节点电压j导纳的虚部。
7.根据权利要求6所述的方法，其特征在于，所述采用粒子群算法计算所述无功优化模型包括：
在判断所述相邻的粒子陷入早熟的情况下，加入自适应微调算子；
采用所述粒子的位置、所述虚拟粒子的位置和所述下一个虚拟粒子的位置基于控制变量原则进行超限检测；
采用潮流计算对每个个体和全局变量进行更新以得出所述无功优化模型；
在判断所述无功优化模型满足所述约束条件的情况下，输出所述无功优化模型；
在判断所述无功优化模型未满足所述约束条件的情况下，再次返回执行利用外推变量调节计算虚拟粒子的位置的步骤，并执行所述粒子群算法。
8.一种用于对光伏发电站进行无功优化的系统，所述系统包括处理器，所述处理器用于执行如权利要求1至7任一所述的方法。
9.一种存储介质，其特征在于，所述存储介质存储有指令，所述指令用于被机器读取以使得所述机器执行如权利要求1至7任一所述的方法。
用于对光伏发电站进行无功优化的方法、系统及存储介质
技术领域
[0001] 本发明涉及光伏发电站调控技术领域，具体地涉及一种用于对光伏发电站进行无功优化的方法、系统及存储介质。
背景技术
[0002] 为应对能源危机，新能源发电站在近年来逐渐兴起。而在众多的新能源发电站中，光伏发电站占据最重要的席位之一。
[0003] 传统的光伏发电站的自动电压控制(AVC)目标仅针对并网点的无功/电压，控制策略通过优先调节逆变器或无功补偿装置的单一方式实现调节，且逆变器无功指令采取等无
功功率/等功率因数等方式计算。对于集电线路长度不均的中大型光伏发电站，传统的自动
电压控制策略会导致站内各节点电压偏差较大，甚至局部节点电压出现越限。
发明内容
[0004] 本发明实施方式的目的是提供一种用于对光伏发电站进行无功优化的方法、系统及存储介质，该方法、系统及存储介质能够解决现有技术中光伏发电站输出无功功率时节
点电压偏差较大的问题。
[0005] 为了实现上述目的，本发明实施方式的一方面提供一种用于对光伏发电站进行无功优化的方法，所述方法包括：
[0006] 接收从发无功/电压指令中获取的无功目标值；
[0007] 采集所述光伏发电站的每个节点电压，其中，所述节点电压包括光伏单元节点电压和并网点电压；
[0008] 对每个节点电压进行分析；
[0009] 在判断每个所述节点电压均位于预设的稳定裕度区间内的情况下，采用光伏逆变器输出无功功率，在所述光伏逆变器输出的无功功率达到最大值的情况下，采用无功补偿
装置输出无功功率；
[0010] 在判断部分节点电压达到所述稳定裕度区间的上限Us+或下限Us‑且所述光伏发电站的无功功率未达到预设的容性/感性无功目标值的情况下，采用所述光伏逆变器输出无
功功率，在所述光伏逆变器输出的无功功率达到限值情况下，再采用无功补偿装置输出无
功功率；
[0011] 在判断达到所述稳定裕度区间的上限Us+或下限Us‑的部分节点电压为所述光伏单元节点电压且所述光伏发电站的无功功率未达到预设的感性/容性无功目标值的情况下，
同时采用所述无功补偿装置和所述光伏逆变器输出无功功率；
[0012] 在判断达到所述稳定裕度区间的上限Us+或下限Us‑的部分节点电压为所述并网点电压且所述光伏发电站的无功功率未达到所述感性/容性无功目标值的情况下，采用所述
无功补偿装置输出无功功率，在所述无功补偿装置输出的无功功率达到限值情况下，再采
用所述光伏逆变器输出无功功率；
[0013] 在判断所述节点电压超出所述稳定裕度区间且未超出预设的极值的情况下，采用所述无功补偿装置输出无功功率，在所述无功补偿装置输出的无功功率达到限值情况下，
再采用所述光伏逆变器输出无功功率；
[0014] 判断输出的所述光伏发电站输出的无功功率是否达到所述无功目标值且每个所述节点电压是否满足所述发无功/电压指令的要求；
[0015] 在判断所述光伏发电站输出的无功功率未达到所述无功目标值和/或至少一个节点电压未满足所述发无功/电压指令的要求的情况下，再次执行采集每个所述节点电压并
进行分析和输出无功功率的操作。
[0016] 可选地，所述方法进一步包括：
[0017] 计算用于控制所述光伏逆变器和/或所述无功补偿装置的无功优化模型；
[0018] 根据所述无功优化模型控制所述光伏逆变器和/或无功补偿装置以输出无功功率。
[0019] 可选地，所述计算用于控制所述光伏逆变器和/或所述无功补偿装置的无功优化模型进一步包括：
[0020] 确定目标函数；
[0021] 确定约束条件；
[0022] 采用粒子群算法计算所述无功优化模型。
[0023] 可选地，所述确定目标函数包括：
[0024] 根据公式(1)计算用于指示所述光伏发电站的并网点实际输出的无功功率与从所述无功/电压指令中获取的与该并网点对应的目标无功功率的差值的第一目标子函数，
[0025] f1＝Q1‑QM， (1)
[0026] 其中，Q1为所述并网点实际输出的无功功率，QM为从所述无功/电压指令中获取的与该并网点对应的目标无功功率；
[0027] 根据公式(2)计算用于指示所述光伏单元节点电压偏差的第二目标子函数，
[0028]
[0029] 其中，Ui为一个光伏单元i的节点电压的幅值，Uav为所述光伏单元节点电压的平均值，SPV为所有所述光伏单元的集合；
[0030] 根据公式(3)计算用于指示所述光伏发电站有功功率损耗的第三目标子函数，
[0031]
[0032] 其中，Ui、θi分别为一个节点电压i的电压幅值和相角，Uj为另一个节点电压j的电压幅值；Gij为一个所述光伏发电站的一个节点电压i和另一个节点电压j导纳的实部；θij是
所述一个节点电压i和所述另一个节点电压j的相角差值；NB为所述光伏发电站的所述节点
电压的数量的总和；SNB为所述光伏发电站的所有所述节点电压的集合；
[0033] 根据公式(4)计算用于指示光伏发电站的无功电源动作的优先级的第四目标子函数，
[0034]
[0035] 其中，a1、a2、a3、a4和a5分别为预设的调节因子，Us+为所述稳定裕度区间的上限，Qsvd为所述无功补偿装置输出的无功功率，Qsvgmin、Qsvgmax分别为所述无功补偿装置输出无功
功率的下限值和上限值，Q1为所述并网点实际输出的无功功率，Qm为从所述无功/电压指令
中获取的与该并网点对应的目标无功功率，Uim为光伏发电站的各个节点电压的最大值，Uav
为所述平均值，Qsvg为所述无功补偿装置输出的无功功率，U1为所述并网点电压，Us‑为所述
稳定裕度区间的下限；
[0036] 根据公式(5)计算所述目标函数，
[0037] minf(x)＝w1f1+w2f2+w3f3+w4f4， (5)
[0038] 其中，f1为所述第一目标子函数，f2为所述第二目标子函数，f3为所述第三目标子函数，f4为第四目标子函数，w1、w2、w3和w4分别为所述第一目标子函数、所述第二目标子函
数、所述第三目标子函数和所述第四目标子函数的权重系数。
[0039] 可选地，所述确定约束条件包括：
[0040] 根据不等式(1)计算所述光伏逆变器的无功调节能力的第一约束，
[0041]
[0042] 其中，QPV，i为一个光伏逆变器i输出的无功功率， 为一个光伏逆变器i输出的无功功率的下限值， 为所述光伏逆变器i输出的无功功率的上限值；
[0043] 根据不等式(2)计算所述节点电压的第二约束，
[0044]
[0045] 其中，Ui为一个节点电压i， 为一个节点电压i的下限值， 为一个节点电压i的上限值，SNB为所有所述节点电压的集合；
[0046] 根据不等式(3)计算所述无功补偿装置的调节能力的第三约束，
[0047]
[0048] 其中，Bc，i为一个无功补偿装置i输出的无功功率， 为一个无功补偿装置i输出的无功功率的下限值， 为一个无功补偿装置i输出的无功功率的上限值，SC为所有
所述无功补偿装置的集合。
[0049] 可选地，所述确定约束条件包括：
[0050] 根据等式(1)确定所述约束条件的等式约束，
[0051]
[0052] 其中：PGi、QGi分别为一个节点电压i所在节点的有功、无功注入；PDi、QDi分别为一个节点电压i所在节点的有功、无功负荷，Bij为所述一个节点电压i和另一个节点电压j导纳的
虚部。
[0053] 可选地，所述采用粒子群算法计算所述无功优化模型包括：
[0054] 初始化粒子的位置和速度；
[0055] 利用外推变量调节计算虚拟粒子的位置；
[0056] 计算出下一虚拟粒子的位置；
[0057] 判断相邻的粒子是否陷入早熟；
[0058] 在判断所述相邻的粒子陷入早熟的情况下，加入自适应微调算子；
[0059] 采用所述粒子的位置、所述虚拟粒子的位置和所述下一个虚拟粒子的位置基于控制变量原则进行超限检测；
[0060] 根据检测结果对每个粒子的速度和位置进行更新；
[0061] 采用潮流计算对每个个体和全局变量进行更新以得出所述无功优化模型；
[0062] 判断所述无功优化模型是否满足所述约束条件；
[0063] 在判断所述无功优化模型满足所述约束条件的情况下，输出所述无功优化模型；
[0064] 在判断所述无功优化模型未满足所述约束条件的情况下，再次返回执行利用外推变量调节计算虚拟粒子的位置的步骤，并执行所述粒子群算法。
[0065] 本发明的另一方面还提供一种用于对光伏发电站进行无功优化的系统，该系统包括控制器，所述控制器用于执行上述所述的方法。
[0066] 本发明的再一方面还提供一种存储介质，所述存储介质存储有指令，所述指令用于被机器读取以使得所述机器执行上述所述的方法。
[0067] 通过过上述技术方案，本发明实施方式提供的用于对光伏发电站进行无功优化的方法、系统及存储介质通过采集光伏发电站的节点电压并加以分析，根据分析的结果控制
光伏发电站的光伏逆变器和无功补偿装置输出无功功率，解决了现有技术中依赖于单一的
无功输出设备而使得节点电压偏差过大的问题，提高了光伏发电站无功优化的效率。
[0068] 本发明实施方式的其它特征和优点将在随后的具体实施方式部分予以详细说明。
附图说明
[0069] 附图是用来提供对本发明实施方式的进一步理解，并且构成说明书的一部分，与下面的具体实施方式一起用于解释本发明实施方式，但并不构成对本发明实施方式的限
制。在附图中：
[0070] 图1是根据本发明的一个实施方式的用于对光伏发电站进行无功优化的方法的流程图；
[0071] 图2是根据本发明的一个实施方式的生成无功优化模型的流程图；
[0072] 图3是根据本发明的一个实施方式的粒子群算法的流程图；
[0073] 图4是根据本发明的一个实施方式的光伏发电站的模型示意图；以及
[0074] 图5所示根据本发明的一个实施方式的不同方式对光伏发电站进行无功优化的结果的曲线图。
具体实施方式
[0075] 以下结合附图对本发明实施方式的具体实施方式进行详细说明。应当理解的是，此处所描述的具体实施方式仅用于说明和解释本发明实施方式，并不用于限制本发明实施
方式。
[0076] 如图1所示是根据本发明的一个实施方式的用于对光伏发电站进行无功优化的方法的流程图。在图1中，该方法可以包括：
[0077] 在步骤S10中，获取调度任务，从调度任务中获取所需发无功/电压指令并计算无功目标值。
[0078] 在步骤S20中，采集光伏发电站的每个节点电压。其中，节点电压包括光伏单元节点电压和并网点电压。另外，采集光伏发电站的节点电压的方式应当是本领域人员所知的，
故此处不详细赘述。
[0079] 在接收到每个节点电压后，可以对每个节点电压进行分析。
[0080] 在步骤S30中，在判断每个节点电压(的幅值)均位于预设的稳定裕度区间内的情况下，采用光伏逆变器输出无功功率，并且在光伏逆变器输出的无功功率达到限值(输出的
无功功率的最大值)的情况下，采用无功补偿装置输出无功功率。
[0081] 在步骤S40中，在判断部分节点电压达到稳定裕度区间的上限Us+或下限Us‑且光伏发电站的无功功率未达到预设的容性/感性无功目标值的情况下，采用光伏逆变器输出无
功功率，在光伏逆变器输出的无功功率达到限值情况下，再采用无功补偿装置输出无功功
率。
[0082] 在判断达到稳定裕度区间的上限Us+或下限Us‑的部分节点电压为光伏单元节点电压且光伏发电站的无功功率未达到预设的感性/容性无功目标值的情况下，同时采用无功
补偿装置和光伏逆变器输出无功功率。
[0083] 在判断达到稳定裕度区间的上限Us+或下限Us‑的部分节点电压为并网点电压且光伏发电站的无功功率未达到感性/容性无功目标值的情况下，采用无功补偿装置输出无功
功率，在无功补偿装置输出的无功功率达到限值情况下，再采用光伏逆变器输出无功功率。
[0084] 在步骤S50中，在判断节点电压超出稳定裕度区间且未超出预设的极值(超出极值的情况下很可能会发生脱网危险事故)的情况下，采用无功补偿装置输出无功功率，在无功
补偿装置输出的无功功率达到限值情况下，再采用光伏逆变器输出无功功率。
[0085] 在步骤S60中，判断输出的光伏发电站(光伏逆变器和/或无功补偿装置)输出的无功功率是否达到无功目标值且每个节点电压是否满足发无功/电压指令的要求。
[0086] 在判断光伏发电站输出的无功功率未达到无功目标值和/或至少一个节点电压未满足发无功/电压指令的要求的情况下，返回执行步骤S20，并根据继续执行该方法的相应
流程。
[0087] 由于光伏逆变器和/或无功补偿装置发无功(输出无功功率)需要根据光伏发电站每个节点多个外部(例如线路的长度)以及内部条件(例如节点电压的幅值)来确定，那么在
本发明的一个实施方式中，该方法可以进一步包括：计算用于控制光伏逆变器和/或无功补
偿装置的无功优化模型，在进一步根据该计算出的无功优化模型控制光伏逆变器和/或无
功补偿装置发无功。
[0088] 在本发明的一个示例中，如图2所示，计算用于控制光伏逆变器和/或无功补偿装置的无功优化模型的方式可以是例如包括以下步骤：
[0089] 在步骤S21中，确定目标函数。在该示例中，确定该目标函数的具体步骤可以是例如：
[0090] 根据公式(1)计算用于指示光伏发电站的并网点实际输出的无功功率与从无功/电压指令中获取的与该并网点对应的目标无功功率的差值的第一目标子函数，
[0091] f1＝Q1‑QM， (1)
[0092] 其中，Q1为并网点实际输出的无功功率，QM为从无功/电压指令中获取的与该并网点对应的目标无功功率；
[0093] 根据公式(2)计算用于指示光伏单元节点电压偏差的第二目标子函数，
[0094]
[0095] 其中，Ui为一个光伏单元i的节点电压的幅值，Uav为光伏单元节点电压的平均值，SPV为所有光伏单元的集合；
[0096] 根据公式(3)计算用于指示光伏发电站有功功率损耗的第三目标子函数，
[0097]
[0098] 其中，Ui、θi分别为一个节点电压i的电压幅值和相角；Uj为另一个节点电压j的电压幅值；Gij为一个光伏发电站的一个节点电压i和另一个节点电压j导纳的实部；θij是一个
节点电压i和另一个节点电压j的相角差值；NB为光伏发电站的节点电压的数量的总和；SNB
为光伏发电站的所有节点电压的集合；
[0099] 根据公式(4)计算用于指示光伏发电站的无功电源动作的优先级的第四目标子函数，
[0100]
[0101] 其中，a1、a2、a3、a4和a5分别为预设的调节因子，Us+为预设的稳定裕度区间的上限，Qsvd为无功补偿装置输出的无功功率，Qsvgmin、Qsvgmax分别为无功补偿装置输出无功功率的下
限值和上限值，Q1为并网点实际输出的无功功率，Qm为从无功/电压指令中获取的与该并网
点对应的目标无功功率，Uim为光伏发电站的各个节点电压的最大值，Uav为平均值，Qsvg为无
功补偿装置输出的无功功率，U1为并网点电压，Us‑为预设的稳定裕度区间的下限；
[0102] 根据公式(5)计算目标函数，
[0103] minf(x)＝w1f1+w2f2+w3f3+w4f4， (5)
[0104] 其中，f1为第一目标子函数，f2为第二目标子函数，f3为第三目标子函数，f4为第四目标子函数，w1、w2、w3和w4分别为第一目标子函数、第二目标子函数、第三目标子函数和第
四目标子函数的权重系数。
[0105] 在步骤S22中，确定约束条件。在该示例中，该约束条件可以包括等式约束(条件)和不等式约束(条件)。那么，确定该等式约束的具体步骤可以是例如：
[0106] 根据等式(1)确定约束条件的等式约束，
[0107]
[0108] 其中：PGi、QGi分别为一个节点电压i所在节点的有功(功率)、无功(功率)注入；PDi、QDi分别为一个节点电压i所在节点的有功、无功负荷，Bij为一个节点电压i和另一个节点电
压j导纳的虚部。
[0109] 确定该不等式约束条件的具体步骤可以是例如：
[0110] 根据不等式(1)计算光伏逆变器的无功调节能力的第一约束，
[0111]
[0112] 其中，QPV，i为一个光伏逆变器i输出的无功功率， 为一个光伏逆变器i输出的无功功率的下限值， 为光伏逆变器i输出的无功功率的上限值；
[0113] 根据不等式(2)计算节点电压的第二约束，
[0114]
[0115] 其中，Ui为一个节点电压i， 为一个节点电压i的下限值， 为一个节点电压i的上限值，SNB为所有节点电压的集合；
[0116] 根据不等式(3)计算无功补偿装置的调节能力的第三约束，
[0117]
[0118] 其中，Bc，i为一个无功补偿装置i输出的无功功率， 为一个无功补偿装置i输出的无功功率的下限值， 为一个无功补偿装置i输出的无功功率的上限值，SC为所有
无功补偿装置的集合。
[0119] 在步骤S23中，采用粒子群算法计算无功优化模型。在该示例中，该粒子群算法的计算步骤可以是例如图3中示出的步骤。在图3中，该粒子群算法可以包括：
[0120] 在步骤S31中，初始化粒子的位置和速度。在该示例中，可以是首先初始化一个规模为N的粒子群，设定每个粒子的各个边界条件，在满足控制变量的约束条件下随机赋予粒
子群中的各粒子i初始位置 和速度
[0121] 在步骤S32中，利用外推变量调节计算虚拟粒子的位置。在该示例中，可以将每个粒子代入目标函数计算出相应的适应值，根据搜索到的最优位置，通过基本粒子群算法中
位置更新式 来更新自己的飞行速度和位置。其中， 为计算
出的粒子i的位置， 为计算出的粒子i的速度。通过粒子i的适应值差异来引导外推方
向，并进一步结合随机数算法，使用公式 在推断出的
粒子i附近计算出虚拟粒子的位置
[0122] 在步骤S33中，计算出下一虚拟粒子的位置。在该示例中，可以是根据步骤S32中的原理进一步推断出下一个虚拟粒子的位置，即可以采用公式(6)表示下一个虚拟粒子的位
置，
[0123]
[0124] 其中， 为下一个虚拟粒子的位置，rand[]为随机数，k为预设的参数因子。
[0125] 在步骤S34中，判断相邻的粒子是否陷入早熟。在该示例中，可以是将粒子i、计算出的粒子i、虚拟粒子、下一个虚拟粒子的多个分量(位置和/或速度)进行比对，判断相邻的
粒子是否出现(位置和/或速度)相近或者相同的情况。
[0126] 在步骤S35中，在判断相邻的粒子陷入早熟的情况下，加入自适应微调算子。在该示例中，可以是在判断相邻的粒子出现(位置和/或速度)相近或相同的情况下，在公式(6)
中加入自适应微调算子ε，即可以采用公式(7)表示该下一虚拟粒子的位置，
[0127]
[0128] (7)
[0129] 其中，ε为加入的自适应微调算子。另外，在判断相邻的粒子没有陷入早熟的情况下，也可以直接执行步骤S36。
[0130] 在步骤S36中，采用粒子的位置、虚拟粒子的位置和下一个虚拟粒子的位置基于控制变量原则进行超限检测。在该示例中，可以是利用预设的目标函数(例如步骤S20中示出
的目标函数)分别计算每个粒子对应的适应值，从而实现对每个粒子的适应度的重新评估，
在保留最优解的情况下，更新每个个体的各个变量(位置和/或速度)，从而实现对整体粒子
群的更新。
[0131] 在步骤S37中，根据检测结果对每个粒子的速度和位置进行更新。该步骤在步骤S36中已经详述，此处不再赘述。
[0132] 在步骤S38中，采用潮流计算对每个个体和全局变量进行更新以得出无功优化模型。该步骤在步骤S36中已经详述，此处不再赘述。
[0133] 在步骤S39中，判断无功优化模型是否满足约束条件。在该示例中，可以是将更新后的粒子群(即无功优化模型)与预设的约束条件(例如步骤S21中示出的约束条件)进行比
对。
[0134] 在步骤S41中，在判断无功优化模型满足约束条件的情况下，输出无功优化模型。在该示例中，在更新后的粒子群满足约束条件的情况下，输出该粒子群(无功优化模型)。
[0135] 在判断无功优化模型未满足约束条件的情况下，再次返回执行步骤S32直到该无功优化模型满足预设的约束条件。
[0136] 本发明的另一方面还提供一种用于对光伏发电站进行无功优化的系统，系统包括处理器。该处理器可以用于执行上述所述方法。
[0137] 在该实施方式中，该处理器可以是例如通用处理器、专用处理器、常规处理器、数字信号处理器(DSP)、多个微处理器、与DSP核心相关联的一个或多个微处理器、控制器、微
控制器、专用集成电路(ASIC)、现场可编程门阵列(FPGA)电路、任何其它类型的集成电路
(IC)、状态机、系统级芯片(SOC)等。
[0138] 本发明的再一方面还提供一种存储介质。该存储介质可以存储有指令，该指令可以用于被机器读取以使得机器执行上述所述的方法。
[0139] 在该实施方式中，该存储介质可以是例如U盘、移动硬盘、只读存储器(ROM，Read‑Only Memory)、随机存取存储器(RAM，Random Access Memory)、磁碟或者光盘等各种可以
存储程序代码的介质。
[0140] 在本发明的一个示例中，以如图4所示出的100MW光伏发电站作为测试对象，采用包括但不限于Matlab7.0的算法软件来编制简化模型、粒子群算法和潮流计算算法。其中，
为了便于算法的简便性，可以将粒子群N的数量设定为20，最大迭代次数Tmax设定为100。
[0141] 在图4中，光伏发电站的装机容量为100MW，共10条集电线路，每回集电线路串联10组箱式逆变器(光伏逆变器)，每组箱式逆变器的容量为1MW(功率因素在‑0.95～0.95之间
可调)，通过0.4kV/35kV箱式变压器、集电线路及5kV/110kV升压变压器并入电网，35kV母线
处安装有1台±20MVar的静止无功发生器14(Static Var Generator，SVG)(即无功补偿装
置)。为便于仿真验证，将每条集电线路(编号为4至13和15至24)上的10台箱式逆变器、10台
箱式变压器等效为1组光伏单元，#1～#10集电线路长度分别约为2.2、2.4、2.6、2.8、3.0、
3.2、3.4、3.6、3.8、4.0km。
[0142] 在图4中预设公共连接点(编号为1)电压为1.0pu(per unit)，从发无功/电压指令中获取的并网点(编号为2)的无功目标值未‑15Mvar，100个光伏单元出力(输出无功功率)
均为8.46MW，预设的稳定裕度区间为[0.968，‑1.602]Un。
[0143] 如图5所示是采用三种不同的方式对图4中示出的光伏发电站进行无功优化的结果。其中，方式1为先采用光伏逆变器输出无功功率，在光伏逆变器的输出达到限值的情况
下，再使用无功补偿装置输出无功功率；方式2为先使用无功补偿装置输出无功功率，在无
功补偿装置的输出达到限值的情况下，再使用光伏逆变器输出无功功率；方式3为采用本发
明提供的用于对光伏发电站进行无功优化的方法对该光伏发电站的无功功率输出进行控
制。
[0144] 在图5中，可以发现，采用方式1的无功优化模型，光伏发电站的各个节点的节点电压相对较高；采用方式2的无功优化模型，光伏发电站的各个节点的节点电压相对较低；而
采用本发明提供的方法的无功优化模型则可以使得每个节点的节点电压的幅值适中。因
此，本发明提供的方法的无功优化效率优于传统的方式1和方式2。
[0145] 另外，表(1)所示为三种方式调控下光伏发电站的各个参数的数值。
[0146] 表(1)
[0147]
[0148] 在表(1)中，每个数据均为标幺值，在计算时节点电压的基准值可以选取设备的额定电压，功率的基准值可以为100MVA。
[0149] 从表(1)中也可以看出，方式1、2和方式3相比，本发明提供的方法使得光伏发电站能够在满足无功目标值的前提下，将有功损耗降低至0.0177MW，并同时能够使得各个节点
电压保持在稳定裕度区间范围内，最后还能够使得无功补偿装置具备一定的无功裕度。因
此，可以看出本发明提供的方法的无功优化效率更高，经济性、环保性也更好。
[0150] 通过过上述技术方案，本发明实施方式提供的用于对光伏发电站进行无功优化的方法、系统及存储介质通过采集光伏发电站的节点电压并加以分析，根据分析的结果控制
光伏发电站的光伏逆变器和无功补偿装置输出无功功率，解决了现有技术中依赖于单一的
无功输出设备而使得节点电压偏差过大的问题，提高了光伏发电站无功优化的效率。
[0151] 以上结合附图详细描述了本发明例的可选实施方式，但是，本发明实施方式并不限于上述实施方式中的具体细节，在本发明实施方式的技术构思范围内，可以对本发明实
施方式的技术方案进行多种简单变型，这些简单变型均属于本发明实施方式的保护范围。
[0152] 另外需要说明的是，在上述具体实施方式中所描述的各个具体技术特征，在不矛盾的情况下，可以通过任何合适的方式进行组合。为了避免不必要的重复，本发明实施方式
对各种可能的组合方式不再另行说明。
[0153] 本领域技术人员可以理解实现上述实施方式方法中的全部或部分步骤是可以通过程序来指令相关的硬件来完成，该程序存储在一个存储介质中，包括若干指令用以使得
一个(可以是单片机，芯片等)或处理器(processor)执行本申请各个实施方式所述方法的
全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(ROM，Read‑Only
Memory)、随机存取存储器(RAM，Random Access Memory)、磁碟或者光盘等各种可以存储程
序代码的介质。
[0154] 此外，本发明实施方式的各种不同的实施方式之间也可以进行任意组合，只要其不违背本发明实施方式的思想，其同样应当视为本发明实施方式所公开的内容。
