基于风速预测的多时段电气互联系统最优潮流计算方法
【专利摘要】本发明公布了一种基于日前风速预测的多时段电气互联系统最优潮流计算方法，适用于电力系统优化控制领域。本发明首先提出了基于变分模态分解和高斯过程回归的风速预测方法，并由此得到了日前预测风速的概率分布曲线。建立了电?气互联系统多时段最优潮流模型，并以总运行费用最小为目标，模型中计及了电力系统和天然气系统的运行约束。采用惩罚成本和备用成本分别描述风电高估和低估造成的影响。算例表明，电力系统与天然气系统存在相互制约关系，综合优化不仅有利于得到全局最优解，更保证了系统的安全可靠性。此外，风电的惩罚成本和备用成本对调度方案有重要影响，为新能源接入背景下系统优化运行提供参考，为调度人员提供决策支持。
【专利说明】
基于风速预测的多时段电气互联系统最优潮流计算方法
技术领域
[0001] 本发明属于电力系统运行分析与控制技术领域，设及一种基于日前风速预测的多 时段电气互联系统最优潮流计算方法。 技术背景
[0002] 近年来，随着经济技术的日益发展，能源枯竭和环境问题加快了天然气产业发展 的步伐。天然气主要优势在于=个方面:相对较低的投资成本、低溫使气体排放W及灵活快 速的响应特性。因而，在可预见的未来，天然气有望在电力系统中大规模应用。美国能源信 息署于2015年颁布了《清洁电力计划》，预计到2035年美国60%的发电机组为燃气机组。燃 气机组建立了电力网与天然气网络的禪合关系，运两个网络之间的强禪合关系要求从去整 体的角度去分析。然而，传统最优潮流计算都是W电力系统为导向，忽略了天然气系统的约 束，由此得到的解可能丧失最优性，甚至可能会违背天然气系统的安全性。因而有必要考虑 考虑电-气互联系统的最优潮流计算。此外，电力系统是一个动态的系统，其运行状态应当 实时得到监控，应当考虑多时段的优化分析。基于上述背景，本发明提出了多时段电-气互 耳关系统最优潮流计算（multi-time combined 邑as and electric optimal power flow, MGEOPF)O
[0003] 另一方面，燃气轮机组的灵活性和快速响应特性可缓解风电等可再生能源的波动 性，运为可再生能源的大规模接入提供了契机。准确的风速预测有利于制定合理的调度计 划，减小调度误差W及旋转备用的配置。基于此，本发明提出了一种日前风速预测方法，W 提高风速预测的准确性。
【发明内容】
[0004] 为克服现有技术上的不足，本发明的目的在于对互联的电力系统与天然气系统进 行联合优化调度，W互联系统的总运行成本最低为优化控制目标，兼顾电力系统与天然气 系统的安全性，同时考虑了大规模风电场的接入，提出了一种日前风速预测的新方法，即基 于日前风速预测的多时段电气互联系统最优潮流计算方法W提高预测精度，减小调度误 差。
[0005] 技术方案:本发明提供一种基于日前风速预测的多时段电气互联系统最优潮流 计算方法，包括W下步骤：
[0006] 步骤1:利用变分模态分解(variational mode decomposition,VMD)将历史风速 序列分解成若干个子序列，W降低原始序列的非平稳性。利用偏相关函数确定预测的输入 变量。利用高斯过程回归（Gaussian process regression,GPR)对子序列进行概率预测，最 后将所有序列预测结果叠加，得到最终的风速预测结果；
[0007] 步骤2:对天然气系统和电力系统的元件进行稳态建模，通过燃气轮机建立两个网 络之间的禪合。W互联系统的综合运行成本最小为目标函数，建立具有确定形式的多时段 电-气互耳关系统最优潮流（multi-time combined gas and electric system optimal power flow,MGEOPF)模型；
[0008] 步骤3:将风速预测误差带来的影响加入目标函数中，引入风电惩罚成本表示因低 估风电场出力造成大量的风能资源浪费，引入风电备用成本表示因高估风电场出力造成的 大量投入旋转备用。采用概率积分的形式表示风电的惩罚成本和备用成本；
[0009] 步骤4:对所建立的模型，采用合适的求解算法进行求解，并对实际风电场进行日 前风速预测实例分析，采用IE邸标准算例对MGEOPF模型进行分析。
[0010] 进一步，步骤1包括W下步骤：
[0011] 步骤IOl=VMD算法是变分问题的求解过程，该算法可分为变分问题的构造和求解。 VMD算法的目的在于将输入信号f分解成有限个子信号（模式），该算法中本征模态函数 (IM巧被重新定义为一个调幅-调频信号，表达为：
[0012] Uk(t)=Ak(t)cos(4)k(t)) (1)式中：Ak(t)为Uk(t)的瞬时幅值。《k(t)为Uk(t)的 瞬时频率，脚―的=妍-的=d脚-/&。Ak( t)及O k( t)相对于相位妍(0来说变化是缓慢的。
[0013] 假定将原始信号分解为K个IMF分量，则对应的约束变分模型表达如下：
[0014]
(2)
[001引式中：{uk(t)} = {山，...，uk}代表分解得到的 K个 IMF 分量；{wk(t)} = {wi，...， W k}表不各分量的频率中屯、。
[0016] 上述约束变分问题可W采用多种求解算法，此处引入二次惩罚项和拉格朗日乘 子入W及惩罚参数a。得到增广Lagrange函数如下：
[0017]
(3)
[0018] 利用交替方向乘子算法求取上述增广Lagrange函数的鞍点，即为约束变分模型的 最优解，从而将原始信号分解为K个窄带IMF分量。
[0019] 步骤102:在利用GH?进行预测前，应当首先确定训练样本各分量的输入变量，可采 用统计工具，即偏相关函数W及偏相关函数图来确定预测的输入变量。
[0020] 具体地，假设Xi为输出变量，如果在滞后k阶处偏相关系数超过95%置信区间，即 在区间，则将xi-k作为其中一个输入变量。值得注意的是，有些情况下 所有的偏相关系数都在95%区间外，此时取Xi-I作为输入变量。
[0021] 步骤103:对于给定的数据集看
舆中：输入数据矩阵XiGRd,输出数 据矩阵yiGR。在给定数据D的有限集合中，f(xW)、f(xW)、…、f(xW)可构成随机变量的一 个集合，且具有联合高斯分布，高斯过程的全部统计特征由均值函数m(x)和协方差函数k (x，x'）组成，即
[0022] f(x)~GP(m(x)，k(x，x'）） （4)
[0023] 将含噪声考虑到观测目标值y中，可建立高斯过程回归问题的一般模型，即
[0024] y = f(x)+e (5)式中：e为独立的高斯白噪声，符合高斯分布，均值为0,方差为62， 即可记作e~N(0,62)。由于噪声e为独立于f(x)的白噪声，当f(x)服从高斯分布，则y同样服 从高斯分布，则其有限观测值联合分布的集合可形成一个高斯过程，即
[0025]
(6)
[00%] 式中：m(x)为均值函数；Sij为Kronecker delta函数，当i = j时，函数Sij = 1;协方 差函数若Pi巧随巧式亲巧，良口
[0027]
(7)
[0028] 式中：I表示NXN的单位矩阵，CU,X)表示NXN的协方差矩阵，KU,X)表示NXN的 核矩阵，称为Gram矩阵，其元素Kij = k(xi，Xj)。
[0029] 根据贝叶斯原理，GP在给定数据D的集合内建立先验函数，在n*个给定测试数据集 於
斗专变为后验分布，则训练数据向量X*的输出向量f和测试数据的输出 向量7之间的联合高斯分布
[0030]
(群
[0031] 由此可W得到主要的GP回归方程，即
[0032]
[0033] (10)
[0034] (11)
[003引其中，预测均值向量哀就是GP回归模型的输出，也即输出向量f*的预测值。
[0036] 进一步，步骤2中模型为：
[0037] 步骤201:天然气系统由气井、管道、加压站、储气罐、阀口等元件组成。天然气从远 距离的气田通过管道传输，再经过配气网络输送至终端用户。对于绝热管道，其首末端节点 分别为m和n，则稳态流量为：
[00：3 引 （12)
[0039] (13)
[0040] 式中：Cmn为管道效率;JImt, JInt为t时段节点压力。
[0041] 由于摩擦力的存在，天然气在管道中流动会有部分能量损失，加压站用于补偿运 部分能量的损耗。调节压力需要消耗额外的能量，运部分能量与流量和加压比成正比。
[0042] 对于加压站k，其流入和流出节点分别为P和q，则其能量消耗为：
[0043] (14)
[0044] (")
[0045] 式中：化t为加压站消耗的功率;故为加压常数;孔为压缩因子;JTkt为等效消耗的流 量;a，e，T为效率转换常数。
[0046] 天然气系统的任一节点应满足流量平衡方程，表示为：
[0047]
(化）
[004引式中:Wgt为气源点注入量;Wdt为气负荷量;Pg,it为第i台发电机的有功出力;祝Pg对 为燃气轮机消耗的天然气量。
[0049] 步骤202: W电-气互联系统联合运行总成本最低为目标函数的MGEOPF模型如下：
[0050] (1)目标函数
[0化1]
(17)
[0052] 式中：T，Nt，Nw，Q S为时段，燃煤机组，风电场，气源点集合;Pw, it为风电场计划出 力;Pav,it为风电场可用风功率。
[0053] 其中，式(17)第一项表示发电成本，通常W发电量的二次函数计算：
[0化4]
(18)
[0055] 式中：ai，bi，ci为发电成本系数。
[0056] 式（17)的第二项为风电成本，该项取决于风电场的所有者，第=项表示天然气成 本，可表示为天然气流量的线性函数：
[0化7]
(19)
[0058] (2)电力《统约巧
[0059] 电力系统多时段约束包括功率平衡约束、发电机出力约束、发电机爬坡约束、节点 电压约束、线路巧率约束，表达如下：
[0060]
[0061 ]式中：Qg,it为第i台发电机无功
出力；Pd,it,Qd,it为节点i电负荷;Vit, 0it为节点i电 压幅值和相角，Qijt =目it-目jt;Gij ,Bij为节点巧日节点j之间的电导和电纳;RUi,RD功发电机 爬坡率;Pi, t为线路1的功率;上标'min '和'max '表示最大和最小值.
[0062] (3)天然气系统约束
[0063] 天然气系统约束包括节点流量平衡方程、气源注入量约束、节点压力约束W及加 压站加压比约束，其表达如下：
[0064]
(21 I
[0065] 式中：K2i，K。，Koi为燃气轮机转换常数；Wgm…，Wgmax为最小和最大气源供应量川…， nmmax为节点压力最小和最大值;^min ,^max为最小和最大加压比。
[0066] 进一步，步骤3包括W下步骤：
[0067] 风能天然具有随机性和波动性，因此，应当在式（17)目标函数中计及风电不确定 性的影响。引入风电惩罚成本用W描述因低估风电带来的风能资源浪费惩罚，表示为：
[006引
（巧
[0069] 式中:k?为惩罚成本系数(低估）;Pr为额定风功率。
[0070] 引入风电备用成本用W描述因高估风电带来的旋转备用配置成本，表示为：
[0071]
(23)
[0072] 巧甲：k?刃爸用成本《数（局巧），av是available的缩与，意思是可用的，此处只是 一种表不方式。
[0073] 进一步，步骤4包括W下步骤：
[0074] 本发明所建立的MGEOPF模型为复杂的非线性规划模型，并且由于式(22)和式(23) 中积分项的引入，难W采用解析算法求解。可W将式(22)和式(23)积分项转化为数值积分， 进而采用非线性规划求解器进行求解。
[0075] 有益效果:本发明采用改进的日前风速预测方法提高预测的精度，并W此为基础， 进行日前调度计划的制定，W天然气和电力系统日运行成本最小为目标函数，兼顾了两个 系统的安全性。模型中同时计及了因风电预测误差带来的惩罚成本和备用成本。本发明适 用于天然气资源和新能源大规模使用的背景下的电力系统分析，为调度人员提供有力的决 策依据。其具有如下优点和技术效果：
[0076] 1)本发明通过求解多时段电-气互联系统最优潮流模型得到电力网与天然气网联 合优化调度方案。实际上，电力网与天然气网存在相互制约关系，对电力网与天然气网的独 立优化将导致优化结果过于乐观，联合优化能为电力系统调度人员W及天然气网调度人员 的正确决策提供依据，确保系统安全运行；
[0077] 2)本发明提出的风速预测方法采用VMD将原始序列分解为若干个子模态，降低了 原始序列的非平稳性，采用GPR方法进行风速预测得到了风速的概率预测结果，提高了预测 的准确性；
[0078] 3)计及了风电出力的不确定性，用惩罚成本表示风电低估成本，用备用成本表示 风电高估成本，所得结果更具有实际意义。
【附图说明】：
[0079] 附图1本发明的VMD分解结果；
[0080] 附图2本发明的风速预测结果图；
[0081 ]附图3本发明的日有功曲线图；
[0082] 附图4本发明的发电机有功出力曲线。
【具体实施方式】
[0083] W下结合附图和实例对本发明的实施作进一步说明，但本发明的实施和包含不限 于此。
[0084] 本实施例的一种基于日前风速预测的多时段电气互联系统最优潮流计算方法，其 首先根据美国可再生能源实验室提供的风速数据进行风速预测，选4号风电场为预测对象， 该区域平均风速为7.668m/s。选用2006年一月份风速数据作为观测样本，每10分钟记录1个 数据点，总共有720个数据点。将每小时内的6个数据取平均值作为该小时的平均风速，便可 得到120个观测点，将前96个点作为训练样本，后24个作为预测样本。附图巧由VMD算法分 解得到的风速子序列，附图2为风速预测结果。通过计算，风速预测的平均平方误差为 0.507，平均绝对误差为0.594。可见本发明所提预测方法具有较高的预测精度。
[0085] 采用修改的IE趾39节点系统和比利时20节点天然气系统进行MGEOPF算例分析。 Ma化ower4.1提供的IE邸39节点系统共有S个分区，10台发电机组，总装机容量7367MW，总 有功负荷6254.23MW，输电线路34条，变压器12台。日有功负荷曲线如附图3。假定区域1的 节点9接入600MW的风电场。风电场切入风速3.5m/s，额定风速15m/s，切出风速25m/s。
[0086] 比利时20节点天然气网有20个节点，21条输气管道，2个加压站，6个气源点，假定 加压站均为电力驱动，气负荷为46.298Mm^d。
[0087] 假定风电高估和低估成本系数均为2000$/MW，则W24小时为调度周期，总优化成 本为6465.403k$，机组出力曲线如附图4。可W发现机组31和机组34都为燃气机组，机组31 始终是满发状态，而机组34出力很小，运是因为机组34在天然气网络中的接入点压力较为 敏感，限制了燃气轮机的出力。
[0088] 选择负荷高峰时段19:00进行分析，该时段各机组出力如附表1，天然气网络的节 点压力W及气源点的供气量如附表2。
[0089] 附表1 19:00发电机出力结果
[0090]
[0091]
[0092]
[0093]
[0094] 可W看出，燃气轮机接入点节点20压力已达到下限，运与上述分析是一致的。该算 例充分说明了电力网络和天然气网络整体分析的必要性，若忽略天然气网络的安全约束， 极有可能导致所得解丧失最优性甚至丧失可行性。
[00%]本实施例公布了一种基于日前风速预测的多时段电气互联系统最优潮流计算方 法，适用于电力系统优化控制领域，其提出了基于变分模态分解和高斯过程回归的风速预 测方法，并由此得到了日前预测风速的概率分布曲线。建立了电-气互联系统多时段最优潮 流模型，并W总运行费用最小为目标，模型中计及了电力系统和天然气系统的运行约束。针 对风电出力的随机性和不确定性，采用惩罚成本和备用成本分别描述风电高估和低估造成 的影响。算例表明，电力系统与天然气系统存在相互制约关系，综合优化不仅有利于得到全 局最优解，更保证了系统的安全可靠性。此外，风电的惩罚成本和备用成本对调度方案有重 要影响。本发明为新能源接入背景下系统优化运行提供参考，为调度人员提供决策支持。
[0096] 本发明采用改进的日前风速预测方法提高预测的精度，并W此为基础，进行日前 调度计划的制定，W天然气和电力系统日运行成本最小为目标函数，兼顾了两个系统的安 全性。模型中同时计及了因风电预测误差带来的惩罚成本和备用成本。本发明适用于天然 气资源和新能源大规模使用的背景下的电力系统分析，为调度人员提供有力的决策依据。
[0097] 此外，本发明通过求解多时段电-气互联系统最优潮流模型得到电力网与天然气 网联合优化调度方案。实际上，电力网与天然气网存在相互制约关系，对电力网与天然气网 的独立优化将导致优化结果过于乐观，联合优化能为电力系统调度人员W及天然气网调度 人员的正确决策提供依据，确保系统安全运行;并提出的风速预测方法采用VMD将原始序列 分解为若干个子模态，降低了原始序列的非平稳性，采用GPR方法进行风速预测得到了风速 的概率预测结果，提高了预测的准确性；
[0098] 本发明计及了风电出力的不确定性，用惩罚成本表示风电低估成本，用备用成本 表示风电高估成本，所得结果更具有实际意义。
[0099] W上显示和描述了本发明的基本原理和主要特征和本发明的优点。本行业的技术 人员应该了解，本发明不受上述实施例的限制，上述实施例和说明书中描述的只是说明本 发明的原理，在不脱离本发明精神和范围的前提下，本发明还会有各种变化和改进，运些变 化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其 等效物界定。
【主权项】
1. 一种基于日前风速预测的多时段电气互联系统最优潮流计算方法，其特征在于包括 W下步骤： 步骤1:利用变分模态分解将历史风速序列分解成若干个子序列，并利用偏相关函数确 定预测的输入变量;然后利用高斯过程回归对子序列进行概率预测，最后将所有序列预测 结果叠加，得到最终的风速预测结果； 步骤2:对天然气系统和电力系统的元件进行稳态建模，通过燃气轮机建立两个网络之 间的禪合；W互联系统的综合运行成本最小为目标函数，建立具有确定形式的多时段电气 互联系统最优潮流模型； 步骤3:将风速预测误差带来的影响加入所述目标函数中，引入风电惩罚成本表示因低 估风电场出力造成的风能资源浪费，引入风电备用成本表示因高估风电场出力造成的投入 旋转备用;采用概率积分的形式表示风电的惩罚成本和备用成本； 步骤4:对所建立的多时段电气互联系统最优潮流模型，采用求解算法进行求解，并对 实际风电场进行日前风速预测实例分析，采用IE邸标准算例对MGEOPF模型进行分析。2. 根据权利要求1所述的一种基于日前风速预测的多时段电气互联系统最优潮流计算 方法，其特征在于:所述步骤1包括W下步骤： 步骤101:变分模态分解算法是变分问题的求解过程，该算法分为变分问题的构造和求 解;所述变分模态分解算法是将输入信号f分解成有限个子信号，该变分模态分解算法中本 征模态函数被重新定义为一个调幅-调频信号Uk(t)，表达为： Uk(t) =Ak(t)cos( Φ??(?)) (1) 式中：k为分解后的第k个子信号;Ak(t)为uk(t)的瞬时幅值；〇k(t)为uk(t)的瞬时频率， 〇>倘=脚怯)=d辦/舶ωk(t)表示为子信号角频率;辦货为相位。 设将原始信号分解为K个IMF分量，则对应的约束变分模型表达如下： 广Λ、. 式中：{11如）} = {山，...，朋}代表分解得到的1(个加。分量；{?如）} = {"1，...，《1(}表 示各分量的频率中屯、;Κ为分解后子信号总数； 步骤102:在利用高斯过程回归进行预测前，首先确定训练样本各分量的输入变量，采 用统计工具，即偏相关函数W及偏相关函数图来确定预测的输入变量； 步骤103:对于给定的数据集合i> = i(A·,，取籠 I，其中：输入数据矩阵XiGRd，输出数据矩 阵 yiER; 在给定数据D的有限集合中，f(χW)、f(χ(2))、…、f(χ(n))构成随机变量的一个集合，且 具有联合高斯分布，高斯过程的全部统计特征由均值函数m(x)和协方差函数k(x，x'）组成， 即： f(X)~GP(m(x)，k(x，x'）） (3) 将含噪声因素加入到观测目标值y中，建立高斯过程回归问题的一般模型，即 y = f(x)+e (4) 式中：ε为独立的高斯白噪声，符合高斯分布，均值为0,方差为护，即可记作?,'~护);由 于噪声ε为独立于f(x)的白噪声，当f(x)服从高斯分布，则y同样服从高斯分布，则y有限观 测值联合分布的集合形成一个高斯过程，即(5) 式中:m(x)为均值函数;Sij为Kronecker de；Ua函数，GP(X)为高斯分布；当i = j时，函数 Sij = l;协方差函数若W矩阵形式表达，即㈱ 式中：1表示NXN的单位矩阵，C(X，X)表示NXN的协方差矩阵，Κ(Χ，Χ)表示NXN的核矩 阵，称为Gram矩阵，其元素1(^ = ^町刮）； 根据贝叶斯原理，GP在给定数据D的集合内建立先验函数，在η*个给定测试数据集的f转变为后验分布，则训练数据向量X*的输出向量f和测试数据的输出向 量哀之间的联合高斯分布其中，预测均值向量京就是GP回归模型的输出，也即输出向量f*的预测值。3. 根据权利要求2所述的一种基于日前风速预测的多时段电气互联系统最优潮流计算 方法，其特征在于:所述步骤101中，对约束变分模型求解算法中引入二次惩罚项和拉格朗 日乘子λ W及惩罚参数α，得到增广Lagrange函数如下：利用交替方向乘子算法求取上述增广Lagrange函数的鞍点，即为约束变分模型的最优 解，从而将原始信号分解为K个窄带IMF分量;j为虚数符号。4. 根据权利要求2所述的一种基于日前风速预测的多时段电气互联系统最优潮流计算 方法，其特征在于：所述步骤102中，设XI为输出变量，如果在滞后k阶处偏相关系数超过 95%置信区间，即在区间[^-1..96/λ/^，1.96Λ游Γ]，则将xi-k作为其中一个输入变量；当所有 的偏相关系数都在95%区间外，此时取xi-i作为输入变量。5. 根据权利要求1所述的一种基于日前风速预测的多时段电气互联系统最优潮流计算 方法计算方法，其特征在于:所述步骤2中模型为： 步骤201:天然气从远距离的气源通过管道传输，再经过配气网络输送至终端用户；对 于绝热管道，设所述绝热管道的首末端节点分别为m和n，则稳态流量为：式中：Cmn为管道效率;JTmt，JTnt为t时段节点压力。 调节压力需要消耗额外的能量，所述额外的能量与流量和加压比成正比； 对于加压站k，设加压站k流入和流出节点分别为P和q，则加压站k能量消耗为：式中：化t为加压站消耗的功率;Bk为加压常数;Zk为压缩因子;Jlkt为等效消耗的流量;曰， 0，Y为效率转换常数;fkt为加压站k在t时段的流量； 天然气系统的任一节点应满足流量平衡方程，表示为：式中:Wgt为气源点注入量;Wdt为气负荷量;Pg,it为第i台发电机的有功出力；的&W)为燃 气轮机消耗的天然气量； 步骤202: W电气互联系统联合运行总成本最低为目标函数的MGEOPF模型如下： (1) 目标函数：式中：T，Nt，Nw，Ω S为时段、燃煤机组、风电场、气源点集合;Pw, it为风电场计划出力； Pav,it为风电场可用风功率;Ci(Pg,it)为第i台常规机组成本函数;Cg(Wg,mt)为风电出力成本 函数;C?(Pav,it-Pw,it)为风电低估惩罚成本;C?(Pw,it-Pav,it)为风电高估备用成本； 其中，式(17)第一项表示发电成本，W发电量的二次函数计算：式中：ai,bi,c功发电成本系数。 式（17)的第二项为风电成本，该项取决于风电场的所有者，第Ξ项表示天然气成本，可 表示为天然气流量的线性函数： Cg(Wg,mt) = gmWg,mt (19) (2) 电力系统约束： 电力系统多时段约束包括功率平衡约束、发电机出力约束、发电机爬坡约束、节点电压 约束、线路功率约束，表达如下：式中：Qg,it为第i台发电机无功出力;Pd,it，Qd,it为节点i电负荷;Vit，0it为节点i电压幅值 和相角，目ijt =目it-目节点巧日节点j之间的电导和电纳;RUi，RD功发电机爬坡率； Pi, t为线路1的功率;上标min和max表示最大和最小值。 (3)天然气系统约束： 天然气系统约束包括节点流量平衡方程、气源注入量约束、节点压力约束W及加压站 加压比约束，其表达如下：式中：Κ2ι，Κ?ι，Κ日功燃气轮机转换常数;为最小和最大气源供应量川mmin，jCax 为节点压力最小和最大值;Rc?in，REmax为最小和最大加压比。6. 根据权利要求1所述的一种基于日前风速预测的多时段电气互联系统最优潮流计算 方法，其特征在于:所述步骤3包括W下步骤： 引入风电惩罚成本用W表示低估风电带来的风能资源浪费惩罚，表示为：式中:kuw为惩罚成本系数(低估）;Pr为额定风功率。 引入风电备用成本用W表示高估风电带来的旋转备用配置成本，表示为：式中:k?为备用成本系数。7. 根据权利要求6所述的一种基于日前风速预测的多时段电气互联系统最优潮流计算 方法，其特征在于:所述步骤4中，建立的电气互联系统最优潮流模型为非线性规划模型，式 (22)和式(23)中积分项的引入，将式(22)和式(23)积分项转化为数值积分，进而采用非线 性规划求解器进行求解。
【文档编号】H02J3/06GK105978037SQ201610628139
【公开日】2016年9月28日
【申请日】2016年8月3日
【发明人】孙国强, 陈霜, 卫志农, 陈 胜, 张思德
【申请人】河海大学