CN104158189B - Electrified railway traction power supply load modeling method based on parameter identification - Google Patents

Electrified railway traction power supply load modeling method based on parameter identification Download PDF

Info

Publication number
CN104158189B
CN104158189B CN201410442604.1A CN201410442604A CN104158189B CN 104158189 B CN104158189 B CN 104158189B CN 201410442604 A CN201410442604 A CN 201410442604A CN 104158189 B CN104158189 B CN 104158189B
Authority
CN
China
Prior art keywords
load
power
electric
model
electric railway
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
Application number
CN201410442604.1A
Other languages
Chinese (zh)
Other versions
CN104158189A (en
Inventor
李庆生
邓朴
陈红坤
王丰元
胡畔
赵庆明
农静
张裕
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
GRID PLANNING RESEARCH CENTER OF GUIZHOU GRID Co
Original Assignee
GRID PLANNING RESEARCH CENTER OF GUIZHOU GRID Co
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by GRID PLANNING RESEARCH CENTER OF GUIZHOU GRID Co filed Critical GRID PLANNING RESEARCH CENTER OF GUIZHOU GRID Co
Priority to CN201410442604.1A priority Critical patent/CN104158189B/en
Publication of CN104158189A publication Critical patent/CN104158189A/en
Application granted granted Critical
Publication of CN104158189B publication Critical patent/CN104158189B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

The invention discloses an electrified railway traction power supply load modeling method based on parameter identification. The method includes the steps that an identification analytical expression which is combined with the characteristics of an electrified railway load locomotive in the running process to determine an electrified railway traction power supply load model is built, a typical electrified railway traction station is selected for testing, measured data of electrified railway traction power supply loads are analyzed and preprocessed, typical parameters in the electrified railway traction power supply load model are identified, and analysis validation is performed on accuracy of the electrified railway traction power supply load model by the adoption of a corresponding example. The impact effect generated when the electrified railway traction power supply loads have access to a power grid can be described better, grid planning and electrified railway development planning are supported, and the degree of coordination between grid planning and electrified railway development planning is increased.

Description

A kind of electric railway traction supply load modeling method based on parameter identification
Technical field
Patent of the present invention relates to a kind of for electrization railway draught load modeling method, particularly relates to a kind of electric railway traction supply load modeling method based on parameter identification.
Background technology
Electric railway is the mode of transportation that current China gives priority to, and it can improve railway transport capacity, improve railway operation, is also conducive to realizing the reasonable distribution of resource simultaneously, cuts operating costs, preserves the ecological environment.State Council approved Mid-term and Long-term Plan of the Railway Network explicitly points out: to the year two thousand twenty, and China railways total kilometrage will reach 120000km, and wherein electric railway is 60000km, and Railway Electric rate is about 60%, and the freight volume proportion born is more than 80%.Along with state-owned railroads 6 speed raisings, train speed is more and more faster, separate unit locomotive power is increasing, the particularly development energetically of high-speed electric railway in recent years, such as Beijing-Tianjin inter-city railway, Shi-Tai passenger line, expensive wide Line for Passenger Transportation, Zheng Xi Line for Passenger Transportation, Shanghai and Nanjing inter-city passenger rail, Shanghai-Hangzhou inter-city passenger rail and Beijing-Shanghai Line for Passenger Transportation is open-minded, a large amount of speed per hour 250,300, the CRH type motor train unit on-line running of 350km/h.
The Large scale construction development of electric railway, one side is had higher requirement to the quality of power supply of utility network and power supply reliability; On the other hand, electrization railway draught load is as typical uneven, the impact load of one, there is running status complicated and changeable and randomness is very large, put into operation and can produce a large amount of harmonic waves, negative-sequence current, inject electric power system by electric railway traction transformer station, certain harm is produced to utility network safety and the quality of power supply.
The shock effect of electric railway traction supply load access electrical network, need planning, design time carry out the research of being correlated with and calculate, to determine its impact on Power Network Transient Stability.At present, relevant, systematic modeling method is less, constrains the coordinated development of electrical network and electric railway to a certain extent.Therefore, for ensureing power grid security, meet electric railway great-leap-forward development requirement, boosting China sustain economic increases fast, need to carry out to this problem the impact load that the deep research of system powers to different electric railway traction to simulate, calculate, ensure safe operation of electric network while ensureing electric railway electricity consumption reliability, improve Electric Power Network Planning and electric railway development plan coordination degree.
Summary of the invention
For above-mentioned prior art, the technical problem that the present invention solves is to provide a kind of electric railway traction supply load modeling method based on parameter identification, can better be described the shock effect of electric railway traction supply load access electrical network, support Electric Power Network Planning and electric railway development plan, raising Electric Power Network Planning and electric railway develop draws a coordination degree.
For solving the problem, a kind of electric railway traction supply load modeling method based on parameter identification of the present invention comprises step:
Step 1, in conjunction with the feature in electric railway load locomotive traveling process, directly assemble on contact line supply arm with direct current machine motor load and synthetic load, determine AC-DC electric locomotive Mathematical Modeling, directly assemble on contact line supply arm with three-phase alternating-current motor load and synthetic load, determine to hand over orthogonal type electric locomotive Mathematical Modeling, and AC-DC electric locomotive Mathematical Modeling is comprehensive with friendship orthogonal type electric locomotive Mathematical Modeling, set up electric railway part throttle characteristics model, in conjunction with the impact characteristics of electric railway traction supply load, determine the model of the impact characteristics of electric railway traction supply load, and insertion switch coefficient obtains by formula
Distinguishing of the electric railway electric locomotive model determined misses analytical expression model, in formula: with for the switching coefficient of three-phase alternating-current motor load and DC motor load, for the meritorious and reactive power that three-phase alternating-current motor absorbs, what absorb for DC motor has power; , for the meritorious and reactive power that synthetic load absorbs, electric railway traction is powered the model of negative impact characteristics;
Its AC-DC electric locomotive Mathematical Modeling is:
In formula; for the electric current of traction electric machine, for voltage transformation coefficient, for the initial active power proportionality coefficient of DC traction motor, for model is initially always gained merit and inactivity, for the indicial response of load, for three-phase side voltage effective value, for state voltage;
Its friendship orthogonal type electric locomotive Mathematical Modeling is:
In formula: U is three-phase side voltage effective value, for steady state voltage, δ is merit angle, for the indicial response of load, for the initial active power of synthetic load, the initial reactive power of synthetic load, X ' is transient state reactance, and E ' is electromotive force after transient state reactance;
Electric railway part throttle characteristics model is:
P i=k 1P j+k 2P z+P s
Q i=k 1Q j+Q s
In formula: with for the switching coefficient of three-phase alternating-current motor load and DC motor load, for the active power that three-phase alternating-current motor absorbs, for the active power that DC motor absorbs, for the meritorious and reactive power that synthetic load absorbs, the reactive power that three-phase alternating-current motor absorbs;
The model of the impact characteristics of electric railway traction supply load is:
In formula: U is three-phase side voltage effective value, for steady state voltage, be respectively and gain merit and the variation coefficient of reactive power impact power with voltage;
Step 2, set up the transformation model between the three-phase side voltage of electric power system and two phase voltages of electric railway traction power supply system, in conjunction with electric railway stating with load, acceleration and the feature of advancing across electric phase-splitting, mathematical modeling is carried out to the power demand initiative feature of electric railway load, and in conjunction with the identification analytical expression model of electric railway electric locomotive model, set up by formula
The electric railway impact load model determined, in formula for the identification analytical expression model of electric railway electric locomotive model, for the power demand initiative feature Mathematical Modeling of electric railway load; for waiting to distinguish parameter by mistake;
The modeling method of the power demand initiative feature Mathematical Modeling of described electric railway load is: the changed power process power of locomotive active demand being equivalent to approximate trapezoid, the Mathematical Modeling of power initiative demand utilization step function represents, the model of foundation is by formula
Determine, in formula: be respectively parameter to be identified, for unit step function; with be respectively the power that in trapezoidal changed power figure, curve level, rising and sloping portion are corresponding; for three superpose on a timeline after compressive stroke power;
Step 3, choose some electric railway traction stations and test, and analyze and preliminary treatment the measured data of electric railway traction supply load, its processing method comprises following steps:
Step 31, choosing some substation bus bars for power supply for locomotive station is measurement point, using one group that measures at this bus place typical and complete impact characteristics curve as analytic target, the change curve of voltage, electric current and apparent power in Measuring Time section,
Step 32, adopts Wavelet Decomposition Technology to carry out filtering and denoising to measured data, completes the pretreatment work of measured data;
Step 4, to the canonical parameter in electric railway traction supply load model, carries out parameter identification, and adopts corresponding example, carries out analysis verification to the accuracy of electric railway traction supply load model, needs the canonical parameter carrying out parameter identification to be: , L and R is respectively inductance and the resistance in traction electric machine loop, and J is the moment of inertia of traction electric machine for the initial active power proportionality coefficient of DC traction motor, for transient state reactance, X is synchronous reactance, for time constant, for the initial active power proportionality coefficient of three-phase alternating-current motor, be respectively and gain merit and the variation coefficient of reactive power impact power with voltage; When carrying out parameter identification to the canonical parameter in electric railway traction supply load model, according to Parameter Sensitivity Analysis result, the improvement PSO algorithm choosing linear decrease weights strategy (LDIW) is distinguished by mistake typical independent parameter.
Beneficial effect of the present invention:
Research at present for the modeling of electric power system impact load is still in the immature stage, not yet forms systematized impact load Modeling Theory.When power system planning, owing to lacking dynamic impulsion load model accurately, static models are often adopted to replace; And the impact load model that actual electric network operational mode calculates employing is mostly static or induction motor static characteristic load model in parallel.On the other hand, for electric railway traction load modeling, current research both domestic and external is not still deeply, the research of its impact characteristics aspect of shorter mention, for this reason, the electric railway traction supply load modeling method that the present invention is based on parameter identification carries out modeling analysis for electric railway impact load, mainly has following feature:
(1) carry out structure and Analysis on Mechanism for electric railway traction power supply system and electric locomotive, summarize the characteristic of electrization railway draught load, do phenomenon analysis and explanation in particular for its impact characteristics;
(2) analyze in conjunction with electric iron impact load actual measurement load curve, the phenomenon that load curve shows is made explanations; Propose to adopt wavelet analysis method to do denoising to actual measurement load curve, the data after denoising are used as identification of Model Parameters;
(3) set up AC-DC respectively towards electrical network PCC point and hand over the part throttle characteristics model of orthogonal type electric locomotive, and by both comprehensively for adding the analytical expression describing electric iron impact load impact characteristics after unified electrified railway traction load characteristic model, form electric iron impact load collective model structure, formulate identification of Model Parameters strategy and implemented.The feasibility of modeling method, the correctness of institute's established model and validity is herein demonstrated by sample calculation analysis;
(4) the initiative feature for impact load makes improvements model, and itself and master mould are compared analysis.
The present invention is by setting up the electric railway traction supply load modeling method based on parameter identification, can better be described the shock effect of electric railway traction supply load access electrical network, support Electric Power Network Planning and electric railway development plan, improve Electric Power Network Planning and electric railway development plan coordination degree, ensure that power grid security, meet electric railway great-leap-forward development requirement.
Accompanying drawing explanation
Fig. 1 is electrization railway draught load isoboles;
Fig. 2 is voltage, electric current, apparent power measured curve figure;
Fig. 3 is the curve chart of curve after denoising;
Fig. 4 is program overall framework figure;
Fig. 5 is impact load model realization flow chart;
Fig. 6 is the simulated power curve chart of load model I;
Fig. 7 is the simulated power curve chart of load model II;
Fig. 8 is the solution of the present invention block diagram.
Embodiment
Be described in further detail the specific embodiment of the present invention below in conjunction with accompanying drawing, as shown in Figure 8, concrete steps are as follows for the solution of the present invention block diagram:
Foundation can describe the electric railway traction supply load model of load impact feature, and in conjunction with the feature in electric railway load locomotive traveling process, the expression formula of electric railway traction supply load model is improved, establishes the identification analytical expression of electric railway traction supply load model.
According to electric railway traction supply load model towards electric power system PCC point feature, set up the transformation model between the three-phase side voltage of electric power system and two phase voltages of electric railway traction power supply system.Traction Station, by traction transformer, becomes two phase voltages from three-phase voltage, the conversion namely between three-phase side voltage and contact line supply arm voltage.
Transformation relation between three-phase side voltage and contact line supply arm voltage is:
Wherein, for contact line supply arm voltage, for traction transformer is to the no-load voltage ratio of supply arm port; for traction transformer port connection angle; for electric power system side three-phase voltage.If c is three-phase side voltage effective value, then contact line supply arm voltage effective value can be obtained by following formula.
Definition for contact line supply arm voltage is to the conversion coefficient of three-phase side voltage effective value, then;
Adopt the load model of " direct current machine synthetic load in parallel " to describe AC-DC electric locomotive part throttle characteristics, direct current machine load and synthetic load are directly assembled on contact line supply arm, sets up AC-DC electric locomotive Mathematical Modeling.
AC-DC electric locomotive load is formed primarily of traction motor, auxiliary equipment and driver's cabin, and the traction electric machine of AC-DC electric locomotive is similar to DC motor.Therefore AC-DC electric locomotive load adopts the model of " platform load is combined in traction electric machine parallel connection " to describe.The state equation of traction electric machine and the active power equation of absorption as follows:
In above formula: for traction electric machine loop terminal voltage, with for traction electric machine factory justice equivalent parameters, with for electric current and the rotating speed of traction electric machine, with be respectively inductance and the resistance in traction electric machine loop; rotation for traction electric machine is used heavy. for mechanical load torque, the quadratic function models of rotating speed is adopted to describe mechanical load torque, for traction motor load rate.
If PCC point busbar voltage is U, be then tied to form vertical just like ShiShimonoseki:
for voltage transformation coefficient, obtained by formula 1.Formula 3 being substituted into formula 2 obtains the traction motor circuit state equation of system three-phase side voltage as excitation:
For reducing identified parameters, improve the identification efficiency to model major parameter and precision, a kind of integrated load model that describes can adopt power function form to represent.
Wherein, reflecting the indicial response of load, is the object of load model parameters identification; determined by steady state operating conditions.
For improving the adaptivity of model, introduce DC traction motor initially meritorious and reactive power proportionality coefficient then have:
Wherein, for model is initially always gained merit and reactive power, therefore, the model structure of AC-DC electric locomotive is made up of formula 4, formula 5 and formula 6.Its power always exports as shown in Equation 7.
Adopt the load model of " three-phase alternating-current motor synthetic load in parallel " to describe and hand over orthogonal type electric locomotive part throttle characteristics, three-phase alternating-current motor load and synthetic load are directly assembled on contact line supply arm, set up and hand over orthogonal type electric locomotive Mathematical Modeling.
Hand over orthogonal type electric locomotive traction drive system primarily of traction transformer, traction convertor and traction electric machine composition, its traction electric machine adopts threephase asynchronous.Locomotive load dynamic part adopts 3 rank electromechanical transient motor models to describe; Comprehensive part adopts the exponential function integrated load model ignoring frequency change impact to describe.
3 rank asynchronous motor state equations:
For improving the adaptivity of model, introduce the initial active power proportionality coefficient of three-phase alternating-current motor then have:
Wherein, for the initial total active power of model, for the initial active power of synthetic load.
Based on above-mentioned analysis, the foundation of electric locomotive model is the transformer station PCC point towards powering for electric railway traction station, and transmission line and tractive power supply system loss consider into by model, so foregoing model can be used for describing the characteristic that electric railway traction power supply system and electric locomotive load change with system voltage.Railway system's electric locomotive is of a great variety, and on different circuits, the type of electric locomotive is not quite similar.Therefore, be necessary aforementioned two class models to integrate, form the electric railway part throttle characteristics model possessing diversified descriptive power.Load model isoboles after integrating as shown in Figure 1.
The switching coefficient of definition three-phase alternating-current motor load and DC motor load is respectively with wherein value 0 or 1, with be asynchronously 0.Following model structure is obtained according to formula 7 and formula 11:
P i=k 1P j+k 2P z+P s
(formula 12)
Q i=k 1Q j+Q s
A, when time, be AC-DC electric locomotive load model;
B, when time, for handing over orthogonal type electric locomotive load model;
C, when time, be AC-DC and friendship orthogonal type locomotive mixing load model.
In formula, for always gaining merit and reactive power; for the meritorious and reactive power that three-phase alternating-current motor absorbs; for the active power that DC motor absorbs; , for the meritorious of comprehensive platform load absorption and reactive power.
Because electric railway load is a kind of typical impact load, in Model in Time Domain, in order to show its impact characteristics, introducing equation and can react actual electric railway changed power process preferably.Therefore a model for the impact characteristics of electric railway traction supply load is described, can by formula
Determine.Wherein, for steady state voltage, be respectively and gain merit and the variation coefficient of reactive power impact power with voltage.
Formula 12 is utilized comprehensively to obtain the collective model of electric railway impact load shown in formula 14 structure with the model of the impact characteristics of electric railway traction supply load.
Owing to being divided into two aspects to the active demand of power in the actual traveling process of electric locomotive:
A, startup and accelerator; Electric locomotive is in startup and accelerator, and great change can occur in traction motor load torque.
B, across electric phase-splitting; Electric locomotive divides phase time with larger speed across electricity, can put on electric power system instantaneously with power rush.
Therefore, a kind of power demand initiative feature to electric railway load is carried out Mathematical Modeling Methods and is, first the power of locomotive active demand is equivalent to the changed power process of approximate trapezoid, is then determined the mathematic(al) structure of this trapezoidal changed power curve by System Discrimination.
Power initiative demand can utilize the Mathematical Modeling of step function to be expressed as follows:
Wherein, be respectively parameter to be identified, for unit step function; with be respectively the power that in trapezoidal changed power figure, curve level, rising and sloping portion are corresponding; for three superpose on a timeline after compressive stroke power.If consider that active power and reactive power are impacted respectively, formula 15 can be utilized to pass through power factor (PF) try to achieve.
The initiative power characteristic of trapezoid model to electric railway load is adopted to be described, can accurately reflect part throttle characteristics and the impact characteristics of load: the process not only sharply rising for electric locomotive load power demand and decline has good description, and fitting effect is desirable within the time period of the small adjustment of its rate of doing work.
The power characteristic that the power characteristic described by formula 14 and formula 15 describe superposes, in order to improve the mathematical description for electric railway impact load.The electric railway impact load model that order improves is shown in formula 14, load model is formula 15 representation models are the then electric railway impact load model structure of 16 expression improvement with the formula.
Wherein, for waiting to distinguish parameter by mistake.Complete model can be obtained by the parameter identification strategy implementing previously to have formulated.
Consider the identifiability of model, according to Parameter Sensitivity Analysis result, choose following parameter as independent parameter to be identified:
In model, other parameter choose representative value, as shown in table 1.
The representative value of table 1 model parameter
Choose typical electric railway traction station to test, and to the measured data analysis of electric railway traction supply load and preliminary treatment.
A kind of measured data analysis and preprocess method of electric railway traction supply load are, to the multiple electric railway traction transformer station in somewhere and measure the bus of side system transformer station, the change curve of voltage, electric current and apparent power in recording period.Using one group that measures at this bus place typical case and also complete impact characteristics curve right as analysis, as shown in Figure 2.In figure, abscissa is sampling number, and sample frequency is power frequency, each sampling interval 0.02s, and amount to duration 48s, ordinate is perunit value, and fiducial value is the value of each electric parameters in 0 moment.
From load composition, in the load data that PCC point records, not only comprise low-pressure side tractive power supply system load, and contain other customer charges a certain proportion of, this makes to contain certain noise in measured signal.From the feature of impact load, a large amount of noises can cover the impact characteristics of impact load to a certain extent.Therefore, in order to eliminate the adverse effect that above-mentioned factor is brought as far as possible, preliminary treatment must be carried out to measured data, to improve the precision of load modeling and parameter identification.
Wavelet analysis technology is adopted to carry out filtering to measured data.Adopt db5 small echo to carry out 5 layers of decomposition to measured data, obtain Fig. 3 measured signal curve, complete the pretreatment work of measured data.
To the canonical parameter in electric railway traction supply load model, carry out parameter identification, and adopt corresponding example, analysis verification is carried out to the accuracy of electric railway traction supply load model.
A kind ofly to the method that the canonical parameter in electric railway traction supply load model carries out parameter identification be, because model number of parameters to be identified is relatively many, in order to improve identification accuracy and Searching efficiency, choose PSO algorithm (population), and adopt linear decrease weights strategy (LDIW) to improve PSO algorithm.
Obtain model parameter by electric railway impact load model shown in programming realization formula 15 in Matlab.Program overall framework as shown in Figure 4.Program is roughly divided into two-layer, is program internal layer in dotted line frame.Skin is the identification of Model Parameters process based on LDIW-PSO, and internal layer is the detailed realization to formula 15.Shown in program internal layer realization flow.
For sake of convenience, claim formula 15 to be model I, formula 16 adopts measured data shown in Fig. 5 to carry out identification to model I and model II for model II, and the simulation result of model I, model II as shown in Figure 6, Figure 7.
(1) model II has better fitting precision relative to model I.Think that main cause is as follows:
Model I take busbar voltage as excitation, by the power response of Solving Equation of State load.The change of load power demand when this model can reflect grid voltage change under different operating mode well, its power demand depends on electrical network busbar voltage completely.Measured power is obtained by transformer station's outlet Current calculation, is tractive power supply system actual power; And survey PCC point voltage except by except tractive power supply system loading effects, also depend on the change of other power loads.That is, both comprise electric locomotive impact characteristics in measured data, comprised again its part throttle characteristics, further comprises the characteristic of other loads a certain amount of.Therefore, model may be caused in this case to be absorbed in local excessive matching state around voltage disturbance, to cause error of fitting.
Model II adds the Analytical Expression of initiative demand power on the basis of model I in the power rush great moment, make model show better impact characteristics.On the one hand, this initiative power does not also rely on system voltage, and it is determined by the operating mode of locomotive operation, and it has effectively broken away from the Identification Errors that other load causes for the impact of system busbar voltage.On the other hand, initiative demand power add the yardstick expanded model and impact characteristics is described, make original model more be absorbed in small changed power and describe, be conducive to reducing value.Therefore, type II remains the descriptive power of original model for part throttle characteristics and impact characteristics while having better matching to impact characteristics.
(2) model I has better generalization ability relative to model II.
Electric railway impact load has very strong randomness, and the amplitude of power rush, time of origin and duration are different, do not possess strict cyclic variation feature.Model II cannot produce the characteristic point of initiative power rush by follow load in real time and accurately.
(3) model I be applicable to required precision not high when do load prediction and electric railway development plan; Model II can be used for revising the identification result of model I, is applicable to analyzing electric railway impact load to the impact of electrical network.
The part throttle characteristics not changing model impact load due to model II describes, and therefore only needs shown 7 parameters of identification table 2.Fitness function value result of calculation 0.0132.
Table 2 Model Independent identified parameters value

Claims (1)

1., based on an electrization railway draught load modeling method for parameter identification, it is characterized in that, comprise the following steps:
Step 1, in conjunction with the feature in electric railway load locomotive traveling process, directly assemble on contact line supply arm with direct current machine motor load and synthetic load, determine AC-DC electric locomotive Mathematical Modeling, directly assemble on contact line supply arm with three-phase alternating-current motor load and synthetic load, determine to hand over orthogonal type electric locomotive Mathematical Modeling, and AC-DC electric locomotive Mathematical Modeling is comprehensive with friendship orthogonal type electric locomotive Mathematical Modeling, set up electric railway part throttle characteristics model, in conjunction with the impact characteristics of electric railway traction supply load, determine the model of the impact characteristics of electric railway traction supply load, and insertion switch coefficient obtains by formula
The identification analytical expression model of the electric railway electric locomotive model determined, in formula: with for the switching coefficient of three-phase alternating-current motor load and DC motor load, for the meritorious and reactive power that three-phase alternating-current motor absorbs, for the active power that DC motor absorbs; for the meritorious and reactive power that synthetic load absorbs, the model of the impact characteristics of electric railway traction supply load;
Its AC-DC electric locomotive Mathematical Modeling is:
In formula: for the electric current of traction electric machine, for voltage transformation coefficient, for the initial active power proportionality coefficient of DC traction motor, for model is initially always gained merit and reactive power, for the indicial response of load, U is three-phase side voltage effective value, for steady state voltage;
Its friendship orthogonal type electric locomotive Mathematical Modeling is:
In formula: U is three-phase side voltage effective value, for steady state voltage, δ is merit angle, for the indicial response of load, for the initial active power of synthetic load, the initial reactive power of synthetic load, for transient state reactance, E ' is electromotive force after transient state reactance;
Electric railway part throttle characteristics model is:
P i=k 1P j+k 2P z+P s
Q i=k 1Q j+Q s
In formula: with for the switching coefficient of three-phase alternating-current motor load and DC motor load, for the active power that three-phase alternating-current motor absorbs, for the active power that DC motor absorbs, for the meritorious and reactive power that synthetic load absorbs, the reactive power that three-phase alternating-current motor absorbs;
The model of the impact characteristics of electric railway traction supply load is:
In formula: U is three-phase side voltage effective value, for steady state voltage, be respectively and gain merit and the variation coefficient of reactive power impact power with voltage;
Step 2, set up the transformation model between the three-phase side voltage of electric power system and two phase voltages of electric railway traction power supply system, in conjunction with electric railway stating with load, acceleration and the feature of advancing across electric phase-splitting, mathematical modeling is carried out to the power demand initiative feature of electric railway load, and in conjunction with the identification analytical expression model of electric railway electric locomotive model, set up by formula
The electric railway impact load model determined, in formula for the identification analytical expression model of electric railway electric locomotive model, for the power demand initiative feature Mathematical Modeling of electric railway load; for parameter to be identified;
The modeling method of the power demand initiative feature Mathematical Modeling of described electric railway load is; The power of locomotive active demand is equivalent to the changed power process of approximate trapezoid, the Mathematical Modeling of power initiative demand utilization step function represents, the model of foundation is by formula
Determine, in formula: be respectively parameter to be identified, for unit step function; with be respectively curve level in trapezoidal changed power figure, rise and power corresponding to sloping portion; for three superpose on a timeline after compressive stroke power;
Step 3, choose some electric railway traction stations and test, and analyze and preliminary treatment the measured data of electric railway traction supply load, its processing method comprises following steps:
Step 31, choosing some substation bus bars for power supply for locomotive station is measurement point, using one group that measures at this bus place typical and complete impact characteristics curve as analytic target, the change curve of voltage, electric current and apparent power in Measuring Time section,
step 32, adopts Wavelet Decomposition Technology to carry out filtering and denoising to measured data, completes the pretreatment work of measured data;
step 4, to the canonical parameter in electric railway traction supply load model, carries out parameter identification, and adopts corresponding example, carries out analysis verification to the accuracy of electric railway traction supply load model, needs the canonical parameter carrying out parameter identification to be: , L and R is respectively inductance and the resistance in traction electric machine loop, and J is the moment of inertia of traction electric machine for the initial active power proportionality coefficient of DC traction motor, for transient state reactance, X is synchronous reactance, for time constant, for the initial active power proportionality coefficient of three-phase alternating-current motor, be respectively and gain merit and the variation coefficient of reactive power impact power with voltage; When carrying out parameter identification to the canonical parameter in electric railway traction supply load model, according to Parameter Sensitivity Analysis result, the improvement PSO algorithm choosing linear decrease weights strategy (LDIW) carries out identification to typical independent parameter.
CN201410442604.1A 2014-09-02 2014-09-02 Electrified railway traction power supply load modeling method based on parameter identification Active CN104158189B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410442604.1A CN104158189B (en) 2014-09-02 2014-09-02 Electrified railway traction power supply load modeling method based on parameter identification

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410442604.1A CN104158189B (en) 2014-09-02 2014-09-02 Electrified railway traction power supply load modeling method based on parameter identification

Publications (2)

Publication Number Publication Date
CN104158189A CN104158189A (en) 2014-11-19
CN104158189B true CN104158189B (en) 2015-07-01

Family

ID=51883632

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410442604.1A Active CN104158189B (en) 2014-09-02 2014-09-02 Electrified railway traction power supply load modeling method based on parameter identification

Country Status (1)

Country Link
CN (1) CN104158189B (en)

Families Citing this family (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2605225C1 (en) * 2015-06-15 2016-12-20 Государственное бюджетное образовательное учреждение высшего образования Нижегородский государственный инженерно-экономический университет (НГИЭУ) Ac electric locomotive simulation device
CN105119278A (en) * 2015-08-26 2015-12-02 河海大学 Special load modeling method
CN106899014B (en) * 2015-12-18 2020-02-18 国网冀北电力有限公司张家口供电公司 Modeling prediction method of electrified railway load based on waveform decomposition
CN105938325B (en) * 2016-04-05 2021-10-29 中国电力科学研究院 System model identification method for AC/DC coordination control
CN106274996B (en) * 2016-08-29 2017-12-01 国网山西省电力公司电力科学研究院 A kind of method that passed through electric locomotive type is identified in traction substation
CN107526858B (en) * 2016-11-07 2023-12-22 北京交通大学 PSCAD/EMTDC-based simulation platform for electric iron traction power supply system
CN107527129B (en) * 2017-06-20 2023-07-18 国网浙江省电力公司温州供电公司 Power grid risk analysis method containing electric iron load
CN109687438B (en) * 2018-12-19 2022-06-24 西南交通大学 Power grid fragile line identification method considering high-speed rail impact load effect
CN109494733B (en) * 2018-12-21 2022-06-07 云南电网有限责任公司电力科学研究院 Identification parameter optimization method and system of power load model
CN113011004B (en) * 2021-02-24 2024-03-15 南京地铁建设有限责任公司 Security check method, device and storage medium
CN113159265B (en) * 2021-03-24 2022-09-09 国网河南省电力公司电力科学研究院 Traction load parameter identification method and system based on SVM-ant colony algorithm
CN113468696B (en) * 2021-08-05 2023-03-24 国核电力规划设计研究院重庆有限公司 Method for selecting section of lead of electric railway traction station
CN115329684B (en) * 2022-10-17 2023-01-24 广东电网有限责任公司中山供电局 Static load model construction method, system, equipment and medium

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101354423A (en) * 2008-09-16 2009-01-28 华北电力大学 System and method for model building of impact load based on actual measurement
CN102122310B (en) * 2011-02-01 2012-08-22 浙江大学 Train diagram-based traction load modeling method
CN102184296B (en) * 2011-05-13 2012-11-28 江西省电力科学研究院 Modelling method of impact load of electrified railway based on actually-measured data

Also Published As

Publication number Publication date
CN104158189A (en) 2014-11-19

Similar Documents

Publication Publication Date Title
CN104158189B (en) Electrified railway traction power supply load modeling method based on parameter identification
CN104201674B (en) Comprehensive load model modeling method considering load low voltage release features
CN102841251B (en) Electric car charging station harmonic wave detection method in microgrid
CN103618492B (en) A kind of Parameter Identification Method of Synchronous Generator based on time-frequency conversion
CN103036498B (en) Synchronous generator practical model parameter examination and identification method based on parameter measure unit (PMU)
CN104362622B (en) The load model of a kind of band small power station unit area power grid and modeling method
Zhang et al. Modeling of regenerative braking energy for electric multiple units passing long downhill section
CN102510263A (en) Method for identifying practical parameters of synchronous generator on basis of load rejection test and numerical difference
CN105628406A (en) High-speed train traction transmission control system fault injection method and simulation system
CN103401238A (en) Method for modeling power load based on measurement-based method
CN102122310A (en) Train diagram-based traction load modeling method
CN102998980A (en) Real-time simulation modeling method for tractive power supply system of electrified railway
Dharmakeerthi et al. Impact of electric vehicle load on power system oscillatory stability
CN107565867A (en) A kind of Parameter Identification Method of Synchronous Generator based on trace sensitivity
Huang et al. Evaluation scheme for EMI of train body voltage fluctuation on the BCU speed sensor measurement
CN102184296B (en) Modelling method of impact load of electrified railway based on actually-measured data
CN106777817B (en) Harmonic analysis modeling method of urban rail transit electric energy quality simulation analysis system
CN101917000B (en) Method for evaluating load of traction substation of electric railway
CN103604988B (en) Elevator feedback electric energy quality and feedback energy-saving efficiency online test method and system
CN1605878A (en) Feeder line singlephase fault and multiphase fault distance measuring method based on wavelet decomposition frequency band feature
CN102055197B (en) Method for establishing controllable serial compensation linearized model
CN105116251A (en) Transformer fault discrimination and protection method
CN107526858B (en) PSCAD/EMTDC-based simulation platform for electric iron traction power supply system
CN204044261U (en) A kind of electrical network resistance characteristic test device frequently
CN105260521A (en) Analytic hierarchy process based load fusion modeling method

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant