Parameter processing method and device for compressed air energy storage device

Abstract

The invention discloses a parameter processing method and device for a compressed air energy storage device, belongs to the technical field of electric digital data processing, and aims at attributes of pressure values and flow values, the relation between the pressure values and the flow values is analyzed in detail, interference variables of all the pressure values in a selected numerical cluster are accurately determined through compliant analysis operation, and the pressure values and the flow values are accurately determined. According to the interference variable, a special interference removing method is constructed for each pressure value, and the method comprises the steps of mechanically determining the cleaning frequency and the unit capacity under different cleaning frequencies, so as to achieve flexible active cleaning, so that the condition that the high pressure value is mistakenly removed can be efficiently avoided, and the cleaning efficiency is improved. In this way, a more accurate cleaning value is obtained, the precision and correctness of future time point pressure value estimation are improved, and accurate detection of the parameter condition at the inlet of the turbine can be achieved.

Description

Technical Field

[0001] The present invention belongs to the technical field of electrical digital data processing, and particularly relates to a parameter processing method and device for a compressed air energy storage device. Background Art

[0002] Compressed-Air Energy Storage (CAES) refers to a method of energy storage in which electrical energy is used to compress air during the low-load period of the power grid, and the compressed air is sealed under high pressure in abandoned mines, sunken undersea gas storage tanks, caves, expired oil and gas wells, or newly built gas storage wells, and the compressed air is released during the peak-load period of the power grid to drive a turbine for power generation.

[0003] During the actual operation of a compressed air energy storage device, as mentioned in the prior art solution with the patent publication number "CN107769233B" and the patent name "A small-capacity compressed air energy storage device and its grid connection method", a pressure sensor is often set at the inlet of its turbine to detect the pressure value at the inlet of the turbine.

[0004] On the other hand, the pressure value sampled by the pressure sensor is often subject to external electromagnetic interference and is prone to generating interference values. Therefore, it is necessary to clean the pressure value to remove the interference value. Currently, the method for cleaning the pressure value is often the Gaussian filtering method. However, since the cleaning method of the Gaussian filtering method defines a fixed unit capacity for all the pressure values to be cleaned, the unit capacity is the number of adjacent sampled pressure values after the sampling time point of the defined pressure value, and the interference removal method of the Gaussian filtering method uses a fixed unit capacity for all the pressure values to be cleaned. Such a fixed processing method often causes misidentification when facing high pressure values at the inlet of the turbine. Because high pressure values are often very similar to interference values, they are often misidentified as interference values and removed. This not only weakens the accuracy and integrity of the pressure value, but also causes a large error between the estimated value of the pressure value at a future time point derived from the estimation model constructed based on such pressure values and the actual situation, and cannot accurately detect the parameter status at the inlet of the turbine. Summary of the Invention

[0005] To solve the defects in the prior art, the present invention proposes a parameter processing method and device for a compressed air energy storage device, effectively avoiding the defects of the Gaussian filtering method in the prior art, which often causes misidentification by using a fixed unit capacity for all the pressure values to be cleaned at the inlet of the turbine of the compressed air energy storage device, weakens the accuracy and integrity of the pressure value, causes a large error between the estimated value of the pressure value at a future time point derived from the estimation model constructed based on such pressure values and the actual situation, and cannot accurately detect the parameter status at the inlet of the turbine.

[0006] The present invention employs the following technical solutions.

[0007] A parameter processing method for a compressed air energy storage device, comprising: A pressure sensor and a flow sensor respectively sample the pressure value and the flow value of the gas at the inlet of the turbine of the compressed air energy storage device, and transmit the pressure value and the flow value to an industrial control computer for processing; The method of transmitting the pressure value and the flow value to an industrial control computer for processing includes: Step 1, receiving the pressure value and the flow value at the inlet of the turbine during a defined period sampled; Step 2, performing several cleanings on all the pressure values, and actively adjusting the unit capacity of the pressure value before each cleaning; Step 3, performing cleaning according to the adjusted unit capacity of all the pressure values, and detecting the operating condition at the inlet of the turbine according to all the pressure values after cleaning.

[0008] Preferably, in Step 1, the pressure value and the flow value are respectively sampled by the pressure sensor and the flow sensor, the starting sampling time points and the sampling speeds of the pressure value and the flow value are the same, the sampling speed is 2 s / time, the sampling time of the pressure value and the flow value is defined as 30 min, and the sampling time is used as the defined period.

[0009] Preferably, in Step 2, the method of actively adjusting the unit capacity of the pressure value includes: Step 2-1, obtaining the final interference amount of each pressure value; Step 2-2, determining whether the pressure value needs to be included in this cleaning according to the level of the final interference amount of each pressure value. If it is necessary, actively adjust the unit capacity of the pressure value by using the final interference amount of each pressure value. If it is not necessary, this pressure value does not perform this cleaning.

[0010] Preferably, Step 2-1 specifically includes: Step 2-1-1, estimating the initial interference amount of the pressure value through the error between the pressure value and the average of its adjacent numerical clusters, and combining the irrelevance between the pressure value and the flow value corresponding to its sampling time point; Step 2-1-2, calculating the final interference amount of each pressure value according to the difference between the initial interference amount of each pressure value and the initial interference amount of the pressure values within its adjacent numerical clusters, and combining the confidence level of the difference.

[0011] Preferably, in Step 2-1-1, the method for estimating the initial interference amount includes: selecting the adjacent numerical cluster of the pressure value, that is, taking the pressure value as the midpoint, and selecting pressure values before and after its sampling time point respectively, and jointly forming the adjacent numerical cluster of the pressure value with the pressure value taken as the midpoint; The operation equation for the starting interference amount of each pressure value is as follows: ; In the equation, is the starting interference amount of the th pressure value, is the value of the th pressure value, is the mean of the near-domain numerical cluster of the th pressure value, is the irrelevance between the th pressure value and the flow value corresponding to the sampling time point of the th pressure value, is to perform standardization on using the Z-score standardization method.

[0012] Preferably, in step 2-1-1, the irrelevance between the th pressure value and the flow value corresponding to the sampling time point of the th pressure value has the following operation equation: ; In the equation, is the sequence code ratio of the th pressure value within its near-domain numerical cluster, is the sequence code ratio of the flow value corresponding to the sampling time point of the th pressure value within its near-domain numerical cluster, is the defined additional value.

[0013] Preferably, step 2-1-1 further includes: estimating the starting interference amount of the flow value through the error between each flow value and the mean of its near-domain numerical cluster, plus the irrelevance between the flow value and the pressure value corresponding to its sampling time point. The operation equation for the starting interference amount of the flow value is as follows: ; In the equation, is the starting interference amount of the th flow value, is the value of the th flow value, is the mean of the near-domain numerical cluster of the th flow value, is the irrelevance between the th flow value and the pressure value corresponding to the sampling time point of the th flow value, is to perform standardization on using the Z-score standardization method.

[0014] Preferably, in step 2-1-2, the operation equation for the final interference amount of the pressure value is: ; In the equation, is the final interference amount of the th pressure value, is the number of pressure values outside the pressure value itself within the near-domain numerical cluster of the th pressure value, is the sequence code of the pressure values outside the pressure value itself within the near-domain numerical cluster of the th pressure value, is the starting interference amount of the th pressure value, is the starting interference amount of the th pressure value outside the pressure value itself within the near-domain numerical cluster of the th pressure value, represents the difference between the starting interference amount of the th pressure value and the starting interference amount of the th pressure value within its near-domain numerical cluster, is the credible amount of , and its operation equation is: ; In the equation, is the outlier amount of the flow value corresponding to the sampling time point of the th pressure value, is the time interval between the th pressure value and the th pressure value outside the pressure value itself within its near-domain numerical cluster.

[0015] Preferably, in step 2-2, to define the critical value of the final interference amount, before each cleaning of each pressure value, compare the final interference amount of the pressure value with the critical value of the final interference amount. If the final interference amount of the pressure value is higher than the critical value of the final interference amount, there is a need to add this cleaning. Conversely, if the final interference amount of the pressure value is not higher than the critical value of the final interference amount, there is no need to add this cleaning for this pressure value.

[0016] Preferably, in step 2-2, the operation equation for actively adjusting the unit capacity is: ; In the equation, is the Euler number, is the current cleaning frequency of the th pressure value, is the th pressure value when adding the During the secondary cleaning, the unit capacity after active adjustment is the unit capacity predefined for all pressure values, is the th pressure value and the last interference amount before the th cleaning, is the th pressure value and the last interference amount before the th cleaning.

[0017] Preferably, in step 3, the method for detecting the operating condition at the inlet of the turbine based on all the pressure values after cleaning includes: Regarding all the pressure values after cleaning as past-period pressure values, then regarding the past-period pressure values as dependent variables and the corresponding sampling time points of the past-period pressure values as independent variables, using the least squares method to obtain the regression equation of the past-period pressure values, and obtaining the estimated values of the pressure values at each sampling time point after the defined period according to the regression equation; Based on the difference between the estimated values and the measured values of the pressure values at each time point after the defined period, determine whether the condition at the inlet of the turbine is abnormal, and thus achieve the detection of the condition at the inlet of the turbine. The method is as follows: Calculate the difference between the measured values and the estimated values of the pressure values at each time point after the defined period; if the difference between the measured value and the estimated value of the pressure value at a certain time point is higher than the predefined difference critical value, the pressure value at that time point is abnormal; if the difference does not exceed the difference critical value, the pressure value at that time point is not abnormal; Sum up the abnormal conditions of the pressure values at each time point after the defined period. After the defined period, starting from the time point when the first abnormal pressure value appears, sum up the number of abnormal pressure values within the defined time; if the ratio of the number of abnormal pressure values to the total number of all pressure values within the defined time exceeds the defined ratio critical value, the condition at the inlet of the turbine is in an abnormal state; if the ratio does not exceed the defined ratio critical value, the condition at the inlet of the turbine is in a normal state.

[0018] A parameter processing device for a compressed air energy storage device, including: A pressure sensor and a flow sensor are provided at the inlet of the turbine of the compressed air energy storage device, and both the pressure sensor and the flow sensor are connected to the industrial control computer; The modules running on the industrial control computer include: A receiving module, which is used to receive the pressure values and flow values at the inlet of the turbine during the defined period transmitted by sampling; An adjustment module, which is used to perform several cleanings on all the pressure values, and actively adjust the unit capacity of the pressure values before each cleaning; A detection module, which is used to perform cleaning according to the unit capacity after adjusting all pressure values, and detect the operating conditions at the inlet of the turbine according to all pressure values after cleaning.

[0019] The beneficial effects of the present invention are that, compared with the prior art, the technical effects of the present invention include: Regarding the attributes of the pressure value and the flow value, carefully analyze the relationship between the pressure value and the flow value. Through compliant analysis operations, accurately determine the interference amount of each pressure value within the selected numerical cluster. According to this interference amount, construct a special anti-interference method for each pressure value, which includes automatically determining its cleaning frequency and the level of the unit capacity at different cleaning frequencies, achieving automatic active cleaning. Through such methods, the situation of high pressure values being wrongly removed can be efficiently avoided, thereby obtaining more accurate cleaning values, improving the accuracy and correctness of the pressure value estimation for future time points, and enabling accurate detection of the parameter conditions at the inlet of the turbine. Description of the Drawings

[0020] Figure 1 is a flowchart of the parameter processing method for the compressed air energy storage device described in the present invention; Figure 2 is a module structure diagram of the industrial control computer described in the present invention. Detailed Embodiments

[0021] To make the objectives, technical solutions, and advantages of the present invention clearer, the following will clearly and completely describe the technical solutions of the present invention with reference to the accompanying drawings in the embodiments of the present invention. The embodiments described herein are only some of the embodiments of the present invention, not all of them. According to the spirit of the present invention, other embodiments obtained by those skilled in the art without creative efforts fall within the protection scope of the present invention.

