Patent Application
20220317646
The present disclosure relates to the field of electrical equipment technologies, and in particular, to a self-adaptive test method for an intelligent prediction algorithm of analog measured values.
A power plant attendant is required to monitor more than 500 analog measurement points with a trip outlet function to monitor a panel. The analog measurement points with the trip outlet function are distributed in different simulation diagrams of a master computer of a monitoring system. It is impossible to find an abnormality in advance by manually tracking a change trend by an attendant. At present, intelligent technologies such as machine learning and the like are increasingly developed, and the intelligent technologies provide important technical means for trend judgment and prediction.
However, it is not negligible that such intelligent technologies rely on fault samples, and a large number of fault samples are required to participate in an algorithm training process to ensure a good implementation effect. Fault sample data is less in the operation and maintenance of a power plant with a mature technology. In addition, there is a lack of a test method for verifying implementation effects of the intelligent technologies. Application effects of the intelligent technologies stay in subjective fuzzy understanding, and even stay in the disposal of accidental events.
In the past, the acquisition of fault samples comes from debugging and troubleshooting after overhaul, and the acquisition of the fault samples is costly. In addition, numerical characteristics of the fault samples cannot completely cover trend characteristics under fault conditions, which have some limitations. Therefore, means for testing performance of an intelligent trend judgment algorithm and a method for providing fault samples for the intelligent trend judgment algorithm are urgently needed at present.
The present disclosure provides a self-adaptive test method for an intelligent prediction algorithm of analog measured values, which may test an implementation effect of an intelligent prediction algorithm of analog measured values of a start-up state of a unit by integrating historical operating conditions, alarm thresholds and switching value signals, so as to automatically perform self-adaptive adjustment according to the historical operating conditions of the unit for testing, so that the obtaining of fault samples and the testing of the implementation effect of the intelligent technology can be solved better, faster and more economically through a computer.
The present disclosure provides a self-adaptive test method for an intelligent prediction algorithm of analog measured values, including the following steps:
(1.1) reading an event recording sequence, an analog measurement point ID and an analog measurement point alarm value from a time sequence event record table, an analog measurement point table and an alarm threshold table;
(1.2) acquiring operation records of a normal operation state of a unit within a statistical cycle to form historical statistics of measured values of an analog measurement point based on switching value signals;
(1.3) calculating simulated measured values of the analog measurement point with time scales by integrating the historical statistics, the analog measurement point alarm value and an analog measurement point current measured value;
(1.4) providing the simulated measured values of the analog measurement point with time scales for the intelligent prediction algorithm of analog measured values of a start-up state of the unit for testing, and calculating sensitivity; and
(1.5) sending an alarm to remind a technician to adjust the algorithm when the sensitivity is less than a threshold.
In the method, the time sequence event record table is a switching value signal set K with a time record, a state record and equipment description set in sequence by a technician, the switching value signal set K including at least a unit starting command signal, a unit steady-state signal and a unit-load to base-load signal; the analog measurement point table is a to-be-tested analog measurement point ID set M set by an attendant; and the alarm threshold table is an analog measurement point first-level alarm set B1 and an analog measurement point second-level alarm set B2.
In the method, the historical statistics of the measured values of the analog measurement point based on the switching value signals is obtained through the following steps:
(3.1) traversing switching value records in the statistical cycle, taking out switching value signals simultaneously satisfying the switching value signal set K in sequence, and storing a time of the switching value signals taken out in a time sequence TL according to the sequence of the switching value signal set K;
(3.2) traversing analog records of the analog measurement point ID set M in the statistical cycle, and taking out a maximum value of measured values of measurement point of the analog measurement point set M with a time scale of the time sequence TL to obtain a measured value set CL of the measurement point; and
(3.3) the measured value set CL of the measurement point being the historical statistics of the measured values of the analog measurement point based on the switching value signals.
In the method, the calculating simulated measured values of the analog measurement point with time scales by integrating the historical statistics, the analog measurement point alarm value and an analog measurement point current measured value involves the following steps:
(4.1) acquiring an average value maxave from the measured value set CL of the measurement point, and acquiring, from the alarm threshold table, a first-level alarm value b1 and a second-level alarm value b2 of a measurement point of the present unit;
(4.2) calculating a stacking slope k1, the stacking slope k1=b1/maxave−1, and calculating a stacking slope k2, the stacking slope k2=b2/maxave−1;
(4.3) setting i=1, a test sample number being yb, and a tolerance d=(k2−k1)/yb;
(4.4) calculating tmp1=k1+(i−1)×d;
(4.5) when i is not greater than yb, storing tmp1 in a vector k(i), i=i+1, performing step (4.4); when i is greater than yb, performing step (4.6); and
(4.6) adding 0.01 to each element in the vector k(i) to obtain a superposition stacking slope vector k;
(4.7) acquiring a time T1 of the unit starting command signal and a time T3 of the unit-load to base-load signal in the switching value signal set K;
(4.8) setting i=1;
(4.9) acquiring an actually measured analog measured value with the time scale as x(t), before the time T1, y(t)=x(t), between T1 and T2, y(t)=x(t)×(1+k(i)), and after the time T3, y(t)=x(t)+(y(T3)−x(T3));
(4.10) when i is not greater than yb, transposing a column vector y to a matrix yy, i=i+1, and performing step (4.9); when i is greater than yb, performing step (4.11); and
(4.11) the matrix yy being the simulated measured values with time scale of the analog measurement point.
In the method, the sensitivity is calculated through the following steps:
(5.1) acquiring the test sample number yb;
(5.2) acquiring the time T1 of the unit starting command signal, the time T2 of the unit steady-state signal and the time T3 of the unit-load to base-load signal in the switching value signal set K, i=1, and f=0;
(5.3) providing the simulated measured values yy(i,:) of the analog measurement point with time scales t for the intelligent prediction algorithm of analog measured values of the start-up state of the unit for testing, and if a testing algorithm finds a fault before the time T2 f=f+1;
(5.4) when i is not greater than yb, i=i+1, performing step (5.3); when i is greater than yb, performing step (5.5); and
(5.5) the sensitivity L=f/yb×100%.
In the method, the threshold δ is 50%.
Compared with the prior art, the present disclosure fills the blank of the engineering field and has the following beneficial effects.
(1) The present disclosure provides a standardized test method for an intelligent prediction algorithm of analog measured values of a start-up state of a unit. The test method according to the present disclosure can completely cover trend characteristics under fault conditions, and break the limitations of the original dependence on actually measured fault samples.
(2) The present disclosure may automatically perform self-adaptive adjustment according to historical operating conditions of the unit to meet test requirements, so that the obtaining of fault samples and the testing of the implementation effect of the intelligent technology can be solved better, faster and more economically through a computer.
(3) The present disclosure quantitatively evaluates the implementation effect of the intelligent technology, provides index support for parameter adjustment of the intelligent technology, the selection of a proper algorithm and implementation means, also realizes the advance of the verification of the implementation effect of the intelligent technology, and prevents losses and uncertain influence caused by the verification through actual engineering.
Specific implementations of the present disclosure are further described below with reference to the accompanying drawings and examples, but implementation and protection of the present disclosure are not limited thereto. It is to be noted that if any of the following processes is not described in detail, they may be realized or understood by those skilled in the art with reference to the prior art.
The present disclosure performs standardization in combination with engineering experience, and provides a test method for testing an intelligent prediction algorithm of analog measured values of a start-up state of unit by integrating historical operating conditions, alarm thresholds and switching value signals. In addition, the present disclosure may automatically perform self-adaptive adjustment according to historical operating conditions of the unit to meet test requirements, so that the obtaining of fault samples and the testing of the implementation effect of the intelligent technology can be solved better, faster and more economically through a computer.
The following is an example analysis on monitoring signals for starting of a power generation condition of a #4 unit in Guangzhou Energy Storage Hydropower Plant from 11:00 to 19:00 on Apr. 10, 2019.
With reference to the flow in
(1.1) An event recording sequence, an analog measurement point ID and an analog measurement point alarm value are read from a time sequence event record table (Table 1 below), an analog measurement point table (Table 2 below) and an alarm threshold table (Table 3 below). The time sequence event record table is a switching value signal set K with a time record, a state record and equipment description set in sequence by a technician. The switching value signal set K includes at least a unit starting command signal, a unit steady-state signal and a unit-load to base-load signal. The analog measurement point table is a to-be-tested analog measurement point ID set M set by an attendant. The alarm threshold table is an analog measurement point first-level alarm set B1 and an analog measurement point second-level alarm set B2.
(1.2) Operation records of a normal operation state of a unit within a statistical cycle are acquired to form historical statistics of measured values of the analog measurement point based on switching value signals. In the present embodiment, the statistical cycle is the past half year.
The historical statistics is obtained through the following steps. Switching value records in the past half year are traversed, switching value signals simultaneously satisfying the switching value signal set K in sequence are taken out, and the time of the switching value signals taken out is stored in a time sequence TL according to the sequence of the switching value signal set K.
Analog records of the analog measurement point ID set M in the past half year are traversed, and a maximum value of measured values of the measurement point of the analog measurement point set M with a time scale of the time sequence TL is taken out to obtain a measured value set CL of the measurement point.
The measured value set CL of the measurement point is the historical statistics of the point measured values of the analog measurement (Table 4 below) based on the switching value signals.
(1.3) Simulated measured values of the analog measurement point with time scales are calculated by integrating the historical statistics, the analog measurement point alarm value and an analog measurement point current measured value, which involves the following specific steps.
(1.3a) An average value maxave is acquired from the measured value set CL of the measurement point, and a first-level alarm value b1 and a second-level alarm value b2 of a measurement point of the present unit are acquired from the alarm threshold table.
(1.3b) A stacking slope k1 is calculated, the stacking slope k1=b1/maxave−1, and a stacking slope k2 is calculated, the stacking slope k2=b2/maxave−1.
In the present embodiment, the calculated stacking slope k1=b1/maxave−1=75/60.7−1=0.2356, and the stacking slope k2=b2/maxave−1=80/60.7−1=0.3180.
(1.3c) i=1, i denotes an ith out-of-limit test sample of the measurement point, a test sample number is yb, and a tolerance d=(k2−k1)/yb.
(1.3d) tmp1=k1+(i−1)×d is calculated, and tmp1 indicates that a measured value of the measurement point can be out of limit by superposing a change speed of tmp1 on the basis of a maximum average measured value of the measurement point.
(1.3e) When i is not greater than yb, tmp1 is stored in a vector k(i), the vector k(i) denotes a superposition change speed of the measured value of the ith out-of-limit test sample of the measurement point, i=i+1, and step (4.4) is performed; when i is greater than yb, step (4.6) is performed.
i=1, tmp1=k1+(i−1)×d=0.2356, k(1)=0.2356;
i=2, tmp1=k1+(i−1)×d=0.2356+1×0.0165, k(2)=0.2521;
i=3, tmp1=k1+(i−1)×d=0.2356+2×0.0165, k(3)=0.2521;
i=4, tmp1=k1+(i−1)×d=0.2356+3×0.0165, k(4)=0.2521;
i=5, tmp1=k1+(i−1)×d=0.2356+4×0.0165, k(5)=0.3015;
k(i)=[0.2356, 0.2521, 0.2521, 0.2521, 0.30151]T
(1.3f) 0.01 is added to each element in the vector k(i) to obtain a superposition stacking slope vector k, i∈[1,yb].
k=k(i)+0.01
k=[0.2456, 0.2621, 0.2785, 0.2950, 0.31151]T
(1.3g) Time T1 of the unit starting command signal and time T3 of the unit-load to base-load signal in the switching value signal set K are acquired. In the present embodiment, the time of the unit starting command signal is T1=2019-04-10 11:13:22, which corresponds to the numerical time 737525.46 in
(1.3h) i=1.
(1.3i) An actually measured analog measured value with the time scale is acquired as x(t); before time T1, y(t)=x(t); between T1 and T2, y(t)=x(t)×(1+k(i)); and after the time T3, y(t)=x(t)+(y(T3)−x(T3)). T2 represents time of the unit steady-state signal. The time of the unit steady-state signal in the present embodiment is T2=2019-04-10 11:17:06, which corresponds to the numerical time 737525.47 in
(1.3j) When i is not greater than yb, a column vector y is transposed to a matrix yy, i=i+1, and step (4.9) is performed; when i is greater than yb, step (4.11) is performed.
(1.3k) The matrix yy(i,:) is the simulated measured values of the analog measurement point with time scales.
When i=1,
y(T3)=75.53; x(T3)=60.7;
before the time T1, y(t)=x(t),
between T1 and T2, y(t)=x(t)×(1+0.2456)
after the time T3, y(t)=x(t)+14.833.
yy(1,:) is a combination of each y(t) after T1 to T3.
When i=2,
y(T3)=76.53; x(T3)=60.70;
before the time T1, y(t)=x(t) prior to,
between T1 and T2, y(t)=x(t)×(1+0.2621)
after the time T3, y(t)=x(t)+15.828.
yy(2,:) is a combination of each y(t) after T1 to T3.
When i=3,
y(T3)=77.52; x(T3)=60.70;
before the time T1, y(t)=x(t),
between T1 and T2, y(t)=x(t)×(1+0.2785)
after the time T3, y(t)=x(t)+16.823.
yy(3,:) is a combination of each y(t) after T1 to T3.
When i=4,
y(T3)=78.52; x(T3)=60.70;
before the time T1, y(t)=x(t)
between T1 and T2, y(t)=x(t)×(1+0.2950)
after the time T3, y(t)=x(t)+17.818.
yy(4,:) is a combination of each y(t) after T1 to T3.
When i=5,
y(T3)=79.51; x(T3)=60.70;
before the time T1, y(t)=x(t),
between T1 and T2, y(t)=x(t)×(1+0 0.3115)
after the time T3, y(t)=x(t)+18.814.
yy(5,:) is a combination of each y(t) after T1 to T3.
(1.4) The simulated measured values of the analog measurement point with time scales are provided for an intelligent prediction algorithm of measured values of a start-up state analog of the unit for testing, and sensitivity is calculated.
As shown in
In the present embodiment, a fault is found prior to the time T2, where f=4 and yb=5, and the sensitivity is L=f/yb×100%=80%.
(1.5) An alarm is sent to remind a technician to adjust the algorithm when the sensitivity is less than a threshold. In the present embodiment, the threshold δ is 50%.
In the present embodiment, the sensitivity L is not less than the threshold, and there is no need to send an alarm to remind the technician to adjust the algorithm.
Therefore, the present disclosure provides a standardized test method for testing an intelligent prediction algorithm of analog measured values of a start-up state of a unit. The testing according to the present disclosure can completely cover trend characteristics under fault conditions, and break the limitations of the original dependence on actually measured fault samples. Self-adaptive adjustment may be automatically performed according to historical operating conditions of the unit so as to meet test requirements, so that the work of obtaining fault samples and testing the implementation effect of the intelligent technology can be solved better, faster and more economically through a computer. At the same time, the implementation effect of the intelligent technology is also quantitatively evaluated, which provides index support for parameter adjustment of the intelligent technology, the selection of a proper algorithm and implementation means, also realizes the advance of the verification of the implementation effect of the intelligent technology, and prevents losses and uncertain influence caused by the verification through actual engineering.
| Number | Date | Country | Kind |
|---|---|---|---|
| 202010457983.7 | May 2020 | CN | national |
| Filing Document | Filing Date | Country | Kind |
|---|---|---|---|
| PCT/CN2020/119659 | 9/30/2020 | WO |
