CN115021256A - 一种大规模交直流输电系统电磁暂态模型自动生成方法 - Google Patents

一种大规模交直流输电系统电磁暂态模型自动生成方法 Download PDF

Info

Publication number
CN115021256A
CN115021256A CN202210942602.3A CN202210942602A CN115021256A CN 115021256 A CN115021256 A CN 115021256A CN 202210942602 A CN202210942602 A CN 202210942602A CN 115021256 A CN115021256 A CN 115021256A
Authority
CN
China
Prior art keywords
class
power
transient
deriving
electromagnetic transient
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
Application number
CN202210942602.3A
Other languages
English (en)
Other versions
CN115021256B (zh
Inventor
刘天琪
何海林
王顺亮
史华勃
马俊鹏
王俊磊
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Sichuan University
Original Assignee
Sichuan University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Sichuan University filed Critical Sichuan University
Priority to CN202210942602.3A priority Critical patent/CN115021256B/zh
Publication of CN115021256A publication Critical patent/CN115021256A/zh
Application granted granted Critical
Publication of CN115021256B publication Critical patent/CN115021256B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J5/00Circuit arrangements for transfer of electric power between ac networks and dc networks
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00Details relating to the application field
    • G06F2113/04Power grid distribution networks
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/10Power transmission or distribution systems management focussing at grid-level, e.g. load flow analysis, node profile computation, meshed network optimisation, active network management or spinning reserve management
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E60/00Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
    • Y02E60/60Arrangements for transfer of electric power between AC networks or generators via a high voltage DC link [HVCD]

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明涉及电力系统仿真技术领域,公开了一种大规模交直流输电系统电磁暂态模型自动生成方法,基于大规模交直流输电系统的机电暂态程序和默认的电磁暂态数据,将包含换流器在内的多种电力设备及其拓扑关系抽象为节点和支路等父类,并派生得到电源、负荷、线路、变压器、换流器、控制器等六个子类,实例化后得到初始的电磁暂态模型,再考虑实际研究对象,采用交直流输电线路的分布参数、变压器的磁路参数和换流器的详细参数等电磁暂态数据对初始的电磁暂态模型进行修正,最终自动生成大规模交直流输电系统的电磁暂态模型,可用于交直流输电线路过电压、变压器磁饱和、交直流谐波耦合、宽频振荡等各种涉及电磁暂态特性问题的研究。

Description

一种大规模交直流输电系统电磁暂态模型自动生成方法
技术领域
本发明涉及电力系统仿真技术领域,具体为一种大规模交直流输电系统电磁暂态模型自动生成方法。
背景技术
“双碳”目标下,以风电、光伏发电等新能源发电为主体的新型电力系统得到了广泛发展,电力系统高度电力电子化和扁平化的特点愈加明显。电力电子设备的大规模接入使得电力系统的动态特性日益复杂,而传统机电暂态分析程序的仿真步长为毫秒级,主要用于潮流计算、短路计算和静态稳定性分析,主要反映系统的工频信息,难以用于采用分布参数的输电线路过电压、交直流谐波耦合与宽频振荡等问题的计算与分析。因此,基于已有的机电暂态数据,或采用机电暂态-电磁暂态混合仿真技术,或采用全电磁暂态仿真技术成为分析电力系统电磁暂态特性的两条主要技术路线:
1)机电暂态-电磁暂态混合仿真技术:
机电暂态-电磁暂态混合仿真技术结合了机电暂态分析程序仿真速度快和电磁暂态分析程序精度高的特点,针对大规模交直流输电系统机电暂态分析过程中新能源厂站、换流站或其他重点关注的区域开展电磁暂态建模分析。该技术主要存在的问题是计算精度严重依赖于机电-电磁暂态接口位置的选择和数据交互算法的设计;机电暂态仿真部分的精度难以满足输电线路过电压、谐波耦合和宽频振荡计算分析等仿真精度要求较高场合的需求,从而导致整个模型的仿真精度受到限制,是一个机电暂态模型向全电磁暂态模型发展的过渡产物;
2)电磁暂态仿真技术:
电磁暂态仿真技术基于交直流输电系统的电压、电流瞬时值,并充分考虑电力电子设备的非线性特征和输电线路的分布参数特性,建立考虑高频动态特性的详细电磁暂态模型。针对现有交直流输电系统的主要参数一般以机电暂态数据格式存在的现状,现有技术多在机电暂态数据的基础上,采用人工或者自动的方式直接在电磁暂态仿真平台搭建仅含工频电路参数的电磁暂态模型。
现有技术方案:
1)措施1:机电暂态-电磁暂态混合仿真技术,可以参考以下文献:
【参考文献1】刘文焯,侯俊贤,汤涌等;考虑不对称故障的机电暂态-电磁暂态混合仿真方法[J],中国电机工程学报,2010.30(13):8-15;
【参考文献2】岳程燕,田芳,周孝信等;电力系统电磁暂态-机电暂态混合仿真接口原理[J],电网技术,2006(01):23-27+88;
【参考文献3】柳勇军,闵勇,梁旭;电力系统数字混合仿真技术综述[J],电网技术,2006(13):38-43。
但该措施计算精度严重依赖于机电-电磁暂态接口位置的选择和数据交互算法的设计;机电暂态模型部分的精度难以满足线路过电压、谐波耦合和宽频振荡计算分析等仿真精度要求较高场合的需求;计算过程中初始化过程繁琐;硬件配置要求高。
2)措施2:电磁暂态仿真技术,可以参考以下文献:
【参考文献4】张民,贺仁睦,孙哲,陶华,石岩;一种机电暂态向电磁暂态的模型转换方案及实现方法[P],北京市:CN101539963B,2011-02-02;
【参考文献5】郭强,朱艺颖,胡涛等;一种自动生成大规模电网电磁暂态仿真模型的方法及系统[P],北京市:CN109766586A,2019-05-17。
但该措施电力设备类型和数目的更新问题未得到有效解决,维护开销较高;仅含系统的工频电路参数信息,没有考虑电力电子设备的非线性特征和输电线路的分布参数特性,难以准确反映交直流输电系统的电磁暂态特性,不能满足操作过电压、谐波耦合和宽频振荡等计算分析需求。
发明内容
针对上述问题,本发明的目的在于提供一种大规模交直流输电系统电磁暂态模型自动生成方法,可用于交直流输电线路过电压、变压器磁饱和、交直流谐波耦合、宽频振荡等各种涉及电磁暂态特性问题的研究。技术方案如下:一种大规模交直流输电系统电磁暂态模型自动生成方法,包括如下步骤:
步骤1:分析机电暂态程序和默认的电磁暂态数据,抽象得到关注区域交直流输电系统的UML类图;并定义顶级父类及其子类,令其具有读数据、写数据、设置坐标的静态方法;
步骤2:依据电力设备节点的数目,抽象得到包括节点类和支路类的次顶级父类;
步骤3:基于节点类,派生得到电源类和负荷类;
步骤4:基于支路类,派生得到交直流输电线路类、变压器类、换流器类和控制器类;至此得到的初始的电磁暂态模型;
步骤5:校核潮流并计算潮流误差,对比机电暂态程序和所述电磁暂态模型的程序,分别计算由机电暂态程序计算得到的有功功率、无功功率和母线电压有效值,并计算N个测量位置的潮流误差;若潮流误差的无穷范数小于误差上限值,则所述电磁暂态模型的程序正确;反之,则排查所述电磁暂态模型的程序,直至满足条件;
步骤6:基于补充的电磁暂态数据,对交直流输电系统的元件类型和参数进行补充和修正;基于包括交直流输电线路的分布参数、换流器的详细控制参数,变压器的磁路参数的详细电磁暂态数据,根据研究对象的不同,选择一项或多项数据来修正所述初始的电磁暂态模型,其余参数仍采用默认的电磁暂态数据;得到最终适用于交直流输电系统的电磁暂态模型。
进一步的,所述步骤2具体为:将节点数目为1的电力设备抽象为节点类,并增加节点名称和节点位置属性;将节点数目大于1的电力设备抽象为支路类,并增加支路名称和支路位置属性。
更进一步的,所述步骤3具体为:
步骤3.1.1:基于节点类,派生得到电源类,并增加包括额定电压、额定频率、额定容量的属性;
步骤3.1.2:基于电源类,派生得到三相电压源,并增加包括阻抗类型、阻抗数值、中心点接地选择、相位大小的电压源属性;
步骤3.1.3:基于电源类,派生得到同步发电机,并增加包括定子电阻数值、直轴不饱和同步电抗、交轴不饱和同步电抗、直轴暂态电抗、交轴暂态电抗、直轴暂态开路时间常数、交轴暂态开路时间常数、直轴次暂态电抗、交轴次暂态电抗、直轴次暂态时间常数、交轴次暂态时间常数的属性;
步骤3.2.1:基于节点类,派生得到负荷类,并增加包括额定容量和额定电压的属性;
步骤3.2.2:基于负荷类,派生得到固定负荷;固定负荷拥有一个设置功率的实例方法SetLoad(double P 0、double Q 0、double V N、double NP、double NQ、double k PF 、double k QF ),满足:
Figure 505793DEST_PATH_IMAGE001
式中,PQ分别为固定负荷吸收的有功功率和无功功率,P 0Q 0分别为额定有功功率和额定无功功率,VV N分别为相电压有效值的实际测量值和额定值,NPNQ分别为有功的电压系数和无功的电压变化系数,Δf为系统的频率变化量,k PFk QF分别为有功和无功的频率变化系数;
通过所述实例方法修改P 0Q 0V NNPNQk PFk QF系数来模拟不同类型的固定负荷;
步骤3.2.3:基于负荷类,派生得到异步电动机负荷,增加包括额定转速、机械阻尼系数、定子电阻、定子电抗、转子电阻、转子电抗、激磁电抗的属性。
更进一步的,所述步骤4具体为:
步骤4.1:基于支路类,派生得到交直流输电线路,增加包括额定频率、正序阻抗、零序阻抗、线路长度的属性;
步骤4.2:基于支路类,派生得到变压器类,增加包括额定频率、额定容量、额定绕组电压、铜耗、铁耗、激磁损耗的属性;
步骤4.3:基于支路类,派生得到换流器类,增加包括换流器个数、每个换流器压降数、桥电流额定值发属性;
步骤4.4:基于支路类,派生得到控制器类,增加包括输入信号、输出信号的属性,以及一个嵌套类ControlSet;
步骤4.5:基于控制器类,派生得到同步发电机的控制器,并设置输入信号、输出信号,以及控制参数;
步骤4.6:基于控制器类,派生得到异步电动机的控制器,并设置它的输入信号、输出信号,以及控制参数;
步骤4.7:基于控制器类,派生得到换流器的控制装置,并设置它的输入信号、输出信号,以及控制参数。
更进一步的,所述步骤5中计算N个测量位置的潮流误差,具体为计算有功功率、无功功率和母线电压有效值的误差系数:
Figure 907956DEST_PATH_IMAGE002
式中,
Figure 232627DEST_PATH_IMAGE003
Figure 757149DEST_PATH_IMAGE004
Figure 918003DEST_PATH_IMAGE005
分别为有功误差、无功误差和母线电压有效值误差;P 1i Q 1i V 1i 分别为由机电暂态程序计算得到的有功功率、无功功率和母线电压有效值;P 2i Q 2i V 2i 分别为由所述电磁暂态模型计算得到的有功功率、无功功率和母线电压有效值;S N 为交直 流输电系统的基准容量,V Ni 为第i个测量位置的额定电压;
若有功功率、无功功率和母线电压有效值的误差系数满足:
Figure 186173DEST_PATH_IMAGE006
则所述电磁暂态模型的程序正确。
与现有技术相比,本发明的有益效果是:本发明提出一种基于面向对象编程的大规模交直流输电系统电磁暂态模型自动生成方法,基于大规模交直流输电系统的机电暂态程序和默认的电磁暂态数据,将包含换流器在内的多种电力设备及其拓扑关系抽象为节点和支路等父类,并派生得到电源、负荷、线路、变压器、换流器、控制器等六个子类,实例化后得到初始的电磁暂态模型,再考虑实际研究对象,采用交直流输电线路的分布参数、变压器的磁路参数和换流器的详细参数等电磁暂态数据对初始的电磁暂态模型进行修正,最终自动生成大规模交直流输电系统的电磁暂态模型,可用于交直流输电线路过电压、变压器磁饱和、交直流谐波耦合、宽频振荡等各种涉及电磁暂态特性问题的研究。
附图说明
图1为118节点交直流输电系统的UML类图。
图2为118节点交直流输电系统电磁暂态模型的整体结构。
图3为系统中的某一组换流器。
图4为系统中的某一控制器。
图5为采用分布参数的某一交直流输电线路。
具体实施方式
下面结合附图和具体实施例对本发明做进一步详细说明。
以118节点交直流输电系统的宽频振荡问题为研究对象,基于该系统的机电暂态程序、换流器的详细参数和输电线路的分布参数,给出该区域交直流输电系统电磁暂态模型自动生成过程。
步骤1、分析机电暂态程序和默认的电磁暂态数据,抽象得到关注区域交直流输电系统的UML类图,如图1所示。图1表示系统中存在的电力设备类及其继承关系,矩形方框表示电力设备类,带箭头的实线表示电力设备类之间的继承关系,其方向为由子类指向父类。
步骤2、命名顶级父类BaseClass,并使其具有ReadData( ),WriteData( )和SetLocation( )等三个静态方法。
步骤3、依据电力设备节点的数目,基于BaseClass派生得到节点类和支路类。
步骤3具体为:
步骤3-1、将节点数目为1的电力设备抽象为节点类,并增加NodeName、NodeLocation等属性;
步骤3-2、将节点数目大于1的电力设备抽象为支路类,并增加BranchName、BranchLocation等属性。
步骤4、基于节点类,派生得到电源类和负荷类。
步骤4-1-1、基于节点类,派生得到电源类,增加BasekV、BaseFreq、S N等属性,以及电源类型选择方法ChooseSource( )。
步骤4-1-2、基于电源类,派生得到三相电压源,增加Config、SignPara、Imped、Moni等结构体来分类存储三相电压源的属性。实例化后,利用机电暂态程序的数据和默认的电磁暂态数据对各个电压源赋值。
步骤4-1-3、基于电源类,派生得到同步发电机,增加Config1、Config2、Interface、Intial、BaseData、EquiData、GeneData、SatuCurve、IntialCon1、IntialCon2、IntialCon3、OutName、OutCtr、Moni等来分类存储同步电动机的属性。实例化后,利用机电暂态程序的数据对各台同步电动机的定子电阻数值、直轴不饱和同步电抗、交轴不饱和同步电抗、直轴暂态电抗、交轴暂态电抗、直轴暂态开路时间常数、交轴暂态开路时间常数、直轴次暂态电抗、交轴次暂态电抗、直轴次暂态时间常数、交轴次暂态时间常数等属性进行赋值,机电暂态数据不具有的数据,则使用默认的电磁暂态数据库中的数据进行赋值。
步骤4-2-1、基于节点类,派生得到负荷类,并增加BasekV、S N等属性,以及一个负荷类型选择方法ChooseLoad( )。
步骤4-2-2、基于负荷类,派生得到固定负荷,增加设置功率的实例方法SetLoad(double P 0、double Q 0、double V N、double NP、double NQ、double k PF 、double k QF ),该实例方法用于设置固定负荷消耗的有功功率和无功功率,拥有额定有功P 0、额定无功Q 0、额定电压V N、有功电压变化系数NP、无功电压变化系数NQ、有功频率变化系数k PF 和无功频率变化系数k QF等七个形参参数,返回值则为固定负荷消耗的实际有功功率和无功功率,并满足:
Figure 365351DEST_PATH_IMAGE001
式中,PQ分别为固定负荷吸收的有功功率和无功功率,P 0Q 0分别为额定有功功率和额定无功功率,VV N分别为相电压有效值的实际测量值和额定值,NPNQ分别为有功的电压系数和无功的电压变化系数,Δf为系统的频率变化量,k PFk QF分别为有功和无功的频率变化系数。
通过上述实例方法Setload(double P 0、double Q 0、double V N、double NP、doubleNQ、double k PF 、double k QF )修改P 0Q 0V NNPNQk PFk QF等系数来模拟不同类型的固定负荷。
实例化后,利用机电暂态程序的数据对各个负荷进行赋值,对于机电暂态数据不具有的数据,则使用默认的电磁暂态数据库中的数据进行赋值。需要特别注意的是,若不考虑负荷的频率特性,则令k PF = 0和k QF = 0;若负荷为无功补偿装置,则令NP = 0,P 0 = 0,NQ= 2。
步骤4-2-3、基于负荷类,派生得到异步电动机负荷,增加GeneData、Opt、SandR_R、SandR_L、MutualSatu、LeakSatu、Moni等结构体来分类存储异步电动机的属性。实例化后,利用机电暂态程序的数据对各台异步电动机的额定转速、机械阻尼系数、定子电阻、定子电抗、转子电阻、转子电抗、激磁电抗等属性进行赋值。对于机电暂态数据不具有的数据,则使用默认的电磁暂态数据库中的数据进行赋值。
步骤5、基于支路类,派生得到线路类、变压器类、换流器类和控制器类。
步骤5-1、基于支路类,派生得到交直流输电线路,增加Config、CopiedData、RXL等结构体来分类存储交直流输电的属性。实例化后,利用机电暂态程序的数据对各条交直流输电线路的Length、BasekV、S NRX LX C等属性进行赋值。对于机电暂态数据不具有的数据,则使用默认的电磁暂态数据库中的数据进行赋值。
步骤5-2、基于支路类,派生得到变压器类,增加Config、WindkV、Satu、Magn等属性。实例化后,利用机电暂态程序的数据对变压器的S N、BaseFreq、X L、CopperLoss、IronLoss等属性进行赋值。对于机电暂态数据不具有的数据,则使用默认的电磁暂态数据库中的数据进行赋值。
步骤5-3、基于支路类,派生得到换流器类(本实施例中,换流器类具体指晶闸管型换流器),增加Config、PLO、ValveData等属性。实例化后,利用机电暂态程序的数据对各个换流器的Type、VoltDrop等属性进行设置。对于机电暂态数据不具有的数据,则使用默认的电磁暂态数据库中的数据进行赋值。
步骤5-4、基于支路类,派生得到控制器类,增加SignalInput、SignalOutput等属性,以及一个控制参数嵌套类CtrlSet。实例化后,利用机电暂态程序的数据对各个控制器的SignalInput、SignalOutput和CtrlPara进行初始化设置。对于机电暂态数据不具有的数据,则使用默认的电磁暂态数据库中的数据进行初始化。
还包括:基于控制器类,派生得到同步发电机的控制器(具体包括励磁机、PSS、原动机和调速器等),并设置输入信号、输出信号,以及控制参数。
基于控制器类,派生得到异步电动机的控制器,并设置它的输入信号、输出信号,以及控制参数。
基于控制器类,派生得到换流器的控制装置,并设置它的输入信号、输出信号,以及控制参数。
步骤6、潮流校核和误差计算。对比机电暂态仿真程序和步骤1-5所得电磁暂态程序分别计算得到的有功功率、无功功率和母线电压有效值,并分别计算有功功率、无功功率和母线电压有效值的误差系数:
Figure 326353DEST_PATH_IMAGE002
其中,
Figure 974504DEST_PATH_IMAGE003
Figure 780786DEST_PATH_IMAGE004
Figure 80049DEST_PATH_IMAGE005
分别为有功误差、无功误差和母线电压有效值误差;P 1i Q 1i V 1i 分别为由机电暂态程序计算得到的有功功率、无功功率和母线电压有效值;P 2i Q 2i V 2i 分别为由所述电磁暂态模型计算得到的有功功率、无功功率和母线电压有效值;S N 为交直 流输电系统的基准容量,V Ni 为第i个测量位置的额定电压。
若有功功率、无功功率和母线电压有效值的误差系数:
Figure 946374DEST_PATH_IMAGE006
则步骤1-5所得电磁暂态程序正确;反之,则需排查程序1-5,直至满足误差系数要求。
具体排查内容:
1)核对自动建模程序有无正确读取全系统电力设备对象的机电暂态数据信息;
2)检查全系统电力设备对象的拓扑连接关系是否与机电暂态数据提供的信息一致。
步骤7、基于补充的换流器详细电磁暂态参数和交直流线路分布参数,对自动生成的交直流输电系统的设备类型和参数进行校正,得到最终适用于交直流输电系统的线路过电压和宽频振荡问题研究的电磁暂态模型,如图2~图5所示。
图2为118节点交直流输电系统电磁暂态模型的整体结构,图中实线为交直流输电线路,中间为实心圆的双圈结构表示发电厂、变电站或者换流站,用于存放除交直流输电线路以外的电力设备实例对象。
图3为交直流输电系统中的一组换流器实例对象,其中左边为整流设备,右边为逆变设备,AM为触发角测量信号,GM关断角测量信号、AO为触发角控制信号、KB为整流设备和逆变设备的开关控制信号。
图4为交直流输电系统的一个名为LCCTRL17_18的控制器实例对象,图中SignalInput为控制器的输入信号端子,SignalOutput为控制器的输出信号端子。
图5为采用分布参数的某一交直流输电线路实例对象,需要指出的是,实例对象采用的是频率相关(相域)模型,其参数设置界面由频率相关(相域)模型设置手册、杆塔参数设置和地平面参数设置界面等部分组成。其中,G1、G2为避雷线,C1、C2、C3为输电线路导体。
图5中,频率相关(相域)模型设置手册详细内容如下:
插补时间设置:开;
曲线拟合起始频:0.5[Hz];
曲线拟合截止频率:1.0E6[Hz];
曲线拟合步数:100;
输入特征导纳极点数:20;
输入特征导纳最大拟合误差:2[%];
输入传输函数每延迟组最大极点数:20;
输入传输函数的最大拟合误差:2[%];
直流校正设置:关闭;
无源性校核:关闭。
图5中,地平面参数设置界面具体内容如下:
电地阻率:100 [ohm*m];
空中互感量积分方法的设置:分析近似(Deri-Semlyen);
地下互感量积分方法的设置:直接数值积分;
空中与地下互感量积分方法的设置:解析近似(Lucca)。
本发明调用现有电磁暂态仿真程序EMTDC提供的分布参数线路模型,参数设置界面如图5所示,且相关参数均使用默认值即可。
综上,本发明将包含换流器在内的多种电力设备及其拓扑关系抽象为节点和支路等父类,并派生得到电源、负荷、线路、变压器、换流器、控制器等六个子类。结合机电暂态程序和默认的电磁暂态数据,得到初始的电磁暂态模型。采用交直流输电线路的分布参数、变压器的磁路参数和换流器详细参数等电磁暂态参数对初始的电磁暂态模型进行修正,得到宽频段的电磁暂态模型。