[0022] As Figure 1 shown, a parameter processing method for a compressed air energy storage device described in the present invention includes: A pressure sensor and a flow sensor respectively sample the pressure value and the flow value of the gas at the inlet of the turbine of the compressed air energy storage device, and transmit the pressure value and the flow value to an industrial control computer for processing; The method of transmitting the pressure value and the flow value to the industrial control computer for processing runs on the industrial control computer and includes: Step 1, receive the pressure value and the flow value at the inlet of the turbine during a defined period sampled and transmitted; The industrial control computer receives the pressure value and the flow value at the inlet of the turbine of the compressed air energy storage device sampled and transmitted. These two parameters present the main parameter conditions at the inlet of the turbine.

[0023] Therefore, by carefully analyzing the relationship between the pressure value and the flow value, for example, when the flow rate increases, the pressure usually also increases. Such a relationship can assist in estimating the conditions at the inlet of the turbine under different operating conditions, thereby accurately performing parameter condition detection. Long-term detection of such pressure values and flow values can estimate the trend of parameter changes at the inlet of the turbine and pre-determine the maintenance plan for the turbine.

[0024] In a preferred but non-limiting embodiment of the present invention, in step 1, the pressure value and the flow value are respectively obtained by sampling through a pressure sensor and a flow sensor. The starting sampling time points and the sampling speed of the pressure value and the flow value are the same, and the sampling speed is 2 s / time. Define the sampling time of the pressure value and the flow value as 30 min, and regard the sampling time of the first time as the defined period. Therefore, a pressure value and a flow value can be obtained at each sampling time point, ensuring the temporal coordination of all the sampled pressure values and flow values, which is conducive to later analysis of the pressure value and the flow value.

[0025] Step 2, perform several cleanings on all the pressure values. Before each cleaning, actively adjust the unit capacity of the pressure value; the unit capacity is the number of pressure values of the adjacent sampling after the sampling time point of the defined pressure value. , performing cleaning on the pressure value means sending the pressure value and the pressure values adjacent to it after the sampling time point into the Gaussian filtering method to perform cleaning on the pressure value. If the number of pressure values adjacent to it after the sampling time point is insufficient , obtain the adjacent pressure values sampled before the sampling time point of the pressure value to make up pressure values, and then send the pressure value and the obtained pressure values into the Gaussian filtering method to perform cleaning on the pressure value. The cleaning method is the Gaussian filtering method. The number of times of several times can be determined according to specific requirements.

[0026] In a preferred but non-limiting embodiment of the present invention, in step 2, the method for actively adjusting the unit capacity of the pressure value includes: Step 2-1, obtain the final interference amount of each pressure value; In a preferred but non-limiting embodiment of the present invention, step 2-1 specifically includes: Step 2-1-1, estimate the initial interference amount of the pressure value through the error between each pressure value and the average of its adjacent numerical clusters, combined with the irrelevance between the pressure value and the flow value corresponding to its sampling time point; In a preferred but non-limiting embodiment of the present invention, in step 2-1-1, the estimation method of the initial interference amount includes: select the adjacent numerical cluster of the pressure value, that is, regard the pressure value as the midpoint, and select A pressure value, together with the pressure value taken as the midpoint, jointly forms a cluster of neighboring numerical values of this pressure value. is a defined value, and its value can be , that is, this pressure value, the pressure values sampled at the previous sampling time points adjacent to the sampling time point of this pressure value, and the pressure values at the subsequent sampling time points adjacent to the sampling time point of this pressure value, a total of two hundred and one pressure values are taken as the cluster of neighboring numerical values of this pressure value; and if there are no pressure values before a pressure value, then the insufficient number is selected after this pressure value. Just like if there are fifty pressure values before this pressure value and there are still fifty insufficient ones, then fifty pressure values are taken successively after this pressure value. A total of one hundred and fifty pressure values are obtained after this pressure value. If the number of pressure values after this pressure value is less than one hundred, similarly, the insufficient number of pressure values is supplemented before this pressure value.

[0027] The calculation equation for the starting interference amount of each pressure value is: ; In the equation, is the starting interference amount of the th pressure value, is the value of the th pressure value, is the mean of the cluster of neighboring numerical values of the th pressure value, is the irrelevance between the th pressure value and the flow value corresponding to the sampling time point of the th pressure value, is to perform standardization on using the Z-score standardization method. is the error between the pressure value and the mean of its cluster of neighboring numerical values.

[0028] In this equation, represents the difference between the th pressure value and the mean of the pressure values of its cluster of neighboring numerical values. The higher is, the more significant the th pressure value is. It is more likely to be a true abnormal value of high pressure caused by changes in the parameters of the turbine inlet itself, or it may also be an interference value. To more precisely identify the specific situation, it is necessary to analyze the The irrelevance between a pressure value and the corresponding flow rate value at its sampling time point is based on the rule that an increase in the gas flow rate at the inlet of the turbine will cause an increase in the gas pressure. Therefore, when the flow rate increases, it often causes an increase in the gas pressure at the inlet of the turbine. If the correlation between a pressure value and the corresponding flow rate value at its sampling time point is lower (higher irrelevance), it indicates that this pressure value is more likely to be an interference value rather than a true abnormal value caused by changes in the parameters of the turbine inlet itself.

[0029] In a preferred but non-limiting embodiment of the present invention, in step 2-1-1, the th pressure value and the th pressure value and the irrelevance between the corresponding flow rate values at their sampling time points The calculation equation is: ;

[0030] In the equation, is the sequence code ratio of the th pressure value within its adjacent numerical cluster, is the sequence code ratio of the flow rate value corresponding to the sampling time point of the th pressure value within its adjacent numerical cluster, is a defined additional value, can be , which is to avoid from being zero, is the difference between the sequence code ratios of the pressure value and the flow rate value within their respective adjacent numerical clusters, which can estimate the irrelevance between these two types of parameters. The higher the difference between the sequence code ratios, the lower the correlation between the pressure value and its corresponding flow rate value, and the higher the irrelevance. This pressure value is more likely to be an interference value. Introducing the additional value can avoid calculation obstacles caused by the same sequence code ratios and improve the accuracy of later estimation. The sequence code ratio is to arrange the values within the corresponding adjacent numerical cluster in ascending order, and then assign corresponding sequence codes to each value in the adjacent numerical cluster starting from 1 in ascending order. The sequence code divided by the number of values in the adjacent numerical cluster is the sequence code ratio of the corresponding value.

