CN104868468A - UPFC optimization configuration method based on life cycle cost - Google Patents

Abstract

The invention discloses a UPFC optimization configuration method based on life cycle cost. Compared with the prior art, a harmony search method serves as a frame, an ant colony system assesses the adaptability of each harmony unit, the problem of LCC based UPFC optimization configuration is effectively solved, the harmony unit created via the method is characterized by diversity, randomness and the like to avoid local optimization, and the application prospects are wide.

Description

UPFC optimal configuration method based on life cycle cost

Technical Field

The invention belongs to the technical field of operation and control of electric power systems, and particularly relates to a UPFC optimal configuration method based on life cycle cost.

Background

With the development of power grid construction, the connection among regional power grids is more and more compact, and under the current power grid network structure, the safe power trading among the regional power grids is difficult. Flexible ac transmission systems (FACTS) are a new technology appearing in recent years, and the latest development of power electronics technology and modern control technology are applied to realize flexible and rapid control over ac transmission system parameters and network structures, so that reasonable distribution of transmission power is realized, power loss and power generation cost are reduced, and system stability and reliability are greatly improved.

The Unified Power Flow Controller (UPFC) is developed for realizing real-time control and dynamic compensation of an alternating current transmission system, can control line impedance, voltage and power angle, has very powerful functions, can realize functions of parallel compensation, series compensation or phase shifter through change of control quantity, can be used for controlling bus voltage and line power flow, improving system dynamic and transient stability, inhibiting system oscillation, and also can quickly convert working state to adapt to emergency state requirements of a power system, and is most representative in FACTS.

The Life Cycle Cost (LCC) theory is used for coordinated and unified planning and management of the life cycle development process of a project, and is widely recognized and applied in power planning decision-making. The theory originally originated in Sweden, entered China at the end of the last 80 th century, and had many successful applications in the fields of planning and designing of power transformation projects, type selection of transformers and power transmission lines, site selection and volume determination of transformer substations, and the like.

At present, the problems of inconvenient management and complex configuration process of a Unified Power Flow Controller (UPFC) are not solved all the time, and the wide popularization of the UPFC is inhibited.

Disclosure of Invention

The invention aims to solve the problems that the existing Unified Power Flow Controller (UPFC) is inconvenient to manage, complex in configuration process and inhibits wide popularization of the UPFC. The UPFC optimal configuration method based on the life cycle cost takes the harmony search algorithm as a framework, adopts the ant colony system to evaluate the fitness of each harmony individual, better solves the UPFC optimal configuration problem, and has good application prospect.

In order to solve the technical problems, the technical scheme adopted by the invention is as follows:

a UPFC optimal configuration method based on life cycle cost is characterized in that: comprises the following steps of (a) carrying out,

step (1), the UPFC device is accessed to a power grid, a UPFC optimal configuration model based on the whole life cycle cost is established according to a steady-state model of the UPFC device, as shown in a formula (1),

optimized object min.LCC (x)

Constraint h (x) 0 (1)

<math> <mrow> <munder> <mi>g</mi> <mo>&OverBar;</mo> </munder> <mo>&le;</mo> <mi>g</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mo>&le;</mo> <mover> <mi>g</mi> <mo>&OverBar;</mo> </mover> </mrow> </math>

LCC (x) is the life cycle cost, and minP_g、Q_RRespectively the active power and the reactive power generated by the generator, theta and V respectively are the phase angle and the amplitude of the node voltage, k_c、Respectively an amplitude control parameter, a phase angle control parameter, Q of the UPFC controllable voltage source_shThe reactive control parameter is a UPFC; h (x) is an equality constraint condition which is a power balance equation of the alternating current system; g (x) is inequality constraint condition, including voltage amplitude and phase angle of AC system, line transmission power constraint, amplitude parameter and phase angle control parameter of controllable voltage source of UPFC,gis the lower limit of the inequality constraint,is the upper limit of the inequality constraint condition;

step (2), acquiring network parameters of the power system;

step (3), setting tone fine-tuning probability of harmony search algorithmTone fine tuning bandwidthHarmony matrix size H, harmony matrix value probability HMCR, relative importance degree alpha of pheromone information, relative importance degree beta of heuristic information, pheromone volatilization coefficient p, constant t and maximum creation times K_maxThe discrete variable to be optimized is the UPFC installation positionUPFC capacityNumber of iterations k_iter0, an initial harmony matrix is randomly generated according to equation (2)

${\mathrm{HM}}_{k_{iter}}=\left[\begin{array}{ccccc}{x}_1^1& x_2^1& ...& x_{n-1}^1& x_n^1\\ x_1^2& x_2^2& ...& x_{n-1}^2& x_n^2\\ \begin{array}{c}.\\ .\\ .\end{array}& \begin{array}{c}.\\ .\\ .\end{array}& \begin{array}{c}.\\ .\\ .\end{array}& \begin{array}{c}.\\ .\\ .\end{array}& \begin{array}{c}.\\ .\\ .\end{array}\\ x_1^{H-1}& x_2^{H-1}& ...& x_{n-1}^{H-1}& x_n^{H-1}\\ x_1^H& x_2^H& ...& x_{n-1}^H& x_n^H\end{array}\right]---\left(2\right)$

Wherein, $x_i^j=x_{i,min}^j+rand*(x_{i,max}^j-x_{i,min}^j),j=1,2,...,H;$ the number of the ith variable to be optimized in the jth population is represented by subscript n, and the number of the population is represented by superscript H;

step (4), the harmony matrix is obtainedEach row vector carries out load flow calculation containing UPFC once, and calculates the adaptive value of each row vector

Step (5), generating a new harmony matrix

Step (6), adjusting the tone fine-tuning probabilityAnd pitch trimming bandwidth

Step (7), new harmony matrix is obtainedEach row vector carries out load flow calculation containing UPFC once, and calculates the adaptive value of each row vectorAnd will harmony matrixAndmerging into harmony matrix of size 2H

Step (8), harmony matrix with size of 2H through ant colony systemPreferably selects the harmony matrix with the size of H

