CN107528323B - Optimal configuration method of dynamic reactive power compensation device - Google Patents

Optimal configuration method of dynamic reactive power compensation device Download PDF

Info

Publication number
CN107528323B
CN107528323B CN201711001724.8A CN201711001724A CN107528323B CN 107528323 B CN107528323 B CN 107528323B CN 201711001724 A CN201711001724 A CN 201711001724A CN 107528323 B CN107528323 B CN 107528323B
Authority
CN
China
Prior art keywords
reactive power
voltage
power compensation
compensation device
fault
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201711001724.8A
Other languages
Chinese (zh)
Other versions
CN107528323A (en
Inventor
陆佳政
朱思国
谭艳军
李波
黄清军
朱远
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
State Grid Corp of China SGCC
State Grid Hunan Electric Power Co Ltd
Disaster Prevention and Mitigation Center of State Grid Hunan Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
State Grid Hunan Electric Power Co Ltd
Disaster Prevention and Mitigation Center of State Grid Hunan Electric Power Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by State Grid Corp of China SGCC, State Grid Hunan Electric Power Co Ltd, Disaster Prevention and Mitigation Center of State Grid Hunan Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN201711001724.8A priority Critical patent/CN107528323B/en
Publication of CN107528323A publication Critical patent/CN107528323A/en
Application granted granted Critical
Publication of CN107528323B publication Critical patent/CN107528323B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/12Circuit arrangements for ac mains or ac distribution networks for adjusting voltage in ac networks by changing a characteristic of the network load
    • H02J3/16Circuit arrangements for ac mains or ac distribution networks for adjusting voltage in ac networks by changing a characteristic of the network load by adjustment of reactive power
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0631Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • G06Q10/06312Adjustment or analysis of established resource schedule, e.g. resource or task levelling, or dynamic rescheduling
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/18Arrangements for adjusting, eliminating or compensating reactive power in networks
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E40/00Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/30Reactive power compensation

Landscapes

  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Engineering & Computer Science (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Power Engineering (AREA)
  • Strategic Management (AREA)
  • Economics (AREA)
  • Game Theory and Decision Science (AREA)
  • Educational Administration (AREA)
  • Development Economics (AREA)
  • Marketing (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Tourism & Hospitality (AREA)
  • Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

The invention discloses an optimal configuration method of a dynamic reactive power compensation device, which comprises the steps of determining a transient voltage instability hidden danger area; selecting a fault set, and simultaneously preliminarily determining candidate access points of the dynamic reactive power compensation device; acquiring the transient response condition of each voltage grade bus in the transient voltage instability hidden danger area; counting the capacity demand range of each candidate access point of the dynamic reactive power compensation device; determining a final optimal dynamic reactive power compensation device access node; recording the voltage sag amplitude and the duration of each node in the transient voltage instability hidden danger area; and calculating the return on investment of each candidate capacity and selecting the final economic optimal dynamic reactive power compensation device capacity. The method of the invention considers the return on investment of the dynamic reactive power compensation device and the result of stable influence on the power grid, thereby exerting the function of the dynamic reactive power compensation device to the maximum extent, realizing smaller calculation amount on the premise of ensuring the result accuracy, and having scientific and reasonable method and obvious effect.

Description

Optimal configuration method of dynamic reactive power compensation device
Technical Field
The invention particularly relates to an optimal configuration method of a dynamic reactive power compensation device.
Background
With the development of national economic technology and the improvement of living standard of people, electric energy becomes essential secondary energy in daily production and life of people, and brings endless convenience to production and life of people.
Considering the inverse distribution of energy structure and load in China, along with the implementation and deepening of the west-east power transmission strategy, a power system has already stepped into a large power grid, ultrahigh-voltage and long-distance transmission era, and five receiving-end power grids of east China, northeast China, Yue Min Qiong and Jingjin Jilu are formed on the national level. The voltage stability problem is increasingly highlighted by the continuous increase of the proportion of external power supplies of large receiving-end power grids, the rapid increase of load levels and the increasing complication of load characteristics. Dynamic reactive power compensation devices represented by static synchronous compensators (STATCOM) are increasingly widely applied to power grids by virtue of the advantages of faster response speed and better performance compared with the traditional reactive power compensation devices. But considering the higher configuration cost, how to optimize the dynamic reactive power compensation configuration scheme to fully play the role of the dynamic reactive power compensation device has great practical significance.
The existing reactive power compensation optimization method based on static voltage stability analysis does not consider the optimal efficacy of the dynamic reactive power compensation device in the fault dynamic process; transient voltage stability analysis optimization schemes based on Impedance Module Margin Index (IMMI), Track Sensitivity Index (TSI) and the like generally have the problem of overlarge calculated amount, and most of the schemes do not consider the huge loss of sensitive load caused by voltage sag with continuously increased occurrence probability in the modern power grid; the existing probability statistical model considering the voltage sag risk also has the problem of large calculation amount, and the influence of different load models in different running states of a power grid on the transient voltage stability analysis result is not considered, so that the final result is incomplete.
Disclosure of Invention
The invention aims to provide an optimal configuration method of a dynamic reactive power compensation device, which can exert the function of the dynamic reactive power compensation device to the maximum extent, is scientific and reasonable and has obvious effect.
The optimal configuration method of the dynamic reactive power compensation device provided by the invention comprises the following steps:
s1, determining a transient voltage instability hidden danger area in a research area according to a power grid structure of the research area and an installed reactive power compensation device;
s2, carrying out simulation analysis on the line load in the transient voltage instability hidden danger area determined in the step S1, outputting a line list with the line load exceeding a rated value by 80.0%, selecting heavy current line faults from the obtained line list to form a fault set, and simultaneously preliminarily determining candidate access points of the dynamic reactive power compensation device in the hidden danger area according to the setting principle of the dynamic reactive power compensation device;
s3, testing the transient response conditions of the buses of each voltage class in the transient voltage instability hidden danger area under all fault conditions in the fault set determined in the step S2 according to the corresponding load proportion models of the induction motors in the research areas under different running states;
s4, according to the transient response conditions of the buses of each voltage class in the transient voltage instability hidden danger area under the fault condition obtained in the step S3, when the candidate access points of the dynamic reactive power compensation device determined in the step S2 are respectively provided with the dynamic reactive power compensation device, a trial and error method is adopted to obtain reactive power compensation capacity required to be arranged under different load models and different fault conditions, and meanwhile, the capacity requirement range of the candidate access points of each dynamic reactive power compensation device is counted;
s5, selecting an access node with the minimum reactive compensation capacity demand range according to the capacity demand range of each candidate access point of the dynamic reactive compensation device obtained in the step S4, and determining the access node as the access node of the final dynamic reactive compensation device;
s6, according to the access node determined in the step S5, selecting a plurality of groups of capacities common in industrial production from the capacity demand range as candidate capacities, and simultaneously recording the voltage sag amplitude and the duration of each node in the transient voltage instability hidden danger area under the condition of each candidate capacity and all faults in a fault set;
and S7, calculating the investment return rate of each candidate capacity according to the voltage sag amplitude value and the duration obtained in the step S6, and selecting the capacity with the highest investment return rate as the final capacity of the dynamic reactive power compensation device.
Step S1, determining a transient voltage instability hidden danger area in the research area, specifically, analyzing a voltage distribution condition in a maximum operation mode in summer and a maximum operation mode in winter and an active power and reactive power transmission condition in each internal region based on a power grid structure of the research area and an installed reactive power compensation device, and thus determining a power grid receiving end of a low voltage area and a large reactive gap as the transient voltage instability hidden danger area.
Step S4, obtaining the reactive compensation capacity by using a trial and error method, specifically performing trial and error on the reactive compensation capacity by using the following rules: the reactive compensation capacity must ensure that after all faults in a fault set occur, the voltage of the bus with the slowest voltage recovery in all buses of 500kV, 220kV, 110kV and 35kV in a transient voltage instability hidden danger area can be raised back to 0.75p.u. within 1s after the faults occur.
Step S7, calculating the return on investment for each candidate capacity, specifically, calculating the return on investment ROI by using the following formula:
Figure BDA0001443520440000031
wherein ROI is the return on investment, R0Annual voltage sag loss cost, R, before installation of dynamic reactive power compensation devicesSThe annual voltage sag loss cost after the reactive power compensation device is installed is shown, and the investment cost for installing the reactive power compensation device is shown as I.
The annual voltage sag loss cost is obtained by replacing the actual sag cost with the sum of annual sag costs of all buses under the fault concentrated accidents (the fault concentrated accidents are all key accidents).
The annual voltage sag loss cost is calculated by adopting the following formula:
Figure BDA0001443520440000041
wherein R is annual voltage sag loss cost, N is annual voltage sag times, N is total number of fault types, and P isiThe ratio of the number of voltage sags caused by the i-th fault to the number of annual sags, EiThe voltage sag loss cost under the i-th fault is obtained.
The investment cost for installing the reactive power compensation device is obtained by adopting the following calculation formula:
I=ρ(γ×f+×q)
wherein I is the investment cost for installing the reactive compensation device; rho is the return on investment; gamma is a mark for judging whether the dynamic reactive power compensation device is installed, and is 1 when installed and 0 when not installed; f is equipment installation cost; unit price of equipment device; and q is the capacity for installing the dynamic reactive power compensation device.
The loss cost E of voltage sag under the i-th faultiSpecifically, the following formula is adopted for calculation:
Figure BDA0001443520440000042
wherein m is the total number of all 500kV, 220kV, 110kV and 35kV buses in a research area, EijTo investigate the loss cost of voltage sag at the jth bus in the area, and EijIs calculated by the formula
Figure BDA0001443520440000043
K is the maximum economic loss value of a single voltage sag, sigma2Is a sensitivity coefficient (the value is 0.3 according to the actual condition), QijThe severity of the voltage sag at the jth bus under the ith type of fault.
Severity of the voltage sag QijThe method specifically comprises the following steps:
1) drawing a coordinate graph by taking the time T as an abscissa and the voltage U as an ordinate; with U being equal to Umax(maximum value of device voltage tolerance amplitude), U ═ Umin(minimum value of device voltage tolerance amplitude), T ═ Tmax(maximum value of device voltage withstand duration) and T ═ Tmin(minimum value of device voltage withstand duration) four straight lines are drawn, dividing the first quadrant of the graph into three regions:
normal zone: the abscissa T is less than or equal to TminOr the ordinate U is greater than or equal to UmaxThe area of (a);
a fault area: the ordinate U is less than or equal to UminAnd the abscissa T is greater than or equal to TmaxThe area of (a);
fuzzy areas: removing the remaining areas of the normal area and the defective area;
the voltage endurance duration is defined as the time required from the occurrence of a fault to the time when the voltage rises back to 0.9p.u., and the voltage sag amplitude is defined as the average value of the effective values of the voltage within the voltage sag duration.
2) The severity of the voltage sag, Q, was calculated as followsij
Definition of UijIs the voltage sag amplitude, T, at the jth bus under the ith faultijFor the voltage withstand duration at the jth bus under a class i fault,
if U isij、TijFalls within the normal region of the graph, then Qij=0;
If U isij、TijFalls within the fault region of the graph, then Qij=1;
If U isij、TijFall into the fuzzy area of the coordinate graph, then
Figure BDA0001443520440000051
The U ismaxIs 0.9 p.u.; u shapeminThe value of (a) is 0.6 p.u.; t ismaxIs 600 ms; t isminIs 20 ms.
The optimal configuration method of the dynamic reactive power compensation device can exert the function of the dynamic reactive power compensation device to the maximum extent, and is scientific, reasonable and remarkable in effect; the influence of different load proportion models corresponding to different operation modes on the transient voltage stability analysis result, huge loss caused by voltage sag with continuously increased occurrence probability in a modern power grid to sensitive loads and the importance of investment return rate under the high-cost background of the dynamic reactive power compensation device are considered, meanwhile, small calculated amount on the premise of ensuring result accuracy is achieved, and the dynamic reactive power compensation device has practicability and superiority from the theoretical and engineering aspects.
Drawings
FIG. 1 is a flow chart of the method of the present invention.
FIG. 2 is a state region diagram illustrating the severity of voltage sag in the method of the present invention.
Detailed Description
FIG. 1 shows a flow chart of the method of the present invention: the optimal configuration method of the dynamic reactive power compensation device provided by the invention comprises the following steps:
s1, determining a transient voltage instability hidden danger area in a research area according to a power grid structure of the research area and an installed reactive power compensation device; the method specifically comprises the steps of analyzing the voltage distribution condition in the summer maximum operation mode and the winter maximum operation mode and the active power and reactive power transmission condition in each internal region based on a power grid structure of a research area and an installed reactive power compensation device, and accordingly determining the power grid receiving end of a low-voltage area and a large reactive power gap as a transient voltage instability hidden danger area.
S2, carrying out simulation analysis on the line loads related to the transient voltage instability hidden danger area determined in the step S1, outputting a line list with the line loads exceeding a rated value by 80.0%, selecting heavy load flow line faults in the line list to form a fault set, and simultaneously primarily determining a candidate access point of the dynamic reactive power compensation device in the hidden danger area according to the setting principle of the dynamic reactive power compensation device;
s3, testing the transient response conditions of the buses of each voltage class in the transient voltage instability hidden danger area under all fault conditions in the fault set determined in the step S2 according to the corresponding load proportion models of the induction motors in the research areas under different running states;
s4, according to the transient response conditions of the buses of each voltage class in the transient voltage instability hidden danger area under the fault condition obtained in the step S3, when the candidate access points of the dynamic reactive power compensation device determined in the step S2 are respectively provided with the dynamic reactive power compensation device, a trial and error method is adopted to obtain reactive power compensation capacity required to be arranged under different load models and different fault conditions, and meanwhile, the capacity requirement range of the candidate access points of each dynamic reactive power compensation device is counted;
when the device is tried, the reactive compensation capacity must ensure that after all faults are concentrated and all faults occur, the voltage of the bus with the slowest recovery in all buses of 500kV, 220kV, 110kV and 35kV in the transient voltage instability hidden danger area can be raised back to 0.75p.u. within 1s after the faults occur.
S5, according to the capacity demand range of each candidate access point of the dynamic reactive power compensation device obtained in the step S4, determining the access node with the minimum reactive power compensation capacity demand range as the access node with the optimal installation effect, and determining the access node as the final dynamic reactive power compensation device;
s6, according to the access node determined in the step S5, selecting a plurality of groups of capacities common in industrial production from the capacity demand range as candidate capacities, and simultaneously recording the voltage sag amplitude and the duration of each node in the transient voltage instability hidden danger area under the condition of each candidate capacity and all faults in a fault set;
s7, calculating the investment return rate of each candidate capacity according to the voltage sag amplitude value and the duration obtained in the step S6, and selecting the capacity with the highest investment return rate as the final capacity of the dynamic reactive power compensation device;
the return on investment is specifically calculated by adopting the following formula:
Figure BDA0001443520440000071
ROI in the formulaTo return on investment, R0Annual voltage sag loss cost, R, before installation of dynamic reactive power compensation devicesSThe annual voltage sag loss cost after the reactive power compensation device is installed is shown, and I is the investment cost for installing the reactive power compensation device; the annual voltage sag loss cost is obtained by replacing the actual sag cost with the sum of annual sag costs of each bus in the case of a fault-concentrated accident (the fault-concentrated accidents are all key accidents), and can be specifically calculated by adopting the following formula:
Figure BDA0001443520440000072
wherein R is annual voltage sag loss cost, N is annual voltage sag times, N is total number of fault types, and P isiThe ratio of the number of voltage sags caused by the i-th fault to the number of annual sags, EiThe voltage sag loss cost under the i-th type fault;
the investment cost for installing the reactive power compensation device is calculated by adopting the following formula:
I=ρ(γ×f+×q)
wherein I is the investment cost for installing the reactive compensation device; rho is the return on investment; gamma is a mark for judging whether the dynamic reactive power compensation device is installed, and is 1 when installed and 0 when not installed; f is equipment installation cost; unit price of equipment device; q is the capacity for installing the dynamic reactive power compensation device;
loss cost of voltage sag under i-th fault EiSpecifically, the following formula is adopted for calculation:
Figure BDA0001443520440000081
wherein m is the total number of all 500kV, 220kV, 110kV and 35kV buses in a research area, EijTo investigate the loss cost of voltage sag at the jth bus in the area, and EijIs calculated by the formula
Figure BDA0001443520440000082
K is the maximum economy of a single voltage sagLoss value, σ2Is a sensitivity coefficient (the value is 0.3 according to the actual condition), QijThe severity of the voltage sag at the jth bus under the ith type of fault.
Severity of voltage sag QijThe method specifically comprises the following steps:
1) drawing a coordinate graph by taking the time T as an abscissa and the voltage U as an ordinate; with U being equal to Umax(maximum value of device voltage tolerance amplitude), U ═ Umin(minimum value of device voltage tolerance amplitude), T ═ Tmax(maximum value of device voltage withstand duration) and T ═ Tmin(minimum value of device voltage withstand duration) four straight lines are drawn, dividing the first quadrant of the graph into three regions:
normal zone: the abscissa T is less than or equal to TminOr the ordinate U is greater than or equal to UmaxThe area of (a);
a fault area: the ordinate U is less than or equal to UminAnd the abscissa T is greater than or equal to TmaxThe area of (a);
fuzzy areas: removing the remaining areas of the normal area and the defective area;
the voltage endurance duration is defined as the time required from the occurrence of a fault to the time when the voltage rises back to 0.9p.u., and the voltage sag amplitude is defined as the average value of the effective values of the voltage within the voltage sag duration.
2) The severity of the voltage sag, Q, was calculated as followsij
Definition of UijIs the voltage sag amplitude, T, at the jth bus under the ith faultijFor the voltage withstand duration at the jth bus under a class i fault,
if U isij、TijFalls within the normal region of the graph, then Qij=0;
If U isij、TijFalls within the fault region of the graph, then Qij=1;
If U isij、TijFall into the fuzzy area of the coordinate graph, then
Figure BDA0001443520440000091
U in the specific implementationmaxIs 0.9 p.u.; u shapeminThe value of (a) is 0.6 p.u.; t ismaxIs 600 ms; t isminIs 20 ms.

