CN112350347B - Power control method, system, device and medium for rail transit vehicle-mounted power grid - Google Patents

Power control method, system, device and medium for rail transit vehicle-mounted power grid Download PDF

Info

Publication number
CN112350347B
CN112350347B CN202011394664.2A CN202011394664A CN112350347B CN 112350347 B CN112350347 B CN 112350347B CN 202011394664 A CN202011394664 A CN 202011394664A CN 112350347 B CN112350347 B CN 112350347B
Authority
CN
China
Prior art keywords
switch state
matrix converter
switch
alternating current
power
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
CN202011394664.2A
Other languages
Chinese (zh)
Other versions
CN112350347A (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.)
Zhengzhou Railway Vocational and Technical College
Original Assignee
Zhengzhou Railway Vocational and Technical College
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 Zhengzhou Railway Vocational and Technical College filed Critical Zhengzhou Railway Vocational and Technical College
Priority to CN202011394664.2A priority Critical patent/CN112350347B/en
Publication of CN112350347A publication Critical patent/CN112350347A/en
Application granted granted Critical
Publication of CN112350347B publication Critical patent/CN112350347B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/24Arrangements for preventing or reducing oscillations of power in networks
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/004Artificial life, i.e. computing arrangements simulating life
    • G06N3/006Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02MAPPARATUS FOR CONVERSION BETWEEN AC AND AC, BETWEEN AC AND DC, OR BETWEEN DC AND DC, AND FOR USE WITH MAINS OR SIMILAR POWER SUPPLY SYSTEMS; CONVERSION OF DC OR AC INPUT POWER INTO SURGE OUTPUT POWER; CONTROL OR REGULATION THEREOF
    • H02M5/00Conversion of ac power input into ac power output, e.g. for change of voltage, for change of frequency, for change of number of phases
    • H02M5/40Conversion of ac power input into ac power output, e.g. for change of voltage, for change of frequency, for change of number of phases with intermediate conversion into dc
    • H02M5/42Conversion of ac power input into ac power output, e.g. for change of voltage, for change of frequency, for change of number of phases with intermediate conversion into dc by static converters
    • H02M5/44Conversion of ac power input into ac power output, e.g. for change of voltage, for change of frequency, for change of number of phases with intermediate conversion into dc by static converters using discharge tubes or semiconductor devices to convert the intermediate dc into ac
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2310/00The network for supplying or distributing electric power characterised by its spatial reach or by the load
    • H02J2310/40The network being an on-board power network, i.e. within a vehicle
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02MAPPARATUS FOR CONVERSION BETWEEN AC AND AC, BETWEEN AC AND DC, OR BETWEEN DC AND DC, AND FOR USE WITH MAINS OR SIMILAR POWER SUPPLY SYSTEMS; CONVERSION OF DC OR AC INPUT POWER INTO SURGE OUTPUT POWER; CONTROL OR REGULATION THEREOF
    • H02M1/00Details of apparatus for conversion
    • H02M1/0003Details of control, feedback or regulation circuits
    • H02M1/0025Arrangements for modifying reference values, feedback values or error values in the control loop of a converter

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Power Engineering (AREA)
  • Theoretical Computer Science (AREA)
  • General Health & Medical Sciences (AREA)
  • Molecular Biology (AREA)
  • Biophysics (AREA)
  • Computational Linguistics (AREA)
  • Data Mining & Analysis (AREA)
  • Evolutionary Computation (AREA)
  • Artificial Intelligence (AREA)
  • Biomedical Technology (AREA)
  • Computing Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Software Systems (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Ac-Ac Conversion (AREA)

Abstract

A power control method, system, device and medium for a rail transit on-board power grid. The invention discloses a power control method of a rail transit vehicle-mounted power grid, which is characterized in that the power control of the conventional rail transit vehicle-mounted power grid is limited by the power control defect of an alternating current-alternating current matrix converter, the global optimal switch state combination of the matrix converter cannot be quickly obtained, and the stability control of the power of the vehicle-mounted power grid cannot be realized.

Description

Power control method, system, device and medium for rail transit vehicle-mounted power grid
Technical Field
The invention relates to the field of power control of a rail transit vehicle-mounted power grid.
Background
With the development of power electronic technology, the requirement for the stability of the vehicle-mounted power grid suitable for rail transit is increasingly urgent, the power supply of the vehicle-mounted power grid can be provided by a large power grid in China and also can be provided by an independent power supply, and the independent power supply can be an independent power supply of various power sources, such as an independent power supply taking solar energy, wind energy and photoelectric complementary forms as an energy source and an independent power supply based on diesel/gasoline engine power generation. The vehicle-mounted power grid of the rail transit needs to be matched with various power supplies to provide stable power for the vehicle, and the power control technology of the vehicle-mounted power grid becomes a research focus.
A matrix converter is needed in a rail transit vehicle-mounted power grid to achieve stable alternating current output, and constant frequency indexes of a power generation system are achieved. In the past decades, ac-ac matrix converters have been used due to their unique characteristics: the input power factor is controllable, the input current harmonic is small, the power density is high, the energy flows in two directions, and the like, and the device is greatly developed, and particularly applied to rail transit.
Since the ac-ac matrix converter has no dc energy storage element, the influence of external disturbances on the converter is very prominent. Therefore, there is a need for a control method, system, device and storage medium with simple principle and good performance, which can improve the noise immunity of the matrix converter to ensure its safe and stable operation.
Generally, the collected voltage and current and a discrete mathematical model of the converter are utilized to circularly calculate the predicted power error, and the optimal vector of the matrix converter is selected by finding the minimum predicted power error to realize the power control of the converter. Because the collected voltage and current information has time delay, the prediction result is influenced, the sampling period is shortened, the prediction precision is favorably improved, but the calculation capability of the chip is higher, therefore, the power prediction error value of the matrix converter is calculated by adopting a particle swarm optimization method in the prior art, although the calculation load of the chip is reduced by the method, the optimal solution cannot be found.
In the power control of the rail transit vehicle-mounted power grid, the ant colony algorithm is introduced into the calculation of the power prediction error of the matrix converter, so that the optimal solution can be found, meanwhile, at the beginning of the ant colony algorithm calculation, the distance value between the switch states is set according to the switch mutual exclusion relation in the alternating current-alternating current matrix converter, the global optimal switch state combination can be rapidly obtained, the obtained optimal switch state combination is applied to the matrix converter, and the vehicle-mounted power grid is ensured to provide stable alternating current. By the power control method, the system, the device and the medium for the rail transit vehicle-mounted power grid, provided by the invention, the global minimum predicted power error can be quickly found, the global optimal switch state combination is obtained, the quick and accurate power control of the alternating current-alternating current matrix converter is realized, the running stability of a rail vehicle which is connected with the output end of the alternating current-alternating current matrix converter and serves as a load is further ensured, and the vehicle-mounted power grid is realized to provide stable alternating current output.
Disclosure of Invention
The invention aims to provide a power control method, a system, a medium and a device for a vehicle-mounted power grid, and the method, the system, the medium and the device can provide stable alternating current for a vehicle and ensure the running stability of the vehicle.
The invention provides a power control method for a rail transit vehicle-mounted power grid, which is characterized by comprising the following steps of 1: initializing an ant colony algorithm according to the number of switches of the AC-AC matrix converter; step 2: determining the distance between the switch states according to the on-off rule of the switch in the AC-AC matrix converter, and placing ants in different switch states; and step 3: the global optimal switch state combination module calculates the predicted power error of the alternating current-alternating current matrix converter every time the ants complete traversal; and 4, step 4: after the maximum iteration times are finished, obtaining the optimal switching state combination of the current AC-AC matrix converter according to the obtained minimum global prediction power error of the AC-AC matrix converter; and 5: and (5) adjusting the alternating current-alternating current matrix converter according to the optimal switching state combination obtained in the step (4) to ensure that the vehicle-mounted power grid provides stable alternating current.
The step 1 specifically comprises the following steps: the number of ants in the ant colony is m, the number of switches in the AC-AC matrix converter is n, and each switch is represented by SW 1 ,SW 2 ,…,SW n Each switch has two states, on and off, expressed as SW 1-1 /SW 1-0 ,SW 2-1 /SW 2-0 ,…,SW n-1 /SW n-0 The ant needs to traverse n switches and 2n switch states, and finally finds the optimal solution of the on-off condition of the n switches.
M is 1.5 times of 2n, namely m =3n, and the distance from the switch state i to the switch state j is d ij (i, j =1,2, \8230;, 2 n), in particular, the distance between the different states of one and the same switch is set to a maximum value, thereby expressing that one and the same switch can only have one state, either on or off, preferably the distance between the different states of one and the same switch is set to infinity.
the concentration of pheromone on a connecting path between the switch state i and the switch state j at the time tc is tau ij (tc). At the initial time, the pheromone concentration on each switch state connecting path is the same and is tau ij (tc)=τ 0 . Then the ants will select the circuit according to a certain probability,
Figure BDA0002814280920000021
representing the probability of ant q transitioning from switch state i to switch state j at time tc.
Figure BDA0002814280920000022
Wherein,
Figure BDA0002814280920000031
representing the desired degree of ant transfer from switch state i to switch state j for a heuristic function,
Figure BDA0002814280920000032
allow com for ants to access the set of switch states, if switch state j belongs to allow com Then, then
Figure BDA0002814280920000033
If the switch state j does not belong to all com Then->
Figure BDA0002814280920000034
Mu represents a pheromone importance factor, and a larger value indicates that the pheromone concentration has a larger effect in metastasis; theta is an important degree factor of the heuristic function, and the larger the value of theta indicates that the heuristic function has larger effect in transfer, ants can be transferred to a short-distance switch state with higher probability;
the step 2 specifically comprises the following steps: and calculating the distance value between the switch states according to the mutual exclusion relation of on-off of the switches.
First, d ij The calculation method of (A) is as follows:
Figure BDA0002814280920000035
Figure BDA0002814280920000036
is all com Represents a set of switch states that have been accessed, γ is a constant;
then, in order to further improve the calculation efficiency and accuracy of the ant colony algorithm, d is calculated ij Making a correction when j is equal to
Figure BDA0002814280920000037
When any switch state is not the same, judging whether the switches respectively corresponding to the switch state i and the switch state j are in a mutual exclusion relationship:
if the two are mutually exclusive, continuously judging whether the current switch state i and the switch state j have the same state indication, namely, both are on or both are off, if the indication states are the same, d ij = ∞; if the indication states are different, it indicates that the switch state i and the switch state j meet the mutual exclusion requirement, and if one switch is on, the two states are necessarily simultaneously existed, so that the distance between the switch state i and the switch state j at the moment is set as the minimum value, which can be 1.
If the switches corresponding to the switch states i and j are not in an exclusive relationship, the distance between the switch states i and j is γ, and γ is a constant greater than 1.
Thus, d ij The final calculation method of (c) is:
Figure BDA0002814280920000041
ants were randomly placed on the switch state.
Initially, allow com 2n-1 elements in the tree, representing all but the ant q-initiated switch state, all over time com The number of elements in (1) is continuously reduced until the number is null, which indicates that all switch states are completely accessed. When the ants release the pheromone, the pheromone on each switch state connecting path gradually disappears, and after all the ants finish one-time circulation, the concentration of the pheromone on each switch state connecting path is updated in real time, which specifically comprises the following steps:
Figure BDA0002814280920000042
in the formula,
Figure BDA0002814280920000043
releasing pheromone concentration for the qth ant on a connection path between a switch state i and a switch state j; delta tau ij The sum of the concentration of pheromones released by all ants on the connection path of the switch state i and the switch state j; the parameter epsilon represents the volatilization degree of the pheromone, and epsilon is more than 0 and less than 1; />
Figure BDA0002814280920000044
Wherein Q is a constant, typically having a value of [20, 2000%]L represents the total amount of pheromone released by ants in one cycle q The length of the path passed by the qth ant is represented by the number of nodes.
The step 3 specifically comprises the following steps:
step 3.1: collecting corresponding voltage and current at the input end of the matrix converter, and transmitting the corresponding voltage and current to the global optimal switch state combination module; step 3.2: the global optimal switching state combination module calculates the next beat of the kth beat sampling point, namely the output power of the kth +1 beat matrix converter according to the power prediction model of the alternating current-alternating current matrix converter; step 3.3: transmitting the voltage and the current of the k +1 th beat output end of the AC-AC matrix converter to a global optimal switch state combination module; step 3.4: the global optimal switch state combination module calculates the actual power of the output end of the alternating current-alternating current matrix converter; step 3.5: and the global optimal switch state combination module calculates the predicted power error of the alternating current-alternating current matrix converter according to the calculation results of the step 3.2 and the step 3.4.
Specifically, the predicted power calculation formula is:
Figure BDA0002814280920000051
wherein it is present>
Figure BDA0002814280920000052
And &>
Figure BDA0002814280920000053
Is the predicted value of the current at the kth +1 th beat at the lower input end of the alpha _ beta coordinate system, and is judged according to the predicted value>
Figure BDA0002814280920000054
And &>
Figure BDA0002814280920000055
Is the k-th beat sampling voltage value of the input end under the alpha-beta coordinate system. The k +1 th beat power at the output end of the AC-AC matrix converter is calculated by the formula
Figure BDA0002814280920000056
Wherein it is present>
Figure BDA0002814280920000057
And &>
Figure BDA0002814280920000058
Is the voltage value of the (k + 1) th beat at the input end under an alpha _ beta coordinate system, and is combined with the circuit>
Figure BDA0002814280920000059
And &>
Figure BDA00028142809200000510
Is the current value of the (k + 1) th beat of the input end under the alpha _ beta coordinate system.
When the AC-AC matrix converter takes a beat of k +1, the predicted value of the output power of the matrix converter at the input end is subtracted from the output power of the matrix converter and then squared to obtain the power prediction error variance err of the AC-AC matrix converter, namely
Figure BDA00028142809200000511
The step 4 specifically comprises the following steps: and setting the objective function as min (err), namely solving the minimum value of the predicted power error, and selecting the optimal switch state combination according to the obtained minimum value of the predicted power error.
The step 5 specifically comprises the following steps: and transmitting the optimal switch state combination to an alternating current-alternating current matrix converter to control the normal operation of the alternating current-alternating current matrix converter, so as to ensure stable alternating current of a rail transit vehicle-mounted power grid.
The power control method for the rail transit vehicle-mounted power grid provided by the first aspect of the invention can ensure that the vehicle-mounted power grid provides stable alternating current to rail vehicles and ensures the running stability of the rail vehicles.
The device can realize the power control method for the rail transit vehicle-mounted power grid provided by the first aspect of the invention, ensure that the vehicle-mounted network provides stable alternating current for rail vehicles and ensure the running stability of the rail vehicles.
Another aspect of the present invention is to provide a computer for power control of a rail transit vehicle-mounted power grid, where the program may implement the power control method for a rail transit vehicle-mounted power grid provided by the first aspect of the present invention, so as to ensure that a vehicle-mounted power grid provides a stable alternating current to a rail vehicle, and ensure the stability of operation of the rail vehicle.
Another aspect of the present invention is to provide a storage medium, where a computer program for executing the power control method for a rail transit vehicle electrical system is stored, so as to ensure that a vehicle electrical system provides stable alternating current to a rail vehicle, and ensure the stability of operation of the rail vehicle.
Has the advantages that: in the power control of the rail transit vehicle-mounted power grid, the improved ant colony algorithm is applied to the control of the matrix converter, the distance factor between the switch states is determined according to the mutual exclusion relation between the switches in the AC-AC matrix converter, then the global optimal solution can be obtained in a short time, the optimal switch state combination of the AC-AC matrix converter is found, the accurate control of the power prediction of the matrix converter is realized, the power control accuracy of the rail transit vehicle-mounted power grid is finally realized, the vehicle-mounted power grid is ensured to provide stable AC, and the running stability of rail vehicles is ensured.
Drawings
FIG. 1: power control structure diagram of rail transit vehicle-mounted power grid
FIG. 2 is a drawing: topological diagram of two-stage AC-AC matrix converter
FIG. 3 is a drawing: computing flow chart for obtaining optimal switch state combination of alternating current-alternating current matrix converter based on ant colony algorithm
Reference numerals
U a 、U b 、U c Is a three input voltage, L f Representing input terminal equivalent inductance, C f Representing input equivalent electricityC, holding; l represents the equivalent inductance of the output end of the AC-AC matrix converter, and R represents the equivalent capacitance of the output end of the AC-AC matrix converter; s ap 、S an 、S bp 、S bn 、S cp 、S cn ,S Ap 、S An 、S Bp 、S Bn 、…、S Xp 、S Xn Switches, S, both of AC-AC matrix converters ap 、S an 、S bp 、S bn 、S cp 、S cn ,S Ap Are composed of two IGBTs with reverse parallel diodes connected according to a common emitter stage, S Ap 、S An 、S Bp 、S Bn 、…、S Xp 、S Xn Each consisting of one IGBT.
Detailed Description
The present invention is further illustrated by the following examples.
At present, the main rail transportation modes comprise high-speed railways, motor train units, heavy-duty freight locomotives, light rails, subway vehicles, monorail trains and magnetic suspension rail transportation, wherein long-distance and high-power main rail transportation generally adopts an alternating current power supply system. The vehicle-mounted power grid of the rail transit needs to be matched with various power supplies to provide stable power for vehicles, the power supply of the vehicle-mounted power grid can be provided by a large power grid in China, and can also be provided by an independent power supply, and the independent power supply can be an independent power supply of various power sources, such as an independent power supply taking a solar energy, wind energy and photoelectric complementary form as an energy source, and an independent power supply based on diesel/gasoline engine power generation. The on-board electrical system requires power control in order to provide a stable alternating current for the rail vehicle.
The rail vehicle obtains energy from a vehicle-mounted power grid, an alternating current-alternating current matrix converter is needed to ensure that the energy is stable and controllable, and if the matrix converter wants to realize power control more quickly and accurately, an optimal combination mode of internal switch states of the matrix converter needs to be obtained.
According to the invention, as shown in figure 1, in a vehicle-mounted power grid, voltage and current values of an input end in sampling beat are collected and transmitted to a global optimal switch state combination module, and the module calculates the predicted power of an output end in the next beat of the sampling beat; meanwhile, the voltage and the current of the next beat of output end of the AC-AC matrix converter are transmitted to a global optimal switch state combination module, and the module calculates the actual power of the output end; the global optimal switch state combination module calculates a power error based on predicted power and actual power, obtains an optimal switch state combination of the alternating current-alternating current matrix converter by solving a minimum power error, adjusts the setting of the matrix converter by using the combination, realizes accurate control of the matrix converter power prediction, further realizes the stability of energy supply of the rail vehicle, and ensures the running stability of the rail vehicle.
Taking T as a discrete period, based on an instantaneous power calculation mode, the discrete mathematical model of the AC-AC matrix converter is as follows:
Figure BDA0002814280920000071
u in 、i in respectively representing the sampled input voltage and current; u. u out 、i out Respectively representing the output voltage and current; l is an output end equivalent inductor, R is an output end equivalent resistor, and C is an output end equivalent capacitor; t is a discrete period, T is a time constant, and k represents the kth beat of samples.
Since the sampling frequency is very high, it can be considered that the voltages at the input terminals of two consecutive times are approximately equal,
Figure BDA0002814280920000072
the predicted power calculation formula is as follows: />
Figure BDA0002814280920000073
Wherein it is present>
Figure BDA0002814280920000074
And &>
Figure BDA0002814280920000075
Is the predicted value of the current at the kth +1 th beat at the lower input end of the alpha _ beta coordinate system, and is judged according to the predicted value>
Figure BDA0002814280920000076
And &>
Figure BDA0002814280920000077
Is the k-th beat sampling voltage value of the input end under the alpha-beta coordinate system. The calculation formula of the k +1 th beat power of the output end of the AC-AC matrix converter is->
Figure BDA0002814280920000078
Wherein +>
Figure BDA0002814280920000079
And &>
Figure BDA00028142809200000710
Is the (k + 1) th beat voltage value at the input end under the alpha _ beta coordinate system, is greater than or equal to>
Figure BDA00028142809200000711
And &>
Figure BDA00028142809200000712
Is the current value of the (k + 1) th beat of the input end under the alpha _ beta coordinate system.
The predicted value of the output power of the matrix converter from the k +1 th beat input end of the AC-AC matrix converter is subtracted from the output power of the matrix converter and then squared to obtain the power prediction error variance err of the AC-AC matrix converter, namely
Figure BDA00028142809200000713
And setting the objective function as min (err), namely solving the minimum error value of the predicted power, obtaining the optimal switch state combination according to the minimum error value, and finally transmitting the optimal switch state combination to the AC-AC matrix converter to control the normal operation of the AC-AC matrix converter.
FIG. 2 shows a topology of a two-stage AC-AC matrix converter, U a 、U b 、U c Is a three input voltage, L f Represents the input terminal equivalent inductance, C f Representing the input terminal equivalent capacitance; l represents AC-AC matrix changeThe equivalent inductance of the output end of the converter, R represents the equivalent capacitance of the output end of the AC-AC matrix converter; s ap 、S an 、S bp 、S bn 、S cp 、S cn ,S Ap 、S An 、S Bp 、S Bn 、…、S Xp 、S Xn Switches, S, both of AC-AC matrix converters ap 、S an 、S bp 、S bn 、S cp 、S cn ,S Ap Composed of two IGBTs with reverse parallel diodes connected in a common emitter stage, S Ap 、S An 、S Bp 、S Bn 、…、S Xp 、S Xn Each consisting of one IGBT. Fig. 2 a-c matrix converter has a common combination of switching states (2) 6 +2 2X ) And (2) carrying out switching state, wherein X represents the last Xth switching group of the inverter stage, and total number of (2X + 6) switches. The switch state combinations are shown in the following table:
serial number S ap S an S bp S bn S cp S cn S Ap S An S Bp S Bn S Xp S Xn
1 0 0 0 0 1 1 0 0 0 0 0 1
2 0 0 1 1 0 0 0 0 1 0 0 0
3 1 1 0 0 0 0 0 1 0 0 1 0
4 1 0 0 1 0 0 0 1 1 1 0 0
5 1 0 0 0 0 1 1 0 0 1 1 0
6 0 0 1 0 0 1 1 0 1 1 1 1
7 0 1 0 0 1 0 1 1 1 0 0 1
8 0 1 1 0 0 0 0 0 0 0 0 0
9 0 0 0 1 1 0 0 0 1 0 1 0
10 0 0 0 0 0 1 0 1 0 1 0 0
0 1 1 1 1 0
2 6 -2 2X -1
In order to quickly obtain a global optimal solution of an objective function min (err), an ant colony algorithm is introduced into the calculation of a predicted power minimum error of the AC-AC matrix converter, and in order to better obtain the global optimal solution, the distance between the switch states is set according to a mutual exclusion relation between the switches, and if the switches have a mutual exclusion requirement, namely two switches are inevitably kept in an open (conducting) state and a closed (disconnecting) state when the AC-AC matrix converter is actually applied, the distance value between the switch states needs to be determined according to the switch states.
For example, the two-stage AC-AC matrix converter of FIG. 2 of the present invention may have S at the inverter stage Ap 、S An Are a group S Bp 、S Bn The control is performed for one group, and the driving pulse meets the principle that one group is turned on and the other group is turned off. Meanwhile, the on-off of the rectifier stage and the inverter stage switching tube also needs to ensure the synchronism.
Fig. 3 shows a calculation flow chart of the ant colony algorithm-based ac-ac matrix converter for obtaining the optimal switch state combination. Firstly, initializing an ant colony algorithm; then, determining the distance between the switch states according to the switch on-off rule; randomly placing ants to any on-off state; selecting the switching state of the next alternating current-alternating current matrix converter for each ant, and calculating the concentration of the pheromone; then judging whether all the switch states are traversed, if the traversal is finished, updating the pheromone table, and if not, continuously selecting the next switch state until all the switch states are traversed; and then, judging whether the maximum iteration times of the ant colony algorithm are finished or not, if so, obtaining a global optimal switch combination state, otherwise, executing the ant colony algorithm again until the maximum iteration times are finished.
The number of ants in the ant colony is m, the number of switches in the AC-AC matrix converter is n, and each switch can be represented as SW 1 ,SW 2 ,…,SW n Each switch has two states, on and off, expressed as SW 1-1 /SW 1-0 ,SW 2-1 /SW 2-0 ,…,SW n-1 /SW n-0 The ant needs to traverse n switches and 2n switch states, and finally finds the optimal solution of the on-off conditions of the n switches. According to the experimental result, m is set to be 1.5 times of 2n, namely m =3n, so that the global optimal solution can be obtained quickly, and local improper convergence is prevented. The distance between the switch state i and the switch state j is set as d ij (i, j =1,2, \8230;, 2 n), in particular, the distance between the different states of one and the same switch is set to a maximum value, thereby expressing that one and the same switch can only have one state, either on or off, preferably the distance between the different states of one and the same switch is set to infinity. the concentration of pheromone on the connection path between the switch state i and the switch state j at the time of tc is tau ij (tc). At the initial time, the pheromone concentration on each switch state connecting path is the same and is tau ij (tc) = τ 0. Then the ants will select the circuit according to a certain probability,
Figure BDA0002814280920000091
representing the probability of ant q transitioning from switch state i to switch state j at time tc.
Figure BDA0002814280920000092
Wherein,
Figure BDA0002814280920000093
representing the desired degree of ant transfer from switch state i to switch state j, as a heuristic function,
Figure BDA0002814280920000094
allow com for ants to access the set of switch states, if switch state j belongs to allow com Then, then
Figure BDA0002814280920000095
If the switch state j does not belong to all com Then->
Figure BDA0002814280920000096
Mu means letterThe larger the value of the pheromone importance factor, the larger the effect of the pheromone concentration in the transfer is; theta is an important degree factor of the heuristic function, the larger the value of theta indicates that the heuristic function has larger effect in transfer, and ants can be transferred to a switch state with short distance with higher probability;
first, d ij The calculation method is as follows:
Figure BDA0002814280920000097
Figure BDA0002814280920000098
is all com Represents a set of switch states that have been accessed, γ is a constant;
then, in order to further improve the calculation efficiency and accuracy of the ant colony algorithm, d is calculated ij Making a correction when j and
Figure BDA0002814280920000101
when any switch state is not the same, judging whether the switches respectively corresponding to the switch state i and the switch state j are in a mutual exclusion relationship:
if the two are mutually exclusive, continuously judging whether the current switch state i and the switch state j have the same state indication, namely, the two states are on or off, and if the indication states are the same, d ij = ∞; if the indication states are different, it indicates that the switch state i and the switch state j meet the mutual exclusion requirement, and if one switch is on, the two states are necessarily simultaneously existed, so that the distance between the switch state i and the switch state j at the moment is set as the minimum value, which can be 1.
If the switches corresponding to the switch states i and j are not mutually exclusive, the distance between the switch states i and j is gamma, and gamma is a constant larger than 1.
Thus, d ij The final calculation method of (a) is as follows:
Figure BDA0002814280920000102
ants were placed on different switch states.
Initially, allow com There are 2n-1 elements, i.e. all other switch states except the one from ant q, all over time com The elements in (b) are continuously reduced until the elements are empty, which indicates that all switch states are completely accessed. When the ants release the pheromone, the pheromone on each switch state connecting path gradually disappears, and after all the ants finish one-time circulation, the concentration of the pheromone on each switch state connecting path is updated in real time, which specifically comprises the following steps:
Figure BDA0002814280920000103
in the formula,
Figure BDA0002814280920000111
releasing pheromone concentration on a connection path of a switch state i and a switch state j for the qth ant; delta tau ij The sum of the concentration of released pheromones of all ants on the connection path of the switch state i and the switch state j; the parameter epsilon represents the volatilization degree of the pheromone, and epsilon is more than 0 and less than 1; />
Figure BDA0002814280920000112
Wherein Q is a constant, typically having a value of [20, 2000%]L represents the total amount of pheromone released by ants in one cycle q The path length of the q-th ant is expressed by the number of nodes.
The maximum iteration number of the whole ant colony algorithm is set to iter _ max, and the value range of the maximum iteration number is usually [100, 600], so that premature convergence is not easy to occur, and a globally optimal solution can be found, preferably, the value of iter _ max is 220.
While one embodiment of the present invention has been described in detail, the description is only a preferred embodiment of the present invention and should not be taken as limiting the scope of the invention. All equivalent changes and modifications made within the scope of the present invention shall fall within the scope of the present invention.

