CN107103163A - A kind of method that current transformer J A simulation models are determined based on Physical Experiment - Google Patents
A kind of method that current transformer J A simulation models are determined based on Physical Experiment Download PDFInfo
- Publication number
- CN107103163A CN107103163A CN201710384368.6A CN201710384368A CN107103163A CN 107103163 A CN107103163 A CN 107103163A CN 201710384368 A CN201710384368 A CN 201710384368A CN 107103163 A CN107103163 A CN 107103163A
- Authority
- CN
- China
- Prior art keywords
- mrow
- current transformer
- msub
- simulation models
- physical experiment
- 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
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 that current transformer J A simulation models are determined based on Physical Experiment, including：Pass through the coiling N on the iron core of current transformer1Circle coil builds the Physical Experiment model of current transformer；By the primary side on-load voltage in Physical Experiment model, it is determined that hysteresis curve to be fitted, the hysteresis curve is the magnetic flux density curve corresponding with the variation relation of magnetic field intensity on iron core；9 initial parameters of predetermined current transformer J A simulation models, the end value of 9 parameters of current transformer J A simulation models is determined according to 9 initial parameters and hysteresis curve iterative fitting to be fitted；Current transformer J A simulation models are determined according to the end value of 9 parameters of current transformer J A simulation models, J A simulation models are used for the measurement process for the current transformer that simulating grid measurement is used.The present invention can build the practical current transformer Physical Experiment model of different types of industry using the relatively low Physical Experiment electric current of the order of magnitude, with very big promotional value.
Description
Technical field
The present invention relates to Relay Protection Technology in Power System field, determined more particularly, to one kind based on Physical Experiment
The method of current transformer J-A simulation models.
Background technology
In power system, the accurate primary current for flowing through equipment can not be measured, it is necessary to which current transformer is according to electricity
The larger primary current of numerical value is converted to the less secondary current of numerical value by magnetic induction principle by certain no-load voltage ratio, is then measured
Reduction obtains primary current.In fact, because this kind of measuring method and current transformer measurement performance difference, it is mutual by electric current
There is certain error in the primary current and actual primary current that the progress of disease reduction of sensor is obtained.Meanwhile, in recent years, in electric power
During system transient modelling, power network there occurs a lot of because CT saturation causes transmission line of electricity, transformer, generator differential
The accident of false protection, has had a strong impact on power supply reliability and power network safety operation.Study transient state process of electric power system pair
CT saturation, so that the effective means that influence is produced on protection is Digital Simulation.
The simulation model of corresponding current transformer is designed for the current transformer that power network measurement is used, electricity can be simulated
Primary current is passed to the actual progress of disease process of secondary current by analogue simulation close to actual current transformer in current transformer
The transient state and steady-state behaviour of change, more accurately to measure the primary current for flowing through equipment.Based on scholar Jiles and Atherton
Propose on ferromagnetic material phenome-nology theory on the basis of, using Nonlinear System of Equations come it is approximate it is various in the case of iron core excitation spy
Property, a kind of current transformer J-A simulation models are proposed, the model is because can preferably reflect the transient characterisitics quilt of current transformer
It is widely recognized as.
But, the order of magnitude of the power network measurement corresponding electric current to be tested of current transformer used is higher, and such as 2500A is existing
It is difficult to reach the so high order of magnitude during its primary side loading current during the BH curve of the practical current transformer of field test, and electric current
Transformer J-A simulation models need to set up based on on-the-spot test BH curve data, therefore due to actual amount of current to be measured
Level is larger to cause the existing method for setting up current transformer J-A simulation models to lack promotional value.And existing current transformer J-A
J-A model parameters in simulation model are typical default, can only represent the situation of certain single kind of current transformer, not have general
Property.And should in ferromagnetic material properties difference summation current transformer type (TPY, P, PR type etc.) difference, simulation software due to existing
The typical default parameters of model be not fully applied to China produce and be used for current transformer that power network measurement uses, it is necessary to according to
Factually border current transformer redefines corresponding current transformer J-A simulation models.
The content of the invention
For the defect of prior art, it is an object of the invention to solve prior art due to actual electric current number to be measured
Magnitude is larger to cause the existing method for setting up current transformer J-A simulation models to lack promotional value, and existing current transformer
J-A model parameters in J-A simulation models are typical default, can only represent the situation of certain single kind of current transformer, not have
The technical problem of versatility.
To achieve the above object, current transformer J-A simulation models are determined based on Physical Experiment the invention provides one kind
Method, including：Pass through the coiling N on the iron core of current transformer1Circle coil builds the Physical Experiment mould of current transformer
Type, N1The specified once electricity of the rated primary current of the current transformer used according to power network measurement and the Physical Experiment model
Stream is determined；By the primary side on-load voltage in the Physical Experiment model, it is determined that hysteresis curve to be fitted, the magnetic hysteresis is returned
Line is the magnetic flux density curve corresponding with the variation relation of magnetic field intensity on the iron core；Predetermined current transformer J-A is emulated
9 parameters and its initial value of model, according to the initial value of 9 parameters and the hysteresis curve iterative fitting to be fitted
The end value of 9 parameters of the current transformer J-A simulation models is determined, 9 parameters are the description 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 models
A simulation models, the current transformer J-A simulation models are used for the measurement for simulating the current transformer that the power network measurement is used
Process.
It should be noted that " Physical Experiment " that the embodiment of the present invention is referred to refers to building electricity under the conditions of conventional finite
Current transformer model, such as available voltage range of Physical Experiment condition are the voltage in 0~800V, but 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
Electric current in 2500A, i.e. practical power systems is up to kA grades.Wherein, the current transformer that power network measurement is used is used for actual
The measurement of power system.
The limited conditions that the embodiment of the present invention is provided based on Physical Experiment, by iron core winding N1Circle coil ensures physics
The added electric current of experiment is constant in the magnetomotive force of iron core compared to the actual current that power network is measured in the magnetomotive force of iron core so that be based on
Physical Experiment condition is set up J-A simulation models and is possibly realized.In addition, the embodiment of the present invention can be with by being fitted J-A simulation models
The feature of many-side description iron core hysteresis curve so that the reliability for the J-A simulation models set up based on Physical Experiment is carried significantly
Rise.
Alternatively, N1Determined by below equation：
Wherein, Ipn_ctThe rated primary current of the current transformer used for power network measurement, In_testFor the Physical Experiment
The rated primary current of model.
Specifically, I is passed throughpn_ctAnd In_testRatio setting, it is ensured that the current transformer model that Physical Experiment is set up
Magnetomotive force is constant compared to the magnetomotive force under the actual measuring condition of power network.
Alternatively, the alternating voltages of different voltage effective values is loaded in the primary side of the Physical Experiment model, it is determined that not
The variation relation of corresponding magnetic flux density and magnetic field intensity with voltage effective value, the alternating voltage correspondence one of every kind of voltage effective value
Plant the variation relation of magnetic flux density and magnetic field intensity；In the corresponding magnetic flux density of the different voltage effective values and magnetic field intensity
A kind of variation relation of saturation is chosen in variation relation as the hysteresis curve to be fitted.
Alternatively, the mistake of the Physical Experiment model primary side is loaded onto in the corresponding alternating voltage of each voltage effective value
Cheng Zhong, gathers the voltage u on multigroup windingpWith 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 intensity H, is obtained especially by below equation：
Wherein, S and l are respectively the core section product and the average length of magnetic path of Physical Experiment model；im1Carried out for primary side
Electric current after leveling processing,Φ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 that primary side is loaded；According to multigroup u under each voltage effective valuepAnd uR
Corresponding multigroup B and H values, it is determined that the variation relation of each corresponding magnetic flux density of voltage effective value and magnetic field intensity.
It should be noted that B curves corresponding with H variation relation can be described as hysteresis curve.
Alternatively, 9 parameters and its initial value of predetermined current transformer J-A simulation models, according at the beginning of 9 parameters
Initial value and the hysteresis curve iterative fitting to be fitted determine 9 parameter values of the current transformer J-A simulation models
End value, including：
(a) control parameter of blending heredity simulated annealing fitting algorithm is determined：Initial temperature Ts, exit temperature Tf, cooling system
Number k and maximum cooling number of times Mmax；
(b) initial value of 9 parameters of the current transformer J-A simulation models is determined, according to the initial of 9 parameters
Value produces initial population G1, and make cooling number of times add up M=0, operation algebraically N=0；
(c) control temperature T=T is resets；
(d) to initial population G1New filial generation G is produced by the genetic manipulation of genetic algorithm2, the new filial generation G2Corresponding current
9 new parameter values of transformer J-A simulation models, determine G2Corresponding B and H variation relation and the B and H to be fitted
Variation relation between evaluation function J (G), more excellent solution is selected according to Metropolis criterions 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 of future generation, N=N+1, and control temperature is reduced to T=T*k, and k is coefficient of temperature drop；
(f) judge whether control temperature reaches exit temperature Tf, if so, then performing step (g), otherwise perform step (d)；
(g) amendment cooling number of times, M=M+1；
(h) judge whether cooling number of times reaches limit value Mmax, if so, then performing step (i), otherwise；Perform step (c)；
(i) output initial population G1The corresponding current transformer J-A simulation models of newest filial generation 9 parameters, as
The end value of 9 parameters of the current transformer J-A simulation models.
Alternatively, the evaluation function J (G) describes new filial generation corresponding B and H variation relation and magnetic hysteresis to be fitted
The similarity of loop line, the smaller then degrees of fitting of J (G) are 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 intensity of Physical Experiment model measurement, HsimulatedIt is strong for the corresponding magnetic field of new filial generation
Degree, n is the sampling u in Physical Experiment modelpAnd uRPoints.
Alternatively, the Metropolis criterions describe the new filial generation G for receiving genetic algorithm generation2Probability, pass through
Below equation is determined：
If Metropolis criterions are new explanation G2Better than G1, i.e. J (G2) ＜ J (G1), then necessarily receive G2；Otherwise only with P
(G2) probability receive G2。
Alternatively, 9 parameters of the current transformer J-A simulation models 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 descriptions is relation between magnetic flux density B (or magnetization M) and magnetic field intensity H.Should
Theory uses 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.Can many-sided description iron core magnetic by above-mentioned 9 parameters
The feature of hysteresis curves so that the reliability for the J-A simulation models set up based on Physical Experiment is greatly promoted.Alternatively, the base
Determine that the method for current transformer J-A simulation models is used to determine the different types of electric current that power network measurement is used in Physical Experiment
The J-A simulation models of transformer；The sectional area S for the different types of current transformer core that the power network measurement is used and average
Length of magnetic path l is different.
In general, by the contemplated above technical scheme of the present invention compared with prior art, with following beneficial effect
Really：
(1) present invention can be mutual using the relatively low practical electric current of the different types of industry of Physical Experiment electric current structure of the order of magnitude
Sensor Physical Experiment model, with very big promotional value.
(2) BH curve that the present invention can be measured based on the Physical Experiment model of the practical type of different industry, is simulated
9 parameters of J-A simulation models, may be directly applied in present digital simulation software, pass through 9 of J-A simulation models
The feature of parameter many-side description iron core hysteresis curve so that the reliability for the J-A simulation models set up based on Physical Experiment is big
Big lifting.
Brief description of the drawings
Fig. 1 provides the method flow that current transformer J-A simulation models are 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 schematic flow sheet；
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 curve result of the test provided in an embodiment of the present invention and J-A simulation models are run
Result schematic diagram；
Fig. 6 is that PR grades of current transformer hysteresis curve result of the test provided in an embodiment of the present invention and J-A simulation models are transported
Row result schematic diagram；
Fig. 7 is that TPY grades of current transformer hysteresis curve result of the test provided in an embodiment of the present invention and J-A simulation models are transported
Row result schematic diagram.
Embodiment
In order to make the purpose , technical scheme and advantage of the present invention be clearer, it is right below in conjunction with drawings and Examples
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 parameters based on Physical Experiment for the embodiment of the present invention and determines method flow
Schematic diagram；Including step 102- steps 104.
In step 101, pass through the coiling N on the iron core of current transformer1Circle coil builds the physics of current transformer
Experimental model, N1The rated primary current of the current transformer used according to power network measurement and the Physical Experiment model it is specified
Primary current is determined.Alternatively, N1Determined by below equation：
Wherein, Ipn_ctThe rated primary current of the current transformer used for power network measurement, 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 relatively low, and the primary side current accordingly brought is relatively low,
The order of magnitude for the rated current that correspondence Physical Experiment adds is relatively low.When the specifications parameter change of iron core, correspondence is different types of
Industrial practicality current transformer physical model.That is, the embodiment of the present invention can utilize the relatively low Physical Experiment of the order of magnitude
The iron core of electric current and different size builds the practical current transformer Physical Experiment model of different types of industry, is pushed away with very big
Wide value.
In step 102, by the primary side on-load voltage in the Physical Experiment model, it is determined that magnetic hysteresis to be fitted is returned
Line, the hysteresis curve is the magnetic flux density curve corresponding with the variation relation of magnetic field intensity on iron core.
In step 103,9 parameters and its initial value of predetermined current transformer J-A simulation models, according to described 9 ginsengs
Several initial values and the hysteresis curve iterative fitting to be fitted determine 9 ginsengs of the current transformer J-A simulation models
Number；9 parameters are the parameter for describing the iron core hysteresis curve.
In step S104, the current transformer J- is determined according to 9 parameters of the current transformer J-A simulation models
A simulation models, the current transformer J-A simulation models are used for the measurement for simulating the current transformer that the power network measurement is used
Process.
The BH curve that the embodiment of the present invention can be measured based on the Physical Experiment model of the practical type of different industry, simulation
Go out 9 parameters of J-A simulation models, may be directly applied in present digital simulation software.
The problem of different types of industry practicality current transformer J-A models can not accurately being obtained for prior art, 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 actual J-A simulation parameters without type current transformer.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, there is provided a kind of current transformer J-A simulation models based on physical test and parameter fitting algorithm
Parameter acquiring method, this method comprises the following steps, and flow chart is as shown in Figure 2.
(1) the current transformer physical model suitable for physical test is built.
(2) basic parameter such as core area S and length of magnetic path l of test current transformer physical model.
(3) test circuit is built, test obtains the hysteresis curve of the current transformer for parameter fitting.
(4) blending heredity simulated annealing fitting algorithm is based on, the J-A that fitting obtains describing current transformer ferromagnetic characteristic is imitated
True mode parameter.
The current transformer physical model further built in (1) suitable for physical test mainly includes following content：
(11) the first winding number of turn is calculated.The rated primary current I of practical current transformerpn_ct(kilo-ampere), dynamic
Model experiment room can not produce kA grades of primary current, in order to ensure that magnetomotive force is constant, and increase is needed 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 turn is
(12) according to the number of turn is calculated, the coil of the number of turn is wound on current transformer core using close by the way of,
Construct the physical model of the current transformer for Physical Experiment.
(3) hysteresis curve of measurement current transform er comprises the following steps.
(31) test circuit is built, the first siding ring two ends that coiling is completed exchange adjustable voltage is connected in parallel to, and seal in
Sample resistance is to measure electric current.
(32) as shown in figure 3, first winding side correspondence primary current, secondary winding side correspondence secondary current.Pass through regulation
The alternating voltage of the different voltage effective values of primary side input, records voltage and sampling of the different voltage swings next time on winding
Ohmically voltage, respectively upAnd uRIf, upCorresponding passage is " passage 1 ", uRCorresponding passage is " passage 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 pair electric current carries out deburring processing.
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 following a cycle of the point, 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 intensity H is calculated by formula (7) and formula (8).
(35) mapped according to above-mentioned data and can obtain the hysteresis curve of current transformer, that is, measure BH curve.
(36) each bar hysteresis curve under different voltage swings is contrasted, the curve for selecting saturation is intended as parameter
The curve of conjunction.
(4) using blending heredity simulated annealing fitting algorithm (genetic algorithm, simulated annealing), it is mutual that fitting obtains electric current
The flow chart of sensor J-A simulation parameters is as shown in figure 4, specifically include following steps.
(a) input B-H arrays and the control parameter of algorithm：Initial temperature Ts, exit temperature Tf, coefficient of temperature drop k and maximum
Cool number of times Mmax, one group of parameter (T can be provided here according to the flexible arrange parameter of required precisions=10, Tf=1, k=
0.96、Mmax=100)；
(b) initialize.The initial value (default value that set initial value reference software is given) of 9 parameters is set, is produced full
The initial population G of sufficient condition1；The related parameter of algorithm：And make cooling number of times add up M=0, operation algebraically N=0；
(c) control temperature T=T is resets；
(d)G1New filial generation G is produced by the genetic manipulation of genetic algorithm2, to G2Progress, which is mapped, obtains BH curve, calculates it
Evaluation function value J (G), shown in one of which evaluation function such as formula (9), selects more excellent solution, with P according to Metropolis criterions
(G2) probability receive new explanation G2Replace G1；
(e) algorithm enters of future generation, N=N+1, and control temperature is reduced to T=T*k；
(f) judge whether to reach exit temperature Tf, if so, then turning (g), otherwise turn (d)；
(g) amendment cooling number of times, M=M+1；
(h) judge whether cooling number of times reaches limit value Mmax, if so, then turning (i), otherwise turn (c)；
(i) optimum results are exported.
Evaluation function J (G), which is described, is fitted the good and bad degree for obtaining parameter, i.e. emulation gained BH curve and original B-H is bent
The fitting similarity of line, expression formula is provided by formula (9), and the smaller then degrees of fitting of J (G) are higher.
Metropolis criterions describe the new filial generation G for receiving genetic algorithm generation2Probability, provided by formula (10).
If J (G2) ＜ J (G1), represent new explanation G2Better than G1, then G is necessarily received2；Otherwise G is only received with certain probability2.
In one optional example, the Current Mutual Inductance based on physical test and blending heredity simulated annealing fitting algorithm that the present invention is provided
Device J-A simulation parameters acquisition methods comprise the following steps：
(1) practical type (P, PR and TPY type) current transformer physical model of difference suitable for physical test is built.
(11) the first winding number of turn is calculated.The rated primary current of practical current transformer is 2500A, the specified electricity in laboratory
Flow for 10A or so, the coil turn of first winding need to be increased according to magnetomotive force equivalence principle.Then need in current transformer once
The side coiling number of turn is
(12) according to the number of turn calculated, the coil of the number of turn is wound on current transformer core, that is, constructing 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 parameters provided in an embodiment of the present invention based on Physical Experiment determine method
Suitable for the current transformer physical model of the practical type of difference of 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
The different types of current transformer that the correspondence power network measurement of current transformer model is used.Therefore the electric current that power network measurement is used is mutual
Sensor also is understood as practical current transformer.Can according to S and l by the current transformer of different practical type be divided into P, PR and
TPY type current transformers, as shown in table 1.
The current transformer basic parameter of table 1
(3) test obtains the hysteresis curve of P, PR and TPY type current transformer, and selects and obtain the deeper song of degree of saturation
Line as parameter fitting curve, it is shown in solid in such as Fig. 5~Fig. 7.
(4) blending heredity simulated annealing fitting algorithm, the J-A models ginseng of fitting description current transformer ferromagnetic characteristic are based on
Number, as shown in table 2, inputs simulation model, the hysteresis curve emulated is as shown in the dotted line in Fig. 5~Fig. 7.
P grades of table 2, 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 curve result of the test provided in an embodiment of the present invention
With J-A simulation model operation result schematic diagrames；Comparison diagram 5 to Fig. 7 experimental result are visible, provide according to embodiments of the present invention
The goodness of fit of the hysteresis curve that the emulation of J-A simulation models is obtained and the hysteresis curve of actual measurement is higher, illustrates present invention implementation
The application effect of the J-A simulation parameters discrimination method based on physical test and blending heredity simulated annealing fitting algorithm that example is provided
Fruit is good, can effectively be fitted the J-A simulation parameters of different practical types, specific application prospect well.
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, it is not used to
The limitation present invention, any modifications, equivalent substitutions and improvements made within the spirit and principles of the invention etc., it all should include
Within protection scope of the present invention.
Claims (9)
1. a kind of method that current transformer J-A simulation models are determined based on Physical Experiment, it is characterised in that including：
Pass through the coiling N on the iron core of current transformer1Circle coil builds the Physical Experiment model of current transformer, N1According to
The rated primary current of the rated primary current for the current transformer that power network measurement is used and the Physical Experiment model is determined；
By the primary side on-load voltage in the Physical Experiment model, it is determined that hysteresis curve to be fitted, the hysteresis curve
For the curve corresponding with the variation relation of magnetic field intensity of the magnetic flux density on the iron core；
9 parameters and its initial value of predetermined current transformer J-A simulation models, according to the initial value of 9 parameters and institute
The end value that hysteresis curve iterative fitting to be fitted determines 9 parameters of the current transformer J-A simulation models is stated, it is described
9 parameters are the parameter for describing the hysteresis curve；
The current transformer J-A emulation mould is determined according to the end value of 9 parameters of the current transformer J-A simulation models
Type, the current transformer J-A simulation models are used for the measurement process for simulating the current transformer that the power network measurement is used.
2. the method according to claim 1 for determining current transformer J-A simulation models, it is characterised in that N1By following
Formula is determined：
<mrow>
<msub>
<mi>N</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mi>I</mi>
<mrow>
<mi>p</mi>
<mi>n</mi>
<mo>_</mo>
<mi>c</mi>
<mi>t</mi>
</mrow>
</msub>
<msub>
<mi>I</mi>
<mrow>
<mi>n</mi>
<mo>_</mo>
<mi>t</mi>
<mi>e</mi>
<mi>s</mi>
<mi>t</mi>
</mrow>
</msub>
</mfrac>
</mrow>
Wherein, Ipn_ctThe rated primary current of the current transformer used for power network measurement, In_testFor the Physical Experiment model
Rated primary current.
3. the method according to claim 1 for determining current transformer J-A simulation models, it is characterised in that described to pass through
In the primary side on-load voltage of the Physical Experiment model, it is determined that hysteresis curve to be fitted, including：
The alternating voltage of different voltage effective values is loaded in the primary side of the Physical Experiment model, different voltage effective values are determined
The variation relation of corresponding magnetic flux density and magnetic field intensity, a kind of alternating voltage magnetic flux density of correspondence of every kind of voltage effective value with
The variation relation of magnetic field intensity；
A kind of change of saturation is chosen in the variation relation of the corresponding magnetic flux density of the different voltage effective values and magnetic field intensity
Change relation is used as the hysteresis curve to be fitted.
4. the method according to claim 1 for determining current transformer J-A simulation models, it is characterised in that described in institute
The primary side for stating Physical Experiment model loads the alternating voltage of different voltage effective values, determines the corresponding magnetic of different voltage effective values
The variation relation of flux density and magnetic field intensity, including：
During the corresponding alternating voltage of each voltage effective value is loaded onto the Physical Experiment model primary side, gather many
Voltage u on the 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 intensity H, is obtained especially by below equation：
<mrow>
<mi>B</mi>
<mo>=</mo>
<mfrac>
<msub>
<mi>&Phi;</mi>
<mn>2</mn>
</msub>
<mi>S</mi>
</mfrac>
</mrow>
<mrow>
<mi>H</mi>
<mo>=</mo>
<mfrac>
<msub>
<mi>N</mi>
<mn>1</mn>
</msub>
<mi>l</mi>
</mfrac>
<msub>
<mi>i</mi>
<mrow>
<mi>m</mi>
<mn>1</mn>
</mrow>
</msub>
</mrow>
Wherein, S and l are respectively the core section product and the average length of magnetic path of Physical Experiment model；im1Leveling is carried out for primary side
Electric current after processing,Φ2For the magnetic of the leveling that passes through in the iron core after steady
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 that primary side is loaded；
According to multigroup u under each voltage effective valuepAnd uRCorresponding multigroup B and H values, it is determined that the corresponding magnetic of each voltage effective value
The variation relation of flux density and magnetic field intensity.
5. the method according to claim 1 for determining current transformer J-A simulation models, it is characterised in that predetermined current
9 parameters and its initial value of transformer J-A simulation models, according to 9 initial parameter values and the magnetic hysteresis to be fitted
Loop line iterative fitting determines the end value of 9 parameter values of the current transformer J-A simulation models, including：
(a) control parameter of blending heredity simulated annealing fitting algorithm is determined：Initial temperature Ts, exit temperature Tf, coefficient of temperature drop k
With maximum cooling number of times Mmax；
(b) initial value of 9 parameters of the current transformer J-A simulation models is determined, is produced according to the initial value of 9 parameters
Raw initial population G1, and make cooling number of times add up M=0, operation algebraically N=0；
(c) control temperature T=T is resets；
(d) to initial population G1New filial generation G is produced by the genetic manipulation of genetic algorithm2, the new filial generation G2Corresponding current mutual inductance
9 parameter values of device J-A simulation models, determine G2Between corresponding B and H variation relation and the hysteresis curve to be fitted
Evaluation function J (G), more excellent solution is selected according to Metropolis criterions based on J (G), with P (G2) probability receive new explanation G2Take
For G1, the P (G2) according to J (G1) and J (G2) determine；
(e) algorithm enters of future generation, N=N+1, and control temperature is reduced to T=T*k, and k is coefficient of temperature drop；
(f) judge whether control temperature reaches exit temperature Tf, if so, then performing step (g), otherwise perform step (d)；
(g) amendment cooling number of times, M=M+1；
(h) judge whether cooling number of times reaches limit value Mmax, if so, then performing step (i), otherwise；Perform step (c)；
(i) output initial population G1The corresponding current transformer J-A simulation models of newest filial generation 9 parameters, be used as the electricity
The end value of 9 parameters of current transformer J-A simulation models.
6. the method according to claim 5 for determining current transformer J-A simulation models, it is characterised in that
It is similar to hysteresis curve to be fitted to H variation relation that the evaluation function J (G) describes the corresponding B of new filial generation
Degree, the smaller then degrees of fitting of J (G) are higher, 9 of the corresponding current transformer J-A simulation models 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 below equation：
<mrow>
<mi>J</mi>
<mrow>
<mo>(</mo>
<mi>G</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msqrt>
<mrow>
<mfrac>
<mn>1</mn>
<mi>n</mi>
</mfrac>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>s</mi>
<mi>a</mi>
<mi>m</mi>
<mi>p</mi>
<mi>l</mi>
<mi>e</mi>
<mi>s</mi>
</mrow>
</munder>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>H</mi>
<mrow>
<mi>m</mi>
<mi>e</mi>
<mi>a</mi>
<mi>s</mi>
<mi>u</mi>
<mi>r</mi>
<mi>e</mi>
<mi>d</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>H</mi>
<mrow>
<mi>s</mi>
<mi>i</mi>
<mi>m</mi>
<mi>u</mi>
<mi>l</mi>
<mi>a</mi>
<mi>t</mi>
<mi>e</mi>
<mi>d</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
</mrow>
Wherein, HmeasuredFor the magnetic field intensity of Physical Experiment model measurement, HsimulatedFor the corresponding magnetic field intensity of new filial generation, n is
Sampling u in Physical Experiment modelpAnd uRPoints.
7. the method according to claim 5 for determining current transformer J-A simulation models, it is characterised in that described
Metropolis criterions describe the new filial generation G for receiving genetic algorithm generation2Probability, determined by below equation：
<mrow>
<mi>P</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>G</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mrow>
<mi>J</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>G</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mo><</mo>
<mi>J</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>G</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mfrac>
<mrow>
<mi>J</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>G</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>J</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>G</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mi>T</mi>
</mfrac>
</mrow>
</msup>
</mtd>
<mtd>
<mrow>
<mi>J</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>G</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>&GreaterEqual;</mo>
<mi>J</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>G</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
2
If Metropolis criterions are new explanation G2Better than G1, i.e. J (G2) ＜ J (G1), then necessarily receive G2；Otherwise only with P (G2)
Probability receives G2。
8. the method according to claim 5 for determining current transformer J-A simulation models, it is characterised in that the electric current
9 parameters of transformer J-A simulation models include：Saturation flux Ms, region pinning coefficient k, region can disturb coupling between coefficient c, domain
Close factor alpha, curve shape coefficient a1、a2、a3, b and correction factor β.
9. the method for the determination current transformer J-A simulation models according to any one of claim 1 to 8, it is characterised in that
It is described to determine that the method for current transformer J-A simulation models is used to determine the different type that power network measurement is used based on Physical Experiment
Current transformer J-A simulation models；The sectional area S for the different types of current transformer core that the power network measurement is used
It 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 true CN107103163A (en) | 2017-08-29 |
| CN107103163B 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) |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107942278A (en) * | 2017-11-16 | 2018-04-20 | 广东电网有限责任公司电力科学研究院 | A kind of simple and easy method and device of calculating current transformer transient Transfer characteristic |
| CN108038078A (en) * | 2017-12-11 | 2018-05-15 | 囯网河北省电力有限公司电力科学研究院 | Calculation method of parameters, device and terminal device in detecting current transformer |
| CN108052017A (en) * | 2017-11-30 | 2018-05-18 | 广东电网有限责任公司电力调度控制中心 | A kind of current transformer Lucas simulation model construction methods |
| CN109308403A (en) * | 2018-11-26 | 2019-02-05 | 荀佳钰 | A method of current transformer J-A simulation model is determined based on Physical Experiment |
| CN109543240A (en) * | 2018-10-30 | 2019-03-29 | 西安理工大学 | A kind of current transformer modeling method based on dynamic area saturation J-A theory |
| CN109884564A (en) * | 2019-03-22 | 2019-06-14 | 华中科技大学 | A kind of magnetic core of transformer magnetic characteristic measurement method and device |
| CN107609298B (en) * | 2017-09-26 | 2020-11-03 | 广东电网有限责任公司电力科学研究院 | Method and device for identifying parameters of Jiles-Atherton model |
| CN112069729A (en) * | 2020-08-28 | 2020-12-11 | 福州大学 | Voltage division method of multi-stage insulation dry type current transformer |
| CN113408232A (en) * | 2021-06-18 | 2021-09-17 | 国网陕西省电力公司电力科学研究院 | Magnetic ring power frequency saturation characteristic simulation method |
| CN117272683A (en) * | 2023-11-13 | 2023-12-22 | 江苏南方永磁科技有限公司 | 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 |
|---|
| 薛菲: "基于蝙蝠算法的启发式智能优化研究与应用", 《中国博士学位论文全文数据库 信息科技辑》 * |
Cited By (17)
| 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 |
| CN107942278A (en) * | 2017-11-16 | 2018-04-20 | 广东电网有限责任公司电力科学研究院 | A kind of simple and easy method and device of calculating current transformer transient Transfer characteristic |
| CN108052017A (en) * | 2017-11-30 | 2018-05-18 | 广东电网有限责任公司电力调度控制中心 | A kind of current transformer Lucas simulation model construction methods |
| CN108038078A (en) * | 2017-12-11 | 2018-05-15 | 囯网河北省电力有限公司电力科学研究院 | Calculation method of parameters, device and terminal device in detecting current transformer |
| 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 |
| CN109543240A (en) * | 2018-10-30 | 2019-03-29 | 西安理工大学 | A kind of current transformer modeling method based on dynamic area saturation J-A theory |
| CN109308403A (en) * | 2018-11-26 | 2019-02-05 | 荀佳钰 | A method of current transformer J-A simulation model is determined based on Physical Experiment |
| 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 |
| CN109884564A (en) * | 2019-03-22 | 2019-06-14 | 华中科技大学 | A kind of magnetic core of transformer magnetic characteristic measurement method and device |
| CN112069729A (en) * | 2020-08-28 | 2020-12-11 | 福州大学 | Voltage division method of multi-stage insulation dry type current transformer |
| CN112069729B (en) * | 2020-08-28 | 2022-04-22 | 福州大学 | Voltage division method of multi-stage insulation dry type current transformer |
| CN113408232A (en) * | 2021-06-18 | 2021-09-17 | 国网陕西省电力公司电力科学研究院 | Magnetic ring power frequency saturation characteristic simulation method |
| CN113408232B (en) * | 2021-06-18 | 2023-09-19 | 国网陕西省电力公司电力科学研究院 | Magnetic ring power frequency saturation characteristic simulation method |
| CN117272683A (en) * | 2023-11-13 | 2023-12-22 | 江苏南方永磁科技有限公司 | Model parameter optimization method and system for magnetic hysteresis effect simulation of ferromagnetic material |
| CN117272683B (en) * | 2023-11-13 | 2024-01-26 | 江苏南方永磁科技有限公司 | Model parameter optimization method and system for magnetic hysteresis effect simulation of ferromagnetic material |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107103163B (en) | 2019-10-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107103163A (en) | A kind of method that current transformer J A simulation models are determined based on Physical Experiment | |
| CN105740569B (en) | Current transformer engineering model method for building up based on transient heavy current measuring technology | |
| Pedra et al. | PSPICE computer model of a nonlinear three-phase three-legged transformer | |
| CN107656221A (en) | A kind of transformer core remanent magnetism evaluation method based on minor loop's slope | |
| CN107317331A (en) | A kind of extra-high voltage transformer axial symmetry D.C. magnetic biasing simulation model | |
| Bhowmick et al. | Online detection of an interturn winding fault in single-phase distribution transformers using a terminal measurement-based modeling technique | |
| CN111475904A (en) | Method for calculating alternating current loss of low-temperature superconducting magnet | |
| Solovev et al. | Current transformers: Transfer functions, frequency response, and static measurement error | |
| Arastou et al. | New method for synchronous generator parameters estimation using load rejection tests data considering operational limitations | |
| CN107958125A (en) | Low frequency model modeling method and system in a kind of three-phase transformer electro-magnetic transient | |
| CN105182261B (en) | Coil inside core field intensity contactless measurement | |
| Kül et al. | Coreloss Estimation via Long Short-Term Memory Model (LSTM) of Dry-Type Transformer based on FEA | |
| CN114355190A (en) | Doubly-fed motor turn-to-turn short circuit fault detection method capable of switching sliding-mode observer | |
| CN110749799B (en) | Extra-high voltage transformer direct current magnetic bias equivalent test method and system | |
| CN109490729B (en) | Residual magnetism assessment method and system based on variable region integration in voltage attenuation process | |
| Duan et al. | Modeling and experimental validation of a dynamic regional saturation JA model for protective current transformer | |
| Zhou et al. | A model considering deep saturation of the iron core for 10 kV potential transformers | |
| Lesniewska et al. | Influence of the number of core air gaps on transient state parameters of TPZ class protective current transformers | |
| Azizian | Nonlinear behavior analysis of split-winding dry-type transformer using a new star model and a coupled field-circuit approach | |
| Chiesa et al. | Hysteretic iron-core inductor for transformer inrush current modeling in EMTP | |
| Carrander | Magnetizing currents in power transformers: Measurements, simulations, and diagnostic methods | |
| Hane | Hysteresis Modeling for Power Magnetic Devices Based on Magnetic Circuit Method | |
| Guo et al. | Modeling of practical protective ct based on lucas model | |
| Yin et al. | Modeling and analysis for practical CT based on transient test and parameter identification | |
| Zhang et al. | Investigation of harmonic and global loss of three-phase transformer based on a permeance capacitance analogy model |
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 |