CN116663433B - Differential optimization method based on differential absorption spectrometer - Google Patents
Differential optimization method based on differential absorption spectrometer Download PDFInfo
- Publication number
- CN116663433B CN116663433B CN202310942708.8A CN202310942708A CN116663433B CN 116663433 B CN116663433 B CN 116663433B CN 202310942708 A CN202310942708 A CN 202310942708A CN 116663433 B CN116663433 B CN 116663433B
- Authority
- CN
- China
- Prior art keywords
- absorption
- differential
- gas
- spectrum
- section
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000010521 absorption reaction Methods 0.000 title claims abstract description 92
- 238000000034 method Methods 0.000 title claims abstract description 78
- 238000005457 optimization Methods 0.000 title claims abstract description 22
- 239000007789 gas Substances 0.000 claims abstract description 81
- 238000001228 spectrum Methods 0.000 claims abstract description 31
- 238000001914 filtration Methods 0.000 claims abstract description 22
- 238000000862 absorption spectrum Methods 0.000 claims abstract description 19
- 230000003287 optical effect Effects 0.000 claims abstract description 17
- 230000008033 biological extinction Effects 0.000 claims abstract description 13
- 230000000694 effects Effects 0.000 claims abstract description 7
- 230000002068 genetic effect Effects 0.000 claims description 13
- 238000005070 sampling Methods 0.000 claims description 9
- 230000003595 spectral effect Effects 0.000 claims description 7
- 238000009499 grossing Methods 0.000 claims description 6
- 238000002835 absorbance Methods 0.000 claims description 3
- 230000006978 adaptation Effects 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 3
- 108090000623 proteins and genes Proteins 0.000 claims description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Chemical compound O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 3
- 238000005259 measurement Methods 0.000 abstract description 2
- 238000005516 engineering process Methods 0.000 description 5
- 230000005540 biological transmission Effects 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 235000013405 beer Nutrition 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 239000012535 impurity Substances 0.000 description 1
- 230000035772 mutation Effects 0.000 description 1
- 230000005070 ripening Effects 0.000 description 1
- 238000010183 spectrum analysis Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/27—Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/12—Computing arrangements based on biological models using genetic models
- G06N3/126—Evolutionary algorithms, e.g. genetic algorithms or genetic programming
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Abstract
The invention discloses a differential optimization method based on a differential absorption spectrometer, which belongs to the technical field of electric data processing and is used for carrying out differential optimization of optical data, and comprises the steps of obtaining wavelength and light intensity information of an absorption spectrum through the differential absorption spectrometer, calculating a concentration value of gas according to the absorption light intensity in the spectrum, filtering out a slow change part with low dependence on the wavelength in the spectrum according to the gas extinction effect by using a Bill law, and removing the slow change part and Rayleigh scattering and Mie scattering extinction coefficients to obtain differential optical density only representing gas molecular absorption; and broadening the absorption spectrum lines of the characteristic wavelength bands corresponding to various gases by using the Fogart function, performing least square fitting on the differential optical density and the reference spectrum, and inverting to obtain the final gas concentration value in the current measurement environment. The method improves and optimizes the traditional differential absorption algorithm, improves the precision for solving various gas concentration values, and can adapt to various complex environments.
Description
Technical Field
The invention discloses a differential optimization method based on a differential absorption spectrometer, and belongs to the technical field of electric data processing.
Background
The general instrument can only detect single gas concentration, can not detect multiple gas concentrations simultaneously, and how to detect various polluted gases in real time and alarm according to local air quality standards is the research focus in the current environment. The differential absorption spectrum technology makes up the defects of the traditional instrument and can detect the concentration of various gases simultaneously; the differential absorption spectrum technology takes solar scattered light or artificially introduced detection light as a light source, obtains the wavelength and the light intensity of the scattered light, and inverts the concentration information of various gases through a differential absorption algorithm after the light source information is obtained, so that the real-time concentration values of the various gases are obtained. The differential absorption spectrometer has the advantages of easy construction of a platform, multiple inversion gas concentration types, multi-platform observation and the like. For the conventional differential absorption technology, the standard absorption section is generally obtained by an ultraviolet-visible band high-resolution standard absorption section in a spectrum database, and the absorption section data is obtained by a gas absorption section of a high-resolution spectrometer at normal temperature, in practical situations, a plurality of different gas concentrations are required to be obtained, and the absorption sections of the gases at different temperatures are different, which leads to an increase of data errors. In addition, in a low-concentration measurement environment, due to the influence of a plurality of environmental conditions such as optical path, concentration, impurities and the like, the information quantity of the concentration of the gas to be measured in the absorption spectrum is small, and a larger error can be generated in the gas concentration value solved by the differential absorption algorithm, and even the error can be more than 10%. Aiming at the defects of the differential absorption algorithm, the traditional algorithm is optimized and improved, so that the method can meet the requirements of multiple conditions and high precision, and the concentration of various gases to be measured can be more accurately obtained under more complex conditions.
Disclosure of Invention
The invention aims to provide a differential optimization method based on a differential absorption spectrometer, which aims to solve the problem of large data error of the differential absorption technology in the prior art.
A differential optimization method based on a differential absorption spectrometer, comprising:
s1, acquiring wavelength and light intensity information of an absorption spectrum by a differential absorption spectrometer, calculating a concentration value of gas from the absorption light intensity in the spectrum, and filtering out a slow change part in the spectrum according to the gas extinction effect by using a Bill law;
s2, separating a fast change part and a slow change part in an absorption spectrum through digital filtering;
s3, removing the slow change part and the Rayleigh scattering and Mie scattering extinction coefficient to obtain the differential optical density only representing the absorption of gas molecules;
S4, broadening absorption spectral lines of the corresponding characteristic wavelength bands of various gases by using the Focus function;
s5, after an effective absorption section is obtained through data processing, the effective absorption section is used as a reference spectrum, and the differential optical density and the reference spectrum are subjected to least square fitting;
s6, inverting the gas concentration, judging, and adding a catastrophe period self-adaptive crossover operation based on a genetic algorithm to optimize;
s7, inverting to obtain a final gas concentration value in the current measuring environment.
S1 comprises the following steps:
the relation between the fast and slow changing parts of the spectrum and the transmitted light intensity is as follows:
;
in the method, in the process of the invention,is the transmitted intensity of ultraviolet light at a certain wavelength, +.>For the incident light intensity of ultraviolet light at a certain wavelength, n is the number of gas species, +.>Representing the broadband structure spectrum, i.e. the slowly varying part of the spectrum,/->For differential absorption cross section, i.e. fast changing part of the spectrum,/->Represents the concentration of the ith species of substance at the transmission path s,/->And->Rayleigh scattering and Mie scattering extinction coefficients, respectively.
S2 comprises the following steps:
the high-pass filtering adopts six-order polynomial fitting to fit the absorption spectrum in the wave band:
;
in the method, in the process of the invention,for the intensity of the light after fitting, +.>Is a coefficient of->Is the wavelength.
S2 comprises the following steps:
the low-pass filtering adopts Savez-Golay smoothing denoising:
;
in the method, in the process of the invention,for smoothed data +.>For smooth coefficient +.>Is the pre-smoothing data.
S3 comprises the following steps:
。
s4 comprises the following steps:
fugget functionThe method comprises the following steps:
;
wherein the method comprises the steps ofRepresenting Doppler spread half-width, < >>Indicating that the gas is inBeam->The absorption line at the point is strong, and the temperature dependence of F is expressed as follows, according to the influence of temperature on the absorption cross section:
;
indicating the reference temperature +.>Strong absorption line and->For the actual temperature +.>For reference temperature->Is the energy of the low energy state of the molecule, +.>Is Planck constant, < >>For the speed of light->Is a boltzmann constant, t represents a variable constant;
x and y are intermediate coefficients:
,/>;
a central beam representing each absorption line, +.>Representing the half width of the pressure broadening.
S4 comprises the following steps:
half-width of Doppler broadeningAnd half width of pressure broadening->The temperature and pressure dependence of (2) are expressed as:
,/>;
in the method, in the process of the invention,taking 1 for symmetrical molecules and 1.5 for asymmetrical molecules in a dimensionless manner; />Indicating the actual atmospheric pressure, +.>Representing a reference spread value, +.>Represents a reference atmospheric pressure, M represents a molecular mass of a gas;
according to the temperature and the atmospheric pressure in the actual situation, comparing the reference temperature and the atmospheric pressure, calculating to obtain the current high-resolution water vapor absorption section, and convoluting the obtained reference section with a spectrometer instrument function curve to obtain the current effective absorption section of each gas.
S5 comprises the following steps:
is provided withIndicate->Seed gas inversion concentration according to ∈>The formula (1) is:
,/>;
will beIs rewritten as a matrix form:
;
in the method, in the process of the invention,represents +.>Component (F)>Line i corresponds to->Component, thereby obtaining->Differential optical Density of seed gas->And effective absorption cross section->Further reversing the concentration of the generated gas>。
S6 comprises the following steps:
s6.1. coding with decimal system, in which each individual comprisesIndividual gene coding value->,/>For the coded gas concentration, each population consisted of 10 random initialized individuals and satisfied:
;
s6.2, selecting fitness functionThe method comprises the following steps:
;
in the method, in the process of the invention,is->Sampling points->For the number of spectral sampling point wavelengths, +.>For the encoded gas concentration, +.>For differential absorption cross section->For the +.>Absorbance values for each sampling point.
S6 comprises the following steps:
s6.3, selecting operation, namely calculating an average value of fitness, and selecting individuals smaller than the average value of the fitness to enter the next generation;
s6.4, adding disaster period self-adaptive crossover operation to define crossover probabilityThe method comprises the following steps:
;
in the method, in the process of the invention,represents the crossover probability->Indicates the number of catastrophic events that have occurred at the present time, < + >>Represents a constant value of less than 1, +.>Represents an integer for adjusting the range of cross probability variation, +.>Indicating that the population is proceeding to->Performing genetic operation;
s6.5, terminating operation, and outputting an individual with the minimum adaptation value as a true value by the optimization algorithm when the solution concentration value of the last individual reaches a required threshold value, wherein the error range of the concentration value finally solved by the optimization algorithm is within the threshold value range.
Compared with the prior art, the invention has the following beneficial effects: the accuracy of solving various gas concentration values is improved by improving and optimizing a traditional differential absorption algorithm, and the method can be suitable for various complex environments;
in order to adapt to more different environments, the standard absorption section can be effectively and accurately obtained at different temperatures, the absorption spectrum lines of the corresponding characteristic wavelength bands of various gases are widened through a Voigt function, and then the obtained reference section is convolved with a spectrometer instrument function curve, so that the current effective absorption section of each gas is finally obtained, and the error is reduced;
in order to adapt to the environment with small concentration and short optical path, a genetic algorithm is integrated into a final data processing process, and the algorithm is optimized by coding, defining a fitness function, selecting, inheriting and terminating, so that the accuracy of data is improved;
in order to reduce the 'early ripening' phenomenon of the genetic algorithm, in the optimization of the genetic algorithm, the catastrophe period self-adaptive crossover operation is added, the global searching capability of population variation is improved by carrying out self-adaptive change on crossover probability, and the adaptability and the accuracy of the algorithm are enhanced.
Drawings
FIG. 1 is a technical flow chart of the present invention;
FIG. 2 is a diagram of a digital filtering process;
FIG. 3 is a schematic diagram of the cross section of the Fogart function for effective absorption.
Detailed Description
For the purpose of making the objects, technical solutions and advantages of the present invention more apparent, the technical solutions in the present invention will be clearly and completely described below, and it is apparent that the described embodiments are some embodiments of the present invention, but not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
A differential optimization method based on a differential absorption spectrometer, comprising:
s1, acquiring wavelength and light intensity information of an absorption spectrum by a differential absorption spectrometer, calculating a concentration value of gas from the absorption light intensity in the spectrum, and filtering a slow change part with low dependence on the wavelength in the spectrum according to the gas extinction effect by using a Bill law;
s2, separating a fast change part and a slow change part in an absorption spectrum through digital filtering;
s3, removing the slow change part and the Rayleigh scattering and Mie scattering extinction coefficient to obtain the differential optical density only representing the absorption of gas molecules;
S4, broadening absorption spectral lines of the corresponding characteristic wavelength bands of various gases by using the Focus function;
s5, after an effective absorption section is obtained through data processing, the effective absorption section is used as a reference spectrum, and the differential optical density and the reference spectrum are subjected to least square fitting;
s6, inverting the gas concentration, judging, and if the gas concentration is measured in an environment with lower gas concentration, measuring the gas concentration with larger error, and optimizing by adding the catastrophe period self-adaptive cross operation based on a genetic algorithm;
s7, inverting to obtain a final gas concentration value in the current measuring environment.
S1 comprises the following steps:
the relation between the fast and slow changing parts of the spectrum and the transmitted light intensity is as follows:
;
in the method, in the process of the invention,is the transmitted intensity of ultraviolet light at a certain wavelength, +.>For the incident light intensity of ultraviolet light at a certain wavelength, n is the number of gas species, +.>Representing the widthWith structural spectrum, i.e. slow varying part of the spectrum,/->For differential absorption cross section, i.e. fast changing part of the spectrum,/->Represents the concentration of the ith species of substance at the transmission path s,/->And->Rayleigh scattering and Mie scattering extinction coefficients, respectively.
S2 comprises the following steps:
the high-pass filtering adopts six-order polynomial fitting to fit the absorption spectrum in the wave band:
;
in the method, in the process of the invention,for the intensity of the light after fitting, +.>Is a coefficient of->Is the wavelength.
S2 comprises the following steps:
the low-pass filtering adopts Savez-Golay smoothing denoising:
;
in the method, in the process of the invention,for smoothed data +.>For smooth coefficient +.>Is the pre-smoothing data.
S3 comprises the following steps:
。
s4 comprises the following steps:
fugget functionThe method comprises the following steps:
;
wherein the method comprises the steps ofRepresenting Doppler spread half-width, < >>Indicating that the gas is in beam->The absorption line at the point is strong, and the temperature dependence of F is expressed as follows, according to the influence of temperature on the absorption cross section:
;
indicating the reference temperature +.>Strong absorption line and->For the actual temperature +.>For referenceTemperature (F)>Is the energy of the low energy state of the molecule, +.>Is Planck constant, < >>For the speed of light->Is a boltzmann constant, t represents a variable constant;
x and y are intermediate coefficients:
,/>;
a central beam representing each absorption line, +.>Representing the half width of the pressure broadening.
S4 comprises the following steps:
half-width of Doppler broadeningAnd half width of pressure broadening->The temperature and pressure dependence of (2) are expressed as:
,/>;
in the method, in the process of the invention,taking 1 for symmetrical molecules and 1.5 for asymmetrical molecules in a dimensionless manner; />Indicating the actual atmospheric pressure, +.>Representing a reference spread value, +.>Represents a reference atmospheric pressure, M represents a molecular mass of a gas;
according to the temperature and the atmospheric pressure in the actual situation, comparing the reference temperature and the atmospheric pressure, calculating to obtain the current high-resolution water vapor absorption section, and convoluting the obtained reference section with a spectrometer instrument function curve to obtain the current effective absorption section of each gas.
S5 comprises the following steps:
is provided withIndicate->Seed gas inversion concentration according to ∈>The formula (1) is:
,/>;
will beIs rewritten as a matrix form:
;
in the method, in the process of the invention,represents +.>Component (F)>Line i corresponds to->Component, thereby obtaining->Differential optical Density of seed gas->And effective absorption cross section->Further reversing the concentration of the generated gas>。
S6 comprises the following steps:
s6.1. coding with decimal system, in which each individual comprisesIndividual gene coding value->,/>For the coded gas concentration, each population consisted of 10 random initialized individuals and satisfied:
;
s6.2, selecting fitness functionThe method comprises the following steps:
;
in the method, in the process of the invention,is->Sampling points->For the number of spectral sampling point wavelengths, +.>For the encoded gas concentration, +.>For differential absorption cross section->For the +.>Absorbance values for each sampling point.
S6 comprises the following steps:
s6.3, selecting an individual with a smaller fitness value to enter the next generation according to the fitness function;
s6.4, adding disaster period self-adaptive crossover operation to define crossover probabilityThe method comprises the following steps:
;
in the method, in the process of the invention,represents the crossover probability->Indicates the number of catastrophic events that have occurred at the present time, < + >>Represents a constant value of less than 1, +.>Represents an integer for adjusting the range of cross probability variation, +.>Indicating that the population is proceeding to->Performing genetic operation;
s6.5, terminating operation, and outputting an individual with the minimum adaptation value as a true value by the optimization algorithm when the solution concentration value of the last individual reaches a required threshold value, wherein the error range of the concentration value finally solved by the optimization algorithm is within the threshold value range.
The technical flow of the invention is shown in figure 1, and comprises the steps of obtaining spectral information, high-pass filtering and low-pass filtering, obtaining differential optical density, obtaining an effective absorption section of gas, fitting by a least square method, then directly outputting a concentration value in a normal concentration environment and a low concentration environment according to different conditions, carrying out a catastrophe genetic algorithm, and outputting an optimal individual meeting a threshold value.
In general, atmospheric scattering comprises Rayleigh scattering and Mie scattering, the Rayleigh scattering and Mie scattering are used as absorption processes when calculated by a differential absorption spectrum technology, and the relationship between the emitted light intensity and the received light intensity is obtained by combining Rayleigh scattering, mie scattering extinction effect and a lambert beer law. The differential absorption algorithm is divided into two parts of rapid change with high wavelength dependence and slow change with low wavelength dependence according to the dependence of the gas extinction effect, the slow change part with low wavelength dependence is needed to be filtered out in the spectrum analysis, the differential optical density only representing the absorption of gas molecules is obtained after the slow change part is filtered out, and the differential optical density and a gas molecule reference absorption section are subjected to least square fitting to obtain the concentration of various gases. Separating the "fast-changing" and "slow-changing" portions of the absorption spectrum typically employs digital filtering, removing the "slow-changing" portions of the spectrum using high-pass digital filtering, such as polynomial regression filtering, and reducing the effects of high-frequency noise using low-pass filtering, such as triangular filtering. The reference absorption section is usually selected by directly using a standard section for processing, the corresponding wavelength is recorded according to the characteristic absorption spectrum peak of different gases in the standard absorption section at a certain wavelength, the data corresponding to the wavelength in the differential absorption section and the data of the standard absorption section are subjected to least square fitting, and the digital filtering process is shown in figure 2.
Based on the gas absorption parameters provided by the spectrum database, under different temperature conditions, characteristic absorption spectrum lines of characteristic wavelength bands of various gases are stretched by utilizing Fogart function fitting, and in the fitting, it is found that various interference gases exist in the characteristic wavelength bands corresponding to certain gases, so that in order to improve the accuracy, firstly, the gas characteristic wavelength Duan Zhonggan is selected, the interference gases exist less, then, all the interference gases existing in the wavelength band are stretched, and then, the absorption cross sections obtained after the stretching participate in the spectrum fitting. Under complex conditions such as low concentration, the genetic algorithm and the traditional differential absorption algorithm are combined and improved, the fitness function is defined according to the gas absorption characteristics, concentration values and the like, an initial data set is generated through operations such as selection, crossover, mutation and the like of the genetic algorithm, iterative optimization is carried out on the data, each generation of the data set is evaluated according to the fitness function, individuals with higher fitness are selected as a population of the next generation, genetic operation is carried out, new individuals are generated, iterative operation of multiple generations is repeated until the individuals with the highest fitness are found, data errors are reduced, and the Fogart function obtains an effective absorption section as shown in a figure 3.
The above embodiments are only for illustrating the technical aspects of the present invention, not for limiting the same, and although the present invention has been described in detail with reference to the foregoing embodiments, it will be understood by those skilled in the art that: the technical solutions described in the foregoing embodiments may be modified or some or all of the technical features may be replaced with other technical solutions, which do not depart from the scope of the technical solutions of the embodiments of the present invention.
Claims (6)
1. The differential optimization method based on the differential absorption spectrometer is characterized by comprising the following steps of:
s1, acquiring wavelength and light intensity information of an absorption spectrum by a differential absorption spectrometer, calculating a concentration value of gas from the absorption light intensity in the spectrum, and filtering out a slow change part in the spectrum according to the gas extinction effect by using a Bill law;
s2, separating a fast change part and a slow change part in an absorption spectrum through digital filtering;
s3, removing the slow change part and the Rayleigh scattering and Mie scattering extinction coefficient to obtain the differential optical density only representing the absorption of gas molecules;
;
In the method, in the process of the invention,is the transmitted intensity of ultraviolet light at a certain wavelength, +.>For the incident light intensity of ultraviolet light at a certain wavelength,nfor the number of gas species>For differential absorption cross section, i.e. fast changing part of the spectrum,/->Represent the firstiSpecies concentration at the delivery path s;
s4, broadening absorption spectral lines of the corresponding characteristic wavelength bands of various gases by using the Focus function;
fugget functionThe method comprises the following steps:
;
wherein the method comprises the steps ofRepresenting Doppler spread half-width, < >>Indicating that the gas is in beam->The absorption line at the point is strong, and the temperature dependence of F is expressed as follows, according to the influence of temperature on the absorption cross section:
;
indicating the reference temperature +.>Strong absorption line and->For the actual temperature +.>For reference temperature->Is the energy of the low energy state of the molecule, +.>Is Planck constant, < >>For the speed of light->Is a boltzmann constant, t represents a variable constant;
xandyis the intermediate coefficient:
,/>;
a central beam representing each absorption line, +.>Representing the half width of the pressure broadening;
half-width of Doppler broadeningAnd half width of pressure broadening->The temperature and pressure dependence of (2) are expressed as:
,/>;
in the method, in the process of the invention,taking 1 for symmetrical molecules and 1.5 for asymmetrical molecules in a dimensionless manner; />Indicating the actual atmospheric pressure, +.>Representing a reference spread value, +.>Indicating that the reference atmospheric pressure is to be used,Mrepresenting the molecular mass of the gas;
comparing the reference temperature with the atmospheric pressure according to the temperature and the atmospheric pressure in the actual situation, calculating to obtain a current high-resolution water vapor absorption section, and convoluting the obtained reference section with a spectrometer instrument function curve to obtain a current effective absorption section of each gas;
s5, after an effective absorption section is obtained through data processing, the effective absorption section is used as a reference spectrum, and the differential optical density and the reference spectrum are subjected to least square fitting;
s6, inverting the gas concentration, and adding catastrophe period self-adaptive crossover operation based on a genetic algorithm to optimize;
s6.1. coding with decimal system, in which each individual comprisesIndividual gene coding value->,/>For the coded gas concentration, each population consisted of 10 random initialized individuals and satisfied:
;
s6.2, selecting fitness functionThe method comprises the following steps:
;
in the method, in the process of the invention,is->Sampling points->For the number of spectral sampling point wavelengths, +.>For the encoded gas concentration, +.>For differential absorption cross section->For the +.>Absorbance values for the individual sampling points;
s7, inverting to obtain a final gas concentration value in the current measuring environment.
2. The differential optimization method based on a differential absorption spectrometer according to claim 1, wherein S1 comprises:
the relation between the fast and slow changing parts of the spectrum and the transmitted light intensity is as follows:
;
in the method, in the process of the invention,representing the broadband structure spectrum, i.e. the slowly varying part of the spectrum,/->And->Rayleigh scattering and Mie scattering extinction coefficients, respectively.
3. The differential optimization method based on a differential absorption spectrometer according to claim 1, wherein S2 comprises:
the high-pass filtering adopts six-order polynomial fitting to fit the absorption spectrum in the wave band:
;
in the method, in the process of the invention,for the intensity of the light after fitting, +.>Is a coefficient of->Is the wavelength.
4. The differential optimization method based on a differential absorption spectrometer according to claim 1, wherein S2 comprises:
the low-pass filtering adopts Savez-Golay smoothing denoising:
;
in the method, in the process of the invention,for smoothed data +.>For smooth coefficient +.>Is the pre-smoothing data.
5. The differential optimization method based on a differential absorption spectrometer according to claim 1, wherein S5 comprises:
is provided withIndicate->Seed gas inversion concentration according to ∈>The formula (1) is:
,/>;
will beIs rewritten as a matrix form:
;
in the method, in the process of the invention,represent the firstiLine 1mEffective absorption cross-section component corresponding to column, +.>First, theiRow corresponding +.>Component, thereby obtaining->Differential optical Density of seed gas->And an effective absorption cross section, thereby reversing the concentration of the process gas +.>。
6. The differential optimization method based on a differential absorption spectrometer according to claim 1, wherein S6 comprises:
s6.3, selecting operation, namely calculating an average value of fitness, and selecting individuals smaller than the average value of the fitness to enter the next generation;
s6.4, adding disaster period self-adaptive crossover operation to define crossover probabilityThe method comprises the following steps:
;
in the method, in the process of the invention,represents the crossover probability->Indicates the number of catastrophic events that have occurred at the present time, < + >>Represents a constant value of less than 1, +.>Represents an integer for adjusting the range of cross probability variation, +.>Indicating that the population is proceeding to->Performing genetic operation;
s6.5, terminating operation, and outputting an individual with a small adaptation value as a true value by the optimization algorithm when the solution concentration value of the last individual reaches a required threshold value, wherein the error range of the concentration value finally solved by the optimization algorithm is within the threshold value range.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310942708.8A CN116663433B (en) | 2023-07-31 | 2023-07-31 | Differential optimization method based on differential absorption spectrometer |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310942708.8A CN116663433B (en) | 2023-07-31 | 2023-07-31 | Differential optimization method based on differential absorption spectrometer |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116663433A CN116663433A (en) | 2023-08-29 |
CN116663433B true CN116663433B (en) | 2023-10-24 |
Family
ID=87722775
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310942708.8A Active CN116663433B (en) | 2023-07-31 | 2023-07-31 | Differential optimization method based on differential absorption spectrometer |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116663433B (en) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103226095A (en) * | 2013-04-08 | 2013-07-31 | 南京国电环保科技有限公司 | Fast calibration method of spectrometer wavelength |
CN106526614A (en) * | 2016-10-25 | 2017-03-22 | 中国科学院合肥物质科学研究院 | Method for optimizing laser radar detection atmospheric composition spectral line analysis |
CN111572729A (en) * | 2020-04-07 | 2020-08-25 | 哈尔滨工程大学 | Thrust distribution method of ship dynamic positioning system based on improved genetic algorithm |
CN115523958A (en) * | 2022-10-15 | 2022-12-27 | 浙江大学 | Gas temperature and concentration synchronous measurement method based on spectrum fast-slow separation principle |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9305257B2 (en) * | 2013-05-20 | 2016-04-05 | International Business Machines Corporation | Adaptive cataclysms in genetic algorithms |
-
2023
- 2023-07-31 CN CN202310942708.8A patent/CN116663433B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103226095A (en) * | 2013-04-08 | 2013-07-31 | 南京国电环保科技有限公司 | Fast calibration method of spectrometer wavelength |
CN106526614A (en) * | 2016-10-25 | 2017-03-22 | 中国科学院合肥物质科学研究院 | Method for optimizing laser radar detection atmospheric composition spectral line analysis |
CN111572729A (en) * | 2020-04-07 | 2020-08-25 | 哈尔滨工程大学 | Thrust distribution method of ship dynamic positioning system based on improved genetic algorithm |
CN115523958A (en) * | 2022-10-15 | 2022-12-27 | 浙江大学 | Gas temperature and concentration synchronous measurement method based on spectrum fast-slow separation principle |
Non-Patent Citations (5)
Title |
---|
LED光源DOAS方法的夜间NO_2气体浓度遥测研究;杨雷 等;光谱学与光谱分析(第05期);全文 * |
Maintaining Healthy Population Diversity Using Adaptive Crossover, Mutation, and Selection;Brian McGinley .etal;IEEE;第15卷(第5期);全文 * |
差分吸收光谱技术在工业污染源烟气排放监测中的应用;孙友文 等;物理学报;第62卷(第01期);第2-4页第2-5节 * |
差分吸收光谱检测H2S气体的反演算法;陈书旺 等;太赫兹科学与电子信息学报;第18卷(第4期);第735页第1节;第736页第2.2节,图3 * |
改进灾变遗传算法在无功优化规划中的应用;林广明 等;电网技术;第34卷(第4期);第129页第2.1节;第130页第2.2节 * |
Also Published As
Publication number | Publication date |
---|---|
CN116663433A (en) | 2023-08-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101105446A (en) | Differential optical absorption spectroscopy air quality detection system | |
CN108918446B (en) | Ultra-low concentration sulfur dioxide ultraviolet difference feature extraction algorithm | |
CN112881321B (en) | Correction method for measuring aerosol extinction coefficient by black carbon instrument | |
CN113049507A (en) | Multi-model fused spectral wavelength selection method | |
CN117349683A (en) | Auto-parts application colour difference anomaly detection system based on spectral data | |
CN116663433B (en) | Differential optimization method based on differential absorption spectrometer | |
CN114384548A (en) | Biological aerosol number concentration profile Raman fluorescence laser radar system and prediction method | |
CN112666104A (en) | DOAS-based gas concentration inversion method | |
CN112964662B (en) | Method for measuring concentration and temperature of high-temperature fuel gas of aircraft engine | |
CN107941758B (en) | Coal quality heat value measuring method based on spectrum pretreatment | |
CN114062306B (en) | Near infrared spectrum data segmentation preprocessing method | |
CN111794744A (en) | Method for monitoring pollution degree of formation water in real time underground | |
CN109521002B (en) | Fuel characteristic measuring method for solid fuel particle flow | |
CN106323888B (en) | Ultra-low emission flue gas measuring method | |
CN112881333B (en) | Near infrared spectrum wavelength screening method based on improved immune genetic algorithm | |
He et al. | Accurate inversion of tropospheric bottom temperature using pure rotational Raman lidar in polluted air condition | |
CN110672552B (en) | Confidence coefficient estimation method for vehicle fuel oil near infrared spectrum detection result | |
CN114720433A (en) | Automatic removing method for useless scattering of fluorescence spectrum | |
CN113702305B (en) | Gas concentration linear measurement method based on self-adaptive differential absorption spectrum technology | |
CN114324206B (en) | Sky condition classification method based on passive multi-axis differential absorption spectrum technology | |
CN110248456B (en) | Real-time automatic analysis method for diagnosing spectrum of low-temperature plasma laser Thomson scattering | |
CN117648537B (en) | Atmospheric pollution real-time monitoring method and system based on hyperspectral technology | |
CN107727597B (en) | Ultraviolet differential concentration inversion and interference suppression method for mixed gas | |
CN111487211B (en) | Incoherent broadband cavity enhanced absorption spectrum fitting waveband selection method | |
CN116008208A (en) | Seawater nitrate concentration characteristic spectrum band selection method |
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 |