Claims (5)

1.一种大规模交直流输电系统电磁暂态模型自动生成方法,其特征在于,包括以下步骤:
步骤1:分析机电暂态程序和默认的电磁暂态数据,抽象得到关注区域交直流输电系统的UML类图;并定义顶级父类及其子类,令其具有读数据、写数据、设置坐标的静态方法;
步骤2:依据电力设备节点的数目,抽象得到包括节点类和支路类的次顶级父类;
步骤3:基于节点类,派生得到电源类和负荷类;
步骤4:基于支路类,派生得到交直流输电线路类、变压器类、换流器类和控制器类;至此得到的初始的电磁暂态模型;
步骤5:校核潮流并计算潮流误差,对比机电暂态程序和所述电磁暂态模型的程序,分别计算由机电暂态程序计算得到的有功功率、无功功率和母线电压有效值,并计算N个测量位置的潮流误差;若潮流误差的无穷范数小于误差上限值,则所述电磁暂态模型的程序正确;反之,则排查所述电磁暂态模型的程序,直至满足条件;
步骤6:基于补充的电磁暂态数据,对交直流输电系统的元件类型和参数进行补充和修正;基于包括交直流输电线路的分布参数、换流器的详细控制参数,变压器的磁路参数的详细电磁暂态数据,根据研究对象的不同,选择一项或多项数据来修正所述初始的电磁暂态模型,其余参数仍采用默认的电磁暂态数据;得到最终适用于交直流输电系统的电磁暂态模型。
2.根据权利要求1所述的一种大规模交直流输电系统电磁暂态模型自动生成方法,其特征在于,所述步骤2具体为:将节点数目为1的电力设备抽象为节点类,并增加节点名称和节点位置属性;将节点数目大于1的电力设备抽象为支路类,并增加支路名称和支路位置属性。
3.根据权利要求1所述的一种大规模交直流输电系统电磁暂态模型自动生成方法,其特征在于,所述步骤3具体为:
步骤3.1.1:基于节点类,派生得到电源类,并增加包括额定电压、额定频率、额定容量的属性;
步骤3.1.2:基于电源类,派生得到三相电压源,并增加包括阻抗类型、阻抗数值、中心点接地选择、相位大小的电压源属性;
步骤3.1.3:基于电源类,派生得到同步发电机,并增加包括定子电阻数值、直轴不饱和同步电抗、交轴不饱和同步电抗、直轴暂态电抗、交轴暂态电抗、直轴暂态开路时间常数、交轴暂态开路时间常数、直轴次暂态电抗、交轴次暂态电抗、直轴次暂态时间常数、交轴次暂态时间常数的属性;
步骤3.2.1:基于节点类,派生得到负荷类,并增加包括额定容量和额定电压的属性;
步骤3.2.2:基于负荷类,派生得到固定负荷;固定负荷拥有一个设置功率的实例方法SetLoad(double P 0、double Q 0、double V N、double NP、double NQ、double k PF 、double k QF ),满足:
Figure 718167DEST_PATH_IMAGE001
式中,PQ分别为固定负荷吸收的有功功率和无功功率,P 0Q 0分别为额定有功功率和额定无功功率,VV N分别为相电压有效值的实际测量值和额定值,NPNQ分别为有功的电压系数和无功的电压变化系数,Δf为系统的频率变化量,k PFk QF分别为有功和无功的频率变化系数;
通过所述实例方法修改P 0Q 0V NNPNQk PFk QF系数来模拟不同类型的固定负荷;
步骤3.2.3:基于负荷类,派生得到异步电动机负荷,增加包括额定转速、机械阻尼系数、定子电阻、定子电抗、转子电阻、转子电抗、激磁电抗的属性。
4.根据权利要求1所述的一种大规模交直流输电系统电磁暂态模型自动生成方法,其特征在于,所述步骤4具体为:
步骤4.1:基于支路类,派生得到交直流输电线路,增加包括额定频率、正序阻抗、零序阻抗、线路长度的属性;
步骤4.2:基于支路类,派生得到变压器类,增加包括额定频率、额定容量、额定绕组电压、铜耗、铁耗、激磁损耗的属性;
步骤4.3:基于支路类,派生得到换流器类,增加包括换流器个数、每个换流器压降数、桥电流额定值发属性;
步骤4.4:基于支路类,派生得到控制器类,增加包括输入信号、输出信号的属性,以及一个嵌套类ControlSet;
步骤4.5:基于控制器类,派生得到同步发电机的控制器,并设置输入信号、输出信号,以及控制参数;
步骤4.6:基于控制器类,派生得到异步电动机的控制器,并设置它的输入信号、输出信号,以及控制参数;
步骤4.7:基于控制器类,派生得到换流器的控制装置,并设置它的输入信号、输出信号,以及控制参数。
5.根据权利要求1所述的一种大规模交直流输电系统电磁暂态模型自动生成方法,其特征在于,所述步骤5中计算N个测量位置的潮流误差,具体为计算有功功率、无功功率和母线电压有效值的误差系数:
Figure 41832DEST_PATH_IMAGE002
式中,
Figure 352727DEST_PATH_IMAGE003
Figure 934887DEST_PATH_IMAGE004
Figure 155784DEST_PATH_IMAGE005
分别为有功误差、无功误差和母线电压有效值误差;P 1i Q 1i V 1i 分 别为由机电暂态程序计算得到的有功功率、无功功率和母线电压有效值;P 2i Q 2i V 2i 分别 为由所述电磁暂态模型计算得到的有功功率、无功功率和母线电压有效值;S N 为交直流输 电系统的基准容量,V Ni 为第i个测量位置的额定电压;
若有功功率、无功功率和母线电压有效值的误差系数满足:
Figure 243826DEST_PATH_IMAGE006
则所述电磁暂态模型的程序正确。
CN202210942602.3A 2022-08-08 2022-08-08 一种大规模交直流输电系统电磁暂态模型自动生成方法 Active CN115021256B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210942602.3A CN115021256B (zh) 2022-08-08 2022-08-08 一种大规模交直流输电系统电磁暂态模型自动生成方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210942602.3A CN115021256B (zh) 2022-08-08 2022-08-08 一种大规模交直流输电系统电磁暂态模型自动生成方法

