The present invention relates to a distributed-generator power-generation-amount estimation apparatus and method.
In recent years, the interconnection amount of distributed generators DG including photovoltaic power generation to an electric power system is increasing for the purpose of protecting the global environment. On the other hand, a distributed generator DG has a risk of sudden change in power generation due to weather such as clouds and wind. Power companies have a need to know the output of the distributed generator DG frequently, such as at intervals of several seconds, and precisely in order to stably supply electric power.
As a conventional technique, Non-Patent Document 1 discloses a technique for estimating output of a distributed generator DG. The document discloses a method for estimating output of a distributed generator DG by using a power-factor method.
In the estimation of the output of a distributed generator DG by the power factor method, by assuming a coordinate plane in which the abscissa is active power P and the ordinate is reactive power Q, putting plots on the PQ coordinate plane according to the active power P and reactive power Q measured in units of banks of a substation or the like while the distributed generator DG is stopped, and performing, for example, linear approximation processing based on these plots, a load characteristic of the electric power system is specified in advance. This load characteristic is a characteristic when there is no distributed generator in the electric power system.
Next, plots are put on the above active power/reactive power coordinate plane according to real time active power P and reactive power Q (measured, for example, at 10-second intervals) measured in units of banks of the substation or the like. Supposing that the distributed generator DG operates at a power factor of 1, the distributed generator DG outputs only active power. Therefore, when it is assuming that the real-time active power measurement value is P1 and the real-time reactive power measurement value is Q1, the output of the distributed generator DG is calculated by Formula (1).
[Mathematical Formula 1]
output of distributed generator DG=active power P0 obtained from load characteristic when reactive power is Q1−measured real time active power P1 (1)
The advantage of the power-factor method is that it is possible to calculate the output of the distributed generator DG at that point each time active power and reactive power is measured in units of banks of a substation or the like, and then it is possible to estimate the output of the distributed generator DG frequently.
Measurement of active power and reactive power in units of banks of a substation means that active power and reactive power in an electric-power supplying target area connected via a feeder (distribution line) to the bank. For this reason, the estimation of the power generation amount means to estimate the power generation amount of the distributed generator in units of areas.
NPL 1: Matsuda et al., “Bunsangata Dengen no Untenjoukyou Suiteisyuhou no Teian” IEEJ Transactions on Power and Energy, Vol. 130 No. 6 Pages 593-602 (2010)
In the method disclosed in Non-Patent Document 1, since the power generation amount of the distributed generator DG is estimated based on the assumption that the power factor of the distributed generator DG is 1, there is a problem that an error is included in the power-generation-amount estimation value of the distributed generator DG in the case where the distributed generator DG is operated at a power factor other than 1.
On the other hand, since it is generally assumed that a power factor of a distributed generator DG is other than 1, it is desirable to calculate the power factor of the distributed generator DG and use it for estimating the output of the distributed generator DG.
In view of the above, the present invention is to provide a distributed-generator power-generation-amount estimation apparatus and method that can estimate output more precisely and easily by reflecting a power factor.
In order to solve the above problem, a configuration described in, for example, claims is adopted. The present invention includes a plurality of means for solving the above problem, and is, for example, a distributed-generator power-generation-amount estimation apparatus that estimates, from active power and reactive power measured at regular intervals by measurement means installed on an electric power system, a power generation amount of a distributed generator interconnected the an electric power system, the distributed-generator power-generation-amount estimation apparatus including power-factor estimation means for estimating, as a ratio of a time change amount of the active power and the reactive power, a power factor of the distributed generator interconnected to a demand side of the electric power system from a measurement point by the measurement means, distributed-generator power-generation-amount estimation means for having, from active power and reactive power measured at a time of no power when the distributed generator interconnected to the electric power system is not generating power, a first characteristic line representing a load characteristic on a coordinate plane based on active power and reactive power, determining, from the measured active power, the measured reactive power, and the power factor of the distributed generator, a second characteristic line representing a distributed-generator characteristic on the coordinate plane, and estimating the difference between the active power and the reactive power at an intersection point of the first characteristic line and the second characteristic line, and the measured active power and the measured reactive power as the power generation amount of the distributed generator, and output means for outputting power-generation-amount estimation information on the distributed generator estimated by the distributed-generator power-generation-amount estimation means.
Furthermore, a distributed-generator power-generation-amount estimation method for estimating, from active power and reactive power measured at regular intervals by measurement means installed on an electric power system, a power generation amount of a distributed generator interconnected to the electric power system, the distributed-generator power-generation-amount estimation method includes estimating, as a ratio of a time change amount of the active power and the reactive power, a power factor of the distributed generator interconnected to a demand side of the electric power system from a measurement point by the measurement means, and obtaining, from active power and reactive power measured at a time of no power when the distributed generator interconnected to the electric power system is not generating power, a first characteristic line representing a load characteristic on a coordinate plane based on active power and reactive power, obtaining, from the measured active power, the measured reactive power, and the power factor of the distributed generator, a second characteristic line representing a distributed-generator characteristic on the coordinate plane, and estimating the difference between the active power and the reactive power at an intersection point of the first characteristic line and the second characteristic line, and the measured active power and the measured reactive power as the power generation amount of the distributed generator.
According to the present invention, it is possible to more accurately and easily estimate a power generation amount of a distributed generator DG as compared with a conventional technique.
Thus, according to an embodiment of the present invention, a business entity operating an electric power system can prepare the setting states of the facilities in the electric power system in appropriate states so as to handle any sudden change in the power generation amount of the distributed generator DG at any time. Accordingly, it is possible to increase the interconnection amount of a distributed generator DG to an electric power system, and to improve the electric power quality more than before when the power generation amount of the distributed generator DG suddenly changes.
An Embodiment of the present invention will be described below with reference to the drawings.
In the input function I, a setting condition used particularly in the calculation function C of the area-unit distributed-generator power-generation-amount estimation apparatus is defined. The input function I includes an estimation-condition input function I1, a normal-range input function I2, and a date/time designation function I3. The estimation-condition input function I1 is used for inputting a setting condition for estimating a power generation amount of a distributed generator DG. The normal-range input function I2 is used for determining whether a power-factor estimation value of the distributed generator DG is an abnormal value. The date/time designation function I3 is used for setting the date and time when data is extracted from a database to verify estimation values of the power factor or the power generation of the distributed generator DG in the past. The input function I will be described later with reference to
The output function O includes a power-generation-amount output function O1, an estimation-process output function O2, a power-factor output function O3, and a history output function O4. The power-generation-amount output function O1 is used for outputting an estimation value of the power generation amount of the distributed generator DG. The estimation-process output function O2 is used for outputting a calculation process of estimating the power generation amount of the distributed generator DG as an image. The power-factor output function O3 is used for outputting a power-factor estimation value of the distributed generator DG. The history output function O4 is used for outputting actual data of the estimation values of the power factor and the power generation amount of the distributed generator DG. The output function O will be described later with reference to
The storage function M is configured, for example, as a database DB, and includes a tidal-current-measurement-value database DB1, a power-factor-control-information database DB2, a load-characteristic database DB3, and an estimation-value database DB4. The tidal-current-measurement-value database DB1 stores a measurement value D1 of a tidal current. The power-factor-control-information database DB2 contains power-factor control information D2. The load-characteristic database DB3 contains load-characteristic data D3. The estimation-value database DB4 stores an estimation value D4 of a calculation result. The power-factor-control-information database DB2 stores a power-factor instruction value when a control function (not shown) of the electric power system controls a power factor of the distributed generator DG in real time.
The calculation function C includes a power-factor estimation function C1, an abnormal-value determination function C2, and a power-generation-amount estimation function C3.
The power-factor estimation function C1 estimates, from the estimation condition received from the estimation-condition input function I1 and a tidal-current measurement value D1 (active power and reactive power) measured in the electric power system and stored in the tidal-current-measurement-value database DB1, a power factor as a whole including a plurality of distributed generators DG on the downstream of the measurement point of the tidal-current measurement value, and outputs the power factor to the abnormal-value determination function C2. When the power generation amount of distributed generators in units of areas is estimated, the tidal-current measurement value D1 (active power and reactive power) is set as the bank installation position of the substation in the electric power system. The power-factor estimation function C1 will be described later with reference to
The abnormal-value determination function C2 determines, based on the information on the power-factor normal range input by the normal-range input function I2, whether the power-factor estimation value output by the power-factor estimation function C1 is within the normal range. Incidentally, the power-factor-control-information database DB2 stores a power-factor instruction value when the control function of the electric power system controls a power factor of the distributed generator DG in real time. When the power-factor control value at the current time is stored in the power-factor-control-information database DB2, the abnormal-value determination function C2 replaces the power-factor estimation value with the power-factor control value, and determines whether the power-factor control value is within the normal range. The power-factor estimation value after the determination is output to the power-generation-amount estimation function C3, the power-factor output function O3, and the estimation-value database DB4. The abnormal-value determination function C2 will be described later with reference to
The power-generation-amount estimation function C3 estimates, from power-factor information output by the abnormal-value determination function C2, load-characteristic information D3 of the load-characteristic database DB3, and tidal-current-measurement-value information D1 of the tidal-current-measurement-value database DB1, the power generation amount (active power and reactive power) of all the distributed generators DG on the downstream of the tidal-current measurement point. The information on the estimated power generation amount is output to the power-generation-amount output function O1 and the estimation-value database DB4. The image of the calculation process is output to an estimation-process output function O3 so that the operator can visually understand the calculation process in the power-generation-amount estimation function C3. The power-generation-amount estimation function C3 will be described later with reference to
Among the information of the power-factor estimation value and the power-generation estimation value of the distributed generator DG recorded in the estimation-value database DB4, the data matching the period output by the date/time designation function I3 is output to the history output function O4.
Hereinafter, the processing of each function will be described in detail.
The estimation-condition input function I1 sets a power-generation-amount estimation cycle and a tidal-current-measurement-value interval for estimating a power factor. The power-generation-amount estimation cycle is a function for setting an interval for executing the estimation processing of the power generation amount of the distributed generator DG interconnecting to the electric power system. In this example, it is assumed that the processing is performed at a 10-second cycle. The tidal-current-measurement-value interval for estimating a power factor is a function for setting a time interval of measurement data, which is tidal-current data (active power and reactive power) measured at another time, from which the power factor of the distributed generator DG is estimated. In this example, it is assumed that the processing is performed at a 10-second cycle.
The normal-range input function I2 designates, as the power-factor estimation value of the distributed generator DG, the normal range in advance. In this example, the normal range is designated as the upper limit of 1.0 and the lower limit of 0.5.
The date/time designation function I3 is a function for setting, when the operator refers to the actual data of the power-factor estimation value of the distributed generator DG and of the power-generation-amount estimation value of the distributed generator DG, the reference period. In the illustrated example, 30 minutes from 12:00 on Aug. 1, 2015 to 12:30 on the same day are designated as the reference period. In
The power-generation-amount output function O1 is a function for outputting power-generation-amount prediction values (active power estimation value Kw and reactive power estimation value KVar) of the distributed generator DG output by the power-generation-amount estimation function C3 as a numerical value or graph.
The estimation-process output function O2 is a function for visually displaying the process in which the power-generation-amount estimation function C3 estimates the power generation amount of the distributed generator DG. In the illustrated example, the process is displayed on the coordinate plane based on the active power and the reactive power.
The power-factor output function O3 is a function for displaying the power-factor estimation value output by the abnormal-value determination function C2 as a numerical value or graph.
The history output function O4 is a function for displaying data of a period designated by the date/time designation function I3 from the past power-factor estimation value of the distributed generator DG and power-generation-amount estimation value of the distributed generator DG that are stored in the estimation-value database DB4. In the illustrated example, the active power, the reactive power, and the power factor for 30 minutes from 12:00 on Aug. 1, 2015 to 12:30 on the same day designated by date/time designation function I3 are displayed as numerical values.
Here, the fluctuation manner of the active power P and the reactive power Q measured in the electric power system is based on that the fluctuation in the load of the consumer interconnected to the system overlaps with the fluctuation in the output of the distributed generator DG. In the present invention, the fluctuation in the active power P and the reactive power Q is caused by that the fluctuation in the load of the customer overlaps with the fluctuation in the output of the distributed generator DG, but that the output of the distributed generator DG generally fluctuates more abruptly than the fluctuation in the load of the consumer is focused on. In other words, if it is limited to a very short time window (several seconds to several tens of seconds), it can be assumed that the change in the load of the consumer is almost 0, and, thus, the change in the active power and the reactive power in this time window is caused by the output of the distributed generator DG. Therefore, by showing the measurement data of 30 seconds from 12:00:00 to 12:00:30 in
[Mathematical Formula 2]
Power-factor estimation value=cos(tan−1(ΔQ/ΔP)) (2)
Where ΔP is the difference in the active power (kW) measured at different times and ΔQ is also the difference in reactive power (kvar). The time intervals for calculating ΔP of the active power and ΔQ of the reactive power are the same, and the timings are also the same. In short, this processing is to determine the power factor from the ratio of the time change amount of the active power and the reactive power.
In the first processing step S11, a tidal-current measurement value at the current time is acquired. In this processing, the measurement values of the active power and the reactive power at the current time are acquired by referring to the tidal-current-measurement-value database DB1.
In the next processing step S12 in
In the processing step S13, a past tidal-current measurement value is acquired. Here, the active power and the reactive power at the time going back by the time length acquired in the processing step S12 from the current time are acquired from the tidal-current-measurement-value database DB1.
In the processing step S14, the active-power change amount ΔP is calculated. ΔP is calculated by Formula (3).
[Mathematical Formula 3]
ΔP=Pt−Pt−Δt (3)
Where, Pt is the active power at the current time, and Pt−Δt is the active power measured at a certain past time.
In the processing step S15, the reactive-power change amount ΔQ is calculated. ΔQ is calculated by Formula (4).
[Mathematical Formula 4]
ΔQ=Qt−Qt−Δt (4)
Where, Qt is the reactive power at the current time, and Qt−Δt is the reactive power measured at a certain past time.
The processing step S16 is a power-factor estimation function. This function calculates the power-factor estimation value using ΔP and ΔQ by Formula (2).
The processing step S17 is output processing of the power-factor estimation value. In this processing, the power-factor estimation value at the current time determined in the processing steps S14, S15, and S16 is output to the abnormal-value determination function C2.
In this embodiment, ΔP is compared with the reference value in processing step S15, but it is also possible to compare ΔQ with the reference value or to compare both ΔP and ΔQ with the reference value.
In the first processing step S21 in
In the processing step S22, a power-factor estimation value is acquired. The power-factor estimation value is acquired from the power-factor estimation function C1.
The processing step S23 is a branching function for determining whether there is power-factor control information at the current time. The power-factor control information referred in this description is information on the power-factor instruction amount in a case where the control system of the electric power system performs control to intentionally lower the power factor of the distributed generator DG in order to stably supply electric power. The power-factor-control-information database DB2 is referred to in order to recognize whether there is the power-factor control information, and the magnitude thereof. Here,
In the processing step S23, the power-factor-control-information database DB2 is referred to in order to confirm whether there is power-factor control information at the current time. When there is power-factor control information at the current time, the processing proceeds to the processing step S24, and the power-factor estimation value received from the power-factor estimation function C1 is replaced with the value recorded in the power-factor-control-information database DB2. On the other hand, when there is no power factor control information at the current time in the power-factor-control-information database DB2, the processing proceeds from the processing step S23 to the processing step S25.
In the processing step S25, it is determined whether the power-factor estimation value is within the normal range. When Formula (5) is satisfied at this function, the power-factor estimation value is determined to be a normal value, or when Formula (5) is not satisfied, the power-factor estimation value is determined to be an abnormal value.
[Mathematical Formula 5]
Lower limit value of power-factor normal range≤power-factor estimation value≤upper limit value of power-factor normal range (5)
Where, the lower limit value and the upper limit value of the power-factor normal range are the setting information of the normal-range input function I2. When it is determined that the power-factor estimation value is in the normal range in the processing step S25, the processing proceeds to the processing step S27. On the other hand, when it is determined that the power-factor estimation value is an abnormal value in the processing step S25, “99” is set as the power-factor estimation value in order to indicate that the power factor is an abnormal value in S26.
In the processing step S27, the processing for outputting the power-factor estimation value is performed. The output destinations of the power-factor estimation value are the power-generation-amount estimation function C3, the power-factor output function O3, and the estimation-value database DB4.
As a technique for obtaining the load-characteristic line L1, the technique described in Non-Patent Document 1 can be applied. By putting plots on the PQ coordinate plane according to the active power P and the reactive power Q measured in units of banks of the substation or the like while the distributed generator DG is stopped and by performing, for example, linear approximation processing based on these plots, the load characteristic of the electric power system can be previously specified. This obtaining method is effective when it can be confirmed that the distributed generator DG is stopped.
In contrast, as an actual operation of the electric power system, if the distributed generator DG only includes photovoltaic power generation, it can be determined that the operation is stopped by obtaining nighttime data, but if the distributed generator DG includes wind power generation, it is necessary to separately confirm nighttime and stoppage of the wind power generation.
As a method for obtaining a line L1 of another load characteristic in consideration of this point, by putting plots on the coordinate plane according to the active power P and the reactive power Q measured in units of banks of the substation or the like in the past without the condition of operation stop, and by performing, for example, linear approximation processing so that the magnitude of the active power P or the reactive power Q is to be maximum in the fourth quadrant (active power P is positive, reactive power Q is negative) in the coordinate plane or in the first quadrant (the active power P and the reactive power Q are positive), the load characteristic of the electric power system can be previously specified.
Next, real-time measurement values of the active power and reactive power are put on the coordinate (indicated as a point N in
Considering the relationship between the point N and the point R, the point R corresponds to the active power and reactive power measured at the measurement point when all the distributed generators DG are not generating power at the current time. On the other hand, the point N is a measurement value, including the power generation effect of the distributed generator DG, of the active power and the reactive power. That is, the divergence width between the point N and the point R is the power generation amount of the distributed generator DG. Thus, the active power component of the power generation amount of the distributed generator DG can be calculated as a value PDG obtained by subtracting the active power PL at the point C from the active power PH at the point R. On the other hand, the reactive power component of the power generation amount of the distributed generator DG can be calculated as the value QDG obtained by subtracting the reactive power QL at the point R from the reactive power QH at the point N.
In the first processing step S31 in
In the processing step S31, according to the acquired load-characteristic information, the load-characteristic line L1 is drawn on the active power-reactive power coordinate plane as described in the description of
In the processing step S32, a measurement value at the current time is acquired. The tidal-current-measurement-value database DB1 is referred to in order to acquire the tidal-current measurement value (active power and reactive power) at the current time.
In the processing step S33, a power-factor estimation value at the current time is acquired. The power-factor estimation value at the current time is acquired from the abnormal-value determination function C2.
The processing step S34 is a branching function for determining whether the power-factor estimation value at the current time is a normal value. That is, when the power-factor estimation value is not “99”, the value is determined to be a normal value, and the processing proceeds to the processing step S35. On the other hand, when the power-factor estimation value is “99”, the value is determined to be an abnormal value, and the processing proceeds to the processing step S38.
In the processing step S35, a power-factor characteristic is generated. As described in the description of
In the processing step S36, on the active power-reactive power coordinate plane, the coordinates of the intersection point R of the load-characteristic line L1 generated in the processing step S31 and the power-factor characteristic line L2 passing through the point N of the tidal-current measurement value generated in the processing step S35 is calculated.
In the processing step S37, the total power generation amount of the distributed generator DG at the current time is estimated. The magnitude of the active power of the total power generation amount of the distributed generator DG is calculated by Formula (6), and the magnitude of the reactive power is calculated by Formula (7).
[Mathematical Formula 6]
Active power generation amount of distributed generator DG=PH−PL (6)
[Mathematical Formula 7]
Reactive power generation amount of distributed generator DG=QH−QL (7)
Where, PH is the active power component (kW) at the coordinates of the intersection point R of the load-characteristic line L1 and the power-factor characteristic line L2, PL is the active power measurement value (kW) at the current time, QH is the reactive power measurement value (kval) at the current time, and QL is the reactive power component (kval) at the coordinates of the intersection point R of the load-characteristic line L1 and the power-factor characteristic line L2.
In the processing step S38, since the power-factor estimation value is an abnormal value, the character “*” meaning “unknown” is set as the power generation amount of the distributed generator DG.
In the processing step S39, the processing for outputting the power-generation-amount estimation value is performed. The output destinations of the power-generation-amount estimation value are the power-generation-amount output function O1 and the estimation-value database DB4.
In the processing step S13, an image of the process of calculating the power-generation-amount estimation value is output. In this processing, an image diagram visualizing the processing from the processing steps S31 to S38 (that is, the image of
According to the present invention, by using the method shown in
Note that, the present invention is not limited to the above embodiment, and includes various modifications. For example, the above embodiment has been described in detail in order for the present invention to be easily understood, and is not necessarily limited to those having all the described configurations.
Number | Date | Country | Kind |
---|---|---|---|
2015-178148 | Sep 2015 | JP | national |
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/JP2016/067945 | 6/16/2016 | WO | 00 |