Technical field
[0001]
　The present invention relates to a system stabilizing controller and power system control system.
Background technique
[0002]
　Power system, many of the generator load, and a power transmission and distribution equipment and control devices. In the steady state, very small disturbance due to load fluctuations and the system switching is applied to the power system. When failure, such as strains of ground fault or short circuit, a large disturbance is applied to the power system. Cause these disturbances, power oscillation having a plurality of oscillation modes component occurs.
[0003]
　If this oscillation modes component is small, or when the attenuation speed, for system oscillations after a disturbance generator is maintained at a sufficiently small range, the power system is in a stable state. However or upset increases, if the attenuation persists or diverge not converge, or the generator is stopped as a result, a wide area power outage in or generated, there is a possibility that the stable operation becomes impossible.
[0004]
　With the expansion of the power system scale, power supply capacity increases, a remote underground, long distance of uneven distribution and transmission lines associated therewith, tends to heavy power flow of proceeds. Since effects of these changes are acting in a direction to lower the stability of the power system, is a concern expansion of cascading effect to large power drop out and a wide area power outage was triggered by the generator step-out during system fault .
[0005]
　Therefore constantly monitors the stability of the system, advance to prevent spread of the wide area power outage, it is necessary system stabilization control method according online to predict.
[0006]
　Recently, a synchronous measurement using GPS (Global Positioning System), are becoming popular phase detector PMU (Phasor Measurement Unit) is to measure the bus voltage phase angle of the power system in substantially real time. By using an on-line measured values, such as system bus voltage measured by the phase detector PMU, it is becoming possible to grasp the power oscillation phenomenon in real-time.
[0007]
　Above background of the original, by operators of power system has been proposed system stabilization control technology for stabilizing monitors the state of the system. For example, Patent Document 1, "plurality of control target power many generators connected a phenomenon becomes unstable to the power system of the electric power system generator group by system fault of the connected power system in the power system stabilizer for stabilizing by blocking a part of the generator at a high speed of at, the system fault occurred before a steady state electric power system of the control target power station at the generator unit and other main and two aircraft system model creating means for calculating a generator constants and the initial state of when equivalently expressed Part 2 machine system model with a constant period of a few minutes the generator group in power within the system, the system an electric quantity detection means for detecting a voltage and current and power of the control target plant failure occurs activation condition, enter the electric quantity detected by the electric quantity detection means the controlled plant after the system fault is removed the main power system from the bus And the external reactance estimating means for calculating the external reactance when viewed, the generator is calculated in two aircraft system model generating means constant and the initial state quantity and the external reactance estimating means external reactance calculated by said electric quantity detection as calculates the upset of the generator groups and other of the main power two equivalent generator the movement of the generator group equivalently represented within the system of the control target power station based on the electric quantity detected by the means operation result and the preset by electrically controlled priority of each generator are electrically controlled generator determining means for determining the electrical control generators required for regulation by the extended like area method, in the electronically controlled generator determining means there is described that the feature "in that a control means for disconnection from the power system selected generator.
CITATION
Patent Literature
[0008]
Patent Document 1: Laid-Open Patent Publication No. 09-046908
Summary of the Invention
Problems that the Invention is to Solve
[0009]
　In recent years, renewable energy represented by solar power and wind power generation is being introduced a large amount to the system. These renewable energy with a steep and pre assumed difficult generator output fluctuation, as a result, there is a possibility that the system characteristic is greatly changed. For example, when it exceeds a certain threshold wind speed, automatically stops the wind power in order to maintain the safety of the device, operation of the cut-out for the power output to zero is generated. Such events have high occurrence frequency in comparison with the short-ground fault generated in the conventional system, when the generator is installed densely in areas adjacent the simultaneous in multiple generators It is expected to be likely to occur.
[0010]
　As the introduction of renewable energy increases, the impact on the system due to sharp output variations events as described above is increased, oscillation modes of the wide area is a concern to destabilize. To cope with this, the stabilization control is necessary based on the high-speed, high-precision determining the stability in the intermediate region of about ten seconds from a few seconds immediately after the output variation events.
[0011]
　Furthermore, because of the fast stabilization control, for each unstable oscillation modes that occur, it is necessary to select a control method stabilization is possible online. In this regard, in the technique described in Patent Document 1, it is difficult to realize a high speed and highly accurate stabilization control.
[0012]
　An object of the present invention from the above is to provide a fast and accurate stabilization control can be realized system stabilizing controller and power system control system.
Means for Solving the Problems
[0013]
　In the present invention in order to solve the above problems, on the basis of the measurement information in multiple measurement point of the power in the system, a system stabilizing control apparatus for controlling a controlled device installed in the power within the system, with the power system a facility information database for storing information of the system equipment which is the system, the control priority determining section for determining the control priority of the control device based on the information of the system equipment, control based on the measurement information and control priority in the plurality measurement point a control target determination section for determining a target, and destabilization time calculation unit for calculating a destabilization time of the electric power system from the measurement information in multiple measurement point, the control of the control device on the basis of information destabilization time and system equipment characterized in that it and a control time determination section for determining a time.
[0014]
　In the present invention also based on the measurement information at a plurality measurement point of the power within the system, a system stabilizing control apparatus for controlling a controlled device installed in the power within the system, the control equipment to be interconnection to the power system including a facility information database storing facility information of the system equipment, a control priority determining section for determining the control priority of the control device based on the facility information, the power system oscillations which is detected from the measurement information upset measuring points and control equipment distance and, a control target determination section for determining a control apparatus to be controlled according to the priority, and the destabilization time calculation unit for calculating a destabilization time of the power system oscillations which is detected from the measurement information, destabilization time and a control time determination section for determining the control time based on the facility information, and executes the control of the control time determined for the determined control equipment.
[0015]
　The present invention, to obtain a measurement information in multiple measurement points on the system stabilizing control device via the first communication means of the power within the system, the second communication to the control device in the system equipment constituting the electric power system a power system control system that provides control commands from the system stabilizing controller via the unit, the system stabilizing control device includes a facility information database for storing information of the system equipment is interconnection to the power grid, grid facilities a control priority determining unit for determining a control priority based control equipment information, and the control target determination section for determining a control target based on the measurement information and control priority in the plurality measurement point, the measurement information at multiple measurement point determining a destabilizing time calculation unit for calculating a destabilization time of the power system, and a control time determination section for determining the control time based on the destabilization time and facility information, and the determined control equipment Characterized in that it comprises an output unit for providing a control command including control time through the second communication means.
Effect of the invention
[0016]
　According to the present invention, it can be a high speed and highly accurate stabilization control can be realized system stabilizing controller and power system control system.
BRIEF DESCRIPTION OF THE DRAWINGS
[0017]
[1] shows an example functional configuration of the grid stability control apparatus according to an embodiment of the present invention.
Diagram illustrating a configuration example of a power system control system according to an embodiment of the present invention; FIG.
[3] the system flow diagram illustrating the stabilization control determining process.
[4] shows the stored contents of the detected motion information table for storing the fluctuation component at each measurement point obtained by shaking detection unit 12 FIG.
[5] multiple measurement point extraction, as well as flow chart showing the upset influence estimation processing means upset influence calculation.
[6] shows the detailed contents of the upset influence point extraction process (process step S11).
It shows an example of FIG. 7 frequency motion information table TB2.
[8] a flowchart showing in detail the distance calculation processing between measurement point in the process step S112.
[9] shows a node information table TB3 facility information database DB.
[10] illustrates the transmission line information table TB4 facility information database DB.
[11] shows a measuring equipment information table TB5 facility information database DB.
It shows an example of FIG. 12 process resulting area upset in the process of step S113 information table TB6.
[13] shows a concrete flow of upset influence calculation process.
[14] shows a concept of upset influence calculation process.
Shows an example of FIG. 15 past motion information table TB7.
It shows a flow chart illustrating the FIG. 16 past motion information calculating process.
[17] shows the flow of the control object determination process.
[18] illustrates the flow of unstable time calculation process.
[19] illustrates the flow of control time determination process.
[20] shows an example of a display screen.
Diagram illustrating a configuration example of FIG. 21 Example 2 system stabilizing control device according to the.
[22] illustrates a flow chart of system stabilization control determining process according to the second embodiment.
Diagram illustrating a configuration example of FIG. 23 Example 3 system stabilizing control device according to the.
Diagram illustrating a configuration example of a power system control system according to [24] Example 3.
FIG 25A] shows a processing flow in the system stabilizing control device 200 according to the third embodiment.
FIG 25B] shows a process flow in the entire system stabilizing control device 300 according to the third embodiment.
DESCRIPTION OF THE INVENTION
[0018]
　Examples of the present invention will be described below in detail with reference to the drawings. Incidentally, only examples below merely exemplary, and are not intended to intended that the invention itself is not limited to the following specific details.
Example 1
[0019]
　Figure 1 is a diagram showing a functional configuration example of a system stabilizing control device according to an embodiment of the present invention.
[0020]
　System stabilizing control device 200 according to this embodiment, the system stabilizing control decision unit 10 includes an information storage unit 30. Among system stabilizing control decision unit 10 is connected to the grid measurement unit 20. System measuring unit 20 has a plurality of points of measurement information for measuring the measuring device 21 (21a, 21b, ··· 21n) is composed of. Information storage unit 30 stores the facility information in the facility information database DB. In the following description, when there is no need to distinguish between the elements by alphabet code, it is possible to omit the alphabet code.
[0021]
　System stabilizing control decision unit 10 includes a measurement information collecting unit 11, and the upset detector 12, a control priority determining unit 13, a control target determination unit 14, a calculating unit 15 destabilization time, control time determination unit including the 16.
[0022]
　Measurement information collecting unit 11, the measurement information of a plurality of points measured in the measuring device 21 of the system measuring unit 20 (system condition amount), received via the information communication network, it collects.
[0023]
　Upset detection unit 12 based on the measurement information collected by the measurement information collecting unit 11, with respect to the measurement information of a specific type (measurement type), for each measurement point, the fluctuation component is a frequency component of upset power to calculate the detection motion information shown.
[0024]
　Control priority determining unit 13, based on the facility information stored in the facility information database DB to calculate the control performance of the plurality of equipment, to determine the control priority for each equipment based on the control performance.
[0025]
　Control target determination unit 14 determines a control object based on the detection motion information calculated in control priority and the upset detector 12 determined in the control priority determining unit 13.
[0026]
　Destabilization time calculating unit 15 based on the detection motion information calculated in the upset detector 12, the system calculates the destabilization time the limit time to destabilize.
[0027]
　Control time determination unit 16 determines the control target based on the destabilization time calculation unit calculated in facility information and destabilization time calculating unit.
[0028]
　Note that system stabilizing control device 200 is equipped with an input operation device such as a monitor 202, a keyboard, a user is capable verify appropriate internal processing contents, and those capable of performing input processing as appropriate required configuration ing.
[0029]
　System measuring unit 20 measures the measurement information in the measurement point in the system, and transmits the measurement information collecting unit 11. Measurement information includes effective power flow, reactive power flow, the system voltage, system current, at least one type in the voltage phase. System measuring unit 20 includes husband plurality of measurement point s arranged measuring device 21a, 21b, · · ·, a 21n. Measuring device 21, for example, PMU is (Phasor Measurement Unit) or the like, performs regular measurements concerning power, and transmits the measurement information of the time series.
[0030]
　The equipment information database DB of the information storage unit 30 stores the facility information indicating the connection information of the system elements that interconnection to the power system. The facility information, the generator being interconnection in the vicinity of the measurement point, the load, the system elements connected to each node in the system such as compensator device information showing the characteristics of the (electrical equipment) (node ​​information D1) , one of the information indicating the impedance or the like of the transmission line that connects each node (transmission line information D2), and information indicating a positional relationship between the latitude and longitude and topologies such as a plurality of measuring points (measuring equipment information D3) It is included. Thus, it is possible to calculate the electrical distance between the measurement point, it is possible to classify the measurement point in the electrical distance. Note that performs a node information D1, the power transmission line information D2, the measuring equipment information D3, 9, 10, separately in detail described by way of example in FIG. 11.
[0031]
　Figure 2 is a diagram illustrating a configuration example of a power system control system according to an embodiment of the present invention.
[0032]
　As shown in FIG. 2, power system control system according to an embodiment of the present invention, a generator G, a node (bus) N (N1, N2, N3), and a substation SS, a compensator device C, comprising a power load LD, and external power system PS, the measuring device 21 (21a, 21b, ··· 21e) and, an information communication network 108, the system stabilizing control device 200. Hereinafter, a generator G, may be referred to as a substation SS, a compensator device C, a power load LD, the measuring device 21 (21a, 21b, ··· 21e) and own line power system including.
[0033]
　Generator G is a generator producing power generation, thermal power generation, hydroelectric power generation, nuclear power, solar power, wind power, biomass power generation, the generator resulting in generating capacity by any of the power generation techniques, including tidal power generation is there.
[0034]
　Generator Ga is thermal power to be installed in the high voltage side of the substation SS in the power system, hydroelectric, a large generator, including nuclear power, including the power generation amount through measuring apparatus 21a and the information communication network 108 sending a system state quantity to the grid stabilization control device 200. Further, the generator Ga receives the control instruction information transmitted from the system stabilizing control device 200 through the measuring device 21a and the information communication network 108, to change the system state variables including the power generation amount in accordance with the control command information.
[0035]
　Generator Gb is photovoltaic installed on the low voltage side of the substation SS in the power grid, wind power generation, a small generator in, including cogeneration, the power generation amount through measuring device 21b and the information communication network 108 the system state quantity comprising transmitting to the system stabilizing control device 200.
[0036]
　Substation SS is installed between the transmission lines of the power within the system, to change the power voltage of the generator Ga is large generator is power than the high voltage side is installed, the power load LD is installed transmitting power to the low voltage side has. The substation SS, power capacitor, compensator device C, such as shunt reactor is connected.
[0037]
　Compensator device C, by changing the reactive power in the power within the system, a device for controlling the voltage distribution of power within the system, power capacitors, shunt reactors, STATCOM (Static Synchronous Compensator: self-commutated reactive power compensator ), SVC (Static Var compensator: static var compensator), and the like. Some of the compensator device C receives the control instruction information transmitted from the system stabilizing control device 200 through the measuring device 21c and the information communication network 108, changes the system state quantity including the power generation amount in accordance with the control command information make.
[0038]
　Power load LD is a facility, including electric motors and lighting equipment that consumes power, representing home, factory, a facility such as a building.
[0039]
　External power system PS is the external power system it is impossible to control from the system stabilizing control device 200 is connected to the own line by tie line.
[0040]
　Measuring device 21 (21a, 21b, ··· 21e) is power generation in the power generator Ga, adjusting the phase of the phase modifying equipment C, power flow value in the transmission lines, a sensor for measuring the system state quantities such as a voltage value, wherein, transmits measured system condition quantity, through the information communication network 108 to the system stabilizing control device 200.
[0041]
　Information communication network 108 is a network capable of transmitting data bidirectionally. Information communication network 108, for example, a wired network or a wireless network, or a combination thereof. Information Communication network 108 may be a so-called Internet, may be a network of dedicated lines.
[0042]
　System stabilizing control device 200 is a device for implementing the system stabilizing control function shown in FIG. System stabilizing control device 200, the measuring device 21 (21a, 21b, ··· 21e) of system condition amount measured in, received through the information communication network 108. Further, the system stabilizing control device 200, the calculated control command information using a systematic state quantity and information accumulated in the interior of the transmitted line, and transmits to the measurement device 21 via the information communication network 108.
[0043]
　System stabilizing control device 200 includes, as its internal structure, the CPU (Central Processing Unit) 201, a display unit 202, a communication unit 203, an input unit 204, a memory 205, a storage device 206, but the bus line 211 It is connected to the. CPU201 executes the calculation program stored in the storage unit 206, generates the calculation and the control signal of the system state, performs such generation of the image data, and displays the image data on the display device 202. Memory 205, image data for display, a memory that temporarily stores a calculation result data of the system state, for example, RAM (Random Access Memory) and the like. Communication means 203, power flow value from the measuring device 21 through the information communication network 108, obtains the system state quantities such as a voltage value. Measurement information collecting unit 11, using the communication unit 203 receives the measurement information of a time series which is measured by a plurality of measuring devices 21.
[0044]
　System stabilizing control device 200 of the user (administrator), the parameters of the various threshold settings or change through a predetermined interface of the input unit 204 can appropriately set the operation of the system stabilizing control device 200 of the own system. The user selects the type of data to be confirmed through a predetermined interface of the input unit 204 can be displayed on the display device 202.
[0045]
　Storage device 206 stores various programs and data. Storage device 206, for example, a HDD (Hard Disk Drive) or a flash memory. Storage device 206, for example, to hold the programs and data, etc. can realize various functions which will be described later. Programs and data stored in the storage unit 206 is executed by being read out to the CPU201 as needed. The storage device 206 enables information storage unit 30 stores various databases (DB).
[0046]
　Figure 3 is a flow diagram showing a system stabilization control determining process. Each processing determination flow in FIG. 3 corresponds to each part of the system stabilizing control decision unit 10 of FIG. Specifically, the process step S10 is upset detector 12 of FIG. 1 in FIG. 3, the processing step S20 is the control priority determining unit, the processing step S30 the control target determination unit 14, the processing step S40 is unstable time calculator 15, the processing step S50 corresponds to the control time determining unit 16. Since Figures 1 and 3 are in the above relationship, in the following description is made with contrasting them.
[0047]
　First, the upset detector 12 of FIG. 1 (the process step S10 in FIG. 3), using the measurement information of the plurality of points measured by the strain measuring unit 20 calculates the detected motion information.
[0048]
　Here, measurement information of a plurality of points measured by the strain measuring unit 20 comprises an effective power flow, reactive power flow, the system voltage, system current, the information of at least one measurement type of the voltage phase. System stabilizing control device 200 may perform the system stabilizing control determination process for a particular one of the measurement types, may be performed system stabilizing control decision process for a plurality of measurement types. The calculation of the detection motion information in the present process, the algorithm used for frequency analysis, i.e., Fourier analysis, wavelet analysis, Puroni method, Matrix Pencil method, HTLS (Hankel Total Least Squares) method at least one of the algorithms as example One to the application.
[0049]
　The following describes motion information calculation algorithm by Puroni method. The motion information calculation algorithm by Puroni method executes the following (1) to (12). Note Upon performing these equations, the number of data N, mode number n of analysis window width is assumed to be set in advance.
[0050]
　The upset detector 12, first input data y is the measurement information of each point (k) (k = 0,1, :, N-1) using the following (1) to (4) run, matrix = a (a 1 , a 2 ,:, a n ) T is calculated. However, the "T" is transposed matrix of superscript, superscript "-1" denotes the inverse matrix.
[0051]
　Here, (1) is a formula for calculating the autoregressive model coefficients vector a, (2) formula (1) defined type autoregressive model coefficients vector a, (3) equation (1) X defining equation of (4) is a defining equation of equation (1) Y. In yet (3), k is time, y is the input measured values ​​(values ​​from when the measuring device 21 for example PMU), n is preconfigured mode order, N is the pre-configured analysis window width (for signal processing, the meaning of the data size Baffuaringu the measure) is.
[0052]
[Number 1]