Publications (2)

Publication Number Publication Date
CN115021256A true CN115021256A (zh) 2022-09-06
CN115021256B CN115021256B (zh) 2022-10-21

Family

ID=83065865

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210942602.3A Active CN115021256B (zh) 2022-08-08 2022-08-08 一种大规模交直流输电系统电磁暂态模型自动生成方法

Country Status (1)

Country Link
CN (1) CN115021256B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116822436A (zh) * 2023-06-30 2023-09-29 四川大学 一种直流输电送端交流系统的振荡风险灵敏度分析方法

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102999668A (zh) * 2012-11-27 2013-03-27 中国电力科学研究院 一种机电暂态数据至电磁暂态数据的自动转换方法
CN103559347A (zh) * 2013-10-31 2014-02-05 浙江大学 一种大规模交直流电力系统电磁暂态仿真模型的构建方法
CN106886616A (zh) * 2015-12-15 2017-06-23 中国电力科学研究院 一种大规模电磁暂态电网仿真的自动分网方法
CN107664954A (zh) * 2017-10-16 2018-02-06 国网四川省电力公司电力科学研究院 一种电力系统电磁暂态可视化仿真模型自动生成方法
CN109638871A (zh) * 2019-01-28 2019-04-16 东北电力大学 考虑风电接入的大规模交直流混联系统主网划分方法
CN109766586A (zh) * 2018-12-14 2019-05-17 中国电力科学研究院有限公司 一种自动生成大规模电网电磁暂态仿真模型的方法及系统
US20210049313A1 (en) * 2019-08-16 2021-02-18 Shanghai Jiao Tong University Electromagnetic transient simulation method for field programmable logic array
CN112434409A (zh) * 2020-11-13 2021-03-02 国家电网有限公司 基于恒定阻抗接口的同步发电机全电磁暂态仿真方法
CN112736896A (zh) * 2020-12-11 2021-04-30 国家电网有限公司 一种适用于全电磁暂态仿真的三阶段初始化方法及系统
CN112883592A (zh) * 2021-03-26 2021-06-01 国网黑龙江省电力有限公司电力科学研究院 一种电网全电磁暂态建模方法

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102999668A (zh) * 2012-11-27 2013-03-27 中国电力科学研究院 一种机电暂态数据至电磁暂态数据的自动转换方法
CN103559347A (zh) * 2013-10-31 2014-02-05 浙江大学 一种大规模交直流电力系统电磁暂态仿真模型的构建方法
CN106886616A (zh) * 2015-12-15 2017-06-23 中国电力科学研究院 一种大规模电磁暂态电网仿真的自动分网方法
CN107664954A (zh) * 2017-10-16 2018-02-06 国网四川省电力公司电力科学研究院 一种电力系统电磁暂态可视化仿真模型自动生成方法
CN109766586A (zh) * 2018-12-14 2019-05-17 中国电力科学研究院有限公司 一种自动生成大规模电网电磁暂态仿真模型的方法及系统
CN109638871A (zh) * 2019-01-28 2019-04-16 东北电力大学 考虑风电接入的大规模交直流混联系统主网划分方法
US20210049313A1 (en) * 2019-08-16 2021-02-18 Shanghai Jiao Tong University Electromagnetic transient simulation method for field programmable logic array
CN112434409A (zh) * 2020-11-13 2021-03-02 国家电网有限公司 基于恒定阻抗接口的同步发电机全电磁暂态仿真方法
CN112736896A (zh) * 2020-12-11 2021-04-30 国家电网有限公司 一种适用于全电磁暂态仿真的三阶段初始化方法及系统
CN112883592A (zh) * 2021-03-26 2021-06-01 国网黑龙江省电力有限公司电力科学研究院 一种电网全电磁暂态建模方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
ZHENDONG TAN 等: ""Automatic Generation and Parameter Verification of Large-scale EMT Simulation Models Based on TSP Projects"", 《2021 POWER SYSTEM AND GREEN ENERGY CONFERENCE (PSGEC)》 *
ZHIHONG LIU 等: ""Automatic Generation and Initialization of EMT Simulation Models for Large-Scale AC-DC Hybrid Power System"", 《 2021 11TH INTERNATIONAL CONFERENCE ON POWER AND ENERGY SYSTEMS (ICPES)》 *
何海林 等: ""用于大规模机电-电磁暂态仿真模型自动转化的分层布局方法"", 《中国电力》 *
程改红 等: "" 大规模交直流电力系统电磁暂态仿真高效建模方法"", 《电网技术》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116822436A (zh) * 2023-06-30 2023-09-29 四川大学 一种直流输电送端交流系统的振荡风险灵敏度分析方法
CN116822436B (zh) * 2023-06-30 2024-02-27 四川大学 一种直流输电送端交流系统的振荡风险灵敏度分析方法

Also Published As

Publication number Publication date
CN115021256B (zh) 2022-10-21

Similar Documents

Publication Publication Date Title
CN108280271B (zh) 基于开关周期平均原理的统一潮流控制器等效建模方法
CN112564175B (zh) 新能源多场站短路比计算方法、装置及存储介质
Jia et al. Fault characteristics and protection adaptability analysis in VSC‐HVDC‐connected offshore wind farm integration system
Liu et al. Impedance modeling of DFIG wind farms with various rotor speeds and frequency coupling
CN115021256B (zh) 一种大规模交直流输电系统电磁暂态模型自动生成方法
Jiang et al. Simplified solid state transformer modeling for real time digital simulator (rtds)
Torquato et al. Practical chart for harmonic resonance assessment of DFIG-Based wind parks
CN113852120A (zh) 一种新能源发电最大并网容量的确定方法及系统
CN109494722B (zh) 电网侧等效阻抗建模方法及系统
CN113224788B (zh) 含有风电接入的交直流混联电网静态电压稳定评估方法
Xu et al. Impedance-based analysis of potential stability risk between grid-forming and grid-following wind turbine systems
Lin et al. Vector-variable based modeling method for large-scale wind farm considering collector system dynamics
Suvorov et al. A novel approach to validation of a positive-sequence modeling of a converter-interfaced generation with hydrogen energy storage in practical power systems
Tang et al. DC-side admittance model of converters with constant power control considering DC-side voltage fluctuation
Im et al. Reactive power control strategy for inverter-based distributed generation system with a programmable limit of the voltage variation at PCC
CN116581757B (zh) 考虑高比例电力电子化设备的负荷模型建模方法及系统
Xiang et al. Equivalent modeling of wind farms considering wake effects and analysis of probability of sub-synchronous oscillations
Zhang et al. Electromechanical Transient Simulation of AC-DC Hybrid Power System Considering MMC Model
Li et al. Parallel impedance-reshaping control strategy for suppressing wide range high-frequency resonance of shunt compensated DFIG system
Zhou et al. Frequency-Coupled Impedance Modeling and Stability Analysis of Multiple Grid-Connected Voltage-Source Converters
Yuan et al. Stability analysis of grid-connected VSC based on impedance modelling
Wang et al. Output Impedance Modeling and Parameter Disturbance Stability Analysis of Photovoltaic Virtual Synchronous Generator
Han et al. A DISCRETE-TIME STATE SPACE MODELLING METHOD APPLYING TO SSO OF RENEWABLE POWER GENERATION
Mo et al. Research on Fault Characteristics of DFIG Wind Farm with Flexible DC Integration
Wang et al. Small-signal dynamics of MMC based DC-grid system

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