CN115329565A - Comprehensive evaluation method and system for complex electromagnetic field environment - Google Patents

Comprehensive evaluation method and system for complex electromagnetic field environment Download PDF

Info

Publication number
CN115329565A
CN115329565A CN202210949666.6A CN202210949666A CN115329565A CN 115329565 A CN115329565 A CN 115329565A CN 202210949666 A CN202210949666 A CN 202210949666A CN 115329565 A CN115329565 A CN 115329565A
Authority
CN
China
Prior art keywords
electromagnetic field
complex
electric energy
energy meter
field signal
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202210949666.6A
Other languages
Chinese (zh)
Other versions
CN115329565B (en
Inventor
王聪
张蓬鹤
赵兵
陈昊
祝毛宁
杨艺宁
宋如楠
杨柳
赵鹏
姜洪浪
王晓东
王爽
赵婷
张玉冠
左嘉
郭清营
姬云涛
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Electric Power Research Institute Co Ltd CEPRI
Original Assignee
China Electric Power Research Institute Co Ltd CEPRI
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China Electric Power Research Institute Co Ltd CEPRI filed Critical China Electric Power Research Institute Co Ltd CEPRI
Priority to CN202210949666.6A priority Critical patent/CN115329565B/en
Publication of CN115329565A publication Critical patent/CN115329565A/en
Application granted granted Critical
Publication of CN115329565B publication Critical patent/CN115329565B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R35/00Testing or calibrating of apparatus covered by the other groups of this subclass
    • G01R35/04Testing or calibrating of apparatus covered by the other groups of this subclass of instruments for measuring time integral of power or current

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Testing Electric Properties And Detecting Electric Faults (AREA)

Abstract

The invention discloses a comprehensive evaluation method and a comprehensive evaluation system for a complex electromagnetic field environment. Wherein, the method comprises the following steps: constructing complex known electromagnetic field signal characteristics based on various known electromagnetic field signals, and constructing known electric energy meter electricity utilization information characteristics and known electric energy meter working state characteristics corresponding to the known electromagnetic field signals; determining a known coupling relation, and constructing a typical characteristic database according to the known coupling relation; when the interference of the complex electromagnetic field environment occurs, constructing complex electromagnetic field signal characteristics based on complex electromagnetic field signals, and constructing complex electric energy meter electricity utilization information characteristics and complex electric energy meter working state characteristics; determining a complex coupling relation; and comparing and screening the complex coupling relation with a typical characteristic database, determining which electromagnetic field interference signal or signals the complex electromagnetic field environment is or, and recording early warning information.

Description

Comprehensive evaluation method and system for complex electromagnetic field environment
Technical Field
The invention relates to the field of complex electromagnetic field monitoring, in particular to a comprehensive evaluation method and a comprehensive evaluation system for complex electromagnetic field environment.
Background
In recent years, the construction of smart grids is accelerated in China, and smart electric energy meters are widely popularized among users. Because the intelligent electric energy meter is easily subjected to electromagnetic interference in the operation process, the electric energy meter needs to be subjected to a performance test for resisting magnetic field interference before leaving a factory, so that the electric energy meter can still normally work in a constant magnetic field environment of 300mT, and the electric energy meter also has a constant magnetic field interference event record in the operation process.
However, the existing constant magnetic field interference cannot distinguish other complex electromagnetic field interference, the recording function of the constant magnetic field interference event is misreported after other electromagnetic field interference occurs, and the evidence obtaining and evidence fixing effects cannot be achieved when some lawless persons use strong electromagnetic fields to perform interference electricity stealing on the electric energy meter.
Disclosure of Invention
According to the invention, a comprehensive evaluation method and a comprehensive evaluation system for a complex electromagnetic field environment are provided, so as to solve the technical problems that the existing constant magnetic field interference cannot distinguish other complex electromagnetic field interference, the recording function of the constant magnetic field interference event is misreported after other electromagnetic field interference occurs, and the evidence obtaining and evidence fixing effects cannot be achieved when some lawbreakers use strong electromagnetic fields to perform interference electricity stealing on an electric energy meter.
According to a first aspect of the present invention, there is provided a comprehensive evaluation method for a complex electromagnetic field environment, comprising:
constructing known electromagnetic field signal characteristics based on various known electromagnetic field signals, and constructing known electric energy meter electricity utilization information characteristics and known electric energy meter working state characteristics corresponding to the known electromagnetic field signals;
determining the known coupling relation of the known electromagnetic field signal characteristics, the known electric energy meter electricity utilization information characteristics and the known electric energy meter working state characteristics, and constructing a typical characteristic database according to the known coupling relation;
when the interference of the complex electromagnetic field environment occurs, constructing complex electromagnetic field signal characteristics based on complex electromagnetic field signals, and constructing complex electric energy meter electricity utilization information characteristics and complex electric energy meter working state characteristics corresponding to the complex electromagnetic field signals;
determining the complex coupling relation of the complex electromagnetic field signal characteristics, the complex electric energy meter electricity utilization information characteristics and the complex electric energy meter working state characteristics;
and comparing and screening the complex coupling relation with the typical characteristic database, determining which electromagnetic field interference signal or signals the complex electromagnetic field environment occurs, and recording early warning information.
Optionally, the known electromagnetic field signal features are constructed based on various known electromagnetic field signals, including:
decomposing the known electromagnetic field signal from a time domain, a frequency domain and an energy domain to obtain the information of the time domain, the frequency domain and the energy domain of the known electromagnetic field signal;
and obtaining various different known electromagnetic field signal characteristics based on the time domain, frequency domain and energy domain information of the known electromagnetic field signals.
Optionally, constructing a known electric energy meter electricity consumption information characteristic and a known electric energy meter working state characteristic corresponding to the known electromagnetic field signal comprises:
analyzing the waveform difference characteristics of adjacent periods of current, voltage, daily freezing and electricity consumption information of the electric energy meter under various interference states of known electromagnetic field signals to obtain electricity consumption information characteristics of the known electric energy meter;
and extracting the metering error and the fault information of the electric energy meter to obtain the working state characteristics of the known electric energy meter.
Optionally, when the interference of the complex electromagnetic field environment occurs, constructing a complex electromagnetic field signal characteristic based on the complex electromagnetic field signal, and constructing a complex electric energy meter electricity consumption information characteristic and a complex electric energy meter working state characteristic corresponding to the complex electromagnetic field signal, including:
when the environment interference of the complex electromagnetic field occurs, decomposing the complex electromagnetic field signal from a time domain, a frequency domain and an energy domain to obtain the information of the time domain, the frequency domain and the energy domain of the complex electromagnetic field signal, and obtaining the characteristics of the complex electromagnetic field signal based on the information of the time domain, the frequency domain and the energy domain of the complex electromagnetic field signal;
and analyzing the waveform difference characteristics of adjacent periods of the current, voltage, daily freezing and electricity consumption information of the electric energy meter under the complex electromagnetic field signal interference state to obtain the electricity consumption information characteristics of the complex electric energy meter, and extracting the metering error and fault information of the electric energy meter to obtain the working state characteristics of the complex electric energy meter.
Optionally, when the complex electromagnetic field environment interference occurs, decomposing the complex electromagnetic field signal from a time domain, a frequency domain and an energy domain to obtain information of the time domain, the frequency domain and the energy domain of the complex electromagnetic field signal, and obtaining the complex electromagnetic field signal characteristics based on the information of the time domain, the frequency domain and the energy domain of the complex electromagnetic field signal, including:
decomposing the data into a group of dynamic modes, setting the number of spatial measurements of each time snapshot as N, the number of time snapshots as M, arranging the time-series data into N multiplied by M snapshots:
Figure 787315DEST_PATH_IMAGE001
in the formula, ti is the starting time of the ith segment, and the vector x is an N-dimensional column vector;
intercepting a1 xM-1 partial matrix in X, and performing singular value decomposition on X1M-1:
Figure 301473DEST_PATH_IMAGE002
wherein, K represents the rank of singular value decomposition of X1M-1;
calculate the K × K projection of a:
Figure 227841DEST_PATH_IMAGE003
Figure 196540DEST_PATH_IMAGE004
decomposing the characteristics of A1:
Figure 815740DEST_PATH_IMAGE005
in the formula, the column of W is a eigenvector, and a diagonal matrix containing corresponding eigenvalue lambada k is 581;
the approximate decomposition for all future times can be expressed as:
Figure 817195DEST_PATH_IMAGE006
where ζ is a spatial coordinate, and wk = ln (λ k)/Δ t initial coefficient value bk (0) is obtained by a pseudo-inverse (equivalent to least square) method:
the first round of decomposition is as follows:
Figure DEST_PATH_IMAGE007
where ψ k (1) denotes the first modality calculated from a complete M snapshots, and XM/2 is defined as:
Figure 422619DEST_PATH_IMAGE008
by performing a second round of decomposition on the second half of the above equation, XM/2 is decomposed into XM/2= X (1) M/2+ X (2) M/2, iterating until the multi-resolution decomposition is completed.
Optionally, the current, voltage, daily freezing and electric quantity power consumption information in the electric energy meter under the complex electromagnetic field signal interference state is subjected to adjacent period waveform difference characteristic analysis to obtain a complex electric energy meter power consumption information characteristic, and an electric meter metering error and fault information are extracted to obtain a complex electric energy meter working state characteristic, including:
calculating the difference value of the characteristic quantities of the adjacent periods, taking the absolute value of the difference value, and taking the absolute value average value as the measurement of the difference of the waveform of the characteristic quantities of the two adjacent periods;
and carrying out adjacent period waveform difference characteristic analysis on the current in the electric energy meter according to the following formula:
setting the sampling point number of the current in each period as N, the period number N, and the sampling values of the current in two adjacent periods as ik-1, ik:
Figure 888236DEST_PATH_IMAGE009
the average of the absolute values of the current differences Δ ik (j) over a period is:
Figure 881600DEST_PATH_IMAGE010
δ k reflects the magnitude of current difference between two adjacent periods, and is normalized to avoid the influence of current values on current deviation:
Figure 370350DEST_PATH_IMAGE011
the waveform similarity is:
Figure 28733DEST_PATH_IMAGE012
optionally, determining the complex coupling relationship of the complex electromagnetic field signal characteristic, the complex electric energy meter electricity consumption information characteristic and the complex electric energy meter working state characteristic includes:
and determining a correlation coefficient calculation formula of the coupling relation of the complex electromagnetic field signal characteristics, the complex electric energy meter electricity consumption information characteristics and the complex electric energy meter working state characteristics as follows:
Figure 286539DEST_PATH_IMAGE013
in the formula, n is sample energy, ρ is a correlation coefficient, and x and y are corresponding elements in two variables.
According to another aspect of the present invention, there is also provided a comprehensive evaluation system for a complex electromagnetic field environment, comprising:
the known characteristic constructing module is used for constructing known electromagnetic field signal characteristics based on various known electromagnetic field signals, and constructing known electric energy meter electricity utilization information characteristics and known electric energy meter working state characteristics corresponding to the known electromagnetic field signals;
the typical database building module is used for determining the known coupling relations of the known electromagnetic field signal characteristics, the known electric energy meter electricity utilization information characteristics and the known electric energy meter working state characteristics, and building a typical characteristic database according to the known coupling relations;
the complex characteristic constructing module is used for constructing complex electromagnetic field signal characteristics based on complex electromagnetic field signals when complex electromagnetic field environment interference occurs, and constructing complex electric energy meter electricity utilization information characteristics and complex electric energy meter working state characteristics corresponding to the complex electromagnetic field signals;
the complex coupling relation determining module is used for determining complex coupling relations among the complex electromagnetic field signal characteristics, the complex electric energy meter electricity utilization information characteristics and the complex electric energy meter working state characteristics;
and the interference signal determining module is used for comparing and screening the complex coupling relation with the typical characteristic database, determining which electromagnetic field interference signal or electromagnetic field interference signals occur in the complex electromagnetic field environment, and recording early warning information.
Optionally, constructing a known feature module comprises:
the known electromagnetic field signal domain information obtaining submodule is used for decomposing the known electromagnetic field signals from a time domain, a frequency domain and an energy domain to obtain the information of the time domain, the frequency domain and the energy domain of the known electromagnetic field signals;
and the known electromagnetic field signal characteristic submodule is used for obtaining various different known electromagnetic field signal characteristics based on the time domain, the frequency domain and the energy domain information of the known electromagnetic field signals.
Optionally, constructing a known feature module comprises:
the sub-module is used for carrying out adjacent period waveform difference characteristic analysis on current, voltage, daily freezing and electricity consumption information of the electric energy meter under various different interference states of known electromagnetic field signals to obtain electricity consumption information characteristics of the known electric energy meter;
and the sub-module for obtaining the working state characteristics of the known electric energy meter is used for extracting the metering error and the fault information of the electric energy meter to obtain the working state characteristics of the known electric energy meter.
Optionally, constructing a typical database module comprising:
the obtained complex electromagnetic field signal characteristic submodule is used for decomposing the complex electromagnetic field signal from a time domain, a frequency domain and an energy domain when the interference of the complex electromagnetic field environment occurs, obtaining the information of the time domain, the frequency domain and the energy domain of the complex electromagnetic field signal, and obtaining the complex electromagnetic field signal characteristic based on the information of the time domain, the frequency domain and the energy domain of the complex electromagnetic field signal;
and the obtained complex electric energy meter characteristic submodule is used for carrying out adjacent period waveform difference characteristic analysis on current, voltage, daily freezing and electric quantity power consumption information in the electric energy meter under the complex electromagnetic field signal interference state to obtain complex electric energy meter power consumption information characteristics, and extracting electric meter metering errors and fault information to obtain complex electric energy meter working state characteristics.
Optionally, obtaining a complex electromagnetic field signal feature submodule includes:
obtaining a complex electromagnetic field signal characteristic unit, which is used for decomposing data into a group of dynamic modes, setting the number of spatial measurements of each time snapshot as N, the number of time snapshots as M, and arranging time sequence data into N multiplied by M snapshots:
Figure 716383DEST_PATH_IMAGE014
in the formula, ti is the starting time of the ith section, and the vector x is an N-dimensional column vector;
intercepting a1 xM-1 partial matrix in X, and performing singular value decomposition on X1M-1:
Figure 754746DEST_PATH_IMAGE015
in the formula, K represents the rank of singular value decomposition of X1M-1;
calculate the K × K projection of a:
Figure 764291DEST_PATH_IMAGE016
Figure 79866DEST_PATH_IMAGE017
decomposing the A1 characteristics:
Figure 415032DEST_PATH_IMAGE018
in the formula, the column of W is a eigenvector, and \581isa diagonal matrix containing corresponding eigenvalues λ k;
the approximate decomposition for all future times can be expressed as:
Figure 878374DEST_PATH_IMAGE019
where ζ is a spatial coordinate, and wk = ln (λ k)/Δ t initial coefficient value bk (0) is obtained by a pseudo-inverse (equivalent to least square) method:
the first round of decomposition is as follows:
Figure 753927DEST_PATH_IMAGE020
where ψ k (1) denotes the first modality calculated from a complete M snapshots, and XM/2 is defined as:
Figure DEST_PATH_IMAGE021
by performing a second round decomposition on the second half of the above equation, XM/2 is decomposed into XM/2= X (1) M/2+ X (2) M/2, iterating until the multi-resolution decomposition is completed.
Optionally, obtaining a complex electric energy meter feature submodule, including:
the characteristic quantity waveform difference determining unit is used for solving the difference of the characteristic quantities of the adjacent periods, taking the absolute value of the difference, and taking the average value of the absolute values as the measurement of the difference of the characteristic quantity waveforms of the adjacent two periods;
the current waveform analyzing unit is used for analyzing the waveform difference characteristics of adjacent cycles of the current in the electric energy meter according to the following formula:
and setting the number of sampling points of current in each period as N, the number of periods N, and the sampling values of current in two adjacent periods as ik-1, ik:
Figure 986325DEST_PATH_IMAGE022
the average of the absolute values of the current differences Δ ik (j) over a period is:
Figure 180808DEST_PATH_IMAGE023
δ k reflects the magnitude of current difference between two adjacent periods, and is normalized to avoid the influence of current values on current deviation:
Figure 131446DEST_PATH_IMAGE024
the waveform similarity is as follows:
Figure 748373DEST_PATH_IMAGE025
optionally, the determining a complex coupling relationship module includes:
and the correlation coefficient calculation formula of the coupling relation of the complex coupling relation used for determining the complex electromagnetic field signal characteristics, the complex electric energy meter electricity utilization information characteristics and the complex electric energy meter working state characteristics is as follows:
Figure 897594DEST_PATH_IMAGE026
in the formula, n is sample energy, ρ is a correlation coefficient, and x and y are corresponding elements in two variables.
Therefore, the signal characteristics of the complex electromagnetic environment are extracted rapidly and comprehensively, the signal characteristics are analyzed, the categories of the complex electromagnetic field environment are divided according to the characteristics, useful interference events are screened and reported, the problem of false alarm of the events is solved, the effectiveness of event recording is ensured, evidence can be solidified in time when artificial abnormal electromagnetic field impact electricity stealing behaviors occur, and the adverse influence of magnetic field interference on the metering performance of the electric energy meter is reduced or avoided.
Drawings
A more complete understanding of exemplary embodiments of the present invention may be had by reference to the following drawings in which:
fig. 1 is a schematic flow chart of a comprehensive evaluation method for a complex electromagnetic field environment according to the embodiment;
FIG. 2 is a schematic diagram of the coupling relationship between the features according to the present embodiment;
FIG. 3 is a schematic diagram of the comprehensive diagnosis, screening and evaluation of complex electromagnetic field environment by analyzing three types of characteristic data according to the embodiment;
fig. 4 is a schematic diagram of a comprehensive evaluation system for a complex electromagnetic field environment according to the embodiment.
Detailed Description
The exemplary embodiments of the present invention will now be described with reference to the accompanying drawings, however, the present invention may be embodied in many different forms and is not limited to the embodiments described herein, which are provided for complete and complete disclosure of the present invention and to fully convey the scope of the present invention to those skilled in the art. The terms used in the exemplary embodiments shown in the drawings are not intended to limit the present invention. In the drawings, the same unit/element is denoted by the same reference numeral.
Unless otherwise defined, terms (including technical and scientific terms) used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs. Further, it will be understood that terms, such as those defined in commonly used dictionaries, should be interpreted as having a meaning that is consistent with their meaning in the context of the relevant art and will not be interpreted in an idealized or overly formal sense.
According to a first aspect of the present invention, there is provided a comprehensive evaluation method 100 for a complex electromagnetic field environment, as shown with reference to fig. 1, the method 100 comprising:
s101, constructing complex known electromagnetic field signal characteristics based on various known electromagnetic field signals, and constructing known electric energy meter electricity utilization information characteristics corresponding to the known electromagnetic field signals, electric energy meter metering data characteristics and known electric energy meter working state characteristics; (ii) a
S102, determining known coupling relations among the known electromagnetic field signal characteristics, the known electric energy meter electricity utilization information characteristics and the known electric energy meter working state characteristics, and constructing a typical characteristic database according to the known coupling relations;
s103, when the interference of the complex electromagnetic field environment occurs, constructing complex electromagnetic field signal characteristics based on complex electromagnetic field signals, and constructing complex electric energy meter electricity utilization information characteristics and complex electric energy meter working state characteristics corresponding to the complex electromagnetic field signals; s104, determining the complex coupling relation among the complex electromagnetic field signal characteristics, the complex electric energy meter electricity utilization information characteristics and the complex electric energy meter working state characteristics;
and S105, comparing and screening the complex coupling relation with the typical characteristic database, determining which electromagnetic field interference signal or signals of the generated complex electromagnetic field environment exist, and recording early warning information.
Specifically, a complex electromagnetic field environment signal is decomposed, the signal is decomposed from a time domain, a frequency domain and an energy domain, the time domain, the frequency domain and the energy domain information of the complex electromagnetic field environment is obtained, and complex electromagnetic field signal characteristics are constructed;
analyzing waveform difference characteristics of adjacent cycles of power consumption information such as current, voltage, daily freezing and electric quantity in the electric energy meter, extracting metering errors, fault information and the like of the electric energy meter, and constructing working state characteristics of the electric energy meter;
extracting different electromagnetic field generation signals such as a constant magnetic field, a power frequency magnetic field, a high-frequency electric field and the like to construct various electromagnetic field signal characteristics.
The coupling relation among the electric energy meter electricity utilization information characteristics, the electric energy meter working state characteristics, the electromagnetic field signal characteristics and the complex electromagnetic field signal characteristics is analyzed, the linear and nonlinear correlation characteristics among the electric energy meter electricity utilization information characteristics, the electric energy meter working state characteristics, the electromagnetic field signal characteristics and the complex electromagnetic field signals are revealed, the redundancy characteristics between the characteristics and the complex electromagnetic field signals are eliminated (refer to fig. 2), and comprehensive diagnosis, screening and evaluation of a complex electromagnetic field environment are carried out by analyzing three types of characteristic data (refer to fig. 3) aiming at three types of characteristic data of electric energy meter electricity utilization information (current, voltage, daily freezing, electric quantity and the like), meter working state (metering error and the like) and typical characteristic quantity of the complex electromagnetic field. The specific method comprises the following steps:
(1) the complex electromagnetic field signals are separated into modal components with different space-time scales according to different resolutions, and the instantaneous frequency and energy characteristics of each mode are collected by utilizing time-frequency distribution to obtain the characteristics of a time domain, a frequency domain and an energy domain.
Decomposing the data into a group of dynamic modes, setting the number of spatial measurements of each time snapshot as N and the number of time snapshots as M, firstly arranging the time series data as N multiplied by M snapshots:
Figure 840142DEST_PATH_IMAGE027
where ti is the start time of the ith segment and the vector x is the N-dimensional column vector.
Then intercepting a1 xM-1 partial matrix in X, and carrying out singular value decomposition on X1M-1:
Figure 12498DEST_PATH_IMAGE028
where, denotes the conjugate transpose, K is the rank of the singular value decomposition of X1M-1.
Calculate the K × K projection of a:
Figure 105219DEST_PATH_IMAGE029
Figure 312209DEST_PATH_IMAGE030
decomposing the A1 characteristics:
Figure 425659DEST_PATH_IMAGE031
in the formula, the column of W is a eigenvector, and \581isa diagonal matrix containing corresponding eigenvalues λ k.
The approximate decomposition for all future times can be expressed as:
Figure 147627DEST_PATH_IMAGE032
where ζ is a spatial coordinate, and wk = ln (λ k)/Δ t initial coefficient value bk (0) is obtained by a pseudo-inverse (equivalent to least square) method.
The first round of decomposition is as follows:
Figure 106356DEST_PATH_IMAGE033
where ψ k (1) indicates that the first modality is calculated from the full M snapshots. XM/2 is defined as:
Figure 167853DEST_PATH_IMAGE034
by performing a second round of decomposition on the second half of the above equation, XM/2 is decomposed into XM/2= X (1) M/2+ X (2) M/2, iterating until the multi-resolution decomposition is completed.
(2) The waveform difference characteristics of adjacent periods of each characteristic quantity of the electric energy meter are analyzed, and the identification of the complex electromagnetic field environment is assisted by analyzing the electricity consumption information such as current, voltage, daily freezing and electric quantity in the electric energy meter.
The difference value of the characteristic quantities of the adjacent periods is firstly solved, then the absolute value of the difference value is obtained, and finally the average value of the absolute value is used for measuring the difference of the characteristic quantity waveforms of the two adjacent periods. The current is taken as an example as follows:
and respectively setting the sampling point number of the current in each period as N, the period number N and the sampling values of the current in two adjacent periods as ik-1,ik.
Figure 373575DEST_PATH_IMAGE035
The average of the absolute values of the current differences Δ ik (j) over a period is:
Figure 520522DEST_PATH_IMAGE036
δ k reflects the magnitude of the current difference between two adjacent periods, and δ k is normalized to avoid the influence of current values on current deviation.
Figure 345259DEST_PATH_IMAGE037
The waveform similarity is as follows:
Figure 526842DEST_PATH_IMAGE038
the change degree of adjacent periods of each characteristic quantity of the electric energy meter can be seen through the waveform similarity, and further the influence of the complex electromagnetic field pulse on the change degree can be reflected.
(3) Even if an abnormal value occurs in the data, because the order of the abnormal value usually does not change obviously, the influence on the correlation coefficient is small, the coupling relation among the electricity consumption information characteristic, the electric energy meter working state characteristic, the electromagnetic field signal characteristic and the complex electromagnetic field signal characteristic of the electric energy meter is obtained, and the correlation coefficient calculation formula of the coupling relation is as follows:
Figure 450935DEST_PATH_IMAGE039
in the formula, n is sample energy, rho is a correlation coefficient, and x and y are corresponding elements in two variables.
By calculating the electricity consumption information characteristics of the electric energy meter, the working state characteristics of the electric energy meter and the correlation coefficient between the electromagnetic field signal characteristics and each complex electromagnetic field signal, the redundant characteristics between each characteristic and the complex electromagnetic field impact scene can be effectively eliminated.
(4) The method comprises the steps of extracting fusion characteristics capable of comprehensively representing three types of data aiming at three types of characteristic data of electric energy meter electricity utilization information (current, voltage, daily freezing, electric quantity and the like), meter working states (metering errors and the like) and typical characteristic quantities of complex electromagnetic fields, constructing a fusion layer of deep characteristics of three types of complex electromagnetic field environments, fully retaining original three types of complex electromagnetic field environment characteristic information by the fusion layer, and carrying out comprehensive diagnosis, screening and evaluation on the complex electromagnetic field environment by analyzing the three types of characteristic data.
Therefore, the signal characteristics of the complex electromagnetic environment are extracted rapidly and comprehensively, the signal characteristics are analyzed, the categories of the complex electromagnetic field environment are divided according to the characteristics, useful interference events are screened and reported, the problem of false alarm of the events is solved, the effectiveness of event recording is ensured, evidence can be solidified in time when artificial abnormal electromagnetic field impact electricity stealing behaviors occur, and the adverse influence of magnetic field interference on the metering performance of the electric energy meter is reduced or avoided.
Optionally, the known electromagnetic field signal characteristics are constructed based on various known electromagnetic field signals, including:
decomposing the known electromagnetic field signal from a time domain, a frequency domain and an energy domain to obtain the information of the time domain, the frequency domain and the energy domain of the known electromagnetic field signal;
and obtaining various different known electromagnetic field signal characteristics based on the time domain, frequency domain and energy domain information of the known electromagnetic field signals.
Optionally, constructing a known electric energy meter electricity information characteristic and a known electric energy meter working state characteristic corresponding to the known electromagnetic field signal comprises:
analyzing the waveform difference characteristics of adjacent periods of current, voltage, daily freezing and electricity consumption information of the electric energy meter under various different interference states of known electromagnetic field signals to obtain the electricity consumption information characteristics of the known electric energy meter;
and extracting the metering error and the fault information of the electric energy meter to obtain the working state characteristics of the known electric energy meter.
Optionally, when the interference of the complex electromagnetic field environment occurs, constructing a complex electromagnetic field signal characteristic based on the complex electromagnetic field signal, and constructing a complex electric energy meter electricity consumption information characteristic and a complex electric energy meter working state characteristic corresponding to the complex electromagnetic field signal, including:
when the environment interference of the complex electromagnetic field occurs, decomposing the complex electromagnetic field signal from a time domain, a frequency domain and an energy domain to obtain the information of the time domain, the frequency domain and the energy domain of the complex electromagnetic field signal, and obtaining the characteristics of the complex electromagnetic field signal based on the information of the time domain, the frequency domain and the energy domain of the complex electromagnetic field signal;
and analyzing the waveform difference characteristics of adjacent periods of the current, voltage, daily freezing and electricity consumption information of the electric energy meter under the complex electromagnetic field signal interference state to obtain the electricity consumption information characteristics of the complex electric energy meter, and extracting the metering error and fault information of the electric energy meter to obtain the working state characteristics of the complex electric energy meter.
Optionally, when the complex electromagnetic field environment interference occurs, decomposing the complex electromagnetic field signal from a time domain, a frequency domain and an energy domain to obtain information of the time domain, the frequency domain and the energy domain of the complex electromagnetic field signal, and obtaining the complex electromagnetic field signal characteristics based on the information of the time domain, the frequency domain and the energy domain of the complex electromagnetic field signal, including:
decomposing data into a group of dynamic modes, setting the number of spatial measurements of each time snapshot as N, the number of time snapshots as M, arranging time sequence data into N × M snapshots:
Figure 22862DEST_PATH_IMAGE040
in the formula, ti is the starting time of the ith section, and the vector x is an N-dimensional column vector;
intercepting a1 xM-1 partial matrix in X, and performing singular value decomposition on X1M-1:
Figure 588973DEST_PATH_IMAGE041
in the formula, K represents the rank of singular value decomposition of X1M-1;
calculate the K × K projection of a:
Figure 625062DEST_PATH_IMAGE042
Figure 782374DEST_PATH_IMAGE043
decomposing the A1 characteristics:
Figure 169493DEST_PATH_IMAGE044
in the formula, the column of W is a eigenvector, and \581isa diagonal matrix containing corresponding eigenvalues λ k;
the approximate decomposition for all future times can be expressed as:
Figure 273715DEST_PATH_IMAGE045
where ζ is a spatial coordinate, and wk = ln (λ k)/Δ t initial coefficient value bk (0) is obtained by a pseudo-inverse (equivalent to least square) method:
the first round of decomposition is as follows:
Figure 590076DEST_PATH_IMAGE020
where ψ k (1) denotes the first modality calculated from a complete M snapshots, and XM/2 is defined as:
Figure 121552DEST_PATH_IMAGE046
by performing a second round of decomposition on the second half of the above equation, XM/2 is decomposed into XM/2= X (1) M/2+ X (2) M/2, iterating until the multi-resolution decomposition is completed.
Optionally, the current, voltage, daily freezing and electric quantity power consumption information in the electric energy meter under the complex electromagnetic field signal interference state is subjected to adjacent period waveform difference characteristic analysis to obtain a complex electric energy meter power consumption information characteristic, and an electric meter metering error and fault information are extracted to obtain a complex electric energy meter working state characteristic, including:
calculating the difference value of the characteristic quantities of the adjacent periods, taking the absolute value of the difference value, and taking the average value of the absolute value as the measurement of the difference of the waveform of the characteristic quantities of the adjacent two periods;
and carrying out adjacent period waveform difference characteristic analysis on the current in the electric energy meter according to the following formula:
and setting the number of sampling points of current in each period as N, the number of periods N, and the sampling values of current in two adjacent periods as ik-1, ik:
Figure 730388DEST_PATH_IMAGE047
the average of the absolute values of the current differences Δ ik (j) over a period is:
Figure 700618DEST_PATH_IMAGE048
δ k reflects the magnitude of the current difference between two adjacent periods, and is normalized to avoid the influence of current values on current deviation:
Figure 711299DEST_PATH_IMAGE049
the waveform similarity is as follows:
Figure 413676DEST_PATH_IMAGE050
optionally, determining the complex coupling relationship of the complex electromagnetic field signal characteristic, the complex electric energy meter power consumption information characteristic and the complex electric energy meter working state characteristic comprises:
and determining a correlation coefficient calculation formula of the coupling relation of the complex electromagnetic field signal characteristics, the complex electric energy meter electricity consumption information characteristics and the complex electric energy meter working state characteristics as follows:
Figure 181912DEST_PATH_IMAGE051
in the formula, n is sample energy, rho is a correlation coefficient, and x and y are corresponding elements in two variables.
Therefore, the signal characteristics of the complex electromagnetic environment are extracted rapidly and comprehensively, the signal characteristics are analyzed, the categories of the complex electromagnetic field environment are divided according to the characteristics, useful interference events are screened and reported, the problem of false alarm of the events is solved, the effectiveness of event recording is ensured, evidence can be solidified in time when artificial abnormal electromagnetic field impact electricity stealing behaviors occur, and the adverse influence of magnetic field interference on the metering performance of the electric energy meter is reduced or avoided.
According to another aspect of the present invention, there is also provided a comprehensive evaluation system 400 for a complex electromagnetic field environment, as shown with reference to fig. 4, the system 400 comprising:
a known characteristic constructing module 410, configured to construct a known electromagnetic field signal characteristic based on various known electromagnetic field signals, and construct a known electric energy meter electricity consumption information characteristic and a known electric energy meter operating state characteristic corresponding to the known electromagnetic field signal;
a typical database constructing module 420, which is used for determining the known coupling relations of the known electromagnetic field signal characteristics, the known electric energy meter electricity utilization information characteristics and the known electric energy meter working state characteristics, and constructing a typical characteristic database according to the known coupling relations;
the complex characteristic constructing module 430 is used for constructing complex electromagnetic field signal characteristics based on complex electromagnetic field signals when complex electromagnetic field environment interference occurs, and constructing complex electric energy meter electricity utilization information characteristics and complex electric energy meter working state characteristics corresponding to the complex electromagnetic field signals;
the complex coupling relation determining module 440 is used for determining the complex coupling relation among the complex electromagnetic field signal characteristics, the complex electric energy meter electricity consumption information characteristics and the complex electric energy meter working state characteristics;
and an interference signal determining module 450, configured to compare and screen the complex coupling relationship with the typical characteristic database, determine which electromagnetic field interference signal or signals the complex electromagnetic field environment is or are, and record early warning information.
Optionally, constructing a known feature module comprises:
the known electromagnetic field signal domain information obtaining submodule is used for decomposing the known electromagnetic field signal from a time domain, a frequency domain and an energy domain to obtain the information of the time domain, the frequency domain and the energy domain of the known electromagnetic field signal;
and the sub-module for obtaining the known electromagnetic field signal characteristics is used for obtaining various different known electromagnetic field signal characteristics based on the time domain, frequency domain and energy domain information of the known electromagnetic field signals.
Optionally, constructing a known feature module comprises:
the sub-module is used for carrying out adjacent period waveform difference characteristic analysis on current, voltage, daily freezing and electricity consumption information of the electric energy meter under various different interference states of known electromagnetic field signals to obtain electricity consumption information characteristics of the known electric energy meter;
and the sub-module for obtaining the working state characteristics of the known electric energy meter is used for extracting the metering error and the fault information of the electric energy meter to obtain the working state characteristics of the known electric energy meter.
Optionally, constructing a typical database module comprising:
the obtained complex electromagnetic field signal characteristic submodule is used for decomposing the complex electromagnetic field signal from a time domain, a frequency domain and an energy domain when the interference of the complex electromagnetic field environment occurs, obtaining the information of the time domain, the frequency domain and the energy domain of the complex electromagnetic field signal, and obtaining the complex electromagnetic field signal characteristic based on the information of the time domain, the frequency domain and the energy domain of the complex electromagnetic field signal;
and the obtained complex electric energy meter characteristic submodule is used for carrying out adjacent period waveform difference characteristic analysis on current, voltage, daily freezing and electric quantity power consumption information in the electric energy meter under the complex electromagnetic field signal interference state to obtain complex electric energy meter power consumption information characteristics, and extracting electric meter metering errors and fault information to obtain complex electric energy meter working state characteristics. Optionally, obtaining a complex electromagnetic field signal feature submodule includes:
obtaining a complex electromagnetic field signal characteristic unit, which is used for decomposing data into a group of dynamic modes, setting the number of spatial measurements of each time snapshot as N, the number of time snapshots as M, and arranging time sequence data into N multiplied by M snapshots:
Figure 893516DEST_PATH_IMAGE052
in the formula, ti is the starting time of the ith section, and the vector x is an N-dimensional column vector;
intercepting a1 xM-1 partial matrix in X, and performing singular value decomposition on X1M-1:
Figure 493124DEST_PATH_IMAGE053
in the formula, K represents the rank of singular value decomposition of X1M-1;
calculate the K × K projection of a:
Figure 428719DEST_PATH_IMAGE054
Figure 12147DEST_PATH_IMAGE055
decomposing the A1 characteristics:
Figure 996284DEST_PATH_IMAGE056
in the formula, the column of W is a eigenvector, and a diagonal matrix containing corresponding eigenvalue lambada k is 581;
the approximate decomposition for all future times can be expressed as:
Figure 902929DEST_PATH_IMAGE057
where ζ is a spatial coordinate, and wk = ln (λ k)/Δ t initial coefficient value bk (0) is obtained by a pseudo-inverse (equivalent to least square) method:
the first round decomposition process is as follows:
Figure 212688DEST_PATH_IMAGE058
where ψ k (1) denotes the first modality calculated from a complete M snapshots, and XM/2 is defined as:
Figure 17833DEST_PATH_IMAGE059
by performing a second round of decomposition on the second half of the above equation, XM/2 is decomposed into XM/2= X (1) M/2+ X (2) M/2, iterating until the multi-resolution decomposition is completed.
Optionally, obtaining a complex electric energy meter feature submodule, including:
the waveform difference determining unit is used for solving the difference of the characteristic quantities of the adjacent periods, taking the absolute value of the difference, and taking the average value of the absolute values as the difference for measuring the waveform of the characteristic quantities of the adjacent two periods;
the current waveform analyzing unit is used for analyzing the waveform difference characteristics of adjacent periods of the current in the electric energy meter according to the following formula:
setting the sampling point number of the current in each period as N, the period number N, and the sampling values of the current in two adjacent periods as ik-1, ik:
Figure 867977DEST_PATH_IMAGE060
the average of the absolute values of the current differences Δ ik (j) over one period is:
Figure 707757DEST_PATH_IMAGE061
δ k reflects the magnitude of current difference between two adjacent periods, and is normalized to avoid the influence of current values on current deviation:
Figure 922838DEST_PATH_IMAGE062
the waveform similarity is as follows:
Figure 887383DEST_PATH_IMAGE063
optionally, the determining a complex coupling relationship module includes:
and the correlation coefficient calculation formula of the coupling relation of the complex coupling relation used for determining the complex electromagnetic field signal characteristics, the complex electric energy meter electricity utilization information characteristics and the complex electric energy meter working state characteristics is as follows:
Figure 744480DEST_PATH_IMAGE064
in the formula, n is sample energy, rho is a correlation coefficient, and x and y are corresponding elements in two variables.
The comprehensive evaluation system 400 for a complex electromagnetic field environment according to an embodiment of the present invention corresponds to the comprehensive evaluation method 100 for a complex electromagnetic field environment according to another embodiment of the present invention, and will not be described herein again.
As will be appreciated by one skilled in the art, embodiments of the present application may be provided as a method, system, or computer program product. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein. The scheme in the embodiment of the application can be implemented by adopting various computer languages, such as object-oriented programming language Java and transliterated scripting language JavaScript.
The present application is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the application. 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.
While the preferred embodiments of the present application have been described, additional variations and modifications in those embodiments may occur to those skilled in the art once they learn of the basic inventive concepts. Therefore, it is intended that the appended claims be interpreted as including the preferred embodiment and all changes and modifications that fall within the scope of the present application.
It will be apparent to those skilled in the art that various changes and modifications may be made in the present application without departing from the spirit and scope of the application. Thus, if such modifications and variations of the present application fall within the scope of the claims of the present application and their equivalents, the present application is intended to include such modifications and variations as well.

Claims (15)

1. A comprehensive evaluation method for a complex electromagnetic field environment is characterized by comprising the following steps:
constructing known electromagnetic field signal characteristics based on various known electromagnetic field signals, and constructing known electric energy meter electricity utilization information characteristics and known electric energy meter working state characteristics corresponding to the known electromagnetic field signals;
determining the known coupling relation of the known electromagnetic field signal characteristics, the known electric energy meter electricity utilization information characteristics and the known electric energy meter working state characteristics, and constructing a typical characteristic database according to the known coupling relation;
when the complex electromagnetic field environment interference occurs, constructing complex electromagnetic field signal characteristics based on complex electromagnetic field signals, and constructing complex electric energy meter electricity utilization information characteristics and complex electric energy meter working state characteristics corresponding to the complex electromagnetic field signals;
determining the complex coupling relation of the complex electromagnetic field signal characteristics, the complex electric energy meter electricity utilization information characteristics and the complex electric energy meter working state characteristics;
and comparing and screening the complex coupling relation with the typical characteristic database, determining which one or more electromagnetic field interference signals the complex electromagnetic field environment generates, and recording early warning information.
2. The method of claim 1, wherein constructing known electromagnetic field signal signatures based on various known electromagnetic field signals comprises:
decomposing the known electromagnetic field signal from a time domain, a frequency domain and an energy domain to obtain the information of the time domain, the frequency domain and the energy domain of the known electromagnetic field signal;
and obtaining various different known electromagnetic field signal characteristics based on the time domain, frequency domain and energy domain information of the known electromagnetic field signals.
3. The method of claim 1, wherein constructing a known power meter power usage information signature and a known power meter operating state signature corresponding to the known electromagnetic field signal comprises:
analyzing the waveform difference characteristics of adjacent periods of current, voltage, daily freezing and electricity consumption information of the electric energy meter under various different interference states of known electromagnetic field signals to obtain the electricity consumption information characteristics of the known electric energy meter;
and extracting the metering error and the fault information of the electric energy meter to obtain the working state characteristics of the known electric energy meter.
4. The method of claim 1, wherein, in the event of a complex electromagnetic field environment interference, constructing a complex electromagnetic field signal characteristic based on a complex electromagnetic field signal, and constructing a complex electric energy meter power consumption information characteristic and a complex electric energy meter operating state characteristic corresponding to the complex electromagnetic field signal, comprises:
when the environment interference of the complex electromagnetic field occurs, decomposing the complex electromagnetic field signal from a time domain, a frequency domain and an energy domain to obtain the information of the time domain, the frequency domain and the energy domain of the complex electromagnetic field signal, and obtaining the characteristics of the complex electromagnetic field signal based on the information of the time domain, the frequency domain and the energy domain of the complex electromagnetic field signal;
and analyzing the waveform difference characteristics of adjacent periods of the current, voltage, daily freezing and electricity consumption information of the electric energy meter under the complex electromagnetic field signal interference state to obtain the electricity consumption information characteristics of the complex electric energy meter, and extracting the metering error and fault information of the electric energy meter to obtain the working state characteristics of the complex electric energy meter.
5. The method of claim 4, wherein when the interference of the complex electromagnetic field environment occurs, decomposing the complex electromagnetic field signal from a time domain, a frequency domain and an energy domain to obtain information of the time domain, the frequency domain and the energy domain of the complex electromagnetic field signal, and obtaining the complex electromagnetic field signal characteristics based on the information of the time domain, the frequency domain and the energy domain of the complex electromagnetic field signal comprises:
decomposing data into a group of dynamic modes, setting the number of spatial measurements of each time snapshot as N, the number of time snapshots as M, arranging time sequence data into N × M snapshots:
Figure 907234DEST_PATH_IMAGE001
in the formula, ti is the starting time of the ith section, and the vector x is an N-dimensional column vector;
intercepting a1 xM-1 partial matrix in X, and performing singular value decomposition on X1M-1:
Figure 280446DEST_PATH_IMAGE002
wherein, K represents the rank of singular value decomposition of X1M-1;
calculate the K × K projection of a:
Figure 675656DEST_PATH_IMAGE003
Figure 755607DEST_PATH_IMAGE004
decomposing the A1 characteristics:
Figure 781332DEST_PATH_IMAGE005
in the formula, the column of W is a eigenvector, and \581isa diagonal matrix containing corresponding eigenvalues λ k;
the approximate decomposition for all future times can be expressed as:
Figure 313944DEST_PATH_IMAGE006
where ζ is a spatial coordinate, and wk = ln (λ k)/Δ t initial coefficient value bk (0) is obtained by a pseudo-inverse method:
the first round of decomposition is as follows:
Figure 591473DEST_PATH_IMAGE007
where ψ k (1) denotes the first modality calculated from a complete M snapshots, and XM/2 is defined as:
Figure 525931DEST_PATH_IMAGE008
by performing a second round of decomposition on the second half of the above equation, XM/2 is decomposed into XM/2= X (1) M/2+ X (2) M/2, iterating until the multi-resolution decomposition is completed.
6. The method of claim 4, wherein the adjacent period waveform difference characteristic analysis is performed on the current, voltage, daily freezing and electricity consumption information of the electric energy meter in the complex electromagnetic field signal interference state to obtain the electricity consumption information characteristic of the complex electric energy meter, and the meter metering error and fault information are extracted to obtain the working state characteristic of the complex electric energy meter, and the method comprises the following steps:
calculating the difference value of the characteristic quantities of the adjacent periods, taking the absolute value of the difference value, and taking the absolute value average value as the measurement of the difference of the waveform of the characteristic quantities of the two adjacent periods;
and carrying out adjacent period waveform difference characteristic analysis on the current in the electric energy meter according to the following formula:
setting the sampling point number of the current in each period as N, the period number N, and the sampling values of the current in two adjacent periods as ik-1, ik:
Figure 722557DEST_PATH_IMAGE009
the average of the absolute values of the current differences Δ ik (j) over a period is:
Figure 476887DEST_PATH_IMAGE010
δ k reflects the magnitude of current difference between two adjacent periods, and is normalized to avoid the influence of current values on current deviation:
Figure 807374DEST_PATH_IMAGE011
the waveform similarity is as follows:
Figure 596338DEST_PATH_IMAGE012
7. the method of claim 1, wherein determining the complex coupling relationship of the complex electromagnetic field signal characteristic, the complex electric energy meter power consumption information characteristic and the complex electric energy meter operating state characteristic comprises:
and determining a correlation coefficient calculation formula of the complex coupling relation of the complex electromagnetic field signal characteristic, the complex electric energy meter power utilization information characteristic and the complex electric energy meter working state characteristic as follows:
Figure 26183DEST_PATH_IMAGE013
in the formula, n is sample energy, ρ is a correlation coefficient, and x and y are corresponding elements in two variables.
8. A comprehensive evaluation system for complex electromagnetic field environments, comprising:
the known characteristic building module is used for building known electromagnetic field signal characteristics based on various known electromagnetic field signals and building known electric energy meter electricity utilization information characteristics and known electric energy meter working state characteristics corresponding to the known electromagnetic field signals;
the typical database building module is used for determining the known coupling relations of the known electromagnetic field signal characteristics, the known electric energy meter electricity utilization information characteristics and the known electric energy meter working state characteristics, and building a typical characteristic database according to the known coupling relations;
the complex characteristic constructing module is used for constructing complex electromagnetic field signal characteristics based on complex electromagnetic field signals when complex electromagnetic field environment interference occurs, and constructing complex electric energy meter electricity utilization information characteristics and complex electric energy meter working state characteristics corresponding to the complex electromagnetic field signals;
the complex coupling relation determining module is used for determining complex coupling relations among the complex electromagnetic field signal characteristics, the complex electric energy meter electricity utilization information characteristics and the complex electric energy meter working state characteristics;
and the interference signal determining module is used for comparing and screening the complex coupling relation with the typical characteristic database, determining which electromagnetic field interference signal or electromagnetic field interference signals occur in the complex electromagnetic field environment, and recording early warning information.
9. The system of claim 8, wherein constructing a known feature module comprises:
the known electromagnetic field signal domain information obtaining submodule is used for decomposing the known electromagnetic field signal from a time domain, a frequency domain and an energy domain to obtain the information of the time domain, the frequency domain and the energy domain of the known electromagnetic field signal;
and the known electromagnetic field signal characteristic submodule is used for obtaining various different known electromagnetic field signal characteristics based on the time domain, the frequency domain and the energy domain information of the known electromagnetic field signals.
10. The system of claim 8, wherein constructing a known feature module comprises:
the sub-module is used for carrying out adjacent period waveform difference characteristic analysis on current, voltage, daily freezing and electricity consumption information of the electric energy meter under various different interference states of known electromagnetic field signals to obtain electricity consumption information characteristics of the known electric energy meter;
and the submodule for obtaining the working state characteristics of the known electric energy meter is used for extracting the metering error and the fault information of the electric energy meter to obtain the working state characteristics of the known electric energy meter.
11. The system of claim 8, wherein building a typical database module comprises:
the obtained complex electromagnetic field signal characteristic submodule is used for decomposing the complex electromagnetic field signal from a time domain, a frequency domain and an energy domain when the interference of the complex electromagnetic field environment occurs, obtaining the information of the time domain, the frequency domain and the energy domain of the complex electromagnetic field signal, and obtaining the complex electromagnetic field signal characteristic based on the information of the time domain, the frequency domain and the energy domain of the complex electromagnetic field signal;
and the obtained characteristic submodule of the complex electric energy meter is used for analyzing the waveform difference characteristics of adjacent periods of current, voltage, daily freezing and electricity consumption information in the electric energy meter under the state of the interference of the complex electromagnetic field signals to obtain the electricity consumption information characteristics of the complex electric energy meter, and extracting the metering error and fault information of the electric energy meter to obtain the working state characteristics of the complex electric energy meter.
12. The system of claim 11, wherein the derive complex electromagnetic field signal characterization submodule comprises:
and obtaining a complex electromagnetic field signal characteristic unit, decomposing the data into a group of dynamic modes, setting the number of spatial measurements of each time snapshot to be N, setting the number of time snapshots to be M, and arranging the time sequence data into N multiplied by M snapshots:
Figure 471070DEST_PATH_IMAGE014
in the formula, ti is the starting time of the ith segment, and the vector x is an N-dimensional column vector;
intercepting a1 xM-1 partial matrix in X, and performing singular value decomposition on X1M-1:
Figure 11773DEST_PATH_IMAGE015
wherein, K represents the rank of singular value decomposition of X1M-1;
calculate the K × K projection of a:
Figure 733873DEST_PATH_IMAGE016
Figure 334618DEST_PATH_IMAGE017
decomposing the characteristics of A1:
Figure 63540DEST_PATH_IMAGE018
in the formula, the column of W is a eigenvector, and \581isa diagonal matrix containing corresponding eigenvalues λ k;
the approximate decomposition for all future times can be expressed as:
Figure 345617DEST_PATH_IMAGE019
where ζ is a spatial coordinate, and wk = ln (λ k)/Δ t initial coefficient value bk (0) is obtained by a pseudo-inverse method:
the first round of decomposition is as follows:
Figure 109173DEST_PATH_IMAGE020
where ψ k (1) denotes the first modality calculated from a complete M snapshots, and XM/2 is defined as:
Figure 208716DEST_PATH_IMAGE021
by performing a second round decomposition on the second half of the above equation, XM/2 is decomposed into XM/2= X (1) M/2+ X (2) M/2, iterating until the multi-resolution decomposition is completed.
13. The system of claim 11, wherein deriving a complex electric energy meter characterization sub-module comprises:
the characteristic quantity waveform difference determining unit is used for solving the difference of the characteristic quantities of the adjacent periods, taking the absolute value of the difference, and taking the average value of the absolute values as the measurement of the difference of the characteristic quantity waveforms of the adjacent two periods;
the current waveform analyzing unit is used for analyzing the waveform difference characteristics of adjacent cycles of the current in the electric energy meter according to the following formula:
and setting the number of sampling points of current in each period as N, the number of periods N, and the sampling values of current in two adjacent periods as ik-1, ik:
Figure 424934DEST_PATH_IMAGE022
the average of the absolute values of the current differences Δ ik (j) over a period is:
Figure 510702DEST_PATH_IMAGE023
δ k reflects the magnitude of current difference between two adjacent periods, and is normalized to avoid the influence of current values on current deviation:
Figure 128765DEST_PATH_IMAGE024
the waveform similarity is:
Figure 71313DEST_PATH_IMAGE025
14. the system of claim 8, wherein the determine complex coupling relationship module comprises:
and the correlation coefficient calculation formula of the coupling relation of the complex coupling relation used for determining the complex electromagnetic field signal characteristics, the complex electric energy meter electricity utilization information characteristics and the complex electric energy meter working state characteristics is as follows:
Figure 587876DEST_PATH_IMAGE026
in the formula, n is sample energy, rho is a correlation coefficient, and x and y are corresponding elements in two variables.
15. A computer-readable storage medium, characterized in that the storage medium stores a computer program for performing the method of any of the preceding claims 1-7.
CN202210949666.6A 2022-08-09 2022-08-09 Comprehensive evaluation method and system for complex electromagnetic field environment Active CN115329565B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210949666.6A CN115329565B (en) 2022-08-09 2022-08-09 Comprehensive evaluation method and system for complex electromagnetic field environment

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210949666.6A CN115329565B (en) 2022-08-09 2022-08-09 Comprehensive evaluation method and system for complex electromagnetic field environment

Publications (2)

Publication Number Publication Date
CN115329565A true CN115329565A (en) 2022-11-11
CN115329565B CN115329565B (en) 2023-04-07

Family

ID=83922225

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210949666.6A Active CN115329565B (en) 2022-08-09 2022-08-09 Comprehensive evaluation method and system for complex electromagnetic field environment

Country Status (1)

Country Link
CN (1) CN115329565B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117435871A (en) * 2023-11-07 2024-01-23 深圳市赛盛技术有限公司 Method, system and medium for reducing electromagnetic interference of equipment based on magnetic field analysis

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5270883A (en) * 1991-08-29 1993-12-14 Mitsubishi Denki Kabushiki Kaisha Magnetic read/write circuit
CN101776740A (en) * 2010-01-19 2010-07-14 云南电网公司 Electric energy metering automation simulation laboratory
CN207780130U (en) * 2017-12-06 2018-08-28 国网重庆市电力公司电力科学研究院 A kind of research system of secondary device
CN110018376A (en) * 2019-05-08 2019-07-16 南华大学 A kind of carbon fibre composite thunder and lightning measurement of electromagnetic environment device and method
CN110097297A (en) * 2019-05-21 2019-08-06 国网湖南省电力有限公司 A kind of various dimensions stealing situation Intellisense method, system, equipment and medium
CN111460900A (en) * 2020-03-04 2020-07-28 中国工程物理研究院应用电子学研究所 Method for quantitative evaluation of equivalence of complex electromagnetic environment construction
CN112948744A (en) * 2019-11-26 2021-06-11 上海交通大学 RFID non-invasive sensing model construction method, device, equipment and medium
CN114814397A (en) * 2021-09-17 2022-07-29 中国电力科学研究院有限公司 Monitoring system and method for transient interference on site of converter station

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5270883A (en) * 1991-08-29 1993-12-14 Mitsubishi Denki Kabushiki Kaisha Magnetic read/write circuit
CN101776740A (en) * 2010-01-19 2010-07-14 云南电网公司 Electric energy metering automation simulation laboratory
CN207780130U (en) * 2017-12-06 2018-08-28 国网重庆市电力公司电力科学研究院 A kind of research system of secondary device
CN110018376A (en) * 2019-05-08 2019-07-16 南华大学 A kind of carbon fibre composite thunder and lightning measurement of electromagnetic environment device and method
CN110097297A (en) * 2019-05-21 2019-08-06 国网湖南省电力有限公司 A kind of various dimensions stealing situation Intellisense method, system, equipment and medium
CN112948744A (en) * 2019-11-26 2021-06-11 上海交通大学 RFID non-invasive sensing model construction method, device, equipment and medium
CN111460900A (en) * 2020-03-04 2020-07-28 中国工程物理研究院应用电子学研究所 Method for quantitative evaluation of equivalence of complex electromagnetic environment construction
CN114814397A (en) * 2021-09-17 2022-07-29 中国电力科学研究院有限公司 Monitoring system and method for transient interference on site of converter station

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
张卫东等: "高压直流换流站瞬态电磁骚扰特征的小样本统计方法" *
王沛等: "变电站内空间电磁场对无线通信单元耦合机理分析" *
范宇清等: "飞行器GNSS接收机静电放电干扰效应分析" *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117435871A (en) * 2023-11-07 2024-01-23 深圳市赛盛技术有限公司 Method, system and medium for reducing electromagnetic interference of equipment based on magnetic field analysis

Also Published As

Publication number Publication date
CN115329565B (en) 2023-04-07

Similar Documents

Publication Publication Date Title
Fang et al. Optimal periodic watermarking schedule for replay attack detection in cyber–physical systems
Baraldi et al. Robust signal reconstruction for condition monitoring of industrial components via a modified Auto Associative Kernel Regression method
Han et al. Wavelet packet based damage identification of beam structures
Biswas et al. Electricity theft pinpointing through correlation analysis of master and individual meter readings
Dervilis et al. On robust regression analysis as a means of exploring environmental and operational conditions for SHM data
Harrou et al. An improved multivariate chart using partial least squares with continuous ranked probability score
Yan et al. Structural damage diagnosis under varying environmental conditions—Part I: A linear analysis
Godoy et al. Relationships between PCA and PLS-regression
CN113518011B (en) Abnormality detection method and apparatus, electronic device, and computer-readable storage medium
Zhao et al. Online fault prognosis with relative deviation analysis and vector autoregressive modeling
JP6361175B2 (en) Abnormality diagnosis apparatus and program
Wang et al. Fault isolation based on residual evaluation and contribution analysis
Yin Data-driven design of fault diagnosis systems
Quiñones-Grueiro et al. An unsupervised approach to leak detection and location in water distribution networks
Rao et al. Sensor fault detection in structural health monitoring using null subspace–based approach
CN115329565B (en) Comprehensive evaluation method and system for complex electromagnetic field environment
JP2000259223A (en) Plant monitoring device
Mahan et al. White Noise Test: detecting autocorrelation and nonstationarities in long time series after ARIMA modeling.
Pan et al. Fault detection with improved principal component pursuit method
Chakour et al. Diagnosis of uncertain nonlinear systems using interval kernel principal components analysis: Application to a weather station
Han et al. Sensor drift detection based on discrete wavelet transform and grey models
Lan et al. Dynamic statistical process monitoring based on generalized canonical variate analysis
CN115018012A (en) Internet of things time sequence anomaly detection method and system under high-dimensional characteristic
US9535917B1 (en) Detection of anomalous utility usage
CN112380073B (en) Fault position detection method and device and readable storage medium

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