[发明专利]一种大规模时序生产模拟加速求解方法

权利要求书

1.一种大规模时序生产模拟加速求解方法，其特征在于包括：

发电机出力转移分布因子矩阵快速计算：在进行未来某一日运行模拟中计算转移分布因子矩阵时，获取当天电力网络拓扑相对于已模拟完毕的第一日网络拓扑的变化，利用矩阵辅助求逆定理公式修正第一日的发电机出力转移分布因子矩阵进而得到当日的发电机出力转移分布因子矩阵；

机组组合模型快速求解；先求解一种场景下无电力网络潮流约束机组组合模型，然后判断支路潮流是否越界，如果越界将所述支路潮流对应的约束添加到约束集中，重新求解机组组合模型，循环该过程直到所有线路潮流均达标，以完成机组组合计算。

2.根据权利要求1所述的大规模时序生产模拟加速求解方法，其特征在于发电机出力转移分布因子矩阵快速计算包括下列步骤：

1)获取运行模拟的第一天对应的节点支路关联矩阵A以及支路导纳矩阵y，由式(1)计算发电机线路潮流转移分布因子G，

G＝yA^T(AyA^T)^-1(1)

2)在运行模拟之前，判断该日对应的电力系统网络拓扑，是否与已计算过的日期的拓扑结构相同，若相同，直接取记录中所述日期对应的G矩阵作为计算结果，若不同，则获取所述日期系统拓扑相对与第一日网络拓扑之间的差别，将新的发电机出力转移分布因子矩阵带入新的优化模型中进行求解。

3.根据权利要求2所述的大规模时序生产模拟加速求解方法，其特征在于所述步骤2)中

若当日线路新增一条支路，则获取关联矢量M_ij和支路导纳y_ij；若当日线路减少一条支路，则相当于新增一条关联矢量相同，导纳值为负的支路，新的转移分布因子矩阵由式(2)计算：

其中，G-yA^TX表示原始系统中所有支路对应的发电机转移分布因子，表示新增支路对应的转移分布因子，

其中，a为新增支路后的节点导纳矩阵，i,j分别表示行和列。

。

4.根据权利要求1所述的大规模时序生产模拟加速求解方法，其特征在于所述机组组合模型快速求解包括下列步骤：

a)在进行每一日的机组组合模型求解时，先求解该日对应的无网络约束机组组合模型,得到最优解所述组合模型参见式(5)：

式中，目标函数f(X)是系统运行总成本最小，决策变量X为日内各时段机组的出力,CX^T≤b约束表示与发电机组相关的运行约束，b表示机组的爬坡速率限制，分别为X^T的下限、上限约束；

b)获取当前发电机线路转移分布因子矩阵约束GX^T≤g，g为线路潮流上限每一行表示一条线路的潮流约束，按线路顺序检查X^T是否满足各条约束，如果不满足，则单独将该条约束添加到原模型中，重新进行求解得到更新解X^T，再次带入GX^T≤g约束中进行检验是否满足约束，循环上述过程直到所有线路潮流均不越界。