Claims (6)

1. A power control method for a rail transit vehicle-mounted power grid is characterized in that,
step 1: initializing an ant colony algorithm according to the number of switches of the AC-AC matrix converter;
step 2: determining the distance between the switch states according to the on-off rule of the switch in the AC-AC matrix converter, and placing ants on different switch states;
and step 3: the global optimal switch state combination module calculates the predicted power error of the alternating current-alternating current matrix converter every time the ants complete traversal;
and 4, step 4: after the maximum iteration times are finished, obtaining the optimal switching state combination of the current AC-AC matrix converter according to the obtained minimum global prediction power error of the AC-AC matrix converter;
and 5: adjusting the alternating current-alternating current matrix converter according to the optimal switching state combination obtained in the step 4 to ensure that the vehicle-mounted power grid provides stable alternating current;
the step 1 specifically comprises the following steps: the number of ants in the ant colony is m, the number of switches in the AC-AC matrix converter is n, and each switch is represented by SW 1 ,SW 2 ,…,SW n Each switch has two states, on and off, expressed as SW 1-1 /SW 1-0 ,SW 2-1 /SW 2-0 ,…,SW n-1 /SW n-0 The ant needs to traverse n switches and 2n switch states, and finally finds the optimal solution of the on-off conditions of the n switches;
the concentration of pheromone on a connecting path between the switch state i and the switch state j at the time tc is tau ij (tc); at the initial time, the pheromone concentration on each switch state connecting path is the same and is tau ij (tc)=τ 0 (ii) a Then the ants will select the circuit according to a certain probability,
Figure FDA0004045651120000011
representing the probability of the ant q transferring from the switch state i to the switch state j at the tc moment;
Figure FDA0004045651120000012
wherein,
Figure FDA0004045651120000013
representing the desired degree of ant transfer from switch state i to switch state j for a heuristic function,
Figure FDA0004045651120000014
allow com for ants to access the set of switch states, if switch state j belongs to allow com Then, then
Figure FDA0004045651120000015
If the switch state j does not belong to all com Then->
Figure FDA0004045651120000016
Mu represents a pheromone importance factor, and a larger value indicates that the pheromone concentration has a larger effect in metastasis; theta is a factor of importance degree of the heuristic function, and the larger the value of theta is, the larger the effect of the heuristic function in transfer is; the distance between the switch state i and the switch state j is set as d ij (i,j=1,2,…,2n);
The step 2 specifically comprises the following steps: d ij The calculation method of (A) is as follows:
Figure FDA0004045651120000021
Figure FDA0004045651120000022
is all com Represents a set of switch states that have been accessed, gamma is a constant greater than 1(ii) a When j and->
Figure FDA0004045651120000023
When any switch state is not the same, judging whether the switches corresponding to the switch state i and the switch state j are in mutual exclusion relation, if not, judging d ij = γ; if the two are mutually exclusive, continuously judging whether the switch state i and the switch state j have the same state indication, if the two are both in an open state or in a closed state, d ij = ∞ if both states indicate that the mutual exclusion requirement is met, i.e. one is open and one is closed, then d ij =1; placing ants on different switch states;
the step 3 specifically comprises the following steps: step 3.1: collecting corresponding voltage and current at the input end of the matrix converter, and transmitting the corresponding voltage and current to the global optimal switch state combination module; step 3.2: the global optimal switch state combination module calculates the next beat of the kth beat sampling point, namely the output power of the kth +1 beat matrix converter according to the power prediction model of the alternating current-alternating current matrix converter; step 3.3: transmitting the voltage and the current of the k +1 th beat output end of the AC-AC matrix converter to a global optimal switch state combination module; step 3.4: the global optimal switch state combination module calculates the actual power of the output end of the alternating current-alternating current matrix converter; step 3.5: the global optimal switch state combination module calculates the predicted power error of the alternating current-alternating current matrix converter according to the calculation results of the step 3.2 and the step 3.4;
step 3 also includes: the predicted power calculation formula is:
Figure FDA0004045651120000024
wherein +>
Figure FDA0004045651120000025
And &>
Figure FDA0004045651120000026
Is the current of the (k + 1) th beat of the input end under the alpha _ beta coordinate systemA predictor value +>
Figure FDA0004045651120000027
And &>
Figure FDA0004045651120000028
Is the kth beat sampling voltage value of the input end under the alpha _ beta coordinate system; the calculation formula of the k +1 th beat power at the output end of the AC-AC matrix converter is->
Figure FDA0004045651120000029
Wherein it is present>
Figure FDA00040456511200000210
And
Figure FDA00040456511200000211
is the voltage value of the (k + 1) th beat at the input end under an alpha _ beta coordinate system, and is combined with the circuit>
Figure FDA00040456511200000212
And &>
Figure FDA00040456511200000213
The current value of the (k + 1) th beat of the input end under an alpha _ beta coordinate system; when the alternating current-alternating current matrix converter beats at the k +1 th time, the input end calculates the difference between the predicted value of the output power of the matrix converter and then squares the difference to obtain the power prediction error variance err, namely ^ er>
Figure FDA0004045651120000031
2. The method of claim 1, further characterized by initially allowing com 2n-1 elements in the ant represent all the switch states except the switch state from which the ant q starts; all over time com The element in (1) is continuously reduced until the element is empty, which indicates the positionThe access is finished in all the switch states; when the ants release the pheromone, the pheromone on each switch state connecting path gradually disappears, and after all the ants finish one-time circulation, the concentration of the pheromone on each switch state connecting path is updated in real time, which specifically comprises the following steps:
Figure FDA0004045651120000032
in the formula,
Figure FDA0004045651120000033
releasing pheromone concentration for the qth ant on a connection path between a switch state i and a switch state j; delta tau ij The sum of the concentration of pheromones released by all ants on the connection path of the switch state i and the switch state j; the parameter epsilon represents the volatilization degree of the pheromone, and epsilon is more than 0 and less than 1;
Figure FDA0004045651120000034
wherein Q is a constant having a value of [20, 2000%]L represents the total amount of pheromone released by ants in one cycle q The length of the path passed by the qth ant is represented by the number of nodes. />
3. The method according to claim 1, wherein step 4 is specifically: setting an objective function as min (err), wherein err represents the power prediction error variance of the AC-AC matrix converter, solving the minimum value of the prediction power error, and selecting the optimal switch state combination according to the obtained minimum value of the prediction power error.
4. A power control system for a rail transit onboard electrical network, which system is capable of implementing the power control method for a rail transit onboard electrical network according to any of claims 1 to 3.
5. A power device for a rail transit vehicle-mounted power grid, which comprises an input parameter module, an alternating current-alternating current matrix converter and a global optimal switch state combination module, and can realize the power control method for the rail transit vehicle-mounted power grid according to any one of claims 1 to 3.
6. A storage medium, on which a computer program is stored, which is capable of implementing the power control method for a rail transit on-board electrical network of any of claims 1-3.
CN202011394664.2A 2020-12-03 2020-12-03 Power control method, system, device and medium for rail transit vehicle-mounted power grid Active CN112350347B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011394664.2A CN112350347B (en) 2020-12-03 2020-12-03 Power control method, system, device and medium for rail transit vehicle-mounted power grid

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011394664.2A CN112350347B (en) 2020-12-03 2020-12-03 Power control method, system, device and medium for rail transit vehicle-mounted power grid