Claims (10)

1. An optimal configuration method of a dynamic reactive power compensation device comprises the following steps:
s1, determining a transient voltage instability hidden danger area in a research area according to a power grid structure of the research area and an installed reactive power compensation device;
s2, carrying out simulation analysis on the line loads related to the transient voltage instability hidden danger area determined in the step S1, outputting a line list with the line loads exceeding a rated value by 80.0%, selecting heavy current line faults from the obtained line list to form a fault set, and simultaneously primarily determining candidate access points of the dynamic reactive power compensation device in the hidden danger area according to the setting principle of the dynamic reactive power compensation device;
s3, testing the transient response conditions of the buses of each voltage class in the transient voltage instability hidden danger area under all fault conditions in the fault set determined in the step S2 according to the corresponding load proportion models of the induction motors in the research areas under different running states;
s4, according to the transient response conditions of the buses of each voltage class in the transient voltage instability hidden danger area under the fault condition obtained in the step S3, when the candidate access points of the dynamic reactive power compensation device determined in the step S2 are respectively provided with the dynamic reactive power compensation device, a trial and error method is adopted to obtain reactive power compensation capacity required to be arranged under different load models and different fault conditions, and meanwhile, the capacity requirement range of the candidate access points of each dynamic reactive power compensation device is counted;
s5, according to the capacity demand range of each candidate access point of the dynamic reactive power compensation device obtained in the step S4, selecting the access node with the minimum reactive power compensation capacity demand range and determining the access node as the access node of the final dynamic reactive power compensation device;
s6, selecting a plurality of groups of candidate capacities from the capacity demand range according to the access node determined in the step S5, and simultaneously recording the voltage sag amplitude and the duration of each node in the transient voltage instability hidden danger area under the condition of each candidate capacity and all faults in the fault set;
and S7, calculating the investment return rate of each candidate capacity according to the voltage sag amplitude value and the duration obtained in the step S6, and selecting the capacity with the highest investment return rate as the final capacity of the dynamic reactive power compensation device.
2. The method according to claim 1, wherein the step S1 is to determine a potential transient voltage instability area in the research area, and specifically, based on the grid structure of the research area and the installed reactive power compensation devices, analyze the voltage distribution and the transmission of active power and reactive power between internal areas in the maximum operation mode in summer and the maximum operation mode in winter, so as to determine the grid receiving end of a low voltage area and a large reactive gap as the potential transient voltage instability area.
3. The method for optimally configuring a dynamic reactive power compensation device according to claim 2, wherein the reactive power compensation capacity is obtained by a trial and error method in step S4, specifically, the following rules are adopted for trial and error of the reactive power compensation capacity: the reactive compensation capacity must ensure that after all faults in a fault set occur, the voltage of the bus with the slowest voltage recovery in all buses of 500kV, 220kV, 110kV and 35kV in a transient voltage instability hidden danger area can be raised back to 0.75p.u. within 1s after the faults occur.
4. The method according to claim 3, wherein the step S7 is to calculate the return on investment for each candidate capacity, specifically, the return on investment ROI is calculated by using the following formula:
Figure FDA0001443520430000021
wherein ROI is the return on investment, R0Annual voltage sag loss cost, R, before installation of dynamic reactive power compensation devicesSThe annual voltage sag loss cost after the reactive power compensation device is installed is shown, and the investment cost for installing the reactive power compensation device is shown as I.
5. The method according to claim 4, wherein the annual voltage sag loss cost is obtained by replacing an actual sag cost with a sum of annual sag costs of each bus in the event of a fault concentration.
6. The optimal configuration method of a dynamic reactive power compensation device according to claim 5, wherein the annual voltage sag loss cost is calculated by using the following formula:
Figure FDA0001443520430000031
wherein R is annual voltage sag loss cost, N is annual voltage sag times, N is total number of fault types, and P isiThe ratio of the number of voltage sags caused by the i-th fault to the number of annual sags, EiThe voltage sag loss cost under the i-th fault is obtained.
7. The method according to claim 6, wherein the investment cost for installing the reactive power compensator is calculated by the following equation:
I=ρ(γ×f+×q)
wherein I is the investment cost for installing the reactive compensation device; rho is the return on investment; gamma is a mark for judging whether the dynamic reactive power compensation device is installed, and is 1 when installed and 0 when not installed; f is equipment installation cost; unit price of equipment device; and q is the capacity for installing the dynamic reactive power compensation device.
8. The method according to claim 7, wherein the i-th fault is a faultLoss of voltage sag EiSpecifically, the following formula is adopted for calculation:
Figure FDA0001443520430000032
wherein m is the total number of all 500kV, 220kV, 110kV and 35kV buses in a research area, EijTo investigate the loss cost of voltage sag at the jth bus in the area, and EijIs calculated by the formula
Figure FDA0001443520430000033
K is the maximum economic loss value of a single voltage sag, sigma2As a coefficient of sensitivity, QijThe severity of the voltage sag at the jth bus in the case of a class i fault.
9. The method according to claim 8, wherein the severity of the voltage sag Q is determined by a method of optimizing the configuration of the dynamic var compensator apparatusijThe method specifically comprises the following steps:
1) drawing a coordinate graph by taking the time T as an abscissa and the voltage U as an ordinate; with U being equal to Umax、u=Umin、t=TmaxAnd T ═ TminFour straight lines are drawn, dividing the first quadrant of the graph into three regions:
normal zone: the abscissa T is less than or equal to TminOr the ordinate U is greater than or equal to UmaxThe area of (a);
a fault area: the ordinate U is less than or equal to UminAnd the abscissa T is greater than or equal to TmaxThe area of (a);
fuzzy areas: removing the remaining areas of the normal area and the defective area;
the voltage endurance duration is defined as the time required from the occurrence of a fault to the time when the voltage rises back to 0.9p.u., and the voltage sag amplitude is defined as the average value of the voltage effective values in the voltage sag duration; the U ismaxDefined as the maximum value of the voltage withstand amplitude of the device, UminDefined as the voltage withstand of the deviceMinimum value of amplitude, TmaxDefined as the maximum value of the voltage withstand duration of the device, TminDefined as the minimum value of the device voltage withstand duration;
2) the severity of the voltage sag, Q, was calculated as followsij
Definition of UijIs the voltage sag amplitude, T, at the jth bus under the ith faultijFor the voltage withstand duration at the jth bus under a class i fault,
if U isij、TijFalls within the normal region of the graph, then Qij=0;
If U isij、TijFalls within the fault region of the graph, then Qij=1;
If U isij、TijFall into the fuzzy area of the coordinate graph, then
Figure FDA0001443520430000041
10. The method according to claim 9, wherein the U is a number of unitsmaxIs 0.9 p.u.; u shapeminThe value of (a) is 0.6 p.u.; t ismaxIs 600 ms; t isminIs 20 ms.
CN201711001724.8A 2017-10-24 2017-10-24 Optimal configuration method of dynamic reactive power compensation device Active CN107528323B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711001724.8A CN107528323B (en) 2017-10-24 2017-10-24 Optimal configuration method of dynamic reactive power compensation device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711001724.8A CN107528323B (en) 2017-10-24 2017-10-24 Optimal configuration method of dynamic reactive power compensation device

