基于分散松弛因子潮流模型的最小发电费用增量获取方法技术

本发明专利技术公开了一种基于分散松弛因子潮流模型的最小发电费用增量获取方法。通过解决只含有一个平衡节点潮流模型的最小发电费用增量模型时，系统总的需求单位变化固定在了平衡节点，所得到的λ，即系统总的需求单位变化时，发电费用最小增量，往往不能反应系统的实际运行。本发明专利技术用原对偶内点法求解提出的模型，根据分散松弛因子的选择，系统的λ是可控的。

【技术实现步骤摘要】

【技术保护点】
一种基于分散松弛因子潮流模型的最小发电费用增量获取方法，其特征在于：包括以下步骤：步骤1：建立分散松弛因子潮流计算模型：式中，x＝[Pτ,U,δ]，fp(x)是有功残差，fq(x)为无功残差，Δδ、ΔU、ΔPτ分别是电压相角、电压幅值、有功不平衡量的修正量，上式中左侧为雅克比矩阵；步骤2：建立分散松弛因子潮流模型的最小发电费用增量计算模型：obj. min.C(Pg)s.t. Pτ＝0Pg‾≤Pg≤Pg‾]]>式中，C(Pg)为目标函数，Pτ＝0为等式约束条件，Pτ表示有功不平衡量；为不等式约束条件，Pg表示发电机的有功出力值，和分别表示发电机的有功出力值的上限和下限；步骤3：获取电力系统的网络参数；步骤4：根据步骤2中建立的含混合直流输电的电力系统最优潮流模型，构造拉格朗日函数如下：L=C(Pg)-yPτ-zT[Pg-l-Pg‾]-wT[Pg+u-Pg‾]-μΣj=1rln(lj)-μΣj=1rln(uj)]]>其中y＝[y1,…,ym]T为等式约束的拉格朗日乘子，z＝[z1,…,zr]T、w＝[w1,…,wr]T为不等式约束的拉格朗日乘子，l＝[l1,…,lr]T、u＝[u1,…,ur]T为不等式约束的松弛变量，μ是扰动因子；步骤5：程序初始化，设置状态量设置初值、拉格朗日乘子初值和罚因子初值、节点次序优化、形成节点导纳矩阵、恢复迭代计数器k'＝1、设置精度要求和最大迭代次数Kmax；步骤6：定义对偶间隙CGap＝lTz‑uTw，计算出CGap的值并判断CGap的值是否小于步骤4中设定的精度要求ε，若小于，则输出计算结果并停止执行后续步骤，若不小于，则继续执行步骤7；步骤7：根据公式μ＝σCGap/2r计算扰动因子μ；步骤8：将Pg代入到步骤1中的模型求解得到电压U、相角δ，并更新Pg步骤9：根据以下方程求解ΔPg、Δy、Δz、Δl、Δu、Δw：H′▿PgPτ▿PgTPτ0ΔPgΔy=LPg′0]]>IL-1Z0IΔzΔl=-L-1LlμLz+▿PgTPgΔPg]]>IU-1W0IΔwΔu=-U-1Luμ-Lw-▿PgTPgΔPg]]>其中：ΔPg、Δy、Δz、Δl、Δu、Δw为Pg、y、z、l、u、w的修正量；步骤10：确定原始变量和对偶变量的迭代步长：αp=0.9995min{min(-lr′Δlr′,Δlr′<0;-ur′Δur′,Δur′<0),1}]]>αd=0.9995min{min(-zr′Δzr′,Δzr′<0;-wr′Δwr′,Δwr′<0),1}]]>步骤11：更新原始变量及拉格朗日乘子；步骤12：判断迭代次数是否大于Kmax，若大于，则退出程序并输出计算不收敛的结果，若不大于，则置迭代次数k'值加1，返回步骤6。...

【技术特征摘要】

【专利技术属性】
技术研发人员：卫志农，黄为民，孙国强，孙永辉，向育鹏，
申请(专利权)人：河海大学，
类型：发明
国别省市：江苏;32