Publications (2)

Publication Number Publication Date
CN112350347A CN112350347A (en) 2021-02-09
CN112350347B true CN112350347B (en) 2023-04-07

Family

ID=74427219

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011394664.2A Active CN112350347B (en) 2020-12-03 2020-12-03 Power control method, system, device and medium for rail transit vehicle-mounted power grid

Country Status (1)

Country Link
CN (1) CN112350347B (en)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107612409A (en) * 2017-09-27 2018-01-19 西安理工大学 Simplify matrix converter model predictive control method with magnetic saturation control
CN109639153A (en) * 2018-12-26 2019-04-16 北京航空航天大学 A kind of model predictive control method of the source Quasi-Z indirect matrix converter
CN109638876A (en) * 2019-03-07 2019-04-16 国网辽宁省电力有限公司大连供电公司 A kind of unit starting sequence optimizing method based on ant group algorithm

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2557294B (en) * 2016-12-05 2022-03-30 Itt Mfg Enterprises Llc Matrix converter control method and system

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107612409A (en) * 2017-09-27 2018-01-19 西安理工大学 Simplify matrix converter model predictive control method with magnetic saturation control
CN109639153A (en) * 2018-12-26 2019-04-16 北京航空航天大学 A kind of model predictive control method of the source Quasi-Z indirect matrix converter
CN109638876A (en) * 2019-03-07 2019-04-16 国网辽宁省电力有限公司大连供电公司 A kind of unit starting sequence optimizing method based on ant group algorithm

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
CRH5动车组牵引系统的多输出间接矩阵变换器的控制研究;王闪闪等;《科技创新与应用》;20201020(第31期);第9-15页 *

