面向大规模交直流系统的宽频机电-电磁混合仿真方法技术方案

本发明专利技术涉及一种面向大规模交直流系统的宽频机电‑电磁混合仿真方法，属于电力系统电磁暂态分析技术领域。该方法首先初始化各个状态变量，接着划分系统为机电子系统和电磁子系统，并分别得到电磁子系统的微分方程和机电子系统的微分方程，接着联立求解机电子系统的微分代数方程，得到机电子系统中各个状态变量的值，将新得到的状态变量带回微分方程迭代求解，直到仿真时间到达设定时间，最终得到含有电力电子单元电力系统的电磁暂态仿真结果。本发明专利技术方法不仅可以满足仿真精度和数值稳定性要求，同时极大地提高了仿真效率，利于工程推广实现。

A Broadband Electromechanical-Electromagnetic Hybrid Simulation Method for Large-Scale AC/DC Systems

【技术实现步骤摘要】
面向大规模交直流系统的宽频机电-电磁混合仿真方法
本专利技术属于电力系统电磁暂态分析
，具体涉及一种面向大规模交直流系统的宽频机电-电磁混合仿真方法。
技术介绍
目前，随着电网规模的扩大，大量以换流器为基础的电力电子装置接入交流电网，换流器与电网之间的相互作用频带越来越宽并且日益复杂化。针对电力系统的仿真，如果单纯采用电磁暂态仿真程序，则仿真规模通常会受到限制，因此电网的大部分区域需要等值简化才能进行全网络的电磁暂态计算，这将导致仿真结果无法准确刻画被等值网络的机电暂态及其电磁暂态特性；如果单纯采用机电暂态程序，则电力电子装置、非线性负荷的电磁暂态特性无法准确描述。因此如何克服现有技术的不足是目前电力系统电磁暂态分析
亟需解决的问题。
技术实现思路
本专利技术的目的是为了解决现有技术的不足，提供一种面向大规模交直流系统的宽频机电-电磁混合仿真方法，该方法基于宽频机电-电磁混合仿真的概念，将需要详细仿真的电力电子装置等根据其母线位置设置在多个不同的电磁子系统中，将其余交流电网设置在机电子系统中进行仿真，对分网后的子系统采用各自独立的电磁暂态程序以及机电暂态程序进行仿真。不同子系统之间宽频相互作用采用本专利技术提出的接口模型进行表征，不仅可以保证仿真精度，同时可以极大地提升分布式混合仿真的仿真效率。为实现上述目的，本专利技术采用的技术方案如下：面向大规模交直流系统的宽频机电-电磁混合仿真方法，其特征在于，包括如下步骤：步骤(1)，全网潮流初始化确定每个状态变量的初始值；步骤(2)，将含有电力电子单元的电力系统划分为一个机电子系统和多个电磁子系统，如图1所示。图1左侧为机电子系统，右侧为电磁子系统，并接有频率相关等值模型，中间WAN表示数据接口，图中1-N表示电磁子系统的编号；步骤(3)，计算电磁子系统的低频段接口等值电路，即戴维南等值电路；通过加流求压法得到第j个电磁子系统的戴维南等值阻抗；戴维南等值电势的计算公式如下：其中，是机电子系统的接口电压幅值和相位参数；δ表示相角，a、b、c表示电力系统三相；ω表示系统角速度；t表示仿真时间；步骤(4)，计算电磁子系统中的高频段接口等值电路，即频率相关等值电路；首先通过建立最小二乘拟合问题，确定接口频率特性拟合后的有理函数；该最小二乘拟合问题对应如下：其中，j为电磁子系统序号，1≤j≤N，s＝j2πf，Yfdne(s)为频率相关节点导纳矩阵，为拟合后频率相关节点导纳矩阵；ρj为第j个电磁子系统的权重系数；利用拟合得到的构建最小二乘问题：ΔAx＝b；其中，ξ1,…,ξl为l个频率采样点；λ1,…,λr为节点导纳矩阵Yfdne(s)已知的r个特征值；求解上述最小二乘问题得到φ1,..φr,c,从而获得频率相关等值模型；所述的频率相关等值模型表示为：其中，将上述频率相关等值模型转换为状态空间形式，如下：其中，xfdne,j为第j个电磁子系统中频率相关等值的状态变量；F3,j为第j个电磁子系统对应的频率相关等值模型转换的状态空间方程；是第j个电磁子系统三相电流；列写第j个电磁子系统对应的微分方程组：其中，为第j个电磁子系统内部的三相状态变量；F1,j、F2,j和F3,j分别对应第j个电磁子系统内部的微分方程、戴维南等值电路对应的微分方程和频率相关等值模型转换的状态空间方程；和分别指第j个电磁子系统三相电流和电压；采用数值积分的方法计算电磁子系统对应的上述微分方程组，得到电磁子系统对应的微分方程对应的电压、电流以及xfdne,j状态变量曲线；步骤(5)，计算机电子系统中的接口电压和电流并采用对称分量法将其转化到120坐标系，即对机电子系统采用诺顿等效，得到其等效导纳计算机电子系统中的等值电流如下：其中，下标t表示为接口处电流或电压，j为电磁子系统序号；该方程为机电子系统等值电流方程；步骤(6)，联立电磁子系统微分方程和机电子系统等值电流方程，求解得到机电子系统中各个节点电压和支路电流的值；步骤(7)，重复步骤(2)-步骤(6)，将新得到的状态变量带回微分方程迭代求解，直到仿真时间到达设定时间，整个过程结束，最终得到含有电力电子单元电力系统的电磁暂态仿真结果，即电压、电流曲线。进一步，优选的是，步骤(1)中，状态变量包括仿真步长、网络拓扑结构和潮流数据。进一步，优选的是，步骤(1)中，所述的初始化的方法具体为：设置电磁暂态子系统仿真步长hemt，设置机电暂态子系统仿真步长hTS，潮流计算得到各节点电压电流初始值。进一步，优选的是，电磁暂态子系统仿真步长hemt取值为10-100μs；机电暂态子系统仿真步长hTS取值为5ms。进一步，优选的是，步骤(2)的具体方法为：将整个电力系统划分为基于机电暂态稳定性的交流电网子系统和基于电磁暂态的直流电网子系统，即机电子系统和电磁子系统。进一步，优选的是，步骤(3)的具体方法为：步骤(3)的具体方法为：求解戴维南等值阻抗，是将网络中所有电源置零，接口母线注入单位电流源Ijl，求得端口电压Ujl，则等值阻抗Zjl＝Ujl/Ijl，j表示第j个电磁子系统，l＝a,b,c表示系统三相。进一步，优选的是，步骤(5)的具体方法为：机电子系统接口处三相电压电流通过傅里叶分析或曲线拟合的方式得到其基波相量值由对称分量法得到其正、负、零序分量：再经由诺顿等效，得到其等效导纳并计算机电子系统中的等值电流得到其接口电路。为了准确快速地分析电力电子装置与交流电网相互作用及其电力系统整体的动态特性(具体包括机电暂态特性和电磁暂态特性)，需要采用具有宽频特性的混合仿真方法进行对整个系统进行仿真。所以，本专利技术提出了一种面向大规模交直流系统的机电-电磁混合仿真方法。该方法需要首先初始化各个状态变量，接着划分系统为机电子系统和电磁子系统，并分别得到电磁子系统的微分方程和机电子系统的微分方程，接着联立求解机电子系统的微分代数方程，得到机电子系统中各个状态变量的值，将新得到的状态变量带回微分方程迭代求解，直到仿真时间到达设定时间，最终得到含有电力电子单元电力系统的电磁暂态仿真结果。低频段指基频附近频段(即频率小于100Hz)，高频段指远高于基频频段(即频率大于100Hz)。本专利技术与现有技术相比，其有益效果为：本专利技术方法在电力系统电磁暂态仿真计算下，由于考虑了机电子系统和电磁子系统的宽频相互作用，可以反映高频电力电子装置与大规模交流电网之间的相互作用。该方法尤其适用于含有模块化多电平模块、高压直流模块等多电力电子开关器件的电力系统电磁暂态仿真计算，具有极大的工程实用价值。本专利技术既能保证大规模交流系统和复杂的电力电子装置的仿真精度，又可以保证整体的仿真效率，满足大规模交直流系统的安全稳定分析的要求。该方法将整个系统划分为一个机电子系统和多个电磁子系统。大规模交流系统划分到机电子系统，其他电力电子装置划分到一个或者多个电磁子系统。两者之间通过频率相关等值(FDNE)模型进行交互，实现全网仿真。该方法不仅可以满足仿真精度和数值稳定性要求，同时极大地提高了仿真效率，利于工程推广实现。附图说明图1是本专利技术方法涉及的网络划分示意图。具体实施方式下面结合实施例对本专利技术作进一步的详细描述。本领域技术人员将会理解，下列实施例仅用于说明本专利技术，而不应视为限定本专利技术的范围。实施例中未注明具本文档来自技高网...