Publications (2)

Publication Number Publication Date
CN107528323A CN107528323A (en) 2017-12-29
CN107528323B true CN107528323B (en) 2020-08-04

Family

ID=60683919

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711001724.8A Active CN107528323B (en) 2017-10-24 2017-10-24 Optimal configuration method of dynamic reactive power compensation device

Country Status (1)

Country Link
CN (1) CN107528323B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109768581B (en) * 2018-12-28 2021-07-06 国网湖南省电力有限公司 Power grid voltage regulation and dynamic reactive support control method for energy storage power station
CN112834867B (en) * 2021-01-06 2022-08-02 南京工程学院 Optimized deployment method of wide-area synchronous intelligent sensor
CN113659581A (en) * 2021-07-20 2021-11-16 广西电网有限责任公司电力科学研究院 Voltage sag state series and parallel compensation equipment management method

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103795068A (en) * 2014-03-05 2014-05-14 广东电网公司电力调度控制中心 Optimal configuration method for high-voltage distribution network dynamic reactive power compensation equipment capacity
CN105846442A (en) * 2016-05-12 2016-08-10 广东电网有限责任公司电力科学研究院 Reactive power compensation optimizing configuration method and reactive power compensation optimizing configuration device
CN106058876A (en) * 2016-06-02 2016-10-26 国网山东省电力公司济南供电公司 Dynamic reactive planning site-selection analysis method and system considering transient voltage stability

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103457279A (en) * 2013-09-23 2013-12-18 广东电网公司电力调度控制中心 Determination method for partition D-STATCOM integrated optimization configuration scheme of large-scale power grid
CN104333019B (en) * 2014-11-27 2016-06-29 宁波华浙电力技术有限公司 10KV overhead transmission line reactive power compensator voltage optimization reconnaissance method
CN105914755B (en) * 2016-04-21 2018-07-24 广州供电局有限公司 Electrical power system dynamic reactive power Optimal Configuration Method and system
CN106532727A (en) * 2016-12-13 2017-03-22 国网北京市电力公司 Configuration method and device of installation points of dynamic reactive power compensator

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103795068A (en) * 2014-03-05 2014-05-14 广东电网公司电力调度控制中心 Optimal configuration method for high-voltage distribution network dynamic reactive power compensation equipment capacity
CN105846442A (en) * 2016-05-12 2016-08-10 广东电网有限责任公司电力科学研究院 Reactive power compensation optimizing configuration method and reactive power compensation optimizing configuration device
CN106058876A (en) * 2016-06-02 2016-10-26 国网山东省电力公司济南供电公司 Dynamic reactive planning site-selection analysis method and system considering transient voltage stability

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
考虑电网特性的动态无功优化配置方法;黄小庆等;《电力自动化设备》;20160930;第127页右栏第2段-第132页左栏第3段 *

