CN113904347B - Parameter optimization method and device for controllable phase shifter additional damping controller - Google Patents
Parameter optimization method and device for controllable phase shifter additional damping controller Download PDFInfo
- Publication number
- CN113904347B CN113904347B CN202111168548.3A CN202111168548A CN113904347B CN 113904347 B CN113904347 B CN 113904347B CN 202111168548 A CN202111168548 A CN 202111168548A CN 113904347 B CN113904347 B CN 113904347B
- Authority
- CN
- China
- Prior art keywords
- objective function
- function value
- universe
- phase shifter
- optimal objective
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000013016 damping Methods 0.000 title claims abstract description 87
- 238000005457 optimization Methods 0.000 title claims abstract description 58
- 238000000034 method Methods 0.000 title claims abstract description 38
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 68
- 230000006870 function Effects 0.000 claims description 138
- 230000008569 process Effects 0.000 claims description 12
- 238000012545 processing Methods 0.000 claims description 12
- 238000004590 computer program Methods 0.000 claims description 8
- 230000007246 mechanism Effects 0.000 claims description 8
- 238000004088 simulation Methods 0.000 claims description 8
- 238000012163 sequencing technique Methods 0.000 claims description 7
- 238000003860 storage Methods 0.000 claims description 5
- 230000010355 oscillation Effects 0.000 abstract description 11
- 230000001629 suppression Effects 0.000 abstract description 4
- 230000005540 biological transmission Effects 0.000 description 11
- 101100499229 Mus musculus Dhrsx gene Proteins 0.000 description 8
- JJWKPURADFRFRB-UHFFFAOYSA-N carbonyl sulfide Chemical compound O=C=S JJWKPURADFRFRB-UHFFFAOYSA-N 0.000 description 8
- 238000010586 diagram Methods 0.000 description 6
- 230000000903 blocking effect Effects 0.000 description 4
- 230000003993 interaction Effects 0.000 description 4
- 230000010363 phase shift Effects 0.000 description 4
- 230000001052 transient effect Effects 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 230000001939 inductive effect Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000012216 screening Methods 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/24—Arrangements for preventing or reducing oscillations of power in networks
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Feedback Control In General (AREA)
Abstract
The invention discloses a parameter optimization method and a device of an additional damping controller of a controllable phase shifter, which are implemented by establishing a power system model comprising the additional damping controller of the controllable phase shifter; acquiring and according to control parameters of an additional damper controller, automatically optimizing parameters under set performance requirements based on a multi-universe optimization algorithm, initializing the multi-universe optimization algorithm, and generating a preset number of universe individuals, wherein each universe individual corresponds to the control parameters with different values; acquiring an optimal objective function value after iteration according to a preset iteration number by calculating a first optimal objective function value in a universe individual; and updating the control parameters of the additional damping controller according to the optimal objective function value. Compared with the prior art, the invention optimizes the control parameters of the controllable phase shifter additional damping controller based on the multi-universe optimization algorithm, thereby realizing the suppression of the oscillation of the power system and improving the stability of the power system.
Description
Technical Field
The invention relates to the technical field of power systems, in particular to a parameter optimization method and device for an additional damping controller of a controllable phase shifter.
Background
Along with the development of national economic construction, the power consumption load is continuously increased, the load rate of the power transmission line is higher and higher, and the power transmission line is limited by environmental factors, economic factors and the like, including newly built power transmission lines and increasingly difficult transformation of old lines, and the construction speed of power transmission facilities is lower than the increase speed of load demands. In addition, the local regional power transmission network has the condition of uneven power flow, different line capacities in certain looped networks have the difference, and overload problems are easy to occur on lines with smaller capacities, so that the looped network power transmission potential cannot be well applied. In order to improve the power transmission capability, a certain measure is needed to balance the trend of the power transmission network, and the current common solutions comprise a scheme of installing series compensation and series reactance, a scheme of adding a FACTS element and the like, wherein the FACTS element becomes a hot spot for research and application nowadays due to the flexible power flow control capability. For example, the unified power flow controller has good application effect in 220KV and 500KV networks, but is difficult to be widely applied due to high manufacturing cost. The controllable phase shifter is also a FACTS element which can be used for adjusting power flow distribution of a power transmission network, is mainly a mechanical-plug type phase shifter at home and abroad at present, but has a simpler structure compared with a unified power flow controller and has good application prospect from the aspect of economy.
At present, the research of the controllable phase shifter is considered to be added with an additional damping controller and a fuzzy logic controller to realize the dynamic adjustment of the phase shift angle of the controllable phase shifter so as to inhibit the oscillation of a system, but most of the controllable phase shifter controls the phase shift angle as a continuous quantity, does not embody the discrete switching characteristic of the gear, and is difficult to accurately reflect the dynamic characteristic of the controllable phase shifter. In addition, the discrete switching characteristic of the controllable phase shifter makes the controllable phase shifter easy to have frequent gear switching, and the controller parameter setting of the controllable phase shifter is important for the controllable phase shifter to exert good dynamic adjustment performance.
Disclosure of Invention
The invention aims to solve the technical problems that: the parameter optimization method and the device for the controllable phase shifter additional damping controller are provided, control parameters of the controllable phase shifter additional damping controller are optimized based on a multi-element universe optimization algorithm, suppression of oscillation of a power system is achieved, and stability of the power system is improved.
In order to solve the technical problems, the invention provides a parameter optimization method of an additional damping controller of a controllable phase shifter, which comprises the following steps:
acquiring and establishing a power system model comprising a controllable phase shifter model according to the controllable phase shifter model; wherein the controllable phase shifter model comprises an additional damping controller;
Initializing and setting a multi-universe optimization algorithm according to control parameters of the additional damping controller to generate a preset number of universe individuals, wherein the preset number of universe individuals respectively correspond to the control parameters with different values;
Calculating the objective function value corresponding to each universe individual according to a first preset formula, and sequencing the objective function values corresponding to all universe individuals to obtain a first optimal objective function value;
performing iteration processing on the first optimal objective function value according to preset iteration times until iteration is finished, and acquiring an optimal objective function value under the current iteration;
And updating the control parameters of the additional damping controller according to the optimal objective function value.
Further, the iterative processing is performed on the first optimal objective function value, specifically:
Acquiring and updating algorithm parameters in the multi-universe optimization algorithm according to the current iteration times and combining a second preset formula;
calculating an updated first optimal objective function value according to a third preset formula by executing a roulette mechanism and combining the updated algorithm parameters, and comparing the updated first optimal objective function value with the first optimal objective function value;
If the updated first optimal objective function value is better than the first optimal objective function value, replacing the updated first optimal objective function value with the first optimal objective function value; if the updated first optimal objective function value is not better than the first optimal objective function value, reserving the first optimal objective function value;
Comparing the current iteration times with a preset maximum iteration times, if the current iteration times are smaller than the preset maximum iteration times, returning to the step of acquiring and updating algorithm parameters in the multi-universe optimization algorithm according to the current iteration times and combining a second preset formula.
Further, the damping controller is disposed at the controllable phase shifter gear adjustment.
Further, the invention also provides a parameter optimizing device of the controllable phase shifter additional damping controller, which comprises: the system comprises a model building module, a setting module, an acquisition module, an iteration module and a first updating module;
The model building module is used for obtaining and building an electric power system model comprising the controllable phase shifter model according to the controllable phase shifter model; wherein the controllable phase shifter model comprises an additional damping controller;
the setting module is used for obtaining and initializing a multi-element universe optimization algorithm according to the control parameters of the additional damping controller to generate a preset number of universe individuals, wherein the preset number of universe individuals respectively correspond to the control parameters with different values;
The acquisition module is used for calculating the objective function value corresponding to each universe individual according to a first preset formula, and sequencing the objective function values corresponding to all universe individuals to acquire a first optimal objective function value;
the iteration module is used for carrying out iteration processing on the first optimal objective function value according to preset iteration times until iteration is finished, and obtaining the optimal objective function value under the current iteration;
the first updating module is used for updating the control parameters of the additional damping controller according to the optimal objective function value.
Further, the iteration module is configured to perform an iteration process on the first optimal objective function value, specifically:
Acquiring and updating algorithm parameters in the multi-universe optimization algorithm according to the current iteration times and combining a second preset formula;
calculating an updated first optimal objective function value according to a third preset formula by executing a roulette mechanism and combining the updated algorithm parameters, and comparing the updated first optimal objective function value with the first optimal objective function value;
If the updated first optimal objective function value is better than the first optimal objective function value, replacing the updated first optimal objective function value with the first optimal objective function value; if the updated first optimal objective function value is not better than the first optimal objective function value, reserving the first optimal objective function value;
Comparing the current iteration times with a preset maximum iteration times, if the current iteration times are smaller than the preset maximum iteration times, returning to the step of acquiring and updating algorithm parameters in the multi-universe optimization algorithm according to the current iteration times and combining a second preset formula.
Further, the damping controller in the model building module is arranged at the gear adjustment position of the controllable phase shifter.
Further, the invention also provides a terminal device, which comprises a processor, a memory and a computer program stored in the memory and configured to be executed by the processor, wherein the parameter optimization method of the controllable phase shifter additional damping controller according to any one of the above is realized when the processor executes the computer program.
Further, the invention also provides a computer readable storage medium, which comprises a stored computer program, wherein the computer program controls equipment where the computer readable storage medium is located to execute the parameter optimization method of the controllable phase shifter additional damping controller according to any one of the above.
Compared with the prior art, the parameter optimization method and device for the additional damping controller of the controllable phase shifter have the following beneficial effects:
Establishing a power system model comprising an additional damping controller of a controllable phase shifter; acquiring and according to the control parameters of the additional damper controller, automatically searching the optimal parameters under the set performance requirements based on a multi-universe optimization algorithm, initializing the multi-universe optimization algorithm, and generating a preset number of universe individuals, wherein the preset number of universe individuals respectively correspond to the control parameters with different values; acquiring the optimal objective function value after iteration according to the preset iteration times by calculating the first optimal objective function value in the preset number of universe individuals; and updating the control parameters of the additional damping controller according to the optimal objective function value. Compared with the prior art, the invention optimizes the control parameters of the controllable phase shifter additional damping controller based on the multi-universe optimization algorithm, thereby realizing the suppression of the oscillation of the power system and improving the stability of the power system.
Drawings
FIG. 1 is a flow chart of an embodiment of a method for optimizing parameters of a controllable phase shifter additional damping controller provided by the present invention;
FIG. 2 is a schematic diagram of an embodiment of a parameter optimization apparatus for a controllable phase shifter with an additional damping controller according to the present invention;
FIG. 3 is a schematic diagram of a controllable phase shifter connected to a four-machine two-area system according to an embodiment of a method for optimizing parameters of an additional damping controller of a controllable phase shifter provided by the present invention;
FIG. 4 is a schematic diagram of a configuration of a controllable phase shifter add-on damping controller according to an embodiment of a method for parameter optimization of a controllable phase shifter add-on damping controller provided by the present invention;
FIG. 5 is a schematic diagram of a phase shifter gear adjustment of a controllable phase shifter with a damping controller under a three-phase fault condition in an embodiment of a method for optimizing parameters of an additional damping controller of a controllable phase shifter provided by the present invention;
FIG. 6 is a schematic diagram showing the comparison of active power waveforms of the phase shifter lines before and after the controllable phase shifter is added to the damping controller under the three-phase fault condition in an embodiment of the parameter optimization method of the controllable phase shifter additional damping controller provided by the invention;
FIG. 7 is a graph showing the comparison of the power angle curves of the generator G4 before and after the controllable phase shifter is added to the damping controller under the three-phase fault condition in an embodiment of the parameter optimization method of the additional damping controller of the controllable phase shifter.
Detailed Description
The following description of the embodiments of the present invention will be made clearly and fully with reference to the accompanying drawings, in which it is evident that the embodiments described are only some, but not all embodiments of the invention. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
Example 1
Referring to fig. 1, fig. 1 is a schematic flow chart of an embodiment of a parameter optimization method of an additional damping controller of a controllable phase shifter, as shown in fig. 1, and the method includes steps 101 to 105, specifically as follows:
step 101: acquiring and establishing a power system model comprising a controllable phase shifter model according to the controllable phase shifter model; wherein the controllable phase shifter model comprises an additional damping controller.
In this embodiment, the obtained controllable phase shifter model is a controllable phase shifter model based on thyristor control, and a power system model including the controllable phase shifter is built according to the controllable phase shifter model based on thyristor control, where the power system model is a four-machine two-area system, as shown in fig. 3, in the drawing, G1, G2, G3 and G4 represent generators for distinguishing four generators in the four-machine two-area system, L7 and L8 represent all active loads and inductive loads and reactive loads, C7 and C8 represent capacitive loads and reactive loads, a straight line with a number of distance units represents a power transmission line, a number represents a line length, TCPST represents the controllable phase shifter, and numbers 1 to 10 represent bus numbers.
In this embodiment, the controllable phase shifter model includes an additional damping controller, and as shown in fig. 4, the input signal of the controllable phase shifter refers to the active power of the line provided by the controllable phase shifter, and the reference signal refers to the active power of the line provided by the controllable phase shifter in a steady state, where the steady state power is the phase shifter line power value before the fault occurs; taking the discrete gear characteristics of the controllable phase shifter into consideration, namely the gear of the controllable phase shifter can only be a negative gear or a positive gear within a certain range, and taking the gear of the controllable phase shifter as output.
As an example in this embodiment, the controllable phase shifter range is an integer range of-13 to +13. In fig. 4, the blocks are 8 blocks from top left to top right to bottom left, and the filtering link, the dead zone link, the blocking link, the lead-lag link, the proportional amplifying link, the clipping link and the last two angular shift links are sequentially referred to in this order. In the view of figure 4 of the drawings,Represents the filter time constant, MIN and MAX represent the dead zone limiting parameter,/>Represents the time constant of the blocking link,And/>Representing the lead time constant and the lag time constant,/>, respectivelyRepresenting the proportional gain,/>、/>Representing upper and lower limits of the limiting link and T in the link of changing the angle into the gear position, representing the gear position number of the unidirectional phase shifter,/>Represents the unidirectional adjustable maximum angle of the phase shifter, and is equal to the previous/>The values are the same, INT (X) stands for the rounding operation on X.
As a preferred scheme in this embodiment, the parameters of part of links in the structure of the additional damping controller for the controllable phase shifter include:、/>、/>、=-/> Unidirectional adjustment of the number of stops/> Initial phase shift angle/>。
Step 102: and initializing and setting a multi-universe optimization algorithm according to the control parameters of the additional damping controller to generate a preset number of universe individuals, wherein the preset number of universe individuals respectively correspond to the control parameters with different values.
In this embodiment, for obtaining control parameters of the additional damping controller of the controllable phase shifter in the electric power system, the ITAE performance index and the index of the number of gear adjustments of the controllable phase shifter are mainly used as the basis for screening the control parameters, and the screened control parameters of the additional damping controller of the controllable phase shifter include time constants of the advance linkThe parameter range is [0, 10], the time constant of the hysteresis link/>The parameter range is [0,0.5] and gain/>The parameter range is [0,0.01].
In this embodiment, the algorithm initialization for the multi-element universe optimization algorithm includes: setting the number M of universe individuals, setting the maximum iteration number N and optimizing parameters as the rangeAnd/>Randomly and initially generating a multi-universe population consisting of M universe individuals in a parameter constraint space: /(I)Wherein the ith cosmic individual is/>The parameter constraint space is: /(I)Wherein/>Respectively express/>Minimum value of/>Respectively express/>Is a maximum value of (a). As a preferred scheme in the embodiment, the maximum iteration times and the number of the multi-element universe population set by the multi-element universe optimization algorithm are respectively/>，。
Step 103: and calculating the objective function value corresponding to each universe individual according to a first preset formula, and sequencing the objective function values corresponding to all universe individuals to obtain a first optimal objective function value.
In this embodiment, after a preset number of universe individuals are generated, calculating an objective function value corresponding to each universe individual according to a first preset formula, where the first preset formula is:
；
Wherein, For simulation time,/>Is the bus voltage deviation value,/>Is the power angle deviation value of the generator,/>For the shift frequency of the phase shifter,/>，/>，/>For each objective function, w1, w2, w 3E (0, 1), and/>。
As a preferred scheme in this embodiment, the voltage deviation value of the bus 8 in fig. 3 and the power angle deviation of the generator G4 are selected to construct an objective function, and the weight coefficients respectively have the values of，/>，/>。
In this embodiment, the calculated objective function values of each universe individual are ranked, mainly based on the multi-element universe optimization algorithm, the quality of each universe individual in the optimizing process is evaluated according to the magnitude of the objective function value, the smaller the objective function value of the universe individual is, the better the universe individual is, and the optimal universe individual and the first optimal objective function value corresponding to the universe individual are recorded.
Step 104: and carrying out iteration processing on the first optimal objective function value according to the preset iteration times until the iteration is ended, and obtaining the optimal objective function value under the current iteration.
In this embodiment, since the algorithm is subjected to iterative processing according to the preset iteration number, during the iterative processing, the algorithm parameters in the multi-element universe optimization algorithm need to be acquired and updated according to the current iteration number in combination with the second preset formula, where the algorithm parameters are updated. The algorithm parameters in the updated multi-element universe optimization algorithm are mainly WEP parameters and TDR parameters; the second preset formula is as follows:
；
；
in the method, in the process of the invention, For the current iteration number,/>For maximum iteration number,/>As an optimal scheme in the embodiment, the values of the multi-element universe optimization algorithm are respectively/>。
In this embodiment, since the objective function value is different for each cosmic individual and the objects in the cosmic individual are transferred, this process follows the roulette mechanism, and the position of the cosmic individual is updated by the position update formula. The location update formula is as follows:
；
Wherein: Normalized objective function for the ith universe,/> For range/>Internal random number,/>The j-th parameter of the kth universe is selected for the roulette mechanism.
In this embodiment, the updated multi-element universe optimization algorithm parameters are combined、/>Acquiring an updated universe population according to a third preset formula, calculating an objective function value corresponding to each universe individual in the updated universe population according to the first preset formula in step 103, and sequencing the calculated objective function values of each universe individual to acquire an updated optimal universe individual and a first optimal objective function value corresponding to the optimal universe individual, wherein the third preset formula is as follows:
When (when) When (1):
；
When (when) When (1):
；
Wherein: The j-th parameter that is the best universe formed at present; /(I) Respectively representing the upper and lower bounds of the variable; /(I)For range/>Is a random number of (a) in the memory.
In this embodiment, the updated first optimal objective function value is compared with the first optimal objective function value obtained in step 103; if the updated first optimal objective function value is better than the first optimal objective function value, replacing the updated first optimal objective function value with the first optimal objective function value; and if the updated first optimal objective function value is not better than the first optimal objective function value, reserving the first optimal objective function value. And meanwhile, comparing the current iteration times with the preset maximum iteration times, if the current iteration times are smaller than the preset maximum iteration times, looping the step 104 until the current iteration times reach the maximum iteration times, exiting the iteration loop, and outputting the optimal universe individual and the corresponding optimal objective function value under the current iteration.
As a preferred solution in this embodiment, the control parameters of the controllable phase shifter additional damping controller may be selected by means of interaction between PSCAD and MATLAB. If a PSCAD model is established by inputting related parameters in PSCAD simulation software, performing time domain simulation on the model to obtain an objective function value of each universe individual, after the subsequent iterative updating, inputting updated parameters into the PSCAD model again to obtain an updated simulation result, and transmitting the ordered optimal universe individual and the corresponding optimal objective function value into MATLAB through information interaction to generate a corresponding image.
Step 105: and updating the control parameters of the additional damping controller according to the optimal objective function value.
In this embodiment, the control parameters of the controllable phase shifter additional damping controller are updated according to the obtained optimal objective function value. As a preferred scheme in this embodiment, the optimized control parameters obtained by calculation are、、/>。
In the embodiment, the damping controller based on the multi-element universe optimization algorithm is added at the gear adjustment position of the controllable phase shifter, so that the low-frequency oscillation of the electric power system with the controllable phase shifter can be effectively restrained.
As an example in this embodiment, as shown in fig. 3, in the four-machine two-area system model connected to the controllable phase shifter, if the bus 8 is provided with a three-phase short-circuit fault, the fault time is [2s,2.5s ], the transient gear output result of the phase shifter connected to the parameter-optimized additional damping controller is shown in fig. 5, and after the parameter-optimized damping controller is added, the gear of the controllable phase shifter is discretely changed to suppress oscillation in the fault transient process.
As an example in this embodiment, in the case of a three-phase fault of the bus 8, the active power curves of the phase shifter lines before and after the controllable phase shifter additional damping controller are additionally installed, as shown in fig. 6; the power angle curves of the generator G4 before and after the additional damping controller of the controllable phase shifter are additionally arranged, as shown in fig. 7, the oscillation can be restrained to a certain extent after the additional damping controller is added under the same fault scene, the power of a circuit and the power angle deviation peak value of the generator are effectively reduced, and the stability of the system is improved.
Referring to fig. 2, fig. 2 is a schematic structural diagram of an embodiment of a parameter optimization device of an additional damping controller for a controllable phase shifter, provided by the present invention, as shown in fig. 2, where the structure includes a model building module 201, a setting module 202, an obtaining module 203, an iteration module 204, and a first updating module 205, specifically as follows:
the model building module 201 is configured to obtain and build a power system model including the controllable phase shifter model according to the controllable phase shifter model; wherein the controllable phase shifter model comprises an additional damping controller.
In this embodiment, the obtained controllable phase shifter model is a controllable phase shifter model based on thyristor control, and a power system model including the controllable phase shifter is built according to the controllable phase shifter model based on thyristor control, where the power system model is a four-machine two-area system, as shown in fig. 3, in the drawing, G1, G2, G3 and G4 represent generators for distinguishing four generators in the four-machine two-area system, L7 and L8 represent all active loads and inductive loads and reactive loads, C7 and C8 represent capacitive loads and reactive loads, a straight line with a number of distance units represents a power transmission line, a number represents a line length, TCPST represents the controllable phase shifter, and numbers 1 to 10 represent bus numbers.
In this embodiment, the controllable phase shifter model includes an additional damping controller, and as shown in fig. 4, the input signal of the controllable phase shifter refers to the active power of the line provided by the controllable phase shifter, and the reference signal refers to the active power of the line provided by the controllable phase shifter in a steady state, where the steady state power is the phase shifter line power value before the fault occurs; taking the discrete gear characteristics of the controllable phase shifter into consideration, namely the gear of the controllable phase shifter can only be a negative gear or a positive gear within a certain range, and taking the gear of the controllable phase shifter as output. As an example in this embodiment, the controllable phase shifter range is an integer range of-13 to +13. In fig. 4, the blocks are 8 blocks from top left to top right to bottom left, and the filtering link, the dead zone link, the blocking link, the lead-lag link, the proportional amplifying link, the clipping link and the last two angular shift links are sequentially referred to in this order. In the view of figure 4 of the drawings,Represents the filter time constant, MIN and MAX represent the dead zone limiting parameter,/>Representing the time constant of the blocking link,/>And/>Representing the lead time constant and the lag time constant,/>, respectivelyRepresenting the proportional gain,/>、/>Representing upper and lower limits of the limiting link and T in the link of changing the angle into the gear position, representing the gear position number of the unidirectional phase shifter,/>Represents the unidirectional adjustable maximum angle of the phase shifter, and is equal to the previous/>The values are the same, INT (X) stands for the rounding operation on X.
As a preferred scheme in this embodiment, the parameters of part of links in the structure of the additional damping controller for the controllable phase shifter include:、/>、/>、/>=- Unidirectional adjustment of the number of stops/> Initial phase shift angle/>。
The setting module 202 is configured to obtain and perform initialization setting on a multi-element universe optimization algorithm according to the control parameters of the additional damping controller, so as to generate a preset number of universe individuals, where the preset number of universe individuals respectively correspond to the control parameters with different values.
In this embodiment, for obtaining control parameters of the additional damping controller of the controllable phase shifter in the electric power system, the ITAE performance index and the index of the number of gear adjustments of the controllable phase shifter are mainly used as the basis for screening the control parameters, and the screened control parameters of the additional damping controller of the controllable phase shifter include time constants of the advance linkThe parameter range is [0, 10], the time constant of the hysteresis link/>The parameter range is [0,0.5] and gain/>The parameter range is [0,0.01].
In this embodiment, the algorithm initialization for the multi-element universe optimization algorithm includes: setting the number M of universe individuals, setting the maximum iteration number N and optimizing parameters as the rangeAnd/>Randomly and initially generating a multi-universe population consisting of M universe individuals in a parameter constraint space: /(I)Wherein the ith cosmic individual is/>The parameter constraint space is: /(I)Wherein/>Respectively express/>Minimum value of/>Respectively express/>Is a maximum value of (a). As a preferred scheme in the embodiment, the maximum iteration times and the number of the multi-element universe population set by the multi-element universe optimization algorithm are respectively/>，。
The obtaining module 203 is configured to calculate an objective function value corresponding to each universe individual according to a first preset formula, and order the objective function values corresponding to all universe individuals to obtain a first optimal objective function value.
In this embodiment, after a preset number of universe individuals are generated, calculating an objective function value corresponding to each universe individual according to a first preset formula, where the first preset formula is:
；
Wherein, For simulation time,/>Is the bus voltage deviation value,/>Is the power angle deviation value of the generator,/>For the shift frequency of the phase shifter,/>，/>，/>For each objective function, w1, w2, w 3E (0, 1), and。
As a preferred scheme in this embodiment, the voltage deviation value of the bus 8 in fig. 3 and the power angle deviation of the generator G4 are selected to construct an objective function, and the weight coefficients respectively have the values of，/>，/>。
In this embodiment, the calculated objective function values of each universe individual are ranked, mainly based on the multi-element universe optimization algorithm, the quality of each universe individual in the optimizing process is evaluated according to the magnitude of the objective function value, the smaller the objective function value of the universe individual is, the better the universe individual is, and the optimal universe individual and the first optimal objective function value corresponding to the universe individual are recorded.
The iteration module 204 is configured to perform an iteration process on the first optimal objective function value according to a preset iteration number until the iteration is ended, and obtain an optimal objective function value under the current iteration.
In this embodiment, since the algorithm is subjected to iterative processing according to the preset iteration number, during the iterative processing, the algorithm parameters in the multi-element universe optimization algorithm need to be acquired and updated according to the current iteration number in combination with the second preset formula, where the algorithm parameters are updated. The algorithm parameters in the updated multi-element universe optimization algorithm are mainly WEP parameters and TDR parameters; the second preset formula is as follows:
；
；
in the method, in the process of the invention, For the current iteration number,/>For maximum iteration number,/>As an optimal scheme in the embodiment, the values of the multi-element universe optimization algorithm are respectively/>。
In this embodiment, since the objective function value is different for each cosmic individual and the objects in the cosmic individual are transferred, this process follows the roulette mechanism, and the position of the cosmic individual is updated by the position update formula. The location update formula is as follows:
；
Wherein: Normalized objective function for the ith universe,/> For range/>Internal random number,/>The j-th parameter of the kth universe is selected for the roulette mechanism.
In this embodiment, the updated multi-element universe optimization algorithm parameters are combined、/>Acquiring an updated universe population according to a third preset formula, calculating an objective function value corresponding to each universe individual in the updated universe population according to the first preset formula in step 103, and sequencing the calculated objective function values of each universe individual to acquire an updated optimal universe individual and a first optimal objective function value corresponding to the optimal universe individual, wherein the third preset formula is as follows:
When (when) When (1):
；
When (when) When (1):
；
Wherein: The j-th parameter that is the best universe formed at present; /(I) Respectively representing the upper and lower bounds of the variable; /(I)For range/>Is a random number of (a) in the memory.
In this embodiment, the updated first optimal objective function value is compared with the first optimal objective function value obtained in step 103; if the updated first optimal objective function value is better than the first optimal objective function value, replacing the updated first optimal objective function value with the first optimal objective function value; and if the updated first optimal objective function value is not better than the first optimal objective function value, reserving the first optimal objective function value. And meanwhile, comparing the current iteration times with the preset maximum iteration times, if the current iteration times are smaller than the preset maximum iteration times, looping the step 104 until the current iteration times reach the maximum iteration times, exiting the iteration loop, and outputting the optimal universe individual and the corresponding optimal objective function value under the current iteration.
As a preferred solution in this embodiment, the control parameters of the controllable phase shifter additional damping controller may be selected by means of interaction between PSCAD and MATLAB. If a PSCAD model is established by inputting related parameters in PSCAD simulation software, performing time domain simulation on the model to obtain an objective function value of each universe individual, after the subsequent iterative updating, inputting updated parameters into the PSCAD model again to obtain an updated simulation result, and transmitting the ordered optimal universe individual and the corresponding optimal objective function value into MATLAB through information interaction to generate a corresponding image.
The first updating module 205 is configured to update the control parameter of the additional damping controller according to the optimal objective function value.
In this embodiment, the control parameters of the controllable phase shifter additional damping controller are updated according to the obtained optimal objective function value. As a preferred scheme in this embodiment, the optimized control parameters obtained by calculation are、/>、/>。
In the embodiment, the damping controller based on the multi-element universe optimization algorithm is added at the gear adjustment position of the controllable phase shifter, so that the low-frequency oscillation of the electric power system with the controllable phase shifter can be effectively restrained.
As an example in this embodiment, as shown in fig. 3, in the four-machine two-area system model connected to the controllable phase shifter, if the bus 8 is provided with a three-phase short-circuit fault, the fault time is [2s,2.5s ], the transient gear output result of the phase shifter connected to the parameter-optimized additional damping controller is shown in fig. 5, and after the parameter-optimized damping controller is added, the gear of the controllable phase shifter is discretely changed to suppress oscillation in the fault transient process.
As an example in this embodiment, in the case of a three-phase fault of the bus 8, the active power curves of the phase shifter lines before and after the controllable phase shifter additional damping controller are additionally installed, as shown in fig. 6; the power angle curves of the generator G4 before and after the additional damping controller of the controllable phase shifter are additionally arranged, as shown in fig. 7, the oscillation can be restrained to a certain extent after the additional damping controller is added under the same fault scene, the power of a circuit and the power angle deviation peak value of the generator are effectively reduced, and the stability of the system is improved.
In summary, the invention relates to a parameter optimization method and a device for an additional damping controller of a controllable phase shifter, which are implemented by establishing a power system model comprising the additional damping controller of the controllable phase shifter; acquiring and according to control parameters of an additional damper controller, automatically optimizing parameters under set performance requirements based on a multi-universe optimization algorithm, initializing the multi-universe optimization algorithm, and generating a preset number of universe individuals, wherein the preset number of universe individuals respectively correspond to the control parameters with different values; acquiring an optimal objective function value after iteration according to a preset iteration number by calculating a first optimal objective function value in a universe individual; and updating the control parameters of the additional damping controller according to the optimal objective function value. Compared with the prior art, the invention optimizes the control parameters of the controllable phase shifter additional damping controller based on the multi-universe optimization algorithm, thereby realizing the suppression of the oscillation of the power system and improving the stability of the power system.
The foregoing is merely a preferred embodiment of the present invention, and it should be noted that modifications and substitutions can be made by those skilled in the art without departing from the technical principles of the present invention, and these modifications and substitutions should also be considered as being within the scope of the present invention.
Claims (6)
1. A method for optimizing parameters of an additional damping controller of a controllable phase shifter, comprising:
Acquiring and establishing a power system model comprising a controllable phase shifter model according to the controllable phase shifter model; the controllable phase shifter model comprises an additional damping controller, and is based on thyristor control, and the power system model is a four-machine two-area system;
Initializing and setting a multi-universe optimization algorithm according to control parameters of the additional damping controller to generate a preset number of universe individuals, wherein the preset number of universe individuals respectively correspond to the control parameters with different values;
Calculating the objective function value corresponding to each universe individual according to a first preset formula, and sequencing the objective function values corresponding to all universe individuals to obtain a first optimal objective function value;
Performing iteration processing on the first optimal objective function value according to preset iteration times until iteration is finished, and obtaining an optimal objective function value under the current iteration, wherein the iteration processing is performed on the first optimal objective function value, specifically: acquiring and updating algorithm parameters in the multi-universe optimization algorithm according to the current iteration times and combining a second preset formula; calculating an updated first optimal objective function value according to a third preset formula by executing a roulette mechanism and combining the updated algorithm parameters, and comparing the updated first optimal objective function value with the first optimal objective function value; if the updated first optimal objective function value is better than the first optimal objective function value, replacing the updated first optimal objective function value with the first optimal objective function value; if the updated first optimal objective function value is not better than the first optimal objective function value, reserving the first optimal objective function value; comparing the current iteration times with a preset maximum iteration times, if the current iteration times are smaller than the preset maximum iteration times, returning to the step of acquiring and updating algorithm parameters in the multi-element universe optimization algorithm according to the current iteration times and combining a second preset formula;
updating the control parameters of the additional damping controller according to the optimal objective function value;
The first preset formula is as follows:
；
Wherein, For simulation time,/>Is the bus voltage deviation value,/>Is the power angle deviation value of the generator,/>For the shift frequency of the phase shifter,/>，/>，/>For each objective function, w1, w2, w 3E (0, 1), and/>；
The second preset formula is as follows:
；
；
in the method, in the process of the invention, For the current iteration number,/>For maximum iteration number,/>Is an algorithm hyper-parameter,/>AndOptimizing algorithm parameters in the algorithm for the updated multi-element universe;
the third preset formula is as follows:
When (when) When (1):
；
When (when) When (1):
；
Wherein: The j-th parameter that is the best universe formed at present; /(I) Respectively representing the upper and lower bounds of the variable; /(I)For range/>Is a random number of (a) in the memory.
2. A method for optimizing parameters of a controllable phase shifter additional damping controller as set forth in claim 1, wherein said additional damping controller is provided at the adjustment of the shift position of the controllable phase shifter.
3. A parameter optimizing apparatus for a controllable phase shifter add-on damping controller, comprising: the system comprises a model building module, a setting module, an acquisition module, an iteration module and a first updating module;
The model building module is used for obtaining and building an electric power system model comprising the controllable phase shifter model according to the controllable phase shifter model; the controllable phase shifter model comprises an additional damping controller, and is based on thyristor control, and the power system model is a four-machine two-area system;
the setting module is used for obtaining and initializing a multi-element universe optimization algorithm according to the control parameters of the additional damping controller to generate a preset number of universe individuals, wherein the preset number of universe individuals respectively correspond to the control parameters with different values;
The acquisition module is used for calculating the objective function value corresponding to each universe individual according to a first preset formula, and sequencing the objective function values corresponding to all universe individuals to acquire a first optimal objective function value;
the iteration module is used for carrying out iteration processing on the first optimal objective function value according to preset iteration times until iteration is finished, and obtaining the optimal objective function value under the current iteration;
The iteration module is used for carrying out iteration processing on the first optimal objective function value, and specifically comprises the following steps: acquiring and updating algorithm parameters in the multi-universe optimization algorithm according to the current iteration times and combining a second preset formula; calculating an updated first optimal objective function value according to a third preset formula by executing a roulette mechanism and combining the updated algorithm parameters, and comparing the updated first optimal objective function value with the first optimal objective function value; if the updated first optimal objective function value is better than the first optimal objective function value, replacing the updated first optimal objective function value with the first optimal objective function value; if the updated first optimal objective function value is not better than the first optimal objective function value, reserving the first optimal objective function value; comparing the current iteration times with a preset maximum iteration times, if the current iteration times are smaller than the preset maximum iteration times, returning to the step of acquiring and updating algorithm parameters in the multi-element universe optimization algorithm according to the current iteration times and combining a second preset formula;
the first updating module is used for updating the control parameters of the additional damping controller according to the optimal objective function value;
The first preset formula is as follows:
；
Wherein, For simulation time,/>Is the bus voltage deviation value,/>Is the power angle deviation value of the generator,/>For the shift frequency of the phase shifter,/>，/>，/>For each objective function, w1, w2, w 3E (0, 1), and/>；
The second preset formula is as follows:
；
；
in the method, in the process of the invention, For the current iteration number,/>For maximum iteration number,/>Is an algorithm hyper-parameter,/>AndOptimizing algorithm parameters in the algorithm for the updated multi-element universe;
the third preset formula is as follows:
When (when) When (1):
；
When (when) When (1):
；
Wherein: The j-th parameter that is the best universe formed at present; /(I) Respectively representing the upper and lower bounds of the variable; /(I)For range/>Is a random number of (a) in the memory.
4. A parameter optimization apparatus for an additional damping controller for a controllable phase shifter as set forth in claim 3 wherein said additional damping controller is disposed in said modeling module at an adjustment of a gear of the controllable phase shifter.
5. A terminal device comprising a processor, a memory and a computer program stored in the memory and configured to be executed by the processor, the processor implementing a method for optimizing parameters of a controllable phase shifter additional damping controller according to any one of claims 1 to 2 when the computer program is executed.
6. A computer readable storage medium, characterized in that the computer readable storage medium comprises a stored computer program, wherein the computer program, when run, controls a device in which the computer readable storage medium is located to perform the parameter optimization method of the controllable phase shifter additional damping controller according to any one of claims 1 to 2.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111168548.3A CN113904347B (en) | 2021-09-30 | 2021-09-30 | Parameter optimization method and device for controllable phase shifter additional damping controller |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111168548.3A CN113904347B (en) | 2021-09-30 | 2021-09-30 | Parameter optimization method and device for controllable phase shifter additional damping controller |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113904347A CN113904347A (en) | 2022-01-07 |
| CN113904347B true CN113904347B (en) | 2024-04-23 |
Family
ID=79190136
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111168548.3A Active CN113904347B (en) | 2021-09-30 | 2021-09-30 | Parameter optimization method and device for controllable phase shifter additional damping controller |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113904347B (en) |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109217289A (en) * | 2018-08-24 | 2019-01-15 | 南京理工大学 | A kind of photovoltaic additional damping controller parameter optimization method based on grey wolf algorithm |
| CN109347097A (en) * | 2018-10-30 | 2019-02-15 | 重庆大学 | A kind of double-fed wind power system sub-synchronous oscillation inhibition strategy based on Modified particle swarm optimization algorithm |
| CN110161849A (en) * | 2019-04-22 | 2019-08-23 | 华能澜沧江水电股份有限公司 | A kind of SVC damping of power oscillation controller parameter optimization method based on improvement population |
| CN110611325A (en) * | 2018-06-15 | 2019-12-24 | 南京理工大学 | Wind power plant subsynchronous oscillation suppression method based on particle swarm optimization |
| CN111079290A (en) * | 2019-12-17 | 2020-04-28 | 中国人民解放军海军航空大学 | Engine model correction method based on improved multivariate cosmic algorithm |
| CN111339635A (en) * | 2020-01-17 | 2020-06-26 | 大连大学 | DNA storage coding optimization method of multivariate cosmic algorithm based on damping factor |
| AU2020101603A4 (en) * | 2020-07-31 | 2020-09-10 | Beijing University Of Civil Engineering And Architecture | A Method of Constructing Wide-area Damping Controller with Iterative Identification for Improving Power System Stability |
| CN112993991A (en) * | 2021-03-09 | 2021-06-18 | 东北电力大学 | Fan two-channel damping control low-frequency oscillation wide-area damping control method |
-
2021
- 2021-09-30 CN CN202111168548.3A patent/CN113904347B/en active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110611325A (en) * | 2018-06-15 | 2019-12-24 | 南京理工大学 | Wind power plant subsynchronous oscillation suppression method based on particle swarm optimization |
| CN109217289A (en) * | 2018-08-24 | 2019-01-15 | 南京理工大学 | A kind of photovoltaic additional damping controller parameter optimization method based on grey wolf algorithm |
| CN109347097A (en) * | 2018-10-30 | 2019-02-15 | 重庆大学 | A kind of double-fed wind power system sub-synchronous oscillation inhibition strategy based on Modified particle swarm optimization algorithm |
| CN110161849A (en) * | 2019-04-22 | 2019-08-23 | 华能澜沧江水电股份有限公司 | A kind of SVC damping of power oscillation controller parameter optimization method based on improvement population |
| CN111079290A (en) * | 2019-12-17 | 2020-04-28 | 中国人民解放军海军航空大学 | Engine model correction method based on improved multivariate cosmic algorithm |
| CN111339635A (en) * | 2020-01-17 | 2020-06-26 | 大连大学 | DNA storage coding optimization method of multivariate cosmic algorithm based on damping factor |
| AU2020101603A4 (en) * | 2020-07-31 | 2020-09-10 | Beijing University Of Civil Engineering And Architecture | A Method of Constructing Wide-area Damping Controller with Iterative Identification for Improving Power System Stability |
| CN112993991A (en) * | 2021-03-09 | 2021-06-18 | 东北电力大学 | Fan two-channel damping control low-frequency oscillation wide-area damping control method |
Non-Patent Citations (1)
| Title |
|---|
| UPFC附加阻尼和PSS控制器参数协调优化方法研究;李旭等;南京工程学院学报(自然科学版);20190315(第01期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113904347A (en) | 2022-01-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Mishra | Neural-network-based adaptive UPFC for improving transient stability performance of power system | |
| Hasanvand et al. | Coordinated design of PSS and TCSC to mitigate interarea oscillations | |
| CN109149620B (en) | Self-energy-storage multi-terminal flexible-straight system control method and system | |
| Hassan et al. | Damping of low-frequency oscillations and improving power system stability via auto-tuned PI stabilizer using Takagi–Sugeno fuzzy logic | |
| Hameed et al. | Reduced rule base self-tuning fuzzy PI controller for TCSC | |
| CN108767869B (en) | Static reactive power compensator voltage adjusting method based on artificial neural network | |
| CN108736519B (en) | Self-adaptive control method and device for virtual synchronous generator of photovoltaic power station | |
| CN114362267B (en) | Distributed coordination optimization method for AC/DC hybrid power distribution network considering multi-objective optimization | |
| CN105629730A (en) | UPFC control method based on nerve network sliding mode control | |
| CN111245032B (en) | Voltage prediction control method considering loss reduction optimization of wind power plant collector line | |
| Shakeel et al. | On-line self-tuning adaptive control of an inverter in a grid-tied micro-grid | |
| CN111555358A (en) | Multi-distributed power supply parallel operation method based on self-adaptive droop control | |
| Kuchibhatla et al. | An elephant herding optimization algorithm-based static switched filter compensation scheme for power quality improvement in smart grid | |
| CN115912405A (en) | Adaptive control strategy for virtual synchronous generator in complex oscillation environment | |
| CN112165095B (en) | Low-voltage distribution network reactive power optimization method based on swarm intelligence | |
| CN114389282A (en) | Affine-based uncertainty modal analysis method | |
| CN113904347B (en) | Parameter optimization method and device for controllable phase shifter additional damping controller | |
| CN109120017A (en) | A kind of Method for Reactive Power Optimization in Power based on improvement particle swarm algorithm | |
| CN106655174A (en) | Comprehensive reconstruction optimization method for power distribution network | |
| Loureiro et al. | ANN-based SVC tuning for voltage and harmonics control in microgrids | |
| CN115912393A (en) | Multi-machine parallel VSG system stability improving method based on RBF neural network | |
| CN114597975A (en) | Power distribution method and control equipment for parallel system | |
| Zhang et al. | Adaptive inter-area power oscillation damping from offshore wind farm and MMC-HVDC using deep reinforcement learning | |
| Meena et al. | Control of DSTATCOM in distribution systems by artificial intelligence based controller | |
| CN106684914A (en) | Adaptive PI control method for improving electric power spring pressure regulating performance |
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 |