[0031] Then, the adjacent numerical cluster of the flow rate value is segmented in the same way as the pressure value, and then the initial interference amount of each flow rate value is estimated in the same way.

[0032] In a preferred but non-limiting embodiment of the present invention, step 2-1-1 further includes: estimating the initial interference amount of the flow rate value through the error between each flow rate value and the mean of its adjacent numerical cluster, plus the irrelevance between the flow rate value and the corresponding pressure value at its sampling time point. The initial interference amount of the flow rate value is proportional to both the error and the irrelevance. The calculation equation for the initial interference amount of the flow rate value is: ;

[0033] In the equation, is the starting interference amount of the th flow value, is the value of the th flow value, is the mean of the near-field numerical cluster of the th flow value, is the irrelevance between the pressure value corresponding to the sampling time point of the th flow value and the th flow value, is to perform standardization on using the Z-score standardization method. is the error between the flow value and the mean of its near-field numerical cluster.

[0034] In this equation, represents the error of the th flow value corresponding to the mean of the flow values in its near-field numerical cluster. The higher is, the more significant the th flow value is, and it is more likely to be a true high-flow outlier caused by changes in the parameters of the turbine inlet itself or an outlier caused by an interference value. To identify the specific situation more precisely, analyze the irrelevance between the th flow value and the pressure value corresponding to its sampling time point , and . If the coherence between a flow value and the pressure value corresponding to its sampling time point is lower (the irrelevance is higher), it indicates that this flow value is more likely to be an interference value rather than a true outlier caused by changes in the parameters of the turbine inlet itself.

[0035] Whether it is the starting interference amount of the pressure value or the flow value above, it all focuses on the significance of the parameter and the coherence between the two types of parameters. The higher the starting interference amount of a parameter, the more likely it is that this parameter is an interference value.

[0036] Step 2-1-2: According to the difference between the starting interference amount of each pressure value and the starting interference amount of the pressure values within its near-field numerical cluster, combined with the confidence level of the difference, calculate the final interference amount of each pressure value.

[0037] ​In the above method, the starting interference amounts of each pressure value have been calculated. However, the starting interference amount only performs analysis based on the near-domain numerical cluster in a short time period and cannot comprehensively represent the global attributes of the pressure value. If the pressure values (near-domain numerical cluster) in a short time period around a pressure value are considered, the global trend in a long time period is often ignored, resulting in inaccurate analysis values. For example, when a pressure value is compared with its near-domain numerical cluster, the error is not small and it is often an outlier interference value. However, when compared with all the pressure values, the global error with all the pressure values is not high and it cannot be determined as an interference value. Therefore, the above starting interference amount needs to be further improved to obtain a more accurate final interference amount for each pressure value.

[0038] To improve the accuracy of estimating the outlier amount of each pressure value, the calculation method of the outlier amount needs to be further improved. According to the method in step 2-1-1, calculate the starting interference amount of all the pressure values within the near-domain numerical cluster of this pressure value. Then, calculate the difference between the starting interference amount of this pressure value and the starting interference amounts of other pressure values in its near-domain numerical cluster. The higher the difference, the more significant the comparison between the outlier amount of this pressure value and the outlier amounts of other values. Also, the confidence level of this difference is analyzed. If the difference between the starting interference amount of a pressure value and the starting interference amounts of other pressure values in its near-domain numerical cluster is higher and the confidence level of this difference is greater, then this pressure value is more likely to be a true interference value.

[0039] Therefore, by combining the differences between each pressure value and the starting interference amounts of other pressure values in its near-domain numerical cluster, and combining the confidence level of this difference, calculate the final interference amount of this pressure value. In a preferred but non-limiting embodiment of the present invention, in step 2-1-2, the calculation equation for the final interference amount of the pressure value is: ; In the equation, is the final interference amount of the th pressure value, is the number of pressure values outside this pressure value itself within the near-domain numerical cluster of the th pressure value, is the serial number of the pressure value outside this pressure value itself within the near-domain numerical cluster of the th pressure value (the serial number of the pressure value outside this pressure value itself within the near-domain numerical cluster of the th pressure value is: arrange the pressure values outside this pressure value itself within the near-domain numerical cluster of the th pressure value in the order of their sampling time points, and then sequentially assign corresponding serial numbers starting from one and incrementing by one to each pressure value according to this order), is the starting interference amount of the th pressure value, is the The starting interference amount of the near-domain numerical cluster of a pressure value outside the pressure value itself at the pressure value, represents the difference between the starting interference amount of the pressure value and the starting interference amount of the pressure value within its near-domain numerical cluster. The higher this difference, the higher the difference between the outlier amount of the pressure value and the starting interference amount of other pressure values in its near domain, indicating that the pressure value is more likely to be an interfering value or an outlier, and its final interference amount is also higher. is the credibility amount of , and its operation equation is: ; In the equation, is the outlier amount of the flow value corresponding to the sampling time point of the pressure value (the outlier amount of this flow value can be the starting interference amount of this flow value), is the time interval between the pressure value and the pressure value outside this pressure value within its near-domain numerical cluster (the time interval between the pressure value and the pressure value outside this pressure value within its near-domain numerical cluster is: the absolute value of the difference between the sampling time point of the pressure value and the sampling time point of this pressure value). The higher the value of the credibility amount, the higher the referenceability of the difference between the starting interference amount of the pressure value and the starting interference amount of the pressure value within its near-domain numerical cluster.

[0040] Here, the following details are carried out on the credibility amount: (1) Regarding the time interval between the pressure value and the pressure value within its near-domain numerical cluster, the lower the time interval, the closer the pressure value and the pressure value within its near-domain numerical cluster in terms of time point, and the difference between the starting interference amounts of the pressure value and the pressure value within its near-domain numerical cluster is more credible. The higher the referability, and the lower the time interval between two pressure values, it indicates that the two pressure values are closer in time point, and the changes between the two pressure values are more coherent and closer. Therefore, if the time interval between two pressure values is not large, the difference in the initial interference amount between the two pressure values is more credible, because such differences represent real numerical differences rather than sudden fluctuations. Therefore, the referability of such differences is better.