【技术保护点】
1.面向大规模交直流系统的宽频机电‑电磁混合仿真方法，其特征在于，包括如下步骤：步骤(1)，全网潮流初始化确定每个状态变量的初始值；步骤(2)，将含有电力电子单元的电力系统划分为一个机电子系统和多个电磁子系统；步骤(3)，计算电磁子系统的低频段接口等值电路，即戴维南等值电路；通过加流求压法得到第j个电磁子系统的戴维南等值阻抗；戴维南等值电势的计算公式如下：

【技术特征摘要】
1.面向大规模交直流系统的宽频机电-电磁混合仿真方法，其特征在于，包括如下步骤：步骤(1)，全网潮流初始化确定每个状态变量的初始值；步骤(2)，将含有电力电子单元的电力系统划分为一个机电子系统和多个电磁子系统；步骤(3)，计算电磁子系统的低频段接口等值电路，即戴维南等值电路；通过加流求压法得到第j个电磁子系统的戴维南等值阻抗；戴维南等值电势的计算公式如下：其中，是机电子系统的接口电压幅值和相位参数；δ表示相角，a、b、c表示电力系统三相；ω表示系统角速度；t表示仿真时间；步骤(4)，计算电磁子系统中的高频段接口等值电路，即频率相关等值电路；首先通过建立最小二乘拟合问题，确定接口频率特性拟合后的有理函数；该最小二乘拟合问题对应如下：其中，j为电磁子系统序号，1≤j≤N，s＝j2πf，Yfdne(s)为频率相关节点导纳矩阵，为拟合后频率相关节点导纳矩阵；ρj为第j个电磁子系统的权重系数；利用拟合得到的构建最小二乘问题：ΔAx＝b；其中，ξ1,…,ξl为l个频率采样点；λ1,…,λr为节点导纳矩阵Yfdne(s)已知的r个特征值；求解上述最小二乘问题得到φ1,..φr,c,从而获得频率相关等值模型；所述的频率相关等值模型表示为：其中，将上述频率相关等值模型转换为状态空间形式，如下：其中，xfdne,j为第j个电磁子系统中频率相关等值的状态变量；F3,j为第j个电磁子系统对应的频率相关等值模型转换的状态空间方程；是第j个电磁子系统三相电流；列写第j个电磁子系统对应的微分方程组：其中，为第j个电磁子系统内部的三相状态变量；F1,j、F2,j和F3,j分别对应第j个电磁子系统内部的微分方程、戴维南等值电路对应的微分方程和频率相关等值模型转换的状态空间方程；和分别指第j个电磁子系统三相电流和电压；采用数值积分的方法计算电磁子系统对应的上述微分方程组，得到电磁子系统对应的微分方程对应的电压、电流以及xfdne,j状态变量曲线；步骤(5)，计算机电子系统中的接口电压和电流并采用对称分量法将其转化到120坐标系，即对机电子系统采用诺顿等效，得到其等效导纳计算机电...

【专利技术属性】
技术研发人员：谢小荣，段荣华，舒德兀，吴琛，刘威，黄伟，程旻，
申请(专利权)人：云南电网有限责任公司，清华大学，
类型：发明
国别省市：云南,53