CN109460575B - 一种基于贝杰龙等值线路的混合仿真接口系统及插值算法 - Google Patents
一种基于贝杰龙等值线路的混合仿真接口系统及插值算法 Download PDFInfo
- Publication number
- CN109460575B CN109460575B CN201811172278.1A CN201811172278A CN109460575B CN 109460575 B CN109460575 B CN 109460575B CN 201811172278 A CN201811172278 A CN 201811172278A CN 109460575 B CN109460575 B CN 109460575B
- Authority
- CN
- China
- Prior art keywords
- electromagnetic
- equivalent
- bergeron
- line
- electromechanical
- 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 67
- 238000004422 calculation algorithm Methods 0.000 title claims abstract description 14
- 238000000034 method Methods 0.000 claims abstract description 10
- 230000008569 process Effects 0.000 claims abstract description 6
- 230000001052 transient effect Effects 0.000 claims description 14
- 238000004458 analytical method Methods 0.000 claims description 7
- 230000005540 biological transmission Effects 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 abstract description 14
- 230000003993 interaction Effects 0.000 abstract description 5
- 238000005516 engineering process Methods 0.000 abstract description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000003990 capacitor Substances 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000002955 isolation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
Classifications
-
- 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
-
- 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
Landscapes
- Engineering & Computer Science (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)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明涉及机电电磁混合仿真技术，具体涉及一种基于贝杰龙等值线路的混合仿真接口系统及插值算法，该值算法包括以下步骤：取混合仿真模型中合适的线路等值成基于贝杰龙等值线路的混合仿真接口模型；获取步骤1所选线路的参数，计算出线路等值阻抗和历史电流源，并赋值给混合仿真接口模型；机电系统在当前时刻t接收电磁系统发送的数据，按照机电系统仿真解算步长∆T解算；电磁系统在当前时刻t接收机电系统t和t‑∆T时刻的历史电流源值进行线性插值；电磁系统按照电磁仿真步长∆t进行解算N个步长，即当电磁系统运算至时刻t=t+∆T时，重复上述解算和交互过程。该法减小了机电系统在电磁解算步长时刻点无法解算引起的解算误差，提高了混合仿真的仿真精度。
Description
技术领域
本发明属于机电电磁混合仿真技术领域，尤其涉及一种基于贝杰龙等值线路的混合仿真接口系统及插值算法。
背景技术
随着大规模新能源、高压直流输电、柔性交流输电、高频电力电子等新型设备接入到电力系统，电力系统变得日趋复杂，区域电网之间的关联性也将大为增加。电网中某个局部条件的改变有可能对整个电网的安全稳定运行产生重要的影响,只孤立考虑局部电网的安全问题已经难以适应现代电力系统发展的要求。
传统的机电暂态仿真程序和电磁暂态仿真程序在分析上述日益复杂的系统特性时,已开始显露出各自的局限性。混合仿真技术综合了机电暂态仿真和电磁暂态仿真各自的优点,对大规模常规电力系统进行机电暂态仿真,对其中重点关注的局部网络或者特定的元件则采用电磁暂态仿真，从而提高了仿真的规模、速度和精度。
发明内容
本发明的目的是提供一种分析庞大复杂电力系统的暂态稳定和动态特性的基于贝杰龙等值线路的混合仿真接口系统及其插值算法。
为实现上述目的，本发明采用的技术方案是：一种基于贝杰龙等值线路的混合仿真接口系统，包括机电系统、电磁系统和基于贝杰龙等值线路的混合仿真接口模型；机电系统包括大规模电网的建模，采用基波、相量、序对系统进行机电暂态分析；电磁系统包括详细电网的建模，采用微分方程进行电磁暂态分析；基于贝杰龙等值线路的混合仿真接口模型包括线路等值阻抗和历史电流源；线路等值阻抗包括线路电阻R和波阻抗Z；历史电流源为关于两端传输延时τ的等值历史电流源。
一种基于贝杰龙等值线路的混合仿真接口系统的插值算法，设机电系统解算步长为ΔT，电磁系统解算步长为Δt，ΔT取值为Δt的整数倍N；包括以下步骤：
步骤1，选取混合仿真模型中合适的线路等值成基于贝杰龙等值线路的混合仿真接口模型；
步骤2，获取步骤1所选线路的参数，计算出线路等值阻抗和历史电流源，并赋值给混合仿真接口模型；
步骤3，机电系统在当前时刻t接收电磁系统发送的数据，按照机电系统仿真解算步长ΔT解算；
步骤4，电磁系统在当前时刻t接收机电系统t和t-ΔT时刻的历史电流源值进行线性插值；
步骤5，电磁系统按照电磁仿真步长Δt进行解算N个步长，即当电磁系统运算至时刻t＝t+ΔT时，重复上述解算和交互过程。
本发明的有益效果：基于贝杰龙等值线路的混合仿真接口模型实现了机电电磁混合仿真，不仅继承了混合仿真的优良特性，而且减小了机电系统在电磁解算步长时刻点无法解算引起的解算误差，利用插值算法，进一步提高了混合仿真的精度。该插值算法对研究庞大复杂电力系统的暂态稳定和动态特性分析具有较强的工程实用价值。
附图说明
图1为本发明一个实施例涉及的基于贝杰龙等值线路的混合仿真接口系统；
图2为本发明一个实施例涉及的基于贝杰龙等值线路的混合仿真接口模型；
图3为本发明一个实施例涉及的基于贝杰龙等值线路的混合仿真接口系统的插值算法。
具体实施方式
下面结合附图对本发明的实施方式进行详细描述。
本实施例采用的技术方案之一是：
一种基于贝杰龙等值线路的混合仿真接口系统，包括机电系统、电磁系统和基于贝杰龙等值线路的混合仿真接口模型。
机电系统包括大规模电网的建模，采用基波、相量、序对系统进行机电暂态的分析。
电磁系统包括详细电网的建模，采用微分方程进行电磁暂态的分析。
基于贝杰龙等值线路的混合仿真接口模型包括线路等值阻抗和历史电流源。线路等值阻抗由线路电阻R和波阻抗Z组成。历史电流源是指关于两端传输延时τ的等值历史电流源。
本实施例采用的技术方案之二是：在机电电磁混合仿真系统中，采用贝杰龙等值线路作为机电电磁混合仿真接口模型，提出了一种线性插值算法。一种基于贝杰龙等值线路的混合仿真接口系统的插值算法。设机电系统解算步长为ΔT，电磁系统解算步长为Δt，ΔT取值为Δt的整数倍N。
具体步骤如下：
S1，选取混合仿真模型中合适的线路等值成基于贝杰龙等值线路的混合仿真接口模型；
S2，获取步骤1所选线路的参数，计算出线路等值阻抗和历史电流源，并赋值给混合仿真接口模型；
S3，机电系统在当前时刻t接收电磁系统发送的数据，按照机电系统仿真解算步长ΔT解算；
S4，电磁系统在当前时刻t接收机电系统t和t-ΔT时刻的历史电流源值进行线性插值；
S5，电磁系统按照电磁仿真步长Δt进行解算N个步长，即当电磁系统运算至时刻t＝t+ΔT时，重复上述解算和交互过程。
如图1所示，为本实施的基于贝杰龙线路的混合仿真接口系统。
左边的代表机电系统，右边的代表电磁系统。机电系统节点m和电磁系统节点k通过贝杰龙线路相互连接。机电系统的节点电压为Vm(t)，节点流入电流为imk(t)。电磁系统的节点电压为Vk(t)，节点流入电流为ikm(t)。贝杰龙线路的线路长度为d,线路电阻为R,线路电感为L。为了引入贝杰龙线路，引入了线路中可能非实际存在的电容C。将线路长度均分成d/Δd段，单位线路长度的电阻为R0,电感为L0,电容为C0。如图1所示，每段线路长度取为Δd，由电阻串联电感，并联电容构成。
如图2所示，为本实施的基于贝杰龙等值线路的混合仿真接口模型。利用该模型，机电系统和电磁系统分别等值成对应的受控历史电流源和线路等值阻抗。具体实施步骤如下：
步骤一，选取混合仿真模型中合适的线路等值成基于贝杰龙等值线路的混合仿真接口模型；
步骤二，获取步骤一所选线路的参数，计算出线路等值阻抗和历史电流源，并赋值给混合仿真接口模型；
步骤三，机电系统在当前时刻t接收电磁系统发送的数据，按照机电系统仿真解算步长ΔT解算；
步骤四，电磁系统在当前时刻t接收机电系统t和t-ΔT时刻的历史电流源值进行线性插值；
步骤五，电磁系统按照电磁仿真步长Δt进行解算N个步长，即当电磁系统运算至时刻t＝t+ΔT时，重复上述解算和交互过程。
如图3所示，本实施例涉及的基于贝杰龙等值线路的混合仿真接口系统的插值算法。具体实施步骤如下：
(1)机电系统在当前时刻t接收电磁系统发送的数据，按照机电系统仿真解算步长ΔT解算；
(2)电磁系统在当前时刻t接收机电系统t和t-ΔT时刻的历史电流源值进行线性插值；
(3)电磁系统按照电磁仿真步长Δt进行解算N个步长，即当电磁系统运算至时刻t＝t+ΔT时，重复上述解算和交互过程。
应当理解的是，本说明书未详细阐述的部分均属于现有技术。
虽然以上结合附图描述了本发明的具体实施方式，但是本领域普通技术人员应当理解，这些仅是举例说明，可以对这些实施方式做出多种变形或修改，而不背离本发明的原理和实质。本发明的范围仅由所附权利要求书限定。
Claims (2)
1.一种基于贝杰龙等值线路的混合仿真接口系统，其特征是，包括机电系统、电磁系统和基于贝杰龙等值线路的混合仿真接口模型，机电系统节点m和电磁系统节点k通过贝杰龙线路相互连接；机电系统包括大规模电网的建模，采用基波、相量、序对系统进行机电暂态分析；电磁系统包括详细电网的建模，采用微分方程进行电磁暂态分析；基于贝杰龙等值线路的混合仿真接口模型包括线路等值阻抗和历史电流源；线路等值阻抗包括线路电阻R和波阻抗Z；历史电流源为关于两端传输延时τ的等值历史电流源。
2.用于权利要求1所述基于贝杰龙等值线路的混合仿真接口系统的插值算法，机电系统和电磁系统分别等值成对应的受控历史电流源和线路等值阻抗；设机电系统解算步长为，电磁系统解算步长为/>，/>取值为/>的整数倍N；其特征是，包括以下步骤：
步骤1，选取混合仿真模型中合适的线路等值成基于贝杰龙等值线路的混合仿真接口模型；
步骤2，获取步骤1所选线路的参数，计算出线路等值阻抗和历史电流源，并赋值给混合仿真接口模型；
步骤3，机电系统在当前时刻t接收电磁系统发送的数据，按照机电系统仿真解算步长解算；
步骤4，电磁系统在当前时刻t接收机电系统t和时刻的历史电流源值进行线性插值；
步骤5，电磁系统按照电磁仿真步长进行解算N个步长，即当电磁系统运算至时刻时，重复上述解算和交互过程。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811172278.1A CN109460575B (zh) | 2018-10-09 | 2018-10-09 | 一种基于贝杰龙等值线路的混合仿真接口系统及插值算法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811172278.1A CN109460575B (zh) | 2018-10-09 | 2018-10-09 | 一种基于贝杰龙等值线路的混合仿真接口系统及插值算法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109460575A CN109460575A (zh) | 2019-03-12 |
| CN109460575B true CN109460575B (zh) | 2023-10-31 |
Family
ID=65607400
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811172278.1A Active CN109460575B (zh) | 2018-10-09 | 2018-10-09 | 一种基于贝杰龙等值线路的混合仿真接口系统及插值算法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109460575B (zh) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113128073B (zh) * | 2021-05-13 | 2022-07-08 | 清鸾科技(成都)有限公司 | 一种多仿真器联合仿真方法和系统 |
| CN115085190B (zh) * | 2022-06-16 | 2024-08-23 | 内蒙古电力(集团)有限责任公司内蒙古电力科学研究院分公司 | 等效贝瑞隆模型的短线路建模方法 |
| CN116738703A (zh) * | 2023-06-07 | 2023-09-12 | 国网江苏省电力有限公司电力科学研究院 | 短线路配电网的仿真方法、装置、存储介质及计算机设备 |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102592003A (zh) * | 2011-12-15 | 2012-07-18 | 天津大学 | 一种用于机电与电磁暂态混合仿真的数据交换方法 |
| CN102708225A (zh) * | 2012-04-13 | 2012-10-03 | 南方电网科学研究院有限责任公司 | 一种交直流大电网电磁暂态仿真的分片调试方法 |
| CN103336460A (zh) * | 2013-05-21 | 2013-10-02 | 中国南方电网有限责任公司电网技术研究中心 | 一种电磁机电混合实时仿真数字化接口的控制方法 |
| CN105205244A (zh) * | 2015-09-14 | 2015-12-30 | 国家电网公司 | 基于机电-电磁混合仿真技术的合环操作仿真系统 |
| CN105738770A (zh) * | 2016-04-01 | 2016-07-06 | 昆明理工大学 | 一种基于故障行波沿线分布特性的线缆混合线路单端行波测距方法 |
| CN106681168A (zh) * | 2015-11-10 | 2017-05-17 | 国网智能电网研究院 | 一种用于mmc换流阀数模混合仿真接口网络的解耦方法 |
| CN106777827A (zh) * | 2017-01-24 | 2017-05-31 | 中国电力科学研究院 | 一种机电‑电磁混合仿真方法及系统 |
| CN107169244A (zh) * | 2017-06-27 | 2017-09-15 | 武汉大学 | 一种机电‑电磁暂态混合仿真接口系统及方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101957398B (zh) * | 2010-09-16 | 2012-11-28 | 河北省电力研究院 | 一种基于机电与电磁暂态混合仿真技术检测计算电网一次时间常数的方法 |
-
2018
- 2018-10-09 CN CN201811172278.1A patent/CN109460575B/zh active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102592003A (zh) * | 2011-12-15 | 2012-07-18 | 天津大学 | 一种用于机电与电磁暂态混合仿真的数据交换方法 |
| CN102708225A (zh) * | 2012-04-13 | 2012-10-03 | 南方电网科学研究院有限责任公司 | 一种交直流大电网电磁暂态仿真的分片调试方法 |
| CN103336460A (zh) * | 2013-05-21 | 2013-10-02 | 中国南方电网有限责任公司电网技术研究中心 | 一种电磁机电混合实时仿真数字化接口的控制方法 |
| CN105205244A (zh) * | 2015-09-14 | 2015-12-30 | 国家电网公司 | 基于机电-电磁混合仿真技术的合环操作仿真系统 |
| CN106681168A (zh) * | 2015-11-10 | 2017-05-17 | 国网智能电网研究院 | 一种用于mmc换流阀数模混合仿真接口网络的解耦方法 |
| CN105738770A (zh) * | 2016-04-01 | 2016-07-06 | 昆明理工大学 | 一种基于故障行波沿线分布特性的线缆混合线路单端行波测距方法 |
| CN106777827A (zh) * | 2017-01-24 | 2017-05-31 | 中国电力科学研究院 | 一种机电‑电磁混合仿真方法及系统 |
| CN107169244A (zh) * | 2017-06-27 | 2017-09-15 | 武汉大学 | 一种机电‑电磁暂态混合仿真接口系统及方法 |
Non-Patent Citations (1)
| Title |
|---|
| 电磁-机电暂态混合仿真接口交互信息限制性分析;陈鹏伟;陶顺;杨洋;肖湘宁;;电工电能新技术(第05期);100-103 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109460575A (zh) | 2019-03-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109460575B (zh) | 一种基于贝杰龙等值线路的混合仿真接口系统及插值算法 | |
| CN106372339B (zh) | 电力电子化电力系统的多速率仿真方法及装置 | |
| CN112100962B (zh) | 一种mmc的电磁暂态仿真方法及系统 | |
| CN110232207B (zh) | 一种基于人工神经网络的电磁干扰滤波器设计方法 | |
| CN109543339B (zh) | 三电平换流器的定导纳建模与实时仿真方法 | |
| CN105850015A (zh) | 具有lc滤波器的电气转换器的控制方法 | |
| CN103336866B (zh) | 一种电磁暂态仿真中含负电阻支路的处理方法 | |
| CN110362937B (zh) | 一种模块化多电平换流器电磁暂态仿真方法及系统 | |
| CN109948185B (zh) | 一种电力系统的解耦仿真方法 | |
| CN113158447B (zh) | 一种大步长移频电磁暂态仿真方法及系统 | |
| CN105260516A (zh) | 一种含开关特性子网络的电磁暂态仿真方法 | |
| CN109614687B (zh) | 两电平桥式换流器的定导纳建模与实时仿真方法 | |
| CN104375876B (zh) | 一种输入量突变情况下的0+误差免疫电磁暂态仿真方法 | |
| CN108959671A (zh) | 半桥和全桥型模块化多电平换流器的实时仿真建模方法 | |
| CN106886617B (zh) | 一种含多vsc的多速率电磁暂态分网方法 | |
| CN110705035A (zh) | 大容量电力电子系统的离散状态事件驱动解耦仿真方法 | |
| CN108917917A (zh) | 一种基于稀疏分解的断路器机械振动信号去噪方法 | |
| CN115618701A (zh) | Mmc的低维导纳电磁暂态建模仿真方法、装置及相关设备 | |
| CN111709209B (zh) | 基于支路指数积分形式的电磁暂态仿真方法及系统 | |
| CN112199914A (zh) | 一种电力电子开关恒导纳模型建立方法及系统 | |
| CN112434411A (zh) | 采用变阶变步长3s-dirk算法的电磁暂态仿真方法 | |
| CN106598911B (zh) | 一种电路网络状态方程的列写方法 | |
| CN112084674B (zh) | 一种基于半隐式松弛的电磁暂态并行仿真方法及系统 | |
| CN114662439B (zh) | 一种vsc换流器的快速仿真方法及系统 | |
| CN116011372B (zh) | 一种基于分数阶电感的双有源桥变换器分数阶建模方法 |
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 |