CN107846024A - 孤岛微电网概率潮流解析计算方法 - Google Patents
孤岛微电网概率潮流解析计算方法 Download PDFInfo
- Publication number
- CN107846024A CN107846024A CN201711315865.7A CN201711315865A CN107846024A CN 107846024 A CN107846024 A CN 107846024A CN 201711315865 A CN201711315865 A CN 201711315865A CN 107846024 A CN107846024 A CN 107846024A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mtd
- mtr
- delta
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/04—Circuit arrangements for ac mains or ac distribution networks for connecting networks of the same frequency but supplied from different sources
- H02J3/06—Controlling transfer of power between connected networks; Controlling sharing of load between connected networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Power Engineering (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明提供一种孤岛微电网概率潮流的解析计算方法。为包含光伏、风电等间歇性DG以及负荷随机波动的孤岛微电网的运行状态分析提供计算方法。现有的孤岛微电网稳态潮流分析方法主要集中在常规潮流计算,仅能获取单一运行点的潮流结果,而关于概率潮流计算的研究也仅限于模拟法,即需要根据大量抽样样本计算统计获得概率潮流结果,计算时间长不利于在线分析,且多数未考虑微电网中DG出力的相关性。本发明方法针对这些局限,提供一种孤岛微电网概率潮流的解析计算方法,能够实现孤岛微电网概率潮流的快速计算,且精度能够符合工程精度要求,具备在线分析应用的可行性。
Description
技术领域
本发明属于孤岛微电网稳态分析领域,尤其涉及一种孤岛微电网概率潮流解析计算方法。
背景技术
微电网(microgrid,MG)为综合利用可再生能源提供了技术手段,根据与主网有无电气联系又可分为并网型与孤岛型,其中孤岛型微电网由于无平衡节点,缺乏主网支撑,通常采用对等控制方式,即由各采取下垂控制策略的分布式电源(distributed generator,DG)、储能装置(energystorage,ES)、静止无功发生器(STATCOM)等共同参与微电网的频率与电压的调节。
孤岛微电网潮流计算是对孤岛微电网规划、运行分析、优化调度的基础,目前已有专利或文献建立了孤岛微电网的潮流模型,并通过信赖域算法、类奔德斯分解方法等提高潮流计算收敛性与收敛速度。诸如光伏、风电这一类DG的出力具有较强的随机性,同时负荷也存在波动,相比于单一运行点的潮流结果,概率潮流(probabilistic load flow,PLF)能够获得各状态变量的概率分布信息,更为全面地反映微电网的运行状态。
求解PLF的方法主要分为模拟法及解析法:模拟法主要是指蒙特卡洛法,其缺陷是计算量过大,通常用于校验概率潮流准确性,不宜在线计算或用于优化调度。而解析法则是基于近似线性化的数学假设,根据不确定量数值特征求取近似解析表达式,仅需要在基准点进行一次潮流计算,计算效率高。目前有文献基于半不变量法求解主从控制模式下的孤岛微电网PLF,计算时将柴油发电机视为平衡节点,储能视为恒功率负荷,与传统电力***概率潮流计算类似,但未考虑对等控制方式下孤岛微电网的运行特点。
发明内容
针对现有技术存在的不足和空白,本发明采用以下技术方案:
一种孤岛微电网概率潮流解析计算方法,其特征在于,包括以下步骤:步骤一:输入原始数据,包括:微电网网架数据、电源和负荷的基本参数以及所对应的概率模型;
步骤二:由已知m个输入随机变量W=[w1,w2,……,wm]的边际概率分布F(W)及相关系数矩阵RW建立满足给定相关性水平的输入变量概率模型;
步骤三:根据各节点DG出力概率模型与负荷概率模型计算各节点功率注入期望值,基于孤岛微电网潮流模型进行基准运行点的稳态潮流计算,获得基准点潮流结果;
步骤四:计算各输入随机变量的各阶半不变量,并将具有相关性的随机变量的各阶半不变量转换为独立随机变量的各阶半不变量;
步骤五:计算节点状态变量、支路潮流变量以及具有下垂控制特性的电源输出功率的各阶半不变量;
步骤六:由Gram-Charlier级数展开获取输出变量的概率密度函数(probabilitydensity function,PDF),及累积分布函数(cumulative distribution function,CDF),输出PLF计算结果。
进一步地,步骤二的具体步骤为:
步骤二A:已知m个输入随机变量W=[w1,w2,……,wm]的边际概率分布F(W),并使用相关系数矩阵RW描述出力相关性,所述相关系数矩阵表达式为:
矩阵中各元素按下式计算:
(2)式中:ρwij为随机变量wi与wj之间的相关系数;Cov(wi,wj)为wi与wj的协方差;σwi与σwj为wi与wj的标准差;
建立满足给定相关性水平的输入变量概率模型的步骤为:
步骤二B:产生满足独立正态分布的样本矩阵S,其维数为m×N,N为样本规模;步骤二C:运用三阶多项式正态变换(third-order polynomial normal transformation,TPNT)以及求解代数方程的方法求取标准正态分布随机变量Y的等效相关系数矩阵RY;
步骤二D:对RY进行Cholesky分解获得其下三角矩阵LY,并由Y=LYS获得具有相关性的服从标准正态分布的样本矩阵Y;
步骤二E:通过等概率转换原则得到满足相关性水平的样本矩阵W,从而获得满足给定相关性水平的输入变量概率模型。
进一步地,步骤三中,所述基于孤岛微电网潮流模型进行基准运行点的稳态潮流计算包括以下步骤:
步骤三A:设定节点功率方程统一表达式:仅设定PQ节点与PV节点两种节点类型,其中PQ节点需列写有功与无功平衡方程,PV节点仅需列写有功平衡方程,方程形式为:
式中:PCi和QCi分别为节点i恒功率电源注入的有功和无功;PDi和QDi分别为节点i具有下垂控制特性的设备注入的有功和无功;PLi和QLi分别为节点i的有功与无功负荷;Pi和Qi分别为节点i注入的总有功功率和无功功率;
其中,所述具有下垂控制特性的设备注入的功率可统一表示为:
(4)、(5)、(6)、(7)式中:PDimax、PDimin、QDimax、QDimin分别为节点i具有下垂控制特性的设备注入的有功和无功的上下限,若无具有下垂控制特性的设备则均取为0;fmax、fmin、Umax、Umin分别为***频率、电压允许上下限;KDfi、KDUi分别为节点i对应的P-f、Q-U下垂系数;f0、f、U0、Ui分别为***频率、电压的空载值与节点i的实际值;
其中,对负荷进行建模可表示为:
(8)式中:PLNi和QLNi分别为节点i在额定工况下的有功与无功负荷;UNi和fN分别为节点i的额定电压与频率,取UNi=1p.u.,fN=1p.u.;Pi和Qi分别为节点i注入的总有功功率和无功功率;Api、Bpi、Cpi、Aqi、Bqi、Cqi分别为节点i的负荷有功与无功功率中恒阻抗型、恒电流型、恒功率型的百分比系数,各满足总和为1;kLpi、kLqi分别为节点i负荷的有功和无功功率的静态频率调节系数;
Ui为节点i***电压的实际值;f为节点i***频率的实际值;
其中,节点注入的有功与无功功率可表示为:
(9)式中:n为节点数目;Gij与Bij分别为节点导纳矩阵的实部与虚部;δij为节点i和节点j的相角差;Ui为节点i***电压的实际值;Uj为节点j***电压的实际值;
步骤三B:牛顿拉夫逊法求解:孤岛微电网中PQ节点与PV节点的数目分别为nPQ与nPV,则总计列写2nPQ+nPV个方程,采用牛顿拉夫逊法求解潮流方程组,其中,修正方程简写为:
(10)式中:ΔP、ΔQ为节点有功与无功不平衡量;Δf、Δδ、ΔU为分别为频率、相角、电压的修正量;J为雅克比矩阵,其分块矩阵分别为:
雅克比矩阵J中各分块矩阵的维数为:
各分块矩阵的表达式为:
(13)式中,各参数由(3)-(9)式中各参数代入计算。
进一步地,步骤四、步骤五的计算步骤为:
步骤四五A:在基准运行点处将潮流方程线性化并忽略2次及以上高次项,可得随机分量的关系为:
(14)式中:ΔW、ΔX、ΔZ分别为节点注入功率、节点状态变量以及支路潮流变量的随机分量;J0为基准点处的雅克比矩阵;X0为基准点处的节点状态量;
当已知W的各阶半不变量,则可根据半不变量的线性性质获得X、Z的各阶半不变量;
步骤四五B:通过对RW进行Cholesky分解获得其下三角矩阵LW,可将m个相关输入随机变量W转换为m个相互独立的变量Wind的组合,转换关系如下:
Wind=LW -1W (15)
相应的J0中与原始m个相关输入随机变量W对应列的元素需要通过LW修正;
步骤四五C:具有下垂控制特性的电源输出功率与***频率及节点电压有关,为获取其输出的概率分布,需要计算具有下垂控制特性的设备输出功率的的半不变量,假定在***运行范围内,对式(2)(4)在基准点线性化,可推得其关系为:
结合式(14)进一步转换为与节点注入功率ΔW的关系后同样可基于半不变量的线性性质由W的各阶半不变量求得具有下垂控制特性的电源输出功率的各阶半不变量。
进一步地,步骤六中获取输出变量的概率密度函数的步骤为:
步骤六A:借助标准正态分布概率密度函数φ以及它的各阶导数的线性组合来逼近最终待求的分布函数,根据已求得状态变量的各阶半不变量为γk,列出由各阶半不变量推导各阶中心矩的公式:
步骤六B:由半不变量γk求得中心矩后,进一步通过中心矩表示Gram-Charlier级数系数,公式为:
步骤六C:将随机变量标准化处理后借助Gram-Charlier级数展开可将随机变量的分布函数表达为由正态随机变量各阶导数组成的级数;
则由Gram-Charlier级数展开获取输出变量的概率密度函数及累积分布函数,输出PLF计算结果。
本发明的主要贡献是提出一种对等控制模式下孤岛微电网概率潮流的解析计算方法。首先综合DG、ES、STATCOM等设备的下垂控制特性建立节点功率方程统一表达式,简化运算与编程;随后在孤岛微电网潮流模型的基础上,扩展出概率潮流模型,运用Nataf变换对微电网中DG的出力相关性进行处理;最后基于半不变量与Gram-Charlier级数展开获得***频率及节点电压等状态变量的完整概率分布信息。
本发明的主要目的提供一种孤岛微电网概率潮流的解析计算方法。为包含光伏、风电等间歇性DG以及负荷随机波动的孤岛微电网的运行状态分析提供计算方法。现有的孤岛微电网稳态潮流分析方法主要集中在常规潮流计算,仅能获取单一运行点的潮流结果,而关于概率潮流计算的研究也仅限于模拟法,即需要根据大量抽样样本计算统计获得概率潮流结果,计算时间长不利于在线分析,且多数未考虑微电网中DG出力的相关性。本发明方法针对这些局限,提供一种孤岛微电网概率潮流的解析计算方法,能够实现孤岛微电网概率潮流的快速计算,且精度能够符合工程精度要求,具备在线分析应用的可行性。
本发明有益效果包括:
1)综合考虑DG、ES、STATCOM等设备的下垂控制特性建立节点功率方程统一表达式,简化了节点分类与编程,在***结构不发生变化情况下,即使发生设备投退的情况,也可不必修改潮流方程与雅克比矩阵结构,只需更新数值。
2)本发明方法的计算结果与蒙特卡洛法计算结果接近,同时在DG容量、出力相关性以及负荷波动的变化下,误差均在可接受范围内,能够满足工程精度要求。
3)本发明方法避免了抽样计算,计算时间短、效率高,与模拟法相比具有显著优势,具备应用于在线计算的可行性。因此在孤岛微电网运行分析、风险预警、优化调度等领域具有应用的可行性。
附图说明
下面结合附图和具体实施方式对本发明进一步详细的说明:
图1为本发明实施例流程示意图;
图2为本发明实施例中算例***网络节点示意图;
图3为本发明实施例中孤岛微电网频率PDF与CDF统计图;
图4为本发明实施例中节点1~18的电压期望值与置信区间统计图;
图5为本发明实施例中节点9STATCOM输出无功的PDF及CDF统计图;
图6为本发明实施例中节点18储能装置输出有功的PDF及CDF统计图;
图7为本发明实施例中调整前后的频率PDF与柴油发电机负载率PDF对比图。
具体实施方式
为让本专利的特征和优点能更明显易懂,下文特举实施例,作详细说明如下:
如图1所示,本发明实施例包括以下步骤:
步骤一:输入原始数据,包括:微电网网架数据、电源和负荷的基本参数以及所对应的概率模型;
步骤二:由已知m个输入随机变量W=[w1,w2,……,wm]的边际概率分布F(W)及相关系数矩阵RW建立满足给定相关性水平的输入变量概率模型;
步骤三:根据各节点DG出力概率模型与负荷概率模型计算各节点功率注入期望值,基于孤岛微电网潮流模型进行基准运行点的稳态潮流计算,获得基准点潮流结果;
步骤四:计算各输入随机变量的各阶半不变量,并将具有相关性的随机变量的各阶半不变量转换为独立随机变量的各阶半不变量;
步骤五:计算节点状态变量、支路潮流变量以及具有下垂控制特性的电源输出功率的各阶半不变量;
步骤六:由Gram-Charlier级数展开获取输出变量的概率密度函数(probabilitydensity function,PDF),及累积分布函数(cumulative distribution function,CDF),输出PLF计算结果。
其中,步骤二的具体步骤为:
步骤二A:已知m个输入随机变量W=[w1,w2,……,wm]的边际概率分布F(W),并使用相关系数矩阵RW描述出力相关性,相关系数矩阵表达式为:
矩阵中各元素按下式计算:
(2)式中:ρwij为随机变量wi与wj之间的相关系数;Cov(wi,wj)为wi与wj的协方差;σwi与σwj为wi与wj的标准差;
建立满足给定相关性水平的输入变量概率模型的步骤为:
步骤二B:产生满足独立正态分布的样本矩阵S,其维数为m×N,N为样本规模;步骤二C:运用三阶多项式正态变换(third-order polynomial normal transformation,TPNT)以及求解代数方程的方法求取标准正态分布随机变量Y的等效相关系数矩阵RY;
步骤二D:对RY进行Cholesky分解获得其下三角矩阵LY,并由Y=LYS获得具有相关性的服从标准正态分布的样本矩阵Y;
步骤二E:通过等概率转换原则得到满足相关性水平的样本矩阵W,从而获得满足给定相关性水平的输入变量概率模型。
其中,步骤三中,基于孤岛微电网潮流模型进行基准运行点的稳态潮流计算包括以下步骤:
步骤三A:设定节点功率方程统一表达式:仅设定PQ节点与PV节点两种节点类型,其中PQ节点需列写有功与无功平衡方程,PV节点仅需列写有功平衡方程,方程形式为:
式中:PCi和QCi分别为节点i恒功率电源注入的有功和无功,例如光伏、风电等采取MPPT控制方式的DG属于此类电源;PDi和QDi分别为节点i具有下垂控制特性的设备注入的有功和无功,例如柴油发电机、燃气轮机、储能、STATCOM等均属于此类电源;PLi和QLi分别为节点i的有功与无功负荷;Pi和Qi分别为节点i注入的总有功功率和无功功率;
其中,具有下垂控制特性的设备注入的功率可统一表示为:
(4)、(5)、(6)、(7)式中:PDimax、PDimin、QDimax、QDimin分别为节点i具有下垂控制特性的设备注入的有功和无功的上下限,不同设备所对应的参数取值及含义见表1,若无具有下垂控制特性的设备则均取为0;fmax、fmin、Umax、Umin分别为***频率、电压允许上下限;KDfi、KDUi分别为节点i对应的P-f、Q-U下垂系数;f0、f、U0、Ui分别为***频率、电压的空载值与节点i的实际值;
其中储能可以工作在充放电两种工况,在不考虑荷电状态的情况下,其PDimin、PDimax取决于储能的最大充放电功率,同样STATCOM可以吸收或发出无功,QDimin、QDimax取决于STATCOM配置容量,根据方程符号及物理含义,充电功率最大值以及STATCOM吸收无功最大值为一负值。
其中,对负荷进行建模时计及负荷的电压和频率静特性,负荷使用恒阻抗、恒电流和恒功率的组合模型来描述,可表示为:
(8)式中:PLNi和QLNi分别为节点i在额定工况下的有功与无功负荷;UNi和fN分别为节点i的额定电压与频率,取UNi=1p.u.,fN=1p.u.;Pi和Qi分别为节点i注入的总有功功率和无功功率;Api、Bpi、Cpi、Aqi、Bqi、Cqi分别为节点i的负荷有功与无功功率中恒阻抗型、恒电流型、恒功率型的百分比系数,各满足总和为1;kLpi、kLqi分别为节点i负荷的有功和无功功率的静态频率调节系数;
Ui为节点i***电压的实际值;f为节点i***频率的实际值;
其中,节点注入的有功与无功功率可表示为:
(9)式中:n为节点数目;Gij与Bij分别为节点导纳矩阵的实部与虚部;δij为节点i和节点j的相角差;Ui为节点i***电压的实际值;Uj为节点j***电压的实际值;
步骤三B:牛顿拉夫逊法求解:孤岛微电网中PQ节点与PV节点的数目分别为nPQ与nPV,则总计列写2nPQ+nPV个方程,采用牛顿拉夫逊法求解潮流方程组,其中,修正方程简写为:
(10)式中:ΔP、ΔQ为节点有功与无功不平衡量;Δf、Δδ、ΔU为分别为频率、相角、电压的修正量;J为雅克比矩阵,其分块矩阵分别为:
雅克比矩阵J中各分块矩阵的维数为:
各分块矩阵的表达式为:
(13)式中,各参数由(3)-(9)式中各参数代入计算。
其中,步骤四、步骤五的计算步骤为:
步骤四五A:在基准运行点处将潮流方程线性化并忽略2次及以上高次项,可得随机分量的关系为:
(14)式中:ΔW、ΔX、ΔZ分别为节点注入功率、节点状态变量以及支路潮流变量的随机分量;J0为基准点处的雅克比矩阵;X0为基准点处的节点状态量;
当已知W的各阶半不变量,则可根据半不变量的线性性质获得X、Z的各阶半不变量;
步骤四五B:通过对RW进行Cholesky分解获得其下三角矩阵LW,可将m个相关输入随机变量W转换为m个相互独立的变量Wind的组合,转换关系如下:
Wind=LW -1W (15)
相应的J0中与原始m个相关输入随机变量W对应列的元素需要通过LW修正;
步骤四五C:具有下垂控制特性的电源输出功率与***频率及节点电压有关,为获取其输出的概率分布,需要计算具有下垂控制特性的设备输出功率的的半不变量,假定在***运行范围内,对式(2)(4)在基准点线性化,可推得其关系为:
结合式(14)进一步转换为与节点注入功率ΔW的关系后同样可基于半不变量的线性性质由W的各阶半不变量求得具有下垂控制特性的电源输出功率的各阶半不变量。
其中,步骤六中获取输出变量的概率密度函数的步骤为:
步骤六A:借助标准正态分布概率密度函数φ以及它的各阶导数的线性组合来逼近最终待求的分布函数,根据已求得状态变量的各阶半不变量为γk,由于Gram-Charlier级数的各项系数与中心矩满足一定的关系,因此先列出由各阶半不变量推导各阶中心矩的公式:
步骤六B:由半不变量γk求得中心矩后,进一步通过中心矩表示Gram-Charlier级数系数,公式为:
步骤六C:将随机变量标准化处理后借助Gram-Charlier级数展开可将随机变量的分布函数表达为由正态随机变量各阶导数组成的级数。
则由Gram-Charlier级数展开获取输出变量的概率密度函数及累积分布函数,输出PLF计算结果。
具体地,本发明实施例运行在配置为Intel Pentium双核处理器,2.0G内存的计算机上,于MATLAB R2013a环境下编写程序实现本发明所提出的方法。以IEEE33节点***作为算例,将节点1与主网断开并配置DG,形成孤岛微电网。
算例***网络节点如图2所示。将***基准容量取1MVA,取孤岛微电网的安全运行范围为fmax=1.004p.u.,fmin=0.996p.u.,Umax=1.05p.u.,Umin=0.95p.u.。各节点恒阻抗、恒电流、恒功率负荷占比统一取为0.3、0.3、0.4,静态频率调节系数取kLpi=2、kLqi=-2。
负荷服从正态分布,负荷期望值见表1,负荷波动的标准差σload取为期望值的10%。
表1算例***负荷数据
在节点1、9、18接入同等容量的光伏电源,假定待分析时段内各光伏电源出力均满足在区间[400kW,600kW]范围的Beta分布,形状参数为α=1.693,β=5.162,各光伏电源之间的相关系数ρDG均取为为0.5。其他接入设备的容量参数见表2,各设备均采取下垂控制策略,下垂系数根据式(3)(5)进行整定。
表2设备参数
以蒙特卡洛法进行10000次抽样计算所得结果作为准确值,分析本发明实施例方法的误差情况,各节点电压与***频率的误差情况见表3。
表3各节点电压与***频率的误差情况
结果表明本发明实施例方法与蒙特卡洛法结果接近,能够满足工程精度的要求。利用本发明实施例方法以及蒙特卡洛法统计求出的的孤岛微电网频率PDF与CDF如图3所示。
本发明实施例方法计算得PDF与CDF和蒙特卡洛法统计得到的结果十分接近,***基本运行在额定频率附近,受光伏及负荷波动性影响可能存在偏离额定频率的情况,但偏离程度不大。
通过PLF计算可获得各节点电压的概率分布信息,节点1~18的电压期望值与置信区间统计图参见图3。
由图4可知,在光伏及负荷存在波动的环境下,PLF能够直观显示孤岛微电网可能得运行范围,并给出具有概率意义的信息,可为孤岛微电网的运行提供参考与决策依据。
此外,验证本发明实施例方法计算采取下垂控制策略的设备输出情况的准确性,其中节点9的STATCOM输出无功以及节点18储能装置输出有功的PDF与CDF参见图5与图6。
由图5和图6可知,本发明实施例方法所获得的设备输出功率的PDF与CDF与蒙特卡洛法所得结果基本接近。考虑节点9的STATCOM配置容量,由于并网点电压低于额定值,因此工作在发出无功的状态,且伴随着节点9电压的不确定性,同样输出无功也呈现不确定性。
而节点18的储能装置基本工作在以较小功率充电或放电的状态,主要由于此时***频率基本在额定值附近波动,且波动范围不大,相对而言储能也在充放电之间波动,且功率同样较小。由此可见本发明实施例方法能够较为准确地估计出具有下垂控制特性设备的可能输出范围,为评估设备运行状态提供参考。
在原算例条件的基础上,分别单独改变算例中光伏装机容量SDG、出力相关性ρDG以及负荷波动标准差σload,分析不同因素对算法误差的影响。其中表中SDG的数值指相对于原算例的光伏装机比例,表中εμmean、εμmax表示期望值的平均与最大偏差绝对值,εσmean、εσmax表示标准差的平均与最大偏差绝对值。不同场景下的误差变化情况见表4。
表4不同场景下的误差情况
可见随着并网光伏容量增加、出力相关性增强以及负荷波动增强,相应的算法误差都会有所增大,主要是由于半不变量法在基准运行点线性化,当输入变量的不确定性增强会引起误差增大。但数据表明误差均处于可接受范围内,都能够满足工程精度要求。
而在计算时间上,本发明实施例方法仅需0.06s即可完成PLF计算,而抽样规模为10000次的蒙特卡洛模拟法需要163.27s,因此本发明实施例方法在允许一定误差的前提下能够在计算时间上具有显著优势,具备应用于在线计算的可行性。
进一步针对上述算例***就本发明PLF方法在孤岛微电网的运行分析中的应用进行举例说明。以柴油发电机组的启停为例,原算例中节点33有4台额定功率500kW的柴油发电机并网运行,***频率与柴油发电机负载率的PDF如图中浅色曲线所示。柴油发电机最低负载要求通常在30%左右,当负载率较低时发电成本较高,经济效益差。由图7可知,柴油发电机此时负载率大致处于30%~45%的范围,负载率较低,考虑进行调整,运用本发明实施例方法计算退出1台柴油发电机后的PLF,调整后***频率与柴油发电机负载率PDF见深色曲线。
通过对比可见,尽管***频率略低于调整前,但仍大致处于0.999p.u.~1.000p.u.之间,频率完全合格,而柴油发电机的负载率均提升到55%左右,较为接近最佳负载状态,经济性较好。通过PLF结果可以判定调整后在***安全性不受到影响的情况下能够提高经济性,因此可以令1台柴油发电机退出运行。
本专利不局限于上述最佳实施方式,任何人在本专利的启示下都可以得出其它各种形式的孤岛微电网概率潮流解析计算方法,凡依本发明申请专利范围所做的均等变化与修饰,皆应属本专利的涵盖范围。
Claims (5)
1.一种孤岛微电网概率潮流解析计算方法,其特征在于,包括以下步骤:
步骤一:输入原始数据,包括:微电网网架数据、电源和负荷的基本参数以及所对应的概率模型;
步骤二:由已知m个输入随机变量W=[w1,w2,……,wm]的边际概率分布F(W)及相关系数矩阵RW建立满足给定相关性水平的输入变量概率模型;
步骤三:根据各节点DG出力概率模型与负荷概率模型计算各节点功率注入期望值,基于孤岛微电网潮流模型进行基准运行点的稳态潮流计算,获得基准点潮流结果;
步骤四:计算各输入随机变量的各阶半不变量,并将具有相关性的随机变量的各阶半不变量转换为独立随机变量的各阶半不变量;
步骤五:计算节点状态变量、支路潮流变量以及具有下垂控制特性的电源输出功率的各阶半不变量;
步骤六:由Gram-Charlier级数展开获取输出变量的概率密度函数(probabilitydensity function,PDF),及累积分布函数(cumulative distribution function,CDF),输出PLF计算结果。
2.根据权利要求1所述的孤岛微电网概率潮流解析计算方法,其特征在于,步骤二的具体步骤为:
步骤二A:已知m个输入随机变量W=[w1,w2,……,wm]的边际概率分布F(W),并使用相关系数矩阵RW描述出力相关性,所述相关系数矩阵表达式为:
矩阵中各元素按下式计算:
<mrow>
<msub>
<mi>&rho;</mi>
<msub>
<mi>w</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
</msub>
<mo>=</mo>
<msub>
<mi>&rho;</mi>
<msub>
<mi>w</mi>
<mrow>
<mi>j</mi>
<mi>i</mi>
</mrow>
</msub>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>C</mi>
<mrow>
<mi>o</mi>
<mi>v</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>w</mi>
<mi>i</mi>
</msub>
<mo>,</mo>
<msub>
<mi>w</mi>
<mi>j</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>&sigma;</mi>
<msub>
<mi>w</mi>
<mi>i</mi>
</msub>
</msub>
<msub>
<mi>&sigma;</mi>
<msub>
<mi>w</mi>
<mi>j</mi>
</msub>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
(2)式中:ρwij为随机变量wi与wj之间的相关系数;Cov(wi,wj)为wi与wj的协方差;σwi与σwj为wi与wj的标准差;
建立满足给定相关性水平的输入变量概率模型的步骤为:
步骤二B:产生满足独立正态分布的样本矩阵S,其维数为m×N,N为样本规模;步骤二C:运用三阶多项式正态变换(third-order polynomial normal transformation,TPNT)以及求解代数方程的方法求取标准正态分布随机变量Y的等效相关系数矩阵RY;
步骤二D:对RY进行Cholesky分解获得其下三角矩阵LY,并由Y=LYS获得具有相关性的服从标准正态分布的样本矩阵Y;
步骤二E:通过等概率转换原则得到满足相关性水平的样本矩阵W,从而获得满足给定相关性水平的输入变量概率模型。
3.根据权利要求2所述的孤岛微电网概率潮流解析计算方法,其特征在于:步骤三中,所述基于孤岛微电网潮流模型进行基准运行点的稳态潮流计算包括以下步骤:
步骤三A:设定节点功率方程统一表达式:仅设定PQ节点与PV节点两种节点类型,其中PQ节点需列写有功与无功平衡方程,PV节点仅需列写有功平衡方程,方程形式为:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&Delta;P</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>C</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>D</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>L</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>P</mi>
<mi>i</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&Delta;Q</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<msub>
<mi>Q</mi>
<mrow>
<mi>C</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>Q</mi>
<mrow>
<mi>D</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>Q</mi>
<mrow>
<mi>L</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>Q</mi>
<mi>i</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
式中:PCi和QCi分别为节点i恒功率电源注入的有功和无功;PDi和QDi分别为节点i具有下垂控制特性的设备注入的有功和无功;PLi和QLi分别为节点i的有功与无功负荷;Pi和Qi分别为节点i注入的总有功功率和无功功率;
其中,所述具有下垂控制特性的设备注入的功率可统一表示为:
<mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>D</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<msub>
<mi>P</mi>
<mrow>
<mi>D</mi>
<mi>i</mi>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
</mtd>
<mtd>
<mrow>
<mi>f</mi>
<mo>&le;</mo>
<msub>
<mi>f</mi>
<mi>min</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>D</mi>
<mi>i</mi>
<mi>min</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>K</mi>
<mrow>
<mi>D</mi>
<mi>f</mi>
<mi>i</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<mi>f</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>f</mi>
<mi>min</mi>
</msub>
<mo><</mo>
<mi>f</mi>
<mo><</mo>
<msub>
<mi>f</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>P</mi>
<mrow>
<mi>D</mi>
<mi>i</mi>
<mi>min</mi>
</mrow>
</msub>
</mtd>
<mtd>
<mrow>
<mi>f</mi>
<mo>&GreaterEqual;</mo>
<msub>
<mi>f</mi>
<mi>max</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>K</mi>
<mrow>
<mi>D</mi>
<mi>f</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>D</mi>
<mi>i</mi>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>D</mi>
<mi>i</mi>
<mi>min</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>f</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>f</mi>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>Q</mi>
<mrow>
<mi>D</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<msub>
<mi>Q</mi>
<mrow>
<mi>D</mi>
<mi>i</mi>
<mi>max</mi>
</mrow>
</msub>
</mtd>
<mtd>
<mrow>
<msub>
<mi>U</mi>
<mi>i</mi>
</msub>
<mo>&le;</mo>
<msub>
<mi>U</mi>
<mi>min</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>Q</mi>
<mrow>
<mi>D</mi>
<mi>i</mi>
<mi>min</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>K</mi>
<mrow>
<mi>D</mi>
<mi>U</mi>
<mi>i</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>U</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<msub>
<mi>U</mi>
<mi>i</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>U</mi>
<mi>min</mi>
</msub>
<mo><</mo>
<msub>
<mi>U</mi>
<mi>i</mi>
</msub>
<mo><</mo>
<msub>
<mi>U</mi>
<mi>max</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>Q</mi>
<mrow>
<mi>D</mi>
<mi>i</mi>
<mi>min</mi>
</mrow>
</msub>
</mtd>
<mtd>
<mrow>
<msub>
<mi>U</mi>
<mi>i</mi>
</msub>
<mo>&GreaterEqual;</mo>
<msub>
<mi>U</mi>
<mi>max</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>K</mi>
<mrow>
<mi>D</mi>
<mi>U</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>Q</mi>
<mrow>
<mi>D</mi>
<mi>i</mi>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>Q</mi>
<mrow>
<mi>D</mi>
<mi>i</mi>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>U</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>U</mi>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
(4)、(5)、(6)、(7)式中:PDimax、PDimin、QDimax、QDimin分别为节点i具有下垂控制特性的设备注入的有功和无功的上下限,若无具有下垂控制特性的设备则均取为0;fmax、fmin、Umax、Umin分别为***频率、电压允许上下限;KDfi、KDUi分别为节点i对应的P-f、Q-U下垂系数;f0、f、U0、Ui分别为***频率、电压的空载值与节点i的实际值;
其中,对负荷进行建模可表示为:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>L</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>L</mi>
<mi>N</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>&lsqb;</mo>
<msub>
<mi>A</mi>
<mrow>
<mi>p</mi>
<mi>i</mi>
</mrow>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>U</mi>
<mi>i</mi>
</msub>
<msub>
<mi>U</mi>
<mrow>
<mi>N</mi>
<mi>i</mi>
</mrow>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>B</mi>
<mrow>
<mi>p</mi>
<mi>i</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>U</mi>
<mi>i</mi>
</msub>
<msub>
<mi>U</mi>
<mrow>
<mi>N</mi>
<mi>i</mi>
</mrow>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>C</mi>
<mrow>
<mi>p</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>&rsqb;</mo>
<mo>&lsqb;</mo>
<mn>1</mn>
<mo>+</mo>
<msub>
<mi>k</mi>
<mrow>
<mi>L</mi>
<mi>p</mi>
<mi>i</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mo>-</mo>
<msub>
<mi>f</mi>
<mi>N</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>Q</mi>
<mrow>
<mi>L</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>Q</mi>
<mrow>
<mi>L</mi>
<mi>N</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>&lsqb;</mo>
<msub>
<mi>A</mi>
<mrow>
<mi>q</mi>
<mi>i</mi>
</mrow>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>U</mi>
<mi>i</mi>
</msub>
<msub>
<mi>U</mi>
<mrow>
<mi>N</mi>
<mi>i</mi>
</mrow>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>B</mi>
<mrow>
<mi>q</mi>
<mi>i</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>U</mi>
<mi>i</mi>
</msub>
<msub>
<mi>U</mi>
<mrow>
<mi>N</mi>
<mi>i</mi>
</mrow>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>C</mi>
<mrow>
<mi>q</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>&rsqb;</mo>
<mo>&lsqb;</mo>
<mn>1</mn>
<mo>+</mo>
<msub>
<mi>k</mi>
<mrow>
<mi>L</mi>
<mi>q</mi>
<mi>i</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mo>-</mo>
<msub>
<mi>f</mi>
<mi>N</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
(8)式中:PLNi和QLNi分别为节点i在额定工况下的有功与无功负荷;UNi和fN分别为节点i的额定电压与频率,取UNi=1p.u.,fN=1p.u.;Pi和Qi分别为节点i注入的总有功功率和无功功率;Api、Bpi、Cpi、Aqi、Bqi、Cqi分别为节点i的负荷有功与无功功率中恒阻抗型、恒电流型、恒功率型的百分比系数,各满足总和为1;kLpi、kLqi分别为节点i负荷的有功和无功功率的静态频率调节系数;Ui为节点i***电压的实际值;f为节点i***频率的实际值;
其中,节点注入的有功与无功功率可表示为:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>P</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<msub>
<mi>U</mi>
<mi>i</mi>
</msub>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>U</mi>
<mi>j</mi>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>G</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msub>
<mi>cos&delta;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>B</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msub>
<mi>sin&delta;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>Q</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<msub>
<mi>U</mi>
<mi>i</mi>
</msub>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>U</mi>
<mi>j</mi>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>G</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msub>
<mi>sin&delta;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>B</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msub>
<mi>cos&delta;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
(9)式中:n为节点数目;Gij与Bij分别为节点导纳矩阵的实部与虚部;δij为节点i和节点j的相角差;Ui为节点i***电压的实际值;Uj为节点j***电压的实际值;
步骤三B:牛顿拉夫逊法求解:孤岛微电网中PQ节点与PV节点的数目分别为nPQ与nPV,则总计列写2nPQ+nPV个方程,采用牛顿拉夫逊法求解潮流方程组,其中,修正方程简写为:
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mi>&Delta;</mi>
<mi>P</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>&Delta;</mi>
<mi>Q</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mi>J</mi>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mi>&Delta;</mi>
<mi>f</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>&Delta;</mi>
<mi>&delta;</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>&Delta;</mi>
<mi>U</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mi>E</mi>
</mtd>
<mtd>
<mi>H</mi>
</mtd>
<mtd>
<mi>N</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>F</mi>
</mtd>
<mtd>
<mi>M</mi>
</mtd>
<mtd>
<mi>L</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mi>&Delta;</mi>
<mi>f</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>&Delta;</mi>
<mi>&delta;</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>&Delta;</mi>
<mi>U</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>10</mn>
<mo>)</mo>
</mrow>
</mrow>
(10)式中:ΔP、ΔQ为节点有功与无功不平衡量;Δf、Δδ、ΔU为分别为频率、相角、电压的修正量;J为雅克比矩阵,其分块矩阵分别为:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mi>E</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&Delta;</mi>
<mi>P</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>f</mi>
</mrow>
</mfrac>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>H</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&Delta;</mi>
<mi>P</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>&delta;</mi>
</mrow>
</mfrac>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>N</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&Delta;</mi>
<mi>P</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>U</mi>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>F</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&Delta;</mi>
<mi>Q</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>f</mi>
</mrow>
</mfrac>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>M</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&Delta;</mi>
<mi>Q</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>&delta;</mi>
</mrow>
</mfrac>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>L</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&Delta;</mi>
<mi>Q</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>U</mi>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>11</mn>
<mo>)</mo>
</mrow>
</mrow>
雅克比矩阵J中各分块矩阵的维数为:
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&Delta;P</mi>
<mrow>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>n</mi>
<mrow>
<mi>P</mi>
<mi>Q</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>n</mi>
<mrow>
<mi>P</mi>
<mi>V</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mo>&times;</mo>
<mn>1</mn>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&Delta;Q</mi>
<mrow>
<msub>
<mi>n</mi>
<mrow>
<mi>P</mi>
<mi>Q</mi>
</mrow>
</msub>
<mo>&times;</mo>
<mn>1</mn>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>E</mi>
<mrow>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>n</mi>
<mrow>
<mi>P</mi>
<mi>Q</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>n</mi>
<mrow>
<mi>P</mi>
<mi>V</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mo>&times;</mo>
<mn>1</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>H</mi>
<mrow>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>n</mi>
<mrow>
<mi>P</mi>
<mi>Q</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>n</mi>
<mrow>
<mi>P</mi>
<mi>V</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mo>&times;</mo>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>n</mi>
<mrow>
<mi>P</mi>
<mi>Q</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>n</mi>
<mrow>
<mi>P</mi>
<mi>V</mi>
</mrow>
</msub>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>N</mi>
<mrow>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>n</mi>
<mrow>
<mi>P</mi>
<mi>Q</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>n</mi>
<mrow>
<mi>P</mi>
<mi>V</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mo>&times;</mo>
<msub>
<mi>n</mi>
<mrow>
<mi>P</mi>
<mi>Q</mi>
</mrow>
</msub>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>F</mi>
<mrow>
<msub>
<mi>n</mi>
<mrow>
<mi>P</mi>
<mi>Q</mi>
</mrow>
</msub>
<mo>&times;</mo>
<mn>1</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>M</mi>
<mrow>
<msub>
<mi>n</mi>
<mrow>
<mi>P</mi>
<mi>Q</mi>
</mrow>
</msub>
<mo>&times;</mo>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>n</mi>
<mrow>
<mi>P</mi>
<mi>Q</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>n</mi>
<mrow>
<mi>P</mi>
<mi>V</mi>
</mrow>
</msub>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>L</mi>
<mrow>
<msub>
<mi>n</mi>
<mrow>
<mi>P</mi>
<mi>Q</mi>
</mrow>
</msub>
<mo>&times;</mo>
<msub>
<mi>n</mi>
<mrow>
<mi>P</mi>
<mi>Q</mi>
</mrow>
</msub>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&Delta;f</mi>
<mrow>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>&times;</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&Delta;&delta;</mi>
<mrow>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>n</mi>
<mrow>
<mi>P</mi>
<mi>Q</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>n</mi>
<mrow>
<mi>P</mi>
<mi>V</mi>
</mrow>
</msub>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
</mrow>
<mo>&times;</mo>
<mn>1</mn>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&Delta;U</mi>
<mrow>
<msub>
<mi>n</mi>
<mrow>
<mi>P</mi>
<mi>Q</mi>
</mrow>
</msub>
<mo>&times;</mo>
<mn>1</mn>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>12</mn>
<mo>)</mo>
</mrow>
</mrow>
各分块矩阵的表达式为:
<mrow>
<mi>E</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&Delta;P</mi>
<mi>i</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>f</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mo>-</mo>
<msub>
<mi>K</mi>
<mrow>
<mi>D</mi>
<mi>f</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>L</mi>
<mi>N</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>&lsqb;</mo>
<msub>
<mi>A</mi>
<mrow>
<mi>p</mi>
<mi>i</mi>
</mrow>
</msub>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>U</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>B</mi>
<mrow>
<mi>p</mi>
<mi>i</mi>
</mrow>
</msub>
<msub>
<mi>U</mi>
<mi>i</mi>
</msub>
<mo>+</mo>
<msub>
<mi>C</mi>
<mrow>
<mi>p</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>&rsqb;</mo>
<msub>
<mi>k</mi>
<mrow>
<mi>L</mi>
<mi>p</mi>
<mi>i</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mi>H</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&Delta;P</mi>
<mi>i</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&delta;</mi>
<mi>j</mi>
</msub>
</mrow>
</mfrac>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>U</mi>
<mi>i</mi>
</msub>
<msub>
<mi>U</mi>
<mi>j</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>G</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msub>
<mi>sin&delta;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>B</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msub>
<mi>cos&delta;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>i</mi>
<mo>&NotEqual;</mo>
<mi>j</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>U</mi>
<mi>i</mi>
<mn>2</mn>
</msubsup>
<msub>
<mi>B</mi>
<mrow>
<mi>i</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>Q</mi>
<mi>i</mi>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mi>j</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
<mrow>
<mi>M</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&Delta;Q</mi>
<mi>i</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&delta;</mi>
<mi>j</mi>
</msub>
</mrow>
</mfrac>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>U</mi>
<mi>i</mi>
</msub>
<msub>
<mi>U</mi>
<mi>j</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>G</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msub>
<mi>cos&delta;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>B</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msub>
<mi>sin&delta;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>i</mi>
<mo>&NotEqual;</mo>
<mi>j</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>U</mi>
<mi>i</mi>
<mn>2</mn>
</msubsup>
<msub>
<mi>G</mi>
<mrow>
<mi>i</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>P</mi>
<mi>i</mi>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mi>j</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
<mrow>
<mi>N</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&Delta;P</mi>
<mi>i</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>U</mi>
<mi>j</mi>
</msub>
</mrow>
</mfrac>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>U</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>G</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msub>
<mi>cos&delta;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>B</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msub>
<mi>sin&delta;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>i</mi>
<mo>&NotEqual;</mo>
<mi>j</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>L</mi>
<mi>N</mi>
<mi>i</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>2</mn>
<msub>
<mi>A</mi>
<mrow>
<mi>p</mi>
<mi>i</mi>
</mrow>
</msub>
<msub>
<mi>U</mi>
<mi>i</mi>
</msub>
<mo>+</mo>
<msub>
<mi>B</mi>
<mrow>
<mi>p</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>&lsqb;</mo>
<mn>1</mn>
<mo>+</mo>
<msub>
<mi>k</mi>
<mrow>
<mi>L</mi>
<mi>p</mi>
<mi>i</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mo>-</mo>
<msub>
<mi>f</mi>
<mi>N</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>-</mo>
<msub>
<mi>U</mi>
<mi>i</mi>
</msub>
<msub>
<mi>G</mi>
<mrow>
<mi>i</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>-</mo>
<mfrac>
<msub>
<mi>P</mi>
<mi>i</mi>
</msub>
<msub>
<mi>U</mi>
<mi>i</mi>
</msub>
</mfrac>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mi>j</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
<mrow>
<mi>L</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&Delta;Q</mi>
<mi>i</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>U</mi>
<mi>j</mi>
</msub>
</mrow>
</mfrac>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>U</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>G</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msub>
<mi>sin&delta;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>B</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msub>
<mi>cos&delta;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>i</mi>
<mo>&NotEqual;</mo>
<mi>j</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>K</mi>
<mrow>
<mi>D</mi>
<mi>U</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>Q</mi>
<mrow>
<mi>L</mi>
<mi>N</mi>
<mi>i</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>2</mn>
<msub>
<mi>A</mi>
<mrow>
<mi>q</mi>
<mi>i</mi>
</mrow>
</msub>
<msub>
<mi>U</mi>
<mi>i</mi>
</msub>
<mo>+</mo>
<msub>
<mi>B</mi>
<mrow>
<mi>q</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>&lsqb;</mo>
<mn>1</mn>
<mo>+</mo>
<msub>
<mi>k</mi>
<mrow>
<mi>L</mi>
<mi>q</mi>
<mi>i</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mo>-</mo>
<msub>
<mi>f</mi>
<mi>N</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>+</mo>
<msub>
<mi>U</mi>
<mi>i</mi>
</msub>
<msub>
<mi>B</mi>
<mrow>
<mi>i</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>-</mo>
<mfrac>
<msub>
<mi>Q</mi>
<mi>i</mi>
</msub>
<msub>
<mi>U</mi>
<mi>i</mi>
</msub>
</mfrac>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mi>j</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>13</mn>
<mo>)</mo>
</mrow>
</mrow>
(13)式中,各参数由(3)-(9)式中各参数代入计算。
4.根据权利要求3所述的孤岛微电网概率潮流解析计算方法,其特征在于:步骤四、步骤五的计算步骤为:
步骤四五A:在基准运行点处将潮流方程线性化并忽略2次及以上高次项,可得随机分量的关系为:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>&Delta;</mi>
<mi>X</mi>
<mo>=</mo>
<msup>
<msub>
<mi>J</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mi>&Delta;</mi>
<mi>W</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>&Delta;</mi>
<mi>Z</mi>
<mo>=</mo>
<msub>
<mi>G</mi>
<mn>0</mn>
</msub>
<msup>
<msub>
<mi>J</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mi>&Delta;</mi>
<mi>W</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>G</mi>
<mn>0</mn>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mo>&part;</mo>
<mi>Z</mi>
<mo>/</mo>
<mo>&part;</mo>
<mi>X</mi>
<mo>)</mo>
</mrow>
<msub>
<mo>|</mo>
<mrow>
<mi>X</mi>
<mo>=</mo>
<msub>
<mi>X</mi>
<mn>0</mn>
</msub>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>14</mn>
<mo>)</mo>
</mrow>
</mrow>
(14)式中:ΔW、ΔX、ΔZ分别为节点注入功率、节点状态变量以及支路潮流变量的随机分量;J0为基准点处的雅克比矩阵;X0为基准点处的节点状态量;
当已知W的各阶半不变量,则可根据半不变量的线性性质获得X、Z的各阶半不变量;
步骤四五B:通过对RW进行Cholesky分解获得其下三角矩阵LW,可将m个相关输入随机变量W转换为m个相互独立的变量Wind的组合,转换关系如下:
Wind=LW -1W (15)
相应的J0中与原始m个相关输入随机变量W对应列的元素需要通过LW修正;
步骤四五C:具有下垂控制特性的电源输出功率与***频率及节点电压有关,为获取其输出的概率分布,需要计算具有下垂控制特性的设备输出功率的的半不变量,假定在***运行范围内,对式(2)(4)在基准点线性化,可推得其关系为:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&Delta;P</mi>
<mrow>
<mi>D</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>K</mi>
<mrow>
<mi>D</mi>
<mi>f</mi>
<mi>i</mi>
</mrow>
</msub>
<mi>&Delta;</mi>
<mi>f</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&Delta;Q</mi>
<mrow>
<mi>D</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>K</mi>
<mrow>
<mi>D</mi>
<mi>U</mi>
<mi>i</mi>
</mrow>
</msub>
<msub>
<mi>&Delta;U</mi>
<mi>i</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>16</mn>
<mo>)</mo>
</mrow>
</mrow>
结合式(14)进一步转换为与节点注入功率ΔW的关系后同样可基于半不变量的线性性质由W的各阶半不变量求得具有下垂控制特性的电源输出功率的各阶半不变量。
5.根据权利要求4所述的孤岛微电网概率潮流解析计算方法,其特征在于:步骤六中获取输出变量的概率密度函数的步骤为:
步骤六A:借助标准正态分布概率密度函数φ以及它的各阶导数的线性组合来逼近最终待求的分布函数,根据已求得状态变量的各阶半不变量为γk,列出由各阶半不变量推导各阶中心矩的公式:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&beta;</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&beta;</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<msub>
<mi>&gamma;</mi>
<mn>2</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&beta;</mi>
<mn>3</mn>
</msub>
<mo>=</mo>
<msub>
<mi>&gamma;</mi>
<mn>3</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&beta;</mi>
<mn>4</mn>
</msub>
<mo>=</mo>
<msub>
<mi>&gamma;</mi>
<mn>4</mn>
</msub>
<mo>+</mo>
<mn>3</mn>
<msup>
<msub>
<mi>&gamma;</mi>
<mn>2</mn>
</msub>
<mn>2</mn>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&beta;</mi>
<mn>5</mn>
</msub>
<mo>=</mo>
<msub>
<mi>&gamma;</mi>
<mn>5</mn>
</msub>
<mo>+</mo>
<mn>10</mn>
<msub>
<mi>&gamma;</mi>
<mn>3</mn>
</msub>
<msub>
<mi>&gamma;</mi>
<mn>2</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&beta;</mi>
<mn>6</mn>
</msub>
<mo>=</mo>
<msub>
<mi>&gamma;</mi>
<mn>6</mn>
</msub>
<mo>+</mo>
<mn>15</mn>
<msub>
<mi>&gamma;</mi>
<mn>4</mn>
</msub>
<msub>
<mi>&gamma;</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<mn>10</mn>
<msup>
<msub>
<mi>&gamma;</mi>
<mn>3</mn>
</msub>
<mn>2</mn>
</msup>
<mo>+</mo>
<mn>15</mn>
<msup>
<msub>
<mi>&gamma;</mi>
<mn>2</mn>
</msub>
<mn>3</mn>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&beta;</mi>
<mn>7</mn>
</msub>
<mo>=</mo>
<msub>
<mi>&gamma;</mi>
<mn>7</mn>
</msub>
<mo>+</mo>
<mn>21</mn>
<msub>
<mi>&gamma;</mi>
<mn>5</mn>
</msub>
<msub>
<mi>&gamma;</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<mn>35</mn>
<msub>
<mi>&gamma;</mi>
<mn>4</mn>
</msub>
<msub>
<mi>&gamma;</mi>
<mn>3</mn>
</msub>
<mo>+</mo>
<mn>105</mn>
<msub>
<mi>&gamma;</mi>
<mn>3</mn>
</msub>
<msup>
<msub>
<mi>&gamma;</mi>
<mn>2</mn>
</msub>
<mn>2</mn>
</msup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>17</mn>
<mo>)</mo>
</mrow>
</mrow>
步骤六B:由半不变量γk求得中心矩后,进一步通过中心矩表示Gram-Charlier级数系数,公式为:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>c</mi>
<mn>0</mn>
</msub>
<mo>=</mo>
<mn>1</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>c</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>c</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>&beta;</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mn>2</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>c</mi>
<mn>3</mn>
</msub>
<mo>=</mo>
<mo>-</mo>
<msub>
<mi>&beta;</mi>
<mn>3</mn>
</msub>
<mo>/</mo>
<mn>6</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>c</mi>
<mn>4</mn>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>&beta;</mi>
<mn>4</mn>
</msub>
<mo>-</mo>
<mn>6</mn>
<msub>
<mi>&beta;</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mn>24</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>c</mi>
<mn>5</mn>
</msub>
<mo>=</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>&beta;</mi>
<mn>5</mn>
</msub>
<mo>-</mo>
<mn>10</mn>
<msub>
<mi>&beta;</mi>
<mn>3</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mn>120</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>c</mi>
<mn>6</mn>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>&beta;</mi>
<mn>6</mn>
</msub>
<mo>-</mo>
<mn>15</mn>
<msub>
<mi>&beta;</mi>
<mn>4</mn>
</msub>
<mo>+</mo>
<mn>45</mn>
<msub>
<mi>&beta;</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<mn>15</mn>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mn>720</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>c</mi>
<mn>7</mn>
</msub>
<mo>=</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>&beta;</mi>
<mn>7</mn>
</msub>
<mo>-</mo>
<mn>21</mn>
<msub>
<mi>&beta;</mi>
<mn>5</mn>
</msub>
<mo>+</mo>
<mn>105</mn>
<msub>
<mi>&beta;</mi>
<mn>3</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mn>5040</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>18</mn>
<mo>)</mo>
</mrow>
</mrow>
步骤六C:将随机变量标准化处理后借助Gram-Charlier级数展开可将随机变量的分布函数表达为由正态随机变量各阶导数组成的级数;
则由Gram-Charlier级数展开获取输出变量的概率密度函数及累积分布函数,输出PLF计算结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711315865.7A CN107846024B (zh) | 2017-12-11 | 2017-12-11 | 孤岛微电网概率潮流解析计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711315865.7A CN107846024B (zh) | 2017-12-11 | 2017-12-11 | 孤岛微电网概率潮流解析计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107846024A true CN107846024A (zh) | 2018-03-27 |
CN107846024B CN107846024B (zh) | 2020-12-25 |
Family
ID=61663761
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711315865.7A Active CN107846024B (zh) | 2017-12-11 | 2017-12-11 | 孤岛微电网概率潮流解析计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107846024B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108964061A (zh) * | 2018-07-23 | 2018-12-07 | 长沙理工大学 | 一种考虑负荷频率电压静特性的含风电交直流电力***概率动态连续潮流新方法 |
CN109347110A (zh) * | 2018-10-29 | 2019-02-15 | 河海大学 | 一种含高比例风电并网的自适应线性化概率潮流计算方法 |
CN109672167A (zh) * | 2019-01-29 | 2019-04-23 | 福州大学 | 基于改进下垂算法的用于光伏低压微电网的功率控制方法 |
CN109950935A (zh) * | 2018-11-08 | 2019-06-28 | 国网甘肃省电力公司电力科学研究院 | 孤岛运行的交直流混联微电网概率潮流方法 |
CN110112790A (zh) * | 2019-06-06 | 2019-08-09 | 南方电网科学研究院有限责任公司 | 有源配电网运行效率的评价方法、装置、设备及存储介质 |
CN110276135A (zh) * | 2019-06-25 | 2019-09-24 | 华北电力大学 | 一种并网停车场的可用容量确定方法、装置及计算设备 |
CN112491056A (zh) * | 2020-11-23 | 2021-03-12 | 广东电网有限责任公司电网规划研究中心 | 一种交直流混联电网的概率潮流计算方法、装置和介质 |
CN113315132A (zh) * | 2021-06-02 | 2021-08-27 | 贵州电网有限责任公司 | 一种含下垂节点的孤岛微电网三相潮流计算方法 |
CN114936454A (zh) * | 2022-05-18 | 2022-08-23 | 东北电力大学 | 一种基于半不变量法的交直流混联***概率潮流计算方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105634008A (zh) * | 2016-01-20 | 2016-06-01 | 华南理工大学 | 考虑风速强度的分布式风力发电并网后概率线损分摊方法 |
US20170317498A1 (en) * | 2016-05-02 | 2017-11-02 | Nec Laboratories America, Inc. | Resiliency Controller for Voltage Regulation in Microgrids |
-
2017
- 2017-12-11 CN CN201711315865.7A patent/CN107846024B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105634008A (zh) * | 2016-01-20 | 2016-06-01 | 华南理工大学 | 考虑风速强度的分布式风力发电并网后概率线损分摊方法 |
US20170317498A1 (en) * | 2016-05-02 | 2017-11-02 | Nec Laboratories America, Inc. | Resiliency Controller for Voltage Regulation in Microgrids |
Non-Patent Citations (1)
Title |
---|
王晗等: "基于改进信赖域算法的孤岛交直流混合微电网潮流计算", 《电力***自动化》 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108964061A (zh) * | 2018-07-23 | 2018-12-07 | 长沙理工大学 | 一种考虑负荷频率电压静特性的含风电交直流电力***概率动态连续潮流新方法 |
CN108964061B (zh) * | 2018-07-23 | 2021-10-08 | 长沙理工大学 | 一种考虑负荷频率电压静特性的含风电交直流电力***概率动态连续潮流计算方法 |
CN109347110A (zh) * | 2018-10-29 | 2019-02-15 | 河海大学 | 一种含高比例风电并网的自适应线性化概率潮流计算方法 |
CN109347110B (zh) * | 2018-10-29 | 2022-01-28 | 河海大学 | 一种含高比例风电并网的自适应线性化概率潮流计算方法 |
CN109950935A (zh) * | 2018-11-08 | 2019-06-28 | 国网甘肃省电力公司电力科学研究院 | 孤岛运行的交直流混联微电网概率潮流方法 |
CN109672167A (zh) * | 2019-01-29 | 2019-04-23 | 福州大学 | 基于改进下垂算法的用于光伏低压微电网的功率控制方法 |
CN110112790A (zh) * | 2019-06-06 | 2019-08-09 | 南方电网科学研究院有限责任公司 | 有源配电网运行效率的评价方法、装置、设备及存储介质 |
CN110112790B (zh) * | 2019-06-06 | 2023-10-20 | 南方电网科学研究院有限责任公司 | 有源配电网运行效率的评价方法、装置、设备及存储介质 |
CN110276135A (zh) * | 2019-06-25 | 2019-09-24 | 华北电力大学 | 一种并网停车场的可用容量确定方法、装置及计算设备 |
CN112491056A (zh) * | 2020-11-23 | 2021-03-12 | 广东电网有限责任公司电网规划研究中心 | 一种交直流混联电网的概率潮流计算方法、装置和介质 |
CN113315132A (zh) * | 2021-06-02 | 2021-08-27 | 贵州电网有限责任公司 | 一种含下垂节点的孤岛微电网三相潮流计算方法 |
CN114936454A (zh) * | 2022-05-18 | 2022-08-23 | 东北电力大学 | 一种基于半不变量法的交直流混联***概率潮流计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107846024B (zh) | 2020-12-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107846024A (zh) | 孤岛微电网概率潮流解析计算方法 | |
Sedghi et al. | Optimal storage planning in active distribution network considering uncertainty of wind power distributed generation | |
Liang et al. | Risk-based uncertainty set optimization method for energy management of hybrid AC/DC microgrids with uncertain renewable generation | |
Eid et al. | Efficient operation of battery energy storage systems, electric-vehicle charging stations and renewable energy sources linked to distribution systems | |
Gupta | Probabilistic optimal reactive power planning with onshore and offshore wind generation, EV, and PV uncertainties | |
CN108306303A (zh) | 一种考虑负荷增长和新能源出力随机的电压稳定评估方法 | |
Roustaee et al. | Multi-objective stochastic operation of multi-microgrids constrained to system reliability and clean energy based on energy management system | |
CN104268345B (zh) | 一种支撑微电网规划设计的全过程暂态仿真方法 | |
CN107069814A (zh) | 配网分布式电源容量布点的模糊机会约束规划方法与*** | |
CN106849112A (zh) | 基于非支配邻域免疫算法的配电网多目标无功优化方法 | |
CN106786977B (zh) | 一种电动汽车充电站的充电调度方法 | |
CN109462257A (zh) | 一种计及多元随机变量电网电压稳定的灵敏度辨识方法 | |
Wang et al. | Two-stage full-data processing for microgrid planning with high penetrations of renewable energy sources | |
Gildenhuys et al. | Optimization of the operational cost and environmental impact of a multi-microgrid system | |
CN106026169A (zh) | 一种基于多微电网并入配电网的分解协调优化方法 | |
CN107122924A (zh) | 一种智能配用电***承载能力评估方法 | |
CN107039981A (zh) | 一种拟直流线性化概率最优潮流计算方法 | |
Xie et al. | Effect analysis of EV optimal charging on DG integration in distribution network | |
Dong et al. | A distributed power transfer limit calculation method for multi-area interconnection power networks | |
Wang et al. | A novel consensus-based optimal control strategy for multi-microgrid systems with battery degradation consideration | |
Zhu et al. | Dynamic multi-objective dispatch considering wind power and electric vehicles with probabilistic characteristics | |
CN111146782B (zh) | 一种主动配电网分层的时变优化追踪方法 | |
Yu et al. | Security region based security-constrained unit commitment | |
Kaur et al. | Design of the ANFIS based optimized frequency control module for an electric vehicle charging station | |
Aljohani | Intelligent Type-2 Fuzzy Logic Controller for Hybrid Microgrid Energy Management with Different Modes of EVs Integration |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |