CN112736896B - 一种适用于全电磁暂态仿真的三阶段初始化方法及系统 - Google Patents
一种适用于全电磁暂态仿真的三阶段初始化方法及系统 Download PDFInfo
- Publication number
- CN112736896B CN112736896B CN202011459124.8A CN202011459124A CN112736896B CN 112736896 B CN112736896 B CN 112736896B CN 202011459124 A CN202011459124 A CN 202011459124A CN 112736896 B CN112736896 B CN 112736896B
- Authority
- CN
- China
- Prior art keywords
- model
- simulation
- direct current
- synchronous generator
- transient simulation
- 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.)
- Active
Links
- 238000004088 simulation Methods 0.000 title claims abstract description 154
- 230000001052 transient effect Effects 0.000 title claims abstract description 109
- 238000011423 initialization method Methods 0.000 title claims abstract description 11
- 238000004364 calculation method Methods 0.000 claims abstract description 55
- 230000005540 biological transmission Effects 0.000 claims abstract description 47
- 238000000034 method Methods 0.000 claims abstract description 34
- 230000001360 synchronised effect Effects 0.000 claims description 62
- 230000001105 regulatory effect Effects 0.000 claims description 28
- 230000009466 transformation Effects 0.000 claims description 20
- 238000005259 measurement Methods 0.000 claims description 14
- 230000008569 process Effects 0.000 claims description 10
- 238000006243 chemical reaction Methods 0.000 claims description 9
- 238000009826 distribution Methods 0.000 claims description 8
- 230000005284 excitation Effects 0.000 claims description 7
- 239000003381 stabilizer Substances 0.000 claims description 7
- 230000006855 networking Effects 0.000 claims description 6
- 238000010586 diagram Methods 0.000 description 10
- 238000004590 computer program Methods 0.000 description 7
- 230000008033 biological extinction Effects 0.000 description 5
- 230000006870 function Effects 0.000 description 4
- 238000012986 modification Methods 0.000 description 4
- 230000004048 modification Effects 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 238000003860 storage Methods 0.000 description 3
- 230000004075 alteration Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 239000003990 capacitor Substances 0.000 description 1
- 230000009194 climbing Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/04—Circuit arrangements for ac mains or ac distribution networks for connecting networks of the same frequency but supplied from different sources
- H02J3/06—Controlling transfer of power between connected networks; Controlling sharing of load between connected networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/04—Power grid distribution networks
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/10—Power 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
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
- Y02E60/60—Arrangements for transfer of electric power between AC networks or generators via a high voltage DC link [HVCD]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Power Engineering (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
技术领域
本发明涉及电力系统建模与仿真技术领域，并且更具体地，涉及一种适用于全电磁暂态仿真的三阶段初始化方法及系统。
背景技术
目前关于电磁仿真及其应用的研究很多，但大多数集中在如何建立电力系统元件的模型及如何提高仿真精度和计算效率方面。全电磁暂态仿真技术应用于大电网稳定性分析主要是为了解决直流输电系统接入后仿真不准确的问题。在交直流混联系统中的全电磁仿真中，模型及系统的初始化在电力系统稳定性分析中是一个必选项，无论是对故障的复现还是电网方式的安排，都需要根据电网不同的运行状态分析电网的稳定性。
在机电暂态仿真中，由于潮流计算与机电暂态的电路在机电暂态程序中，由于潮流计算与机电暂态的电路部分都采用相同的基波向量模型，初始化并不是难题。机电暂态程序采用潮流结果进行初始化，在仿真开始时就会对潮流结果进行校验，存在电压偏差时会报错并退出计算。机电暂态中主要对控制系统的动态进行详细建模，潮流中将控制特性改为了PQ、PV、Vθ节点，控制器的限幅则变为了Qmax、Qmin等限制量，会产生限幅的冲突。除此之外，没有其他初始化问题。
电磁暂态仿真的初始化方法主要分为三种方式：相量解方法、电磁暂态法、潮流解方法。全电磁暂态仿真由于仿真规模较大，向量解方法和电磁暂态法的适用性并不强，一般采用基于潮流解的仿真初始化方法进行求解。但是，直流输电系统由于潮流模型与电磁模型的不一致，其初始化面临着较大的困难，初始化时间更长。
发明内容
针对上述问题本发明提出了一种适用于全电磁暂态仿真的三阶段初始化方法，包括：
针对直流输电系统模型，读入模型参数，并进行潮流计算；
通过潮流计算确定机电暂态仿真模型数据，并将机电暂态仿真模型数据转换为全电暂态仿真模型数据；
针对全电暂态仿真模型数据，进入三阶段仿真，包括：
第一阶段，闭锁直流触发系统，确定交流侧两端电压源的幅值和角度为潮流计算结果，根据潮流计算结果对交流侧同步发电机进行充电和初始化；
第二阶段，利用潮流计算结果中交流侧两端电压源的角度锁定直流控制系统，启动直流电路以及触发系统，对直流电路的预定元件进行充电，并调整交直流接口处的有功功率和无功功率误差；
第三阶段，将直流输电系统的直流控制系统进行反向初始化，并将控制系统放开，使直流输电系统进入全电磁暂态的常规仿真。
可选的，将机电暂态仿真模型数据转换为全电暂态仿真模型数据，准换后全电暂态仿真模型数据中发电机模型采用
采用dq变换模型；全电暂态仿真模型数据中线路模型，包括：网络内部采用PI模型，分网联络线采用分布参数模型；全电暂态仿真模型数据中负荷模型采用恒阻抗模型，三相等效模型使用对称的星型接线方式。
可选的，第一阶段的同步发电机的充电过程，包括：
将同步发电机转换为电压源模型，对电压源模型进行预设时间的仿真；
将电压源模型转换为同步发电机模型，电磁方程部分放开，转动方程部分锁定，采用恒转速模式进行计算；
将同步发电机模型的转动方程部分放开，同步发电机模型部分转化为正常模式仿真，获取测量值；
利用同步发电机模型的测量值，对除同步发电机模型外的模型进行初始化，并将控制系统接入，控制交流系统进行仿真；
除同步发电机模型外的模型，包括：原动机模型、调速控制系统模型、励磁系统模型和电力系统稳定器PSS模型。
可选的，调整交直流接口处的有功功率和无功功率误差，包括：
根据电磁暂态中元件的参数，调整全电暂态仿真模型数据中的参数；
锁定控制系统中的定电压控制环节；
调节逆变侧变压器的变比，仿真计算直流输电系统初始化完成后的直流电压，并调整直流电压与潮流计算结果一致；
调节整流侧变压器的变比，仿真计算直流输电系统初始化完成后的整流侧触发角，并调整整流侧触发角与潮流结果一致；
放开定直流电压控制环节，进行仿真，检查两侧交流有功功率和无功功率与潮流结果的偏差，若偏差符合预设阈值，则完成初始化调整。
可选的，第三阶段，根据直流输电系统控制方式进行不同的初始化计算，逆变侧的熄弧角控制方式包括实测型和预测型；
所述实测型控制利用PI进行闭环控制；
所述预测型通过熄弧角对逆变侧进行触发控制。
本发明还提出了一种适用于全电磁暂态仿真的三阶段初始化系统，包括：
潮流计算单元，针对直流输电系统模型，读入模型参数，并进行潮流计算；
模型转换单元，通过潮流计算确定机电暂态仿真模型数据，并将机电暂态仿真模型数据转换为全电暂态仿真模型数据；
仿真单元，针对全电暂态仿真模型数据，进入三阶段仿真，包括：
第一阶段，闭锁直流触发系统，确定交流侧两端电压源的幅值和角度为潮流计算结果，根据潮流计算结果对交流侧同步发电机进行充电和初始化；
第二阶段，利用潮流计算结果中交流侧两端电压源的角度锁定直流控制系统，启动直流电路以及触发系统，对直流电路的预定元件进行充电，并调整交直流接口处的有功功率和无功功率误差；
第三阶段，将直流输电系统的直流控制系统进行反向初始化，并将控制系统放开，使直流输电系统进入全电磁暂态的常规仿真。
可选的，将机电暂态仿真模型数据转换为全电暂态仿真模型数据，准换后全电暂态仿真模型数据中发电机模型采用
采用dq变换模型；全电暂态仿真模型数据中线路模型，包括：网络内部采用PI模型，分网联络线采用分布参数模型；全电暂态仿真模型数据中负荷模型采用恒阻抗模型，三相等效模型使用对称的星型接线方式。
可选的，第一阶段的同步发电机的充电过程，包括：
将同步发电机转换为电压源模型，对电压源模型进行预设时间的仿真；
将电压源模型转换为同步发电机模型，电磁方程部分放开，转动方程部分锁定，采用恒转速模式进行计算；
将同步发电机模型的转动方程部分放开，同步发电机模型部分转化为正常模式仿真，获取测量值；
利用同步发电机模型的测量值，对除同步发电机模型外的模型进行初始化，并将控制系统接入，控制交流系统进行仿真；
除同步发电机模型外的模型，包括：原动机模型、调速控制系统模型、励磁系统模型和电力系统稳定器PSS模型。
可选的，调整交直流接口处的有功功率和无功功率误差，包括：
根据电磁暂态中元件的参数，调整全电暂态仿真模型数据中的参数；
锁定控制系统中的定电压控制环节；
调节逆变侧变压器的变比，仿真计算直流输电系统初始化完成后的直流电压，并调整直流电压与潮流计算结果一致；
调节整流侧变压器的变比，仿真计算直流输电系统初始化完成后的整流侧触发角，并调整整流侧触发角与潮流结果一致；
放开定直流电压控制环节，进行仿真，检查两侧交流有功功率和无功功率与潮流结果的偏差，若偏差符合预设阈值，则完成初始化调整。
可选的，第三阶段，根据直流输电系统控制方式进行不同的初始化计算，逆变侧的熄弧角控制方式包括实测型和预测型；
所述实测型控制利用PI进行闭环控制；
所述预测型通过熄弧角对逆变侧进行触发控制。
本发明的具有较好的可操作性，适用于大多数直流输电系统的初始化，也可以适应大多数的电磁暂态仿真软件，具有更好的计算适应性。
附图说明
图1为本发明一种适用于全电磁暂态仿真的三阶段初始化方法的流程图；
图2为本发明一种适用于全电磁暂态仿真的三阶段初始化方法的三阶段初始化示意图；
图3为本发明一种适用于全电磁暂态仿真的三阶段初始化方法常规直流潮流模型图；
图4为本发明一种适用于全电磁暂态仿真的三阶段初始化系统的结构图。
具体实施方式
现在参考附图介绍本发明的示例性实施方式，然而，本发明可以用许多不同的形式来实施，并且不局限于此处描述的实施例，提供这些实施例是为了详尽地且完全地公开本发明，并且向所属技术领域的技术人员充分传达本发明的范围。对于表示在附图中的示例性实施方式中的术语并不是对本发明的限定。在附图中，相同的单元/元件使用相同的附图标记。
除非另有说明，此处使用的术语(包括科技术语)对所属技术领域的技术人员具有通常的理解含义。另外，可以理解的是，以通常使用的词典限定的术语，应当被理解为与其相关领域的语境具有一致的含义，而不应该被理解为理想化的或过于正式的意义。
本发明提出了一种适用于全电磁暂态仿真的三阶段初始化方法，如图1所示，包括：
针对直流输电系统模型，读入模型参数，并进行潮流计算；
通过潮流计算确定机电暂态仿真模型数据，并将机电暂态仿真模型数据转换为全电暂态仿真模型数据；
针对全电暂态仿真模型数据，进入三阶段仿真，包括：
第一阶段，闭锁直流触发系统，确定交流侧两端电压源的幅值和角度为潮流计算结果，根据潮流计算结果对交流侧同步发电机进行充电和初始化；
第二阶段，利用潮流计算结果中交流侧两端电压源的角度锁定直流控制系统，启动直流电路以及触发系统，对直流电路的预定元件进行充电，并调整交直流接口处的有功功率和无功功率误差；
第三阶段，将直流输电系统的直流控制系统进行反向初始化，并将控制系统放开，使直流输电系统进入全电磁暂态的常规仿真。
将机电暂态仿真模型数据转换为全电暂态仿真模型数据，准换后全电暂态仿真模型数据中发电机模型采用
采用dq变换模型；全电暂态仿真模型数据中线路模型，包括：网络内部采用PI模型，分网联络线采用分布参数模型；全电暂态仿真模型数据中负荷模型采用恒阻抗模型，三相等效模型使用对称的星型接线方式。
第一阶段的同步发电机的充电过程，包括：
将同步发电机转换为电压源模型，对电压源模型进行预设时间的仿真；
将电压源模型转换为同步发电机模型，电磁方程部分放开，转动方程部分锁定，采用恒转速模式进行计算；
将同步发电机模型的转动方程部分放开，同步发电机模型部分转化为正常模式仿真，获取测量值；
利用同步发电机模型的测量值，对除同步发电机模型外的模型进行初始化，并将控制系统接入，控制交流系统进行仿真；
除同步发电机模型外的模型，包括：原动机模型、调速控制系统模型、励磁系统模型和电力系统稳定器PSS模型。
调整交直流接口处的有功功率和无功功率误差，包括：
根据电磁暂态中元件的参数，调整全电暂态仿真模型数据中的参数；
锁定控制系统中的定电压控制环节；
调节逆变侧变压器的变比，仿真计算直流输电系统初始化完成后的直流电压，并调整直流电压与潮流计算结果一致；
调节整流侧变压器的变比，仿真计算直流输电系统初始化完成后的整流侧触发角，并调整整流侧触发角与潮流结果一致；
放开定直流电压控制环节，进行仿真，检查两侧交流有功功率和无功功率与潮流结果的偏差，若偏差符合预设阈值，则完成初始化调整。
第三阶段，根据直流输电系统控制方式进行不同的初始化计算，逆变侧的熄弧角控制方式包括实测型和预测型；
实测型控制利用PI进行闭环控制；
预测型通过熄弧角对逆变侧进行触发控制。
下面结合实施例对本发明进行进一步的说明：
本发明中仿真示意图如图2所示，包括：
步骤1：读入网络中的母线、线路、变压器、发电机、负荷等模型参数，直流输电系统的模型如图3所示，具体的潮流计算公式为(1)和(2)，利用牛顿拉夫逊方法进行大电网的潮流计算。
步骤2：将机电暂态仿真模型数据转换为全电磁暂态仿真数据，其中机电暂态模型简化处理为对称的三相模型。
步骤3：第一阶段仿真。闭锁直流触发系统，交流侧两端电压源的幅值和角度为潮流计算结果，对交流侧元件进行充电和初始化计算。
步骤4：第二阶段仿真。利用潮流计算结果中两侧的触发角锁定直流控制系统，启动直流电路以及触发系统，对直流线路部分的元件进行充电，并进行调整使交直流接口处的有功功率和无功功率误差尽量减小。
步骤5：第三阶段仿真。将直流输电系统的控制系统进行反向初始化，并将直流系统的控制系统放开，进入全电磁暂态的正常仿真。
步骤2中，转换后发电机采用dq变换模型，线路分为两种情况，网络内部采用PI模型，分网联络线采用分布参数模型，负荷采用简化的恒阻抗模型，三相等效采用对称的星型接线方式；
步骤3中的同步发电机充电过程包含以下步骤：
步骤301：同步发电机转换为电压源模型，仿真时间为0.2s，其中电压爬坡时间为0.1s，即前0.1s电压源的电压从0线性增加到潮流计算解。
步骤302：将电压源模型改为同步发电机模型，电磁方程部分放开，转动方程部分锁定，采用恒转速模式计算0.15s，在此过程中同步发电机的转速偏差保持为0，即无偏差。
步骤303：将同步发电机的转动方程部分放开，同步发电机模型部分转化为正常模式，仿真0.15s。
步骤304：利用同步发电机的测量值对同步发电机的其他模型进行初始化，包括原动机模型、调速控制系统、励磁系统、电力系统稳定器(PSS)模型，将控制系统接入，整个交流系统进行正常仿真。
步骤4的调整包含以下步骤：
步骤401：按照电磁暂态中元件的参数，调整潮流模型中的参数，特别时变压器的参数、线路的参数、滤波器电容等的参数。
步骤402：锁定控制系统中的定电压控制环节，保证在初始化过程中，定电压控制部分不起作用。
步骤403：调节逆变侧变压器的变比，仿真计算并观察直流系统初始化完成后的直流电压，反复调整直到与潮流结果一致。
步骤404：调节整流侧变压器的变比，仿真计算并观察直流系统初始化完成后的整流侧触发角，反复调整直到与潮流结果一致。
步骤405：放开定直流电压控制环节，进行仿真，检查两侧交流有功功率和无功功率与潮流结果是否一致，若偏差较大，则是潮流模型中的参数存在错误，重新执行步骤403，检查模型参数中的错误，若偏差较小，则完成初始化调整。
步骤5需要根据直流输电系统控制方式进行不同的初始化计算方式。逆变侧的熄弧角控制方式分为实测型和预测型两种：实测型控制利用PI进行闭环控制；预测型控制通过计算得到的熄弧角对逆变侧进行触发控制，其实质是一开环控制，熄弧角的计算公式(3)。为了得到与潮流结果匹配的结果，对于预测型熄弧角控制系统，需要设置合适的dx参数才能获得合适的熄弧角，该参数由公式(3)计算得到。
本发明还提出了一种适用于全电磁暂态仿真的三阶段初始化系统200，如图4所示，包括：
潮流计算单元201，针对直流输电系统模型，读入模型参数，并进行潮流计算；
模型转换单元202，通过潮流计算确定机电暂态仿真模型数据，并将机电暂态仿真模型数据转换为全电暂态仿真模型数据；
仿真单元203，针对全电暂态仿真模型数据，进入三阶段仿真，包括：
第一阶段，闭锁直流触发系统，确定交流侧两端电压源的幅值和角度为潮流计算结果，根据潮流计算结果对交流侧同步发电机进行充电和初始化；
第二阶段，利用潮流计算结果中交流侧两端电压源的角度锁定直流控制系统，启动直流电路以及触发系统，对直流电路的预定元件进行充电，并调整交直流接口处的有功功率和无功功率误差；
第三阶段，将直流输电系统的直流控制系统进行反向初始化，并将控制系统放开，使直流输电系统进入全电磁暂态的常规仿真。
将机电暂态仿真模型数据转换为全电暂态仿真模型数据，准换后全电暂态仿真模型数据中发电机模型采用
采用dq变换模型；全电暂态仿真模型数据中线路模型，包括：网络内部采用PI模型，分网联络线采用分布参数模型；全电暂态仿真模型数据中负荷模型采用恒阻抗模型，三相等效模型使用对称的星型接线方式。
第一阶段的同步发电机的充电过程，包括：
将同步发电机转换为电压源模型，对电压源模型进行预设时间的仿真；
将电压源模型转换为同步发电机模型，电磁方程部分放开，转动方程部分锁定，采用恒转速模式进行计算；
将同步发电机模型的转动方程部分放开，同步发电机模型部分转化为正常模式仿真，获取测量值；
利用同步发电机模型的测量值，对除同步发电机模型外的模型进行初始化，并将控制系统接入，控制交流系统进行仿真；
除同步发电机模型外的模型，包括：原动机模型、调速控制系统模型、励磁系统模型和电力系统稳定器PSS模型。
调整交直流接口处的有功功率和无功功率误差，包括：
根据电磁暂态中元件的参数，调整全电暂态仿真模型数据中的参数；
锁定控制系统中的定电压控制环节；
调节逆变侧变压器的变比，仿真计算直流输电系统初始化完成后的直流电压，并调整直流电压与潮流计算结果一致；
调节整流侧变压器的变比，仿真计算直流输电系统初始化完成后的整流侧触发角，并调整整流侧触发角与潮流结果一致；
放开定直流电压控制环节，进行仿真，检查两侧交流有功功率和无功功率与潮流结果的偏差，若偏差符合预设阈值，则完成初始化调整。
第三阶段，根据直流输电系统控制方式进行不同的初始化计算，逆变侧的熄弧角控制方式包括实测型和预测型；
实测型控制利用PI进行闭环控制；
预测型通过熄弧角对逆变侧进行触发控制。
本发明的具有较好的可操作性，适用于大多数直流输电系统的初始化，也可以适应大多数的电磁暂态仿真软件，具有更好的计算适应性。
本领域内的技术人员应明白，本申请的实施例可提供为方法、系统、或计算机程序产品。因此，本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。本申请实施例中的方案可以采用各种计算机语言实现，例如，面向对象的程序设计语言Java和直译式脚本语言JavaScript等。
本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
尽管已描述了本申请的优选实施例，但本领域内的技术人员一旦得知了基本创造性概念，则可对这些实施例作出另外的变更和修改。所以，所附权利要求意欲解释为包括优选实施例以及落入本申请范围的所有变更和修改。
显然，本领域的技术人员可以对本申请进行各种改动和变型而不脱离本申请的精神和范围。这样，倘若本申请的这些修改和变型属于本申请权利要求及其等同技术的范围之内，则本申请也意图包含这些改动和变型在内。
Claims (10)
1.一种适用于全电磁暂态仿真的三阶段初始化方法，所述方法包括：
针对直流输电系统模型，读入模型参数，并进行潮流计算；
通过潮流计算确定机电暂态仿真模型数据，并将机电暂态仿真模型数据转换为全电磁暂态仿真模型数据；
针对全电磁暂态仿真模型数据，进入三阶段仿真，包括：
第一阶段，闭锁直流触发系统，确定交流侧两端电压源的幅值和角度为潮流计算结果，根据潮流计算结果对交流侧同步发电机进行充电和初始化；
第二阶段，利用潮流计算结果中交流侧两端电压源的角度锁定直流控制系统，启动直流电路以及触发系统，对直流电路的预定元件进行充电，并调整交直流接口处的有功功率误差和无功功率误差；
第三阶段，将直流输电系统的直流控制系统进行反向初始化，并将控制系统放开，使直流输电系统进入全电磁暂态的常规仿真。
2.根据权利要求1所述的方法，所述将机电暂态仿真模型数据转换为全电磁暂态仿真模型数据，转换后全电磁暂态仿真模型数据中发电机模型采用dq变换模型；全电磁暂态仿真模型数据中线路模型，包括：网络内部采用PI模型，分网联络线采用分布参数模型；全电磁暂态仿真模型数据中负荷模型采用恒阻抗模型，三相等效模型使用对称的星型接线方式。
3.根据权利要求1所述的方法，所述第一阶段的同步发电机的充电过程，包括：
将同步发电机转换为电压源模型，对电压源模型进行预设时间的仿真；
将电压源模型转换为同步发电机模型，电磁方程部分放开，转动方程部分锁定，采用恒转速模式进行计算；
将同步发电机模型的转动方程部分放开，同步发电机模型部分转化为正常模式仿真，获取测量值；
利用同步发电机模型的测量值，对除同步发电机模型外的模型进行初始化，并将控制系统接入，控制交流系统进行仿真；
除同步发电机模型外的模型，包括：原动机模型、调速控制系统模型、励磁系统模型和电力系统稳定器PSS模型。
4.根据权利要求1所述的方法，所述调整交直流接口处的有功功率和无功功率误差，包括：
根据电磁暂态中元件的参数，调整全电磁暂态仿真模型数据中的参数；
锁定控制系统中的定电压控制环节；
调节逆变侧变压器的变比，仿真计算直流输电系统初始化完成后的直流电压，并调整直流电压与潮流计算结果一致；
调节整流侧变压器的变比，仿真计算直流输电系统初始化完成后的整流侧触发角，并调整整流侧触发角与潮流结果一致；
放开定直流电压控制环节，进行仿真，检查两侧交流有功功率和无功功率与潮流结果的偏差，若偏差符合预设阈值，则完成初始化调整。
5.根据权利要求1所述的方法，所述第三阶段，根据直流输电系统控制方式进行不同的初始化计算，逆变侧的熄弧角控制方式包括实测型和预测型；
所述实测型控制利用PI进行闭环控制；
所述预测型通过熄弧角对逆变侧进行触发控制。
6.一种适用于全电磁暂态仿真的三阶段初始化系统，所述系统包括：
潮流计算单元，针对直流输电系统模型，读入模型参数，并进行潮流计算；
模型转换单元，通过潮流计算确定机电暂态仿真模型数据，并将机电暂态仿真模型数据转换为全电磁暂态仿真模型数据；
仿真单元，针对全电磁暂态仿真模型数据，进入三阶段仿真，包括：
第一阶段，闭锁直流触发系统，确定交流侧两端电压源的幅值和角度为潮流计算结果，根据潮流计算结果对交流侧同步发电机进行充电和初始化；
第二阶段，利用潮流计算结果中交流侧两端电压源的角度锁定直流控制系统，启动直流电路以及触发系统，对直流电路的预定元件进行充电，并调整交直流接口处的有功功率误差和无功功率误差；
第三阶段，将直流输电系统的直流控制系统进行反向初始化，并将控制系统放开，使直流输电系统进入全电磁暂态的常规仿真。
7.根据权利要求6所述的系统，所述将机电暂态仿真模型数据转换为全电磁暂态仿真模型数据，转换后全电磁暂态仿真模型数据中发电机模型采用dq变换模型；全电磁暂态仿真模型数据中线路模型，包括：网络内部采用PI模型，分网联络线采用分布参数模型；全电磁暂态仿真模型数据中负荷模型采用恒阻抗模型，三相等效模型使用对称的星型接线方式。
8.根据权利要求6所述的系统，所述第一阶段的同步发电机的充电过程，包括：
将同步发电机转换为电压源模型，对电压源模型进行预设时间的仿真；
将电压源模型转换为同步发电机模型，电磁方程部分放开，转动方程部分锁定，采用恒转速模式进行计算；
将同步发电机模型的转动方程部分放开，同步发电机模型部分转化为正常模式仿真，获取测量值；
利用同步发电机模型的测量值，对除同步发电机模型外的模型进行初始化，并将控制系统接入，控制交流系统进行仿真；
除同步发电机模型外的模型，包括：原动机模型、调速控制系统模型、励磁系统模型和电力系统稳定器PSS模型。
9.根据权利要求6所述的系统，所述调整交直流接口处的有功功率和无功功率误差，包括：
根据电磁暂态中元件的参数，调整全电磁暂态仿真模型数据中的参数；
锁定控制系统中的定电压控制环节；
调节逆变侧变压器的变比，仿真计算直流输电系统初始化完成后的直流电压，并调整直流电压与潮流计算结果一致；
调节整流侧变压器的变比，仿真计算直流输电系统初始化完成后的整流侧触发角，并调整整流侧触发角与潮流结果一致；
放开定直流电压控制环节，进行仿真，检查两侧交流有功功率和无功功率与潮流结果的偏差，若偏差符合预设阈值，则完成初始化调整。
10.根据权利要求6所述的系统，所述第三阶段，根据直流输电系统控制方式进行不同的初始化计算，逆变侧的熄弧角控制方式包括实测型和预测型；
所述实测型控制利用PI进行闭环控制；
所述预测型通过熄弧角对逆变侧进行触发控制。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011459124.8A CN112736896B (zh) | 2020-12-11 | 2020-12-11 | 一种适用于全电磁暂态仿真的三阶段初始化方法及系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011459124.8A CN112736896B (zh) | 2020-12-11 | 2020-12-11 | 一种适用于全电磁暂态仿真的三阶段初始化方法及系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112736896A CN112736896A (zh) | 2021-04-30 |
| CN112736896B true CN112736896B (zh) | 2023-11-03 |
Family
ID=75599543
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011459124.8A Active CN112736896B (zh) | 2020-12-11 | 2020-12-11 | 一种适用于全电磁暂态仿真的三阶段初始化方法及系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112736896B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113872229B (zh) * | 2021-09-15 | 2024-02-09 | 国网冀北电力有限公司电力科学研究院 | 一种柔性直流输电操作模拟方法及装置 |
| CN115021256B (zh) * | 2022-08-08 | 2022-10-21 | 四川大学 | 一种大规模交直流输电系统电磁暂态模型自动生成方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105322540A (zh) * | 2015-11-12 | 2016-02-10 | 中国电力科学研究院 | 一种交直流大电网电磁暂态模型的稳态运行方式建立方法 |
| CN109004638A (zh) * | 2018-06-29 | 2018-12-14 | 中国电力科学研究院有限公司 | 一种柔性直流工程机电-电磁混合仿真方法及系统 |
| CN109586275A (zh) * | 2018-10-18 | 2019-04-05 | 清华大学 | 交直流混联系统电磁暂态仿真的快速启动方法与装置 |
| CN109766586A (zh) * | 2018-12-14 | 2019-05-17 | 中国电力科学研究院有限公司 | 一种自动生成大规模电网电磁暂态仿真模型的方法及系统 |
-
2020
- 2020-12-11 CN CN202011459124.8A patent/CN112736896B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105322540A (zh) * | 2015-11-12 | 2016-02-10 | 中国电力科学研究院 | 一种交直流大电网电磁暂态模型的稳态运行方式建立方法 |
| CN109004638A (zh) * | 2018-06-29 | 2018-12-14 | 中国电力科学研究院有限公司 | 一种柔性直流工程机电-电磁混合仿真方法及系统 |
| CN109586275A (zh) * | 2018-10-18 | 2019-04-05 | 清华大学 | 交直流混联系统电磁暂态仿真的快速启动方法与装置 |
| CN109766586A (zh) * | 2018-12-14 | 2019-05-17 | 中国电力科学研究院有限公司 | 一种自动生成大规模电网电磁暂态仿真模型的方法及系统 |
Non-Patent Citations (2)
| Title |
|---|
| Electromechanical Transient and Electromagnetic Transient Hybrid Modeling and Simulation of Back to Back MMC-HVDC Project;Yunlong Fan 等;2018 2nd IEEE Conference on Energy Internet and Energy System Integration (EI2);第1-6页 * |
| MMC-UPFC 电磁–机电混合仿真技术研究;叶小晖 等;电网技术;第43卷(第4期);第1122-1129页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112736896A (zh) | 2021-04-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Kalcon et al. | Small-signal stability analysis of multi-terminal VSC-based DC transmission systems | |
| CN113300383B (zh) | 一种机电暂态建模方法、系统、设备及存储介质 | |
| CN112736896B (zh) | 一种适用于全电磁暂态仿真的三阶段初始化方法及系统 | |
| Mueller et al. | An efficient method of determining operating points of droop-controlled microgrids | |
| Zhu et al. | Reduced-order dynamic model of modular multilevel converter in long time scale and its application in power system low-frequency oscillation analysis | |
| Arunprasanth et al. | Generalized frequency-domain controller tuning procedure for VSC systems | |
| CN106021768A (zh) | 含分布式电源接入的配电网简化建模方法 | |
| CN108110783A (zh) | 一种同步调相机和lcc-hvdc之间的协调控制方法和装置 | |
| Lin et al. | Average modelling of medium frequency DC–DC converters in dynamic studies | |
| Baek et al. | Nonlinear parameter optimization of FACTS controller via real-time digital simulator | |
| Kyesswa et al. | A Matlab-based dynamic simulation module for power system transients analysis in the eASiMOV framework | |
| Kenyon et al. | Criticality of inverter controller order in power system dynamic studies–Case study: Maui island | |
| Almeida et al. | A software tool for the determination of dynamic equivalents of power systems | |
| Abdlrahem et al. | Fixed-order loop shaping robust controller design for parametric models to damp inter-area oscillations | |
| CN113241736A (zh) | 一种用于新能源送出线路保护适应性分析的方法及系统 | |
| Guerreiro et al. | Novel estimation framework for short-circuit current contribution of type IV wind turbines at transient and steady-state of the faults | |
| CN116861576A (zh) | 一种短路故障下双馈风机峰值电流计算方法、装置及介质 | |
| CN107994598B (zh) | 一种交直流协调控制方法及系统 | |
| Kaberere et al. | Comparative analysis and numerical validation of industrial-grade power system simulation tools: Application to small-signal stability | |
| CN108964122B (zh) | 一种风电机组电气模型频率响应特性的验证方法和装置 | |
| Huang et al. | A review on applications of holomorphic embedding methods | |
| Llerins et al. | Modelling and Simulation of Power Systems with Grid-Connected Converters in OpenModelica | |
| Jalili-Marandi et al. | Model-in-the-Loop real-time simulation in phasor domain | |
| Kabiri et al. | DigSILENT modelling of power electronic converters for distributed generation networks | |
| Valdivia et al. | Low-frequency reduced-order modeling approach and implementation of grid emulation in hardware-in-the-loop platforms |
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 |