[0053]
[Number 2]

[0054]
[Number 3]

[0055]
[Formula 4]

[0056]
　In next upset detector 12, by executing the equation (5) to (7) below, to calculate the eigenvalue lambda. Here (5) is (1) autoregressive model coefficients calculated by the formula (a 1 , a 2 , · ·, a n characteristic equation for calculating the eigenvalues using a), (6) the solution Z and eigenvalues relational expression lambda, (7) equation is attenuated eigenvalue lambda factor sigma, a relational expression similar frequency omega.
[0057]
　Specifically first, the (1) to (4) autoregressive model coefficients calculated by formula (a 1 , a 2 , · ·, a n ) to be substituted into the following equation (5). (5) formula, Z i is an n-order polynomial to the solution (characteristic equation), n pieces of Z i , sigma i , f i are obtained as a solution, respectively.
[0058]
　Thereafter, the upset detector 12, (5) the solution Z of formula characteristic equation, (6) and (7) Z is represented by the formula of the form i substitutes, eigenvalue lambda i , attenuation coefficient (attenuation factor ) sigma i [1 / sec], upset frequency f i [Hz] = omega i to calculate a / 2 [pi. Here Δt represents the sampling period.
[0059]
[Formula 5]

[0060]
[Number 6]

[0061]
[Number 7]

[0062]
　In next upset detector 12, the input data y (k) (k = 0,1 , ··, N-1), and the calculated Z i (i = 1, 2, · ·, n) using, the following (8) to (11), matrix = b (B 1 , B 2 , · ·, B n ) T is calculated. Where (8) calculates the complex amplitude vector b formula (9) is (8) defining formula of vector b, (10) the definition formula of the solution Z of (8), ( 11) is a defining equation of equation (8) of w. Note (8) elements in formula is the complex amplitude B of each oscillation modes, and y is the measurement value in (11).
[0063]
[Number 8]

[0064]
[Number 9]

[0065]
[Formula 10]

[0066]
[Number 11]

[0067]
　Finally, in shaking detection unit 12, (8) to (11) the complex amplitude B calculated from equation (B 1 , B 2 ,:, B n ) absolute value R from the formula (12) i , polarization angle Shita I seek. At this time, the absolute value R i is the amplitude, polarization angle theta i means the initial phase. Note equation (12), the complex amplitude B and the amplitude R, which is a relational expression of the phase theta.
[0068]
[Number 12]

[0069]
　By the above processing, upset frequency, damping coefficients, amplitude, detecting motion information including the initial phase is obtained. Detecting motion information is by including upset frequency, damping coefficients, and one of the amplitude, may represent fluctuation component at each measurement point. Further, by using such a detection motion information, it is possible to predict the future fluctuation component.
[0070]
　Figure 4 is a view showing stored contents of the detected motion information table TB1 storing the fluctuation component at each measurement point obtained by the above equation (1) (12) of the process.
[0071]
　Detection motion information table TB1 in FIG. 4 has an entry for each detection motion information. The entry corresponding to one detection motion information, and 100 (upset name in the drawing) detects upset identification information for identifying the detected motion information, motion characteristics 101 showing characteristics of upset represented by the detection motion information including the door. Detection upset identification information 100 includes a name 100a of the measuring point corresponding to the input data, and a detection motion information number (No.) 100b indicating the mode of the detection motion information.
[0072]
　If the detected motion information using a Puroni method described above is calculated, per type of the measurement information of one measurement point, detecting motion information of the same number as the mode number n is calculated. Accordingly, in each point of the p number of points (p ≧ 1), when using the measurement information of q type (q ≧ 1), n ​​× p × q pieces of sensing motion information is calculated. Upset characteristics 101, upset frequency 102 calculated damping coefficient 103 includes an amplitude L, initial phase 105.
[0073]
　In this way, the upset detector 12 (process step S10 in FIG. 3), but calculates the detected motion information, which can be utilized as a further processed information as follows. An example of processing information is upset influence is shown by a dotted line as S10A processing step of determining a perturbation influence in FIG.
[0074]
　In the processing step S10A upset detector 12, a detecting motion information obtained in the process step S10, based on the data of the equipment information database DB stored in the information storage unit 30, associated with each fluctuation component to be judged identical to the plurality extracting the measurement point, and calculates the upset influence representative of the magnitude of the influence of the fluctuation component. By determining the system stabilizing control using the agitation impact, since a wide large influence range is upset, it can be selected preferentially control for stabilizing the high impact upset to Unstable lineages.
[0075]
　A series of processing of the processing step S10A upset detector 12 (more measurement point extraction, as well as upset influence calculation) will be described in detail with reference to FIG. 13 from FIG.
[0076]
　First, FIG. 5 is a flowchart showing the upset influence estimation process that represents a series of processing of the processing step S10A (multiple measurement point extraction, as well as upset influence calculation). The flowchart of FIG. 5, is executed in the upset detector 12 of FIG. 1, which is obtained by embodying the process of the processing step S10A of FIG. In the first processing step S11 in FIG. 5 performs upset influence point extracting process of extracting a plurality measurement point associated with each fluctuation component to be determined the same as. In the next process step S12, it calculates the upset influence representative of the magnitude of the influence of the fluctuation component. The detailed processing contents of the processing step S11 is shown in the process flowchart of FIG. 6, the detailed processing contents of the processing step S12 shown in the flowchart of FIG.
[0077]
　Process flowchart of FIG. 6 shows the detailed contents of the upset influence point extraction process (process step S11). In the first processing step S111 in this case, to classify a plurality of detection motion information calculated by the upset detector 12 in upset frequency, it stores the result in the information storage section 30 as the frequency motion information table TB2. Note that although the plurality of detecting motion information calculated by the upset detection unit 12 utilizes those already stored in the information storage unit 30 as the detecting motion information table TB1 described above.
[0078]
　Figure 7 shows an example of a frequency motion information table TB2. The table TB2 of Fig. 7 as is apparent from the comparison with the table TB1 of Fig. 4, the table TB2 is obtained by replacing the described sequence a table TB1 in Fig. 4 in terms of upset frequency 102.
[0079]
　To describe in more detail the frequency motion information table TB2 in FIG. 7, which includes an entry for each frequency motion information. The entry corresponding to one frequency motion information, having at least upset the frequency 102 of the fluctuation component corresponding to the frequency motion information, and the detected upset identification information 100 (sway name in the drawing), the upset frequency close to the upset frequency and a sway characteristics 101 showing characteristics of one detection motion information. Detection upset identification information 100, detected upset identification information 100a of the detection motion information is the same as 100b. Upset characteristics 101 includes a damping coefficient 103, and amplitude 104, and an initial phase 105.
[0080]
　In the example of this frequency motion information table TB2, the frequency motion information having upset frequency close to 0.18, No. in the measurement point PMU1 And one of the detected motion information, No. in the measurement point PMU2 It is included and one of the detected motion information. Similarly, the frequency motion information having upset frequency close to 0.65, No. in the measurement point PMU1 Detecting motion information of 2 are included. Hereinafter, it may be referred to as measurement points included in one frequency motion information and the measurement point group candidates.
[0081]
　Note Upon this classification, the upset detector 12, using the frequency width shown the difference in upset frequency defined in advance, the detection motion information difference upset frequency is within the frequency range, the frequency motion information of the same operating frequency it may be classified. In the example of detecting motion information table TB1 of the foregoing, No. of measurement point PMU2 Upset frequency of the detected motion information of 1 is 0.19, No. in the measurement point PMU1 1 upset frequency of the detected motion information is a 0.18, which are classified into the same frequency motion information. By perturbation frequency classifies detected motion information within the frequency range on the same frequency motion information, each detection motion information, in a case that contains the errors due to noise, the exact distribution of the system of the upset It can be estimated.
[0082]
　Here, the frequency width delimiting the upset frequencies to be included in the same class may be determined based on the number of measurement information used as input data. As an example, the frequency width Δf of the sway frequency included in the same frequency motion information to (13). (13) is an equation representing a method of determining the frequency width Δf based on the sampling theorem. Here, the number of input data of the measurement information using the N in the calculation of the detection motion information, a constant specified sampling frequency of the measurement information? Fs, alpha in advance. In general, the frequency resolution in the frequency analysis is high when the sampling frequency of the input data is smaller, since higher when the number of data is large, the equation (13), the frequency width Δf of the sway frequency included in the same frequency motion information properly it can be selected.
[0083]
[Formula 13]

[0084]
　Then, shaking detection unit 12 performs measurement point distance calculation process for calculating the distance between the measurement point on the basis of the data of the equipment information database DB stored in the information storage unit 30 in the process step S112 in FIG. 6.
[0085]
　Figure 8 is a flow chart illustrating further details of the distance calculation processing between measurement point in the process step S112. In the first processing step S114 in FIG. 8, in order from the measurement point of the obtained measurement information by measuring the information collection unit 11, selects the two measurement points.
[0086]
　In the next processing step S115, the measurement point distance between the two measurement point selected is calculated based on the stored data to the facility information database DB stored in the information storage unit 30. Here and in reference to the equipment information database DB of the data includes a node information D1 as described above, the transmission line information D2, and a measuring equipment information D3.
[0087]
　9, among the data of the equipment information database DB, which is a diagram illustrating an example of a node information table TB3 for storing node information D1.
[0088]
　Node information table TB3 for storing node information D1, for each node in the power system, a node name 200 for identifying the node, for example, the position information 201 indicating the position of the node in terms of longitude and latitude, the node generator which is connected to, including power load, and a characteristic 202 of the compensator device. In the example of this figure, the generator 202a is interconnection to the node A and node B. The node A has thermal generator interconnection rated capacity name is 100 called G1, wind power generator is interconnection to the node B is 200 rated capacity named G2.
[0089]
　In the example of this figure, the power load 202b are interconnection to the node AA and BB. The node AA and housing the rated capacity of the named L1 is 1000 and interconnection, factory rated capacity of the named L2 is 2000 is interconnection to the node BB.
[0090]
　In the example of this figure, nodes AAA and BBB node compensator device 202C of is interconnection. Rated capacity named D1 to the node AAA is 10 SC (Static Condenser: capacitor power) has to interconnection, the node BBB rated capacity named D2 is 20 ShR (Shunt Reactor: min road reactor) are interconnection.
[0091]
　10, among the data in the equipment information database DB, which is a diagram showing an example of a transmission line information table TB4 for storing power transmission line information D2.
[0092]
　Transmission line information D2, for each transmission line in a power system, a power line name 300 for identifying a transmission line, both ends node 301 that indicates the nodes at both ends of the transmission line, the positive-phase resistance 302a of the transmission line, positive-phase reactance 302b, and a positive phase capacitance 302C. In the illustrated example, the transmission line a is present between the node A and the node B, the transmission line positive-phase resistance 0.01 a, positive-phase reactance 0.2, a positive-phase capacitance 0.1 is there. The transmission line b is present between node B and node C, the positive-phase resistance 0.02 of the transmission line b, the positive-phase reactance 0.5, positive-phase capacitance is 0.2.
[0093]
　11, among the data of the equipment information database DB, which is a diagram showing an example of a measurement facility information table TB5 for storing measuring equipment information D3.
[0094]
　Measuring equipment information table TB5 includes measuring equipment information, for each point of measurement within a power system, two measurement point 400a, identification information 400b, and a measurement point distance 403 between them. Here measurement point and indicates the point where the measuring device 21 is connected. In the example of table TB5 in this figure, measurement point 400a, 400b identification information includes measurement point name 401a indicating the measuring device 21, and 401b, the node name 402a included in the node information, and 402b. Measurement point of measurement point name PMU1 is connected to the node A, measurement point of measurement point name PMU2 is connected to node B.
[0095]
　Also, measurement point the distance between the electrical distance between the two measurement points, two generators which are interconnection in the vicinity of the measurement point, are calculated based on the system elements such as power load, phase modifiers equipment is that value. For example, the distance D between the measurement point is calculated by the following equation (14). Note (14) is a definition equation of the measurement point distance D.
[0096]
[Number 14]

[0097]
　Here, D is showing an electrical distance which is proportional to the impedance between the two measurement points. n is from one of the two measurement points, indicating the generator to be interconnection within the electrical distance specified by a certain threshold, the power load, the number of strains elements such as tone-phase equipment. The C i from one of the two measurement points, indicating the rated capacity of the i-th line elements interconnection within the electrical distance prescribed by certain threshold. Arufa, Beta I denotes the coefficient. By using the electric distance between the two measurement points, it is possible to determine whether the fluctuation component of the two measurement points are the same.
[0098]
　Figure the indicated data structure of these facilities information database DB is an example, the information storage unit 30 may store the more detailed facility information. For example, the node information D1 as the information of the node thermal generator to interconnection may include the power generation mobility characteristics such as the governor constants. The node information D1 the generator by natural energy such as wind power machine as the information of the node that interconnection may comprise a stochastic generation amount variation characteristic calculated by statistical analysis using the past history.
[0099]
　Further, in the calculation of the measurement point distance D in the measuring facility information D3, upset detection unit 12, the coefficient β of the equation (1) i and may be determined based on the power generation mobility characteristics described above. Further, upset detection unit 12, the impedance and power generators, power load, without using the information of the system elements, such as phase modifiers equipment, using the latitude and longitude information contained in the node information, and calculates the distance between the measurement point , it may be the value as measurement point distance D.
[0100]
　Upset detector 12 in the next processing step S116 in FIG. 8 with respect measurement point of the obtained measurement information by measuring the information collection unit 11 determines whether the calculated distance between the measurement point for all combinations of the two points . If it is determined that the calculated distance between the measurement point for all combinations of two points (Yes), upset detector 12 terminates the distance calculation processing between measurement point. If there is a combination of the measurement point is not calculated distance measurement point (No), upset detection unit 12, for another 2 measurement point processing step S114, and repeats the process step S115. By the above process, it calculates the distance between the measurement point for the combination of all the second measurement point.
[0101]
　Return to the upset impact point extraction process of FIG. 6, to resume the explanation. The process flow of FIG. 6, after the measurement point distance calculation (processing step S112), upset detector 12, a frequency motion information obtained in the processing step S111 are classified by the distance between the measurement point, regions upset result storing in the information storage section 30 as information (processing step S113). At this time, agitation detection unit 12 for each frequency motion information (upset frequency), using the distance between the measurement point, Ward's method, by classifying the frequency motion information using a clustering method represented by the k-means method it may be.
[0102]
　Figure 12 is a diagram showing an example of a region motion information table TB6 obtained by processing of the processing step S113. As is clear from looking at the area motion information table TB6 in FIG 12, which is the frequency motion information table TB2 in FIG. 7, in consideration of the measurement point name 401 measuring equipment information table TB5 in Figure 11, summarizes the upset number 500 those were.
[0103]
　It will be described in detail area motion information table TB6 in FIG. Area motion information table TB6 contains an entry for each classified by the distance between the measurement point region motion information. The entry corresponding to one region motion information, and upset number 500 indicating a fluctuation component represented by the region motion information (No.), The upset frequency 102 of the frequency motion information corresponding to the area motion information, the area including the measurement point name 102 that indicates the measurement point of frequency motion information corresponding to motion information, the upset characteristics 101 showing characteristics of the fluctuation component. Upset characteristics 101 includes similarly to upset characteristics of the frequency motion information (FIG. 7), and the attenuation coefficient 103, and amplitude 104, and an initial phase 105.
[0104]
　In the example of this figure, the detection motion information of upset frequency 0.18 observed at measurement point PMU1, and is also detected motion information of the upset frequency 0.18 observed at measurement point PMU3, in different areas motion information It has been classified. This measurement point distance between the measurement point PMU1 and measurement point PMU3 is, than the distance between other measurement point, because relatively large. Thus, the frequency motion information is classified by using the storage information of the facility information database DB.
[0105]
　If the difference between the upset frequency fluctuation component detected by the plurality of measurement point is within a predetermined range, based on the storage information of the facility information database DB, between the plurality of measurement points electrical distance by a plurality of It classifies measurement point, by extracting the measurement point belonging to each category as a measurement point group can identify areas where the fluctuation component is present.
[0106]
　In upsetting influence estimation process in FIG. 5, after the processing step S12, processing steps S11, performs upset influence calculation process for calculating the upset impact with area motion information of FIG. 12. Figure 13 shows a concrete flow of upset influence calculation process.
[0107]
　First, in the process step S121 of upset detector 12, each organized regions upset information as shown in FIG. 12, it extracts a measurement point included in the region motion information. For example, when explained in the example of the region motion information table TB6 in FIG 12, upset detection unit 12, measurement refers to the "1" for the first to upset No. 500 was detected upset frequency 102 the fluctuation component of the "0.18" as a point 401, "PMU1", to extract, such as "PMU3". Similarly, upset detector 12 refers to the "2" for upset No. 500, upset frequency 102 extracts the like "PMU1" as measurement point 401 detects the fluctuation component of "0.65". Similarly upset detector 12 refers to the "3" for the upset No. 500, upset frequency 102 extracts the like "PMU2" as measurement point 401 detects the fluctuation component of "0.18".
[0108]
　Next, in the process step S122 of upset detector 12, the upset characteristics corresponding to the measurement point extracted in the processing step S114 in FIG. 8, calculates the upset influence of each region motion information. By calculating the upset impact may represent the magnitude of the effect of each fluctuation component is given to the stability of the power system.
[0109]
　With reference to FIG. 14, illustrating the concept of upset influence calculation process. In Figure 14, the horizontal axis represents time and the vertical axis measurement information x ijk and upset influence J obtained from the measurement information ijk shows. In this example, the power variation shows a case in which the damped oscillation with time. It will be described in more detail Figure 14.
[0110]
　In FIG. 14, J ijk was upset impact area motion information j and measurement type k at measurement point i, x ijk and measurement information of the measurement type k with area motion information j in measurement point i, t1, t2 (t1 t2> t1, upset impact J past region motion information from the analysis reference time ijk is calculated. Similarly, setting t2> 0> t1, analysis upset impact J of the current region motion information including reference time ijk is calculated, t2> t1> 0 and If set, future area motion information from the analysis reference time upset impact J of ijk is calculated.
[0114]
　For example, t2> by t1> 0 and sets, to predict the time variation of the future fluctuation component in a certain measurement point, agitation impact J based on the predicted fluctuation component ijk is calculated and its upset influence J ijk , it is possible to determine whether the measurement point in the future, become unstable. Upset influence J ijk amplitude B in ijk and damping coefficient sigma ijk By using future perturbations influence J ijk can be predicted.
[0115]
　(15) upset influence J by the definition of formula ijk by calculating the amplitude B ijk is large and the attenuation coefficient sigma ijk is poor (attenuation coefficient sigma ijk is large) upset influence region motion information increases. Here, instead of upsetting degree of impact of (15), the amplitude B ijk or attenuation coefficient sigma ijk may define the value of itself as upset influence.
[0116]
　Further, in the process of shaking detection unit 12 refers to the past region motion information, from the past region motion information, select a past region motion information corresponding to the status of the measurement point group in the current area motion information and, upset impact J considering the duration of the fluctuation component of the past area motion information that has been selected ijk may be calculated.
[0117]
　For example, motion detection unit 12 in the past region motion information, if the period of fluctuation component having the same upset frequency exists continuously is period threshold value greater than or equal to the predetermined upset influence J ijk set large If period fluctuation component having the same upset frequency exists continuously is less than the period threshold value, upset influence J ijk may be set small. Thus upset influence J ijk by establishing, to eliminate sway due to the influence of noise included in the momentarily system fluctuation, the fluctuation component present constantly can be monitored preferentially.
[0118]
　Figure 15 is a diagram showing an example of the past motion information table TB7. Past motion information table TB7 of Figure 15 As is apparent from the comparison with the area motion information table TB6 in FIG 12, the past motion information table TB7 are those organized further upset impact and case information details imparting to it is.
[0119]
　For the past motion information table TB7, it will be described in more detail. Past motion information table TB7 has an entry for each case. Entry one case includes the case details 600 illustrating examples, at least one region motion information corresponding to the case. Case details 600, time 601, date or day of week 602, month or season 603 includes an event type 604, and the like. Area motion information is upset number 500, upset frequency 102, upset influence 605, the measurement point 401, the attenuation coefficients 103, the amplitude 104, including the initial phase 105.
[0120]
　Upset detector 12 classifies the historical area motion information for each case, the information obtained by statistically processing, is stored as the past motion information table TB7 in the information storage unit 30. Upset detection unit 12, the time of measurement information used in the calculation of the detected motion information, date, day of the week, month, season, based on the information of the event type, to set the example.
[0121]
　Upset detector 12, by setting the case with the time 601, it is possible to reflect the effect of a variation in tidal conditions in the system of day the statistical processing area motion information. Change of tide situation in this case, are due to the life behavior in the day, with an example of the working hours, which are separated by morning and evening attendance clock-out time. In the example of FIG. 15, Case No. In 1 8:00 to 10:00 is, Case No. In the two-12:00 to 14:00 is set.
[0122]
　Upset detector 12, by setting the case with the date or day of the week 602, it is possible to reflect the effect of a variation in tidal conditions in the system in a week statistical processing region motion information. Change of tide situation in this case, are due to the life behavior of each day of the week to be an example weekdays and holidays. In the example of this figure, case No. In 1 weekends is, Case No. In February - gold has been set.
[0123]
　Upset detector 12, by setting the case with monthly or seasonal 603 may reflect the effect of a variation in tidal conditions in the system in one year statistical processing region motion information. Change of tide situation in this case, cooling demand in summer, are due to the life behavior of the month or every season to an example of the heating demand in winter. In the example of this figure, case No. March in 1 - May is, case No. 2 In July - September has been set.
[0124]
　Upset detector 12, by setting the case with the event type 604, may reflect system fault, the effect of a variation in tidal conditions in the system due to events such as system switch to the statistical processing of the area motion information. In the example of this figure, case No. 1 In no event, Case No. 1LG of 2, node AA has been set.
[0125]
　Although not illustrated in FIG. 15, upset detector 12, in addition to the setting method of the above example, power flow amount of pre-configured node or transmission lines to be monitored, the power generation amount of the generator, sunny, cloudy - weather information, including a section or temperature information of rain, set the case with the classification item, it may be carried out statistical processing of the area motion information.
[0126]
　Here, the description upset detection unit 12 about the past motion information calculating process of calculating a past motion information table TB7 in Figure 15. Figure 16 is an example of a flowchart showing the past motion information calculating process. Note the process flow of FIG. 16 will be performed positioned downstream of the upset influence estimation processing flow of FIG.
[0127]
　In Figure 16, motion detector 12 first stores the measurement information of a plurality time points in the process step S161.
[0128]
　Then, shaking detection unit 12, using the measurement information stored in the processing step S162, it calculates the area motion information and perturbations impact for multiple time points. At this time, agitation detection unit 12 in accordance with step of the process steps S11 ~ S12 of FIG. 5 described above, calculates the area motion information and upset influence.
[0129]
　Then, shaking detection unit 12 sets the classification of the case in the process step S163. At this time, agitation detection unit 12 as the classification items, the time shown in the past motion information, date, day, month, season, using at least one of information of an event type. Upset detector 12 for each of the categories to be used to determine start and end of the period, the classification name of the event type, or the minimum and maximum values ​​of the numerical information, sets a classification category.
[0130]
　Next, in the process step S164 of upset detector 12, based on the classification of the cases set in the process step S161, it calculates the area motion information for each case, the agitation impact. At this time, agitation detection unit 12, region motion information in the processing step S162, classifies an analysis time calculated sway impact every case, information of the area motion information and perturbations impact for several times, which is classified into the same case a representative value calculated by performing statistical processing, the region motion information of each case, and upset impact. Upset detector 12, as a representative value, region motion information, the average value for each upset impact, the intermediate value is used.
[0131]
　The above is the past motion information calculation process. In still upset influence calculation processing, motion detection unit 12, by referring to the past motion information calculated by the past motion information calculation processing to extract the closest examples to the situation of the analysis time, agitation impact of the extracted Case it may be applied. Upset detector 12, by using a statistical processing result of a past region motion information, exactly the fluctuation component that eliminates upset by momentary effect of noise included in the system oscillations in the analysis time, preferentially monitored it can be selected.
[0132]
　Note After the processing step S122 in FIG. 13, upset detector 12 in the process step S123, in accordance with (16), upset influence J region motion information j in measurement point i ij is calculated. (16) is a definition equation of the upset influence J in consideration of multiple types. Here, beta equation (16) k is a coefficient set for the measurement information type k. Beta K may be as the reciprocal of the nominal value of each measurement type k. Upset detector 12, by using the nominal value of each measurement type, the upset impact for measurement information of a plurality of types, can be calculated by normalizing a ratio to rated. Further, in order to extract the perturbation effect of a particular measurement type, the coefficient β corresponding to the measurement type k k may be set larger than the value of the coefficient for the other types.
[0133]
[Number 16]

The scope of the claims
[Claim 1]
　Based on the measurement information at a plurality measurement point of the power within the system, a system stabilizing control apparatus for controlling a controlled device installed in the power within the system,
　stores the information of the system equipment is interconnection to the power system equipment and information database, a control priority determining section for determining the control priority of the control device based on the system equipment information, the control target determination for determining a control target based on the measurement information and the control priority of the plurality measurement point parts and, said plurality destabilization time to calculate a destabilization time of the power system from the measurement information in the measurement point calculator, the control time of said control device based on the destabilization time and the system equipment information system stabilizing control device characterized by and a control time determination section for determining.
[Claim 2]
　A system stabilizing control device according to claim 1,
　wherein the destabilization time calculation unit, agitation frequency of sway of the frequency components of the power system, by using the damping coefficient, and at least one of the amplitude non system stabilizing control device and calculates the stabilization time.
[Claim 3]
　A system stabilizing control device according to claim 2,
　wherein the information of the system equipment, the control means of the system equipment including the control device, the controllable amount, and characterized in that it comprises at least one of operation time system stabilizing control unit for.
[Claim 4]
　A system stabilizing control device according to claim 3,
　wherein the information of the system equipment, and latitude and longitude of each of the plurality of measurement points, the topology between the plurality of measurement points, the plurality of measurement and impedance between points, the plurality of the characteristics of the power devices that interconnection in the vicinity of each of the measurement point, include at least one, the system stabilizing control device, characterized in that.
[Claim 5]
　A system stabilizing control device according to claim 4,
　wherein the facility information database stores the past motion information, the shaking detection unit, the past motion information corresponding to the situation of the measurement point group select, based on the duration of the frequency components in the past motion information of the selected, the system stabilizing control apparatus characterized by and calculates the upset impact.
[Claim 6]
　A system stabilizing control device according to claim 1,
　wherein the control target determination unit, between the electrical distance or said two points between said control equipment interconnection point and the measurement point where instability upset is observed system stabilizing control apparatus characterized by determining a control target using an impedance.
[Claim 7]
　Based on the measurement information at a plurality measurement point of the power within the system, a system stabilizing control apparatus for controlling a controlled device installed in the power within the system,
　the system equipment including the control device to be interconnection to the power system a facility information database storing facility information, and a control priority determining section for determining the control priority of the control device based on the facility information, and upset the measurement point of the power system oscillations which is detected from the measurement information of the controlled device distance, a control target determination section for determining the control apparatus to be controlled according to the priority, and the destabilization time calculation unit for calculating a destabilization time of the power system oscillations which is detected from the measurement information, wherein line depreciation, characterized in that it executes control of the destabilization time and said a control time determination section for determining the control time based on the facility information, control time determined for determined control equipment Of the control device.
[8.]
　Obtained system stabilizing control device measurement information in multiple measurement point of the power in the system through the first communication means, via the second communication means to the control device of the system features which constitute the electric power system a power system control system for providing a control command from the system stabilizing control device,
　the system stabilizing control device includes a facility information database for storing information of the system equipment is interconnection to the power system, the system equipment wherein the control priority determining section for determining the control priority of the control device, a control target determination section for determining a control target based on the measurement information and control priority in the plurality measurement point, the measurement information in multiple measurement point based on the information a destabilization time calculation unit for calculating a destabilization time of the electric power system from a control time determination section for determining the control time based on the destabilization time and the facility information, the determined Power system control system, characterized in that it comprises an output unit for the control command give via the second communication means comprising a control device and determined control time that.
[Claim 9]
　A power system control system of claim 8,
　upset comparing the system stabilizing controller for comparing the assumed perturbation calculation unit that calculates an assumed upset based on the facility information, the measurement information and the assumed perturbation and a section, power system control system, wherein the control object determination section for determining a control target based on the upset comparison result and the control priority.
[Claim 10]
　A power system control system of claim 8,
　wherein the system is connected to the output of the stabilizer control device, obtaining information including the control time determined with the control device to which the determined total system stabilizing control device the provided,
　該全system stabilizing control device, the received control device and the additional control decision section based on the control time determining additional control devices and additional control time, the additional control devices and additional control time the system features power system control system, characterized in that it comprises an output unit for transmitting to the control equipment.