CN107103163B - A method of current transformer J-A simulation model is determined based on Physical Experiment - Google Patents
A method of current transformer J-A simulation model is determined based on Physical Experiment Download PDFInfo
- Publication number
- CN107103163B CN107103163B CN201710384368.6A CN201710384368A CN107103163B CN 107103163 B CN107103163 B CN 107103163B CN 201710384368 A CN201710384368 A CN 201710384368A CN 107103163 B CN107103163 B CN 107103163B
- Authority
- CN
- China
- Prior art keywords
- current transformer
- simulation model
- physical experiment
- model
- current
- 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
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
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)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The present invention discloses a kind of method for determining current transformer J-A simulation model based on Physical Experiment, comprising: passes through the coiling N on the iron core of current transformer1Circle coil constructs the Physical Experiment model of current transformer；By the primary side on-load voltage in Physical Experiment model, determine that hysteresis loop to be fitted, the hysteresis loop are the magnetic flux density curve corresponding with the variation relation of magnetic field strength on iron core；9 initial parameters of predetermined current mutual inductor J-A simulation model determine the end value of 9 parameters of current transformer J-A simulation model according to 9 initial parameters and hysteresis loop iterative fitting to be fitted；Current transformer J-A simulation model, the measurement process for the current transformer that J-A simulation model is used for simulating grid measurement are determined according to the end value of 9 parameters of current transformer J-A simulation model.The present invention can construct the different types of practical current transformer Physical Experiment model of industry using the lower Physical Experiment electric current of the order of magnitude, have very big promotional value.
Description
Technical field
The present invention relates to Relay Protection Technology in Power System fields, are determined more particularly, to a kind of based on Physical Experiment
The method of current transformer J-A simulation model.
Background technique
In the power system, the accurate primary current for flowing through equipment can not measure, and need current transformer according to electricity
The biggish primary current of numerical value is converted to the lesser secondary current of numerical value by certain no-load voltage ratio by magnetic induction principle, is then measured
Reduction obtains primary current.In fact, because of this kind of measurement method and current transformer measurement performance difference, it is mutual by electric current
There are certain errors with practical primary current for the primary current that the progress of disease reduction of sensor obtains.Meanwhile in recent years, in electric power
During system transient modelling, power grid has occurred a lot of since CT saturation causes transmission line of electricity, transformer, generator differential
The accident of false protection, has seriously affected power supply reliability and power network safety operation.Study transient state process of electric power system pair
CT saturation, so that the effective means having an impact to protection is Digital Simulation.
The simulation model that corresponding current transformer is designed for the current transformer that power grid measurement uses, can simulate electricity
Primary current is passed by analogue simulation close to actual current mutual inductor to the actual progress of disease process of secondary current in current transformer
The transient state and steady-state performance of change, more accurately to measure the primary current for flowing through equipment.Based on scholar Jiles and Atherton
It is proposed about ferromagnetic material phenome-nology theory on the basis of, using Nonlinear System of Equations come in approximate various situations iron core excitation spy
Property, propose a kind of current transformer J-A simulation model, transient characterisitics quilt of the model because can preferably reflect current transformer
It is widely recognized as.
But the order of magnitude of the corresponding electric current to be tested of current transformer used in power grid measurement is higher, it is existing such as 2500A
When testing the BH curve of practical current transformer, when primary side loading current, is difficult up to the so high order of magnitude, and electric current
Mutual inductor J-A simulation model needs to establish based on on-the-spot test BH curve data, therefore due to practical amount of current to be measured
Grade is larger to cause the existing method for establishing current transformer J-A simulation model to lack promotional value.And existing current transformer J-A
The case where J-A model parameter in simulation model is typical default, can only represent certain single kind of current transformer, does not have general
Property.And due to there are ferromagnetic material properties difference and current transformer type (TPY, P, PR type etc.) difference, it should in simulation software
The typical default parameters of model is not fully suitable for China and produces and be used for the current transformer that power grid measurement uses, need according to
Factually border current transformer redefines corresponding current transformer J-A simulation model.
Summary of the invention
In view of the drawbacks of the prior art, it is an object of the invention to solve the prior art due to practical electric current number to be measured
Magnitude is larger to cause the existing method for establishing current transformer J-A simulation model to lack promotional value, and existing current transformer
The case where J-A model parameter in J-A simulation model is typical default, can only represent certain single kind of current transformer, does not have
The technical issues of versatility.
To achieve the above object, the present invention provides one kind determines current transformer J-A simulation model based on Physical Experiment
Method, comprising: pass through the coiling N on the iron core of current transformer1Circle coil constructs the Physical Experiment mould of current transformer
Type, N1The rated primary current of the current transformer used according to power grid measurement and the specified primary electricity of the Physical Experiment model
Stream determines；By the primary side on-load voltage in the Physical Experiment model, hysteresis loop to be fitted is determined, the magnetic hysteresis is returned
Line is the magnetic flux density curve corresponding with the variation relation of magnetic field strength on the iron core；Predetermined current mutual inductor J-A emulation
9 parameters and its initial value of model, according to the initial value of 9 parameters and the hysteresis loop iterative fitting to be fitted
Determine that the end value of 9 parameters of the current transformer J-A simulation model, 9 parameters are to describe the iron core magnetic hysteresis
The parameter of loop line；The current transformer J- is determined according to the end value of 9 parameters of the current transformer J-A simulation model
A simulation model, the current transformer J-A simulation model are used to simulate the measurement for the current transformer that the power grid measurement uses
Process.
It should be noted that " Physical Experiment " that the embodiment of the present invention refers to refers to building electricity under the conditions of conventional finite
Current transformer model, such as the available voltage range of Physical Experiment condition are 0~800V, but the voltage in practical power systems
Up to 500kV.Available current range is 0~10A under the conditions of Physical Experiment, and the electric current in practical power systems is reachable
2500A, i.e., the electric current in practical power systems is up to kA grades.Wherein, the current transformer that power grid measurement uses is used for reality
The measurement of electric system.
The limited conditions that the embodiment of the present invention is provided based on Physical Experiment, by iron core winding N1Circle coil guarantees physics
The added electric current of experiment is constant in the magnetomotive force of iron core compared to the actual current of power grid measurement in the magnetomotive force of iron core, so that being based on
Physical Experiment condition is established J-A simulation model and is possibly realized.In addition, the embodiment of the present invention can be with by fitting J-A simulation model
The feature of various aspects description iron core hysteresis loop, so that the reliability of the J-A simulation model based on Physical Experiment foundation mentions significantly
It rises.
Optionally, N1It is determined by following formula:
Wherein, Ipn_ctFor the rated primary current for the current transformer that power grid measurement uses, In_testFor the Physical Experiment
The rated primary current of model.
Specifically, pass through Ipn_ctAnd In_testRatio setting, ensure that the current transformer model that Physical Experiment is established
Magnetomotive force is constant compared to magnetomotive force under the conditions of power grid actual measurement.
Optionally, it in the alternating voltage of the primary side load different voltages virtual value of the Physical Experiment model, determines not
The variation relation of corresponding magnetic flux density and magnetic field strength with voltage effective value, the alternating voltage corresponding one of every kind of voltage effective value
The variation relation of kind magnetic flux density and magnetic field strength；In the corresponding magnetic flux density of the different voltages virtual value and magnetic field strength
A kind of variation relation of saturation is chosen in variation relation as the hysteresis loop to be fitted.
Optionally, the mistake of the Physical Experiment model primary side is loaded onto the corresponding alternating voltage of each voltage effective value
Cheng Zhong acquires the voltage u on winding described in multiple groupspWith the voltage u on sampling resistorR, the sampling resistor and the winding string
Connection；Determine every group of upAnd uRCorresponding magnetic flux density B and magnetic field strength H, obtains especially by following formula:
Wherein, S and l is respectively the core section product and the average length of magnetic path of Physical Experiment model；im1For primary side progress
Leveling treated electric current,Φ2For the leveling that passes through in the iron core it is steady after
Magnetic flux；Φ2=Φ1-(max(Φ1)+min(Φ1))/2, Φ1To remove the magnetic after the biasing magnetic flux passed through in the iron core
It is logical, Φ1=Φ-Φpz, Φ is that integration method measures the initial magnetic flux passed through in the iron core,ΦpzFor biasing
Magnetic flux,T is the AC voltage cycle of primary side load；According to multiple groups u under each voltage effective valuepAnd uR
Corresponding multiple groups B and H value, determine the variation relation of each voltage effective value corresponding magnetic flux density and magnetic field strength.
It should be noted that B curve corresponding with the variation relation of H can be described as hysteresis loop.
Optionally, 9 parameters and its initial value of predetermined current mutual inductor J-A simulation model, according at the beginning of 9 parameters
Initial value and the hysteresis loop iterative fitting to be fitted determine 9 parameter values of the current transformer J-A simulation model
End value, comprising:
(a) control parameter of blending heredity simulated annealing fitting algorithm: initial temperature T is determineds, exit temperature Tf, cooling system
Number k and maximum cooling number Mmax；
(b) initial value for determining 9 parameters of the current transformer J-A simulation model, according to the initial of 9 parameters
Value generates initial population G1, and cooling number is enabled to add up M=0, run algebra N=0；
(c) resetting control temperature T=Ts；
(d) to initial population G1New filial generation G is generated by the genetic manipulation of genetic algorithm2, the new filial generation G2Corresponding current
9 new parameter values of mutual inductor J-A simulation model, determine G2The variation relation of corresponding B and H and the B and H to be fitted
Variation relation between evaluation function J (G), more excellent solution is selected according to Metropolis criterion based on J (G), with P (G2) it is general
Rate receives new explanation G2Replace G1, the P (G2) according to J (G1) and J (G2) determine；
(e) algorithm enters the next generation, N=N+1, and control temperature is reduced to T=T*k, and k is coefficient of temperature drop；
(f) judge to control whether temperature reaches exit temperature Tf, if so, (g) is thened follow the steps, it is no to then follow the steps (d)；
(g) amendment cooling number, M=M+1；
(h) judge whether cooling number reaches limit value Mmax, if so, thening follow the steps (i), otherwise；It executes step (c)；
(i) initial population G is exported1The corresponding current transformer J-A simulation model of newest filial generation 9 parameters, as
The end value of 9 parameters of the current transformer J-A simulation model.
Optionally, the evaluation function J (G) describe new filial generation corresponding B and H variation relation and magnetic hysteresis to be fitted
The similarity of loop line, the smaller then degree of fitting of J (G) is higher, the corresponding current transformer J-A emulation of the higher new filial generation of explanation of degree of fitting
9 parameters of model are more accurate, and J (G) is determined especially by following publicity:
Wherein, HmeasuredFor the magnetic field strength of Physical Experiment model measurement, HsimulatedIt is strong for the corresponding magnetic field of new filial generation
Degree, n are the sampling u in Physical Experiment modelpAnd uRPoints.
Optionally, the Metropolis criterion describes the new filial generation G for receiving genetic algorithm generation2Probability, pass through
Following formula determines:
If Metropolis criterion is new explanation G2Better than G1, i.e. J (G2) < J (G1), then centainly receive G2；Otherwise only with P
(G2) probability receive G2。
Optionally, 9 parameters of the current transformer J-A simulation model include: saturation flux Ms, region pinning coefficient
K, region can disturb coefficient of coup α between coefficient c, domain, curve shape coefficient a1、a2、a3, b and correction factor β.
Specifically, J-A theoretical description is relationship between magnetic flux density B (or magnetization M) and magnetic field strength H.It should
Theory saturation flux Ms, region pinning coefficient k, region can disturb coefficient of coup α between coefficient c, domain, curve shape coefficient a1、a2、
a3, 9 parameters such as b and correction factor β curvilinear characteristic is described.It can many-sided description iron core magnetic by above-mentioned 9 parameters
The feature of hysteresis curves, so that the reliability of the J-A simulation model based on Physical Experiment foundation greatly promotes.Optionally, the base
It is used to determine the different types of electric current that power grid measurement uses in the method that Physical Experiment determines current transformer J-A simulation model
The J-A simulation model of mutual inductor；The sectional area S for the different types of current transformer core that the power grid measurement uses and average
Length of magnetic path l is different.
In general, through the invention it is contemplated above technical scheme is compared with the prior art, have below beneficial to effect
Fruit:
(1) it is mutual to can use the different types of practical electric current of industry of the lower Physical Experiment electric current building of the order of magnitude by the present invention
Sensor Physical Experiment model has very big promotional value.
(2) BH curve that the present invention can be measured based on the Physical Experiment model of the different practical types of industry, simulates
9 parameters of J-A simulation model,
It may be directly applied in present digital simulation software, described by 9 parameter various aspects of J-A simulation model
The feature of iron core hysteresis loop, so that the reliability of the J-A simulation model based on Physical Experiment foundation greatly promotes.
Detailed description of the invention
Fig. 1 provides the method flow that current transformer J-A simulation model is determined based on Physical Experiment for the embodiment of the present invention
Schematic diagram；
Fig. 2 is the current transformer J-A simulation parameters acquisition side provided in an embodiment of the present invention based on Physical Experiment
Method flow diagram；
Fig. 3 is Physical Experiment precircuit schematic diagram provided in an embodiment of the present invention；
Fig. 4 is blending heredity simulated annealing fitting algorithm flow chart provided in an embodiment of the present invention；
Fig. 5 is that P grades of current transformer hysteresis loop test result provided in an embodiment of the present invention and J-A simulation model are run
Result schematic diagram；
Fig. 6 is that PR grades of current transformer hysteresis loop test result provided in an embodiment of the present invention and J-A simulation model are transported
Row result schematic diagram；
Fig. 7 is that TPY grades of current transformer hysteresis loop test result provided in an embodiment of the present invention and J-A simulation model are transported
Row result schematic diagram.
Specific embodiment
In order to make the objectives, technical solutions, and advantages of the present invention clearer, with reference to the accompanying drawings and embodiments, right
The present invention is further elaborated.It should be appreciated that the specific embodiments described herein are merely illustrative of the present invention, and
It is not used in the restriction present invention.
Fig. 1 provides the current transformer J-A model parameter based on Physical Experiment for the embodiment of the present invention and determines method flow
Schematic diagram；Including step 102- step 104.
In step 101, pass through the coiling N on the iron core of current transformer1Circle coil constructs the physics of current transformer
Experimental model, N1The rated primary current of the current transformer used according to power grid measurement and the Physical Experiment model it is specified
Primary current determines.Optionally, N1It is determined by following formula:
Wherein, Ipn_ctFor the rated primary current for the current transformer that power grid measurement uses, In_testFor the Physical Experiment
The rated primary current of model.
Specifically, the order of magnitude for the primary side voltage that Physical Experiment adds is lower, and corresponding bring primary side current is relatively low,
The order of magnitude for the rated current that corresponding Physical Experiment adds is relatively low.When the variation of the specifications parameter of iron core, correspondence is different types of
The practical current transformer physical model of industry.That is, the embodiment of the present invention can use the lower Physical Experiment of the order of magnitude
The iron core of electric current and different size constructs the different types of practical current transformer Physical Experiment model of industry, has very big push away
Wide value.
In step 102, by the primary side on-load voltage in the Physical Experiment model, determine that magnetic hysteresis to be fitted is returned
Line, the hysteresis loop are the magnetic flux density curve corresponding with the variation relation of magnetic field strength on iron core.
In step 103,9 parameters and its initial value of predetermined current mutual inductor J-A simulation model are joined according to described 9
Several initial values and the hysteresis loop iterative fitting to be fitted determine 9 ginsengs of the current transformer J-A simulation model
Number；9 parameters are to describe the parameter of the iron core hysteresis loop.
In step S104, the current transformer J- is determined according to 9 parameters of the current transformer J-A simulation model
A simulation model, the current transformer J-A simulation model are used to simulate the measurement for the current transformer that the power grid measurement uses
Process.
The BH curve that the embodiment of the present invention can be measured based on the Physical Experiment model of the different practical types of industry, simulation
9 parameters of J-A simulation model out, may be directly applied in present digital simulation software.
Aiming at the problem that prior art cannot accurately obtain different types of industry practical current transformer J-A model, this
Invention provides a kind of current transformer J-A simulation parameters acquisition methods based on physical test and parameter fitting algorithm,
It is intended to obtain the practical J-A simulation parameters for not having to type current mutual inductor.The parameter of acquisition is used directly for mathematics point
Analysis or magnetic Transient State Simulation Software (Power Systems Computer Aided Design, PSCAD) Digital Simulation are structure
The simulation model of on-site current transformer is built, analysis field accident plays the role of great.
According to the present invention, a kind of current transformer J-A simulation model based on physical test and parameter fitting algorithm is provided
Parameter acquiring method, this approach includes the following steps, and flow chart is as shown in Figure 2.
(1) building is suitable for the current transformer physical model of physical test.
(2) basic parameters such as the core area S and length of magnetic path l of current transformer physical model are tested.
(3) test circuit built is taken, test obtains the hysteresis loop of the current transformer for parameter fitting.
(4) it is based on blending heredity simulated annealing fitting algorithm, the J-A that fitting obtains description current transformer ferromagnetic characteristic is imitative
True mode parameter.
The current transformer physical model that further building is suitable for physical test in (1) mainly includes following content:
(11) first winding the number of turns is calculated.The rated primary current I of practical current transformerpn_ct(kilo-ampere), dynamic
Model experiment room cannot generate kA grades of primary current, in order to guarantee that magnetomotive force is constant, need to increase according to magnetomotive force equivalence principle
The coil turn of first winding.If the rated current I that physical test addsn_test, then need in current transformer primary side coiling
The number of turns is
(12) according to the number of turns is calculated, the coil of the number of turns is wound on current transformer core by the way of using close,
Construct the physical model of the current transformer for Physical Experiment.
(3) hysteresis loop of measurement current transform er includes the following steps.
(31) test circuit built is taken, the first siding ring both ends that coiling is completed are connected in parallel to exchange adjustable voltage, and seal in
Sample resistance is to measure electric current.
(32) as shown in figure 3, first winding side corresponds to primary current, secondary winding side corresponds to secondary current.Pass through adjusting
The alternating voltage of the different voltages virtual value of primary side input, the voltage on winding and the sampling next time of record different voltages size
Ohmically voltage, respectively upAnd uRIf upCorresponding channel is " channel 1 ", uRCorresponding channel is " channel 2 ".
(33) exciting current original value is calculated, as shown in formula (1)；Current leveling is carried out, as shown in formula (2)；Using smooth
Function carries out deburring processing to electric current.
Iron core magnetic flux original value is calculated using off-line voltage integration method, as shown in formula (3)；Magnetic flux leveling is then carried out, such as
Shown in formula (4)~(6), ΦpzFor magnetic flux bias, magnetic flux bias is the virtual value in the point future a cycle, and the moment becomes
Change.
im1=im-(max(im)+min(im))/2 (2)
Φ1=Φ-Φpz (5)
Φ2=Φ1-(max(Φ1)+min(Φ1))/2 (6)
(34) magnetic flux density B and magnetic field strength H is calculated by formula (7) and formula (8).
(35) hysteresis loop of current transformer, i.e. measurement BH curve can be obtained according to the mapping of above-mentioned data.
(36) each hysteresis loop under different voltages size is compared, it is quasi- as parameter to select the curve that one has been saturated
The curve of conjunction.
(4) blending heredity simulated annealing fitting algorithm (genetic algorithm, simulated annealing) is used, it is mutual that fitting obtains electric current
The flow chart of sensor J-A simulation parameters is as shown in figure 4, specifically comprise the following steps.
(a) control parameter of B-H array and algorithm: initial temperature T is inputteds, exit temperature Tf, coefficient of temperature drop k and maximum
Cool down number Mmax, one group of parameter (T can be provided here according to required precision flexible setting parameters=10, Tf=1, k=
0.96、Mmax=100)；
(b) it initializes.The initial value (set initial value reference software to default value) for setting 9 parameters generates full
The initial population G of sufficient condition1；The relevant parameter of algorithm: and cooling number is enabled to add up M=0, run algebra N=0；
(c) resetting control temperature T=Ts；
(d)G1New filial generation G is generated by the genetic manipulation of genetic algorithm2, to G2It is mapped to obtain BH curve, calculates it
Evaluation function value J (G) shown in one of evaluation function such as formula (9), selects more excellent solution according to Metropolis criterion, with P
(G2) probability receive new explanation G2Replace G1；
(e) algorithm enters the next generation, N=N+1, and control temperature is reduced to T=T*k；
(f) judge whether to reach exit temperature Tf, if so, turning (g), otherwise turn (d)；
(g) amendment cooling number, M=M+1；
(h) judge whether cooling number reaches limit value Mmax, if so, turning (i), otherwise turn (c)；
(i) optimum results are exported.
Evaluation function J (G) describes the superiority and inferiority degree for being fitted to obtain parameter, i.e., emulation gained BH curve and original B-H are bent
The fitting similarity of line, expression formula are provided by formula (9), and the smaller then degree of fitting of J (G) is higher.
Metropolis criterion describes the new filial generation G for receiving genetic algorithm generation2Probability, provided by formula (10).
If J (G2) < J (G1), indicate new explanation G2Better than G1, then centainly receive G2；Otherwise G is only received with certain probability2.?
In one optional example, the Current Mutual Inductance provided by the invention based on physical test and blending heredity simulated annealing fitting algorithm
Device J-A simulation parameters acquisition methods include the following steps:
(1) building is suitable for practical type (P, PR and TPY type) the current transformer physical model of difference of physical test.
(11) first winding the number of turns is calculated.The rated primary current of practical current transformer is 2500A, the specified electricity in laboratory
Stream is 10A or so, need to increase the coil turn of first winding according to magnetomotive force equivalence principle.It then needs primary in current transformer
Side coiling the number of turns is
(12) according to calculated the number of turns, the coil of the number of turns is wound on current transformer core, that is, constructs and is used for
The physical model of the current transformer of Physical Experiment.
(2) the core area S and length of magnetic path l of current transformer physical model are tested.
Specifically, the current transformer J-A model parameter provided in an embodiment of the present invention based on Physical Experiment determines method
The current transformer physical model of the practical type of difference suitable for physical test.Wherein, the Current Mutual Inductance of different practical types
The difference of device physical model is that their core section product S is different with average length of magnetic path l.Wherein, different practical types
Current transformer model corresponds to the different types of current transformer that power grid measurement uses.Therefore the electric current that power grid measurement uses is mutual
Sensor also is understood as practical current transformer.Can according to S and l by the current transformer of Bu Tong practical type be divided into P, PR and
TPY type current transformer, as shown in table 1.
1 current transformer basic parameter of table
(3) test obtains the hysteresis loop of P, PR and TPY type current transformer, and selects to obtain the deeper song of degree of saturation
Curve of the line as parameter fitting, as shown in solid in Fig. 5~Fig. 7.
(4) blending heredity simulated annealing fitting algorithm, the J-A model ginseng of fitting description current transformer ferromagnetic characteristic are based on
Number inputs simulation model, the hysteresis loop emulated is as shown in the dotted line in Fig. 5~Fig. 7 as shown in table 2.
2 P grades of table, PR grades, TPY grades of current transformer J-A model parameters
Fig. 5-Fig. 7 is respectively P, PR and TPY grade of current transformer hysteresis loop test result provided in an embodiment of the present invention
With J-A simulation model operation result schematic diagram；Comparison diagram 5 to Fig. 7 experimental result as it can be seen that providing according to embodiments of the present invention
The goodness of fit of the hysteresis loop of hysteresis loop and actual measurement that J-A simulation model emulates is higher, illustrates implementation of the present invention
What example provided is imitated based on the application of the J-A simulation parameters discrimination method of physical test and blending heredity simulated annealing fitting algorithm
Fruit is good, can effectively be fitted the J-A simulation parameters of different practical types, specific good application prospect.
As it will be easily appreciated by one skilled in the art that the foregoing is merely illustrative of the preferred embodiments of the present invention, not to
The limitation present invention, any modifications, equivalent substitutions and improvements made within the spirit and principles of the present invention should all include
Within protection scope of the present invention.
Claims (7)
1. a kind of method for determining current transformer J-A simulation model based on Physical Experiment characterized by comprising
Pass through the coiling N on the iron core of current transformer1Circle coil constructs the Physical Experiment model of current transformer, N1According to
The rated primary current of the rated primary current for the current transformer that power grid measurement uses and the Physical Experiment model determines；
By the primary side on-load voltage in the Physical Experiment model, hysteresis loop to be fitted is determined, it is described to be fitted
Hysteresis loop is the magnetic flux density curve corresponding with the variation relation of magnetic field strength on the iron core；
9 parameters and its initial value of predetermined current mutual inductor J-A simulation model, according to the initial value of 9 parameters and institute
The end value that hysteresis loop iterative fitting to be fitted determines 9 parameters of the current transformer J-A simulation model is stated, it is described
9 parameters are the parameter of the description hysteresis loop to be fitted；9 parameter packets of the current transformer J-A simulation model
It includes: saturation flux Ms, region pinning coefficient k, region can disturb coefficient of coup α between coefficient c, domain, curve shape coefficient a1、a2、a3、b
And correction factor β；
9 parameters and its initial value of predetermined current mutual inductor J-A simulation model, according to 9 initial parameter values and described
Hysteresis loop iterative fitting to be fitted determines the end value of 9 parameter values of the current transformer J-A simulation model, packet
It includes:
(a) control parameter of blending heredity simulated annealing fitting algorithm: initial temperature T is determineds, exit temperature Tf, coefficient of temperature drop k
With maximum cooling number Mmax；
(b) initial value for determining 9 parameters of the current transformer J-A simulation model is produced according to the initial value of 9 parameters
Raw initial population G1, and cooling number is enabled to add up M=0, run algebra N=0；
(c) resetting control temperature T=Ts；
(d) to initial population G1New filial generation G is generated by the genetic manipulation of genetic algorithm2, the new filial generation G2Corresponding current mutual inductance
9 parameter values of device J-A simulation model, determine G2Between the variation relation of corresponding B and H and the hysteresis loop to be fitted
Evaluation function J (G), more excellent solution is selected according to Metropolis criterion based on J (G), with P (G2) probability receive new explanation G2It takes
For G1, the P (G2) according to J (G1) and J (G2) determine；
(e) genetic algorithm enters the next generation, N=N+1, and control temperature is reduced to T=T*k, and k is coefficient of temperature drop；
(f) judge to control whether temperature reaches exit temperature Tf, if so, (g) is thened follow the steps, it is no to then follow the steps (d)；
(g) amendment cooling number, M=M+1；
(h) judge whether cooling number reaches limit value Mmax, if so, thening follow the steps (i), otherwise；It executes step (c)；
(i) initial population G is exported1The corresponding current transformer J-A simulation model of newest filial generation 9 parameters, as the electricity
The end value of 9 parameters of current transformer J-A simulation model；
The current transformer J-A emulation mould is determined according to the end value of 9 parameters of the current transformer J-A simulation model
Type, the current transformer J-A simulation model are used to simulate the measurement process for the current transformer that the power grid measurement uses.
2. the method according to claim 1 for determining current transformer J-A simulation model based on Physical Experiment, feature exist
In N1It is determined by following formula:
Wherein, Ipn_ctFor the rated primary current for the current transformer that power grid measurement uses, In_testFor the Physical Experiment model
Rated primary current.
3. the method according to claim 1 for determining current transformer J-A simulation model based on Physical Experiment, feature exist
In, by the primary side on-load voltage in the Physical Experiment model, determining hysteresis loop to be fitted, comprising:
In the alternating voltage of the primary side load different voltages virtual value of the Physical Experiment model, different voltages virtual value is determined
The variation relation of corresponding magnetic flux density and magnetic field strength, the alternating voltage of every kind of voltage effective value correspond to a kind of magnetic flux density with
The variation relation of magnetic field strength；
A kind of change of saturation is chosen in the variation relation of the corresponding magnetic flux density of the different voltages virtual value and magnetic field strength
Change relationship is as the hysteresis loop to be fitted.
4. the method according to claim 3 for determining current transformer J-A simulation model based on Physical Experiment, feature exist
In the alternating voltage of the primary side load different voltages virtual value in the Physical Experiment model determines that different voltages have
The variation relation of valid value corresponding magnetic flux density and magnetic field strength, comprising:
During the corresponding alternating voltage of each voltage effective value is loaded onto the Physical Experiment model primary side, acquire more
Voltage u on group windingpWith the voltage u on sampling resistorR, the sampling resistor and the windings in series；
Determine every group of upAnd uRCorresponding magnetic flux density B and magnetic field strength H, obtains especially by following formula:
Wherein, S and l is respectively the core section product and the average length of magnetic path of Physical Experiment model；im1It is leveled for primary side
Treated electric current,Φ2Magnetic after being leveled surely for what is passed through in the iron core
It is logical；Φ2=Φ1-(max(Φ1)+min(Φ1))/2, Φ1To remove the magnetic flux after the biasing magnetic flux passed through in the iron core, Φ1
=Φ-Φpz, Φ is that integration method measures the initial magnetic flux passed through in the iron core,ΦpzTo bias magnetic flux,T is the AC voltage cycle of primary side load；
According to multiple groups u under each voltage effective valuepAnd uRCorresponding multiple groups B and H value, determine the corresponding magnetic of each voltage effective value
The variation relation of flux density and magnetic field strength.
5. the method according to claim 4 for determining current transformer J-A simulation model based on Physical Experiment, feature exist
In,
It is similar to hysteresis loop to be fitted to the variation relation of H that the evaluation function J (G) describes the corresponding B of new filial generation
Degree, the smaller then degree of fitting of J (G) is higher, and 9 of the corresponding current transformer J-A simulation model of the higher new filial generation of explanation of degree of fitting
The numerical value of parameter is more accurate, and J (G) is determined especially by following formula:
Wherein, HmeasuredFor the magnetic field strength of Physical Experiment model measurement, HsimulatedFor the corresponding magnetic field strength of new filial generation, n is
Sampling u in Physical Experiment modelpAnd uRPoints.
6. the method according to claim 4 for determining current transformer J-A simulation model based on Physical Experiment, feature exist
In the Metropolis criterion describes the new filial generation G for receiving genetic algorithm generation2Probability, pass through following formula determine:
If Metropolis criterion is new explanation G2Better than G1, i.e. J (G2) < J (G1), then centainly receive G2；Otherwise only with P (G2)
Probability receives G2。
7. the side according to any one of claims 1 to 6 for determining current transformer J-A simulation model based on Physical Experiment
Method, which is characterized in that described to determine the method for current transformer J-A simulation model for determining that power grid measures based on Physical Experiment
The J-A simulation model of the different types of current transformer used；The different types of Current Mutual Inductance that the power grid measurement uses
The sectional area S of device iron core is different with average length of magnetic path l.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710384368.6A CN107103163B (en) | 2017-05-26 | 2017-05-26 | A method of current transformer J-A simulation model is determined based on Physical Experiment |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710384368.6A CN107103163B (en) | 2017-05-26 | 2017-05-26 | A method of current transformer J-A simulation model is determined based on Physical Experiment |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107103163A CN107103163A (en) | 2017-08-29 |
| CN107103163B true CN107103163B (en) | 2019-10-25 |
Family
ID=59669098
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710384368.6A Active CN107103163B (en) | 2017-05-26 | 2017-05-26 | A method of current transformer J-A simulation model is determined based on Physical Experiment |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107103163B (en) |
Families Citing this family (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107609298B (en) * | 2017-09-26 | 2020-11-03 | 广东电网有限责任公司电力科学研究院 | Method and device for identifying parameters of Jiles-Atherton model |
| CN107942278B (en) * | 2017-11-16 | 2020-08-14 | 广东电网有限责任公司电力科学研究院 | Simple method and device for calculating transient transmission characteristic of current transformer |
| CN108052017A (en) * | 2017-11-30 | 2018-05-18 | 广东电网有限责任公司电力调度控制中心 | A kind of current transformer Lucas simulation model construction methods |
| CN108038078B (en) * | 2017-12-11 | 2021-03-16 | 囯网河北省电力有限公司电力科学研究院 | Parameter calculation method and device in current transformer detection and terminal equipment |
| CN109543240B (en) * | 2018-10-30 | 2022-09-16 | 西安理工大学 | Current transformer modeling method based on dynamic region saturation J-A theory |
| CN109308403B (en) * | 2018-11-26 | 2022-11-04 | 荀佳钰 | Method for determining J-A simulation model of current transformer based on physical experiment |
| CN109884564B (en) * | 2019-03-22 | 2020-02-14 | 华中科技大学 | Method and device for measuring magnetic characteristics of transformer magnetic core |
| CN112069729B (en) * | 2020-08-28 | 2022-04-22 | 福州大学 | Voltage division method of multi-stage insulation dry type current transformer |
| CN113408232B (en) * | 2021-06-18 | 2023-09-19 | 国网陕西省电力公司电力科学研究院 | Magnetic ring power frequency saturation characteristic simulation method |
| CN117272683B (en) * | 2023-11-13 | 2024-01-26 | 江苏南方永磁科技有限公司 | Model parameter optimization method and system for magnetic hysteresis effect simulation of ferromagnetic material |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN202600121U (en) * | 2012-04-14 | 2012-12-12 | 河北工业大学 | Device for obtaining hysteresis loop of dc biasing of transformer iron-core material |
| CN102967842A (en) * | 2012-10-24 | 2013-03-13 | 重庆大学 | Method for on-line diagnosing gradually-changing fault of electronic current transformers |
| CN103279625A (en) * | 2013-06-15 | 2013-09-04 | 国家电网公司 | Equivalent circuit model building method of transformer iron cores in different structures |
| CN104200055A (en) * | 2014-09-26 | 2014-12-10 | 国家电网公司 | Excitation inrush current simulation method and device for ultra-high voltage transformer under condition of no-load closing |
| CN104991994A (en) * | 2015-06-09 | 2015-10-21 | 国网天津市电力公司 | Transformer field current simulation method based on J-A magnetic hysteresis model |
| CN105302975A (en) * | 2015-11-09 | 2016-02-03 | 国网福建省电力有限公司 | Harmonic transmission modeling method of electromagnetic current transformer |
| CN105740569A (en) * | 2016-02-24 | 2016-07-06 | 国家电网公司 | Current transformer engineering model building method based on transient large current testing technology |
-
2017
- 2017-05-26 CN CN201710384368.6A patent/CN107103163B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN202600121U (en) * | 2012-04-14 | 2012-12-12 | 河北工业大学 | Device for obtaining hysteresis loop of dc biasing of transformer iron-core material |
| CN102967842A (en) * | 2012-10-24 | 2013-03-13 | 重庆大学 | Method for on-line diagnosing gradually-changing fault of electronic current transformers |
| CN103279625A (en) * | 2013-06-15 | 2013-09-04 | 国家电网公司 | Equivalent circuit model building method of transformer iron cores in different structures |
| CN104200055A (en) * | 2014-09-26 | 2014-12-10 | 国家电网公司 | Excitation inrush current simulation method and device for ultra-high voltage transformer under condition of no-load closing |
| CN104991994A (en) * | 2015-06-09 | 2015-10-21 | 国网天津市电力公司 | Transformer field current simulation method based on J-A magnetic hysteresis model |
| CN105302975A (en) * | 2015-11-09 | 2016-02-03 | 国网福建省电力有限公司 | Harmonic transmission modeling method of electromagnetic current transformer |
| CN105740569A (en) * | 2016-02-24 | 2016-07-06 | 国家电网公司 | Current transformer engineering model building method based on transient large current testing technology |
Non-Patent Citations (1)
| Title |
|---|
| 基于蝙蝠算法的启发式智能优化研究与应用;薛菲;《中国博士学位论文全文数据库 信息科技辑》;20170215;正文第4章 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107103163A (en) | 2017-08-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107103163B (en) | A method of current transformer J-A simulation model is determined based on Physical Experiment | |
| CN107656221A (en) | A kind of transformer core remanent magnetism evaluation method based on minor loop's slope | |
| CN105740569B (en) | Current transformer engineering model method for building up based on transient heavy current measuring technology | |
| Chiesa et al. | Transformer model for inrush current calculations: Simulations, measurements and sensitivity analysis | |
| Guerra et al. | Current transformer model | |
| Høidalen et al. | Developments in the hybrid transformer model–Core modeling and optimization | |
| CN111475904A (en) | Method for calculating alternating current loss of low-temperature superconducting magnet | |
| CN109444776A (en) | Three-phase full-bridge converter iron core remanent magnetism measuring method, system and storage medium | |
| Rezaei-Zare et al. | An accurate current transformer model based on Preisach theory for the analysis of electromagnetic transients | |
| CN105182261B (en) | Coil inside core field intensity contactless measurement | |
| Solovev et al. | Current transformers: Transfer functions, frequency response, and static measurement error | |
| CN112131778A (en) | Transformer residual magnetism assessment method based on particle swarm optimization | |
| CN110399677A (en) | Based on the transformer DC magnetic bias simulation method for improving J-A formula under bias state | |
| Yahiou et al. | Modified method for transformer magnetizing characteristic computation and point‐on‐wave control switching for inrush current mitigation | |
| CN109446707A (en) | A kind of Y/ Δ transformer winding turn-to-turn short circuit vibration acceleration calculation method | |
| CN110289595B (en) | Zero sequence overcurrent protection setting method of high-impedance transformer and related device | |
| CN109490729B (en) | Residual magnetism assessment method and system based on variable region integration in voltage attenuation process | |
| Albert et al. | Measurement based transformer modelling approach | |
| Duan et al. | Modeling and experimental validation of a dynamic regional saturation JA model for protective current transformer | |
| CN110749799A (en) | Extra-high voltage transformer direct current magnetic bias equivalent test method and system | |
| Chiesa et al. | Hysteretic iron-core inductor for transformer inrush current modeling in EMTP | |
| Sato et al. | Study on an accurate iron loss calculation method considering the non-uniformity of the magnetic flux density | |
| Zhou et al. | A model considering deep saturation of the iron core for 10 kV potential transformers | |
| Bui et al. | Characterization and modeling of a current transformer working under thermal stress | |
| Azizian | Nonlinear behavior analysis of split-winding dry-type transformer using a new star model and a coupled field-circuit approach |
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 |