CN109543240B - 一种基于动态区域饱和j-a理论的电流互感器建模方法 - Google Patents
一种基于动态区域饱和j-a理论的电流互感器建模方法 Download PDFInfo
- Publication number
- CN109543240B CN109543240B CN201811277513.1A CN201811277513A CN109543240B CN 109543240 B CN109543240 B CN 109543240B CN 201811277513 A CN201811277513 A CN 201811277513A CN 109543240 B CN109543240 B CN 109543240B
- Authority
- CN
- China
- Prior art keywords
- current
- saturation
- dynamic
- magnetic
- theory
- 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
- 238000000034 method Methods 0.000 title claims abstract description 50
- 230000005415 magnetization Effects 0.000 claims abstract description 22
- 238000004088 simulation Methods 0.000 claims abstract description 20
- 230000005284 excitation Effects 0.000 claims abstract description 13
- 238000004364 calculation method Methods 0.000 claims abstract description 12
- 230000005291 magnetic effect Effects 0.000 claims description 63
- 230000000694 effects Effects 0.000 claims description 20
- 239000011162 core material Substances 0.000 claims description 14
- 230000008569 process Effects 0.000 claims description 14
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical group [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 claims description 12
- 230000003068 static effect Effects 0.000 claims description 12
- 230000008878 coupling Effects 0.000 claims description 11
- 238000010168 coupling process Methods 0.000 claims description 11
- 238000005859 coupling reaction Methods 0.000 claims description 11
- 230000009466 transformation Effects 0.000 claims description 8
- 230000035699 permeability Effects 0.000 claims description 7
- 229920006395 saturated elastomer Polymers 0.000 claims description 7
- 230000005381 magnetic domain Effects 0.000 claims description 6
- 230000002441 reversible effect Effects 0.000 claims description 5
- 238000013016 damping Methods 0.000 claims description 4
- 230000004907 flux Effects 0.000 claims description 3
- 230000001052 transient effect Effects 0.000 description 25
- 238000012360 testing method Methods 0.000 description 24
- 238000004458 analytical method Methods 0.000 description 4
- 230000008859 change Effects 0.000 description 4
- 238000002307 isotope ratio mass spectrometry Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000011156 evaluation Methods 0.000 description 3
- 230000009471 action Effects 0.000 description 2
- 238000005094 computer simulation Methods 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 230000007257 malfunction Effects 0.000 description 2
- 238000004804 winding Methods 0.000 description 2
- 101100499229 Mus musculus Dhrsx gene Proteins 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000011157 data evaluation Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000005562 fading Methods 0.000 description 1
- 230000005294 ferromagnetic effect Effects 0.000 description 1
- 238000003475 lamination Methods 0.000 description 1
- 238000012821 model calculation Methods 0.000 description 1
- 230000009022 nonlinear effect Effects 0.000 description 1
- 230000001681 protective effect Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/36—Circuit design at the analogue level
- G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
-
- 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)
- Computer Hardware Design (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Testing Of Short-Circuits, Discontinuities, Leakage, Or Incorrect Line Connections (AREA)
Abstract
本发明公开了一种基于动态区域饱和J‑A理论的电流互感器建模方法，具体包括如下步骤：步骤1，获取待建模CT的基本参数，步骤2，设置仿真环境参数，步骤3，根据dM/dH求解微分方程得dM/dt；步骤4，进行区域饱和判断，并根据判断结果进行动态量dMf(t)/dH的更新；步骤5，更新励磁电流，并计算当前二次侧电流i2；步骤6，标记当前计算的计及动态损耗后的磁化强度值Mf(i)和励磁电流值im(i)分别为下一次计算的初始值，继续求解二次侧电流，直至完成仿真计算。本发明通过引入区域饱和的概念，改进电流互感器的动态J‑A模型，给出了提出的动态区域饱和J‑A电流互感器模型的非迭代格式的数值解。
Description
技术领域
本发明属于电流互感器的建模分技术领域，涉及一种基于动态区域饱和J-A理论的电流互感器建模方法。
背景技术
保护用电流互感器(P CT)在配合继电保护系统的安全可靠运行中具有重要的作用。然而，由于铁芯材料的非线性影响，P CT很容易暂态饱和，这可能导致二次侧继电保护的故障或误操作。随着技术的发展，对CT暂态特性的研究逐渐从线性分析和物理实验转变为数字仿真。线性等效分析显然不能解释饱和，而广泛的物理实验的手段是非常昂贵的。因此，以瞬时值的形式准确地对CT(电流互感器)进行建模是当务之急。
目前，J-A理论已经被一些专业软件，例如PSCAD/EMTDC、ATP，用于CT模型的模拟。相比于其他磁滞理论，J-A理论更能直观地描述微观能量过程与宏观磁滞现象，是铁磁元件建模分析的重要理论。一般来说，J-A理论模型由于在频域中的不同应用而包括静态和动态模型。静态J-A模型仅考虑了磁滞现象；动态J-A模型增加的动态过程可适用于数百千Hz变换磁场的讨论，但其中的过剩损耗难以计算，当前过剩损耗的线性化计算模式在衰减暂态情况下精度仍不理想。在宽范围的电流水平测试下，利用J-A理论的CT模型在结果准确性上存在着明显的差异，尚不足以准确地显示CT的饱和特性，特别是在暂态模式下对保护效果的分析影响较大。综上所述，现有的基于J-A理论的电流互感器建模方法对暂态饱和的阐述，效果不佳，且不具有普适性。因此，需要进一步研究适用于CT暂态饱和的建模方法，来准确分析配合保护准确动作的相关问题。
发明内容
本发明的目的是提供一种基于动态区域饱和J-A理论的电流互感器建模方法，该方法通过引入区域饱和的概念，改进电流互感器的动态J-A模型，给出了提出的动态区域饱和J-A电流互感器模型的非迭代格式的数值解。本发明所采用的技术方案是，一种基于动态区域饱和J-A理论的电流互感器建模方法，具体包括如下步骤：
步骤1，获取待建模CT的基本参数，包括：变比、二次侧负载、二次侧漏抗以及铁心的内径、外径、截面积；通过离线参数辨识的方式，获得电流互感器J-A理论的五个基本参数，分别为：形状参数a，磁畴壁耦合系数α，可逆运动系数c，损耗系数k及铁心材料的饱和磁化强度Ms；
步骤2，设置仿真环境参数，包括：初始剩磁值M(i-1)和初始励磁电流im(i-1)、仿真步长以及仿真时长；
步骤3，利用已知的初始励磁电流im(i-1)更新当前的磁场强度H(i)，再利用Weiss耦合关系，得到仅考虑磁滞效应的磁化率微分表达式dM/dH，根据dM/dH求解微分方程得dM/dt；
步骤4，进行区域饱和判断，并根据判断结果进行动态量dMf(t)/dH的更新；
步骤5，更新励磁电流，并计算当前二次侧电流i2；
步骤6，标记当前计算的计及动态损耗后的磁化强度值Mf(i)和励磁电流值im(i)分别为下一次计算的初始值：M(i-1)和im(i-1)，然后开始输入新的一次侧电流，继续求解二次侧电流，直至完成仿真计算。
本发明的特点还在于，
步骤3的具体过程如下：
步骤3.1，通过如下公式(1)计算磁场强度H；
其中，l是环路积分的等效磁链长度，单位m；N是CT变比，im是励磁电流；
步骤3.2，利用已知的励磁电流im(i-1)更新当前的磁场强度H(i)，基于如下公式(2)所描述的Weiss耦合关系，通过初始剩磁值M(i-1)和更新后磁场强度H(i)来对有效磁场强度He做进一步的更新，得到当前的有效磁场强度He(i)；
He＝H+α·M (2)；
其中，He是有效磁场强度；α是磁畴壁耦合系数；
步骤3.3，采用如下公式(3)所描述的Langevin微分表达式求得无磁滞效应的磁化率dMan/dH(i)；
其中，Ms是铁心材料的饱和磁化强度；a是形状参数；
步骤3.4，结合静态J-A理论推导出的磁化率微分表达式，如下公式(4)所示，即可得到仅考虑磁滞效应的磁化率微分表达式dM/dH：
其中，c是描述畴壁形变能损的系数；
δ的取值如下公式(5)所示：
其中，sign的定义为：
根据dM/dH求解微分方程得dM/dt。
步骤4的具体过程为：
区域饱和的判断，仅考虑两段式区域饱和的划分，即：重饱和区域和轻饱和区域，考虑动态作用，公式(4)的静态J-A理论推导的磁化率公式要进一步更新，通过如下公式(7)进行修正：
式中，λ是衰减常数，ωn是自然频率，Mf是计及动态损耗后的磁化强度值；
当
满足时，判断为弱磁饱和区，否则为强磁饱和区域。
工作于强磁饱和区时，有dMf/dt＝dM/dt
工作于弱饱和区时，有dMf/dt＝C×dM/dt，C为常数。
步骤5的具体过程如下：
步骤5.1，采用如下公式(8)求解励磁电流的微分增量：
式中，A是铁心的截面积，单位m2；μ0是真空磁导率，μ0＝4π×10-7T·m/A；L2为CT的漏感和负载电感之和，R2为CT的漏阻和负载电阻之和；
步骤5.2，结合励磁电流初始值，再利用如下公式(9)，完成已知一次侧电流i1实时计算二次侧电流i2计算的CT数值模型
im＝i1/N-i2 (9)。
本发明的有益效果是，本发明所建立的电流互感器J-A动态区域饱和模型在动态J-A模型的基础上改进，既考虑铁心磁滞饱和的损耗，又针对暂态饱和的特征考虑了动态区域饱和下磁化过程，并且定步长的时序数值求解可以实现CT暂态饱和模型的在线实时模拟，易于和专业软件进行嵌套或联合仿真。
附图说明
图1是本发明本发明一种基于动态区域饱和J-A理论的电流互感器建模方法中采用的CT大通流动态模拟试验电路图；
图2是本发明本发明一种基于动态区域饱和J-A理论的电流互感器建模方法的实施例中暂态6kA、1/7额定负载的测试结果示意图；
图3是本发明本发明一种基于动态区域饱和J-A理论的电流互感器建模方法的实施例中暂态48kA、额定负载的模型测试结果示意图。
具体实施方式
下面结合附图和具体实施方式对本发明进行详细说明。
本发明一种基于动态区域饱和J-A理论的电流互感器建模方法，具体包括如下过程：
步骤1，获取建模所需要的基本参数。一个电流互感器模型对应三部分实物：铁心、绕组和负载。对系统保护工作而言，CT暂态饱和所工作的频带下，绕组的趋肤、邻近效应对于二次侧变换的小电流而言可以忽略。而建模的核心即为精确地描述铁心内的磁场情况。因此，首先需知道待建模CT的基本参数(变比、二次侧负载、二次侧漏阻抗以及铁心的内径、外径、截面积等几何尺寸信息)，并通过离线参数辨识的方式获得电流互感器J-A理论的5个基本参数(形状参数a，磁畴壁耦合系数α，可逆运动系数c，损耗系数k及铁心材料的饱和磁化强度Ms)；
步骤2，设置仿真环境参数。包括：仿真步长以及仿真时长。由于所建立的CT暂态饱和模型为数值求解方式，还需给定仿真初始情况，包括初始剩磁值M(i-1)和励磁电流初始值im(i-1)；
步骤3，首先通过公式(1)得到磁场强度H，并利用已知的励磁电流im(i-1)更新当前的磁场强度H(i)。同时考虑到公式(2)所描述的Weiss耦合关系，再通过初始剩磁值M(i-1)和更新后磁场强度H(i)来对有效磁场强度He做进一步的更新得到当前的He(i)。进而用公式(3)所描述的Langevin微分表达式求得无磁滞效应的磁化率dMan/dH(i)。结合静态J-A理论推导出的磁化率微分表达式(式(4))，即可完成初步的仅考虑磁滞效应的磁化率微分表达式dM/dH；
其中，l是环路积分的等效磁链长度，单位m；N是CT变比，im是励磁电流。
He＝H+α·M (2)；
式中，He是有效磁场强度；α是磁畴壁耦合系数。
式中，Ms是铁心材料的饱和磁化强度；a是形状参数；
其中，c是描述畴壁形变能损的系数；
δ的取值如下公式(5)所示：
其中，sign的定义为：
可根据dM/dH求解微分方程得dM/dt。
步骤4，区域饱和的判断和动态量dMf(t)/dH的更新。本发明为增强在线实时计算的效率，仅考虑了两段式区域饱和的划分，即：将动态磁化量分为两种状态：重饱和区域和轻饱和区域。考虑动态作用，公式(4)的静态J-A理论推导的磁化率公式需要进一步更新。采用如下公式(7)进行动态修正：
式中，λ是衰减常数，ωn是自然频率，Mf是计及动态损耗后的磁化强度值。新增两个参数后，此类动态J-A模型可以被求解计算出来。
当
满足时，判断为弱磁饱和区，否则为强磁饱和区域。
工作于强磁饱和区时，有dMf/dt＝dM/dt；
工作于弱饱和区时，有dMf/dt＝C×dM/dt，C为常数，一般可取经验值为2。
步骤5，更新励磁情况求解二次测电流i2。首先按照公式(8)求解励磁电流的微分增量。
式中，A是铁心的截面积，单位m2；μ0是真空磁导率，μ0＝4π×10-7T·m/A；L2为CT的漏感和负载电感之和，R2为CT的漏阻和负载电阻之和。
结合励磁电流初始值，再利用公式(9)，完成已知一次侧电流i1实时计算二次侧电流i2计算的CT数值模型。
im＝i1/N-i2 (9)
步骤6，标记当前计算的计及动态损耗后的磁化强度值Mf(i)和励磁电流值im(i)为下一次计算初始值：M(i-1)和im(i-1)。然后开始输入新的一次侧电流，继续求解二次侧电流，直至完成仿真计算。
本发明一种基于动态区域饱和J-A理论的电流互感器建模方法的工作原理为：引入区域饱和的思想对电流互感器动态J-A模型进行改进，着重强调铁心在暂态条件下的过剩损耗问题，给出了所提出的动态区域饱和J-A电流互感器模型无迭代格式的数值求解计算方法，使其能够在线实时仿真电流互感器的二次侧输出结果。有关动态区域的判断依据推到如下：
根据电力系统暂态故障的分析可知，电流互感器铁心内磁场变化的频率和衰减直接受一、二次侧等效故障回路的时间常数、以及短路电流水平决定。系统中高频量噪声的信号能量较弱，加之电流互感器铁心的叠片工艺，使铁心内的涡流作用已经不明显，故将动态J-A磁化的微分式表达简化为：
此时，通解解析式为
其中，C1是关于dM(H)/dt的变量，C2为常数。
式中，λ是衰减常数，ωn是自然频率。
如果不考虑自然衰减和与之对应的自然频率问题，且当前工作在弱饱和区时，平均磁导率较高，可认为
此时等效的作用相当于dMf(t)/dt倍增一个常量。而当铁心在重饱和区时，平均磁导率变化较慢，此时λ与方向参数、磁导率和磁场微增量呈缓慢的动态变化，此时等效的作用相当于dMf(t)/dt相加了一个与磁化变化方向一致常量。
为降低微分方程求解的复杂度，修正的式(10)按倍增常数和相加常数完成。其中，关于dMf/dH属于轻、重饱和的判断，按dMf(t)/dt的瞬时值以及J-A和电流互感器参数确定。
可以先假设磁化过程无磁滞，即M-H关系满足Langevin方程：
即
再根据可逆磁化量Mrev的概念，即Mrev＝c(Man—M)。该定义是由计算磁畴能损忽略三次高阶项简化得来，即计及磁滞的可逆磁化部分与实际偏小，因此对其做磁化率描述后有
若以原有效场He作的Langevin方程的区间缩放，根据方程(14)，易得
与式(15)整理后有如下不等式成立
根据He＝H+αM的关系，因α恒定且小于1，故对上式再进行缩放，最终判断强、弱磁饱和场的判断阈值选择为：当
满足时，判断为弱磁饱和区；反之为强磁饱和。
由此，J-A动态磁化模型中，将以二区域饱和态对动态磁化率进行实时修正，并求取此时的dMf/dH，用于电流互感器铁心励磁电流的计算，还原CT二次侧电流。
实施例1
以图1所示的CT大通流动态模拟试验的现场测试为例，来说明和比对本发明的正确性和优越性。其中试验对象为P级电流互感器，型号10P20，1200/1，额定负载30VA，其铁心尺寸信息为：r内径×r外径×A＝0.345×0.495m×(7.5×2.5)cm2，试验实际负载有1Ω、1/7额定负载、额定负载，进行了暂态通流6～48kA的若干组大电流试验，获得了大量实际暂态饱和的录波数据后，用以验证和测试建模方法的正确性。同时与静态J-A建模方法和动态J-A建模方法进行性能对照，以体现本发明建模方法的优越性。
建模方法测试的步骤如下：
1)记录所测试CT的基本参数，并将离线辨识的J-A模型参数(a、c、k、α、Ms)在静态J-A建模方法、动态J-A建模方法以及本发明中的动态区域饱和J-A建模方法中，设置为同组参数；
2)将录波试验数据统一设置为仿真步长50μs/点，总仿真时间与试验通流时间一致；
3)将一次侧电流作为输入，分别利用三种方法对若干组试验数据进行测试，得到各自所计算的二次侧电流，并与实际录波的二次侧电流进行比对，以验证建模方法的正确性；
4)定量说明建模结果与试验数据间的精度还原效果，用用均方差IRMS和波形相似度Rr两个指标来评价：
式中：num是评价数据点的总个数；Itest是试验数据的二次侧电流；Imodel是模型计算的二次侧电流。
图2和图3分别展示为试验6kA和48kA暂态电流水平下，10P20,1200/1的电流互感器1/7额定负载和额定负载的计算结果。
表1是相同CT在暂态4中通流水平下配合3种负载，分别在静态J-A建模方法、动态J-A建模方法以及本发明中的动态区域饱和J-A建模方法中与二次侧试验数据评价的结果统计。
表1电流互感器数值模型的结果评价
注：加粗下划线部分为对应试验场景下的最佳指标.
根据表1的数据评价，当负载较小时，本发明所建立的CT模型在传变6kA暂态电流(二次侧折合5A)的电流平均方差为0.72A；传变12kA暂态电流(二次侧折合10A)的IRMS为0.97A；传变24kA暂态电流(二次侧折合20A)的IRMS为1.71A；试验最高的暂态48kA(二次侧折合40A)下IRMS为1.76A，比静态模型的误差小30％，相比动态模型的小了5.9％。本发明所建立的模型淡化了高频场才重点刻画动态损耗问题，对弱磁场下的暂态饱和变化过程进行了磁滞能损修正，使之最大程度地符合实际。相比于动态J-A模型，改进的动态区域饱和模型在6～12kA低通流下的轻饱和工作区的提升效果也十分明显，少数不足于动态模型的也都较为接近；在48kA以上的深度饱和情况，本发明模型的精度也保持较高水平，能够反映实际情况。
本发明提出的动态区域饱和模型能够在还原宽范围的暂态通流试验数据测试中，能够取得较好的效果。另该测试里各模型的数值计算是直接对原始数据在采样率50kHz下求解的，数值计算方式能够适应所有数据集中夹杂的试验噪声，所设计的模型计算方式也较稳定。
Claims (4)
1.一种基于动态区域饱和J-A理论的电流互感器建模方法，其特征在于：具体包括如下步骤：
步骤1，获取待建模CT的基本参数，包括：变比、二次侧负载、二次侧漏阻抗以及铁心的内径、外径、截面积；
通过离线参数辨识的方式，获得电流互感器J-A理论的五个基本参数，分别为：形状参数a，磁畴壁耦合系数α，可逆运动系数c，损耗系数k及铁心材料的饱和磁化强度Ms；
步骤2，设置仿真环境参数，包括：初始剩磁值M(i-1)和初始励磁电流im(i-1)、仿真步长以及仿真时长；
步骤3，利用已知的初始励磁电流im(i-1)更新当前的磁场强度H(i)，再利用Weiss耦合关系，得到仅考虑磁滞效应的磁化率微分表达式dM/dH，根据dM/dH求解微分方程得dM/dt；
步骤4，进行区域饱和判断，并根据判断结果进行动态量dMf(t)/dH的更新；
步骤5，更新励磁电流，并计算当前二次侧电流i2；
步骤6，标记当前计算的计及动态损耗后的磁化强度值Mf(i)和励磁电流值im(i)分别为下一次计算的初始值：M(i-1)和im(i-1)，然后开始输入新的一次侧电流，继续求解二次侧电流，直至完成仿真计算。
2.根据权利要求1所述的一种基于动态区域饱和J-A理论的电流互感器建模方法，其特征在于：所述步骤3的具体过程如下：
步骤3.1，通过如下公式(1)计算磁场强度H；
其中，l是环路积分的等效磁链长度，单位m；N是CT变比，im是励磁电流；
步骤3.2，利用已知的励磁电流im(i-1)更新当前的磁场强度H(i)，基于如下公式(2)所描述的Weiss耦合关系，通过初始剩磁值M(i-1)和更新后磁场强度H(i)来对有效磁场强度He做进一步的更新，得到当前的有效磁场强度He(i)；
He＝H+α·M (2)；
其中，He是有效磁场强度；α是磁畴壁耦合系数；
步骤3.3，采用如下公式(3)所描述的Langevin微分表达式求得无磁滞效应的磁化率dMan/dH(i)；
其中，Ms是铁心材料的饱和磁化强度；a是形状参数；
步骤3.4，结合静态J-A理论推导出的磁化率微分表达式，如下公式(4)所示，即可得到仅考虑磁滞效应的磁化率微分表达式dM/dH：
其中，c是描述畴壁形变能损的系数；
δ的取值如下公式(5)所示：
其中，sign的定义为：
根据dM/dH求解微分方程得dM/dt。
3.根据权利要求2所述的一种基于动态区域饱和J-A理论的电流互感器建模方法，其特征在于：所述步骤4中进行区别饱和判断的过程为：
将动态磁化量分为两种状态：重饱和区域和轻饱和区域，考虑动态作用，公式(4)的静态J-A理论推导的磁化率公式要进一步更新，通过如下公式(7)进行修正：
式中，λ是衰减常数，ωn是自然频率，Mf是计及动态损耗的磁化强度值；
当
满足时，判断为弱磁饱和区，否则为强磁饱和区域；
工作于强磁饱和区时，有dMf/dt＝dM/dt；
工作于弱饱和区时，有dMf/dt＝C×dM/dt，C为常数。
4.根据权利要求3所述的一种基于动态区域饱和J-A理论的电流互感器建模方法，其特征在于：所述步骤5的具体过程如下：
步骤5.1，采用如下公式(8)求解励磁电流的微分增量：
式中，A是铁心的截面积，单位m2；μ0是真空磁导率，μ0＝4π×10-7T·m/A；L2为CT的漏感和负载电感之和，R2为CT的漏阻和负载电阻之和；
步骤5.2，结合励磁电流初始值，再利用公式(9)，完成已知一次侧电流i1实时计算二次侧电流i2计算的CT数值模型
im＝i1/N-i2 (9)。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811277513.1A CN109543240B (zh) | 2018-10-30 | 2018-10-30 | 一种基于动态区域饱和j-a理论的电流互感器建模方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811277513.1A CN109543240B (zh) | 2018-10-30 | 2018-10-30 | 一种基于动态区域饱和j-a理论的电流互感器建模方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109543240A CN109543240A (zh) | 2019-03-29 |
| CN109543240B true CN109543240B (zh) | 2022-09-16 |
Family
ID=65845491
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811277513.1A Active CN109543240B (zh) | 2018-10-30 | 2018-10-30 | 一种基于动态区域饱和j-a理论的电流互感器建模方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109543240B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113792417B (zh) * | 2021-08-23 | 2023-10-20 | 鞍钢集团矿业有限公司 | 一种基于给矿参数与智能算法的强磁机电流优化方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0882990A2 (de) * | 1997-06-04 | 1998-12-09 | Siemens Aktiengesellschaft | Verfahren und Vorrichtung zur Erkennung und Korrektur eines gesättigten Stromverlaufs eines Stromwandlers |
| GB0409767D0 (en) * | 2004-04-30 | 2004-06-09 | Dunlop Aerospace Ltd | Method of calculating a hysteresis charcteristics and hysteresis charcteristic modelling apparatus |
| CN107015187A (zh) * | 2017-06-14 | 2017-08-04 | 国家电网公司 | 一种保护用电流互感器j‑a模型参数检验方法 |
| CN107103163A (zh) * | 2017-05-26 | 2017-08-29 | 华中科技大学 | 一种基于物理实验确定电流互感器j‑a仿真模型的方法 |
-
2018
- 2018-10-30 CN CN201811277513.1A patent/CN109543240B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0882990A2 (de) * | 1997-06-04 | 1998-12-09 | Siemens Aktiengesellschaft | Verfahren und Vorrichtung zur Erkennung und Korrektur eines gesättigten Stromverlaufs eines Stromwandlers |
| GB0409767D0 (en) * | 2004-04-30 | 2004-06-09 | Dunlop Aerospace Ltd | Method of calculating a hysteresis charcteristics and hysteresis charcteristic modelling apparatus |
| CN107103163A (zh) * | 2017-05-26 | 2017-08-29 | 华中科技大学 | 一种基于物理实验确定电流互感器j‑a仿真模型的方法 |
| CN107015187A (zh) * | 2017-06-14 | 2017-08-04 | 国家电网公司 | 一种保护用电流互感器j‑a模型参数检验方法 |
Non-Patent Citations (2)
| Title |
|---|
| 基于J-A动态磁滞模型的电流互感器谐波变换建模及实验验证;陈利翔等;《电气技术》;20160815(第08期);全文 * |
| 电流互感器J-A模型参数辨识及大通流动模试验;雷阳等;《中国电机工程学报》;20160930;全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109543240A (zh) | 2019-03-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105740569B (zh) | 基于暂态大电流测试技术的电流互感器工程模型建立方法 | |
| CN112364490B (zh) | 一种铁磁材料磁滞特性的模型参数辨识方法 | |
| CN106257462A (zh) | 一种变压器内部短路故障非线性仿真的方法 | |
| CN109444776A (zh) | 三相三柱变压器铁芯剩磁测算方法、系统及存储介质 | |
| Rezaei-Zare et al. | An accurate current transformer model based on Preisach theory for the analysis of electromagnetic transients | |
| CN109884402B (zh) | 一种三维非对称结构高温超导磁体交流损耗的获取方法 | |
| CN110399677A (zh) | 基于偏磁状态下改进j-a公式的变压器直流偏磁仿真模拟方法 | |
| CN103149478B (zh) | 一种直流偏磁条件下电流互感器工作点检测方法 | |
| Saeed et al. | Improved inductance calculation in variable power inductors by adjustment of the reluctance model through magnetic path analysis | |
| CN101975892B (zh) | 一种换流阀用饱和电抗器的损耗计算方法 | |
| CN109460609B (zh) | 一种变压器绕组匝间短路反向环流计算方法 | |
| CN102708295B (zh) | 一种电工钢片偏磁特性的分析方法 | |
| CN109543240B (zh) | 一种基于动态区域饱和j-a理论的电流互感器建模方法 | |
| CN109145447A (zh) | 一种可控磁饱和电抗器建模方法 | |
| Tian et al. | Neural network model for magnetization characteristics of ferromagnetic materials | |
| CN110289595B (zh) | 一种高阻抗变压器的零序过流保护整定方法及相关装置 | |
| CN103258124A (zh) | 一种高压变电站电磁设备工频磁场评估方法 | |
| Mork et al. | Parameter estimation methods for five-limb magnetic core model | |
| CN113255285B (zh) | 基于ja模型与对偶原理的单相四柱式变压器的建模方法 | |
| Zhao et al. | Study on dynamic hysteretic and loss properties of silicon steel sheet under hybrid harmonic and DC bias excitation | |
| CN116822287A (zh) | 换流变压器匝间短路多耦合场模型降阶方法及系统 | |
| Lin et al. | Losses calculation in transformer tie plate using the finite element method | |
| CN115238554A (zh) | 一种基于多物理场双向耦合精细化建模的分析方法 | |
| CN113239649B (zh) | 一种变压器的建模方法 | |
| Ma et al. | Modeling and Multi-Objective Optimization Design of High-Speed on/off Valve 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 | ||
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20240723 Address after: 052300 No. 0369, fangbei East Street, Xinji City, Hebei Province Patentee after: HEBEI SHENKE ELECTRIC POWER CO.,LTD. Country or region after: China Address before: 710048 No. 5 Jinhua South Road, Shaanxi, Xi'an Patentee before: XI'AN University OF TECHNOLOGY Country or region before: China |