CN108593557B - Remote measurement error compensation method based on TE-ANN-AWF (transverse electric field analysis) -based mobile pollution source - Google Patents

Remote measurement error compensation method based on TE-ANN-AWF (transverse electric field analysis) -based mobile pollution source Download PDF

Info

Publication number
CN108593557B
CN108593557B CN201810205703.6A CN201810205703A CN108593557B CN 108593557 B CN108593557 B CN 108593557B CN 201810205703 A CN201810205703 A CN 201810205703A CN 108593557 B CN108593557 B CN 108593557B
Authority
CN
China
Prior art keywords
weight
awf
interference
ann
fusion
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201810205703.6A
Other languages
Chinese (zh)
Other versions
CN108593557A (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.)
Hangzhou Dianzi University
Original Assignee
Hangzhou Dianzi University
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 Hangzhou Dianzi University filed Critical Hangzhou Dianzi University
Priority to CN201810205703.6A priority Critical patent/CN108593557B/en
Publication of CN108593557A publication Critical patent/CN108593557A/en
Application granted granted Critical
Publication of CN108593557B publication Critical patent/CN108593557B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/18Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/08Learning methods
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N2021/1793Remote sensing
    • G01N2021/1795Atmospheric mapping of gases

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computational Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Software Systems (AREA)
  • Health & Medical Sciences (AREA)
  • General Engineering & Computer Science (AREA)
  • Pathology (AREA)
  • Evolutionary Biology (AREA)
  • Molecular Biology (AREA)
  • Computing Systems (AREA)
  • Computational Linguistics (AREA)
  • Biophysics (AREA)
  • Biomedical Technology (AREA)
  • Artificial Intelligence (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Immunology (AREA)
  • Evolutionary Computation (AREA)
  • Biochemistry (AREA)
  • Analytical Chemistry (AREA)
  • Operations Research (AREA)
  • Probability & Statistics with Applications (AREA)
  • Chemical & Material Sciences (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention discloses a remote measurement error compensation method based on a TE-ANN-AWF mobile pollution source. A virtual observation method is provided to realize the multi-element deconstruction of a unit observation sequence, the compensation of a single interference channel virtual observation sequence is realized through an ANN error prediction model, and then the compensated multi-element virtual observation sequence is subjected to fusion reconstruction by adopting a multi-element self-adaptive weighting fusion method. Aiming at the problem of weight convergence in the fusion algorithm, an exponential forgetting method is introduced into the model, and the advantages of good weight estimation capability of TE and self-adaptive adjustment of the weight of AWF are combined, so that the dynamic performance of the error compensation process is improved.

Description

Remote measurement error compensation method based on TE-ANN-AWF (transverse electric field analysis) -based mobile pollution source
Technical Field
The invention designs a remote sensing detection error compensation method for a mobile pollution source based on TE-ANN-AWF, which belongs to the technical field of error compensation of remote sensing detection instruments for the mobile pollution source.
Background
The mobile pollution source refers to a source that discharges air pollutants without a stationary device, such as a motor vehicle, a mobile construction machine, a ship, an airplane, and the like that discharge exhaust gas during movement. The pollution of mobile sources, particularly heavy diesel trucks, commercial gasoline vehicles, old automobiles and engineering machinery vehicles, has become a main source of air pollution and is an important cause of urban ultrafine particulate matters and photochemical smog pollution. The world health organization in 2012 has promoted diesel exhaust from suspected carcinogens to clear carcinogens. In order to effectively treat the mobile pollution source, especially to identify and control the high pollution discharge motor vehicle, the emission of the mobile pollution source needs to be effectively detected in real time.
At present, the detection technology of the mobile pollution source is various, such as reflecting the pollution degree of the atmosphere through plant growth or determining the concentration of the pollutant through chemical detection, but the methods are difficult to perform real-time online detection on the mobile pollution source. Remote sensing is a fast and efficient method of detecting mobile sources of pollution. The method is based on the precise analysis of the optical absorption spectrum of the gas, and the concentration of the gas is inverted according to the absorption properties of the gas components in the environment in ultraviolet, visible and infrared spectrum bands. In 1988, the first Remote Sensing monitoring system (RSD) using non-dispersive infrared technology (NDIR) was developed by the university of denver, usa, and was capable of simultaneously detecting emission concentrations of CO2, CO, and HC in the exhaust gas of a motor vehicle. In the 90 s of the 20 th century, the university of denver developed a remote sensing monitoring system for measuring the concentration of NO by using an ultraviolet absorption technology (Ultra Violet, UV), and the problem of water vapor absorption during NDIR measurement was overcome. Subsequently, MD-LaserTech developed NO and HC remote sensing systems based on the ultraviolet Differential technology (UV-DOAS). Due to the fact that the emitting wavelength band of a Tunable Diode Laser (TDL) is narrow, the wavelength band which is not interfered can be selected for measuring different measured gas components, and real-time monitoring requirements of high time resolution, high sensitivity, high selectivity and the like are met. Then, the massachusetts institute of technology developed a telemetry system based on the infrared laser differential absorption spectroscopy (tdlas) (tunable diode laser absorption spectrometer) for detecting the excessive emission of mobile pollution sources such as motor vehicles. Currently, remote sensing detection technology has become the mainstream technical means for real-time online detection of gas concentration.
The remote sensing measurement method has high automation degree, and a large number of vehicles passing through the road section can be measured as long as the equipment is erected beside the road. Meanwhile, the method has small influence on traffic. The university of purjohn morse, liriot, uk, TSI corporation, and others, have developed roadside monitoring systems that can detect pollutants including CO, NOx, and others, using technologies based on TDLAS, UV-DOAS, and others. However, the measuring method has limited use conditions and is greatly interfered by the environment, such as the ambient temperature, the humidity, the air pressure, the wind speed and the like, which have important influence on the detection data. And the urban geographic environment is mostly complex and changeable, such as: the urban canyon effect causes a sudden change in airflow. Due to the complexity of the external environment, the German Frounhf physical measurement technology is additionally provided with the measurement of meteorological parameters such as wind speed, wind direction, temperature and humidity, air pressure and the like on the basis of detecting the pollutant parameters, and the pollutant detection result is corrected.
Disclosure of Invention
Aiming at the problem that a remote sensing detection method of a mobile pollution source is easily interfered by an external environment, the invention provides a new error compensation model TE-ANN-AWF by combining three methods of TE transfer entropy, ANN artificial neural network and AWF self-adaptive weighting fusion. In the model, the causal correlation between interference and measurement is quantitatively analyzed by using TE transfer entropy, and a determination method of non-significant causal relationship is introduced. The idea of virtual observation is put forward in the model to realize the multi-element deconstruction of the unit observation sequence, and then the multi-element virtual observation sequence is reconstructed by a multi-element fusion method. The method for introducing the index forgetting into the model combines the advantages of the good weight estimation of the TE and the self-adaptive adjustment of the weight of the AWF, and improves the dynamic performance of the error compensation process.
The technical scheme of the invention is as follows:
the method comprises the following steps: measuring samples under different interference effects are obtained through an environment simulation experiment; based on the experimental measurement sample, performing interference correlation analysis through TE transfer entropy so as to determine a measurement error source and measure the unbalance degree between multiple interferences; the directivity of the TE transfer entropy is utilized to lead out the quantization standard and the judgment method of the non-obvious causal relationship;
step two: acquiring a training sample set of a single interference channel by adopting an environment simulation smoke box experiment platform, and establishing a measurement error prediction model of each interference by a neural network ANN method;
step three: realizing the multivariate deconstruction of the unit observation sequence by a virtual observation method, carrying out error compensation on the deconstructed multivariate virtual observation sequence by an ANN error prediction model under different interferences, and then carrying out fusion reconstruction on the compensated multivariate virtual observation sequence by adopting a multivariate self-adaptive weighting fusion method;
aiming at the problem of weight convergence in the fusion algorithm, an exponential forgetting method is introduced into the model, and the weight estimation capability of TE and the advantage of weight self-adaptive adjustment of the multivariate self-adaptive weighting fusion method are combined, so that the dynamic performance of the error compensation process is improved.
In the first step, TE transfer entropy is introduced to carry out correlation causal analysis on the interference of remote sensing measurement aiming at the characteristic that the interference of an external environment can be detected, and the directivity of the transfer entropy is utilized to introduce a quantitative standard and a judgment method of an insignificant causal relationship;
assuming that the interferences are mutually independent to a certain extent, a measurement sequence under the change of single interference is obtained through a simulation experiment platform, so that the transfer entropy TE from the temperature interference to the measured value is calculatedT->COThe transfer entropy TE of the humidity disturbance to the measured valueH->COTransfer entropy TE of barometric disturbance to measured valueP->COTransfer entropy TE of wind speed disturbance factor to measured valueW->COIn which the maximum backward transfer entropy TE is taken0As a balance of non-causal relationsA quantity standard;
TE0=max{TECO->T、TECO->H、TECO->P、TECO->W} (1)。
in the third step, the fusion reconstruction of the compensated multivariate virtual observation sequence by adopting the multivariate adaptive weighting fusion method specifically comprises the following steps: carrying out self-adaptive fusion on the multivariate observation values of all interference channels under the criterion of minimum mean square error; three methods of TE transfer entropy, ANN artificial neural network and AWF self-adaptive weighting fusion are closely combined with each other through a virtual observation method and an exponential forgetting mechanism.
In the third step, an exponential forgetting method is introduced into the model to combine the advantages of good weight estimation capability of TE and adaptive adjustment of AWF weight, specifically
In the TE-ANN model, although the transfer entropy has good weight estimation capability, the weight cannot be adjusted according to errors so as to gradually converge, so that the distribution of the errors also keeps fluctuating and cannot converge; the AWF adaptive weighting fusion algorithm has the advantage that the weight can be gradually converged under the criterion of minimum mean square error without any prior knowledge; in order to combine the good weight estimation advantage of TE with the self-adaptive weight adjustment advantage of AWF, a forgetting mechanism is introduced;
the obtained confidence weight K of the transmission entropy estimation and the optimal weighting factor W according to the minimum mean square error criterion*Selecting a weighting factor of βnGet K and W through fusion
Figure GDA0002532541900000041
As shown in the following formula;
Figure GDA0002532541900000042
wherein n represents the fusion of the nth observation,
Figure GDA0002532541900000043
a weighting factor, K, representing a confidence weight, KnTo representThe nth observation value is used for transmitting the confidence weight value of the entropy estimation,
Figure GDA0002532541900000044
a weighting factor, W, representing an optimal weighting factor Wn *And representing the weight of the AWF after the nth observation value is fused.
Considering the characteristics of stable estimation of the transmission entropy and the self-adaptive adjustment and convergence of the AWF weight, the estimated value of the initial dynamic process of the weight is mainly focused on KnWhile the steady state process should be focused on Wn *And gradually converge to Wn *In order to embody the above-mentioned characteristics, the weighting factor is βnThe following characteristics need to be satisfied:
Figure GDA0002532541900000045
i represents a natural number; (3)
in order to satisfy the above conditions, the following function is constructed;
dn=(1-b)/(1-a·bn),n=1,2,3… (4)
wherein b is a forgetting factor, a is an attenuation factor, and 0< b <1< a; the weighting coefficient is obtained from equation (3):
Figure GDA0002532541900000046
thus, the weight of each channel is obtained, as shown in the formula:
Figure GDA0002532541900000047
compared with the prior art, the invention has the advantages that:
(1) the invention provides a new error compensation model TE-ANN-AWF by combining three methods of TE transfer entropy, ANN artificial neural network and AWF self-adaptive weighting fusion. Compared with the traditional error compensation method, the method does not need to know a priori multiple interference measurement samples, can effectively compensate the measurement error caused by the external environment interference, and improves the applicability and the anti-interference capability of remote sensing measurement.
(2) According to the characteristic that external interference can be detected, the TE transfer entropy is introduced to carry out correlation causal analysis on the interference of remote sensing measurement, and the directivity of the transfer entropy is utilized to introduce the quantization standard and the judgment method of the non-obvious causal relationship.
(3) The invention introduces an adaptive weighted fusion algorithm AWF, and adaptively fuses the observed values of each interference channel under the criterion of minimum mean square error, thereby ensuring the long-term stability of error compensation.
(4) In order to further improve the dynamic performance of the error compensation process, the invention combines the advantages of the good weight estimation of TE and the self-adaptive adjustment of the weight of AWF by introducing a forgetting mechanism.
Drawings
FIG. 1 is a flow chart of the method of the present invention;
FIG. 2 is a diagram illustrating a comparison of transmission entropy for different interference factors according to the present invention;
FIG. 3 is a diagram of a prediction model of a telemetry error neural network under single interference according to the present invention;
FIG. 4 is a CO measurement error distribution diagram under temperature-concentration according to the present invention;
FIG. 5 is a diagram of the CO measurement error distribution under wind speed-concentration according to the present invention;
FIG. 6 is a diagram showing the error distribution of CO measurement under air pressure-concentration in accordance with the present invention;
FIG. 7 is a diagram illustrating reconstruction of a virtual observation sequence according to the present invention;
FIG. 8 is a diagram of a TE-ANN-AWF error compensation model of the present invention;
Detailed Description
In order to make the technical innovation point realized by the invention easy to understand, the following describes the realization mode of the invention in detail with reference to fig. 1, and the specific steps are as follows:
the method comprises the following steps: and obtaining measurement samples under different interference effects through an environment simulation experiment, and then preprocessing and normalizing the measurement samples. And performing interference correlation analysis through TE transfer entropy based on experimental measurement samples, thereby determining a measurement error source and measuring the degree of imbalance among multiple interferences. And the directivity of the TE transfer entropy is utilized to lead out the quantization standard and the judgment method of the non-significant causal relationship.
Let Xn and Yn be two sequences of environmental interference changes and telemetric observations with discrete states of Xn and Yn at time n, and Xn and Yn may be approximated as steady state markov processes of order k and order l, respectively, then the entropy of transfer from Yn and Xn is defined as follows:
Figure GDA0002532541900000051
wherein, TY→XDenotes the entropy of Transfer of Y to X (Transfer Entroy), un=(xn+1,xn,yn (l)),p(un) Represents a state xn+1And the sequence xn (k),yn (l)Probability of simultaneous occurrence; p (x)n+1|xn (k),yn (l)) Indicating that at time n, x is knownn (k),yn (l)On the premise of (1), xn+1The conditional probability of (a); p (x)n+1|xn (k)) Denotes xn (k)Given the known premise of xn+1When x is a conditional probability ofnWhen the state at a certain moment is completely determined by the historical state of the state, the transfer entropy is zero.
The remote sensing detection is based on the optical absorption principle to finish the detection of the mobile pollution source sewage gas, and the optical signal has the problems of absorption, scattering, beam deflection, beam diffusion and the like. Detection is mostly carried out in an outdoor environment, the measuring process is easily influenced by factors such as temperature, humidity, air pressure, wind speed, wind direction and dust, and the problem is a complex nonlinear dynamics problem under the combined action of a plurality of factors. The transfer entropy can quantitatively measure the linear and nonlinear relations among system variables, has good anti-noise capability, and is often applied to depicting dynamic nonlinear characteristics in a complex system. Therefore, the transfer entropy is selected in the compensation model to perform correlation analysis on the measurement data, which helps to trace back the source of error generation and compensate the error more reasonably and accurately.
In order to verify the effectiveness of the compensation model, CO is selected as a target detection object, the measurement environment of temperature, humidity, air pressure and wind speed changes is simulated through an environment simulation experiment platform, and the target object is measured by adopting remote measuring equipment. Transfer entropy is adopted in the compensation model to measure the imbalance degree between multiple interferences. Assuming that the interferences are mutually independent to a certain extent, a measurement sequence under the change of single interference is obtained through a simulation experiment platform, so that the transfer entropy TE from the temperature interference to the measured value is calculatedT->CO, transfer entropy TE of humidity disturbance to measured valueH->CO, transfer entropy TE of barometric disturbance to measured valueP->CO, transfer entropy TE of wind speed disturbance factors to measured valuesW->CO, with maximum reverse transfer entropy TE0 taken as a measure of non-causal relationship.
TE0=max{TECO->T、TECO->H、TECO->P、TECO->W}
As can be seen from FIG. 2, the transfer entropy TE of wind speed, temperature, barometric pressure to CO measurementW->CO、TET->CO、TEP->COAre clearly greater than the transfer entropy TE of humidity to CO measurementH->COAnd the transfer entropy TE of humidity to CO measurementH->COAnd reverse transfer entropy TE0The difference is not large. From the information theory perspective, the information contained in the measurement sequence can be obviously explained from the wind speed, temperature and air pressure interference sequence, and the explaining part can be obtained from the humidity sequence to be smaller. Therefore, the wind speed, the temperature and the air pressure have obvious causal association on the measurement result, namely, the three environmental factors have larger interference on the measurement, and the proportion of the three factors in the measurement error is larger. From this it can be seen the transfer entropy TE of the humidity to the measurementH->COAnd reverse transfer entropy TE0Very close, reflecting that there is no obvious causal relationship of humidity to CO measurements. The effect of humidity may thus not be taken into account in the compensation of the CO measurement.
Step two: acquiring a training sample set of a single interference channel by adopting an environment simulation smoke box experiment platform, and establishing a measurement error prediction model of each interference by a neural network ANN method;
according to the graph shown in FIG. 3, the prediction model of the neural network for the CO gas telemetry error under temperature interference is shown. The telemetering error is not only influenced by interference factors, but also the absolute value of the error is influenced by the real concentration of the gas to be detected. Therefore, the input data of the neural network is the standard gas concentration n and the temperature t, and the output data is the measurement error e. And (4) selecting 7 hidden layer numbers according to a hidden layer point selection method of the forward network and the effect of actually processing data. Through smog chamber environmental simulation platform, fixed other factors regulate and control the incasement temperature to this training sample under obtaining single temperature interference. Similarly, a neural network prediction model under wind speed interference and air pressure interference can be established. Fig. 4-6 are error distributions for a test sample under single temperature, single air pressure, and single wind speed disturbances, respectively.
After three prediction models of temperature, wind speed and air pressure are established, error prediction can be respectively carried out on the measured values under certain temperature, wind speed and air pressure, and then the three prediction errors are summed to obtain a final compensation value.
Step three: according to the TE-ANN-AWF model structure diagram shown in FIG. 8, the transfer entropy TE, the neural network ANN and the adaptive weighting fusion AWF are combined for final estimation of the disturbance measurement value through virtual observation deconstruction and an exponential forgetting mechanism.
1) Firstly, realizing the multivariate deconstruction of a unit observation sequence by a virtual observation method, carrying out corresponding error compensation on the deconstructed multivariate virtual observation sequence by an ANN error prediction model under different interferences, and then carrying out fusion reconstruction on the compensated multivariate virtual observation sequence by adopting a multivariate self-adaptive weighting fusion method AWF.
Aiming at an ANN error prediction model, the real concentration of the gas to be detected directly influences the error distribution. In fact, the actual concentration of the gas to be measured cannot be known, and only approximate estimation can be performed through a historical sequence, which also causes inaccurate predicted measurement errors. Therefore, the error compensation values for the respective channels are directly added, and the effect is not obvious. According to fig. 7, in order to apply the method of data fusion to the error compensation model, the observation sequence needs to be deconstructed. For this purpose, the concept of virtual observation is first introduced.
For example, the remote sensing detection of CO is interfered by temperature, air pressure and wind speed in the environment, and if the mutual coupling of the three interference factors on the measurement result is small, the error component can be regarded as additive noise. Thus, the measured value can be regarded as represented by the following formula, where Y is the measured value, yr is the true value, Wt, Ww, Wp are the environmental interference noise, and θ is the measurement random noise.
Y=yr+Wt+Ww+Wp
Wherein, Wt, Ww, Wp are the measurement error of single interference, which can be predicted by the single interference ANN model.
The measured value at a single temperature disturbance is Yt, and it can be considered as a value excluding both of Ww and Wp, as shown in the following equation.
Yt=yr+Wt+θ=Y-Ww-Wp
It is clear that the measured value Yt at temperature single interference does not exist in practice. In fact, it is a reconstruction of the original observations and is therefore referred to as a "virtual observation". Similarly, the measured value Yp under the single disturbance of the air pressure and the measured value Yw under the single disturbance of the wind speed can be obtained as shown in the following formula.
Yw=yr+Ww+θ=Y-Wt-Wp
Yp=yr+Wp+θ=Y-Wt-Ww
After the observation is deconstructed by the concept of virtual observation, reconstruction is also needed. Wind speed has a greater influence on the measurement results than temperature, barometric pressure, especially when wind speed is greater than 5 m/s. Obviously, the confidence of the measured value at this time should be minimal. To express the magnitude of this confidence, a transitive entropy is introduced to express the confidence of the three virtual observations. It can be known from the transmission entropy analysis that the larger the transmission entropy of the interference is, the stronger the causal relationship of the interference measurement is, and the stronger the interference is on the measurement result. The confidence of the virtual observations is thus inversely proportional to the transfer entropy. Meanwhile, in order to keep the total weight of 1, the weights of the three virtual observations whose confidence levels are different need to be calculated, and the confidence levels of the three virtual observations can be obtained as follows.
Figure GDA0002532541900000081
Figure GDA0002532541900000082
Figure GDA0002532541900000083
Wherein TET、TEP、TEWCan be estimated according to the principle of transfer entropy, and KT、KP、KWSatisfying the following equation.
Kt+Kp+Kw=1
And reconstructing by using the estimated confidence degrees of different virtual observations as the weight values of the virtual observations. The measurement results after reconstruction are shown in the following formula.
Figure GDA0002532541900000084
2) Secondly, aiming at the problem of weight convergence in the fusion algorithm, an exponential forgetting method is introduced to combine the advantages of good weight estimation capability of TE and self-adaptive adjustment of the weight of AWF, so that the dynamic performance of the error compensation process is improved.
In the TE-ANN model, although the transfer entropy has a good capability of weight estimation, the weight cannot be adjusted according to the error so as to gradually converge, so that the distribution of the error also keeps fluctuating and cannot converge. The AWF adaptive weighting fusion algorithm has the advantage that the weight can be gradually converged under the criterion of minimum mean square error without any prior knowledge. In order to combine the advantages of good weight estimation of TE and the advantages of adaptive adjustment of AWF weight, a forgetting mechanism is introduced.
Pre-entropy of the obtained transferEstimated confidence weight K and optimal weighting factor W according to minimum mean square error criterion*Selecting a weighting coefficient { β n }, and adding K to W*The fusion was carried out as shown in the following formula.
Figure GDA0002532541900000091
Where n represents the fusion of the nth observation.
Considering the characteristics of the smooth estimation of the transmission entropy and the adaptive adjustment and convergence of the weight of the AWF, the estimated value of the initial dynamic process of the weight is mainly referred to Kn, and the steady-state process is mainly referred to Wn and gradually converges to Wn. In order to embody the above-described characteristics, the weighting coefficients { β n } need to satisfy the following characteristics:
Figure GDA0002532541900000092
in order to satisfy the above conditions, the following function is constructed.
dn=(1-b)/(1-a·bn),n=1,2,3…
Wherein b is a forgetting factor, a is an attenuation factor, and 0< b <1< a. By
Figure GDA0002532541900000093
The available weighting coefficients:
βn K=dn,βn W=1-dn
therefore, the weight of each channel can be obtained, as shown in the formula:
Figure GDA0002532541900000094
wherein Wn ^ satisfies the following formula, and the obtained new weight Wn ^ is used as a new weight to replace the original Wn ^.
Figure GDA0002532541900000095
The above examples are provided only for the purpose of describing the present invention, and are not intended to limit the scope of the present invention. The scope of the invention is defined by the appended claims. Various equivalent substitutions and modifications can be made without departing from the spirit and principles of the invention, and are intended to be within the scope of the invention.

