Disclosure of Invention
Aiming at the defects in the prior art: due to the fact that information is refreshed caused by massive information alarm and flooding alarm, workload of operation and maintenance personnel is caused, important information is omitted, and information quantity and processing efficiency of system information expression are reduced; the invention provides a compression display method and a compression display system for information redundancy jitter characteristics, and provides huge adjustment brought by a large amount of redundant event information to data persistence storage and offline analysis.
In order to achieve the purpose of the invention, the invention adopts the technical scheme that:
on one hand, the compression display method of the information redundancy jitter characteristic comprises the following steps:
s1, acquiring alarm information jitter characteristic parameters based on the SCADA alarm system in a set acquisition period;
s2, constructing an alarm information jitter rate calculation model according to the alarm information jitter characteristic parameters in the step S1;
s3, judging whether the alarm information jitter rate obtained by the alarm information jitter rate calculation model established in the step S2 exceeds a preset threshold value, if so, entering a step S4, otherwise, normally displaying the alarm information based on the SCADA alarm system;
s4, selecting a time window by utilizing a statistical big data theorem;
and S5, executing and constructing a compression display and decompression strategy based on the time window selected in the step S4, and displaying the alarm information.
Further, step S2 specifically includes the following sub-steps:
s21, calculating the instantaneous frequency value and alarm arrival rate of the alarm information jitter characteristic parameter according to the information jitter characteristic parameter in the sampling period in the step S1;
s22, calculating a jitter rate normalization attenuation index according to the time attenuation factor;
s23, calculating a jitter rate index according to the instantaneous frequency value in the step S21 and the jitter rate normalized attenuation index in the step S22, wherein the jitter rate index is expressed as:
wherein the content of the first and second substances,
for the jitter rate index, the convolution operation, freq (n) is the instantaneous frequency value of the nth sample point in the T period,
normalizing the attenuation index for the jitter rate;
and S24, calculating the jitter rate of the alarm information according to the jitter rate index in the step S23 and the alarm arrival rate in the step S21 to obtain an alarm information jitter rate calculation model.
The further scheme has the following beneficial effects:
the weight of the influence of the alarm reaching each interval on the current is differentiated, and the weight of the interval with long time is smaller than the weight of the interval with short time according to the principle of time attenuation, so that the interval density distribution is convenient to differentiate.
Further, the calculation formula of the alarm arrival rate in step S21 is as follows:
wherein R isTIs the alarm information arrival rate within the interval of T time periods, StimeFor the length of the period of arrival of the alarm information in the T period, ACTThe alarm information arrival number in the T period.
The further scheme has the following beneficial effects:
the main idea of this solution is to take into account not only the influence of the previous time interval, but also the influence of the average rate of arrival of the actual alarm. The influence of local extreme values generated by unbalanced alarm interval distribution is avoided.
Further, step S22 is specifically:
s221, calculating a decay exponent according to the time decay factor, which is expressed as:
wherein D isT,nThe attenuation index of the nth sampling point in the T period is shown, e is a constant, tau is a time index, and N is the total number of alarm sampling points in the SCADA-based system;
s222, calculating a normalization coefficient according to the attenuation exponent in step S221, and expressing as:
wherein A isTIs a normalized coefficient;
s223, calculating a jitter rate normalized attenuation index according to the attenuation index in step S221 and the normalized coefficient in step S222, which is expressed as:
wherein the content of the first and second substances,
the decay index is normalized for jitter rate.
The further scheme has the following beneficial effects:
the normalization used in this scheme is to make the coefficients as weighting coefficients, which need to be normalized as a function of probability density.
Further, the calculation formula of the information jitter rate in step S24 is represented as:
where ACIndex (T) is the information jitter rate, RTIs the alarm arrival rate.
The further scheme has the following beneficial effects:
the scheme integrates two event frequency algorithms of interval density reciprocal (frequency) and alarm arrival rate. And the advantages of two algorithms for reaching the speed can be well balanced by adopting a geometric mean form to synthesize the two algorithms, and the corresponding disadvantages are avoided.
Further, step S5 specifically includes the following sub-steps:
step S5 specifically includes the following substeps:
s51, extracting key features in the alarm information jitter feature parameters in the step S1 in the time window selected in the step S4, and performing information compression according to the key features to construct a group of composite alarm information vector groups;
s52, inquiring whether the composite alarm information vector group in the step S51 has repeated redundant information through a preset alarm information characteristic parameter cache table, if so, combining to generate a new alarm information vector group, and entering the step S53; if not, normally displaying the composite alarm information vector group;
and S53, setting repeated shaking identification state indication bits in the new alarm information vector group generated by combination, and displaying the alarm information.
In another aspect, a system for compressed presentation of information redundancy dithering features includes:
the information acquisition device is used for acquiring the jitter characteristic parameters of the alarm information based on the SCADA alarm system in a set acquisition period;
the alarm information jitter rate processing device is used for calculating the alarm information jitter rate according to the alarm information jitter characteristic parameters, judging whether the alarm information jitter rate exceeds a preset threshold value or not, and carrying out corresponding processing according to the judgment result;
time window selection means for selecting a time window using the statistical big data theorem
And the compression display and decompression strategy device is used for performing compression display and decompression operation on the alarm information based on the SCADA alarm system.
The invention has the following beneficial effects:
collecting alarm information jitter characteristic parameters based on an SCADA alarm system in a set sampling period, analyzing the alarm information jitter characteristic parameters, constructing an alarm information jitter rate calculation model based on the collected information to obtain an alarm information jitter rate, judging whether the alarm information jitter rate exceeds a preset threshold value, if so, selecting a time window by utilizing a statistical big data theorem, executing and constructing a compression display and decompression strategy, and displaying the alarm information, otherwise, directly displaying the alarm information based on the SCADA alarm system; the method has good inhibition effect on flooding alarm and low-efficiency redundant alarm, solves the problem of screen refreshing and screen flickering caused by mass alarm information, improves the efficiency and comfort of a personnel monitoring system from the perspective of human-computer communication engineering, reduces the working strength of operation and maintenance personnel, can fully improve the information expression and storage efficiency of an SCADA alarm log system, improves the loading capacity and the information layout utilization rate of information in the same time-space range, improves the probability distribution of the alarm information to be more balanced, namely meets the requirement of the maximum entropy, and maximizes the alarm information amount (information entropy).
Detailed Description
The following description of the embodiments of the present invention is provided to facilitate the understanding of the present invention by those skilled in the art, but it should be understood that the present invention is not limited to the scope of the embodiments, and it will be apparent to those skilled in the art that various changes may be made without departing from the spirit and scope of the invention as defined and defined in the appended claims, and all matters produced by the invention using the inventive concept are protected.
In order to solve the problems:
1. the information is 'refreshed' when the alarm is given for massive information and the alarm is given for flooding, so that operation and maintenance personnel cannot authenticate and carefully check specific information, and thus important information is omitted;
2. certain 'ptering is not stopped' event information does not stop refreshing the screen, so that important information can quickly move out of the sight of personnel, and the information quantity and efficiency of important information expression of the system are reduced;
3. the massive redundant event information also brings great challenges to data persistence storage and offline analysis and search work;
the invention provides a method and a system for compressing and displaying information redundancy jitter characteristics.
As shown in fig. 1, a method for displaying information redundancy jitter characteristics by compression includes the following steps S1-S5:
s1, acquiring alarm information jitter characteristic parameters based on the SCADA alarm system in a set acquisition period;
s2, constructing an alarm information jitter rate calculation model according to the alarm information jitter characteristic parameters in the step S1;
in practice, the alarm information will be full of a large amount of alarm information with individual events jittering and repeating, and according to the information maximum entropy principle: assuming that an alarm event system is composed of N events, because these events have the characteristic of randomness, an event can randomly occur multiple times (redundant repetition) at different time instants, so the information quantity that the event system can express can be described according to the definition of information entropy, and the frequency (probability) of occurrence of different events should be guaranteed to be as same as possible by knowing the maximum entropy of information transmission, so compressing the probability of occurrence of these events or occupying the proportion of the layout, and maximizing the entropy of the alarm information is the main optimization target of the embodiment of the present invention.
In practice, the alarm information jitter rate calculation model should satisfy the following principle:
1. the number (Alarm Count AC) of the arriving Alarm information in unit time is reflected, namely the number of the actually arriving alarms in the predefined unit time;
2. the interval length and the interval distribution density among all information are embodied, namely under the condition that the principle 1 is the same, different alarm arrival intervals are distinguished, and the arrival intervals are denser than sparse intervals;
3. the influence of the natural attenuation factor (decay index DI) on each interval is to distinguish the time positions where the arrival interval density occurs under the same condition of the principle 2, wherein the influence of the interval density occurring later is larger than that occurring earlier, and the influence weight of the same interval is attenuated along with the time.
According to the three principles, firstly, an alarm information jitter rate calculation model needs to be constructed to obtain an alarm information jitter rate, namely the alarm information jitter rate is the arrival rate of repeated alarm information and can also be expressed as the number of the arrival of the repeated alarm information in unit time.
As shown in fig. 2, in this embodiment, step S2 specifically includes the following sub-steps:
s21, calculating the instantaneous frequency value and alarm arrival rate of the alarm information jitter characteristic parameter according to the information jitter characteristic parameter in the sampling period in the step S1;
in this embodiment, the calculation formula of the alarm arrival rate in step S21 is represented as:
wherein R isTFor alarm information arrival rate within a time interval T, StimeFor the length of the period of arrival of the alarm information in the T period, ACTThe alarm information arrival number in the T period.
In practice, the instantaneous frequency value is calculated as:
(unit: number/second), wherein Δ T
nIs the time interval difference, expressed as: delta T
n=t
n+1-t
n n∈[0,N-1]And N is not less than 2, wherein t
n+1Is the acquisition time of the (n + 1) th sampling point, t
nAnd N is the total number of the alarm sampling points in the SCADA system at the acquisition time of the nth sampling point.
S22, calculating a jitter rate normalization attenuation index according to the time attenuation factor;
as shown in fig. 3, in this embodiment, step S22 specifically includes:
s221, calculating a decay exponent according to the time decay factor, which is expressed as:
wherein D is
T,nIs the decay exponent of the nth sampling point in the T period, the decay exponent D
T,nA probability density coefficient and a time constant tau, e is a constant, tau is a time index for representing the decay speed, when the time T-n is tau, the decay speed is equal to the original value
Multiplying, wherein N is the total number of the alarm sampling points in the SCADA-based system;
s222, calculating a normalization coefficient according to the attenuation exponent in step S221, and expressing as:
wherein A isTIs a normalized coefficient;
s223, calculating a jitter rate normalized attenuation index according to the attenuation index in step S221 and the normalized coefficient in step S222, which is expressed as:
wherein the content of the first and second substances,
normalizing the decay index for jitter rate, normalizing the decay index for jitter rate
Is also a probability density function, satisfies
While being a monotonically decreasing sequence, i.e.:
in practice, the amount of the liquid to be used,
is a normalized coefficient having a probability density property.
S23, calculating a jitter rate index from the instantaneous frequency value in step S21 and the normalized attenuation index in step S22, as:
wherein the content of the first and second substances,
is the jitter rate index (unit: one/second), is the convolution operation, freq (n) is the instantaneous frequency value of the nth sampling point in the T period;
in practice, the calculated instantaneous frequency value is convolved with the normalized attenuation index, and then the convolution structure is normalized, which can also be expressed as:
and S24, calculating the jitter rate of the alarm information according to the jitter rate index in the step S23 and the alarm arrival rate in the step S21 to obtain an alarm information jitter rate calculation model.
In this embodiment, the calculation formula of the information jitter rate in step S24 is represented as:
where ACIndex (T) is the information jitter rate, RTIs the alarm arrival rate.
In practice, the geometric mean of the actual rate and the decay interval factor may be taken in the embodiment of the present invention, and the specific operation result is shown in table 2.
Table 2 calculation example 1
According to the calculation results shown in table 2, it can be calculated that:
comprises the following steps:
is 0.21 (one/second)
Alarm arrival rate R
TComprises the following steps:
(one/second)
The alarm arrival rate acidex is calculated as follows:
(one/second)
Wherein the jitter rate index is analyzed according to the influence of different density distributions, as shown in table 3,
TABLE 3 calculation of different density distributions versus jitter rate index
The density distribution pairs of different intervals can be seen from the table
The calculations have a significant effect, and the effect of high density intervals occurring later is greater,
consistent with the desired results defined by the embodiments of the present invention.
S3, judging whether the alarm information jitter rate obtained by the alarm information jitter rate calculation model established in the step S2 exceeds a preset threshold value, if so, entering a step S4, otherwise, normally displaying the alarm information based on the SCADA alarm system;
in practice, the alarm information jitter rate is detected in real time based on the SCADA alarm system, when the alarm information jitter rate meets a preset threshold value, a corresponding compression display and decompression strategy is executed, otherwise, the processing is carried out according to the original processing mode.
S4, selecting a time window by utilizing a statistical big data theorem;
in practice, the alarm jitter rate according to the ANSI/ISA-18.2-2016 recommendations for alarm systems in the united states should be within a reasonably acceptable range for maximum performance of the SCADA system; it is specified in section 16.4 of the standard that the number of important alarms per operator terminal should meet the requirements of table 4, and exceeding this recommended standard may result in the omission of important alarms or overload of personnel.
TABLE 4 average alarm arrival Rate
According to the alarm arrival rate given by the American ANSI/ISA-18.2-2016, the time range is usually selected to be 10 minutes or 600 seconds, if the alarm frequency of the system is faster or slower, the window size can be adjusted to be larger or smaller according to the actual situation of the system, for example, 1 hour or 1 minute can be selected as the time window; the number of alarms in a larger sliding time window is more usually, the statistical distribution characteristic of data tends to actual data according to the principle of the statistical majority theorem, the influence of accidental rules is reduced, and meanwhile, N-1/N tends to 1 under the condition that the number of N is certain large because the interval of N alarms is N-1, so that the R calculated at intervals can be made to tend to 1T、ΦTAnd the statistical mean arrival rate converges to obtain an appropriate time window.
And S5, executing and constructing a compression display and decompression strategy based on the time window selected in the step S4, and displaying the alarm information.
As shown in fig. 4, in this embodiment, step S5 specifically includes the following sub-steps:
s51, extracting key features in the alarm information jitter feature parameters in the step S1 in the time window selected in the step S4, and performing information compression according to the key features to construct a group of composite alarm information vector groups;
in practice, the alarm event information is usually an information vector described by using a multidimensional vector, and generally consists of an event ID (number), a timestamp, key information (key) of an event, detailed information description, a state, a level, and other attributes of the event, where the key information is used as a unique and necessary feature of point index information for the alarm point, and information compression is performed on a redundant event in the alarm information through the key information to obtain a set of composite alarm information vector groups, where the redundant event has the following features:
1. the event key information is the same, and is usually an event occurring at one information point;
2. the same or related state characteristics, such as action and reset of one point and out-of-limit reset of analog quantity;
3. the high jitter rate generates a large amount of information in a short time, and the information is basically the same and repeated except the state information such as a few times;
it can also be understood that: the idea of 'removing the same and different things' is adopted, the same repeated alarms are 'aggregated' into a piece of composite alarm information by a similar method of extracting 'a formula', the composite number of alarm events can be embodied, meanwhile, the difference part is compressed and displayed in a dimensionality reduction mode, and the repeated information is eliminated to the maximum extent.
S52, inquiring whether the composite alarm information vector group in the step S51 has repeated redundant information through a preset alarm information characteristic parameter cache table, if so, combining to generate a new alarm information vector group, and entering the step S53, otherwise, normally displaying the composite alarm information vector group;
in practice, an alarm information characteristic parameter cache table is constructed by an SCADA (supervisory control and data acquisition) based alarm system according to the existing alarm information characteristic parameters, for example, if the current alarm is N, the latest repeated N-1 information is found, then the difference of the two pieces of information is compared, as shown in part A in FIG. 5, the common part of the two pieces of information is proposed and merged, then a vector table composed of difference items is generated to describe different parts, as shown in part B in FIG. 5, a new alarm information vector group is generated by merging, the difference information described by the new alarm information vector group is the content of the latest arriving alarm information, therefore, whenever the inquiry is carried out by presetting the alarm information characteristic parameter cache table, the merged repeated redundant information into the latest alarm vector group is displayed as the latest alarm information in the SCADA based alarm system, and if no redundant information exists, the normal display is directly performed.
And S53, setting repeated shaking identification state indication bits in the new alarm information vector group generated by combination, and displaying the alarm information.
In practice, as shown in part C in fig. 5, repeated jitter flag state indication bits are set in a new alarm information vector group and generated to distinguish from ordinary alarm information, and when all compressed information of an alarm needs to be checked (decompressed information), the alarm information of different parts can be restored and displayed by performing double-click operation on the alarm information with the repeated jitter flag state indication bits in the SCADA alarm system.
In this embodiment, the following experiment is performed to illustrate a method, as shown in table 5, for a set of possible arriving alarm information, the repeated jitter alarm information of the same key information key is selected here, the information is an information vector described by dimensions such as < serial number, time, key information key, information, status >, and the like, where the same point has multiple alarm vectors at different times and in different statuses.
Table 5 display layout information before compression (total 3 records)
However, since the key information key, description, etc. are the same in different vectors due to different events generated by one information key point, the common part can be extracted as the grouped group header content, and then the information of different sequence numbers, events, levels, states, etc. is repackaged into a new group of vectors, and the group of information vector displays the layout occupying one original alarm display, as shown in table 6, and at the same time, a flag bit for compression and de-duplication is added in the merged information vector to display, the letter R represents that there is duplicated information in the information vector, the number in parentheses is the number of duplicated information, here is 2, which represents that the information contains 2 pieces of the same type information in addition to the information, and the state information is the last recorded state.
TABLE 6 compressed display layout information (reduced to 1 record)
The new similar alarm information is always displayed on the outermost (new) side of the vector group, so that an operator can know that not only the newly generated alarm information exists, but also similar redundant repeated alarm data is contained, the original vector dimension is reduced from 6 dimensions to 3 dimensions, the dimension reduction compression processing of the alarm data from high dimension to low dimension is carried out, the redundant information is eliminated, as shown in table 7, the new vector group can be seen to be organized into a tree structure, the tree root is the public (same) information of the information, and the tree leaf node is a new dimension reduction vector group.
TABLE 7 Structure of compressed display layout information
As shown in fig. 6, a compressed display system based on the redundant jitter characteristics of SCADA alarm system includes:
the information acquisition device is used for acquiring the jitter characteristic parameters of the alarm information based on the SCADA alarm system in a set acquisition period;
the alarm information jitter rate processing device is used for calculating the alarm information jitter rate according to the alarm information jitter characteristic parameters, judging whether the alarm information jitter rate exceeds a preset threshold value or not, and carrying out corresponding processing according to the judgment result;
time window selection means for selecting a time window using the statistical big data theorem
And the compression display and decompression strategy device is used for performing compression display and decompression operation on the alarm information based on the SCADA alarm system.
The compression display system for the information redundancy jitter characteristics provided by the embodiment of the invention has the beneficial effects of the compression display method for the information redundancy jitter characteristics.
The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
The principle and the implementation mode of the invention are explained by applying specific embodiments in the invention, and the description of the embodiments is only used for helping to understand the method and the core idea of the invention; meanwhile, for a person skilled in the art, according to the idea of the present invention, there may be variations in the specific embodiments and the application scope, and in summary, the content of the present specification should not be construed as a limitation to the present invention.
It will be appreciated by those of ordinary skill in the art that the embodiments described herein are intended to assist the reader in understanding the principles of the invention and are to be construed as being without limitation to such specifically recited embodiments and examples. Those skilled in the art, having the benefit of this disclosure, may effect numerous modifications thereto and changes may be made without departing from the scope of the invention in its aspects.