Step (9), judging whether the iteration times are larger than the maximum creation times K_maxIf the UPFC installation position n is greater than the UPFC installation position n, exiting and outputting the UPFC configuration parameters_pUPFC controllable voltage source amplitude parameter k_cUPFC controllable voltage source phase angle parameterUPFC reactive power control parameter Q_sh(ii) a If not more than the maximum creation times K_maxThen, the number of iterations k is set_iterAnd (5) adding 1 to the value and returning to the step.

The UPFC optimal configuration method based on the life cycle cost is characterized in that: and (2) acquiring network parameters of the power system, including bus serial number, name, negative active power, load reactive power, compensation capacitor, branch number of the power transmission line, serial resistance, serial reactance, parallel conductance, parallel susceptance, transformer transformation ratio and impedance, generator active power output, reactive upper and lower limits and economic parameters.

The invention has the beneficial effects that: compared with the prior art, the UPFC optimal configuration method based on the life cycle cost takes the harmony search algorithm as a frame, adopts the ant colony system to evaluate the fitness of each harmony individual, and is used for solving the UPFC optimal configuration problem based on the LCC.

Drawings

FIG. 1 is a flow chart of a UPFC optimal configuration method based on full lifecycle cost in accordance with the present invention.

Fig. 2 is a schematic diagram of the equivalent injected active power and reactive power of the UPFC of the present invention at node i.

Detailed Description

The invention will be further described with reference to the accompanying drawings.

Compared with the prior art, the UPFC optimal configuration method based on the life cycle cost adopts the harmony search algorithm as a frame, adopts the ant colony system to evaluate the fitness of each harmony individual, is used for solving the UPFC optimal configuration problem based on the LCC, avoids falling into local optimization due to the characteristics of diversity, randomness and the like of the harmony individual during the creation of the algorithm, better solves the UPFC optimal configuration problem based on the LCC, and particularly comprises the following steps,

step (1), the UPFC device is accessed to a power grid, a UPFC optimal configuration model based on the whole life cycle cost is established according to a steady-state model of the UPFC device, as shown in a formula (1),

optimized object min.LCC (x)

Constraint h (x) 0 (1)

<math> <mrow> <munder> <mi>g</mi> <mo>&OverBar;</mo> </munder> <mo>&le;</mo> <mi>g</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mo>&le;</mo> <mover> <mi>g</mi> <mo>&OverBar;</mo> </mover> </mrow> </math>

Where LCC (x) is the full life cycle cost, min. LCC (x) is the minimum full life cycle cost,P_g、Q_Rrespectively the active power and the reactive power generated by the generator, theta and V respectively are the phase angle and the amplitude of the node voltage, k_c、Respectively an amplitude control parameter, a phase angle control parameter, Q of the UPFC controllable voltage source_shThe reactive control parameter is a UPFC; h (x) is an equality constraint condition which is a power balance equation of the alternating current system; g (x) is an inequality constraint condition, which comprises the voltage amplitude and the phase angle of the alternating current system, the transmission power constraint of the line, the amplitude parameter and the phase angle control parameter of the controllable voltage source of the UPFC;gis the lower limit of the inequality constraint,is the upper limit of the inequality constraint.

As shown in FIG. 2, Δ P_ij,jΔQ_ijRespectively equivalent injected active power and reactive power, delta P, of UPFC at node i_ji,jΔQ_jiActive power and reactive power are injected equivalently at node j for the UPFC,the voltage phasors at nodes i and j respectively,is the voltage phasor of the UPFC controllable voltage source,is the current phasor, g, of a UPFC controllable current source_ij、b_ijRespectively, the conductance of the line between nodes i, j andsusceptance, B is the ground admittance of the line,

the basic equation of the UPFC under the per unit system is as follows,

the invention divides the node into common node and UPFC node according to whether the node of the AC system is connected with the UPFC device, because the UPFC device is connected on the common node, the voltage amplitude U of the corresponding control and state variable at the common node_iAnd phase angle theta_iOn the basis, a UPFC variable k is added_c、Q_shWherein k is_cIs the amplitude parameter of the UPFC controllable voltage source,For UPFC controllable voltage source phase angle parameter, Q_shFor the UPFC reactive control parameter, for the UPFC node, the power balance equation is as follows,