Claims (4)

1. The method for compensating the telemetering error of the mobile pollution source based on the TE-ANN-AWF specifically comprises the following steps:
the method comprises the following steps: measuring samples under different interference effects are obtained through an environment simulation experiment; based on the experimental measurement sample, performing interference correlation analysis through TE transfer entropy so as to determine a measurement error source and measure the unbalance degree between multiple interferences; the directivity of the TE transfer entropy is utilized to lead out the quantization standard and the judgment method of the non-obvious causal relationship; ignoring interference factors which have non-significant causal relationship with the measurement result in the subsequent step;
step two: acquiring a training sample set of a single interference channel by adopting an environment simulation smoke box experiment platform, and establishing a measurement error prediction model of each interference by a neural network ANN method;
step three: realizing the multivariate deconstruction of the unit observation sequence by a virtual observation method, carrying out error compensation on the deconstructed multivariate virtual observation sequence by an ANN error prediction model under different interferences, and then carrying out fusion reconstruction on the compensated multivariate virtual observation sequence by adopting a multivariate self-adaptive weighting fusion method;
aiming at weight convergence in the fusion algorithm, an exponential forgetting method is introduced into the model to combine the weight pre-estimation capability of TE and the weight self-adaptive adjustment of the multivariate self-adaptive weighting fusion method, so that the dynamic performance of the error compensation process is improved.
2. The TE-ANN-AWF-based mobile pollution source telemetry error compensation method of claim 1, wherein: in the first step, TE transfer entropy is introduced to carry out correlation causal analysis on the interference of remote sensing measurement aiming at the characteristic that the interference of an external environment can be detected, and the directivity of the transfer entropy is utilized to introduce a quantitative standard and a judgment method of an insignificant causal relationship;
assuming that the interferences are mutually independent to a certain extent, a measurement sequence under the change of single interference is obtained through a simulation experiment platform, so that the transfer entropy TE from the temperature interference to the measured value is calculatedT->COThe transfer entropy TE of the humidity disturbance to the measured valueH->COTransfer entropy TE of barometric disturbance to measured valueP->COTransfer entropy TE of wind speed disturbance factor to measured valueW->COIn which the maximum backward transfer entropy TE is taken0As a measure of non-causal relationships;
TE0=max{TECO->T、TECO->H、TECO->P、TECO->W} (1)。
3. the TE-ANN-AWF-based mobile pollution source telemetry error compensation method of claim 1, wherein: in the third step, the fusion reconstruction of the compensated multivariate virtual observation sequence by adopting the multivariate adaptive weighting fusion method specifically comprises the following steps: carrying out self-adaptive fusion on the multivariate observation values of all interference channels under the criterion of minimum mean square error; three methods of TE transfer entropy, ANN artificial neural network and AWF self-adaptive weighting fusion are closely combined with each other through a virtual observation method and an exponential forgetting mechanism.
4. The TE-ANN-AWF-based mobile pollution source telemetry error compensation method of claim 1, wherein: in the third step, an exponential forgetting method is introduced into the model to combine the TE weight estimation capability and the AWF weight adaptive adjustment, and the method specifically comprises the following steps:
the obtained confidence weight K of the transmission entropy estimation and the optimal weighting factor W according to the minimum mean square error criterion*Selecting a weighting factor of βnH, mixing K and W*Carrying out fusion to obtain
Figure FDA0002532541890000021
As shown in the following formula;
Figure FDA0002532541890000022
wherein n represents the fusion of the nth observation,
Figure FDA0002532541890000023
a weighting factor, K, representing a confidence weight, KnRepresenting the confidence weight of the transmission entropy estimation of the nth observed value,
Figure FDA0002532541890000024
represents the optimal weighting factor W*Weighting coefficient of Wn *Representing the weight of the AWF after the nth observation value is fused;
considering the characteristics of stable estimation of the transmission entropy and the self-adaptive adjustment and convergence of the AWF weight, the estimated value of the initial dynamic process of the weight is mainly focused on KnWhile the steady state process should be focused on Wn *And gradually converge to Wn *In order to embody the above-mentioned characteristics, the weighting factor is βnThe following characteristics need to be satisfied:
βi=βi-1b;0<b<1;
Figure FDA0002532541890000025
i represents a natural number; (3)
in order to satisfy the above conditions, the following function is constructed;
dn=(1-b)/(1-a·bn),n=1,2,3…(4)
wherein b is a forgetting factor, a is an attenuation factor, and 0< b <1< a; the weighting coefficient is obtained from equation (3):
Figure FDA0002532541890000031
thus, the weight of each channel is obtained, as shown in the formula:
Figure FDA0002532541890000032
CN201810205703.6A 2018-03-13 2018-03-13 Remote measurement error compensation method based on TE-ANN-AWF (transverse electric field analysis) -based mobile pollution source Active CN108593557B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810205703.6A CN108593557B (en) 2018-03-13 2018-03-13 Remote measurement error compensation method based on TE-ANN-AWF (transverse electric field analysis) -based mobile pollution source

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810205703.6A CN108593557B (en) 2018-03-13 2018-03-13 Remote measurement error compensation method based on TE-ANN-AWF (transverse electric field analysis) -based mobile pollution source

