CN105243232A - 将场分析集成到电网络求解的电磁暂态仿真方法及系统 - Google Patents

将场分析集成到电网络求解的电磁暂态仿真方法及系统 Download PDF

Info

Publication number
CN105243232A
CN105243232A CN201510740268.3A CN201510740268A CN105243232A CN 105243232 A CN105243232 A CN 105243232A CN 201510740268 A CN201510740268 A CN 201510740268A CN 105243232 A CN105243232 A CN 105243232A
Authority
CN
China
Prior art keywords
equipment
model
finite element
analysis
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.)
Granted
Application number
CN201510740268.3A
Other languages
English (en)
Other versions
CN105243232B (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.)
State Grid Corp of China SGCC
Electric Power Research Institute of State Grid Hubei Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
Electric Power Research Institute of State Grid Hubei Electric Power Co Ltd
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 State Grid Corp of China SGCC, Electric Power Research Institute of State Grid Hubei Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN201510740268.3A priority Critical patent/CN105243232B/zh
Publication of CN105243232A publication Critical patent/CN105243232A/zh
Application granted granted Critical
Publication of CN105243232B publication Critical patent/CN105243232B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • 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

Landscapes

  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明公开了一种将场分析集成到电网络求解的电磁暂态仿真方法及系统,其特征在于:引入设备特性方程,建立表征设备运行健康状态的模型和相应判断设备健康状态的参考值;建立交直流区域大电网的电磁暂态仿真模型,计算结果;将计算结果与参考值比较,判断设备运行是否正常;建立通讯,通过开发的专用程序,将仿真结果加载到有限元模型中;基于场分析,通过有限元方法计算设备的详细电气特性;专家系统判断设备的运行状态,确定设备正常与否,必要时发出警告提示安排检修;对模型进行校验、修正。本发明实现了在系统进行电网络分析的同时,对设备的运行状态进行实时或离线的分析,有效指导电网设备状态评价、检修及故障分析。

Description

