CN118090078A - Leakage online monitoring method for closed circulation water cooling system - Google Patents
Leakage online monitoring method for closed circulation water cooling system Download PDFInfo
- Publication number
- CN118090078A CN118090078A CN202410472879.3A CN202410472879A CN118090078A CN 118090078 A CN118090078 A CN 118090078A CN 202410472879 A CN202410472879 A CN 202410472879A CN 118090078 A CN118090078 A CN 118090078A
- Authority
- CN
- China
- Prior art keywords
- pressure
- data point
- temperature
- scatter diagram
- pressure data
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 46
- 238000001816 cooling Methods 0.000 title claims abstract description 45
- 238000000034 method Methods 0.000 title claims abstract description 40
- 238000012544 monitoring process Methods 0.000 title claims abstract description 25
- 238000010586 diagram Methods 0.000 claims abstract description 139
- YHXISWVBGDMDLQ-UHFFFAOYSA-N moclobemide Chemical compound C1=CC(Cl)=CC=C1C(=O)NCCN1CCOCC1 YHXISWVBGDMDLQ-UHFFFAOYSA-N 0.000 claims abstract 4
- 101100134058 Caenorhabditis elegans nth-1 gene Proteins 0.000 claims description 13
- 238000010606 normalization Methods 0.000 claims description 10
- 230000008569 process Effects 0.000 claims description 4
- 230000008859 change Effects 0.000 abstract description 11
- 238000012545 processing Methods 0.000 abstract description 4
- 230000000875 corresponding effect Effects 0.000 description 9
- 230000000694 effects Effects 0.000 description 3
- 230000002159 abnormal effect Effects 0.000 description 2
- 230000005856 abnormality Effects 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000000498 cooling water Substances 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000035772 mutation Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000003313 weakening effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M3/00—Investigating fluid-tightness of structures
- G01M3/02—Investigating fluid-tightness of structures by using fluid or vacuum
- G01M3/26—Investigating fluid-tightness of structures by using fluid or vacuum by measuring rate of loss or gain of fluid, e.g. by pressure-responsive devices, by flow detectors
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Examining Or Testing Airtightness (AREA)
Abstract
The invention relates to the technical field of data processing, in particular to an online leakage monitoring method for a closed circulation water cooling system, which comprises the following steps: the method comprises the steps of obtaining a pressure data sequence and a temperature data sequence of a closed circulation water cooling system, respectively constructing a pressure scatter diagram and a temperature scatter diagram, obtaining the primary noise degree of each pressure data point according to the pressure value of each pressure data point in the pressure scatter diagram and the pressure value of the surrounding pressure data points, combining the correlation of each pressure data point and the temperature data point to obtain the final noise degree of each pressure data point in the pressure scatter diagram, correcting the order of MA items in an ARIMA model, predicting the pressure data sequence, and judging whether leakage occurs in the closed circulation water cooling system. According to the invention, through analyzing the change characteristics of the pressure data and the temperature data, the influence of noise data points is weakened, the order of MA items is corrected, and the accuracy of on-line leakage monitoring of the closed circulation water cooling system is improved.
Description
Technical Field
The invention relates to the technical field of data processing, in particular to an online leakage monitoring method for a closed circulation water cooling system.
Background
The closed circulation water cooling system leakage on-line monitoring is an important technical application, and can effectively ensure the normal operation of the water cooling system and the safety and stability of equipment. Once the closed circulation water cooling system leaks, equipment damage or work abnormality can be caused, the system leakage problem can be found in time through online monitoring and prediction, and loss and maintenance cost are reduced. The existing method for predicting the pressure of the closed circulation water cooling system is an ARIMA model, has the advantages of strong interpretability, wide application range, good prediction capability, and capability of effectively capturing the characteristics of time sequence data and providing accurate prediction results.
The prior art has the problem that various external interference factors such as temperature change, mechanical vibration, electromagnetic interference and the like can exist in the system operation environment in the process of collecting the pressure of the closed circulation water cooling system, and the factors can influence the measurement of the pressure sensor, so that the collected data is interfered to generate noise. The existing ARIMA model obtains the order of MA items by analyzing ACF and PACF graphs of time series data to determine the order of MA items. And the noise data of the pressure may show similar characteristics with the real abnormal data of the pressure in the graph, thereby reducing the accuracy of on-line leakage monitoring of the closed circulation water cooling system.
Disclosure of Invention
The invention provides an online leakage monitoring method for a closed circulation water cooling system, which aims to solve the existing problems.
The invention discloses an on-line leakage monitoring method for a closed circulation water cooling system, which adopts the following technical scheme:
The embodiment of the invention provides an on-line leakage monitoring method for a closed circulation water cooling system, which comprises the following steps:
Acquiring a pressure data sequence and a temperature data sequence of a closed circulation water cooling system, and respectively constructing a pressure scatter diagram and a temperature scatter diagram according to the pressure data sequence and the temperature data sequence;
Obtaining the primary noise degree of each pressure data point in the pressure scatter diagram according to the pressure value of each pressure data point in the pressure scatter diagram and the pressure value of the surrounding pressure data point of each pressure data point; obtaining the primary noise degree of each temperature data point in the temperature scatter diagram according to the acquisition mode of the primary noise degree of each pressure data point in the pressure scatter diagram;
According to the primary noise degree of each pressure data point in the pressure scatter diagram and the primary noise degree of each temperature data point in the temperature scatter diagram, obtaining the correlation between each pressure data point in the pressure scatter diagram and each temperature data point in the temperature scatter diagram;
Obtaining the final noise degree of each pressure data point in the pressure scatter diagram according to the correlation between each pressure data point in the pressure scatter diagram and each temperature data point in the temperature scatter diagram and the primary noise degree of each pressure data point in the pressure scatter diagram;
correcting the order of the MA item in the ARIMA model according to the final noise degree of each pressure data point in the pressure scatter diagram to obtain the corrected MA item order;
According to the modified MA item order, predicting the pressure data sequence by using an ARIMA model to obtain a pressure data predicted value; and judging whether the closed circulation water cooling system leaks or not according to the predicted value of the pressure data.
Further, the method respectively constructs a pressure scatter diagram and a temperature scatter diagram according to the pressure data sequence and the temperature data sequence, and comprises the following specific steps:
Each pressure data in the pressure data sequence corresponds to one time, the time is taken as a horizontal axis, the pressure data is taken as a vertical axis, and a pressure scatter diagram of the pressure data sequence is constructed; and each temperature data in the temperature data sequence corresponds to one time, the time is taken as a horizontal axis, and the temperature data is taken as a vertical axis, so that a temperature scatter diagram of the temperature data sequence is constructed.
Further, the step of obtaining the primary noise level of each pressure data point in the pressure scatter diagram according to the pressure value of each pressure data point in the pressure scatter diagram and the pressure value of the surrounding pressure data point of each pressure data point comprises the following specific steps:
The method comprises the steps of obtaining a reference data segment of each pressure data point in a pressure scatter diagram, and obtaining a specific calculation formula corresponding to the primary noise degree of each pressure data point according to the pressure value of each pressure data point in the pressure scatter diagram and the pressure value of each pressure data point in the reference data segment of each pressure data point, wherein the specific calculation formula comprises the following steps:
In the method, in the process of the invention, Representing the primary noise level of the nth pressure data point in the pressure scatter plot,/>Representing the pressure value of the nth pressure data point in the pressure scatter plot,/>Mean value of pressure values representing all pressure data points in the reference data segment of the nth pressure data point in the pressure scatter plot,/>Cosine value representing minimum angle between the line connecting the nth and (n-1) th pressure data points and the line connecting the nth and (n+1) th pressure data points in pressure scatter diagram,/>Representing the number of pressure data points in the reference data segment of the nth pressure data point in the pressure scatter plot,/>Absolute value representing the difference in pressure values of the nth and nth-1 pressure data points in the pressure scatter plot,/>Absolute value representing the difference in pressure values of the u-th and u-1-th pressure data points in the reference data segment of the nth pressure data point in the pressure scatter plot,/>As an exponential function based on natural constants,/>As a linear normalization function,/>As a function of absolute value.
Further, the acquiring the reference data segment of each pressure data point in the pressure scatter diagram includes the following specific steps:
In the pressure scatter diagram, a sequence of any one pressure data point and K nearest pressure data points is recorded as a reference data segment of the any one pressure data point, Is a preset quantity threshold.
Further, the step of obtaining the correlation between each pressure data point in the pressure scatter diagram and each temperature data point in the temperature scatter diagram according to the primary noise degree of each pressure data point in the pressure scatter diagram and the primary noise degree of each temperature data point in the temperature scatter diagram comprises the following specific steps:
acquiring a reference data segment of each temperature data point in the temperature scatter diagram;
According to the pressure scatter diagram and the difference of the primary noise degree of the pressure data points and the temperature data points in the temperature scatter diagram, obtaining a first characteristic of each pressure data point in the pressure scatter diagram;
And obtaining the correlation between each pressure data point in the pressure scatter diagram and each temperature data point in the temperature scatter diagram according to the first characteristic of each pressure data point in the pressure scatter diagram and the difference between the pressure data point and the reference data segment of the temperature data point in the pressure scatter diagram and the temperature scatter diagram.
Further, the acquiring the reference data segment of each temperature data point in the temperature scatter diagram includes the following specific steps:
and obtaining the reference data segment of each temperature data point in the temperature scatter diagram according to the acquisition mode of the reference data segment of each pressure data point in the pressure scatter diagram.
Further, according to the first feature of each pressure data point in the pressure scatter diagram and the differences between the pressure data point and the reference data segment of the temperature data point in the pressure scatter diagram and the temperature scatter diagram, the correlation between each pressure data point in the pressure scatter diagram and each temperature data point in the temperature scatter diagram is obtained, and the corresponding specific calculation formula is as follows:
In the middle of Representing the correlation of the nth pressure data point in the pressure scatter plot with the nth temperature data point in the temperature scatter plot,/>Representing the primary noise level of the nth pressure data point in the pressure scatter plot,/>Representing a primary noise level of an nth temperature data point in the temperature scatter plot; /(I)Representing the slope of the nth and nth-1 pressure data point line in the pressure scattergram,/>Representing the slope of the nth and nth-1 temperature data point links in the temperature scattergram; /(I)Representing a number of pressure data points in a reference data segment of an nth pressure data point in the pressure scatter plot; /(I)Slope of the line between the (u) th and (u-1) th pressure data points in the reference data segment representing the (n) th pressure data point in the pressure scatter plot,/>A slope of a line connecting the (u) th and (u-1) th temperature data points in a reference data segment representing the (n) th temperature data point of the temperature scattergram; /(I)Representing the percentile of the nth pressure data point in its reference data segment in a pressure scatter plot,/>Representing the percentile of the nth temperature data point in its reference data segment in a temperature scatter plot,/>For a preset first constant,/>As an exponential function based on natural constants,/>As an absolute value function,/>Is the first feature of the nth pressure data point in the pressure scatter plot.
Further, the final noise degree of each pressure data point in the pressure scatter diagram is obtained according to the correlation between each pressure data point in the pressure scatter diagram and each temperature data point in the temperature scatter diagram and the primary noise degree of each pressure data point in the pressure scatter diagram, and the corresponding specific calculation formula is as follows:
In the middle of Representing the final noise level of the nth pressure data point in the pressure scatter plot,/>Representing the primary noise level of the nth pressure data point in the pressure scatter plot,/>Representing the correlation of the nth pressure data point in the pressure scatter plot with the nth temperature data point in the temperature scatter plot,/>As a linear normalization function,/>Is a preset second constant.
Further, the steps of the MA term in the ARIMA model are corrected according to the final noise degree of each pressure data point in the pressure scatter diagram, and the corrected steps of the MA term are obtained, and the corresponding specific calculation formula is as follows:
In the middle of Representing the modified MA item order,/>Representing the number of pressure data points in the pressure scatter plot,/>Representing the final noise level of the nth pressure data point in the pressure scatter plot,/>Representing a preset MA item order,/>Representing a rounding down, c is a preset third constant.
Further, the method for judging whether the closed circulation water cooling system leaks according to the predicted value of the pressure data comprises the following specific steps:
and when the predicted value of the pressure data is smaller than a preset judging threshold value, judging that the closed circulation water cooling system leaks.
The technical scheme of the invention has the beneficial effects that:
The method comprises the steps of obtaining a pressure data sequence and a temperature data sequence of a closed circulation water cooling system, respectively constructing a pressure scatter diagram and a temperature scatter diagram, accurately reflecting the change trend of data in an image, obtaining the primary noise degree of each pressure data point in the pressure scatter diagram according to the pressure value of each pressure data point in the pressure scatter diagram and the pressure value of surrounding pressure data points, combining the correlation of each pressure data point and the temperature data point, obtaining the final noise degree of each pressure data point in the pressure scatter diagram, correcting the order of MA items in an ARIMA model, weakening the influence of noise data, better describing the short-term fluctuation and random fluctuation of the pressure data, improving the prediction precision of the model, predicting the pressure data sequence, judging whether the closed circulation water cooling system leaks, providing more accurate prediction and analysis, enabling the system to be more convenient to operate and regulate, and improving the efficiency and stability of the system.
Drawings
In order to more clearly illustrate the embodiments of the invention or the technical solutions in the prior art, the drawings that are required in the embodiments or the description of the prior art will be briefly described, it being obvious that the drawings in the following description are only some embodiments of the invention, and that other drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.
FIG. 1 is a flow chart of the steps of a method for monitoring leakage of a closed circulation water cooling system on line;
FIG. 2 is a graph of pressure data in the present embodiment;
fig. 3 is a graph of temperature data in the present embodiment.
Detailed Description
In order to further describe the technical means and effects adopted by the invention to achieve the preset aim, the following is a detailed description of specific embodiments, structures, features and effects of a closed circulation water cooling system leakage on-line monitoring method according to the invention with reference to the accompanying drawings and preferred embodiments. In the following description, different "one embodiment" or "another embodiment" means that the embodiments are not necessarily the same. Furthermore, the particular features, structures, or characteristics of one or more embodiments may be combined in any suitable manner.
Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs.
The invention provides a specific scheme of an on-line leakage monitoring method for a closed circulation water cooling system, which is specifically described below with reference to the accompanying drawings.
Referring to fig. 1, a flow chart of steps of a leakage online monitoring method for a closed circulation water cooling system according to an embodiment of the invention is shown, and the method includes the following steps:
step S001: and acquiring a pressure data sequence and a temperature data sequence of the closed circulation water cooling system, and respectively constructing a pressure scatter diagram and a temperature scatter diagram according to the pressure data sequence and the temperature data sequence.
When the pressure data and the temperature data of the closed circulation water cooling system are collected, the closed circulation water cooling system needs to be ensured to be carried out in the same time period, and specifically, the pressure data at each moment can be collected by using a pressure sensor to obtain a pressure data sequence, and the temperature data at each moment can be collected by using a temperature sensor to obtain a temperature data sequence. The specific acquisition duration is thirty minutes, and the acquisition frequency is one second, which is described as an example, but other values may be set in other embodiments, and the embodiment is not limited thereto. Constructing a pressure scatter diagram of a pressure data sequence by taking time as a horizontal axis and pressure data as a vertical axis; and constructing a temperature scatter diagram of the temperature data sequence by taking time as a horizontal axis and temperature data as a vertical axis. Fig. 2 is a graph of pressure data, and fig. 3 is a graph of temperature data. What needs to be described is: in fig. 2, the horizontal axis represents time, the vertical axis represents pressure data, and the unit of pressure data is bar; in fig. 3, the horizontal axis represents time, the vertical axis represents temperature data, and the unit of temperature data is degrees celsius.
Step S002: obtaining the primary noise degree of each pressure data point in the pressure scatter diagram according to the pressure value of each pressure data point in the pressure scatter diagram and the pressure value of the surrounding pressure data point of each pressure data point; and obtaining the primary noise degree of each temperature data point in the temperature scatter diagram according to the acquisition mode of the primary noise degree of each pressure data point in the pressure scatter diagram.
The numerical representation of each pressure data point in the graph and the change characteristics of the surrounding data are analyzed to obtain the primary noise degree. Wherein the more abnormal the numerical performance of each data point, the stronger the abrupt change that occurs at each pressure data point, the greater the primary noise level of that pressure data point.
In the pressure scatter diagram, a sequence of any one pressure data point and K nearest pressure data points is recorded as a reference data segment of the any one pressure data point, thereby obtaining a reference data segment of each pressure data point in the pressure scatter diagram,For the preset number threshold, in this embodiment/>50, Which is described as an example, other values may be set in other embodiments, and the present example is not limited thereto. According to the logical AND operation described above, the primary noise level for each pressure data point is calculated:
In the method, in the process of the invention, Representing the primary noise level of the nth pressure data point in the pressure scatter plot,/>Representing the pressure value of the nth pressure data point in the pressure scatter plot,/>Mean value of pressure values representing all pressure data points in the reference data segment of the nth pressure data point in the pressure scatter plot,/>Cosine value representing minimum angle between the line connecting the nth and (n-1) th pressure data points and the line connecting the nth and (n+1) th pressure data points in pressure scatter diagram,/>Representing the number of pressure data points in the reference data segment of the nth pressure data point in the pressure scatter plot. /(I)Absolute value representing the difference in pressure values of the nth and nth-1 pressure data points in the pressure scatter plot,/>Absolute value representing the difference in pressure values of the u-th and u-1-th pressure data points in the reference data segment of the nth pressure data point in the pressure scatter plot,/>The present embodiment uses/>, as an exponential function based on natural constantsThe inverse proportion relation and normalization processing are presented, and an implementer can set an inverse proportion function and a normalization function according to actual conditions; /(I)Normalizing the data values to/>, as a linear normalization functionWithin the interval,/>As a function of absolute value.
In the formulaThe relative difference of the nth pressure data point in the pressure scatter plot with respect to the average pressure value in the reference data segment is shown, and the greater the difference value is, the greater the possibility of abrupt change in the pressure data point is, the greater the primary noise level is. /(I)The larger the value is, the sharper the included angle is, the larger the difference between the values of the pressure data points and the values of the data points which are adjacent to the left and right sides is, and the larger the probability of mutation on the corresponding pressure data points is, and the larger the primary noise degree is. /(I)The average value representing the absolute value of the difference between the pressure values of the nth and the n-1 th pressure data points in the pressure scatter plot and the absolute value of the difference between the nth and the u-1 th pressure data points in the reference data segment of the nth pressure data point, the greater this value is indicative of the greater the likelihood of a sudden change occurring in the nth pressure data point, the greater its primary noise level.
Obtaining the reference data segment of each temperature data point in the temperature scatter diagram according to the acquisition mode of the reference data segment of each pressure data point in the pressure scatter diagram; and obtaining the primary noise degree of each temperature data point in the temperature scatter diagram according to the acquisition mode of the primary noise degree of each pressure data point in the pressure scatter diagram.
Step S003: and obtaining the correlation between each pressure data point in the pressure scatter diagram and each temperature data point in the temperature scatter diagram according to the primary noise degree of each pressure data point in the pressure scatter diagram and the primary noise degree of each temperature data point in the temperature scatter diagram.
By analyzing the pressure value representation of each pressure data point and its surrounding data change characteristics through the above steps, the primary noise level of each pressure data point is obtained. However, analysis of only the pressure value representation of each pressure data point and the varying characteristics of its surrounding data points does not accurately distinguish between noise data and anomaly data for the pressure. In a closed-cycle water system, when there is a leak in the system, the pressure will drop, while the temperature in the system will rise due to the reduced flow of cooling water. So that the pressure data and the temperature data have a certain negative correlation in the system. Therefore, we can analyze the correlation between the pressure data and the temperature data, and correct the primary noise level of each pressure data point in the pressure scatter plot to obtain the final noise level of each pressure data point.
Calculating a correlation between each pressure data point and a temperature data point:
In the middle of Representing the correlation of the nth pressure data point in the pressure scatter plot with the nth temperature data point in the temperature scatter plot,/>Representing the primary noise level of the nth pressure data point in the pressure scatter plot,/>Representing the primary noise level of the nth temperature data point in the temperature scatter plot. /(I)Representing the slope of the nth and nth-1 pressure data point line in the pressure scattergram,/>The slope of the nth and nth-1 temperature data point lines in the temperature scattergram is shown. /(I)Representing the number of pressure data points in the reference data segment of the nth pressure data point in the pressure scatter plot. /(I)Slope of the line between the (u) th and (u-1) th pressure data points in the reference data segment representing the (n) th pressure data point in the pressure scatter plot,/>The slope of the line connecting the (u) th and (u-1) th temperature data points in the reference data segment of the nth temperature data point of the temperature scattergram is represented. /(I)Representing the percentile of the nth pressure data point in its reference data segment in a pressure scatter plot,/>Representing the percentile of the nth temperature data point in its reference data segment in a temperature scatter plot,/>For a preset first constant, the denominator is avoided to be 0, in this embodiment/>This is described as an example, but other values may be set in other embodiments, and the present example is not limited thereto. What needs to be described is: the percentile is a well known technique, e.g. the percentile of 3 in the sequences 1,2,3,4,5 is 50%. The data points in the reference data segment are ordered by time,/>The present embodiment uses/>, as an exponential function based on natural constantsThe inverse proportion relation and normalization processing are presented, and an implementer can set an inverse proportion function and a normalization function according to actual conditions; /(I)As an absolute value function,/>Is the first feature of the nth pressure data point in the pressure scatter plot.
In the formulaThe difference value between the primary noise degree of the nth pressure data point in the pressure scatter diagram and the nth temperature data point in the temperature scatter diagram is shown, and the smaller the difference value is, the smaller the numerical performance of the two data at the corresponding moment of the nth pressure data point and the smaller the difference of the data change characteristics is, and the stronger the correlation of the two data is. /(I)The ratio of the slope of the nth and nth-1 pressure data point lines in the pressure scatter diagram to the slope of the nth and nth-1 temperature data point lines in the temperature scatter diagram is shown, and the relative magnitude relation of the slope and the slope is shown. /(I)An average value representing a ratio of a slope of the nth and nth-1 pressure data point lines in the pressure scattergram to a slope of the nth and nth-1 temperature data point lines in the reference data segment of the nth temperature data point of the temperature scattergram. /(I)The smaller the value, the more strongly the pressure data is correlated with the temperature data change at that time. /(I)The larger the difference between the relative magnitudes of the pressure data and the temperature data in the respective reference data segments at the time corresponding to the nth pressure data point, the stronger the negative correlation therebetween.
Step S004: and obtaining the final noise degree of each pressure data point in the pressure scatter diagram according to the correlation of each pressure data point in the pressure scatter diagram and each temperature data point in the temperature scatter diagram and the primary noise degree of each pressure data point in the pressure scatter diagram.
After the correlation between the pressure data and the temperature data at the corresponding time of each pressure data point is obtained, the primary noise degree of each pressure data point is combined for analysis, and the final noise degree of each pressure data point is obtained. The concrete calculation mode is as follows:
In the middle of Representing the final noise level of the nth pressure data point in the pressure scatter plot. /(I)Representing the primary noise level of the nth pressure data point in the pressure scatter plot, the greater the value will correspond to the greater the final noise level. /(I)The correlation of the nth pressure data point in the pressure scatter plot and the nth temperature data point in the temperature scatter plot is shown, and the larger the value is, the more the change characteristic between the pressure data and the temperature data accords with the characteristics in the system, and the smaller the corresponding final noise degree is. End use/>Normalize the final noise level for each pressure data point to/>Between,/>For a preset second constant, in this embodiment/>2, Which is described by way of example, other values may be set in other embodiments, and the example is not limited thereto,/>Normalizing the data values to/>, as a linear normalization functionWithin the interval.
Step S005: and correcting the order of the MA item in the ARIMA model according to the final noise degree of each pressure data point in the pressure scatter diagram to obtain the corrected MA item order.
By the analysis of the above steps, the final noise level of each pressure data point in the pressure scatter plot is obtained. The final noise level of each pressure data point in the pressure data is then combined with the MA term order obtained by the prior art to obtain an accurate MA term order, and this is taken as an example to describe an initial value 2 of the MA term order preset in this embodiment, and other values may be set in other embodiments, which is not limited in this embodiment. In the ARIMA model, more noise data means more random and unpredictable fluctuations in the time series, which can have an impact on the observations, resulting in more complex data and more fluctuations. To better capture the effects of these noise data, a higher order MA model needs to be used to improve the adaptation of the model. According to the logic, the order of the corrected MA item is obtained:
In the middle of Representing the modified MA term order. /(I)Representing the number of pressure data points in the pressure scatter plot,/>Representing the final noise level of the nth pressure data point in the pressure scatter plot,/>Representing the average of the final noise levels for all pressure data points in the pressure scatter plot. /(I)Representing a preset MA item order,/>The third constant is a predetermined value, and c is 0.5 in this embodiment, which is described as an example, but other values may be set in other embodiments, and this embodiment is not limited thereto.
What needs to be described is: the ARIMA (Autoregressive Integrated Moving Average) model is a commonly used time series prediction model for analyzing and predicting the trend and periodicity of time series data. The ARIMA model consists of three parts: autoregressive (AR), differential (I), and Moving Average (MA), which are well known techniques, and specific methods are not described herein. The MA term order is a parameter in the ARIMA model.
Step S006: according to the modified MA item order, predicting the pressure data sequence by using an ARIMA model to obtain a pressure data predicted value; and judging whether the closed circulation water cooling system leaks or not according to the predicted value of the pressure data.
According to the modified MA term order, the ARIMA model is used to predict the pressure data sequence to obtain a pressure data predicted value, and the preset judgment threshold value in this embodiment is 1, which is described by way of example, other values may be set in other embodiments, and this embodiment is not limited, and when the pressure data predicted value is smaller than the preset judgment threshold value, it is determined that leakage occurs in the closed circulation water cooling system.
The present invention has been completed.
The above description is only of the preferred embodiments of the present invention and is not intended to limit the invention, but any modifications, equivalent substitutions, improvements, etc. within the principles of the present invention should be included in the scope of the present invention.
Claims (10)
1. The on-line leakage monitoring method for the closed circulation water cooling system is characterized by comprising the following steps of:
Acquiring a pressure data sequence and a temperature data sequence of a closed circulation water cooling system, and respectively constructing a pressure scatter diagram and a temperature scatter diagram according to the pressure data sequence and the temperature data sequence;
Obtaining the primary noise degree of each pressure data point in the pressure scatter diagram according to the pressure value of each pressure data point in the pressure scatter diagram and the pressure value of the surrounding pressure data point of each pressure data point; obtaining the primary noise degree of each temperature data point in the temperature scatter diagram according to the acquisition mode of the primary noise degree of each pressure data point in the pressure scatter diagram;
According to the primary noise degree of each pressure data point in the pressure scatter diagram and the primary noise degree of each temperature data point in the temperature scatter diagram, obtaining the correlation between each pressure data point in the pressure scatter diagram and each temperature data point in the temperature scatter diagram;
Obtaining the final noise degree of each pressure data point in the pressure scatter diagram according to the correlation between each pressure data point in the pressure scatter diagram and each temperature data point in the temperature scatter diagram and the primary noise degree of each pressure data point in the pressure scatter diagram;
correcting the order of the MA item in the ARIMA model according to the final noise degree of each pressure data point in the pressure scatter diagram to obtain the corrected MA item order;
According to the modified MA item order, predicting the pressure data sequence by using an ARIMA model to obtain a pressure data predicted value; and judging whether the closed circulation water cooling system leaks or not according to the predicted value of the pressure data.
2. The method for on-line monitoring leakage of a closed circulation water cooling system according to claim 1, wherein the method comprises the following steps of:
Each pressure data in the pressure data sequence corresponds to one time, the time is taken as a horizontal axis, the pressure data is taken as a vertical axis, and a pressure scatter diagram of the pressure data sequence is constructed; and each temperature data in the temperature data sequence corresponds to one time, the time is taken as a horizontal axis, and the temperature data is taken as a vertical axis, so that a temperature scatter diagram of the temperature data sequence is constructed.
3. The method for on-line monitoring leakage of a closed circulation water cooling system according to claim 1, wherein the step of obtaining the primary noise level of each pressure data point in the pressure scatter diagram according to the pressure value of each pressure data point in the pressure scatter diagram and the pressure value of the surrounding pressure data point of each pressure data point comprises the following specific steps:
The method comprises the steps of obtaining a reference data segment of each pressure data point in a pressure scatter diagram, and obtaining a specific calculation formula corresponding to the primary noise degree of each pressure data point according to the pressure value of each pressure data point in the pressure scatter diagram and the pressure value of each pressure data point in the reference data segment of each pressure data point, wherein the specific calculation formula comprises the following steps:
In the method, in the process of the invention, Representing the primary noise level of the nth pressure data point in the pressure scatter plot,/>Representing the pressure value of the nth pressure data point in the pressure scatter plot,/>Mean value of pressure values representing all pressure data points in the reference data segment of the nth pressure data point in the pressure scatter plot,/>Cosine value representing minimum angle between the line connecting the nth and (n-1) th pressure data points and the line connecting the nth and (n+1) th pressure data points in pressure scatter diagram,/>Representing the number of pressure data points in the reference data segment of the nth pressure data point in the pressure scatter plot,/>Absolute value representing the difference in pressure values of the nth and nth-1 pressure data points in the pressure scatter plot,/>Absolute value representing the difference in pressure values of the u-th and u-1-th pressure data points in the reference data segment of the nth pressure data point in the pressure scatter plot,/>As an exponential function based on natural constants,/>As a linear normalization function,/>As a function of absolute value.
4. The method for on-line monitoring leakage of a closed circulation water cooling system according to claim 3, wherein the step of obtaining the reference data segment of each pressure data point in the pressure scatter diagram comprises the following specific steps:
In the pressure scatter diagram, a sequence of any one pressure data point and K nearest pressure data points is recorded as a reference data segment of the any one pressure data point, Is a preset quantity threshold.
5. The method for on-line monitoring leakage of a closed circulation water cooling system according to claim 1, wherein the step of obtaining the correlation between each pressure data point in the pressure scatter diagram and each temperature data point in the temperature scatter diagram according to the primary noise level of each pressure data point in the pressure scatter diagram and the primary noise level of each temperature data point in the temperature scatter diagram comprises the following specific steps:
acquiring a reference data segment of each temperature data point in the temperature scatter diagram;
According to the pressure scatter diagram and the difference of the primary noise degree of the pressure data points and the temperature data points in the temperature scatter diagram, obtaining a first characteristic of each pressure data point in the pressure scatter diagram;
And obtaining the correlation between each pressure data point in the pressure scatter diagram and each temperature data point in the temperature scatter diagram according to the first characteristic of each pressure data point in the pressure scatter diagram and the difference between the pressure data point and the reference data segment of the temperature data point in the pressure scatter diagram and the temperature scatter diagram.
6. The method for on-line monitoring leakage of a closed circulation water cooling system according to claim 5, wherein the step of obtaining the reference data segment of each temperature data point in the temperature scatter diagram comprises the following specific steps:
and obtaining the reference data segment of each temperature data point in the temperature scatter diagram according to the acquisition mode of the reference data segment of each pressure data point in the pressure scatter diagram.
7. The method for on-line monitoring leakage of a closed circulation water cooling system according to claim 5, wherein the correlation between each pressure data point in the pressure scatter diagram and each temperature data point in the temperature scatter diagram is obtained according to the first characteristic of each pressure data point in the pressure scatter diagram and the difference between the pressure data point in the pressure scatter diagram and the reference data segment of the temperature data point in the temperature scatter diagram, and the corresponding specific calculation formula is as follows:
In the middle of Representing the correlation of the nth pressure data point in the pressure scatter plot with the nth temperature data point in the temperature scatter plot,/>Representing the primary noise level of the nth pressure data point in the pressure scatter plot,/>Representing a primary noise level of an nth temperature data point in the temperature scatter plot; /(I)Representing the slope of the nth and nth-1 pressure data point line in the pressure scattergram,Representing the slope of the nth and nth-1 temperature data point links in the temperature scattergram; /(I)Representing a number of pressure data points in a reference data segment of an nth pressure data point in the pressure scatter plot; /(I)Slope of the line between the (u) th and (u-1) th pressure data points in the reference data segment representing the (n) th pressure data point in the pressure scatter plot,/>A slope of a line connecting the (u) th and (u-1) th temperature data points in a reference data segment representing the (n) th temperature data point of the temperature scattergram; /(I)Representing the percentile of the nth pressure data point in its reference data segment in a pressure scatter plot,/>Representing the percentile of the nth temperature data point in its reference data segment in a temperature scatter plot,/>For a preset first constant,/>As an exponential function based on natural constants,/>As an absolute value function,/>Is the first feature of the nth pressure data point in the pressure scatter plot.
8. The method for on-line monitoring leakage of a closed circulation water cooling system according to claim 1, wherein the final noise level of each pressure data point in the pressure scatter diagram is obtained according to the correlation between each pressure data point in the pressure scatter diagram and each temperature data point in the temperature scatter diagram and the primary noise level of each pressure data point in the pressure scatter diagram, and the specific corresponding calculation formula is as follows:
In the middle of Representing the final noise level of the nth pressure data point in the pressure scatter plot,/>Representing the primary noise level of the nth pressure data point in the pressure scatter plot,/>Representing the correlation of the nth pressure data point in the pressure scatter plot with the nth temperature data point in the temperature scatter plot,/>As a linear normalization function,/>Is a preset second constant.
9. The method for online monitoring leakage of a closed circulation water cooling system according to claim 1, wherein the steps of MA items in the ARIMA model are corrected according to the final noise degree of each pressure data point in the pressure scatter diagram, and the corrected steps of MA items are obtained, and the corresponding specific calculation formula is as follows:
In the middle of Representing the modified MA item order,/>Representing the number of pressure data points in the pressure scatter plot,/>Representing the final noise level of the nth pressure data point in the pressure scatter plot,/>Representing a preset MA item order,/>Representing a rounding down, c is a preset third constant.
10. The method for on-line monitoring leakage of a closed circulation water cooling system according to claim 1, wherein the step of judging whether leakage occurs in the closed circulation water cooling system according to the predicted value of the pressure data comprises the following specific steps:
and when the predicted value of the pressure data is smaller than a preset judging threshold value, judging that the closed circulation water cooling system leaks.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410472879.3A CN118090078B (en) | 2024-04-19 | 2024-04-19 | Leakage online monitoring method for closed circulation water cooling system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410472879.3A CN118090078B (en) | 2024-04-19 | 2024-04-19 | Leakage online monitoring method for closed circulation water cooling system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN118090078A true CN118090078A (en) | 2024-05-28 |
| CN118090078B CN118090078B (en) | 2024-06-25 |
Family
ID=91154893
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202410472879.3A Active CN118090078B (en) | 2024-04-19 | 2024-04-19 | Leakage online monitoring method for closed circulation water cooling system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN118090078B (en) |
Citations (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2016101690A1 (en) * | 2014-12-22 | 2016-06-30 | 国家电网公司 | Time sequence analysis-based state monitoring data cleaning method for power transmission and transformation device |
| US20190050515A1 (en) * | 2018-06-27 | 2019-02-14 | Intel Corporation | Analog functional safety with anomaly detection |
| KR101956715B1 (en) * | 2017-09-01 | 2019-03-11 | 군산대학교산학협력단 | Wind direction prediction method and apparatus for yaw control of wind turbines |
| CN210603662U (en) * | 2019-11-02 | 2020-05-22 | 山东艾克索仑电气有限公司 | Temperature measurement structure for permanent magnet electric roller bearing |
| CN115099143A (en) * | 2022-06-22 | 2022-09-23 | 青岛理工大学 | Uncertainty reasoning-based driving all-electric actuator early fault diagnosis method |
| CN116320042A (en) * | 2023-05-16 | 2023-06-23 | 陕西思极科技有限公司 | Internet of things terminal monitoring control system for edge calculation |
| JP2023113280A (en) * | 2022-02-03 | 2023-08-16 | 大成建設株式会社 | Air leakage detection method |
| CN116738353A (en) * | 2023-08-15 | 2023-09-12 | 安拓思纳米技术（苏州）有限公司 | Pharmaceutical workshop air filter element performance detection method based on data analysis |
| KR20230137021A (en) * | 2022-03-21 | 2023-10-04 | 주식회사 플로워크연구소 | detecting method and system of pipe abnormality using artificial intelligence |
| CN116881745A (en) * | 2023-09-07 | 2023-10-13 | 深圳市特安工业科技有限公司 | Pressure transmitter abnormality monitoring method based on big data |
| CN117167903A (en) * | 2023-11-03 | 2023-12-05 | 江苏中安建设集团有限公司 | Artificial intelligence-based foreign matter fault detection method for heating ventilation equipment |
| CN117725792A (en) * | 2023-12-18 | 2024-03-19 | 清华大学合肥公共安全研究院 | Building safety monitoring and early warning method based on time sequence and model correction |
| CN117870943A (en) * | 2024-01-22 | 2024-04-12 | 中国三峡建工(集团)有限公司 | Multi-sensor-based data optimization acquisition system in grouting process |
-
2024
- 2024-04-19 CN CN202410472879.3A patent/CN118090078B/en active Active
Patent Citations (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2016101690A1 (en) * | 2014-12-22 | 2016-06-30 | 国家电网公司 | Time sequence analysis-based state monitoring data cleaning method for power transmission and transformation device |
| KR101956715B1 (en) * | 2017-09-01 | 2019-03-11 | 군산대학교산학협력단 | Wind direction prediction method and apparatus for yaw control of wind turbines |
| US20190050515A1 (en) * | 2018-06-27 | 2019-02-14 | Intel Corporation | Analog functional safety with anomaly detection |
| CN210603662U (en) * | 2019-11-02 | 2020-05-22 | 山东艾克索仑电气有限公司 | Temperature measurement structure for permanent magnet electric roller bearing |
| JP2023113280A (en) * | 2022-02-03 | 2023-08-16 | 大成建設株式会社 | Air leakage detection method |
| KR20230137021A (en) * | 2022-03-21 | 2023-10-04 | 주식회사 플로워크연구소 | detecting method and system of pipe abnormality using artificial intelligence |
| CN115099143A (en) * | 2022-06-22 | 2022-09-23 | 青岛理工大学 | Uncertainty reasoning-based driving all-electric actuator early fault diagnosis method |
| CN116320042A (en) * | 2023-05-16 | 2023-06-23 | 陕西思极科技有限公司 | Internet of things terminal monitoring control system for edge calculation |
| CN116738353A (en) * | 2023-08-15 | 2023-09-12 | 安拓思纳米技术（苏州）有限公司 | Pharmaceutical workshop air filter element performance detection method based on data analysis |
| CN116881745A (en) * | 2023-09-07 | 2023-10-13 | 深圳市特安工业科技有限公司 | Pressure transmitter abnormality monitoring method based on big data |
| CN117167903A (en) * | 2023-11-03 | 2023-12-05 | 江苏中安建设集团有限公司 | Artificial intelligence-based foreign matter fault detection method for heating ventilation equipment |
| CN117725792A (en) * | 2023-12-18 | 2024-03-19 | 清华大学合肥公共安全研究院 | Building safety monitoring and early warning method based on time sequence and model correction |
| CN117870943A (en) * | 2024-01-22 | 2024-04-12 | 中国三峡建工(集团)有限公司 | Multi-sensor-based data optimization acquisition system in grouting process |
Non-Patent Citations (6)
| Title |
|---|
| 文富勇;: "基于时间序列的地下水封洞库变形监测模型研究", 人民长江, no. 24, 28 December 2015 (2015-12-28) * |
| 朱少民 等: "基于ARIMA和LSTM组合模型的核电厂主泵状态预测", 核动力工程, 30 April 2022 (2022-04-30) * |
| 杨娟;戚建江;沈毅;: "ARIMA模型在杭州市中小学生咳嗽症状监测中的应用", 生物数学学报, no. 03, 15 September 2011 (2011-09-15) * |
| 王立华;孙少通;单恒;梁佳玮;: "基于ARIMA预测模型和混合补偿算法的医用耗材检漏仪应用研究", 山东科技大学学报(自然科学版), no. 03, 17 May 2019 (2019-05-17) * |
| 解建仓;王;雷社平;李想;吕正祥;: "基于ARIMA模型的大坝安全监测数据分析与预测", 人民黄河, no. 10, 10 October 2018 (2018-10-10) * |
| 谭海旺 等: "基于XGBoost-LSTM组合模型的光伏发电功率预测", 太阳能学报, 31 August 2022 (2022-08-31) * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN118090078B (en) | 2024-06-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111275288B (en) | XGBoost-based multidimensional data anomaly detection method and device | |
| CN108681633B (en) | Condensate pump fault early warning method based on state parameters | |
| Li et al. | A novel diagnostic and prognostic framework for incipient fault detection and remaining service life prediction with application to industrial rotating machines | |
| CN110648480B (en) | Single variable alarm system and method based on change rate | |
| CN102736546A (en) | State monitoring device of complex electromechanical system for flow industry and method | |
| CN109298633A (en) | Chemical production process fault monitoring method based on adaptive piecemeal Non-negative Matrix Factorization | |
| CN111855219A (en) | Diesel engine lubricating oil entering security parameter prediction method based on grey theory | |
| CN111881574A (en) | Wind turbine generator key component reliability modeling method based on distribution function optimization | |
| CN112414694A (en) | Equipment multistage abnormal state identification method and device based on multivariate state estimation technology | |
| CN117786371B (en) | Temperature monitoring data optimization prediction analysis method and system | |
| CN117195007A (en) | Heat exchanger performance prediction method and system | |
| CN118395091B (en) | Data efficient acquisition system based on everything interconnection | |
| Perry | Identifying the time of polynomial drift in the mean of autocorrelated processes | |
| CN117029968A (en) | Traffic data diagnosis method, system, storage medium and electronic equipment | |
| CN113887119A (en) | River water quality prediction method based on SARIMA-LSTM | |
| CN118211082A (en) | Oil filter element residual life prediction method and system based on data analysis | |
| Pan et al. | Multi-sensor signal fusion method for rolling bearing based on the standard relative mean–variance value and random weighting algorithm | |
| Yao et al. | A framework via impulses-oriented Gini index and extreme value distribution for rolling bearing dynamic fault alarm and identification | |
| CN118090078B (en) | Leakage online monitoring method for closed circulation water cooling system | |
| CN113469408A (en) | Running state trend prediction method and system for phase modulator | |
| CN116641941A (en) | Hydraulic system early fault dynamic detection method based on typical variable analysis | |
| CN117470529A (en) | Static friction fault detection method and system for valve of process control system | |
| CN114112390B (en) | Nonlinear complex system early fault diagnosis method | |
| CN110516323B (en) | Pressure pipeline damage prediction method based on time sequence analysis | |
| CN113283157A (en) | System, method, terminal and medium for predicting life cycle of intelligent stamping press part |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |