CN111259538B - Vibration alarm method, storage medium and electronic equipment - Google Patents
Vibration alarm method, storage medium and electronic equipment Download PDFInfo
- Publication number
- CN111259538B CN111259538B CN202010038460.9A CN202010038460A CN111259538B CN 111259538 B CN111259538 B CN 111259538B CN 202010038460 A CN202010038460 A CN 202010038460A CN 111259538 B CN111259538 B CN 111259538B
- Authority
- CN
- China
- Prior art keywords
- value
- equivalent
- vibration acceleration
- sequence
- determining
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 59
- 230000001133 acceleration Effects 0.000 claims abstract description 208
- 238000011161 development Methods 0.000 claims abstract description 47
- 238000012544 monitoring process Methods 0.000 claims description 20
- 238000012216 screening Methods 0.000 claims description 14
- 238000004590 computer program Methods 0.000 claims description 6
- 238000013461 design Methods 0.000 abstract description 10
- 238000010586 diagram Methods 0.000 description 18
- 230000002354 daily effect Effects 0.000 description 12
- 238000011156 evaluation Methods 0.000 description 10
- 230000008569 process Effects 0.000 description 8
- 238000004364 calculation method Methods 0.000 description 6
- 238000004891 communication Methods 0.000 description 6
- 238000012545 processing Methods 0.000 description 4
- 101100083446 Danio rerio plekhh1 gene Proteins 0.000 description 3
- 238000013016 damping Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000003203 everyday effect Effects 0.000 description 2
- 238000012423 maintenance Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000010485 coping Effects 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 230000008439 repair process Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
Landscapes
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
The embodiment of the invention provides a vibration alarm method, a storage medium and electronic equipment, wherein the corresponding equivalent value can be determined according to vibration acceleration data of measuring points, and a dynamic safety threshold value is determined according to the equivalent values accumulated one by one; determining an overrun probability value of the vibration acceleration according to the maximum vibration acceleration data of the measuring points, and further determining an overrun probability development trend coefficient according to the overrun probability value of the vibration acceleration; and determining an alarm state through the vibration acceleration equivalent value of the measuring point, the dynamic threshold value operation result, the effective value development trend coefficient and the overrun probability development trend coefficient. The vibration alarm method solves the problems that false alarm or untimely alarm is easy to cause when a design permission value is used as an alarm threshold value.
Description
Technical Field
The invention relates to the field of safety monitoring, in particular to a vibration alarm method, a storage medium and electronic equipment.
Background
The vibration safety assessment of the fan tower is a key problem of operation safety monitoring of the fan tower. The existing fan tower vibration safety assessment method is mainly used for roughly obtaining the tower vibration safety level based on comparison between the vibration monitoring parameter value and the design allowable value. Such an evaluation method causes a larger error in the evaluation result due to the following reasons: firstly, the designed safety allowable value is calculated mainly based on a simulation technology, and the error is larger; secondly, though the structural design of the same type of tower barrels is consistent, the construction errors, the foundation soil layers and the material characteristics are difficult to ensure to be completely consistent, so that the vibration resistance of the tower barrels with the same structure are different. Third, the damping of vibration resistance of the tower can gradually occur in the operation process, and if the design safety allowable value is always used as the safety evaluation threshold value, the damping is very unscientific.
Disclosure of Invention
In view of the above, an object of the present invention is to provide a vibration evaluation method, a storage medium, and an electronic device, which solve the problem that false alarm or untimely alarm is easily caused by using a design permission value as an alarm threshold.
In a first aspect, an embodiment of the present invention provides a vibration alarm method, including:
acquiring vibration acceleration data of at least one measuring point, and determining a vibration acceleration data sequence of the measuring point in a first period;
determining an equivalent value to be detected corresponding to a first period according to the vibration acceleration data sequence;
determining a sequence of equivalent values to be detected in a second period according to the equivalent values to be detected corresponding to the plurality of first periods;
dividing a to-be-measured value interval by taking the minimum value in the to-be-measured equivalent value sequence as a starting point and the maximum value as an end point, and dividing the to-be-measured value interval into a plurality of first intervals;
calculating the number of elements in the equivalent value sequence to be measured falling into each first interval;
calculating a first percentage, wherein the first percentage is the percentage of the number of elements falling into each first interval in an equivalent value sequence to be detected to the total number of elements in the equivalent value set;
screening at least one second interval of which the first percentage exceeds a preset percentage;
screening out a second interval with the maximum endpoint value as a safety threshold interval;
determining a dynamic threshold according to the starting point value and the ending point value of the safety threshold interval; and
and determining an alarm state according to the current vibration acceleration data and the dynamic threshold value.
Preferably, determining the equivalent value to be measured corresponding to the first period according to the vibration acceleration data sequence in the first period includes:
reading a vibration acceleration data sequence to be detected in a first period, and searching out the minimum value and the maximum value in the vibration acceleration data sequence to be detected;
determining a vibration acceleration value interval to be measured by taking the minimum value as a starting point and the maximum value as an end point, and dividing the vibration acceleration value interval to be measured into a plurality of third intervals;
calculating the number of the vibration acceleration data between the minimum value and the maximum value of the vibration acceleration data falling into each third interval;
screening out a third interval with the largest vibration acceleration data as an equivalent interval; and
and determining equivalent vibration acceleration data according to the starting point value and the end point value of the equivalent interval.
Preferably, the vibration acceleration data is vibration acceleration peak value or vibration acceleration effective value; and
and the equivalent vibration acceleration data are equivalent vibration acceleration peak value and equivalent vibration acceleration effective value.
Preferably, determining the sequence of equivalent values to be measured for the second period according to the equivalent values to be measured corresponding to the plurality of first periods includes:
accumulating the effective values of the equivalent vibration acceleration calculated in the first time period one by one to obtain a second time period equivalent vibration acceleration effective value set sequence; and
the dynamic threshold is equivalent vibration acceleration peak-to-peak dynamic threshold and equivalent vibration acceleration effective value dynamic threshold.
Preferably, determining the alarm state from the current vibration acceleration data and the dynamic threshold value comprises:
determining an equivalent value development trend coefficient according to the equivalent value sequence to be detected corresponding to each measuring point;
determining an overrun probability development trend coefficient according to the vibration acceleration data sequence and the dynamic threshold value;
determining a first security level according to the equivalent value to be detected of each current measuring point and the corresponding dynamic threshold value;
determining a second security level according to the development trend coefficient;
determining a third security level according to the overrun probability development trend coefficient; and
and determining the alarm state according to the first security level, the second security level and the third security level.
Preferably, determining the equivalent value development trend coefficient according to the equivalent value sequence to be measured corresponding to each measuring point includes:
determining an equivalent value overrun probability value of the first time period according to the equivalent value sequence of the first time period;
obtaining a set of equivalent value overrun probability values of the second time period according to the equivalent value overrun probability values of the first time period; and
and fitting elements in the set of equivalent value overrun probability values of the second period by using a least square method to generate a linear line segment, and calculating the slope of the line segment, wherein the slope of the line segment is the overrun probability development trend coefficient.
Preferably, determining the equivalent overrun probability value for the first time period from the equivalent sequence of the first time period comprises:
according to a set of equivalent values to be measured in a third period, calculating the maximum value of the equivalent values to be measured in the set of equivalent values to be measured in the third period;
obtaining an equivalent value sequence of the first period according to the maximum value of the equivalent value to be detected;
reading a comparison value, wherein the number of elements in the equivalent value sequence of the first period is taken as the comparison value;
determining a first number according to the comparison value and the dynamic threshold value; and
and determining an equivalent value overrun probability value of the first period according to the first quantity and the comparison value.
Preferably, the equivalent value sequence is vibration acceleration peak value and vibration acceleration effective value;
the equivalent value overrun probability value is a vibration acceleration peak value overrun probability value and a vibration acceleration effective value overrun probability value;
the overrun probability development trend coefficient is a vibration acceleration peak value overrun probability development trend coefficient and a vibration acceleration equivalent value overrun probability development trend coefficient; and
and determining the overrun probability development trend coefficient according to the overrun probability development trend coefficient of the vibration acceleration peak value and the overrun probability development trend coefficient of the vibration acceleration effective value.
In a second aspect, a computer readable storage medium is provided for storing computer program instructions which, when executed by a processor, implement a method as described above.
In a third aspect, there is provided an electronic device comprising a memory and a processor, wherein the memory is for storing one or more computer program instructions, wherein the one or more computer program instructions are executed by the processor to implement the method as described above.
Thus, on the one hand, the current dynamic threshold is updated daily according to the vibration parameters, and a proper vibration safety threshold interval is automatically calculated based on historical vibration data of the fan. In this embodiment, the current safety threshold of the monitoring system is automatically corrected as the monitoring data increases daily, compared to the mode in which the vibration design allowable value is fixed as the safety threshold. The vibration safety threshold interval calculated through the embodiment is gradually increased along with the monitored sample data, and the vibration safety threshold interval is more and more close to the current practical vibration allowable value of the fan tower, so that the stability of vibration alarm data is improved. On the other hand, the embodiment of the invention calculates the safety evaluation state based on the daily development trend coefficient and the overrun probability coefficient of the vibration acceleration. Compared with the technical scheme that the design allowable value is used and compared with the safety threshold value in the prior art, the technical scheme that the compared result is directly used as the alarm state can solve the problems of inaccurate alarm and false alarm of the current fan vibration alarm method, obviously improve the accuracy of fan tower vibration safety early warning, further reduce the fan tower collapse accident and reduce the daily maintenance cost.
Drawings
The above and other objects, features and advantages of the present invention will become more apparent from the following description of embodiments of the present invention with reference to the accompanying drawings, in which:
FIG. 1 is a schematic diagram of a vibration alarm system of a blower in accordance with an embodiment of the present invention;
FIG. 2 is a schematic diagram of an alarm monitoring system for a blower in accordance with an embodiment of the present invention;
FIG. 3 is a flow chart of a vibration alert method of an embodiment of the present invention;
FIG. 4 is a schematic diagram of a process of reading vibration acceleration data by a computer according to an embodiment of the present invention;
FIG. 5 is a data flow diagram of an equivalent vibratory acceleration peak-to-peak calculation process according to an alternative embodiment of the invention;
FIG. 6 is a flow chart of equivalent acceleration data processing according to an embodiment of the present invention;
FIG. 7 is a schematic diagram of a process of screening equivalent vibration acceleration data according to an embodiment of the present invention;
FIG. 8 is a schematic diagram of equivalent vibration acceleration data processed during a second period of time in accordance with an embodiment of the present invention;
FIG. 9 is a data flow diagram of a current dynamic threshold calculation method according to an embodiment of the present invention;
FIG. 10 is a diagram illustrating screening of security zones according to an embodiment of the present invention;
FIG. 11 is a flow chart of an embodiment of the present invention for determining an alarm state;
FIG. 12 is a flow chart of determining an equivalent value trend coefficient according to an embodiment of the present invention;
FIG. 13 is a flow chart of determining overrun probability trend coefficients according to an embodiment of the present invention;
fig. 14 is a block diagram of the electronic device of the embodiment of the present invention.
Detailed Description
The present invention is described below based on examples, but the present invention is not limited to only these examples. In the following detailed description of the present invention, certain specific details are set forth in detail. The present invention will be fully understood by those skilled in the art without the details described herein. Well-known methods, procedures, flows, components and circuits have not been described in detail so as not to obscure the nature of the invention.
Moreover, those of ordinary skill in the art will appreciate that the drawings are provided herein for illustrative purposes and that the drawings are not necessarily drawn to scale.
Unless the context clearly requires otherwise, throughout the description and the claims, the words "comprise", "comprising", and the like are to be construed in an inclusive sense as opposed to an exclusive or exhaustive sense; that is, it is the meaning of "including but not limited to".
In the description of the present invention, it should be understood that the terms "first," "second," and the like are used for descriptive purposes only and are not to be construed as indicating or implying relative importance. Furthermore, in the description of the present invention, unless otherwise indicated, the meaning of "a plurality" is two or more.
The following describes the technical scheme in the embodiment of the present invention in detail with reference to the drawings.
FIG. 1 is a schematic diagram of a vibration alarm system of a blower in accordance with an embodiment of the present invention. As shown in fig. 1, the vibration sensors 1 to n are mounted on the tower of the blower, and the collected data may be displacement, velocity or acceleration according to the type of the sensor. The sensor is connected to the data acquisition module 20 by a cable. The data acquisition module 20 is connected with the data communication module 21 through a bus. The sensors 1-n, the data acquisition module 20 and the data communication module 21 are installed in the fan tower. The fan vibration monitoring module 30 operates on a computer (which may be an industrial computer or a server) in an industrial control room and is in communication connection with the data communication module through a cable, an optical cable or a wireless signal.
FIG. 2 is a schematic diagram of an alarm monitoring system for a blower in accordance with an embodiment of the present invention. As illustrated in FIG. 2, the blower vibration monitoring module 30 performs the function of alarm monitoring. In the prior art, the safety threshold 30 is a constant that is fixed and calculated based on the design allowable value of the fan. For example, in the prior art, if the set safety threshold 30 is low, the vibration parameter 31 is compared with the low safety threshold 30, so that false alarm information 34 may be generated, and after the manager finds the alarm information 34, the manager needs to go to repair before the time, and a large amount of false alarm information 34 may cause the manager to feel tired of coping with the situation, and affect the normal management work of the fan. In contrast, in the prior art, if the safety threshold 30 set by the fan is higher, if the fan is in a dangerous state at this time, since the vibration parameter 31 is smaller than the higher safety threshold 30, no alarm information 34 is generated, a serious accident is missed, and the fan is more seriously damaged due to the continuous operation of the fan, which affects the service life of the fan, and even generates a serious accident. In addition, because the wind power station has a severe environment, important parts are located at a high altitude of tens of meters, and the wind speed change is irregular, so that different wind power sets have own vibration characteristics, and the conventional fan vibration detection technology is easy to cause the problem of false alarm or missing alarm. In the embodiment of the invention, the safety threshold 30 is obtained by daily updating of historical operation data of the fan, and the alarm monitor 33 further calculates the alarm information 34 by the dynamic safety threshold 30 obtained by daily updating, so that more accurate alarm information 34 can be given to the manager.
Fig. 3 is a flow chart of a vibration alert method of an embodiment of the present invention. As shown in fig. 3, the vibration alarm method of the present embodiment includes the steps of:
and step S100, obtaining vibration acceleration data of at least one measuring point to determine a vibration acceleration data sequence of the measuring point in a first period. In an alternative embodiment, a plurality of acceleration sensors are installed on the fan tower, and vibration acceleration peak value and vibration acceleration effective value of the measuring point are obtained. The vibration acceleration peak-to-peak value is collected and stored by the data collection module 20 to form a vibration acceleration peak-to-peak value sequence. In the same way, a sequence of vibration acceleration effective values can be obtained. And a computer positioned in the industrial control room reads the vibration acceleration peak-to-peak value sequence and the vibration acceleration effective value sequence through a data network.
And 200, determining an equivalent value to be detected corresponding to the first period according to the vibration acceleration data sequence. In an alternative embodiment, the vibration acceleration data sequence is a vibration acceleration peak-to-peak data sequence and a vibration acceleration effective value data sequence. The first period is a period length of one acquisition cycle per day. The step of determining the equivalent value to be measured corresponding to the first period by the vibration acceleration data sequence in the first period is to convert the vibration acceleration peak value data sequence into an equivalent vibration acceleration peak value of the same day according to the collected vibration acceleration peak value data sequence, and convert the vibration acceleration peak value data sequence into an equivalent vibration acceleration peak value. In the same method, the vibration acceleration effective value data sequence is converted into an equivalent vibration acceleration effective value. A preferred embodiment of this embodiment is to start acquisition from 0:00 of the day, to stop acquisition at 23:59, and to designate the acquisition period of 0:00 to 23:29 of the day as the first period.
FIG. 4 is a schematic diagram of a process of reading vibration acceleration data by a computer according to an embodiment of the present invention. As shown in fig. 4, the computer reads the vibration acceleration data at 1 minute acquisition intervals. For example 0 during the day: at time 00, the acquired vibration acceleration data is A1; vibration acceleration data acquired at the time of 0:01 is A2; similarly, at time 23:59, the vibration acceleration data collected is A1440. Thus, 1440 vibration acceleration data may be acquired during a first time period. In an alternative embodiment, the vibration acceleration data is vibration acceleration peak-to-peak data.
FIG. 5 is a data flow diagram of an equivalent vibration acceleration peak-to-peak calculation process of an alternative embodiment. As shown in fig. 5, 1440 vibration acceleration peak-to-peak data are sent to the equivalent value calculation module to obtain an equivalent vibration acceleration peak-to-peak value on the same day. In the same way, the vibration acceleration data is vibration acceleration effective value data. An equivalent vibration acceleration effective value of the same day can be obtained.
FIG. 6 is a flow chart of equivalent acceleration data processing according to an embodiment of the present invention. As shown in fig. 6, the step S200 may include the steps of:
step S210, a vibration acceleration data sequence to be detected in a first period is read, and the minimum value and the maximum value in the vibration acceleration data sequence to be detected are searched.
And the computer reads the vibration acceleration data sequence to be measured in the first period. When the first period is one day, step S210 is to read 1440 vibration acceleration data of the same day. In the following description, 1440 vibration acceleration data are taken as an example of reading vibration acceleration data once every minute and collecting vibration acceleration data every day. It is easy to understand that a person skilled in the art can set different sampling intervals and time period lengths corresponding to the vibration acceleration data sequences according to the needs of the application scenario and the processing capability of the configured hardware device.
Specifically, the minimum value and the maximum value in the vibration acceleration data sequence to be detected are searched. And marking the minimum value in the vibration acceleration data sequence to be measured as min, and marking the maximum value as max.
And S220, determining a vibration acceleration value interval to be measured by taking the minimum value as a starting point and the maximum value as an end point, and dividing the vibration acceleration value interval to be measured into a plurality of third intervals.
Specifically, the section of the vibration acceleration data is divided with the minimum value min as a start point and the maximum value max as an end point. The vibration acceleration data section is divided into a plurality of third sections according to the specified section interval.
In a preferred embodiment, the vibration acceleration data interval is divided into 5 third intervals.
The interval of the third interval is:
the 5 third intervals are respectively:
[min，min+Δ]，
[min+Δ，min+2Δ]，
[min+2Δ，min+3Δ]，
[min+3Δ，min+4Δ]，
[min+4Δ，min+5Δ]
step 230, calculating the number of the vibration acceleration data between the minimum value and the maximum value of the vibration acceleration data falling into each third interval; and calculating the interval with the maximum number of vibration acceleration peak-to-peak value sequences.
Fig. 7 is a schematic diagram of a process of screening equivalent vibration acceleration data according to an embodiment of the present invention. As shown in fig. 7, each third section is first named according to the section interval, for example, "third section-1", "third section-2", and so on. And determining the starting point and the ending point of each third section, wherein two data columns, such as a section starting point and a section ending point, respectively correspond to the starting point value and the ending point value of each third section. For example, the number of data falling into the third section "third section-1" is 288, the number of data falling into the third section "2" is 300, the number of data falling into the third section "3" is 260, the number of data falling into the third section "4" is 330, and the number of data falling into the third section "5" is 262.
Step 240, screening out the third interval with the largest falling vibration acceleration data as the equivalent interval
As shown in fig. 7, the number of data falling into the "third section-4" is 330, and it can be found by comparison that the number of data falling into it is the largest, and therefore, it is determined that the "third section-4" is the section that "falls into the largest number". Then "third interval-4" is the screened equivalent interval.
And step 250, determining equivalent vibration acceleration data according to the starting point value and the end point value of the equivalent interval.
Taking the equivalent intermediate value as the equivalent vibration acceleration peak value of the current day. The calculation steps are as follows: intermediate value of equivalent section= (section start+section end)/2. As shown in fig. 7, the intermediate value of the equivalent interval= (min1+3Δ1+min1+4Δ1)/2.
Step 200 is performed daily to generate equivalent vibration acceleration data for a day, according to the method shown in fig. 6. When the read vibration acceleration data is a vibration acceleration peak value, the equivalent vibration acceleration peak value of the same day is obtained; when the read vibration acceleration data is a vibration acceleration effective value, the equivalent vibration acceleration effective value of the same day is obtained. In a preferred embodiment, step 200 is performed once a day, and Y days later, Y equivalent vibration acceleration data are generated, which are a plurality of equivalent values to be measured.
Step 300, determining a sequence of equivalent values to be measured in the second period according to the equivalent values to be measured corresponding to the first periods.
The equivalent value to be measured corresponding to the first period is obtained through the step 200, and a plurality of equivalent vibration acceleration data obtained through the step 200 form a sequence of equivalent values to be measured in the second period. In a preferred embodiment, the sequence of equivalent values to be measured for the second period refers to a plurality of equivalent vibration acceleration data obtained since the start of monitoring.
Fig. 8 is a schematic diagram of equivalent vibration acceleration data processed in the second period of time according to the embodiment of the present invention. As shown in fig. 8, an equivalent vibration acceleration data 42 is generated every day since the start of monitoring, and the vibration acceleration data 42 is calculated by step 200. For example, on the first day, an equivalent value eq1 is calculated; the next day, an equivalent value eq2 is calculated; similarly, on day Y, an equivalent value eqi (eqi, i=1, 2,3 … …, Y) is calculated, and after several days, Y equivalent vibration acceleration data from the start of monitoring are generated. And starting to monitor Y equivalent vibration acceleration data, namely a sequence of equivalent values to be detected in the second period.
FIG. 9 is a data flow diagram of a current dynamic threshold calculation method according to an embodiment of the present invention. As shown in fig. 9, multiple equivalent values result in a current dynamic threshold. The process of calculating the current dynamic threshold is described in detail in steps 400 through 900.
And 400, dividing a to-be-measured value interval by taking the minimum value in the to-be-measured equivalent value sequence of the second period as a starting point and the maximum value as an end point, and dividing the to-be-measured value interval into a plurality of first intervals.
As shown in fig. 9, the daily equivalent vibration acceleration data 41 from the start of monitoring is read. In a preferred embodiment, the equivalent vibration acceleration data 41 is an equivalent vibration acceleration peak-to-peak value or an equivalent vibration acceleration effective value. First, the vibration acceleration data is an equivalent vibration acceleration peak-to-peak value, and will be described in detail. And screening out the maximum value and the minimum value of the peak-to-peak value sequence of the equivalent vibration acceleration. The maximum value in the equivalent vibration acceleration peak-to-peak value sequence is max1, and the minimum value in the equivalent vibration acceleration peak-to-peak value sequence is min1.
The section of the vibration acceleration peak-to-peak value is divided by taking the minimum value min1 as a starting point and taking the maximum value max1 as an end point.
10 sections are re-divided according to the sections of the vibration acceleration peak value. In a preferred embodiment, the interval between the minimum value min1 and the maximum value max1 is uniformly divided into 10 first intervals, and each interval is:
the 10 first intervals are respectively:
[min1，min1+Δ1]，
[min1+Δ1，min1+2Δ1]，
[min1+2Δ1，min1+3Δ1]，
[min1+3Δ1，min1+4Δ1]，
[min1+4Δ1，min1+5Δ1]，
[min1+5Δ1，min1+6Δ1]，
[min1+6Δ1，min1+7Δ1]，
[min+17Δ1，min1+8Δ1]，
[min1+8Δ1，min1+9Δ1]，
[min1+9Δ1，min+10Δ1]，
in the same method, when the vibration acceleration data is an equivalent vibration acceleration effective value, 10 sections can be divided again for the vibration acceleration effective value.
Step 500, calculating the number of elements in the equivalent value sequence to be measured falling into each first interval. Fig. 10 is a schematic diagram of screening a security interval according to an embodiment of the present invention. As shown in fig. 10, the number of data falling into the interval of the measurement points is counted first, for example, in the second period and 300 equivalent values are generated, and in the data line with the sequence number of 1, the number of data falling into the interval is 30.
Step 600, calculating a first percentage, where the first percentage is a percentage of the number of elements in the equivalent value sequence to be measured falling into each first interval to the total number of elements in the equivalent value sequence.
As shown in fig. 10, in the data line with the sequence number of 1, the number of data falling is 30, and the first percentage is 30/300=10%.
Step 700, screening out at least one second interval with the first percentage exceeding the predetermined percentage.
As shown in fig. 10, a first percentage in each data line with sequence numbers 1 to 10 is calculated, and the first percentage of each data line is compared with a predetermined percentage. If the first percentage is greater than 10%, the section consisting of the section start point and the section end point in the line of data is marked as the second section. For example, the second column of data row number 3 is "first interval-3", whose first percentage is 10%, and which does not exceed the first percentage, is not the second interval; the second column of data line number 4 is "first section-4", the first percentage of which is 11%, exceeding the first percentage, so that the section consisting of the start point of this data line and the end point of the data line is marked as the second section. By the above-mentioned method, the data lines where the first section-4, the first section-6, the first section-7 and the first section-9 are located are the second sections corresponding to the sections formed by the data start point and the data end point.
Step 800, screening out the second interval with the largest endpoint value as the safety threshold interval.
As shown in fig. 10, all the data rows with the largest end points of the intervals in the second interval are screened out. For example, the data line with the sequence number 9 is marked as the second section, and it is the data line in the second section in which the section end value is the largest. Thus, a data line with a sequence number of 9 is selected, and the start point and end point of the section constitute a safety threshold section.
And step 900, determining a dynamic threshold according to the starting point value and the ending point value of the safety threshold interval.
And calculating the intermediate value as the current dynamic threshold value of the vibration equivalent acceleration data according to the intermediate value of the safety threshold value interval. If the equivalent value to be measured uses an equivalent vibration acceleration peak value, obtaining a current dynamic threshold value of the equivalent vibration acceleration peak value; in the same method, the effective value of the equivalent vibration acceleration is taken as the equivalent value to be measured, and the current dynamic threshold value of the effective value of the equivalent vibration acceleration is obtained.
Step 1000, determining an alarm state according to the current vibration acceleration data and the dynamic threshold value. FIG. 11 is a flow chart of an embodiment of the present invention for determining an alarm condition.
As shown in fig. 11, step 1000 includes the following sub-steps:
step 1100, determining an equivalent value development trend coefficient according to the equivalent value sequence to be measured corresponding to each measuring point. Fig. 12 is a flowchart for determining an equivalent value development trend coefficient according to an embodiment of the present invention. As shown in fig. 12, step 1100 includes the following two sub-steps:
at step 1110, a linear segment is fit using least squares.
According to the time sequence, the vibration acceleration peak value after starting monitoring is formed into a time sequence, the data in the time sequence is fitted according to a least square method, a linear line segment can be obtained, and the slope coefficient of the linear line segment obtained through fitting is calculated.
In the same method, the vibration acceleration effective value from the beginning of monitoring is formed into a time sequence, the data in the sequence is fitted according to a least square method, a linear line segment can be obtained, and the slope coefficient of the linear line segment obtained by fitting is calculated.
In step 1120, the slope coefficient of the line segment is the vibration acceleration peak-to-peak trend coefficient and the vibration acceleration effective value super trend coefficient of the analyzed measurement point, respectively.
And step 1200, determining an overrun probability development trend coefficient according to the vibration acceleration data sequence and the dynamic threshold value.
Fig. 13 is a flowchart for determining an overrun probability trend coefficient according to an embodiment of the present invention. As shown in fig. 13, step 1200 includes the following three sub-steps:
step 1210 determines an equivalent overrun probability value for the second time period based on the equivalent sequence of the first time period, and in a preferred embodiment, the equivalent of the first time period is maximum vibration acceleration data for each minute of the first time period, which may be a maximum vibration acceleration peak-to-peak value or a maximum vibration acceleration effective value.
Step 1210 includes the steps of: according to the set of equivalent values to be measured in the third period, calculating the maximum value of the equivalent values to be measured in the set of equivalent values to be measured in the third period; obtaining an equivalent value sequence of the first period according to the maximum value of the equivalent value to be detected; reading a comparison value, wherein the number of elements in the equivalent value sequence of the first period is taken as the comparison value; determining a first number according to the comparison value and the dynamic threshold value; and determining an equivalent value overrun probability value of the first period according to the first quantity and the comparison value.
The vibration acceleration peak-to-peak value will be described below as an example. Specifically, according to a set of equivalent values to be measured in a third period, calculating the maximum value of the equivalent values to be measured in the set of equivalent values to be measured in the third period; the third time period is the time length of one minute, and the maximum vibration acceleration peak-to-peak value of each minute is calculated. 1440 data are formed during the day. Obtaining an equivalent value sequence of the first period according to the maximum value of the equivalent value to be detected; and obtaining a maximum vibration acceleration peak-to-peak value sequence of the first time period. The comparison value is read, the number of elements in the equivalent value sequence of the first period is taken as the comparison value, specifically 1440 maximum vibration acceleration peak values are read in the first period, and the 1440 maximum vibration acceleration peak values are the comparison value. Determining a first number according to the comparison value and the dynamic threshold value; specifically, the number of maximum vibration acceleration peaks exceeding the dynamic threshold (derived in step 900) is a first number. And determining the vibration acceleration peak-to-peak overrun probability value of the first period according to the first quantity and the comparison value (1440). Specifically, the ratio of the first number to 1440 is the vibration acceleration peak-to-peak overrun probability value of the first period.
And similarly, calculating the overrun probability value of the vibration acceleration effective value in the first period.
Step 1220 obtains a set of equivalent overrun probability values for the second time period based on the equivalent overrun probability values for the first time period.
And forming a set by using the equivalent vibration acceleration peak value overrun probability values of all the first time periods, or forming a set by using the equivalent vibration acceleration effective value overrun probability values of all the first time periods.
Step 1230 uses least square method to fit the elements in the set of equivalent value overrun probability values of the second period to generate a linear line segment, and calculates the slope of the line segment, where the slope of the line segment is the overrun probability development trend coefficient.
In an optional embodiment, the equivalent value overrun probability set of the second period is an equivalent vibration acceleration peak-to-peak overrun probability value set or an equivalent vibration acceleration effective value overrun probability value set, elements in the set are fitted to generate a linear line segment, and the slope of the line segment is calculated to obtain a vibration acceleration peak-to-peak overrun probability development trend coefficient or a vibration acceleration effective value overrun probability development trend coefficient respectively. And calculating the over-limit probability development trend coefficient of the peak value of the vibration acceleration peak or the average value of the over-limit probability development trend coefficient of the vibration acceleration effective value as the over-limit probability development trend coefficient.
Step 1300, determining a first security level according to the equivalent value to be detected of each current measuring point and the corresponding dynamic threshold value;
and determining a first security level according to the equivalent value to be detected of each current measuring point and the corresponding dynamic threshold value. The first security level determination method is as follows: taking the quotient of the peak value of the equivalent vibration acceleration and the current dynamic threshold value of the peak value of the equivalent vibration acceleration as a first proportional coefficient. When the first coefficient is smaller than 1, the safety level is equal to 0, when the first coefficient is larger than or equal to 1 and smaller than 2, the safety level is equal to 1, when the first coefficient is larger than or equal to 2 and smaller than 4, the safety level is equal to 2, and when the first coefficient is larger than or equal to 4, the safety level is equal to 3.
Another method of first security level determination is: and taking the quotient of the effective value of the equivalent vibration acceleration and the current dynamic threshold value of the effective value of the equivalent vibration acceleration as a second proportionality coefficient. When the second coefficient is smaller than 1, the safety level is equal to 0, when the second coefficient is larger than or equal to 1 and smaller than 2, the safety level is equal to 1, when the second coefficient is larger than or equal to 2 and smaller than 3, the safety level is equal to 2, and when the second coefficient is larger than or equal to 3, the safety level is equal to 3.
And step 1400, determining a third security level according to the overrun probability development trend coefficient.
The third security level determination method is as follows: taking the overrun probability development trend coefficient as a fifth proportion coefficient; when the fifth coefficient is less than or equal to 0.1, the safety level is equal to 0, when the fifth coefficient is greater than or equal to 0.1 and less than 0.5, the safety level is equal to 1, when the fifth coefficient is greater than or equal to 0.5 and less than 1, the safety level is equal to 2, and when the fifth coefficient is greater than or equal to 1, the safety level is equal to 3.
Step 1500, determining the alarm state according to the first security level, the second security level and the third security level. And determining a third security level according to the overrun probability development trend coefficient.
And determining a second security level according to the development trend coefficient. The second security level determination method is as follows: taking the vibration acceleration peak value development trend coefficient as a third proportion coefficient; and when the third coefficient is smaller than or equal to 0.1, the safety level is equal to 0, when the third coefficient is larger than or equal to 0.1 and smaller than 0.5, the safety level is equal to 1, when the third coefficient is larger than or equal to 0.5 and smaller than 1, the safety level is equal to 2, and when the third coefficient is larger than or equal to 1, the safety level is equal to 3.
Another method of second security level determination is: taking the development trend coefficient of the vibration acceleration effective value as a fourth proportion coefficient; and when the fourth coefficient is smaller than or equal to 0.1, the safety level is equal to 0, when the fourth coefficient is larger than or equal to 0.1 and smaller than 0.5, the safety level is equal to 1, when the fourth coefficient is larger than or equal to 0.5 and smaller than 1, the safety level is equal to 2, and when the fourth coefficient is larger than or equal to 1, the safety level is equal to 3.
And determining the alarm state according to the first security level, the second security level and the third security level. The security level is equal to 0, the evaluation state is the first state, and the information sent to the manager is 'security'; the security level is equal to 1, the evaluation state is the second state, and the information sent to the manager is 'attention'; the security level is 2, the evaluation state is the third state, and the information sent to the manager is a warning; the security level is 3, the evaluation state is the fourth state, and the information sent to the manager is an alarm.
Fig. 14 is a schematic diagram of an electronic device according to an embodiment of the invention. As shown in fig. 14, in this embodiment, the electronic device may be a server or a terminal, and the terminal may be an intelligent device such as a mobile phone, a computer, a tablet computer, or the like. As shown, the electronic device includes: at least one processor 141; a memory 140 communicatively coupled to the at least one processor; and a communication component 142 communicatively coupled to the storage medium, the communication component 142 receiving and transmitting data under control of the processor; the memory 140 stores instructions that may be executed by the at least one processor 141, and the instructions are executed by the at least one processor 141 to implement the task allocation method according to the embodiment of the present invention.
In particular, the memory 140 is used as a non-volatile computer readable storage medium for storing non-volatile software programs, non-volatile computer executable programs, and modules. The processor 141 executes various functional applications of the device and data processing, i.e., implements the task allocation method described above, by running nonvolatile software programs, instructions, and modules stored in the memory.
Memory 140 may include a storage program area that may store an operating system, at least one application program required for functionality, and a storage data area; the storage data area may store a list of options, etc. In addition, memory 140 may include high-speed random access memory, and may also include non-volatile memory, such as at least one magnetic disk storage device, flash memory device, or other non-volatile solid-state storage device. In some embodiments, memory 140 optionally includes memory remotely located relative to processor 141, which may be connected to an external device via a network. Examples of such networks include, but are not limited to, the internet, intranets, local area networks, mobile communication networks, and combinations thereof.
One or more modules are stored in memory 140 that, when executed by one or more processors 141, perform the task allocation method of any of the method embodiments described above.
The product may perform the method disclosed in the embodiment of the present application, and have corresponding functional modules and beneficial effects of the performing method, and technical details not described in detail in the embodiment of the present application may be referred to the method disclosed in the embodiment of the present application.
The invention also relates to a computer readable storage medium for storing a computer readable program for causing a computer to perform some or all of the above-described method embodiments.
That is, it will be understood by those skilled in the art that all or part of the steps in implementing the methods of the embodiments described above may be implemented by a program, which is stored in a storage medium and includes several instructions for causing a device (which may be a single-chip microcomputer, a chip or the like) or a processor (processor) to perform all or part of the steps in the methods of the embodiments described herein. And the aforementioned storage medium includes: a U-disk, a removable hard disk, a Read-Only Memory (ROM), a random access Memory (RAM, random Access Memory), a magnetic disk, or an optical disk, or other various media capable of storing program codes.
Thus, on the one hand, the current dynamic threshold is updated daily according to the vibration parameters, and a proper vibration safety threshold interval is automatically calculated based on historical vibration data of the fan. In this embodiment, the current safety threshold of the monitoring system is automatically corrected as the monitoring data increases daily, compared to the mode in which the vibration design allowable value is fixed as the safety threshold. The vibration safety threshold interval calculated through the embodiment is gradually increased along with the monitored sample data, and the vibration safety threshold interval is more and more close to the current practical vibration allowable value of the fan tower, so that the stability of vibration alarm data is improved. On the other hand, the embodiment of the invention calculates the safety evaluation state based on the daily development trend coefficient and the overrun probability coefficient of the vibration acceleration. Compared with the technical scheme that the design allowable value is used and compared with the safety threshold value in the prior art, the technical scheme that the compared result is directly used as the alarm state can solve the problems of inaccurate alarm and false alarm of the current fan vibration alarm method, obviously improve the accuracy of fan tower vibration safety early warning, further reduce the fan tower collapse accident and reduce the daily maintenance cost.
The above description is only of the preferred embodiments of the present invention and is not intended to limit the present invention, and various modifications and variations may be made to the present invention by those skilled in the art. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present invention should be included in the protection scope of the present invention.
Claims (7)
1. A vibration alarm method, comprising:
acquiring vibration acceleration data of at least one measuring point, and determining a vibration acceleration data sequence of the measuring point in a first period;
determining an equivalent value to be detected corresponding to a first period according to the vibration acceleration data sequence;
determining a sequence of equivalent values to be detected in a second period according to the equivalent values to be detected corresponding to the plurality of first periods;
dividing a to-be-measured value interval by taking the minimum value in the to-be-measured equivalent value sequence as a starting point and the maximum value as an end point, and dividing the to-be-measured value interval into a plurality of first intervals;
calculating the number of elements in the equivalent value sequence to be measured falling into each first interval;
calculating a first percentage, wherein the first percentage is the percentage of the number of elements in an equivalent value sequence to be detected falling into each first interval to the total number of elements in the equivalent value sequence;
screening at least one first area with the first percentage exceeding a preset percentage as a second interval;
screening out a second interval with the maximum endpoint value as a safety threshold interval;
determining a dynamic threshold according to the starting point value and the ending point value of the safety threshold interval; and
determining an alarm state according to the current vibration acceleration data and the dynamic threshold value;
wherein, determining the equivalent value to be measured corresponding to the first period according to the vibration acceleration data sequence includes:
reading a vibration acceleration data sequence to be detected in a first period, and searching out the minimum value and the maximum value in the vibration acceleration data sequence to be detected;
determining a vibration acceleration value interval to be measured by taking the minimum value as a starting point and the maximum value as an end point, and dividing the vibration acceleration value interval to be measured into a plurality of third intervals;
calculating the number of the vibration acceleration data between the minimum value and the maximum value of the vibration acceleration data falling into each third interval;
screening out a third interval with the largest vibration acceleration data as an equivalent interval; and
and determining equivalent vibration acceleration data according to the starting point value and the end point value of the equivalent interval.
2. The vibration alert method according to claim 1, wherein the vibration acceleration data is a vibration acceleration peak value or a vibration acceleration effective value; and
and the equivalent vibration acceleration data is equivalent vibration acceleration peak value or equivalent vibration acceleration effective value.
3. The vibration alert method according to claim 2, wherein determining the sequence of equivalent values to be measured for the second period from the plurality of equivalent values to be measured for the first period comprises:
accumulating the effective values of the equivalent vibration acceleration calculated in the first time period one by one to obtain a second time period equivalent vibration acceleration effective value set sequence; and
the dynamic threshold is an equivalent vibration acceleration effective value dynamic threshold.
4. The vibration alert method according to claim 1, wherein determining an alert state based on current vibration acceleration data and the dynamic threshold value comprises:
determining an equivalent value development trend coefficient according to the equivalent value sequence to be detected corresponding to each measuring point;
determining an overrun probability development trend coefficient according to the vibration acceleration data sequence and the dynamic threshold value;
determining a first security level according to the equivalent value to be detected of each current measuring point and the corresponding dynamic threshold value;
determining a second security level according to the development trend coefficient;
determining a third security level according to the overrun probability development trend coefficient; and
determining the alarm state according to the first security level, the second security level and the third security level;
wherein, determining the equivalent value development trend coefficient according to the equivalent value sequence to be measured corresponding to each measuring point comprises:
forming a time sequence of vibration acceleration peaks from the beginning of monitoring according to the time sequence, fitting the data in the time sequence according to a least square method to obtain a first linear line segment, and calculating the slope coefficient of the first linear line segment;
forming a time sequence of vibration acceleration effective values from the beginning of monitoring, fitting the data in the time sequence according to a least square method to obtain a second linear line segment, and calculating the slope coefficient of the second linear line segment;
taking the slope coefficient of the first linear line segment as a vibration acceleration peak value development trend coefficient and taking the slope coefficient of the second linear line segment as a vibration acceleration effective value super development trend coefficient;
wherein determining the overrun probability development trend coefficient according to the vibration acceleration data sequence and the dynamic threshold value comprises:
determining an equivalent value overrun probability value of the first time period according to the equivalent value sequence of the first time period;
obtaining a set of equivalent value overrun probability values of the second time period according to the equivalent value overrun probability values of the first time period; and
fitting elements in the set of equivalent value overrun probability values of the second period by using a least square method to generate a linear line segment, calculating the slope of the line segment, and taking the slope of the line segment as an overrun probability development trend coefficient;
wherein determining the equivalent overrun probability value for the first time period based on the equivalent sequence of the first time period comprises:
according to a set of equivalent values to be measured in a third period, calculating the maximum value of the equivalent values to be measured in the set of equivalent values to be measured in the third period;
obtaining an equivalent value sequence of the first period according to the maximum value of the equivalent value to be detected;
reading a comparison value, and taking the value of an element in the equivalent value sequence of the first time period as the comparison value;
determining a number of comparison values exceeding the dynamic threshold as a first number; and
and determining an equivalent value overrun probability value of the first time period according to the first quantity and the quantity of elements in the equivalent value sequence of the first time period.
5. The vibration alert method according to claim 4, wherein the equivalent value sequence is a vibration acceleration peak-to-peak value and a vibration acceleration effective value;
the equivalent value overrun probability value is a vibration acceleration peak value overrun probability value and a vibration acceleration effective value overrun probability value;
the overrun probability development trend coefficient is a vibration acceleration peak value overrun probability development trend coefficient and a vibration acceleration equivalent value overrun probability development trend coefficient; and
and determining the overrun probability development trend coefficient according to the overrun probability development trend coefficient of the vibration acceleration peak value and the overrun probability development trend coefficient of the vibration acceleration effective value.
6. A computer readable storage medium storing computer program instructions which, when executed by a processor, implement the method of any one of claims 1-5.
7. An electronic device comprising a memory and a processor, wherein the memory is configured to store one or more computer program instructions, wherein the one or more computer program instructions are executed by the processor to implement the method of any of claims 1-5.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010038460.9A CN111259538B (en) | 2020-01-14 | 2020-01-14 | Vibration alarm method, storage medium and electronic equipment |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010038460.9A CN111259538B (en) | 2020-01-14 | 2020-01-14 | Vibration alarm method, storage medium and electronic equipment |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111259538A CN111259538A (en) | 2020-06-09 |
| CN111259538B true CN111259538B (en) | 2023-08-01 |
Family
ID=70945257
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010038460.9A Active CN111259538B (en) | 2020-01-14 | 2020-01-14 | Vibration alarm method, storage medium and electronic equipment |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111259538B (en) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111636509A (en) * | 2020-07-17 | 2020-09-08 | 上海三一重机股份有限公司 | Cab Vibration Reduction System and Excavator |
| CN112040351B (en) * | 2020-08-25 | 2022-04-26 | 西安因联信息科技有限公司 | Mechanical equipment vibration alarm information pushing method |
| CN116150895A (en) * | 2021-11-19 | 2023-05-23 | 株洲中车时代电气股份有限公司 | Converter vibration early warning method, device, equipment and medium |
| CN114357567B (en) * | 2021-12-09 | 2024-12-03 | 中国建筑第二工程局有限公司 | A wind vibration monitoring method, system, storage medium and computer equipment based on BIM |
| CN114239661A (en) * | 2021-12-21 | 2022-03-25 | 佳讯飞鸿（北京）智能科技研究院有限公司 | Real-time early warning method and device, electronic equipment and medium |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN206555715U (en) * | 2017-02-27 | 2017-10-13 | 安顺 | A kind of oil pipeline warning device |
| CN107633661A (en) * | 2017-08-28 | 2018-01-26 | 国家电网公司 | Pump-storage generator runout alarm method and device |
| CN109901537A (en) * | 2019-03-18 | 2019-06-18 | 北京大通惠德科技有限公司 | Mechanical equipment method for monitoring operation states and system for edge calculations side |
| CN110648480A (en) * | 2019-10-12 | 2020-01-03 | 山东科技大学 | A rate-of-change-based single-variable alarm system and method |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US10508974B2 (en) * | 2016-02-01 | 2019-12-17 | Computational Systems, Inc. | Storing analytical machine data based on change in scalar machine data indicating alert condition |
-
2020
- 2020-01-14 CN CN202010038460.9A patent/CN111259538B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN206555715U (en) * | 2017-02-27 | 2017-10-13 | 安顺 | A kind of oil pipeline warning device |
| CN107633661A (en) * | 2017-08-28 | 2018-01-26 | 国家电网公司 | Pump-storage generator runout alarm method and device |
| CN109901537A (en) * | 2019-03-18 | 2019-06-18 | 北京大通惠德科技有限公司 | Mechanical equipment method for monitoring operation states and system for edge calculations side |
| CN110648480A (en) * | 2019-10-12 | 2020-01-03 | 山东科技大学 | A rate-of-change-based single-variable alarm system and method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111259538A (en) | 2020-06-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111259538B (en) | Vibration alarm method, storage medium and electronic equipment | |
| CN112348237B (en) | Abnormal trend detection method for dynamic drilling data | |
| CN107066831B (en) | Regional comprehensive environment evaluation method, device and system | |
| US10508974B2 (en) | Storing analytical machine data based on change in scalar machine data indicating alert condition | |
| CN109213654A (en) | A kind of method for detecting abnormality and device | |
| CN116308305B (en) | Bridge health monitoring data management system | |
| CN108535759B (en) | Remote radiation protection monitoring method, device and system for million-kilowatt nuclear power station | |
| CN114446020B (en) | Linkage early warning management method, system, storage medium and equipment | |
| CN106896779A (en) | Maintenance period forecasting system and maintenance period prediction meanss | |
| JP6972868B2 (en) | Safety judgment program, safety judgment method and safety judgment device | |
| CN113006872A (en) | Monitoring device and correction method for key branch air volume of mine ventilation network | |
| CN118051863B (en) | Health data acquisition system and method based on digital metering technology | |
| CN119511394A (en) | Method and system for tunnel stratum monitoring and early warning based on cosmic ray muons | |
| CN108180935A (en) | Fault detection method and device of sensor | |
| CN119509466A (en) | A dynamic environmental sedimentation observation method | |
| CN117978628B (en) | Communication control method and system based on intelligent park | |
| KR101507903B1 (en) | Apparatus and method for estimating observational reference threshold for the event registration of the driving risk index | |
| CN118091294B (en) | Method for detecting resistor | |
| CN113487141A (en) | Bearing state evaluation method based on multi-source information clustering fusion | |
| CN114565330B (en) | Health degree evaluation method, system, equipment and storage medium of precipitation observation equipment | |
| CN110020000B (en) | Method and device for determining abnormal wind speed data | |
| CN117147959A (en) | Electricity larceny prevention analysis method and device | |
| CN111898094B (en) | Method and device for processing foundation pit monitoring data, electronic equipment and storage medium | |
| KR102394464B1 (en) | A SYSTEM AND METHOD FOR CALCULATING THE COMPREHENSIVE CIVIL COMPLAINT INDEX AND CONTROLING THE CRITICAL CIVIL COMPLIAINT POINTS USING IoT TIME SERIES NOISE, VIBRATION DATA AND REGULATORY STANDARDS AND THE CIVIL COMPLAINT INFORMATION | |
| JP6777228B2 (en) | Control devices, communication systems, and control methods |
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 |