将场分析集成到电网络求解的电磁暂态仿真方法及系统
技术领域
本发明涉及电力设备电磁暂态仿真领域,具体涉及一种将场分析集成到电网络求解的电磁暂态仿真方法及系统。
背景技术
目前电力设备特性分析方法主要有基于场分析的有限元分析方法和基于电网络的电磁暂态分析方法。
基于场分析的有限元分析方法,能够准确仿真电气设备在加载电流、电压的电、磁、温度、形变特性,对于各种复杂模型和边界条件、不均匀或非线性材料的情况,都具有较高的处理能力,能够直接反映设备的物理本质和运行工况。但该方法需要建立复杂的物理模型,运行时间较长,当考虑的物理过程较复杂时,计算量大,程序调试较慢且不易与实时的电网的一次设备网络建立联系,只能针对某个孤立具体的设备深入分析。常见的有限元分析软件主要有ANSYS,COMSOL,JMAG,FLUX等。
基于电网络分析的电磁暂态方法,目前常见的国内外广泛使用的仿真软件主要有:国际公认的电磁暂态仿真标准程序EMTP,加拿大曼尼托巴直流研究中心开发的PSCAD/EMTDC,中国电科院开发的EMTPE和PSASP,美国电力技术公司开发的PSS/E,加拿大RTDS公司开发的RTDS以及加拿大魁北克TET公司开发的HYPERSIM等。其中,EMTDC、RTDS、HYPERSIM等均是基于EMTP内核开发,与EMTP有相同的底层运算原理,只是软件的应用侧重不同。对于全球最普遍使用的EMTP软件,该方法主要是面向基于电网络求解的电磁暂态分析,研究电力系统稳态潮流及发生故障和过电压时,一次与二次系统的暂态过程。但该方法由于是基于电路计算,虽能得到通过设备的电压、电流、功率,但无法考虑电力设备内部的具体电磁特性,难以据此准确判断设备的运行状态。
目前我国电力检修工作多采用状态检修,该方法是以电力设备当前的实际工作状况为依据,通过状态监测、分析手段,掌控设备的现实运行状态,对故障或隐患的早期征兆、故障部位、故障严重程度及发展趋势做出判断,从而有效指导设备检修与更换,大大降低了检修的成本,提高了检修的精确度。但该方法需要对设备的具体运行状态及电磁特性有较清楚的了解。
为了解决上述问题,本专利提出了将场分析的有限元方法与基于电网络分析的电磁暂态方法相结合,对电力设备的具体特性进行仿真,在求解电网络的同时对设备的特性进行分析,有效指导电力设备的状态检修。
发明内容
针对上述现有技术上的不足,本发明提供一种将场分析集成到电网络求解的电磁暂态仿真方法。
一种将场分析集成到电网络求解的电磁暂态仿真方法,其特征在于:所述方法包括以下步骤:
步骤1):预先对电磁暂态仿真的设备模型进行二次开发,引入设备特性方程,建立表征设备运行健康状态的模型和相应判断设备健康状态的参考值;
步骤2):建立交直流区域大电网的电磁暂态仿真模型,基于电路分析求解系统网络,得到电压、电流、功率的计算结果;
步骤3):将电压、电流计算结果与设备运行健康状态模型中预设的设备健康运行参考值比较,判断设备运行是否正常;
步骤4):若计算结果超过参考值范围或需要进行设备深度分析,则建立通讯,通过开发的专用程序,将上述仿真结果加载到有限元模型中;
步骤5):基于场分析,通过有限元方法计算设备的详细电气特性;
步骤6):专家系统根据场分析计算结果判断设备的运行状态,确定设备正常与否,必要时发出警告提示安排检修;
步骤7):将场分析结果结合现场设备状态,反馈到步骤1)所述的设备运行健康状态模型,对模型进行校验、修正。
所述步骤1)到步骤3)为电力系统网络实时仿真,进行硬件在环(HIL)测试;步骤4)到步骤7)为针对具体电力设备的离线仿真,基于场分析。
实时仿真基于电路导纳矩阵求解,主要使用Hypersim软件进行电网仿真,EMTP-RV软件对其进行局部校验及补充;离线仿真基于有限元方法针对设备进行电磁场、温度场、应力及形变场分析,使用但不限于ANSYS,COMSOL,JMAG,FLUX有限元软件。
所述步骤1)中,表征设备运行健康状态的模型是指,在原Hypersim软件中基于EMTP算法的元件模型基础上进行二次开发,对设备模型加入逻辑算法,将反映电气设备电磁、热等特性的方程通过逻辑算法实现,并添加到原EMTP模型中,由计算出的设备的电压、电流、功率,通过方程间接反映出设备的运行特性,如变压器铁芯损耗、变压器线圈温升,从而快速判断设备健康状态,表征设备运行健康状态的模型及相应正常工作参考值的确定通过离线有限元程序结合现场测量参数校验。
所述表征设备运行健康状态的变压器铁芯损耗、变压器线圈温升模型,其主要特性按照下述方法计算、表征,并建立正常工作状态下设备特性的参考值:
变压器铁芯损耗使用下述方法计算:
1 T ∫ 0 T dP I r o n l o s s e s ( t ) d t = k h B m 2 fk f + 1 T ∫ 0 T { σd 2 1 12 [ d B d t ( t ) ] 2 + k e [ d B d t ( t ) ] 3 2 } k f d t - - - ( 1 )
其中:Pironloss是铁芯损耗,Bm是铁心区域最大磁通密度,f是电流频率,σ为铁心材料的电导率,d为铁心厚度,T为周期,kh为涡流损耗系数,ke为附加损耗系数,kf为硅钢片叠片系数,kh,ke,kf由离线有限元分析结果结合现场设备监测数据得到,Bm采用有限元计算得到铁心区域最大磁通密度Bm,然后将其与原、副边电流建立非线性映射关系,通过变压器原、副边电流作为状态量反映铁心区域最大磁通密度Bm,按照步骤2)计算原、副边电流,通过原、副边电流作为状态量反映变压器铁芯损耗,然后建立正常工作状态变压器铁芯损耗参考值作为设备正常与否的判据;
与上述铁心计算方法类似,所述表征设备运行健康状态的模型中,变压器线圈温升与原、副边电流也建立映射关系,温升按照下述方法计算:
Q l o s s = I 2 R t = c m Δ T Q e x c h a n g e = t h A ( Δ T + 40 - T o i l ) = Q l o s s - - - ( 2 )
Δ T = I 2 R h A + T o i l - 40 - - - ( 3 )
其中:ΔT是变压器线圈温升,I是原边电流,R是原边电阻,h是变压器线圈的对流换热系数,A是变压器线圈与油的接触面积,Toil是变压器油温,40代表“标准环境温度”,据此用原边电流的数值表征变压器线圈温升,建立正常工作状态变压器线圈温升参考值及其对应的电流值,将其作为设备正常与否的判据;
类似的,所述表征设备运行健康状态的模型中,按照步骤2)中计算出断路器处暂态过电压,建立正常工作状态断路器处能承受的暂态过电压参考值,将其作为设备正常与否的判据。
所述步骤4)中,专用程序使用MATLAB软件编写,程序将Hypersim软件计算的电压、电流结果从ScopeView中导出,转换成*.mat或*.txt格式,对于直流状态量,专用程序通过逻辑过滤掉信号尖峰,再对信号积分取平均值得到结果;对于交流状态量,专用程序对信号进行滤波,保留信号基波及直流分量两部分作为结果,最后,将结果作为载荷加载到有限元模型中,通过运行MATLAB专用程序可实现自动加载、调用、完成有限元分析。
所述步骤5)中,有限元分析模型针对系统各电气主设备预先建好,离线运行,为保证计算速度,主要使用2D有限元模型进行频域、暂态仿真,3D有限元模型仅开展频域分析。
一种将场分析集成到电网络求解的电磁暂态仿真系统,其特征在于:包括进行有限元分析的有限元分析计算机以及对电网络进行建模仿真的电网仿真机,电网仿真机对电网络进行仿真计算并将仿真计算结果数据通过通信系统输入到数据处理系统中,数据处理系统对仿真计算结构数据进行处理后通过通信系统输送到有限元计算机内进行基于场分析的有限元计算。
所述有限元分析计算机和电网仿真机分别连接有显示器装置以及输入设备。
所述电网仿真机为安装Hypersim软件的计算机
本发明的技术效果体现在:
本发明针对电力设备电磁暂态仿真中存在的问题,提出了一种将场分析集成到电网络求解的电力设备电磁暂态仿真的方法。该方法结合了基于场分析的有限元仿真和基于电网络电磁暂态仿真方法,能够表征电力设备健康运行状态模型,实现电力设备的实时和离线仿真,并在离线状态通过有限元反向校验设备健康状态模型,并识别设备早期的故障征兆,对设备故障部位,故障严重程度及发展趋势做出判断,为电力设备的检修提供方便,保证电力设备安全、稳定、健康运行。
附图说明
图1是将场分析集成到电网络求解的电磁暂态仿真方法的流程图;
图2是建立表征设备运行健康状态模型的流程图;
图3是有限元分析流程图;
图4是将场分析集成到电网络求解的电磁暂态仿真系统结构图;
图5是仿真得到的变压器原边的电压电流及其励磁电流波形界面图;
图6是变压器铁心的磁密云图。
具体实施方式
以下将结合本发明中的附图,对本发明中的技术方案进行完整地描述。
一种将场分析集成到电网络求解的电磁暂态仿真系统,其特征在于:包括进行有限元分析的有限元分析计算机1以及对电网络进行建模仿真的电网仿真机4,电网仿真机4对电网络进行仿真计算并将仿真计算结果数据通过通信系统2输入到数据处理系统3中,数据处理系统3对仿真计算结构数据进行处理后通过通信系统2输送到有限元计算机1内进行基于场分析的有限元计算。
所述有限元分析计算机1和电网仿真机4分别连接有显示器装置以及输入设备。
所述电网仿真机4为安装Hypersim软件的计算机
本发明实施例提供一种将场分析集成到电网络求解的高压直流输电系统直流偏磁下的变压器设备电磁暂态仿真方法,实现方式包括以下步骤:
步骤1:预先建立表征变压器运行健康状态的模型并设置变压器正常运行状态下不同电流或电压值对应的励磁电流、铁心损耗、电磁场、热场等参考值。
为建立能够反映变压器特性的精确模型,本发明的实施方式是结合逻辑算法对Hypersim软件模型库中基于EMTP算法的变压器元件模型进行二次开发,图2为流程图,具体实现方式如下:
步骤1-1:确定变压器设备在Hypersim软件库中的原EMTP模型,该模型可以直接从Hypersim软件库中拖出使用;
步骤1-2:引入变压器设备的电流、电压值与电磁场、热场、形变特性的关系方程,如表征变压器的铁芯损耗,按照下述公式计算:
1 T ∫ 0 T dP I r o n l o s s e s ( t ) d t = k h B m 2 fk f + 1 T ∫ 0 T { σd 2 1 12 [ d B d t ( t ) ] 2 + k e [ d B d t ( t ) ] 3 2 } k f d t
其中:Pironloss是铁芯损耗,Bm是铁心区域最大磁通密度,f是电流频率,σ为铁心材料的电导率,d为铁心厚度,T为周期,kh为涡流损耗系数,ke为附加损耗系数,kf为硅钢片叠片系数。
步骤1-3:通过设备特性方程计算得到变压器不同电流、电压值对应的电磁场、热场、形变参考值,建立电流、电压与电气设备具体特性(如铁芯损耗,温升等)的联系;
例如步骤1-2中要计算变压器的铁芯损耗,kh,ke,kf由离线有限元分析结果结合现场设备监测数据得到,Bm则根据采样到的变压器原、副边电流值基于有限元法计算对应的铁心区域最大磁通密度Bm,从而得到了变压器原、副边电流与最大磁通密度Bm(或铁芯损耗)的对应关系。
步骤1-4:根据变压器正常运行时特性(如温升、铁芯损耗等)的参数范围,对应到步骤1-3所述的电流、电压与设备特性的关系,确定反映设备特性正常的电压、电流参考值,以此作为设备正常与否的判据。比如型号为ODFS-334MVA/500kV的变压器,其原边额定相电压为317.54kV,原边额定电流为1051.8A,励磁电流为0.037%,0.39A,额定空载损耗(铁耗)72.4kW;
步骤1-5:将步骤1-4、步骤1-3得到的设备正常工作的判据及电流、电压与设备特性的关系集成到步骤1-1所述设备的原EMTP模型中,得到表征设备运行健康状态的模型。
步骤2:实时仿真,进行网络计算:
在Hypersim软件中编辑模式下对该变压器所在的大电网包括整流侧换流站、直流输电线路、逆变侧换流站进行建模,所要分析的电力设备采用步骤1中二次开发的模型;在仿真控制模式下选择在线仿真模式,通过Hypersim软件对设备所在区域电网进行仿真,建立节点导纳矩阵,采用节点导纳方程求解网络变量。Hypersim自动将任务生成C代码分配给各并行处理器,实时仿真计算得到目标电网设备节点的电压、电流值。在变压器工作在原边电流为479.97A,原边相电压为302.26kV时(如图5所示),从变压器中性点涌入三相变压器共5A的直流偏磁电流,此时通过所述表征设备运行健康状态的模型算出各相变压器励磁电流畸变,峰值达到7.7A,如图5所示,铁心损耗109kW。
步骤3:初次判断电力设备是否正常
将步骤2中得到的计算电流或电压值与步骤1中预先设置的设备健康运行下的参考参数值进行比较:如果在设备正常运行的参考值范围之内,则判断设备为正常工作状态,无需进行进一步分析与检修;如果超过设备正常运行的参考值范围,则需要进一步分析研究设备在该计算电流、电压值对应的电磁场、热场以及形变详细特性,再次判断设备是否运行正常。本实施例中变压器原边电压电流均在正常范围内,励磁电流远远超过正常值,铁心损耗也较额定值有所增加,需要进一步分析研究变压器在该电流、电压值下的详细特性。
步骤4:若设备电压、电流参数超过参考值或需要进行设备深度分析,则将上述仿真结果加载到有限元模型中,进行离线分析。
离线仿真的目的是对超过设备正常运行参考值范围的情况,进一步判断电力设备是否需要安排检修。仿真需要在有限元分析软件中进行,使用但不限于ANSYS,COMSOL,JMAG,FLUX软件。
仿真首先是建立实时仿真结果与有限元仿真软件之间的通讯,在此使用MATLAB软件编写专用程序。专用程序将Hypersim软件计算的电压、电流结果从ScopeView中导出,转换成*.mat或*.txt格式,并将参数自动导入到有限元程序计算,并导出有限元计算结果。对于直流状态量,专用程序通过逻辑过滤掉信号尖峰,再对信号积分取平均值得到结果;对于交流状态量,专用程序对信号进行滤波,保留信号基波及直流分量两部分叠加作为结果。
步骤5:对电力主设备进行有限元仿真分析。图3为有限元分析流程,具体实时方式如下步骤:
步骤5-1:根据不同厂家电力设备的几何尺寸,并遵循有限元软件中建模规则,预先建立好系统各电气主设备的有限元模型。有限元建模可以综合考虑仿真精度和设备结构,对设备的物理结构进行适当的简化处理,对仿真结果影响不大的模型可以忽略不计。对于不同的工况的仿真,只要使用的电力设备型号相同,则不需要对有限元模型进行反复修改,只需通过步骤4所述专用程序改变需要导入的参数。
步骤5-2:分配材料属性,例如设备材料主要电磁、热、结构特性参数:导热系数、材料密度、磁导率、电导率,相对介电参数等,这些也是在建模之后已经设置好。
步骤5-3:划分网格,施加载荷和边界条件。本发明中施加载荷的方式也是通过调用步骤4所述的基于MATLAB开发的专用程序,运行该程序,可以将转换为*.mat或*.txt格式的实时仿真电流、电压值作为载荷自动调用并加载到有限元模型中。边界条件的设置按照设备现场运行工况设置。
步骤5-4:根据需要添加需要进行的求解类型:频域分析或是暂态分析;计算电磁特性还是热特性等,进行求解。
步骤5-5:观察设备的电磁场或热场分析结果并导出。根据步骤1-4的电磁暂态仿真结果,变压器原边电流为479.97A,原边相电压为302.26kV,从变压器中性点涌入三相变压器共5A的直流偏磁电流,算得变压器铁心的磁密云图如图6所示。
步骤6:专家根据仿真云图及其他相关有限元后处理结果再次判断电力设备是否健康,确定设备正常与否,必要时发出警告提示安排检修。图6所示在直流偏磁发生时,铁心最大磁密达到1.74T,超过1.7T,已趋于饱和,需要加强监视。
步骤7:校验、修正设备运行健康状态模型
该实施方式是对实时仿真计算得到的电流、电压值,导入有限元仿真软件计算得到的电磁场、热场与步骤1中预先设置的参考值进行对比,对比较结果有差异的设备的健康状态模型进行修改,以校验步骤1设置的参考值及对应关系的正确性。
本发明针对电力设备电磁暂态仿真中存在的问题,提出了一种将场分析集成到电网络求解的电力设备电磁暂态仿真的方法。该方法结合了基于场分析的有限元仿真和基于电网络电磁暂态仿真方法,能够表征电力设备健康运行状态模型,实现电力设备的实时和离线仿真,并在离线状态通过有限元反向校验设备健康状态模型,并识别设备早期的故障征兆,对设备故障部位,故障严重程度及发展趋势做出判断,为电力设备的检修提供方便,保证电力设备安全、稳定、健康运行。