[0041] (2)Regarding the initial interference amount of the flow value corresponding to the sampling time point of the th pressure value, the lower this initial interference amount, it indicates that the calculated value is more accurate, and the referability of is stronger. If the initial interference amount of the flow value corresponding to the sampling time point of the pressure value is not high, it indicates that the change of the flow value around this sampling time point is very stable. Therefore, the calculated is more accurate. In such a situation, it can well represent the real outlier situation, rather than the non-real outlier situation formed by the fluctuation of the flow value itself. Therefore, the referability of such differences is better.

[0042] Since the near-domain numerical clusters of each pressure value contain the pressure value itself and its adjacent pressure values before and after it, therefore, when calculating the initial interference amount of all pressure values in the near-domain numerical cluster of this pressure value, essentially based on a wider pressure value interval, when analyzing the initial interference amount of all pressure values in the near-domain numerical cluster of a pressure value, more corresponding pressure values are involved, which can make the finally calculated interference amount better represent the global attributes of the pressure value, and the calculation of the finally interference amount involves the outlier difference and its credibility amount, which is conducive to more accurately identifying the interference value and improving the credibility of the later estimation.

[0043] Step 2-2: Determine whether this pressure value needs to be included in this cleaning according to the level of the final interference amount of each pressure value. If there is a need, actively adjust the unit capacity of this pressure value using the final interference amount of each pressure value. If there is no need, this pressure value will not perform this cleaning.

[0044] Before adjusting the unit capacity of each pressure value, because the interference amounts of the approximate interference values represented by different pressure values are different, and the cleaning frequencies they need to perform are also different. In step 2-2, the cleaning frequency of each pressure value is accurately determined at the beginning.

[0045] In the preferred but non-limiting embodiment of the present invention, in step 2-2, a critical value of the final interference amount needs to be defined. Here, the critical value can be defined as , this value can be determined according to specific requirements. Before each cleaning, for each pressure value, compare it with the last interference amount of this pressure value and the critical value of the last interference amount according to this rule. If the last interference amount of this pressure value is higher than the critical value of the last interference amount, it indicates that the interference value of this pressure value has not been completely removed, the cleaning function is not good, and there is still a need to add this cleaning to improve the accuracy of the pressure value. On the contrary, if the last interference amount of this pressure value is not higher than the critical value of the last interference amount, it indicates that the cleaning function of this pressure value has achieved the goal and can represent the true situation of this pressure value. At this time, this pressure value does not need to be added to this cleaning to avoid over-treatment.

[0046] Before each cleaning, for each pressure value, according to this compliance confirmation rule, obtain its most appropriate cleaning frequency to ensure that when efficiently removing the interference value, the original attributes of the pressure value are maintained as much as possible.

[0047] Then, before each cleaning, for all pressure values that need to be added to this cleaning, actively adjust their unit capacity according to the level of their last interference amount to achieve a flexible and active cleaning process.

[0048] When a pressure value undergoes the first cleaning, if its last interference amount is not low, it also indicates that the probability of being interfered by the interference value is not small. In this case, to more efficiently remove the interference value, a relatively large unit capacity is required because a relatively large unit capacity can introduce a larger amount of pressure values for cleaning, thereby enhancing the function of eliminating the interference value and improving the accuracy and reliability of the pressure value. When a pressure value is in several cleaning cycles, according to the rule that the higher the last interference amount, the higher the required unit capacity. In addition, the cleaning function of this pressure value after the previous cleaning also needs to be considered. If the cleaning function of this pressure value after the previous cleaning is not good, for example, after the previous cleaning, the last interference amount of this pressure value is still not low and there is no significant change compared with before, it indicates that the previous cleaning did not completely remove the interference value. Then, a higher unit capacity is required for more detailed cleaning during this cleaning. Through such methods, the true effect of each cleaning can be accurately estimated, providing a basis for adjusting the level of the unit capacity of the next cleaning.

[0049] According to the above rule analysis, actively adjust the level of the unit capacity of each pressure value at different cleaning frequencies to ensure that during the cleaning process, not only can the appropriate unit capacity be flexibly selected according to the interference value attribute of the pressure value to efficiently remove the interference value, but also the true change attribute of the pressure value can be well maintained, providing pressure value conditions with excellent accuracy for the subsequent pressure value analysis.

[0050] In a preferred but non-limiting embodiment of the present invention, in step 2-2, the operation equation for actively adjusting the unit capacity is: ; In the equation, is the Euler number, is the current cleaning frequency of the th pressure value, is the main parameter for determining which type of operation method to use. Different cleaning frequencies correspond to different operation rules, showing the flexible matching of the pressure value processing process. is the th pressure value when adding the th cleaning. The actively adjusted unit capacity is the main purpose of this operation. Whether it is accurately determined will determine the cleaning function. is the unit capacity predefined for all pressure values. Its value can be ten. It gives a starting capacity for the operation, ensuring that all values have a common starting unit capacity without the influence of other factors. is the th pressure value. The last interference amount before adding the th cleaning. This represents the interference value situation of the current pressure value when the cleaning is about to be performed. It is the main basis for determining the adjustment of the unit capacity level. is the th pressure value. The last interference amount before adding the th cleaning. By comparing it with the current interference amount, the cleaning function of the previous cleaning can be estimated. and have the same operation method, only the cleaning frequency is different. has its corresponding numerical meaning in the equation and is used as the weight for adjusting the unit capacity level according to the change of the interference amount.

[0051] In this equation, represents the cleaning function of the th pressure value when adding the th cleaning. If the cleaning function of the th cleaning is better, should be higher. When performing the th cleaning, the unit capacity should be reduced. If the cleaning function of the th cleaning is not good, then should be lower. When performing the th cleaning, the unit capacity should be increased. Therefore, through construct the inverse relationship between the cleaning function of the th cleaning and the level of the unit capacity during the th cleaning.

