The present invention relates to a prediction device, a prediction method, and a program.
Priority is claimed on Japanese Patent Application No. 2018-156568, filed Aug. 23, 2018, the content of which is incorporated herein by reference.
In the related art, a plant or a mechanical device may use a prediction model to guide monitoring or the like. For example, process data is collected and physical models of a device and the like provided in the plant or a statistical technique is used to construct a prediction model. Then, a value which is a norm for a process amount is obtained from the constructed prediction model, and the value is used to perform monitoring, control, and determination of abnormality. Patent Literature 1 discloses a technique where a plurality of data sets that are a part selected from various process data are acquired, each of the data sets is used to construct a prediction model, and an integrated value of prediction values calculated by a plurality of the constructed prediction models is used to monitor a plant or the like.
[PTL 1] Japanese Unexamined Patent Application Publication No. 2015-127914
In general, there is an error in the accuracy of prediction by the prediction model, and there is a possibility that the prediction value is shifted to an unsafe side by the amount of the error. For this reason, when the prediction value based on the prediction model is used as it is to monitor or control the plant, there is a possibility of leading to an undesired result.
In Patent Literature 1, among the plurality of prediction models, a large weight is applied to the prediction value calculated from the prediction model having a small error and a small weight is applied to the prediction value calculated from the prediction model having a large error to calculate a weighted mean of the prediction values to integrate the prediction values, so that the influence of the errors of the prediction models is reduced. However, for example, when a plurality of appropriate prediction models cannot be constructed or the like, the method described in Patent Literature 1 cannot be used.
The present invention provides a prediction device, a prediction method, and a program capable of solving the above-described problem.
According to one aspect of the present invention, there is provided a prediction device including: a data collection unit that collects process data of a device; a prediction model construction unit that constructs a prediction model having a predetermined input variable of first process data as an input value and having a predetermined output variable of the process data as an output value, and an error calculation model which calculates a prediction error of the prediction model, based on the first process data collected by the data collection unit; and a prediction unit that outputs a corrected prediction value which is obtained by correcting a prediction value of the output variable with the prediction error calculated based on the error calculation model, the prediction value being calculated based on the input variable of second process data collected by the data collection unit and on the prediction model.
According to one aspect of the present invention, the prediction unit adds or subtracts the prediction error to or from the prediction value to correct the prediction value such that the corrected prediction value is not safer or is less efficient than the prediction value before correction.
According to one aspect of the present invention, the prediction device further includes a state monitoring unit that compares the process data with a predetermined threshold value to determine whether or not the process data is abnormal; and an operation-amount determination unit that calculates an operation amount which improves the corrected prediction value when the state monitoring unit determines that the process data is abnormal.
According to one aspect of the present invention, the prediction device further includes a first output unit that outputs the operation amount, which is calculated by the operation-amount determination unit, to a control device of the device.
According to one aspect of the present invention, the prediction device further includes a second output unit that displays the corrected prediction value and a graph, which visualizes the prediction model, in a superimposed manner.
According to one aspect of the present invention, there is provided a prediction method including: a step of collecting process data of a device; a step of constructing a prediction model having a predetermined input variable of first process data as an input value and having a predetermined output variable of the process data as an output value, and an error calculation model which calculates a prediction error of the prediction model, based on the first process data collected in the step of collecting the process data; a step of collecting second process data of an evaluation target; and a step of outputting a corrected prediction value that is obtained by correcting a prediction value of the output variable with the prediction error calculated based on the error calculation model, the prediction value being calculated based on the input variable of the collected second process data and on the prediction model.
According to one aspect of the present invention, there is provided a program that causes a computer to function as means for collecting process data of a device; means for constructing a prediction model having a predetermined input variable of first process data as an input value and having a predetermined output variable of the process data as an output value, and an error calculation model which calculates a prediction error of the prediction model, based on the first process data collected in a step of collecting the process data; means for collecting second process data of an evaluation target; and means for outputting a corrected prediction value that is obtained by correcting a prediction value of the output variable with the prediction error calculated based on the error calculation model, the prediction value being calculated based on the input variable of the collected second process data and on the prediction model.
According to the prediction device, the prediction method, and the program, the prediction value based on the influence of the prediction error of the prediction model can be output.
Hereinafter, a prediction device according to a first embodiment of the present invention will be described with reference to
The plant illustrated in
The prediction device 30 acquires current various process data from the gas turbine 10 to predict an operating state of the gas turbine 10 based on the acquired process data and a prediction model. For example, a value predicted by the prediction device 30 may be the value of process data, which represents an operating state of the gas turbine 10 for the future ahead of a predetermined time, or may be an estimate value that estimates a value that is not directly measurable. Here, the process data is, for example, measurement data such as temperature and pressure measured by sensors provided at places in the gas turbine 10 and the generator 15. The measurement data includes physical property data of the fuel gas and atmosphere air that are taken into the gas turbine 10 to be used for actual operation, and measurement data of operating environments such as atmospheric temperature and humidity. The measurement data includes identification information of the sensors, measurement values, measurement times, and the like. The process data includes control values (opening degree command values of the fuel flow regulation valves 16A to 16C) that are generated by the device 20 to control the gas turbine 10. The process data includes values converted from the acquired process data or values calculated from a plurality of process data. The prediction device 30 of the present embodiment can output a prediction value corrected to be on a safer side in consideration of a prediction error of the prediction model, whereas a general prediction device outputs a prediction value based on the prediction model. Next, the prediction device 30 will be described.
As illustrated in
The data collection unit 31 collects process data from the plant or a mechanical device that is a monitoring target.
The data storage 32 stores the process data, which is collected by the data collection unit 31, in the storage unit 37.
The data extraction unit 33 extracts data, which is required to construct a prediction model, from the process data collected by the data collection unit 31. For example, the data extraction unit 33 extracts data of the types required to construct the prediction model or extracts values in the required range (removes outliers and the like). In the case of a prediction model that predicts the combustion vibration of the combustor 12, the data of the types required to construct the prediction model is, for example, vibration data obtained by measuring the vibration of combustion air inside the combustor 12 (or data obtained by analyzing the frequency of vibration data by fast Fourier analysis), the opening degree command values of the fuel flow regulation valves 16A to 16C, the inlet temperature of the turbine 13, the angle of the IGV 17, and the like.
The prediction model construction unit 34 constructs a prediction model that predicts an operating state of the plant or the mechanical device by a statistical technique such as multiple regression analysis or Gaussian process regression, machine learning such as random forest, or deep learning such as a neural network. The prediction model construction unit 34 constructs an error calculation model that calculates an error of the constructed prediction model (variation or uncertainty of prediction). For example, the prediction model construction unit 34 constructs a prediction model and an error calculation model that have the values of predetermined input variables of the process data extracted by the data extraction unit 33, as input values, and have the value of a predetermined output variable as an output value to learn a relationship between both the input and output variables. The input variables are, for example, the opening degree command values of the fuel flow regulation valves 16A to 16C, the inlet temperature of the turbine 13, the degree of the IGV 17, atmospheric temperature, atmospheric humidity, an output of the gas turbine 10, vehicle compartment pressure, and the like. The output variable is, for example, the level of combustion vibration of the combustion air inside the combustor 12, emissions such as NOx and CO, a performance index such as output efficiency, or the like. For example, the prediction model construction unit 34 constructs a prediction model (function or the like) for combustion vibration which has the values of predetermined input variables (opening degree command values of the fuel flow regulation valves 16A to 16C, the inlet temperature of the turbine 13, and the degree of the IGV 17) of the process data extracted by the data extraction unit 33, as input values, and has the value of a predetermined output variable (vibration data of combustion vibration) as an output value to define a relationship between both the input and output variables. The error calculation model is a calculation formula that calculates a statistic such as a difference between the process data given as teacher data and the root mean square of the prediction value. In the case of the prediction model using regression analysis, a confidence interval of the prediction value may be used, and in the case of the prediction model using Gaussian process regression, an error directly obtained by a Gaussian process regression technique may be used. Examples of the prediction model and the error calculation model will be described later with reference to
The prediction unit 35 predicts an output value based on a prediction error for the predetermined output variables, based on the input variables of the process data collected by the data collection unit 31, the prediction model, and the error calculation model. At this time, the prediction unit 35 corrects the value of the output variable, which is predicted by the prediction model, with the value of the prediction error for the value of the output variable to generate a final prediction value. More specifically, the prediction unit 35 adds or subtracts the prediction error to or from the prediction value such that the corrected value is not safer than the value before correction or the corrected value is less efficient than the value before correction. The prediction unit 35 outputs the prediction value after addition or after subtraction (after correction) as the final prediction value. In such a manner, the prediction unit 35 uses the prediction error to obtain the prediction value on a safe side in terms of equipment protection or contract. For example, when the output variable is combustion vibration or the amount of emission of NOx or CO, the prediction unit 35 adds a prediction error to a prediction value to correct the prediction value in an increasing direction. When the output variable is a variable related to efficiency, the prediction unit 35 subtracts a prediction error from a prediction value to correct the prediction value in a decreasing direction to thus calculate a final prediction value.
The output unit 36 outputs a prediction result.
The storage unit 37 stores the process data, the prediction model, the error calculation model, and the like.
Here, the prediction model and the error calculation model will be described.
In the case of multiple regression analysis, a prediction value y is expressed by an equation using a plurality of explanatory variables x1, x2, etc. When the prediction model using single regression is considered for convenience of description, the prediction value y can be obtained by the following equation using an explanatory variable x.
y=α+βx (1)
At this time, a variation (error) σe{circumflex over ( )} of the prediction value y is estimated by the following equation (2).
Here, the hat ({circumflex over ( )}) means an estimate value, n means the data count, and i means the data number.
Therefore, the prediction value (mean value) and the distribution of variations (errors) thereof are as illustrated in
The same idea applies when a multivariate regression spline, a feedforward neural network, or the like is used as the prediction model.
The horizontal axis of
When the prediction model construction unit 34 constructs a prediction model by random forest regression, the prediction model construction unit 34 calculates the prediction model and an error calculation model which are visualized by a prediction value Y (line 6b) having a step shape and lines 6a and 6c, which illustrate a variation (error) of 2σ centered on Y, for an explanatory variable X1 as illustrated in
When the prediction model construction unit 34 constructs a prediction model by Gaussian process regression, the prediction model construction unit 34 can calculate lines 7a and 7c illustrating variations with respect to a line 7b illustrating the prediction model. In the case of Gaussian process regression, an error which differs according to the magnitude of the explanatory variable X1 can be calculated as illustrated.
In the case of Gaussian process regression, a distribution f(x) of a response surface can be obtained from data D (aggregation of sets of the explanatory variable x and an output y) as expressed by the following equation (3).
p(f(x)|D)=N(kt(K+σ2IN)−1y,
K
0
−k
t(K+σ2IN)−1k) (3)
Here, when σ is a variance of observation noise, σp is a variance of the prior distribution of a prediction target, and θ is a scaling parameter, p(y|x, σ2), K0, k, and K(x, x′) are as follows.
p(y|x,σ2)=N(y|f(x),σ2) (4)
K
0
=K(x,x),k=(K(x,x1), . . . ,K(x,xN))t (5)
At this time, for example, the prediction value y (line 7b) and the prediction value y±2σ (y+2σ is the line 7a and y−2π is the line 7c) are as illustrated in
As described above, various prediction models and prediction errors illustrated in
The lines C1 to C3 illustrated in
Next, the flow of a prediction model construction process of the present embodiment will be described.
First, the data collection unit 31 acquires process data including the values of an input variable and an output variable required to construct a prediction model (step S11). Next, the data storage 32 stores the acquired process data in the storage unit 37 (step S12). Next, the data extraction unit 33 extracts and reads out process data, which is required for a predetermined prediction model, from the storage unit 37 to output the extracted process data to the prediction model construction unit 34 (step S13). The prediction model construction unit 34 sets the input variable and the output variable from the extracted process data (step S14). The prediction model construction unit 34 uses a technique such as multiple regression analysis, random forest regression, Gaussian process regression, and a neural network to construct, for example, the prediction model and the error calculation model illustrated in
First, the data collection unit 31 acquires process data of an evaluation target including a predetermined input variable (step S21). The data storage 32 stores the acquired process data in the storage unit 37. Next, the data extraction unit 33 extracts and reads out process data, which is required to calculate a prediction value, from the storage unit 37. Next, the prediction unit 35 reads out a predetermined prediction model which predicts a prediction value of the evaluation target, and a predetermined error calculation model from the storage unit 37. The prediction unit 35 inputs the process data into the prediction model to calculate the prediction value (step S22). The prediction unit 35 inputs the prediction value and the process data into the error calculation model to calculate a prediction error (step S23). When the technique of constructing the prediction model is Gaussian process regression, the process data is input into the prediction model, so that the prediction value corresponding to the value of the process data and the prediction error can be obtained at the same time. The prediction unit 35 adds or subtracts the prediction error to or from the prediction value to calculate a final prediction value (step S24). At this time, the prediction unit 35 may output the prediction value before correction and the prediction error in addition to the final prediction value.
According to this embodiment, a model which calculates a prediction error together with a prediction value can be constructed from a plurality of process data (learning data). Process data of an evaluation target is input into a constructed model, and thus even when the influence of a prediction error of the prediction model is considered to the maximum extent, a corrected prediction value (final prediction value) which allows safe operation of the plant or the like can be obtained. In addition, there is no need for constructing a plurality of prediction models in order to obtain one prediction value.
A prediction device 30A of a second embodiment is a so-called operation guidance device that determines whether or not a prediction value predicted by the prediction unit 35 is within a predetermined allowable range, and when the prediction value is within the allowable range, provides an operation amount to converge the prediction value within the allowable range or information which guides determination of such an operation amount.
Among configurations according to the second embodiment of the present invention, the same functional units as those forming the prediction device 30 according to the first embodiment of the present invention are denoted by the same reference signs, and description thereof will be omitted. The prediction device 30A according to the second embodiment includes a state monitoring unit 38 and an operation-amount determination unit 39 in addition to the configuration of the first embodiment.
The state monitoring unit 38 monitors process data. Specifically, the state monitoring unit 38 compares the process data with a threshold value set for each of the process data to determine that there is abnormality, when the process data deviates from the threshold value. The threshold value used for the determination may be set based on a prediction model constructed by the prediction model construction unit 34. The state monitoring unit 38 may perform threshold determination on a final prediction value as a monitoring target, the final prediction value being predicted by the prediction unit 35 based on the process data.
When the state monitoring unit 38 determines that there is abnormality, the operation-amount determination unit 39 determines an operation amount or a control value of the plant or the mechanical device to avoid the abnormality. For example, when the level of combustion vibration is high, the operation-amount determination unit 39 determines an operation amount that lowers the level of combustion vibration (for example, how much the opening degree of the fuel flow regulation valve 16A is reduced or increased, or the like). For example, when the amount of emission of NOx or CO is large, the operation-amount determination unit 39 determines an operation amount that reduces the amount of emission thereof. For example, when the output efficiency of the gas turbine 10 is low, the operation-amount determination unit 39 determines an operation amount that improves the output efficiency. As will be described later, when an operation amount is determined, a prediction model which has the operation amount or process data related to the operation amount (for example, the flow rate of fuel supplied from the main system when the operation amount is the opening degree of the fuel flow regulation valve 16A) as an explanatory variable (input variable) can be used.
The output unit 36 outputs the operation amount, which is determined by the operation-amount determination unit 39, to the device 20. Alternatively, the output unit 36 displays the operation amount, which is determined by the operation-amount determination unit 39, on a display of the prediction device 30A or the like.
First, the data collection unit 31 acquires process data of an evaluation target including a predetermined input variable (step S31). The data collection unit 31 outputs the process data to the state monitoring unit 38. The state monitoring unit 38 compares each of a plurality of the process data with a corresponding threshold value (step S32). When there is process data deviating from the threshold value (step S33: Yes), the state monitoring unit 38 notifies the operation-amount determination unit 39 of a detection of abnormality. The operation-amount determination unit 39 acquires the process data (process data acquired in step S31) including an input variable for which the abnormality is detected, to determine a safe operation amount (step S34). For example, the operation-amount determination unit 39 instructs the prediction unit 35 to output a final prediction value. The prediction unit 35 inputs the process data, which is acquired in step S31, into an input model to output the prediction value. Here, reference will be made to
When the process data is determined to be within the threshold value in step S33, the process is repeated from step S31 for the next process data.
A process in which the prediction device 30A displays guidance information that guides determination of an operation amount which improves the operating state, instead of the operation amount, will be described with reference to
According to the present embodiment, in addition to the effects of the first embodiment, the plant or the mechanical device can be stably operated by the operation amount that is determined by the operation-amount determination unit 39 based on the prediction value on a safe side based on uncertainty of the prediction model, the prediction value being predicted by the prediction unit 35, or guidance information output by the output unit 36.
A computer 900 is, for example, a personal computer (PC) or a server terminal device including a CPU 901, a main storage device 902, an auxiliary storage device 903, an input and output interface 904, and a communication interface 905. The prediction devices 30 and 30A described above are implemented in the computer 900. Then, the operation of each of the processing units described above is stored in the auxiliary storage device 903 in the form of a program. The CPU 901 reads out the program from the auxiliary storage device 903 to expand the program in the main storage device 902, and to then execute the above processes according to the program. The CPU 901 secures a storage area, which corresponds to the storage unit 37, in the main storage device 902 according to the program. The CPU 901 secures a storage area, which stores data being processed, in the auxiliary storage device 903 according to the program.
In at least one embodiment, the auxiliary storage device 903 is one example of a non-transitory medium. Other examples of the non-transitory medium include a magnetic disk, a magneto-optical disk, a CD-ROM, a DVD-ROM, a semiconductor memory, and the like that are connected via the input and output interface 904. When the program is delivered to the computer 900 by a communication line, the computer 900 which receives the delivery may expand the program in the main storage device 902 to execute the above processes. The program may realize a part of the above-described functions. Further, the program may be a so-called differential file (differential program) that realizes the above-described functions in combination with another program already stored in the auxiliary storage device 903.
In addition, well-known components can be appropriately replaced with the components in the embodiments without departing from the spirit of the present invention. The technical scope of the present invention is not limited to the embodiments, and various modifications can be made without departing from the spirit of the present invention. The output unit 36 is one example of a first output unit and a second output unit.
According to the prediction device, the prediction method, and the program, the prediction value based on the influence of the prediction error of the prediction model can be output.
Number | Date | Country | Kind |
---|---|---|---|
2018-156568 | Aug 2018 | JP | national |
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/JP2019/029006 | 7/24/2019 | WO | 00 |