Publications (2)

Publication Number Publication Date
CN108593557A CN108593557A (en) 2018-09-28
CN108593557B true CN108593557B (en) 2020-08-11

Family

ID=63626189

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810205703.6A Active CN108593557B (en) 2018-03-13 2018-03-13 Remote measurement error compensation method based on TE-ANN-AWF (transverse electric field analysis) -based mobile pollution source

Country Status (1)

Country Link
CN (1) CN108593557B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109374532B (en) * 2018-12-17 2021-03-30 杭州电子科技大学 Pollution emission telemetering error compensation method based on transfer entropy and self-adaptive fusion
CN109492830B (en) * 2018-12-17 2021-08-31 杭州电子科技大学 Mobile pollution source emission concentration prediction method based on time-space deep learning
CN112825157B (en) * 2019-11-20 2022-10-04 天津大学 Gasification gas production prediction method, device, equipment and storage medium

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5137616A (en) * 1991-04-04 1992-08-11 Surface Combustion, Inc. Gas analysis system for furnaces and the like
CN105445344A (en) * 2015-12-30 2016-03-30 桂林电子科技大学 Temperature compensation method of system for detecting heavy metals in water environment
CN106878375A (en) * 2016-12-22 2017-06-20 中国民航大学 A kind of cockpit pollutant monitoring method based on distribution combination sensor network
CN107170219A (en) * 2017-04-24 2017-09-15 杭州电子科技大学 A kind of mobile pollution source on-line monitoring system and method
CN107300550A (en) * 2017-06-21 2017-10-27 南京大学 A kind of method based on BP neural network model prediction atmosphere heavy metal concentration
CN107423467A (en) * 2017-04-18 2017-12-01 武汉科技大学 A kind of three-dimensional pollution sources localization method close to lake bank

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20110010318A1 (en) * 2007-08-17 2011-01-13 Institutt For Energiteknikk System and method for empirical ensemble- based virtual sensing
FR2994272B1 (en) * 2012-08-06 2014-09-05 Commissariat Energie Atomique METHOD OF LOCATING ORIGIN OF GAS FLOWS IN A GEOGRAPHICAL SPACE, INVOLVING A SELECTION OF MEASUREMENTS
US10112618B2 (en) * 2015-10-11 2018-10-30 Rimalu Technologies, Inc. Traffic pollution indicator

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5137616A (en) * 1991-04-04 1992-08-11 Surface Combustion, Inc. Gas analysis system for furnaces and the like
CN105445344A (en) * 2015-12-30 2016-03-30 桂林电子科技大学 Temperature compensation method of system for detecting heavy metals in water environment
CN106878375A (en) * 2016-12-22 2017-06-20 中国民航大学 A kind of cockpit pollutant monitoring method based on distribution combination sensor network
CN107423467A (en) * 2017-04-18 2017-12-01 武汉科技大学 A kind of three-dimensional pollution sources localization method close to lake bank
CN107170219A (en) * 2017-04-24 2017-09-15 杭州电子科技大学 A kind of mobile pollution source on-line monitoring system and method
CN107300550A (en) * 2017-06-21 2017-10-27 南京大学 A kind of method based on BP neural network model prediction atmosphere heavy metal concentration

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
Finding the Direction of Disturbance Propagation in a Chemical Process Using Transfer Entropy;Margret Bauer et al.;《TRANSACTIONS ON CONTROL SYSTEMS TECHNOLOGY》;20070131;第12-21页 *
On-road Heavy-duty Vehicle Emissions Monitoring System;Gary A. Bishop et al.;《Environ. Sci. Technol》;20150121;第1639-1645页 *
Recurrent Air Quality Predictor Based on Meteorology- and Pollution-Related Factors;Ke Gu et al.;《TRANSACTIONS ON INDUSTRIAL INFORMATICS》;20180115;第3946-3955页 *
Roadside automobile emission monitoring with Peltier-cooled diode laser spectrometer;J. F. Kastner et al.;《Proc. SPIE》;19970522;第103-109页 *
应用TDLAS探测气体的误差分析与补偿研究;矫晓敏, 尚丽平;《传感器与微系统》;20081231;第70-72,76页 *

