CN107221935A - 基于置信度变换的交直流混联输电系统区间潮流方法 - Google Patents
基于置信度变换的交直流混联输电系统区间潮流方法 Download PDFInfo
- Publication number
- CN107221935A CN107221935A CN201710437922.2A CN201710437922A CN107221935A CN 107221935 A CN107221935 A CN 107221935A CN 201710437922 A CN201710437922 A CN 201710437922A CN 107221935 A CN107221935 A CN 107221935A
- Authority
- CN
- China
- Prior art keywords
- mrow
- mtd
- msub
- mtr
- msup
- 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
-
- 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)
- Power Engineering (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
一种基于置信度变换的交直流混联输电系统区间潮流方法,包括如下步骤:(1)以交直流区间参数形式刻画系统中的不确定变量,建立交直流混联系统的区间潮流模型;(2)对置信度变换方法进行适应性推导,将区间优化模型转换为确定性数学模型,实现含区间不确定变量的约束转换;(3)采用预测校正原始对偶内点方法分别对转换后的确定性模型进行求解,通过修正方程计算仿射方向,得到仿射扰动因子后回代得到校正方向,进而得到修正量。本发明中采用的PDT方法对不确定性问题转换效率较高,可灵活方便地将区间不确定变量进行确定性转换;区间优化的结果可为决策者对交直流混联输电系统进行合理的风险控制和规划提供重要参考。
Description
技术领域
本发明涉及一种基于置信度变换的交直流混联输电系统区间潮流(IOPF)方法,属电力输送技术领域。
背景技术
高压直流输电不仅在功率控制方面具有快速性和灵活性,而且在远距离跨区域、大容量输送电能上优势明显,已成为现代电网互联的重要方式。2020年国家电网将建成“三华”特高压同步电网骨干网架,建成“四纵七横”特高压目标网架和21回特高压直流,特高压交直流混联系统在高压输电中扮演着越来越重要的角色,直流系统逐渐成为影响电力系统规划与运行的关键因素之一。
另外,由于设备故障或检修,系统参数的测量、传输等误差以及近年来新能源发电的规模性并网使得电力系统的分析与计算需要面对更多不确定性因素。为了帮助相关人员更全面地获取系统的运行信息,捕捉系统的薄弱环节及潜在危机,对交直流混联系统进行不确定性分析计算显得尤为必要。
近年来,一些学者对交直流混联系统的不确定性分析进行了诸多有益探索。比如,采用实用化概率的方法进行稳定性分析,完善了系统和节点的风险指标。但概率数学的方法必须以大量的历史数据为基础,以构造精确的概率分布函数;一种适于求解含新能源的交直流混联系统的区间最优潮流问题,并将区间问题转化为两个确定性问题求解,取得了较好的计算效果。其采用的区间数学只需区间变量的外延参数,建模方便,对历史数据的依赖度较少。PDT(Probability Degree Transformation)是一种将区间优化模型转换为确定性数学模型求解的转换方法,叫作置信度变换。
发明内容
本发明的目的是,为了解决交直流混联系统的不确定性问题,提出一种基于置信度变换的交直流混联输电系统区间潮流方法
一种基于置信度变换的交直流混联输电系统区间潮流方法,其特征在于,所述方法包括如下步骤:
(1)以交直流区间参数形式刻画系统中的不确定变量,建立交直流混联系统的区间潮流模型;
(2)对置信度变换(PDT)的方法进行适应性推导,将区间优化模型转换为确定性数学模型,实现含区间不确定变量的约束转换;
(3)采用预测校正原始对偶内点方法分别对转换后的确定性模型进行求解,通过修正方程计算仿射方向,得到仿射扰动因子后回代得到校正方向,进而得到修正量。
本发明步骤(1)中所述交直流混联系统的区间潮流模型如下:
目标函数如下式:
式中,目标函数f(x)为系统损耗;SG为发电机节点集合;PGi为节点i发电机有功功率PDi为节点i负荷的有功功率。
交流系统潮流方程等式约束如下式:
式中,Vi为节点i的电压幅值;δi为节点i的电压相角;IAC为与交流节点i所联接支路的对应节点集合;QRi为节点i无功电源无功功率;Yij为节点导纳矩阵第i行j列元素的幅值,αij为对应的相角,δij=δi-δj-αij;i,j∈SB,SB为节点集合。
直流系统潮流方程等式约束如下式:
式中,k∈Sd,Sd为直流节点集合;IDC为与直流系统节点i所联接支路对应的节点集合;Idk和Idk分别表示直流节点k的电流和电压;为直流节点k的换流器功率因数角;正、负号分别对应于整流器和逆变器。
换流器方程及直流网络方程等式约束如下式:
式中,θdk为直流节点k的换流器控制角;Xck为直流节点k的换流器换相电抗;Kγ表示由换流器换相电抗产生的换相系数;gdkj为消去联络节点后的直流网络节点电导矩阵Gd的元素;Sdk为换流站类型,“+1”表示整流测,“-1”表示逆变侧。
系统安全约束和直流变量约束如下式:
式中,SPV为PV(发电机)节点集合;SPQ为PQ节点集合;和(·)分别代表变量(·)的上、下限。
本发明步骤(1)中所述交直流系统区间参数表示如下:
与常规OPF模型不同,式(1)-(5)的交直流混联系统OPF模型考虑变量的不确定性,并以区间数的形式表示,考虑负荷功率以服从区间分布的随机变量表示如下:
式中,和分别为有功负荷变量区间分布的上下界;和分别为无功负荷变量区间分布的上下界。将式(6)代入式(1)-(5)的交直流混联系统OPF模型,即转换成交直流混联系统的IOPF模型,一个包含区间不确定变量的非线性区间优化问题。
本发明步骤(2)中所述含区间不确定变量的约束转换如下:
区间优化方法只需区间变量的外延参数,无需构造精确的概率分布模型,本发明假设区间变量[X]±=[X-,X+]服从均匀分布,x<[X±]的置信度定义为:
引入置信度λ,λ∈[0,1],则P(x<[X±])≥λ可转换为:
x≤X-λ+X+(1-λ) (8)
同理,P(x>[X±])≥λ可转换为:
x≥X+λ+X-(1-λ) (9)
当λ=0时,式(8)表达为x<X+,表示决策者只关心区间上边界,即完全乐观;当λ=1时,式(8)等效为x<X-,表示决策者更倾向于下边界取值,即完全悲观,即λ的取值体现了决策者对不确定变量的谨慎程度。由式(7)-式(9)可以看出:置信度λ是将区间变量的边界信息转换为确定性变量的比例参数,是根据具体应用场景而设定的,以实现不同电网运行场景的切换。
区间不确定变量转换方法思路为:利用区间置信度公式,使区间约束条件成立的概率满足给定的区间置信度水平,从而实现区间约束条件的确定性转化。对于一般含区间变量的非线性优化问题:
式中:[X-,X+]为不等式约束区间;[Y-,Y+]为等式约束区间。引入区间不等式约束变量的置信度λX=[λX1,λX2,...,λXm]T,m为含区间变量不等式约束数。根据式(7)-式(9)的转换方式,式(10)中的区间不等式约束可转换为:
同理,由于区间等式约束可转化为两个不等式约束,引入的置信度λY=[λY1,λY2,...,λYl]T,式(10)中的区间等式约束采用如下转换方式:
式中,n为含区间变量等式约束数。
本发明步骤(3)中所述预测校正原始对偶内点法基于原对偶内点法理论,通过修正方程计算仿射方向,得到仿射扰动因子后回代得到校正方向,进而得到修正量。
1、原对偶内点法
OPF计算模型可简化为非线性优化形式:
式中,f(x)为目标函数;h(x)为等式约束函数;g(x)为不等式约束函数,g和分别为其下、上限。
原对偶内点法求解该非线性规划问题的步骤如下:
首先,引入松弛变量,并将模型转换为拉格朗日函数形式:
式中,y,z,w均为拉格朗日乘子;μ为引入的障碍常数;l,u为不等式约束的松弛变量;r为不等式个数。
由一阶最优性条件,导出KKT方程组:
其中,L,U,Z,W为r×r全对角矩阵,由上式的最后两个方程可以解得μ=(lTz-uTw)/(2r),定义对偶间隙Gap=lTz-uTw。为提高算法收敛效率,引入中心参数σ,障碍常数μ=σGap/(2r)。
牛顿–拉夫逊法处理式(15)的非线性方程组,将其线性化,可转换为如下自治常微分方程组:
式中, 和分别为h(x)和g(x)海森矩阵的线性叠加。
2、预测校正原始对偶内点法
通过预测步求出仿射方向修正量,然后利用其估计互补方程泰勒展开式的2阶项,求出校正步修正量。算法在每次迭代中只增加1次回代运算,可有效减少迭代次数,提高收敛速度。预测校正步骤如下:
一、预测步
1)设定中心参数σ=0;
2)求解式(16),得到仿射方向Δxaff、Δyaff、Δzaff、Δlaff、Δuaff、Δwaff;
3)确定迭代步长和
4)计算仿射方向的互补间隙Gapaff;
5)根据当前点的对偶间隙与仿射方向下的对偶间隙间的关系动态确定中心参数值:
δ=(Gapaff/Gap)3 (19)
式中,对偶间隙Gap=lTz-uTw。
二、校正步
6)对互补松弛条件进行校正
7)L'x校正为:
本发明的有益效果是:本发明中采用的PDT方法对不确定性问题转换效率较高,可灵活方便地将区间不确定变量进行确定性转换;区间优化的结果可为决策者对交直流混联输电系统进行合理的风险控制和规划提供重要参考,两者结合应用对于推动交直流混联系统的发展具有重要意义。
附图说明
图1为交直流系统的示意图;
图2为不同λ条件下的计算结果;
图3为节点电压幅值计算结果;
图4为各发电机优化结果。
具体实施方式
下面结合本发明实施例,对本发明进行详细的描述。
本实施例一种基于置信度变换的交直流混联输电系统区间潮流方法,包括如下步骤:
1)初始化。读入交直流系统参数,给定原始与对偶变量初值,初始化迭代次数iter等;
2)确定发明内容中式(6)负荷功率的区间分布情况,并将其代入发明内容中式(1)-(5)的交直流混联最优潮流模型;
3)采用发明内容中式(11)-(12)对IOPF模型进行置信度变换,设定置信度值,将原区间问题进行确定性转换;
4)应用预测校正原始对偶内点法求解转换后的确定性模型;
5)判断对偶间隙Gap≤ε是否成立,若不满足,转到步骤6),若满足,输出计算结果;
6)迭代次数iter=iter+1,若iter小于最大迭代次数,转步骤4),否则终止计算,调试程序。
以下为以IEEE30系统为基础,基于置信度变换的IOPF问题的具体事例。
将6-7交流支路使用直流线路替代,形成交直流混联输电系统。换流站变压器电压比调节范围为10%,节点6为整流侧,节点7为逆变侧;整流侧与逆变侧换相电阻Xc1=Xc2=0.015;换流器功率因数取0.96。
表1和图2给出了IEEE-30节点系统负荷5种不同置信度情况下的计算结果。可以看出:负荷功率预测区间的置信度越大,负荷水平越低,则发电机总有功出力越高,系统的网损也随着置信度λ的增加而减小。
表1不同λ条件下的计算结果
置信度 | 0 | 0.25 | 0.50 | 0.75 | 1.00 |
目标函数/MW | 2.52 | 2.36 | 2.18 | 2.08 | 1.95 |
总有功出力/MW | 314.26 | 299.93 | 285.57 | 271.31 | 257.01 |
由于置信度λ是考虑了调度运行决策者对系统的不确定因素的乐观程度,即区间负荷功率的预测值。表3的结果表明,置信度转换的方法能够将不确定因素融入到本文的区间优化模型中,该方法对于交直流混联输电系统进行合理的风险控制以及规划具有重要的指导意义。
由发明内容中式(12)可知,λ=1时,hl(x)=Yl -对应于本文的负荷区间下限值,即对负荷功率预测较为悲观;同理,当λ=0时,对应于负荷区间上限值,即对负荷功率预测较为乐观。λ=1或λ=0是较为极端情况,接下来将基于一般情况(λ取0.25)分析。
表2系统IDOPF计算结果
图3和表2分别为λ=0.25时IEEE-30节点交直流混联系统各节点电压幅值和IOPF的计算结果。可以看出,IOPF得到了负荷区间不确定性影响下优化目标的区间,其给出的是系统资源协调利用的两个极端情况;结合置信度转换的方法,可根据实际负荷情况对置信度设置值进行调整,实现不同负荷条件下系统的运行情况的计算。
图4为IEEE-30节点系统各个发电机出力优化结果。可以看出,通过置信度转换后的发电机出力结果均处于区间优化法求得的区间边界内,置信度转换的求解结果考虑了决策者对系统的风险信息和实际的调控要求信息等,给出了在区间计算范围及明确的调度结果。同时,区间优化法的电压幅值区间上、下边界结果可为运行人员进行安全评估提供重要支撑。
上面对本发明的较佳实施方式作了详细说明,但是本发明并不限于上述实施方式,在本领域的普通技术人员所具备的知识范围内,还可以在不脱离本发明宗旨的前提下做出各种变化。
Claims (5)
1.一种基于置信度变换的交直流混联输电系统区间潮流方法,其特征在于,所述方法包括如下步骤:
(1)以交直流区间参数形式刻画系统中的不确定变量,建立交直流混联系统的区间潮流模型;
(2)对置信度变换方法进行适应性推导,将区间优化模型转换为确定性数学模型,实现含区间不确定变量的约束转换;
(3)采用预测校正原始对偶内点法分别对转换后的确定性模型进行求解,通过修正方程计算仿射方向,得到仿射扰动因子后回代得到校正方向,进而得到修正量。
2.根据权利要求1所述的基于置信度变换的交直流混联输电系统区间潮流方法,其特征在于,所述交直流混联系统的区间潮流模型如下:
目标函数为:
式中,目标函数f(x)为系统损耗;SG为发电机节点集合;PGi为节点i发电机有功功率PDi为节点i负荷的有功功率;
交流系统潮流方程等式约束为:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>V</mi>
<mi>i</mi>
</msub>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>&Element;</mo>
<msub>
<mi>I</mi>
<mrow>
<mi>A</mi>
<mi>C</mi>
</mrow>
</msub>
</mrow>
</munder>
<msub>
<mi>Y</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msub>
<mi>V</mi>
<mi>j</mi>
</msub>
<msub>
<mi>cos&delta;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>G</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>D</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>V</mi>
<mi>i</mi>
</msub>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>&Element;</mo>
<msub>
<mi>I</mi>
<mrow>
<mi>A</mi>
<mi>C</mi>
</mrow>
</msub>
</mrow>
</munder>
<msub>
<mi>Y</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msub>
<mi>V</mi>
<mi>j</mi>
</msub>
<msub>
<mi>sin&delta;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>Q</mi>
<mrow>
<mi>R</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>Q</mi>
<mrow>
<mi>D</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>;</mo>
</mrow>
式中,Vi为节点i的电压幅值;δi为节点i的电压相角;IAC为与交流节点i所联接支路的对应节点集合;QRi为节点i无功电源无功功率;Yij为节点导纳矩阵第i行j列元素的幅值,αij为对应的相角,δij=δi-δj-αij;i,j∈SB,SB为节点集合;
直流系统潮流方程等式约束为:
式中,k∈Sd,Sd为直流节点集合;IDC为与直流系统节点i所联接支路对应的节点集合;Idk和Idk分别表示直流节点k的电流和电压;为直流节点k的换流器功率因数角;正、负号分别对应于整流器和逆变器;
换流器方程及直流网络方程等式约束为:
式中,为直流节点k的换流器控制角;Xck为直流节点k的换流器换相电抗;Kγ表示由换流器换相电抗产生的换相系数;gdkj为消去联络节点后的直流网络节点电导矩阵Gd的元素;Sdk为换流站类型,“+1”表示整流测,“-1”表示逆变侧;
系统安全约束和直流变量约束为:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<munder>
<msub>
<mi>P</mi>
<mrow>
<mi>G</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>&OverBar;</mo>
</munder>
<mo>&le;</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>G</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>&le;</mo>
<mover>
<msub>
<mi>P</mi>
<mrow>
<mi>G</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>&OverBar;</mo>
</mover>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>&ForAll;</mo>
<mi>i</mi>
<mo>&Element;</mo>
<msub>
<mi>S</mi>
<mrow>
<mi>P</mi>
<mi>G</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<munder>
<mi>Q</mi>
<mo>&OverBar;</mo>
</munder>
<mrow>
<mi>R</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>&le;</mo>
<msub>
<mi>Q</mi>
<mrow>
<mi>R</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>&le;</mo>
<msub>
<mover>
<mi>Q</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>R</mi>
<mi>i</mi>
</mrow>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>&ForAll;</mo>
<mi>i</mi>
<mo>&Element;</mo>
<msub>
<mi>S</mi>
<mrow>
<mi>P</mi>
<mi>V</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<munder>
<msub>
<mi>V</mi>
<mi>i</mi>
</msub>
<mo>&OverBar;</mo>
</munder>
<mo>&le;</mo>
<msub>
<mi>V</mi>
<mi>i</mi>
</msub>
<mo>&le;</mo>
<mover>
<msub>
<mi>V</mi>
<mi>i</mi>
</msub>
<mo>&OverBar;</mo>
</mover>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>&ForAll;</mo>
<mi>i</mi>
<mo>&Element;</mo>
<msub>
<mi>S</mi>
<mrow>
<mi>P</mi>
<mi>Q</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<munder>
<msub>
<mi>V</mi>
<mrow>
<mi>d</mi>
<mi>k</mi>
</mrow>
</msub>
<mo>&OverBar;</mo>
</munder>
<mo>&le;</mo>
<msub>
<mi>V</mi>
<mrow>
<mi>d</mi>
<mi>k</mi>
</mrow>
</msub>
<mo>&le;</mo>
<mover>
<msub>
<mi>V</mi>
<mrow>
<mi>d</mi>
<mi>k</mi>
</mrow>
</msub>
<mo>&OverBar;</mo>
</mover>
</mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<munder>
<msub>
<mi>I</mi>
<mrow>
<mi>d</mi>
<mi>k</mi>
</mrow>
</msub>
<mo>&OverBar;</mo>
</munder>
<mo>&le;</mo>
<msub>
<mi>I</mi>
<mrow>
<mi>d</mi>
<mi>k</mi>
</mrow>
</msub>
<mo>&le;</mo>
<mover>
<msub>
<mi>I</mi>
<mrow>
<mi>d</mi>
<mi>k</mi>
</mrow>
</msub>
<mo>&OverBar;</mo>
</mover>
</mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<munder>
<msub>
<mi>K</mi>
<mrow>
<mi>d</mi>
<mi>k</mi>
</mrow>
</msub>
<mo>&OverBar;</mo>
</munder>
<mo>&le;</mo>
<msub>
<mi>K</mi>
<mrow>
<mi>d</mi>
<mi>k</mi>
</mrow>
</msub>
<mo>&le;</mo>
<mover>
<msub>
<mi>K</mi>
<mrow>
<mi>d</mi>
<mi>k</mi>
</mrow>
</msub>
<mo>&OverBar;</mo>
</mover>
</mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>;</mo>
</mrow>
式中,SPV为PV(发电机)节点集合;SPQ为PQ节点集合;和(·)分别代表变量(·)的上、下限。
3.根据权利要求1所述的基于置信度变换的交直流混联输电系统区间潮流方法,其特征在于,所述交直流系统区间参数:
与常规OPF模型不同,所述交直流混联系统的区间潮流模型考虑变量的不确定性,并以区间数的形式表示,考虑负荷功率以服从区间分布的随机变量表示如下:
式中,和分别为有功负荷变量区间分布的上下界;和分别为无功负荷变量区间分布的上下界;将上式代入交直流混联系统的区间潮流模型各式,即转换成交直流混联系统的IOPF模型,一个包含区间不确定变量的非线性区间优化问题。
4.根据权利要求1所述的基于置信度变换的交直流混联输电系统区间潮流方法,其特征在于,所述含区间不确定变量的约束转换如下:
区间优化方法只需区间变量的外延参数,无需构造精确的概率分布模型,假设区间变量[X]±=[X-,X+]服从均匀分布,x<[X±]的置信度定义为:
<mrow>
<mi>P</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo><</mo>
<mo>&lsqb;</mo>
<msup>
<mi>X</mi>
<mo>&PlusMinus;</mo>
</msup>
<mo>&rsqb;</mo>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mn>1</mn>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>x</mi>
<mo>&le;</mo>
<msup>
<mi>X</mi>
<mo>-</mo>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<msup>
<mi>X</mi>
<mo>+</mo>
</msup>
<mo>-</mo>
<mi>x</mi>
</mrow>
<mrow>
<msup>
<mi>X</mi>
<mo>+</mo>
</msup>
<mo>-</mo>
<msup>
<mi>X</mi>
<mo>-</mo>
</msup>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<msup>
<mi>X</mi>
<mo>-</mo>
</msup>
<mo>&le;</mo>
<mi>x</mi>
<mo>&le;</mo>
<msup>
<mi>X</mi>
<mo>+</mo>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>0</mn>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>x</mi>
<mo>&GreaterEqual;</mo>
<msup>
<mi>X</mi>
<mo>+</mo>
</msup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
引入置信度λ,λ∈[0,1],则P(x<[X±])≥λ可转换为:
x≤X-λ+X+(1-λ) (2)
同理,P(x>[X±])≥λ可转换为:
x≥X+λ+X-(1-λ) (3)
当λ=0时,式(2)表达为x<X+,表示决策者只关心区间上边界,即完全乐观;当λ=1时,式(2)等效为x<X-,表示决策者更倾向于下边界取值,即完全悲观,即λ的取值体现了决策者对不确定变量的谨慎程度;由式(1)-式(3)可以看出:置信度λ是将区间变量的边界信息转换为确定性变量的比例参数,是根据具体应用场景而设定的,以实现不同电网运行场景的切换;
区间不确定变量转换方法思路为:利用区间置信度公式,使区间约束条件成立的概率满足给定的区间置信度水平,从而实现区间约束条件的确定性转化;对于一般含区间变量的非线性优化问题:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>F</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>s</mi>
<mo>.</mo>
<mi>t</mi>
<mo>.</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>&le;</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow>
<mi>e</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow>
<mi>g</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>&le;</mo>
<mo>&lsqb;</mo>
<msup>
<mi>X</mi>
<mo>-</mo>
</msup>
<mo>,</mo>
<msup>
<mi>X</mi>
<mo>+</mo>
</msup>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow>
<mi>h</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mo>&lsqb;</mo>
<msup>
<mi>Y</mi>
<mo>-</mo>
</msup>
<mo>,</mo>
<msup>
<mi>Y</mi>
<mo>+</mo>
</msup>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
式中:[X-,X+]为不等式约束区间;[Y-,Y+]为等式约束区间;引入区间不等式约束变量的置信度λX=[λX1,λX2,...,λXm]T,m为含区间变量不等式约束数;
根据式(1)-式(3)的转换方式,式(4)中的区间不等式约束可转换为:
<mrow>
<msub>
<mi>g</mi>
<mi>k</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>&le;</mo>
<msubsup>
<mi>X</mi>
<mi>k</mi>
<mo>-</mo>
</msubsup>
<msub>
<mi>&lambda;</mi>
<mrow>
<mi>X</mi>
<mi>k</mi>
</mrow>
</msub>
<mo>+</mo>
<msubsup>
<mi>X</mi>
<mi>k</mi>
<mo>+</mo>
</msubsup>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>&lambda;</mi>
<mrow>
<mi>X</mi>
<mi>k</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
<mo>:</mo>
<mi>m</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
同理,由于区间等式约束可转化为两个不等式约束,引入的置信度λY=[λY1,λY2,...,λYl]T,式(4)中的区间等式约束采用如下转换方式:
<mrow>
<msub>
<mi>h</mi>
<mi>l</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msubsup>
<mi>Y</mi>
<mi>l</mi>
<mo>-</mo>
</msubsup>
<msub>
<mi>&lambda;</mi>
<mrow>
<mi>Y</mi>
<mi>l</mi>
</mrow>
</msub>
<mo>+</mo>
<msubsup>
<mi>Y</mi>
<mi>l</mi>
<mo>+</mo>
</msubsup>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>&lambda;</mi>
<mrow>
<mi>Y</mi>
<mi>l</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mi>l</mi>
<mo>=</mo>
<mn>1</mn>
<mo>:</mo>
<mi>n</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
式中,n为含区间变量等式约束数。
5.根据权利要求1所述的基于置信度变换的交直流混联输电系统区间潮流方法,其特征在于,所述预测校正原始对偶内点法通过预测步求出仿射方向修正量,然后利用其估计互补方程泰勒展开式的二阶项,求出校正步修正量;算法在每次迭代中只增加1次回代运算,有效减少迭代次数,提高收敛速度;预测校正步骤如下:
一、预测步
(1)设定中心参数σ=0;
(2)求解下式常微分方程组:
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msup>
<mi>H</mi>
<mo>&prime;</mo>
</msup>
</mtd>
<mtd>
<mrow>
<msub>
<mo>&dtri;</mo>
<mi>x</mi>
</msub>
<mi>h</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mo>&dtri;</mo>
<mi>x</mi>
</msub>
<mi>g</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mo>&dtri;</mo>
<mi>x</mi>
</msub>
<mi>g</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msubsup>
<mo>&dtri;</mo>
<mi>x</mi>
<mi>T</mi>
</msubsup>
<mi>h</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msubsup>
<mo>&dtri;</mo>
<mi>x</mi>
<mi>T</mi>
</msubsup>
<mi>g</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mi>I</mi>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msubsup>
<mo>&dtri;</mo>
<mi>x</mi>
<mi>T</mi>
</msubsup>
<mi>h</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mi>I</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mi>L</mi>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mi>Z</mi>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mi>U</mi>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mi>W</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msup>
<mi>&Delta;x</mi>
<mrow>
<mi>a</mi>
<mi>f</mi>
<mi>f</mi>
</mrow>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msup>
<mi>&Delta;y</mi>
<mrow>
<mi>a</mi>
<mi>f</mi>
<mi>f</mi>
</mrow>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msup>
<mi>&Delta;z</mi>
<mrow>
<mi>a</mi>
<mi>f</mi>
<mi>f</mi>
</mrow>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msup>
<mi>&Delta;w</mi>
<mrow>
<mi>a</mi>
<mi>f</mi>
<mi>f</mi>
</mrow>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msup>
<mi>&Delta;l</mi>
<mrow>
<mi>a</mi>
<mi>f</mi>
<mi>f</mi>
</mrow>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msup>
<mi>&Delta;u</mi>
<mrow>
<mi>a</mi>
<mi>f</mi>
<mi>f</mi>
</mrow>
</msup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>L</mi>
<mi>x</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>L</mi>
<mi>y</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>L</mi>
<mi>z</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>L</mi>
<mi>w</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>L</mi>
<mi>l</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>L</mi>
<mi>u</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>;</mo>
</mrow>
式中,f(x)为目标函数;h(x)为等式约束函数;g(x)为不等式约束函数 和分别为h(x)和g(x)海森矩阵的线性叠加;y、z、w均为拉格朗日乘子;μ为障碍常数;l、u为不等式约束的松弛变量;r为不等式个数;
Lx,Ly,Lz,Lw,Ll和Lu为由一阶最优性条件得到的KKT方程组:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>L</mi>
<mi>y</mi>
</msub>
<mo>=</mo>
<mi>h</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>L</mi>
<mi>z</mi>
</msub>
<mo>=</mo>
<mi>g</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>l</mi>
<mo>-</mo>
<munder>
<mi>g</mi>
<mo>&OverBar;</mo>
</munder>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>L</mi>
<mi>w</mi>
</msub>
<mo>=</mo>
<mi>g</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>u</mi>
<mo>-</mo>
<mover>
<mi>g</mi>
<mo>&OverBar;</mo>
</mover>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>L</mi>
<mi>l</mi>
</msub>
<mo>=</mo>
<mi>L</mi>
<mi>Z</mi>
<mi>e</mi>
<mo>-</mo>
<mi>&mu;</mi>
<mi>e</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>L</mi>
<mi>u</mi>
</msub>
<mo>=</mo>
<mi>U</mi>
<mi>W</mi>
<mi>e</mi>
<mo>+</mo>
<mi>&mu;</mi>
<mi>e</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>L</mi>
<mi>x</mi>
</msub>
<mo>=</mo>
<msub>
<mo>&dtri;</mo>
<mi>x</mi>
</msub>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msubsup>
<mo>&dtri;</mo>
<mi>x</mi>
<mi>T</mi>
</msubsup>
<mi>h</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mi>y</mi>
<mo>-</mo>
<msubsup>
<mo>&dtri;</mo>
<mi>x</mi>
<mi>T</mi>
</msubsup>
<mi>g</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>+</mo>
<mi>w</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mo>&dtri;</mo>
<mi>x</mi>
</msub>
<mi>g</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<msup>
<mi>L</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>(</mo>
<mrow>
<msubsup>
<mi>L</mi>
<mi>l</mi>
<mi>&mu;</mi>
</msubsup>
<mo>+</mo>
<msub>
<mi>ZL</mi>
<mi>z</mi>
</msub>
</mrow>
<mo>)</mo>
<mo>+</mo>
<msup>
<mi>U</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>(</mo>
<mrow>
<msubsup>
<mi>L</mi>
<mi>u</mi>
<mi>&mu;</mi>
</msubsup>
<mo>+</mo>
<msub>
<mi>WL</mi>
<mi>w</mi>
</msub>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>;</mo>
</mrow>
3
其中,L,U,Z,W为r×r全对角矩阵
由常微分方程组得到仿射方向Δxaff、Δyaff、Δzaff、Δlaff、Δuaff、Δwaff;
式中,(·)aff为(·)的仿射量;
(3)确定迭代步长和
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>&alpha;</mi>
<mi>P</mi>
<mrow>
<mi>a</mi>
<mi>f</mi>
<mi>f</mi>
</mrow>
</msubsup>
<mo>=</mo>
<mn>0.9995</mn>
<mi>min</mi>
<mo>{</mo>
<mfrac>
<mrow>
<mo>-</mo>
<msubsup>
<mi>l</mi>
<mi>i</mi>
<mrow>
<mi>a</mi>
<mi>f</mi>
<mi>f</mi>
</mrow>
</msubsup>
</mrow>
<mrow>
<msubsup>
<mi>&Delta;l</mi>
<mi>i</mi>
<mrow>
<mi>a</mi>
<mi>f</mi>
<mi>f</mi>
</mrow>
</msubsup>
</mrow>
</mfrac>
<mo>,</mo>
<msubsup>
<mi>&Delta;l</mi>
<mi>i</mi>
<mrow>
<mi>a</mi>
<mi>f</mi>
<mi>f</mi>
</mrow>
</msubsup>
<mo><</mo>
<mn>0</mn>
<mo>;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<mo>-</mo>
<msubsup>
<mi>u</mi>
<mi>i</mi>
<mrow>
<mi>a</mi>
<mi>f</mi>
<mi>f</mi>
</mrow>
</msubsup>
</mrow>
<mrow>
<msubsup>
<mi>&Delta;u</mi>
<mi>i</mi>
<mrow>
<mi>a</mi>
<mi>f</mi>
<mi>f</mi>
</mrow>
</msubsup>
</mrow>
</mfrac>
<mo>,</mo>
<msubsup>
<mi>&Delta;u</mi>
<mi>i</mi>
<mrow>
<mi>a</mi>
<mi>f</mi>
<mi>f</mi>
</mrow>
</msubsup>
<mo><</mo>
<mn>0</mn>
<mo>;</mo>
<mn>1</mn>
<mo>}</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>&alpha;</mi>
<mi>D</mi>
<mrow>
<mi>a</mi>
<mi>f</mi>
<mi>f</mi>
</mrow>
</msubsup>
<mo>=</mo>
<mn>0.9995</mn>
<mi>min</mi>
<mo>{</mo>
<mfrac>
<mrow>
<mo>-</mo>
<msubsup>
<mi>z</mi>
<mi>i</mi>
<mrow>
<mi>a</mi>
<mi>f</mi>
<mi>f</mi>
</mrow>
</msubsup>
</mrow>
<mrow>
<msubsup>
<mi>&Delta;z</mi>
<mi>i</mi>
<mrow>
<mi>a</mi>
<mi>f</mi>
<mi>f</mi>
</mrow>
</msubsup>
</mrow>
</mfrac>
<mo>,</mo>
<msubsup>
<mi>&Delta;z</mi>
<mi>i</mi>
<mrow>
<mi>a</mi>
<mi>f</mi>
<mi>f</mi>
</mrow>
</msubsup>
<mo><</mo>
<mn>0</mn>
<mo>;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<mo>-</mo>
<msubsup>
<mi>w</mi>
<mi>i</mi>
<mrow>
<mi>a</mi>
<mi>f</mi>
<mi>f</mi>
</mrow>
</msubsup>
</mrow>
<mrow>
<msubsup>
<mi>&Delta;w</mi>
<mi>i</mi>
<mrow>
<mi>a</mi>
<mi>f</mi>
<mi>f</mi>
</mrow>
</msubsup>
</mrow>
</mfrac>
<mo>,</mo>
<msubsup>
<mi>&Delta;w</mi>
<mi>i</mi>
<mrow>
<mi>a</mi>
<mi>f</mi>
<mi>f</mi>
</mrow>
</msubsup>
<mo><</mo>
<mn>0</mn>
<mo>;</mo>
<mn>1</mn>
<mo>}</mo>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>;</mo>
</mrow>
(4)计算仿射方向的互补间隙Gapaff:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msup>
<mi>Gap</mi>
<mrow>
<mi>a</mi>
<mi>f</mi>
<mi>f</mi>
</mrow>
</msup>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>+</mo>
<msubsup>
<mi>&alpha;</mi>
<mi>P</mi>
<mrow>
<mi>a</mi>
<mi>f</mi>
<mi>f</mi>
</mrow>
</msubsup>
<msup>
<mi>&Delta;l</mi>
<mrow>
<mi>a</mi>
<mi>f</mi>
<mi>f</mi>
</mrow>
</msup>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>+</mo>
<msubsup>
<mi>&alpha;</mi>
<mi>D</mi>
<mrow>
<mi>a</mi>
<mi>f</mi>
<mi>f</mi>
</mrow>
</msubsup>
<msup>
<mi>&Delta;z</mi>
<mrow>
<mi>a</mi>
<mi>f</mi>
<mi>f</mi>
</mrow>
</msup>
<mo>)</mo>
</mrow>
<mo>-</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>(</mo>
<mi>u</mi>
<mo>+</mo>
<msubsup>
<mi>&alpha;</mi>
<mi>P</mi>
<mrow>
<mi>a</mi>
<mi>f</mi>
<mi>f</mi>
</mrow>
</msubsup>
<msup>
<mi>&Delta;u</mi>
<mrow>
<mi>a</mi>
<mi>f</mi>
<mi>f</mi>
</mrow>
</msup>
<mo>)</mo>
<mo>(</mo>
<mi>w</mi>
<mo>+</mo>
<msubsup>
<mi>&alpha;</mi>
<mi>D</mi>
<mrow>
<mi>a</mi>
<mi>f</mi>
<mi>f</mi>
</mrow>
</msubsup>
<msup>
<mi>&Delta;w</mi>
<mrow>
<mi>a</mi>
<mi>f</mi>
<mi>f</mi>
</mrow>
</msup>
<mo>)</mo>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>;</mo>
</mrow>
(5)根据当前点的对偶间隙与仿射方向下的对偶间隙间的关系动态确定中心参数值:δ=(Gapaff/Gap)3;
式中,对偶间隙Gap=lTz-uTw;
二、校正步
(6)对互补松弛条件进行校正
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>Z</mi>
<mi>&Delta;</mi>
<mi>l</mi>
<mo>+</mo>
<mi>L</mi>
<mi>&Delta;</mi>
<mi>z</mi>
<mo>=</mo>
<mo>-</mo>
<msubsup>
<mi>L</mi>
<mi>l</mi>
<mi>&mu;</mi>
</msubsup>
<mo>-</mo>
<msup>
<mi>&Delta;Z</mi>
<mrow>
<mi>a</mi>
<mi>f</mi>
<mi>f</mi>
</mrow>
</msup>
<msup>
<mi>&Delta;l</mi>
<mrow>
<mi>a</mi>
<mi>f</mi>
<mi>f</mi>
</mrow>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>W</mi>
<mi>&Delta;</mi>
<mi>u</mi>
<mo>+</mo>
<mi>U</mi>
<mi>&Delta;</mi>
<mi>z</mi>
<mo>=</mo>
<mo>-</mo>
<msubsup>
<mi>L</mi>
<mi>u</mi>
<mi>&mu;</mi>
</msubsup>
<mo>-</mo>
<msup>
<mi>&Delta;W</mi>
<mrow>
<mi>a</mi>
<mi>f</mi>
<mi>f</mi>
</mrow>
</msup>
<msup>
<mi>&Delta;u</mi>
<mrow>
<mi>a</mi>
<mi>f</mi>
<mi>f</mi>
</mrow>
</msup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>;</mo>
</mrow>
(7)L'x校正为:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>L</mi>
<mi>x</mi>
<mrow>
<mo>&prime;</mo>
<mo>&prime;</mo>
</mrow>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>L</mi>
<mi>x</mi>
<mo>&prime;</mo>
</msubsup>
<mo>+</mo>
<msub>
<mo>&dtri;</mo>
<mi>x</mi>
</msub>
<mi>g</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>(</mo>
<msup>
<mi>L</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<msup>
<mi>&Delta;Z</mi>
<mrow>
<mi>a</mi>
<mi>f</mi>
<mi>f</mi>
</mrow>
</msup>
<msup>
<mi>&Delta;l</mi>
<mrow>
<mi>a</mi>
<mi>f</mi>
<mi>f</mi>
</mrow>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msup>
<mi>U</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<msup>
<mi>&Delta;W</mi>
<mrow>
<mi>a</mi>
<mi>f</mi>
<mi>f</mi>
</mrow>
</msup>
<msup>
<mi>&Delta;u</mi>
<mrow>
<mi>a</mi>
<mi>f</mi>
<mi>f</mi>
</mrow>
</msup>
<mo>)</mo>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>.</mo>
</mrow>
4
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710437922.2A CN107221935B (zh) | 2017-06-12 | 2017-06-12 | 基于置信度变换的交直流混联输电系统区间潮流方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710437922.2A CN107221935B (zh) | 2017-06-12 | 2017-06-12 | 基于置信度变换的交直流混联输电系统区间潮流方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107221935A true CN107221935A (zh) | 2017-09-29 |
CN107221935B CN107221935B (zh) | 2019-12-31 |
Family
ID=59948370
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710437922.2A Active CN107221935B (zh) | 2017-06-12 | 2017-06-12 | 基于置信度变换的交直流混联输电系统区间潮流方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107221935B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110097263A (zh) * | 2019-04-18 | 2019-08-06 | 新奥数能科技有限公司 | 综合能源系统的设备调控方法及装置 |
CN111799799A (zh) * | 2020-07-13 | 2020-10-20 | 福州大学 | 一种基于区间泰勒展开法的交直流混合配电网区间潮流计算方法 |
CN112039058A (zh) * | 2020-08-18 | 2020-12-04 | 国网山东省电力公司临沂供电公司 | 基于风功率预测区间的机组组合方法、系统、介质及设备 |
CN115034675A (zh) * | 2022-06-28 | 2022-09-09 | 贵州电网有限责任公司 | 一种基于全纯函数嵌入的电力系统直流最优潮流评估方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103795091A (zh) * | 2014-01-14 | 2014-05-14 | 河海大学 | 基于预测校正内点法的混合直流输电系统最优潮流方法 |
CN105552915A (zh) * | 2016-02-25 | 2016-05-04 | 国网福建省电力有限公司 | 一种考虑静态安全约束的含mmc-hvdc交直流混合电网优化潮流计算方法 |
-
2017
- 2017-06-12 CN CN201710437922.2A patent/CN107221935B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103795091A (zh) * | 2014-01-14 | 2014-05-14 | 河海大学 | 基于预测校正内点法的混合直流输电系统最优潮流方法 |
CN105552915A (zh) * | 2016-02-25 | 2016-05-04 | 国网福建省电力有限公司 | 一种考虑静态安全约束的含mmc-hvdc交直流混合电网优化潮流计算方法 |
Non-Patent Citations (3)
Title |
---|
JIANG C.等: "A nonlinear interval number programming method for uncertain optimization problems", 《EUROPEAN JOURNAL OF OPERATIONAL RESEARCH》 * |
卫志农 等: "含VSC-HVDC 的交直流系统内点法最优潮流计算", 《中国电机工程学报》 * |
鲍海波 等: "考虑新能源的交直流系统区间最优潮流", 《中国电机工程学报》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110097263A (zh) * | 2019-04-18 | 2019-08-06 | 新奥数能科技有限公司 | 综合能源系统的设备调控方法及装置 |
CN111799799A (zh) * | 2020-07-13 | 2020-10-20 | 福州大学 | 一种基于区间泰勒展开法的交直流混合配电网区间潮流计算方法 |
CN111799799B (zh) * | 2020-07-13 | 2022-03-08 | 福州大学 | 一种基于区间泰勒展开法的交直流混合配电网区间潮流计算方法 |
CN112039058A (zh) * | 2020-08-18 | 2020-12-04 | 国网山东省电力公司临沂供电公司 | 基于风功率预测区间的机组组合方法、系统、介质及设备 |
CN115034675A (zh) * | 2022-06-28 | 2022-09-09 | 贵州电网有限责任公司 | 一种基于全纯函数嵌入的电力系统直流最优潮流评估方法 |
CN115034675B (zh) * | 2022-06-28 | 2024-07-09 | 贵州电网有限责任公司 | 一种基于全纯函数嵌入的电力系统直流最优潮流评估方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107221935B (zh) | 2019-12-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101860037B (zh) | 一种高压直流输电系统网侧谐波电流的确定方法 | |
CN107221935A (zh) | 基于置信度变换的交直流混联输电系统区间潮流方法 | |
CN105633948B (zh) | 一种分布式能源接入电力系统随机模糊潮流算法 | |
Wu et al. | Coordinated optimal power flow for integrated active distribution network and virtual power plants using decentralized algorithm | |
CN106503839A (zh) | 一种海上风电场环形集电网络分层规划方法 | |
CN105846456A (zh) | 一种交直流互联电网风、火协调动态经济调度优化方法 | |
CN107134799A (zh) | 一种基于短路比的特高压电网换流器无功配置方法 | |
CN106505633B (zh) | 一种风光接入容量确定方法及装置 | |
CN104915724B (zh) | 交直流并联输电通道功率优化分配方法和系统 | |
CN111900732A (zh) | 考虑风电不确定性及无功潮流与电压约束的风电集群接入方式规划方法 | |
CN105811407A (zh) | 一种基于分布式牛顿法的微电网一次调频控制方法 | |
Aziz et al. | Admm-based distributed optimization of hybrid mtdc-ac grid for determining smooth operation point | |
CN105224760A (zh) | 一种基于风电场的vsc-hvdc并网系统可靠性计算方法 | |
CN105552915A (zh) | 一种考虑静态安全约束的含mmc-hvdc交直流混合电网优化潮流计算方法 | |
CN109598377B (zh) | 一种基于故障约束的交直流混合配电网鲁棒规划方法 | |
CN109274116A (zh) | 一种交直流混合微电网连续潮流计算方法 | |
CN105958496A (zh) | 基于连续潮流的含vsc系统最大输电能力计算方法 | |
CN105048468A (zh) | 基于分布式计算的输配电网一体化电压稳定评估方法 | |
CN106936152A (zh) | 考虑换流站损耗特性的交直流系统电压无功协调控制方法 | |
Zhai et al. | Distributed optimal power flow for VSC-MTDC meshed AC/DC grids using ALADIN | |
Gan et al. | Coordinated planning of large-scale wind farm integration system and transmission network | |
CN106208099A (zh) | 一种基于二层规划的电力系统无功优化方法及其应用 | |
CN110994697B (zh) | 含光储联合体的交直流配电网优化运行控制方法及系统 | |
Li et al. | Decentralized distributionally robust coordinated dispatch of multiarea power systems considering voltage security | |
CN113158573A (zh) | 一种基于深度学习的小水电群区域配网最优潮流获取方法 |
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 |