<math> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <mrow> <msub> <mi>&Delta;P</mi> <mrow> <mi>u</mi> <mi>p</mi> <mi>f</mi> <mi>c</mi> <mi>k</mi> </mrow> </msub> <mo>=</mo> <msubsup> <mi>P</mi> <mrow> <mi>u</mi> <mi>p</mi> <mi>f</mi> <mi>c</mi> <mi>k</mi> </mrow> <mi>s</mi> </msubsup> <mo>-</mo> <msub> <mi>U</mi> <mrow> <mi>u</mi> <mi>p</mi> <mi>f</mi> <mi>c</mi> <mi>k</mi> </mrow> </msub> <munder> <mo>&Sigma;</mo> <mrow> <mi>j</mi> <mo>&Element;</mo> <mi>J</mi> </mrow> </munder> <mrow> <msub> <mi>U</mi> <mi>j</mi> </msub> <mrow> <mo>(</mo> <mrow> <msub> <mi>G</mi> <mrow> <mi>k</mi> <mi>j</mi> </mrow> </msub> <msub> <mi>cos&theta;</mi> <mrow> <mi>k</mi> <mi>j</mi> </mrow> </msub> <mo>+</mo> <msub> <mi>B</mi> <mrow> <mi>k</mi> <mi>j</mi> </mrow> </msub> <msub> <mi>sin&theta;</mi> <mrow> <mi>k</mi> <mi>j</mi> </mrow> </msub> </mrow> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>&Delta;P</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>&Delta;Q</mi> <mrow> <mi>u</mi> <mi>p</mi> <mi>f</mi> <mi>c</mi> <mi>k</mi> </mrow> </msub> <mo>=</mo> <msubsup> <mi>Q</mi> <mrow> <mi>u</mi> <mi>p</mi> <mi>f</mi> <mi>c</mi> <mi>k</mi> </mrow> <mi>s</mi> </msubsup> <mo>-</mo> <msub> <mi>U</mi> <mrow> <mi>u</mi> <mi>p</mi> <mi>f</mi> <mi>c</mi> <mi>k</mi> </mrow> </msub> <munder> <mo>&Sigma;</mo> <mrow> <mi>j</mi> <mo>&Element;</mo> <mi>J</mi> </mrow> </munder> <mrow> <msub> <mi>U</mi> <mi>j</mi> </msub> <mrow> <mo>(</mo> <mrow> <msub> <mi>G</mi> <mrow> <mi>k</mi> <mi>j</mi> </mrow> </msub> <msub> <mi>sin&theta;</mi> <mrow> <mi>k</mi> <mi>j</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>B</mi> <mrow> <mi>k</mi> <mi>j</mi> </mrow> </msub> <msub> <mi>cos&theta;</mi> <mrow> <mi>k</mi> <mi>j</mi> </mrow> </msub> </mrow> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>&Delta;Q</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>&Delta;P</mi> <mrow> <mi>u</mi> <mi>p</mi> <mi>f</mi> <mi>ct</mi> </mrow> </msub> <mo>=</mo> <msubsup> <mi>P</mi> <mrow> <mi>u</mi> <mi>p</mi> <mi>fct</mi> </mrow> <mi>s</mi> </msubsup> <mo>-</mo> <msub> <mi>U</mi> <mrow> <mi>u</mi> <mi>p</mi> <mi>fct</mi> </mrow> </msub> <munder> <mo>&Sigma;</mo> <mrow> <msup> <mi>j</mi> <mo>&prime;</mo> </msup> <mo>&Element;</mo> <msup> <mi>J</mi> <mo>&prime;</mo> </msup> </mrow> </munder> <mrow> <msub> <mi>U</mi> <msup> <mi>j</mi> <mo>&prime;</mo> </msup> </msub> <mrow> <mo>(</mo> <mrow> <msub> <mi>G</mi> <mrow> <msup> <mi>tj</mi> <mo>&prime;</mo> </msup> </mrow> </msub> <msub> <mi>cos&theta;</mi> <mrow> <msup> <mi>tj</mi> <mo>&prime;</mo> </msup> </mrow> </msub> <mo>+</mo> <msub> <mi>B</mi> <mrow> <msup> <mi>tj</mi> <mo>&prime;</mo> </msup> </mrow> </msub> <msub> <mi>sin&theta;</mi> <mrow> <msup> <mi>tj</mi> <mo>&prime;</mo> </msup> </mrow> </msub> </mrow> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>&Delta;P</mi> <mrow> <mi>j</mi> <mi>i</mi> </mrow> </msub> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>&Delta;Q</mi> <mrow> <mi>u</mi> <mi>p</mi> <mi>fct</mi> </mrow> </msub> <mo>=</mo> <msubsup> <mi>Q</mi> <mrow> <mi>u</mi> <mi>p</mi> <mi>fct</mi> </mrow> <mi>s</mi> </msubsup> <mo>-</mo> <msub> <mi>U</mi> <mrow> <mi>u</mi> <mi>p</mi> <mi>fct</mi> </mrow> </msub> <munder> <mo>&Sigma;</mo> <mrow> <msup> <mi>j</mi> <mo>&prime;</mo> </msup> <mo>&Element;</mo> <msup> <mi>J</mi> <mo>&prime;</mo> </msup> </mrow> </munder> <mrow> <msub> <mi>U</mi> <msup> <mi>j</mi> <mo>&prime;</mo> </msup> </msub> <mrow> <mo>(</mo> <mrow> <msub> <mi>G</mi> <mrow> <msup> <mi>tj</mi> <mo>&prime;</mo> </msup> </mrow> </msub> <msub> <mi>sin&theta;</mi> <mrow> <msup> <mi>tj</mi> <mo>&prime;</mo> </msup> </mrow> </msub> <mo>+</mo> <msub> <mi>B</mi> <mrow> <msup> <mi>tj</mi> <mo>&prime;</mo> </msup> </mrow> </msub> <msub> <mi>cos&theta;</mi> <mrow> <msup> <mi>tj</mi> <mo>&prime;</mo> </msup> </mrow> </msub> </mrow> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>&Delta;Q</mi> <mrow> <mi>j</mi> <mi>i</mi> </mrow> </msub> </mrow> </mrow> </mtd> </mtr> </mtable> </mfenced> </math>

the subscript k represents a UPFC installed at the i end of the node ij of the branch, and the subscript t represents a UPFC installed at the j end of the node ij of the branch; delta P_upfck、ΔQ_upfckUnbalance amounts of i active power and reactive power of a node of a UPFC (unified power flow controller) arranged at an i end of an ij node of a branch circuit respectively; delta P_upfct、ΔQ_upfctUnbalance amounts of j active power and reactive power of a node of a UPFC (unified power flow controller) arranged at an i end of an ij node of a branch circuit respectively;respectively installing UPFC for the i end of a node i of a branch ij, and injecting active power and reactive power into the node i;respectively installing UPFC for the i end of a node of a branch ij, and injecting active power and reactive power into the node j; u shape_upfckSetting the voltage amplitude of the alternating current node of the kth UPFC; u shape_upfctSetting the voltage amplitude of the alternating current node with the t-th UPFC; j represents all nodes connected with the alternating current node provided with the kth UPFC, and J represents the jth alternating current node connected with the alternating current node provided with the kth UPFC; u shape_jThe voltage amplitude of a jth alternating current node connected with an alternating current node provided with a kth UPFC; theta_kjIs the voltage angle difference between the alternating current node provided with the kth UPFC and the jth alternating current node connected with the kth UPFC; g_kj、B_kjRespectively setting the conductance and susceptance between the k-th UPFC AC node and the j-th AC node connected with the k-th UPFC AC node; j ' represents all nodes connected with the alternating current node provided with the t-th UPFC, and J ' represents the J ' th alternating current node connected with the alternating current node provided with the t-th UPFC; u shape_j'The voltage amplitude of the jth' alternating current node connected with the alternating current node provided with the tth UPFC; theta_tj'Is the voltage angle difference between the alternating current node provided with the t-th UPFC and the j' th alternating current node connected with the t-th UPFC; g_tj'、B_tj'Respectively setting the conductance and susceptance between the t-th UPFC AC node and the j' th AC node connected with the t-th UPFC AC node;

the composition of the UPFC life cycle cost here is as follows:

initial investment cost C_I

The initial investment cost, i.e., the cost of the infrastructure, typically includes the purchase cost of the equipment, construction costs, installation costs, and other dynamic costs, and may be generally expressed as

<math> <mrow> <msub> <mi>C</mi> <mi>I</mi> </msub> <mo>=</mo> <munderover> <mo>&Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>m</mi> </munderover> <msub> <mi>p</mi> <mi>i</mi> </msub> </mrow> </math>

Wherein m represents the total number of devices; p is a radical of_iRepresenting the initial investment in the ith equipment. The initial input cost of the power transmission system mainly considers the cost C of the transformer substation_subCost of transmission line C_cabAnd installation cost C_insCompensation equipment cost C_comAnd a platform construction cost C that may be generated according to a site of construction_rigAnd land cost C_landIs that is

C_I＝C_sub+C_cab+C_ins+C_com+C_rig+C_land

Annual maintenance cost C_Mt

The annual maintenance cost mainly refers to the expenses of materials and labor required by operation and maintenance in the whole life cycle of the equipment, and mainly comprises equipment operation expenses, maintenance expenses, operation and maintenance personnel expenses and the like. Because the maintenance period and the cost of the power equipment are relatively stable, the maintenance period and the cost can be generally estimated according to historical average maintenance conditions

C_Mt＝∑N_mt·C_jmt

In the formula, N_mtThe annual average number of overhauls of class j equipment, C_jmtThe average maintenance cost. For the condition that historical overhaul data is difficult to collect, the conversion according to the initial investment cost is often selected

<math> <mrow> <msub> <mi>C</mi> <mrow> <mi>M</mi> <mi>t</mi> </mrow> </msub> <mo>=</mo> <mover> <mo>&Sigma;</mo> <mi>m</mi> </mover> <msub> <mi>f</mi> <mi>m</mi> </msub> <mo>&CenterDot;</mo> <msub> <mi>p</mi> <mi>i</mi> </msub> </mrow> </math>

In the formula (f)_mConverting the operation and maintenance into coefficients; p is a radical of_iRepresenting the initial investment in the ith equipment.

(iii) abandonment cost C_D

The waste cost is a residual value that can be recovered when the power equipment is returned to waste after the end of the life cycle, and can be expressed by the following formula

<math> <mrow> <msub> <mi>C</mi> <mi>D</mi> </msub> <mo>=</mo> <munder> <mo>&Sigma;</mo> <mrow> <mi>i</mi> <mo>&Element;</mo> <mi>M</mi> </mrow> </munder> <mrow> <mo>(</mo> <msub> <mi>C</mi> <mrow> <mi>i</mi> <mi>e</mi> <mi>d</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>C</mi> <mrow> <mi>i</mi> <mi>e</mi> <mi>r</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> </math>

In the formula, C_iedCost for disposal of the equipment i, C_ierThe residual value of the equipment is obtained by conversion by adopting a conversion coefficient according to the original values of different equipment;

year operation cost C_Ot

The annual operating cost is mainly the operating losses of the transmission system per unit time, including substation losses and line losses, and can be expressed as

<math> <mrow> <msub> <mi>C</mi> <mrow> <mi>O</mi> <mi>t</mi> </mrow> </msub> <mo>=</mo> <mover> <mi>&Sigma;</mi> <mi>k</mi> </mover> <mi>&beta;</mi> <mo>&CenterDot;</mo> <msub> <mi>W</mi> <mi>i</mi> </msub> <mo>&CenterDot;</mo> <msub> <mi>T</mi> <mi>i</mi> </msub> <mo>&CenterDot;</mo> <mi>u</mi> <mo>+</mo> <mover> <mi>&Sigma;</mi> <mi>h</mi> </mover> <mi>&rho;</mi> <mo>&CenterDot;</mo> <msub> <mi>W</mi> <mi>j</mi> </msub> <mo>&CenterDot;</mo> <msub> <mi>T</mi> <mrow> <mi>cos</mi> <mtext></mtext> <mi>t</mi> </mrow> </msub> <mo>&CenterDot;</mo> <mi>u</mi> </mrow> </math>

In the formula, k is the total number of lines; beta represents the line loss rate; w_iThe active power transmitted by the ith line is represented; t is_iRepresenting the annual running time of the ith line; u is the average electricity selling price; h represents the total number of substations; rho is the loss of the transformer substation; w_jThe output of the transformer substation is active; t is_costThe maximum annual loss time of the transformer substation;

fault cost of five years C_Ft

The annual fault cost refers to the economic loss caused by the fault, including the power failure loss to the user and the economic loss caused by the fault of the power department, which is mainly related to the occurrence time, duration, frequency of power failure and user type, and the more common calculation method can be represented by the following formula

<math> <mrow> <msub> <mi>C</mi> <mrow> <mi>F</mi> <mi>t</mi> </mrow> </msub> <mo>=</mo> <mover> <mi>&Sigma;</mi> <mi>m</mi> </mover> <msub> <mi>&lambda;</mi> <mi>i</mi> </msub> <mo>&CenterDot;</mo> <mover> <mi>&Sigma;</mi> <mi>h</mi> </mover> <msub> <mi>t</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>&CenterDot;</mo> <mi>R</mi> <mrow> <mo>(</mo> <msub> <mi>t</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> <mo>&CenterDot;</mo> <msub> <mi>W</mi> <mi>j</mi> </msub> </mrow> </math>

Wherein m represents the total number of devices; h represents the total number of substations; lambda [ alpha ]_iIndicating the failure rate of the device i; t is t_ijRepresenting the power failure time of the substation j caused by the failure of the equipment i; w_jThe output of the transformer substation is active; r (t)_ijJ) indicates that substation j corresponds to blackout time t_ijWhen the cost and the specific failure power failure time can not be directly obtained, the cost is usually estimated through the power generation ratio or the electricity selling price and the unavailability rate, and the failure cost is reduced to be

<math> <mrow> <msub> <mi>C</mi> <mrow> <mi>F</mi> <mi>t</mi> </mrow> </msub> <mo>=</mo> <mover> <mi>&Sigma;</mi> <mi>h</mi> </mover> <msub> <mi>W</mi> <mi>j</mi> </msub> <mo>&CenterDot;</mo> <msub> <mi>T</mi> <mi>t</mi> </msub> <mo>&CenterDot;</mo> <mi>&eta;</mi> <mo>&CenterDot;</mo> <mi>u</mi> </mrow> </math>

Where η is the system unavailability, T_tThe annual running time of the system;

LCC of UPFC access scheme

An LCC evaluation model is established according to the initial input cost, the maintenance cost, the abandonment cost, the operation cost and the fault cost, and the expression converted into the current value is

<math> <mrow> <mi>L</mi> <mi>C</mi> <mi>C</mi> <mo>=</mo> <msub> <mi>C</mi> <mi>I</mi> </msub> <mo>+</mo> <munderover> <mo>&Sigma;</mo> <mrow> <mi>t</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>t</mi> <mo>=</mo> <mi>n</mi> </mrow> </munderover> <mfrac> <mrow> <msup> <mrow> <mo>(</mo> <mn>1</mn> <mo>+</mo> <mi>r</mi> <mo>)</mo> </mrow> <mi>n</mi> </msup> <mo>-</mo> <mn>1</mn> </mrow> <mrow> <mi>r</mi> <msup> <mrow> <mo>(</mo> <mn>1</mn> <mo>+</mo> <mi>r</mi> <mo>)</mo> </mrow> <mi>n</mi> </msup> </mrow> </mfrac> <mrow> <mo>(</mo> <msub> <mi>C</mi> <mrow> <mi>O</mi> <mi>t</mi> </mrow> </msub> <mo>+</mo> <msub> <mi>C</mi> <mrow> <mi>M</mi> <mi>t</mi> </mrow> </msub> <mo>+</mo> <msub> <mi>C</mi> <mrow> <mi>F</mi> <mi>t</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>+</mo> <mfrac> <msub> <mi>C</mi> <mi>D</mi> </msub> <msup> <mrow> <mo>(</mo> <mn>1</mn> <mo>+</mo> <mi>r</mi> <mo>)</mo> </mrow> <mi>n</mi> </msup> </mfrac> </mrow> </math>

In the formula, C_I、C_DFor a one-time cost, C_Ot、C_Mt、C_FtIs the average cost in years; r is the discount rate, t is the engineering life cycle, and n is the year;

step (2), acquiring network parameters of the power system; the method comprises the following steps: the number, name, negative active power, load reactive power, compensation capacitance of the bus, the number of a branch of the power transmission line, the number of a head end node and a tail end node, series resistance, series reactance, parallel conductance, parallel susceptance, transformer transformation ratio and impedance, the upper and lower limits of active power output and reactive power of the generator and economic parameters;

step (3) setting the tone fine-tuning probability of the HASTone fine tuning bandwidthHarmony matrix size H, harmony matrix value probability HMCR, relative importance degree alpha of pheromone information, relative importance degree beta of heuristic information, pheromone volatilization coefficient p, constant t and maximum creation times K_maxThe discrete variable to be optimized is the UPFC installation positionUPFC capacityNumber of iterations k_iterWhen the sum of the initial harmony matrix is 0, the initial harmony matrix is randomly generated according to the following formula

${\mathrm{HM}}_{k_{iter}}=\left[\begin{array}{ccccc}{x}_1^1& x_2^1& ...& x_{n-1}^1& x_n^1\\ x_1^2& x_2^2& ...& x_{n-1}^2& x_n^2\\ \begin{array}{c}.\\ .\\ .\end{array}& \begin{array}{c}.\\ .\\ .\end{array}& \begin{array}{c}.\\ .\\ .\end{array}& \begin{array}{c}.\\ .\\ .\end{array}& \begin{array}{c}.\\ .\\ .\end{array}\\ x_1^{H-1}& x_2^{H-1}& ...& x_{n-1}^{H-1}& x_n^{H-1}\\ x_1^H& x_2^H& ...& x_{n-1}^H& x_n^H\end{array}\right]---\left(2\right)$

Wherein, $x_i^j=x_{i,min}^j+rand*(x_{i,max}^j-x_{i,min}^j),j=1,2,...,H;$ the number of the ith variable to be optimized in the jth population is represented by subscript n, and the number of the population is represented by superscript H;

step (4), according to the formula (3), the harmony matrix is alignedOne power flow calculation with the UPFC is performed per row vector,

wherein, Delta U, Delta_upfc、ΔU_upfcRespectively including the voltage phase angle and the voltage amplitude of the common node, the voltage phase angle and the voltage amplitude of the UPFC node, the active residual error delta P and the reactive residual error delta Q of the common node, and the active residual error delta P of the UPFC node_upfcAnd a reactive residual error delta Q_upfcThe following were used:

<math> <mrow> <msub> <mi>&Delta;P</mi> <mrow> <mi>u</mi> <mi>p</mi> <mi>f</mi> <mi>c</mi> </mrow> </msub> <mo>=</mo> <msubsup> <mi>P</mi> <mrow> <mi>u</mi> <mi>p</mi> <mi>f</mi> <mi>c</mi> </mrow> <mi>s</mi> </msubsup> <mo>-</mo> <msub> <mi>U</mi> <mrow> <mi>u</mi> <mi>p</mi> <mi>f</mi> <mi>c</mi> </mrow> </msub> <munder> <mo>&Sigma;</mo> <mrow> <mi>j</mi> <mo>&Element;</mo> <mi>J</mi> </mrow> </munder> <msub> <mi>U</mi> <mi>j</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>G</mi> <mrow> <mi>k</mi> <mi>j</mi> </mrow> </msub> <msub> <mi>cos&theta;</mi> <mrow> <mi>k</mi> <mi>j</mi> </mrow> </msub> <mo>+</mo> <msub> <mi>B</mi> <mrow> <mi>k</mi> <mi>j</mi> </mrow> </msub> <msub> <mi>sin&theta;</mi> <mrow> <mi>k</mi> <mi>j</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>&Delta;P</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> </mrow> </math>

<math> <mrow> <msub> <mi>&Delta;Q</mi> <mrow> <mi>u</mi> <mi>p</mi> <mi>f</mi> <mi>c</mi> </mrow> </msub> <mo>=</mo> <msubsup> <mi>Q</mi> <mrow> <mi>u</mi> <mi>p</mi> <mi>f</mi> <mi>c</mi> </mrow> <mi>s</mi> </msubsup> <mo>-</mo> <msub> <mi>U</mi> <mrow> <mi>u</mi> <mi>p</mi> <mi>f</mi> <mi>c</mi> </mrow> </msub> <munder> <mo>&Sigma;</mo> <mrow> <mi>j</mi> <mo>&Element;</mo> <mi>J</mi> </mrow> </munder> <msub> <mi>U</mi> <mi>j</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>G</mi> <mrow> <mi>k</mi> <mi>j</mi> </mrow> </msub> <msub> <mi>sin&theta;</mi> <mrow> <mi>k</mi> <mi>j</mi> </mrow> </msub> <mo>+</mo> <msub> <mi>B</mi> <mrow> <mi>k</mi> <mi>j</mi> </mrow> </msub> <msub> <mi>cos&theta;</mi> <mrow> <mi>k</mi> <mi>j</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>&Delta;Q</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> </mrow> </math>

the Jacobian matrix of the UPFC nodes is:

<math> <mfenced open = '' close = ''> <mtable> <mtr> <mtd> <mrow> <mi>H</mi> <mo>=</mo> <mo>-</mo> <mfrac> <mrow> <mo>&part;</mo> <msub> <mi>&Delta;P</mi> <mrow> <mi>u</mi> <mi>p</mi> <mi>f</mi> <mi>c</mi> </mrow> </msub> </mrow> <mrow> <mo>&part;</mo> <msub> <mi>&delta;</mi> <mrow> <mi>u</mi> <mi>p</mi> <mi>f</mi> <mi>c</mi> </mrow> </msub> </mrow> </mfrac> </mrow> </mtd> <mtd> <mrow> <mi>J</mi> <mo>=</mo> <mo>-</mo> <mfrac> <mrow> <mo>&part;</mo> <msub> <mi>P</mi> <mrow> <mi>u</mi> <mi>p</mi> <mi>f</mi> <mi>c</mi> </mrow> </msub> </mrow> <mrow> <mo>&part;</mo> <msub> <mi>U</mi> <mrow> <mi>u</mi> <mi>p</mi> <mi>f</mi> <mi>c</mi> </mrow> </msub> </mrow> </mfrac> </mrow> </mtd> </mtr> </mtable> </mfenced> </math>

<math> <mfenced open = '' close = ''> <mtable> <mtr> <mtd> <mrow> <mi>K</mi> <mo>=</mo> <mo>-</mo> <mfrac> <mrow> <mo>&part;</mo> <msub> <mi>&Delta;Q</mi> <mrow> <mi>u</mi> <mi>p</mi> <mi>f</mi> <mi>c</mi> </mrow> </msub> </mrow> <mrow> <mo>&part;</mo> <msub> <mi>&delta;</mi> <mrow> <mi>u</mi> <mi>p</mi> <mi>f</mi> <mi>c</mi> </mrow> </msub> </mrow> </mfrac> </mrow> </mtd> <mtd> <mrow> <mi>L</mi> <mo>=</mo> <mo>-</mo> <mfrac> <mrow> <mo>&part;</mo> <msub> <mi>&Delta;Q</mi> <mrow> <mi>u</mi> <mi>p</mi> <mi>f</mi> <mi>c</mi> </mrow> </msub> </mrow> <mrow> <mo>&part;</mo> <msub> <mi>U</mi> <mrow> <mi>u</mi> <mi>p</mi> <mi>f</mi> <mi>c</mi> </mrow> </msub> </mrow> </mfrac> </mrow> </mtd> </mtr> </mtable> </mfenced> </math>

according to the load flow calculation result, the adaptive value of each row vector is obtained

Step (5), a new harmony matrix is generated according to the following formula

Novel harmony x'_i＝(x′₁,x′₂,…,x′_n) Newly solved first variable x'₁Has a probability of HMCR selected fromInAny value of (a), having a probability of 1-HMCR selected fromAny other value is generated as follows:

<math> <mrow> <msubsup> <mi>x</mi> <mi>i</mi> <mo>&prime;</mo> </msubsup> <mo>=</mo> <mfenced open = '{' close = ''> <mtable> <mtr> <mtd> <mrow> <msubsup> <mi>x</mi> <mi>i</mi> <mo>&prime;</mo> </msubsup> <mo>&Element;</mo> <mrow> <mo>(</mo> <msubsup> <mi>x</mi> <mi>i</mi> <mn>1</mn> </msubsup> <mo>,</mo> <msubsup> <mi>x</mi> <mi>i</mi> <mn>2</mn> </msubsup> <mo>,</mo> <mo>...</mo> <mo>,</mo> <msubsup> <mi>x</mi> <mi>i</mi> <mi>H</mi> </msubsup> <mo>)</mo> </mrow> <mo>,</mo> </mrow> </mtd> <mtd> <mrow> <mi>i</mi> <mi>f</mi> </mrow> </mtd> <mtd> <mrow> <mi>r</mi> <mi>a</mi> <mi>n</mi> <mi>d</mi> <mo>&lt;</mo> <mi>H</mi> <mi>M</mi> <mi>C</mi> <mi>R</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msubsup> <mi>x</mi> <mi>i</mi> <mo>&prime;</mo> </msubsup> <mo>&Element;</mo> <msub> <mi>X</mi> <mi>i</mi> </msub> <mo>,</mo> </mrow> </mtd> <mtd> <mrow> <mi>o</mi> <mi>t</mi> <mi>h</mi> <mi>e</mi> <mi>r</mi> <mi>w</mi> <mi>i</mi> <mi>s</mi> <mi>e</mi> </mrow> </mtd> <mtd> <mrow></mrow> </mtd> </mtr> </mtable> </mfenced> </mrow> </math>

second, if new harmony x'_iFrom the harmony matrixTo fine tune it, the operation is as follows:

<math> <mrow> <msubsup> <mi>x</mi> <mi>i</mi> <mo>&prime;</mo> </msubsup> <mo>=</mo> <mfenced open = '{' close = ''> <mtable> <mtr> <mtd> <mrow> <msubsup> <mi>x</mi> <mi>i</mi> <mo>&prime;</mo> </msubsup> <mo>+</mo> <mi>r</mi> <mi>a</mi> <mi>n</mi> <mi>d</mi> <mo>*</mo> <msub> <mi>BW</mi> <msub> <mi>k</mi> <mrow> <mi>i</mi> <mi>t</mi> <mi>e</mi> <mi>r</mi> </mrow> </msub> </msub> <mo>,</mo> </mrow> </mtd> <mtd> <mrow> <mi>i</mi> <mi>f</mi> </mrow> </mtd> <mtd> <mrow> <mi>r</mi> <mi>a</mi> <mi>n</mi> <mi>d</mi> <mo>&lt;</mo> <msub> <mi>PAR</mi> <msub> <mi>k</mi> <mrow> <mi>i</mi> <mi>t</mi> <mi>e</mi> <mi>r</mi> </mrow> </msub> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msubsup> <mi>x</mi> <mi>i</mi> <mo>&prime;</mo> </msubsup> <mo>,</mo> </mrow> </mtd> <mtd> <mrow> <mi>o</mi> <mi>t</mi> <mi>h</mi> <mi>e</mi> <mi>r</mi> <mi>w</mi> <mi>i</mi> <mi>s</mi> <mi>e</mi> </mrow> </mtd> <mtd> <mrow></mrow> </mtd> </mtr> </mtable> </mfenced> </mrow> </math>

then <math> <mrow> <msub> <mi>HM</mi> <mrow> <msub> <mi>k</mi> <mrow> <mi>i</mi> <mi>t</mi> <mi>e</mi> <mi>r</mi> </mrow> </msub> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>=</mo> <mfenced open = '[' close = ']'> <mtable> <mtr> <mtd> <msubsup> <mi>x</mi> <mn>1</mn> <mo>&prime;</mo> </msubsup> </mtd> </mtr> <mtr> <mtd> <msubsup> <mi>x</mi> <mn>2</mn> <mo>&prime;</mo> </msubsup> </mtd> </mtr> <mtr> <mtd> <mtable> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> </mtable> </mtd> </mtr> <mtr> <mtd> <msubsup> <mi>x</mi> <mi>H</mi> <mo>&prime;</mo> </msubsup> </mtd> </mtr> </mtable> </mfenced> <mo>;</mo> </mrow> </math>

Step (6), the fine tuning probability of the tone is adjusted according to the following formulaAnd pitch trimming bandwidth ${\mathrm{BW}}_{k_{iter}+1}:$

${\mathrm{PAR}}_{k_{iter}+1}={\mathrm{PAR}}_{min}+\frac{{\mathrm{PAR}}_{max}-{\mathrm{PAR}}_{min}}{IM}*k_{iter}$

${\mathrm{BW}}_{k_{iter}+1}={\mathrm{BW}}_{max}\exp (\frac{\ln ({\mathrm{BW}}_{min}/{\mathrm{BW}}_{max})}{IM}*k_{iter})$

Step (7), new harmony matrixEach row vector is subjected to load flow calculation containing UPFC once, and the adaptive value of each row vector is obtained according to the load flow calculation resultAnd will harmony matrixAndcombined to size 2H

Step (8) of using the ant colony system to perform harmony matrix with the size of 2HPreferably selects the harmony matrix with the size of HThe method comprises the following specific steps:

pheromone solubility update: <math> <mrow> <msubsup> <mi>&psi;</mi> <mi>i</mi> <mrow> <msub> <mi>k</mi> <mrow> <mi>i</mi> <mi>t</mi> <mi>e</mi> <mi>r</mi> </mrow> </msub> <mo>+</mo> <mn>1</mn> </mrow> </msubsup> <mo>=</mo> <msubsup> <mi>p&psi;</mi> <mi>i</mi> <msub> <mi>k</mi> <mrow> <mi>i</mi> <mi>t</mi> <mi>e</mi> <mi>r</mi> </mrow> </msub> </msubsup> <mo>+</mo> <msubsup> <mi>&Delta;&psi;</mi> <mi>i</mi> <msub> <mi>k</mi> <mrow> <mi>i</mi> <mi>t</mi> <mi>e</mi> <mi>r</mi> </mrow> </msub> </msubsup> <mo>;</mo> </mrow> </math>

wherein <math> <mrow> <msubsup> <mi>&Delta;&psi;</mi> <mi>i</mi> <msub> <mi>k</mi> <mrow> <mi>i</mi> <mi>t</mi> <mi>e</mi> <mi>r</mi> </mrow> </msub> </msubsup> <mo>=</mo> <mfenced open = '{' close = ''> <mtable> <mtr> <mtd> <mrow> <mfrac> <mi>t</mi> <mrow> <msub> <mi>ed</mi> <mi>i</mi> </msub> </mrow> </mfrac> <mo>,</mo> </mrow> </mtd> <mtd> <mrow> <mi>i</mi> <mi>f</mi> </mrow> </mtd> <mtd> <mrow> <msub> <mi>LCC</mi> <mrow> <mi>k</mi> <mi>i</mi> <mi>t</mi> <mi>e</mi> <mi>r</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>&gt;</mo> <msub> <mi>LCC</mi> <mrow> <mi>k</mi> <mi>i</mi> <mi>t</mi> <mi>e</mi> <mi>r</mi> </mrow> </msub> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mn>0</mn> <mo>,</mo> </mrow> </mtd> <mtd> <mrow> <mi>o</mi> <mi>t</mi> <mi>h</mi> <mi>e</mi> <mi>r</mi> <mi>w</mi> <mi>i</mi> <mi>s</mi> <mi>e</mi> </mrow> </mtd> <mtd> <mrow></mrow> </mtd> </mtr> </mtable> </mfenced> <mo>,</mo> </mrow> </math> ed_iIs thatAndthe Euler distance of each row vector;

determining the size of the probability state transition according to the determined pheromone solubility

<math> <mrow> <msubsup> <mi>PST</mi> <mi>i</mi> <mrow> <msub> <mi>k</mi> <mrow> <mi>i</mi> <mi>t</mi> <mi>e</mi> <mi>r</mi> </mrow> </msub> <mo>+</mo> <mn>1</mn> </mrow> </msubsup> <mo>=</mo> <mfrac> <mrow> <msup> <mrow> <mo>&lsqb;</mo> <msubsup> <mi>&psi;</mi> <mi>i</mi> <mrow> <msub> <mi>k</mi> <mrow> <mi>i</mi> <mi>t</mi> <mi>e</mi> <mi>r</mi> </mrow> </msub> <mo>+</mo> <mn>1</mn> </mrow> </msubsup> <mo>&rsqb;</mo> </mrow> <mi>&alpha;</mi> </msup> <msup> <mrow> <mo>&lsqb;</mo> <msub> <mi>LCC</mi> <mrow> <msub> <mi>k</mi> <mrow> <mi>i</mi> <mi>t</mi> <mi>e</mi> <mi>r</mi> </mrow> </msub> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>&rsqb;</mo> </mrow> <mi>&beta;</mi> </msup> </mrow> <mrow> <msubsup> <mo>&Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>H</mi> </msubsup> <mrow> <msup> <mrow> <mo>&lsqb;</mo> <msubsup> <mi>&psi;</mi> <mi>i</mi> <mrow> <msub> <mi>k</mi> <mrow> <mi>i</mi> <mi>t</mi> <mi>e</mi> <mi>r</mi> </mrow> </msub> <mo>+</mo> <mn>1</mn> </mrow> </msubsup> <mo>&rsqb;</mo> </mrow> <mi>&alpha;</mi> </msup> <msup> <mrow> <mo>&lsqb;</mo> <msub> <mi>LCC</mi> <mrow> <msub> <mi>k</mi> <mrow> <mi>i</mi> <mi>t</mi> <mi>e</mi> <mi>r</mi> </mrow> </msub> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>&rsqb;</mo> </mrow> <mi>&beta;</mi> </msup> </mrow> </mrow> </mfrac> </mrow> </math>

According toOptimizing harmony matrix with size H

Step (9), judging whether the iteration times are larger than the maximum creation times K_maxIf the calculated result is larger than the UPFC configuration parameter, exiting and outputting the result of calculating unconvergence to obtain the UPFC configuration parameter; if not more than the maximum creation times K_maxThen, the number of iterations k is set_iterAnd (5) adding 1 to the value and returning to the step.

The foregoing shows and describes the general principles and broad features of the present invention and advantages thereof. It will be understood by those skilled in the art that the present invention is not limited to the embodiments described above, which are described in the specification and illustrated only to illustrate the principle of the present invention, but that various changes and modifications may be made therein without departing from the spirit and scope of the present invention, which fall within the scope of the invention as claimed. The scope of the invention is defined by the appended claims and equivalents thereof.

Claims (2)

1. The UPFC optimal configuration method based on the life cycle cost is characterized by comprising the following steps: comprises the following steps of (a) carrying out,

step (1), the UPFC device is accessed to a power grid, a UPFC optimal configuration model based on the whole life cycle cost is established according to a steady-state model of the UPFC device, as shown in a formula (1),

optimized object min.LCC (x)

Constraint h (x) 0 (1)

LCC (x) is the life cycle cost, and minP_g、Q_RRespectively the active power and the reactive power generated by the generator, theta and V respectively are the phase angle and the amplitude of the node voltage, k_c、Respectively an amplitude control parameter, a phase angle control parameter, Q of the UPFC controllable voltage source_shThe reactive control parameter is a UPFC; h (x) is an equality constraint condition which is a power balance equation of the alternating current system; g (x) is inequality constraint condition, including voltage amplitude and phase angle of AC system, line transmission power constraint, amplitude parameter and phase angle control parameter of controllable voltage source of UPFC,gis the lower limit of the inequality constraint,is the upper limit of the inequality constraint condition;

step (2), acquiring network parameters of the power system;

step (3), setting tone fine-tuning probability of harmony search algorithmTone fine tuning bandwidthHarmony matrix size H, harmony matrix value probability HMCR, relative importance degree alpha of pheromone information, relative importance degree beta of heuristic information, pheromone volatilization coefficient p, constant t and maximum creation times K_maxThe discrete variable to be optimized is the UPFC installation positionUPFC capacityNumber of iterations k_iter0, an initial harmony matrix is randomly generated according to equation (2)

Wherein, the number of the ith variable to be optimized in the jth population is represented by subscript n, and the number of the population is represented by superscript H;

step (4), the harmony matrix is obtainedEach row vector carries out load flow calculation containing UPFC once, and calculates the adaptive value of each row vector

Step (5), generating a new harmony matrix

Step (6), adjusting the tone fine-tuning probabilityAnd pitch trimming bandwidth

Step (7)New harmony matrixEach row vector carries out load flow calculation containing UPFC once, and calculates the adaptive value of each row vectorAnd will harmony matrixAndmerging into harmony matrix of size 2H

Step (8), harmony matrix with size of 2H through ant colony systemPreferably selects the harmony matrix with the size of H

Step (9), judging whether the iteration times are larger than the maximum creation times K_maxIf the UPFC installation position n is greater than the UPFC installation position n, exiting and outputting the UPFC configuration parameters_pUPFC controllable voltage source amplitude parameter k_cUPFC controllable voltage source phase angle parameterUPFC reactive power control parameter Q_sh(ii) a If not more than the maximum creation times K_maxThen, the number of iterations k is set_iterAnd (5) adding 1 to the value and returning to the step.

2. The full life cycle cost based UPFC optimal configuration method of claim 1, wherein: and (2) acquiring network parameters of the power system, including bus serial number, name, negative active power, load reactive power, compensation capacitor, branch number of the power transmission line, serial resistance, serial reactance, parallel conductance, parallel susceptance, transformer transformation ratio and impedance, generator active power output, reactive upper and lower limits and economic parameters.