Also Published As

Publication number Publication date
CN108593557A (en) 2018-09-28

Similar Documents

Publication Publication Date Title
US20180321208A1 (en) Determining the net emissions of air pollutants
CN108593557B (en) Remote measurement error compensation method based on TE-ANN-AWF (transverse electric field analysis) -based mobile pollution source
Arhami et al. Predicting hourly air pollutant levels using artificial neural networks coupled with uncertainty analysis by Monte Carlo simulations
Irie et al. Eight-component retrievals from ground-based MAX-DOAS observations
CN106920007B (en) PM based on second-order self-organizing fuzzy neural network2.5Intelligent prediction method
Timmreck et al. The interactive stratospheric aerosol model intercomparison project (ISA-MIP): motivation and experimental design
Shi et al. An inversion method for estimating strong point carbon dioxide emissions using a differential absorption Lidar
CN109374532B (en) Pollution emission telemetering error compensation method based on transfer entropy and self-adaptive fusion
CN111579504A (en) Atmospheric pollution component vertical distribution inversion method based on optical remote sensing
JP7365415B2 (en) Methods for calibrating gas sensors
Liguori et al. Estimation of the minimum measurement time interval in acoustic noise
CN113758890A (en) Gas concentration calculation method, device, equipment and storage medium
Watne et al. Tackling data quality when using low-cost air quality sensors in citizen science projects
De Vito et al. Dynamic multivariate regression for on-field calibration of high speed air quality chemical multi-sensor systems
Šavli et al. The prospects for increasing the horizontal resolution of the Aeolus horizontal line‐of‐sight wind profiles
Reshma Analysis and prediction of air quality
Ma et al. On-line wavenumber optimization for a ground-based CH4-DIAL
Shichkin et al. Training algorithms for artificial neural networks for time series forecasting of greenhouse gas concentrations
CN110322015B (en) Vehicle inspection data generation method
CN114814092A (en) IP index measuring method based on BP neural network
Li et al. Back-propagation neural network-based modelling for soil heavy metal
CN114485800A (en) Remote quality control method suitable for gas multi-parameter mobile monitor
CN114565136A (en) Air quality prediction optimization method based on generation countermeasure network
Barbes et al. The use of artificial neural network (ANN) for prediction of some airborne pollutants concentration in urban areas
Xi et al. An error compensation method for remote sensing measurement of mobile source emissions

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