Also Published As

Publication number Publication date
CN107528323A (en) 2017-12-29

Similar Documents

Publication Publication Date Title
CN108418255B (en) Planning method and system suitable for extra-high voltage direct current outgoing power grid containing high-permeability new energy
CN103414197B (en) A kind of reactive-load compensation equipment collocation method of large-scale wind power delivery system
CN103296685A (en) SVC (static var compensator) compensation strategy optimizing method
CN107528323B (en) Optimal configuration method of dynamic reactive power compensation device
CN107181253B (en) Power grid planning method based on power grid dynamic reliability probability index
CN110728065A (en) Voltage sag fault probability estimation and loss estimation method
CN102646155B (en) A kind of load modeling method being applicable to calculation of short-circuit current
CN116054224A (en) Quick evaluation method for commutation failure of multi-feed direct-current transmission system
Du et al. Hosting capacity assessment in distribution networks considering wind–photovoltaic–load temporal characteristics
CN111162541B (en) Dynamic partitioning method for voltage control of power system
CN110661268B (en) Dynamic reactive power compensation requirement determining method and system for new energy collection station
CN111626591A (en) Risk assessment method for flexible platform area
CN103400309A (en) Method for evaluating influence of ultrahigh-voltage access on stability of provincial power grid
Guo et al. Research on Dispatching and Control Strategy of UHVDC Blocking Fault in Sending-end Power Grid Based on K-means Clustering
CN109560561A (en) Active distribution network asymmetrical three-phase operation state emulation mode, system and storage medium
CN104143039A (en) System impedance calculating method
ZHANG et al. Analysis of voltage stability weak nodes location and optimization measures in DC multi-infeed receiving-end power grid
CN110445128B (en) Sensitivity-based real-time evaluation method for new energy consumption capability of regional power grid
CN112906172B (en) Onshore grid-connected point optimal selection method and system for large-scale offshore wind farm
Meng et al. Voltage stability research of receiving-end network based on real-time classification load model
Shan et al. Mechanism-Aware Distribution Network Reconfiguration with Transient Voltage Stability Enhancement
CN111600393B (en) Method for reducing voltage measurement data of transformer substation in different voltage classes and identification method
CN114243758B (en) Overcurrent protection self-adaptive setting method based on neural network model fitting
Ma et al. A real-time data quality analysis method and application for power systems based on balance calculation
Zhang et al. Voltage safety risk assessment of low and medium voltage distribution networks with high proportions of distributed photovoltaics

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