Claims (10)

1.一种将场分析集成到电网络求解的电磁暂态仿真方法,其特征在于:所述方法包括以下步骤:
步骤1):预先对电磁暂态仿真的设备模型进行二次开发,引入设备特性方程,建立表征设备运行健康状态的模型和相应判断设备健康状态的参考值;
步骤2):建立交直流区域大电网的电磁暂态仿真模型,基于电路分析求解系统网络,得到电压、电流、功率的计算结果;
步骤3):将电压、电流计算结果与设备运行健康状态模型中预设的设备健康运行参考值比较,判断设备运行是否正常;
步骤4):若计算结果超过参考值范围或需要进行设备深度分析,则建立通讯,通过开发的专用程序,将上述仿真结果加载到有限元模型中;
步骤5):基于场分析,通过有限元方法计算设备的详细电气特性;
步骤6):专家系统根据场分析计算结果判断设备的运行状态,确定设备正常与否,必要时发出警告提示安排检修;
步骤7):将场分析结果结合现场设备状态,反馈到步骤1)所述的设备运行健康状态模型,对模型进行校验、修正。
2.根据权利要求1所述的将场分析集成到电网络求解的电磁暂态仿真方法,其特征在于:所述步骤1)到步骤3)为电力系统网络实时仿真,进行硬件在环(HIL)测试;步骤4)到步骤7)为针对具体电力设备的离线仿真,基于场分析。
3.根据权利要求2所述的将场分析集成到电网络求解的电磁暂态仿真方法,其特征在于:实时仿真基于电路导纳矩阵求解,主要使用Hypersim软件进行电网仿真,EMTP-RV软件对其进行局部校验及补充;离线仿真基于有限元方法针对设备进行电磁场及温度场分析,使用但不限于ANSYS,COMSOL,JMAG,FLUX有限元软件。
4.根据权利要求1所述的将场分析集成到电网络求解的电磁暂态仿真方法,其特征在于:所述步骤1)中,表征设备运行健康状态的模型是指,在原Hypersim软件中基于EMTP算法的元件模型基础上进行二次开发,对设备模型加入逻辑算法,将反映电气设备电磁、热等特性的方程通过逻辑算法实现,并添加到原EMTP模型中,由计算出的设备的电压、电流、功率,通过方程间接反映出设备的运行特性,如变压器铁芯损耗、变压器线圈温升,从而快速判断设备健康状态,表征设备运行健康状态的模型及相应正常工作参考值的确定通过离线有限元程序结合现场测量参数校验。
5.根据权利要求4所述的将场分析集成到电网络求解的电磁暂态仿真方法,其特征在于:所述表征设备运行健康状态的变压器铁芯损耗、变压器线圈温升模型,其主要特性按照下述方法计算、表征,并建立正常工作状态下设备特性的参考值:
变压器铁芯损耗使用下述方法计算:
1 T ∫ 0 T dP I r o n l o s s ( t ) d t = k h B m 2 fk f + 1 T ∫ 0 T { σd 2 1 12 [ d B d t ( t ) ] 2 + k e [ d B d t ( t ) ] 3 2 } k f d t - - - ( 1 )
其中:Pironloss是铁芯损耗,Bm是铁心区域最大磁通密度,f是电流频率,σ为铁心材料的电导率,d为铁心厚度,T为周期,kh为涡流损耗系数,ke为附加损耗系数,kf为硅钢片叠片系数,kh,ke,kf由离线有限元分析结果结合现场设备监测数据得到,Bm采用有限元计算得到铁心区域最大磁通密度Bm,然后将其与原、副边电流建立非线性映射关系,通过变压器原、副边电流作为状态量反映铁心区域最大磁通密度Bm,按照步骤2)计算原、副边电流,通过原、副边电流作为状态量反映变压器铁芯损耗,然后建立正常工作状态变压器铁芯损耗参考值作为设备正常与否的判据;
与上述铁心计算方法类似,所述表征设备运行健康状态的模型中,变压器线圈温升与原、副边电流也建立映射关系,温升按照下述方法计算:
Q l o s s = I 2 R t = c m Δ T Q e x c h a n g e = t h A ( Δ T + 40 - T o i l ) = Q l o s s - - - ( 2 )
Δ T = I 2 R h A + T o i l - 40 - - - ( 3 )
其中:ΔT是变压器线圈温升,I是原边电流,R是原边电阻,h是变压器线圈的对流换热系数,A是变压器线圈与油的接触面积,Toil是变压器油温,40代表“标准环境温度”,据此用原边电流的数值表征变压器线圈温升,建立正常工作状态变压器线圈温升参考值及其对应的电流值,将其作为设备正常与否的判据;
类似的,所述表征设备运行健康状态的模型中,按照步骤2)中计算出断路器处暂态过电压,建立正常工作状态断路器处能承受的暂态过电压参考值,将其作为设备正常与否的判据。
6.根据权利要求1所述的将场分析集成到电网络求解的电磁暂态仿真方法,其特征在于:所述步骤4)中,专用程序使用MATLAB软件编写,程序将Hypersim软件计算的电压、电流结果从ScopeView中导出,转换成*.mat或*.txt格式,对于直流状态量,专用程序通过逻辑过滤掉信号尖峰,再对信号积分取平均值得到结果;对于交流状态量,专用程序对信号进行滤波,保留信号基波及直流分量两部分作为结果,最后,将结果作为载荷加载到有限元模型中,通过运行MATLAB专用程序可实现自动加载、调用、完成有限元分析。
7.根据权利要求1所述的将场分析集成到电网络求解的电磁暂态仿真方法,其特征在于:所述步骤5)中,有限元分析模型针对系统各电气主设备预先建好,离线运行,为保证计算速度,主要使用2D有限元模型进行频域、暂态仿真,3D有限元模型仅开展频域分析。
8.一种利用权利要求1所述方法进行电磁暂态仿真的系统,其特征在于:包括进行有限元分析的有限元分析计算机(1)以及对电网络进行建模仿真的电网仿真机(4),电网仿真机(4)对电网络进行仿真计算并将仿真计算结果数据通过通信系统(2)输入到数据处理系统(3)中,数据处理系统(3)对仿真计算结构数据进行处理后通过通信系统(2)输送到有限元计算机(1)内进行基于场分析的有限元计算。
9.根据权利要求8所述的将场分析集成到电网络求解的电磁暂态仿真系统,其特征在于:所述有限元分析计算机(1)和电网仿真机(4)分别连接有显示器装置以及输入设备。
10.根据权利要求8所述的将场分析集成到电网络求解的电磁暂态仿真系统,其特征在于:所述电网仿真机(4)为安装Hypersim软件的计算机。
CN201510740268.3A 2015-11-03 2015-11-03 将场分析集成到电网络求解的电磁暂态仿真方法及系统 Active CN105243232B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510740268.3A CN105243232B (zh) 2015-11-03 2015-11-03 将场分析集成到电网络求解的电磁暂态仿真方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510740268.3A CN105243232B (zh) 2015-11-03 2015-11-03 将场分析集成到电网络求解的电磁暂态仿真方法及系统

