1.一种安全约束的电力系统运行模拟模型快速生成与求解方法，定义输电线路、电缆、变压器以及连接两个母线的输电设备为“支路”；定义电力系统中的所有母线为“节点”；定义K为系统支路数，N为节点系统节点数，M为系统发电机数；
其特征在于，该方法包括以下步骤：
1)获得在运行模拟时间范围内所有涉及到所有支路及节点信息，根据支路与节点的连接关系以及各支路电抗，计算电力系统考虑所有支路的原始节点阻抗矩阵、负荷转移分布因子原始矩阵以及发电机转移分布因子原始矩阵；
2)依次开始进行逐日运行模拟计算，获取当前模拟日的支路开断状态，根据支路开断状态，将负荷转移分布因子原始矩阵与发电机转移分布因子原始矩阵进行修正后得到当前模拟日考虑支路开断状态后的负荷转移分布因子矩阵与发电机转移分布因子矩阵；
3)根据无安全约束的电力系统机组组合优化模型采用混合整数规划优化求解技术进行求解得到各发电机组各时段出力，对于每一时段，根据发电机各时段出力以及节点负荷，计算各支路潮流，根据各支路潮流以及支路潮流极限判断各支路过载情况；根据过载线路对应的负荷转移分布因子以及发电机转移分布因子矩阵的行，生成机组组合模型中的安全约束，加入无安全约束的电力系统机组组合优化模型中利用混合整数规划优化求解技术再次进行求解，得到发电机各时段出力后再次判断是否有支路过载；通过多次迭代直至没有支路潮流过载，得到考虑安全约束的各机组各时段出力，完成当日运行模拟并跳转至下一日进行运行模拟，依次对待模拟日进行运行模拟得到全年考虑安全约束的运行模拟结果。
2.如权利要求1所述方法，其特征在于，所述步骤1)具体包括：
1.1)根据各支路电抗参数，形成支路导纳矩阵：
式中，y为支路导纳矩阵，xl为支路l的阻抗，l＝1、2…..K；
1.2)根据所有支路与节点连接关系，形成节点-支路关联矩阵，以及发电机-节点关联矩阵；首先，对各支路形成支路-节点关联矢量Hl，该矢量共有N个元素，对于支路l，在Hl中其起始节点i对应的元素置1，在终止节点j对应元素置-1：
利用支路-节点关联矢量形成节点-支路关联矩阵：
A＝[H1 H2 … HK]；
将系统中连接发电机组容量最大的节点选定系统的松弛节点，将松弛节点对应行删去，得到降阶节点-支路关联矩阵
然后，对各发电机建立关联矢量Rm，该矢量共有N个元素，对于机组m，m＝1、2….M，设该机组所在节点为i，在Rm中对应的第i个元素置1：
利用发电机-节点关联矢量形成发电机-节点关联矩阵：
AG＝[R1 R2 … RM]；
根据选定的系统松弛节点，将松弛节点对应行删去，得到降阶节点-支路关联矩阵
1.3)根据支路导纳矩阵以及节点-支路关联矩阵形成发电机转移分布因子原始矩阵以及负荷转移分布因子原始矩阵，根据下式计算原始节点阻抗矩阵X、负荷转移分布因子原始降阶矩阵 以及发电机转移分布因子原始降阶矩阵
对于该松弛节点 分别在 以及 对应的列插入一行元素均为0的向量，得到负荷转移分布因子原始矩阵GG以及发电机转移分布因子原始矩阵GD。
3.如权利要求1所述方法，其特征在于，所述步骤2)具体包括：
2.1)定义d为当前模拟日序号，定义当日断开的支路序号的集合为
Ω:{b1,b2,...,bB}，定义系统当前断开的支路为b1，定义b1阻抗为 b1起始节点为i，终止节点为j；
2.2)计算b1的中间变量
上式中，Xii，Xjj，Xij，Xji分别为原始节点阻抗矩阵X的对应元素；
2.3)计算考虑支路b1断开后的节点阻抗矩阵的修正矩阵X′，根据修正矩阵X′以及原始节点阻抗矩阵X计算考虑支路b1断开后的节点阻抗矩阵
其中，X的下标表示X在原始节点阻抗矩阵X中元素的行号与列号，例如：Xi1为原始节点阻抗矩阵X第i行第1列元素，Xj1为X第j行第1列元素，Xi2为X第i行第2列元素，Xj2为X第j行第2列元素，XiN为X第i行第N列元素，XiN为X第i行第N列元素；
2.4)根据负荷转移分布因子原始降阶矩阵 以及原始节点阻抗矩阵X计算考虑支路b1断开后的负荷转移分布因子的修正矩阵 根据修正矩阵 以及负荷转移分布因子原始降阶矩阵 计算考虑支路b1断开后的负荷转移分布因子降阶矩阵
其中，g1i为负荷转移分布因子原始降阶矩阵 第1行第i列元素，g1j为 第1行第j列元素，g2i为 第2行第i列元素，g2j为 第2行第j列元素，gLi为 第K行第i列元素，gLj为 第K行第j列元素；Xi1为原始节点阻抗矩阵X第i行第1列元素，Xj1为X第j行第1列元素，Xi2为X第i行第2列元素，Xj2为X第j行第2列元素，XiN为X第i行第N列元素，XiN为X第i行第N列元素；
2.5)根据考虑支路b1断开后的节点阻抗降阶矩阵计算 计算考虑支路b1断开后的发电机转移分布因子降阶矩阵
2.6)对于选定的松弛节点 分别在 以及 对应的列插入一行元素均为0的
向量，得到考虑支路b1断开后负荷转移分布因子矩阵 以及发电机转移分布因子矩阵
2.7)将 分别作为新的原始节点阻抗矩阵、节点转移分布因子原始降
阶矩阵、发电机转移分布因子原始降阶矩阵，计算下一条支路断开后系统的节点阻抗矩阵、节点转移分布因子降阶矩阵、发电机转移分布因子降阶矩阵，继续步骤2.2)～2.7)直至处理完d日所有断开的线路，得到第d日的发电机转移分布因子矩阵 负荷转移分布因子矩阵 生成第d日机组组合的安全约束：
其中，P为各类型机组各时段出力组成的向量，为机组组合优化模型的决策变量， 与分别为支路的潮流上下限，L为系统各节点负荷组成的向量。
4.如权利要求1所述方法，其特征在于，所述步骤3)具体包括：
3.1)定义系统起作用安全约束集合 表示空集，定义计算循环变量s＝0，计算第d日不包含安全约束 机组组合模型，得到各机组出力X(0)，根据X(0)
以及各节点负荷向量L，计算支路潮流：
根据F(0)、 以及 判断各支路过载情况，记录过载的支路集合为Θ(1)，若集合Θ(1)(1) (1)
为空集，则跳转至步骤3.3)，若集合Θ 不为空集，则更新起作用安全约束集合Θ＝Θ ，定义 为矩阵 的第l行元素，定义 为矩阵 的第l行元素，定义 为向量
的第l个元素，定义 为向量 的第l行元素，计算第d日包含起作用安全约束机组组合模型，起作用约束表示为：
求解包含安全约束集合Θ的机组组合模型，获得第一次迭代中考虑安全约束集合Θ(1)
的机组组合模型的最优解P ；
(s)
3.2)s＝s+1，进入下一次迭代，对于第s次迭代，根据P 以及各节点负荷向量L，计算支路潮流：
(s) (s)
根据F 、 以及 判断各支路过载情况，记录过载的支路集合为Θ ，若集合为空(s-1) (s)
集，则跳转至步骤3.3)，若集合不为空集，则更新起作用安全约束集合Θ＝Θ ∪Θ ，计算第d日包含起作用安全约束机组组合模型，起作用约束表示为：
(s+1)
求解包含安全约束集合Θ的机组组合模型，获得P ，重新实施步骤3.2)；
(s)
3.3)将P 作为考虑系统安全约束后机组组合的最优解，d＝d+1，跳转至步骤2.1)计算第d+1日的运行模拟。