[0052] In this equation, mainly involves the th pressure value at the The cleaning function of the second cleaning and the last interference amount of the pressure value at the current time, that is: When it is not small, it indicates that the last interference amount of the pressure value before adding the second cleaning is not low, that is, the interference amplitude of the current pressure value by the interference value is very large. A certain unit capacity is required to effectively eliminate the interference value and reduce the effect of the interference value.

[0053] And if is not high, then is not small, indicating that the cleaning function of the previous cleaning is not good, and the last interference amount of the pressure value after the previous cleaning has not been significantly reduced. At this time, is not small, which indicates that during the second cleaning, the unit capacity needs to be increased to more efficiently remove the interference value; and if is not small, then is not high, indicating that the cleaning function of the previous cleaning is good, and the last interference amount of the pressure value after the previous cleaning has been significantly reduced. At this time, even if is not small, but will not be too high, which means that in the case of good cleaning function in the previous cleaning, there is no need to increase the unit capacity too much to avoid removing too many values that represent the true changing attributes of the pressure value.

[0054] When is not high, it indicates that the last interference amount of the pressure value before adding the second cleaning is not high. In terms of maintaining the correctness of the pressure value, it is expected that the unit capacity is not too high to avoid excessive cleaning of the pressure value. And if is not high, then is not small, indicating that the cleaning function of the previous cleaning is not good, and the last interference amount of the pressure value after the previous cleaning has not been significantly reduced. However, because is not high, at this time will not be too high, which means that in the case of good cleaning function in the previous cleaning, there is no need to increase the unit capacity too much to avoid removing too many values that represent the true changing attributes of the pressure value; and if is not small, then is not high, indicating that the cleaning function of the previous cleaning is good, and the last interference amount of the pressure value after the previous cleaning has been significantly reduced, and is not high. At this time The value will be lower, which better indicates that under the condition of low current interference and good previous cleaning function, a relatively low unit capacity should be maintained to greatly maintain the true variation property of the pressure value.

[0055] Through the method of actively adjusting the unit capacity in step 2-2, before each cleaning, according to the variation property of the pressure value, such as the interference amount and its variation trend, the unit capacity of each pressure value can be adjusted flexibly. Such flexible adjustment enables the subsequent cleaning process to approach the true situation of the pressure value, greatly improving the cleaning accuracy. For example, in the detection of the parameter situation at the inlet of the turbine, the pressure values sampled at different times are often interfered by interference values of different amplitudes and types. Through such active unit capacity adjustment, appropriate anti-interference methods can be applied according to different pressure value situations, enabling the pressure value after cleaning to more accurately represent the true parameter situation at the inlet of the turbine. Based on such pressure values, more accurate detection of the parameter situation at the inlet of the turbine can be carried out, which is conducive to efficiently monitoring the inlet of the turbine.

[0056] Step 3: Perform cleaning according to the unit capacity adjusted for all pressure values, and perform condition detection on the inlet of the turbine according to the overall pressure value after cleaning.

[0057] Use the Gaussian filtering method to perform cleaning on all pressure values added to each cleaning. During each cleaning, perform cleaning and anti-interference according to the unit capacity adjusted for each pressure value. This can not only efficiently remove interference values but also maintain the trend and main properties of the pressure values, and finally obtain the pressure values after anti-interference. The pressure values after anti-interference contain the true values of the variation of the parameters of the inlet of the turbine itself. Based on this true value, accurate detection of the parameter situation at the inlet of the turbine can be carried out in the later stage.

[0058] In a preferred but non-limiting embodiment of the present invention, in step 3, the method for performing condition detection on the inlet of the turbine according to the overall pressure value after cleaning includes: Taking the overall pressure value after cleaning as the past period pressure value, then taking the past period pressure value as the dependent variable and the sampling time point corresponding to the past period pressure value as the independent variable, and using the least squares method to obtain the regression equation of the past period pressure value. According to the regression equation, obtain the estimated values of the pressure values at each sampling time point after the defined period; define the next sampling time point after the defined period as , The next sampling time point of .

[0059] According to the difference between the estimated values and the measured values (the measured values are the pressure values sampled by the pressure sensor) of the pressure values at each time point after the defined period, determine whether the situation at the inlet of the turbine is abnormal, and thus complete the condition detection of the inlet of the turbine. The method is as follows: Calculate the difference between the measured value and the estimated value of the pressure at each time point after the defined period (this difference is the absolute value of the quantity obtained by subtracting the estimated value from the measured value); if the difference between the measured value and the estimated value of the pressure at a certain time point is higher than the predefined difference critical value, the pressure at this time point is abnormal; if the difference is not higher than the difference critical value, the pressure at this time point is not abnormal; the difference critical value between the measured value and the estimated value of the pressure defined in this application can be If the difference between the estimated value and the measured value of the pressure at the time point is higher than it is determined that the pressure at the time point is abnormal. Conversely, the pressure at the time point is not abnormal. For the next sampling time point facing is Regarding the time point and all the pressure values sampled previously at this time point as the new past pressure values, repeat steps 1 to 3 again. Take the obtained interference-removed past pressure values as the dependent variable and the corresponding sampling time points of the past pressure values as the independent variable, use the least squares method to obtain the regression equation of the past pressure values, and obtain the estimated value of the pressure at the time point according to the regression equation, and then confirm whether the pressure at the time point is abnormal; for the next sampling time point facing the time point Regarding the time point and all the pressure values sampled previously at this time point as the new past pressure values, repeat steps 1 to 3 again. Take the obtained interference-removed past pressure values as the dependent variable and the corresponding sampling time points of the past pressure values as the independent variable, use the least squares method to obtain the regression equation of the past pressure values, and obtain the estimated value of the pressure at the time point according to the regression equation, and then confirm whether the pressure at the time point is abnormal, and continue to loop accordingly to achieve efficient detection of whether the pressure value is abnormal.

[0060] Total the abnormal conditions of the pressure values at each time point after the defined period. After the defined period, starting from the time point when the first abnormal pressure value appears, total the number of abnormal pressure values within the defined time; if the ratio of the number of abnormal pressure values within the total number of all pressure values within the defined time is higher than the defined ratio critical value, the condition at the inlet of the turbine is in an abnormal state. If the ratio is not higher than the defined ratio critical value, the condition at the inlet of the turbine is in a normal state. The defined time in this application can be 60s, and the defined ratio critical value can be 20%. The defined time and the defined ratio critical value can be set according to specific requirements. For example, after Starting from this point in time, continuously sample all the pressure values within 60 s. During these 60 s, if the number of abnormal pressure values (the absolute value of the quantity obtained by subtracting the estimated value from the measured value) is higher than the pressure value), and the number of such values exceeds 20% of the total number of all pressure values within these 60 s, it is determined that the condition at the inlet of the turbine is in an abnormal condition, and there is a risk of failure at the inlet of the turbine. It is advisable to use corresponding methods. For example, check the inlet of the turbine and deal with the corresponding risks. If the magnetic mesh filter at the inlet of the turbine is blocked, perform cleaning and other treatment methods. If the number of abnormal pressure values within these 60 s does not exceed 20% of the total number of all pressure values within these 60 s, it is determined that the condition at the inlet of the turbine is in a normal condition. The normal condition means that the inlet of the turbine is operating correctly. When the condition at the inlet of the turbine is in an abnormal condition, the abnormal condition message can also be displayed on the display screen of the industrial control computer to inform the staff to deal with the corresponding risks.

[0061] For any point in time after the defined period, it is possible to continuously sample the pressure values for 60 s starting from this point in time, and estimate whether the condition at the inlet of the turbine is in an abnormal condition according to the same method, so as to achieve efficient detection of the parameter condition at the inlet of the turbine.

[0062] Therefore, by performing steps 1 to 3 in this application, the final interference amount of each pressure value can be estimated more accurately, improving the accuracy of abnormal detection. The calculation of the final interference amount involves outlier discrimination and its confidence level, which is conducive to more accurately identifying interference values and improving the credibility of estimation. By performing efficient estimation on the pressure values after removing interference and performing abnormal condition estimation based on the difference between the estimated value and the measured value, the parameter changes at the inlet of the turbine can be detected efficiently, and thus the risks that may occur at the inlet of the turbine can be detected and handled efficiently.

[0063] As Figure 2 shown, a parameter processing device for a compressed air energy storage device according to the present invention includes: A pressure sensor and a flow sensor are provided at the inlet of the turbine of the compressed air energy storage device. The pressure sensor and the flow sensor are both connected to the industrial control computer; the pressure sensor and the flow sensor are used to respectively sample the pressure value and the flow value of the gas at the inlet of the turbine of the compressed air energy storage device, and transmit the pressure value and the flow value to the industrial control computer for processing; The modules running on the industrial control computer include: A receiving module, which is used to receive the pressure value and the flow value at the inlet of the turbine during the defined period sampled and transmitted; An adjustment module, which is used to perform several cleanings on all the pressure values, and actively adjust the unit capacity of the pressure value before each cleaning; A detection module, which is used to perform cleaning according to the unit capacity adjusted according to the overall pressure value, and perform condition detection on the inlet of the turbine according to the overall pressure value after cleaning.

[0064] The beneficial effects of the present invention are that, compared with the prior art, the technical effects of the present invention include: Regarding the attributes of the pressure value and the flow value, carefully analyze the relationship between the pressure value and the flow value. Through compliant analysis operations, accurately determine the interference amount of each pressure value within the selected numerical cluster. According to this interference amount, construct a special anti-interference method for each pressure value, which includes automatically determining its cleaning frequency and the level of the unit capacity at different cleaning frequencies, achieving automatic active cleaning. Through such methods, it is possible to efficiently avoid the situation where high pressure values are wrongly removed, thereby obtaining more accurate cleaning values, improving the accuracy and correctness of the future pressure value estimation, and achieving accurate detection of the inlet parameter conditions of the turbine.

[0065] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention and not to limit them. Although the present invention has been described in detail with reference to the above embodiments, those skilled in the art should understand that: still can modify the specific implementation manners of the present invention or make equivalent replacements, and any modification or equivalent replacement without departing from the spirit and scope of the present invention shall be covered by the protection scope of the claims of the present invention.

Claims

1. A parameter processing method for a compressed air energy storage device, characterized in that: include: The pressure sensor and the flow sensor respectively sample the pressure value and flow value of the gas at the inlet of the turbine of the compressed air energy storage device, and transmit the pressure value and the flow value to the industrial computer for processing; The method of transmitting the pressure value and flow value to the industrial computer for processing includes: Step 1, receiving the pressure value and flow value at the turbine inlet in a defined period of time transmitted by sampling; Step 2, performing a plurality of cleaning operations on the entire pressure value, and actively adjusting the unit capacity of the pressure value before each cleaning operation; Step 3, perform cleaning according to the unit capacity adjusted according to the overall pressure value, and perform status detection on the turbine inlet according to the overall pressure value after cleaning.

2. The parameter processing method for a compressed air energy storage device according to claim 1, characterized in that: In step 1, the pressure value and flow value are sampled by the pressure sensor and the flow sensor respectively. The starting sampling time point and sampling speed of the pressure value and the flow value are the same, and the sampling speed is 2s / time. The sampling time of the pressure value and the flow value is defined as 30min, and the sampling time is used as the defined period.

3. The parameter processing method for a compressed air energy storage device according to claim 2, characterized in that: In step 2, the method for actively adjusting the unit capacity of the pressure value includes: Step 2-1, obtaining the final interference amount of each pressure value; Step 2-2, based on the final interference amount of each pressure value, determine whether the pressure value needs to be included in this cleaning. If so, use the final interference amount of each pressure value to actively adjust the unit capacity of the pressure value. If not, the pressure value will not be cleaned.

4. The parameter processing method for a compressed air energy storage device according to claim 3, characterized in that: Step 2-1 specifically includes: Step 2-1-1, estimating the initial interference amount of the pressure value through the error between each pressure value and the mean of its near-domain value cluster, combined with the independence between the pressure value and the flow value corresponding to the sampling time point; Step 2-1-2, based on the difference between the initial interference amount of each pressure value and the initial interference amount of the pressure values ​​in its nearby value cluster, combined with the credibility of the difference, the final interference amount of each pressure value is calculated.

5. The parameter processing method for a compressed air energy storage device according to claim 4, characterized in that: In step 2-1-1, the method for estimating the initial interference includes: firstly selecting the near-domain value cluster of the pressure value, that is, taking the pressure value as the midpoint, and selecting the near-domain value cluster before and after the sampling time point. The pressure value and the pressure value taken as the midpoint together form a near-domain value cluster of the pressure value; The calculation equation of the initial interference quantity of each pressure value is: ； In the equation, It is The initial disturbance of the pressure value, It is The value of the pressure value, It is The mean of the cluster of values ​​in the vicinity of the pressure value, It is The pressure value and The independence of the flow values ​​corresponding to the sampling time points of the pressure values, The Z-score standardization method is used to Implement standardization; In step 2-1-1, The pressure value and The independency of the flow values ​​corresponding to the sampling time points of the pressure values The operational equation is: ； In the equation, It is The ratio of the sequence codes of the pressure values ​​in their nearby value clusters, It is The ratio of the sequence codes of the flow value corresponding to the sampling time point of each pressure value in its near-domain value cluster, It is a defined additional value.

6. The parameter processing method for a compressed air energy storage device according to claim 5, characterized in that: Step 2-1-1 also includes: estimating the initial interference amount of the flow value by the error between each flow value and the mean of its near-domain value cluster, plus the independence between the flow value and the pressure value corresponding to the sampling time point, the calculation equation of the initial interference amount of the flow value is: ； In the equation, It is The initial disturbance of the flow value, It is The value of the flow value, It is The mean of the cluster of values ​​in the neighborhood of the flow value, It is The flow value and The independence of the corresponding pressure values ​​at the sampling time of each flow value, The Z-score standardization method is used to Implement standardization; In step 2-1-2, the calculation equation of the final interference value of the pressure value is: ； In the equation, It is The final disturbance of the pressure value, It is The number of pressure values ​​in the cluster of nearby values ​​of a pressure value other than the pressure value itself, It is The sequence code of the pressure value other than the pressure value itself in the cluster of local values ​​of the pressure value. It is The initial disturbance of the pressure value, It is The cluster of values ​​near the pressure value is the value outside the pressure value itself. The initial disturbance of the pressure value, Representative The initial disturbance of the pressure value and the value of the cluster of its nearby values The difference in the initial interference amount of the pressure value is yes The credible quantity, its operation equation is: ； In the equation, It is The outlier of the flow value corresponding to the sampling time point of the pressure value, It is The pressure value and the number of its neighboring value clusters outside the pressure value The time interval between pressure values.

7. The parameter processing method for a compressed air energy storage device according to claim 6, characterized in that: In step 2-2, the critical value of the final interference amount is defined. Before each cleaning, each pressure value is compared based on the final interference amount of the pressure value and the critical value of the final interference amount. If the final interference amount of the pressure value is higher than the critical value of the final interference amount, it is necessary to add it to this cleaning. Conversely, if the final interference amount of the pressure value is not higher than the critical value of the final interference amount, it is not necessary to add this pressure value to this cleaning.

8. The parameter processing method for a compressed air energy storage device according to claim 7, characterized in that: In step 2-2, the operational equation for actively adjusting the unit capacity is: ； In the equation, is the Euler number, It is The cleaning frequency at the current pressure value, It is The pressure value is added During the first cleaning, the unit capacity is automatically adjusted. It is a unit capacity predefined for the overall pressure value. It is The pressure value is added The final disturbance amount before cleaning, It is The pressure value is added The final disturbance amount before the first cleaning.

9. The parameter processing method for a compressed air energy storage device according to claim 8, characterized in that: In step 3, the method for performing condition detection on the inlet of the turbine according to the overall pressure value after cleaning includes: The overall pressure value after cleaning is regarded as the past pressure value, and then the past pressure value is regarded as the dependent variable and the sampling time point corresponding to the past pressure value is regarded as the independent variable, and the regression equation of the past pressure value is obtained by using the least square method, and the estimated value of the pressure value at each sampling time point after the defined period is obtained according to the regression equation; According to the difference between the estimated value and the measured value of the pressure value at each time point after the defined period, it is confirmed whether the condition at the inlet of the turbine is abnormal, and the condition detection at the inlet of the turbine is achieved accordingly, and the method is as follows: Calculate the difference between the measured value and the estimated value of the pressure value at each time point after the defined period; if the difference between the measured value and the estimated value of the pressure value at a time point is higher than a predefined difference threshold, the pressure value at the time point is abnormal; if the difference is not higher than the difference threshold, the pressure value at the time point is not abnormal; A total is performed to determine whether the pressure value at each time point after the defined time period is in an abnormal condition. After the defined time period, starting from the time point when the first abnormal pressure value is generated, the number of abnormal pressure values ​​in the defined time is totaled; if the ratio of the number of abnormal pressure values ​​to the total number of all pressure values ​​in the defined time is higher than a defined ratio critical value, the condition at the turbine inlet is in an abnormal condition; if the ratio is not higher than the defined ratio critical value, the condition at the turbine inlet is in a normal condition.

10. A parameter processing device for a compressed air energy storage device, characterized in that: include: A pressure sensor and a flow sensor are provided at the inlet of the turbine of the compressed air energy storage device, and both the pressure sensor and the flow sensor are connected to the industrial computer; The modules running on the industrial computer include: A receiving module, which is used to receive the pressure value and flow value at the inlet of the turbine in a defined period of time transmitted by sampling; A regulating module, which is used to perform a number of cleanings on the entire pressure value and actively adjust the unit capacity of the pressure value before each cleaning; The detection module is used to perform cleaning according to the unit capacity adjusted according to the overall pressure value, and perform status detection on the inlet of the turbine according to the overall pressure value after cleaning.

Images