Publications (2)

Publication Number Publication Date
CN105243232A true CN105243232A (zh) 2016-01-13
CN105243232B CN105243232B (zh) 2018-05-22

Family

ID=55040880

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510740268.3A Active CN105243232B (zh) 2015-11-03 2015-11-03 将场分析集成到电网络求解的电磁暂态仿真方法及系统

Country Status (1)

Country Link
CN (1) CN105243232B (zh)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105956332A (zh) * 2016-05-31 2016-09-21 河北工业大学 基于专家系统的变压器电磁优化设计方法
CN108595902A (zh) * 2018-07-16 2018-09-28 山东大学 计及谐波影响的直线电机铁耗计算方法及模型
CN111090265A (zh) * 2019-03-28 2020-05-01 惠州市奇思创想科技有限公司 一种基于大数据的工业监控系统及其工作流程
CN111339624A (zh) * 2018-12-19 2020-06-26 国网新疆电力有限公司 基于psasp和emtp/atp短路电流直流分量计算方法
CN112329284A (zh) * 2020-10-08 2021-02-05 华中科技大学 一种同步电机瞬态参数动态时变特性辨识方法
CN114638140A (zh) * 2022-05-19 2022-06-17 国网江西省电力有限公司电力科学研究院 变压器在直流偏磁暂态过程下短期允许运行时长计算方法
CN114970259A (zh) * 2022-05-20 2022-08-30 中国南方电网有限责任公司超高压输电公司广州局 换流变压器建模方法、装置、设备、存储介质和程序产品
CN112329284B (zh) * 2020-10-08 2024-10-22 华中科技大学 一种同步电机瞬态参数动态时变特性辨识方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070282588A1 (en) * 2006-05-31 2007-12-06 Powersmiths International, Inc. Electric power system training simulation software and process for electric power system training
KR20110109553A (ko) * 2010-03-31 2011-10-06 주식회사 하이닉스반도체 전압 인가 장치의 전압 인가 방법
CN104361236A (zh) * 2014-11-17 2015-02-18 国家电网公司 电力设备健康状态的评估方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070282588A1 (en) * 2006-05-31 2007-12-06 Powersmiths International, Inc. Electric power system training simulation software and process for electric power system training
KR20110109553A (ko) * 2010-03-31 2011-10-06 주식회사 하이닉스반도체 전압 인가 장치의 전압 인가 방법
CN104361236A (zh) * 2014-11-17 2015-02-18 国家电网公司 电力设备健康状态的评估方法

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105956332A (zh) * 2016-05-31 2016-09-21 河北工业大学 基于专家系统的变压器电磁优化设计方法
CN105956332B (zh) * 2016-05-31 2019-12-10 河北工业大学 基于专家系统的变压器电磁优化设计方法
CN108595902A (zh) * 2018-07-16 2018-09-28 山东大学 计及谐波影响的直线电机铁耗计算方法及模型
CN111339624A (zh) * 2018-12-19 2020-06-26 国网新疆电力有限公司 基于psasp和emtp/atp短路电流直流分量计算方法
CN111339624B (zh) * 2018-12-19 2023-10-20 国网新疆电力有限公司 基于psasp和emtp/atp短路电流直流分量计算方法
CN111090265A (zh) * 2019-03-28 2020-05-01 惠州市奇思创想科技有限公司 一种基于大数据的工业监控系统及其工作流程
CN112329284A (zh) * 2020-10-08 2021-02-05 华中科技大学 一种同步电机瞬态参数动态时变特性辨识方法
CN112329284B (zh) * 2020-10-08 2024-10-22 华中科技大学 一种同步电机瞬态参数动态时变特性辨识方法
CN114638140A (zh) * 2022-05-19 2022-06-17 国网江西省电力有限公司电力科学研究院 变压器在直流偏磁暂态过程下短期允许运行时长计算方法
CN114970259A (zh) * 2022-05-20 2022-08-30 中国南方电网有限责任公司超高压输电公司广州局 换流变压器建模方法、装置、设备、存储介质和程序产品
CN114970259B (zh) * 2022-05-20 2024-05-31 中国南方电网有限责任公司超高压输电公司广州局 换流变压器建模方法、装置、设备、存储介质和程序产品