Also Published As

Publication number Publication date
CN112350347A (en) 2021-02-09

Similar Documents

Publication Publication Date Title
CN108964469B (en) Full-bridge double LLC resonant converter with parallel-series structure
CN110768531A (en) Bidirectional high-frequency isolation type DC/DC module
CN109149922B (en) Power factor correction circuit and alternating current charger for electric automobile using same
CN112350347B (en) Power control method, system, device and medium for rail transit vehicle-mounted power grid
Brovont et al. DM and CM modeling of non-isolated buck converters for EMI filter design
Lima et al. A bidirectional isolated integrated AC–DC converter based on an interleaved 3-Level T-type power converters
CN105226986A (en) A kind of inverter and control method thereof eliminating the pulsation of input side secondary power
Farrakhov et al. Quasi-resonant buck converter for high power application
Kopacz et al. Three-Level Interleaved Non-isolated DC/DC Converter as a Battery Interface in an EV Charging System with Bipolar DC-Link.
Hu et al. A transformer‐less boost DC‐DC converter with high gain and continuous input current for fuel cell vehicle
CN210536515U (en) Fuel cell automobile power system simulation platform
Safayatullah et al. A Three-Port DC-DC-DC Converter based on Dual Active Bridge Series Resonant Topology for Electric Vehicle DC Fast Charging Applications
Hossain et al. Power stage design of a synchronous buck converter for battery charger application
CN112821788B (en) Hybrid charging control device and control method
Zhang et al. Analysis of an LLC–DAB integrated DC–DC converter topology with low loss and low cost
Wang et al. Theoretical comparison, Real-time emulation and Experiment validation of DC/DC Converter for Fuel Cell Electric Vehicle
Karami et al. Optimal Control of Triple Active Bridge Based on Deep Machine Learning Techniques
Ramanath et al. Equivalent Modeling, Design and Analysis of Integrated Magnetics Ćuk Converter
Dinis et al. Maximizing energy transfer in wireless power transfer systems using maximum power point tracking for in-motion ev and phev charging
CN112350586B (en) Ant colony algorithm-based alternating current-alternating current matrix converter power control method, system, medium and device
Lu et al. An integrated charging equalizer based on LLC resonant DC-DC converter for series-connected battery string in plug-in hybrid electric vehicle
Wang et al. Impact of boost converter switching frequency on optimal operation of fuel cell systems
Lin Implementation of a soft switching converter with series DC–DC circuits and single transformer
CN221328827U (en) Control circuit for reducing current distortion by utilizing voltage phase lag
Mishra et al. A 10-kW Active Front End Rectifier Fed Dual Active Bridge Converter For EV Charging

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