Also Published As

Publication number Publication date
CN105243232B (zh) 2018-05-22

Similar Documents

Publication Publication Date Title
CN105243232A (zh) 将场分析集成到电网络求解的电磁暂态仿真方法及系统
Bialek et al. Benchmarking and validation of cascading failure analysis tools
CN102332709B (zh) 一种等效高压直流输电直流保护数字仿真方法及系统
CN102707628B (zh) 一种用于电网安全稳定控制的实时仿真试验研究方法
CN103076548B (zh) 一种用表面电导率和泄漏电流预测闪络电压的方法
Xu et al. Robust transient stability-constrained optimal power flow with uncertain dynamic loads
CN103954925B (zh) 一种基于rtds实时仿真的故障录波器动态测试方法
CN103872681A (zh) 一种基于主配网一体化的在线实时合环方法
CN102184209A (zh) 一种基于cim电网模型接口的仿真数据接入方法
CN106354898A (zh) 一种基于总应变能量密度的焊缝疲劳寿命计算方法
CN103532137B (zh) 一种三相四线低压配电网的状态估计方法
Jalili-Marandi et al. A real-time dynamic simulation tool for transmission and distribution power systems
CN105205244B (zh) 基于机电-电磁混合仿真技术的合环操作仿真系统
CN104716646B (zh) 一种基于注入电流的节点耦合度分析方法
CN104901306A (zh) 一种考虑连锁故障的电网运行安全裕度计算方法
CN103699723B (zh) 一种发电厂机组动力系统模型校核方法
CN101867225A (zh) 变电站综合自动化及继电保护系统一体化调试方法
CN102445642A (zh) 一种电力设备金具电晕试验方法
CN103324858A (zh) 配电网三相潮流状态估计方法
CN103678798A (zh) 一种用于含分布式电源配电网的电磁暂态实时仿真方法
CN117639086A (zh) 一种配电台区的拓扑识别功能评价方法及装置
CN115912615A (zh) 一种配电网状态评估方法及系统
CN108233356B (zh) 一种光伏逆变器控制器一致性评估方法及其评估平台
CN103823143A (zh) 一种特高压增容后变电站设备安全性分析的仿真方法
Ravikumar et al. Efficient modeling of hil multi-grid system for scalability & concurrency in cps security testbed

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant