Cite

INTRODUCTION

With the rapid development of economic and people living standard unceasing improvement the popularity of the development of heavy industry and transportation, with a large number of harmful substances into the environment atmosphere, changed the composition of normal air, worsen air quality[1]. Health problems, which caused by the atmospheric pollution has caused governments and people of great importance to it. Therefore, monitoring SO2 and NOx become one of the important subject of environmental protection[2]. For the accurate monitoring of the real-time quality, ecological environment and pollution to the environment, to provide accurate basis for the environmental protection departments at various levels supervision, management and environment decision, urgent need a large number of modern environmental monitoring instruments.

Through to the present situation of flue gas analyzer in the domestic and foreign research and analysis, portable flue gas on the domestic market at present the main monitoring method is electrochemical analysis instrument, infrared gas analyzer and differential absorption spectrum. Electrochemical analysis instrument each sensor can only be measured a chemical composition in flue gas, if need to measure six ingredients you need six sensors. The chemical battery principle led to zero drift in the development of the system and cross interference, sensor lifetime is also difficult to solve the problem[3]. There is no absorption peak to NO2 in the infrared spectrum analyzer, cannot monitor NO2, but can only detect the NO content, through the assumptions to NO2 concentrations. At the same time as infrared spectrum analyzer adopts wheel filter method, each time only one gas measurement, can not measure a variety of gases at the same time[4]. Differential absorption measured gas to uv and visible light waves, on the basis of the differential absorption spectrum by the strength of the differential absorption spectrum inversion gas concentrations, in the traditional algorithm of differential absorption, the application of least square method of gas concentration for calculating optimal estimation equation is complex, and easy in pathological conditions[5]. Accordingly based on Lambert Beer law and absorbance of dual stack, a mutually recursive iteration is presented more gas concentration decoding algorithm, the process of the algorithm is presented, and verify the effectiveness of the algorithm. This system using the Ocean Optic micro spectrometer combining with embedded technology company has developed a set of portable flue gas analyzer[6]. This instrument is to overcome the disadvantages of the electrochemical sensor, had been achieved and sensor life problems, and at the same time for accurate measurement of gases and no cross interference.

WORKING PRINCIPLE OF THE SYSTEM

Portable flue gas analyzer system structure consists of three parts, part data collection, data processing part and the human-computer interaction part. Data acquisition part by sampling probe, flue gas sampling pump, an ultraviolet light and miniature spectrum analyzer, collecting gas mainly accomplished by ultraviolet light and miniature spectrum analyzer spectrum before and after the data is converted into digital signal transmitted to the data processing part. Data processing part of ARM microcontroller CPU core module and peripheral circuit, main is to electrical signals were collected by using corresponding algorithm to calculate the gas concentration of all kinds of gas in the flue gas composition. The human-computer interaction part consists of LCD display and keyboard, mainly to complete the real-time display of gas concentration, and through the keyboard to control the storage of data and related parameters Settings.

THE BASIC PRINCIPLE OF LOOP ITERATION METHOD
Lambert-Beerlaw

The mathematical model of Lambert-Beer’s law [7-8] saidA=λγ(1/T)=Kβχ

In the formula: A is absorbance; T is the transmittance, project the light intensity is than the incident light intensity; C is light-absorbing substance concentration; B is the absorption layer thickness. Its physical meaning is: when a beam of parallel light through a vertical uniform a material suction light scattering, the absorbance and the concentration and absorption material suction light is directly proportional to the thickness. Absorbance [16] A binary additive think: if the binary and multicomponent mixture components have absorbed A wave number, the total absorbance at the wave number is equal to the arithmetic of absorbance and at all levels.

The Mathematical model

The logarithmic of Lambert-Beer’s law said:A=lgRλDλSλDλ=Kc

In the formula: A is absorbance; Rλ is photon number (zero point); Sλ is the number of through the photon; λ is selected for a wavelength; K is constant (associated with the length of wavelength and absorption tube); C is a single gas concentration; Dλ is dark spectrum photon number (related to the integration time can take a fixed constant).

Two kinds of mixed gas Lambert-Beer absorbance expression can be represented as:A1=lgRλ1DλSλ1DλA2=lgRλ2DλSλ2Dλ

According to the law of superposition of binary:A=A1+A2lgRλDλSλDλ=lgRλ1DλSλ1Dλ+lgRλ2DλSλ2Dλ

i.e.RλDλSλDλ=Rλ1DλSλ1Dλ×Rλ2DλSλ2Dλ

By formula (7) can be know in a particular incident light wavelength lambda, the mixed gas of photon number and subtract its dark photons through the photon number, equal to the selected wavelength under different gas on the number of incident photons and through the photon number minus the dark photon number of the product.

In the actual gas absorption model of gas absorption and scattering all play a role in the formation of spectrum, so the experiment of gas absorption spectra contain two parts. Part of gas molecules on the absorption of photons, ultraviolet band is mainly caused by the electron transition of the gas molecules absorbed. Another part of the scattering is gas or smoke without selective absorption caused by light attenuation. Because this system adopts the removable detection method, gas extraction, you first after a pretreatment device for flue gas filtration, drying, cooling and other processing. Filtration process can be 0. Lum or above level particle filter out completely, and flue gas molecules of the diameter of the order of magnitude are below 1 nm. Instrument of the selected work between the wavelength of 190 nm - 290nm, greater than the molecules or residual particle diameter, so that the gases to be detected is mainly gas molecules Rayleigh scattering.

The actual operation process we with nitrogen gas as zero, the spectral curve is called the zero gas, nitrogen gas molecules, although about 190 nm - 290 nm uv absorption, but also produces Rayleigh scattering. That is zero gas line is light after zero gas scattering spectral lines. By spectrogram, found that in the absence of absorption bands, such as nitric oxide line outside the three absorption peak is coincidence with zero gas lines, namely in the absence of gas absorption bands of two gases scattering light is the same. So the Rayleigh scattering of gas absorption interference can be ruled out by the method of using the zero gas line as a reference.

According to the Lambert Beer-law, exclude mie scattering and Rayleigh scattering in gas detection, then the device applicable gas absorption model becomes:ln(I(λ)I0(λ))=L(σi(λ)Ci)

Which represents the light source through I0(λ) zero gas scattering by spectrometer in wavelength λ after the received light intensity, σi (λ) represents the ith kind of gases at wavelength λ absorption cross section, L is the length of the pool by absorption, and Ci is the ith the concentration of the gas. Which gas is equal to the total absorbance and of each gas absorption degree.

CIRCULAR ITERATION STEPS

Circular iteration method using various gases in 190 ~ 290 nm band characteristic absorption peak, combined with the absorbance of dual stack, a feature in some gas absorption point assuming that other gas no absorption, launched the initial concentration of the gas, then switch to another characteristic absorption point, the gas absorbed photon number from the measured total absorbed photon number subtracting, get another gas initial concentration, and so on to get the initial concentration of each gas. Then return to the first, the characteristics of the gas absorption point from the total number of photons absorbed minus every other gas absorption of the photon number, again to get the first gas concentration of an iteration, and so on to get other gas concentrations of an iteration, the repeated iteration until the concentration difference between two times smaller than a certain value, the concentration of each gas to the gas concentration accurately. Algorithm steps are as follows.

First step: to solve the first gas in the gas mixture of initial concentration C10, and selects a characteristic wavelength λ1, under the specific wavelength of a gas in the gas mixture has obvious characteristic absorption peak, and other gas absorption peak is small at this point, is to read the value of the absorbed photon number Sλ1 and solve RλDλSλDλ is obtained by absorbance look-up table and the concentration of the gas as the initial concentration of mixed gas in the gas.

Step 2: solving the second gas in the gas mixture of initial concentration: select the second characteristic wavelength in the wavelength of λ2 second kind of gas has obvious absorption peak, and other gas absorption peak is weak, then read the absorbed photon number Sλ2 considered under the band is only two kinds of mixed gas, according to formula (7) to calculate the second gas absorbance, look-up table against solving calculation, the second gas concentrations, as the second kind of initial concentration C20 of gases.

Step 3: calculating of other gases in the mixed gases, the initial concentration: the mth characteristic wavelengths selected λm, the wavelength of the first m gas has obvious absorption peak, read the wavelength of absorption of the photon number Sλm by formula (7) and look-up table type inversion calculation of this kind of gas concentration as the initial concentration Cm0 of the gas.

Step 4: iterative inversion to calculate the first gas concentration of first order recursive: will the desires of all kinds of gas concentration in formula (7), again read λ1 of Sλ1 reverse the first gas concentration C11 of first order recursive.

Step 5: repeat the second and third step to calculate gas m first order recursive concentration C1m.

Step 6: calculation of gas concentration of the adjacent two iterative error; Calculating first-order differential iteration of each gas is shown in the following formula,Δ1=|ComCm1|

Select Δ1 maximum of first-order differential iteration as the iteration error, i.eΔ1=max{Δm1,m=1,2,,M}

Step 7: repeat the fourth, fifth and sixth step, until the iteration error is less than the given value namely.Δ1=max{Δm1,m=1,2,,M}

Step 8: termination of the algorithm, it will be the last time the concentration as a final concentration of gases.

THE EXPERIMENTAL RESULTS AND ANALYSIS

The obtained by experiments, So2No2, the NO and NH3 four gas absorption lines as shown in figure 1, in the range of 190 nm to 190 nm, four types of gas are interfered with each other. At wavelength of 1 (273.33nnl) So2 and No2 are absorbed. Also at wavelength of 2 (231.33nm) So2 and No2 are absorbed. Only at wavelength of 3 (225.88nm) NO, So2 and No2 are absorbed. At wavelength of 4 (208.23nm) NH3, No2 and So2 are absorbed. The So2 and No2 have absorbed at each wavelength point, For NH3 and NO, and, as long as the concentration of So2 and No2 first came out, and then the corresponding concentration on wavelength 3 or 4 of the absorbance minus the can get their absorbance. So should first of all, the concentration of So2 and No2, then the concentration of NH3 and NO.

Figure 1.

Mixed gas absorption spectrum

Calculate the concentration of So2 and No2

For the concentration of NH3 and NO are the premise of know the concentration of So2 and No2. However, So2 and No2 throughout the working wave band are interfered with each other, need an algorithm to eliminate their interference. Here, on the basis of the look-up table using loop iteration method and calculation method to determine the concentration of So2 and No2.

As shown in figure 2, At the wavelength of 1 and 2 at the So2 and No2 were absorbed, they differ greatly in the absorption intensity at two wavelengths, the absorbance of the absorbance of So2 and No2 at 273.33nm is far less than the same concentration of No2 in 231.33nm, and the So2 absorbance at 231.33nm is far less than the same concentration of So2 in the absorbance at 273.33nm. Then we can be absorbed in 231.33nm No2 called the main absorption of No2, No2 in the absorption of 273.33nm is called So2 interference, So2 in the absorption at 273.33nm is the main absorption, in 231.33 absorption interference on No2 main absorption. Iterative method is through circulation calculation to gradually eliminate the interference of No2 on So2, the true concentration approaching So2 and No2.

Figure 2.

The absorption spectrum of 2

Cyclic iteration method to overcome the conflict there understand equations caused no solution of the problem, and this method can be easily implemented by software programming. Statistical chart concentrations of So2 and No2 calculated with iteration number as shown in the figure 3-4. The abscissa is the number of iterations, the longitudinal concentration coordinates So2 and No2 (ppm).

Figure 3.

The results of calculation curves of iteration number and So2 concentration

Figure 4.

The results of calculation curves of iteration number and No2 concentration

According to the running results and statistics can be seen in figure, the process of iteration is a successive approximation process, the calculation results of So2 with the increase in the number of iterations decreases gradually, and tends to be stable, but No2 is exactly the opposite, results with increasing iteration times increase. Cyclic iterative method in two iterations later effect is not obvious, which is to say as long as two iterations. This can be calculated through No2 and So2 concentrations.

Calculate the concentration of NO and NH3

Because the NO and NH3 absorption in the selected wavelengths of 3 and a wavelength of 4 does not interfere with each other, so it can be separated by the concentration of total absorbance and the So2 and NO : calculation. Hypothesis has been obtained in the mixed gas of No2 and So1 for C1 concentration and C2 concentration, the NO and NH3 for C3 and C4, then at the wavelength of 3 according to the superposition of absorbance of the available,A3=AA1A2=lg(I0I)A1A2

A1 is the C1 concentration of No2 absorbance at 225.8nm, A2 is the C2 concentration of So2 absorbance at 225.8nm, A is a mixture of gases in the 225.88nm total absorbance at 225.88nm, A3 is NO of total absorbance. I0 is the spectral intensity of 225.88nm through zero gas, I is the intensity of 225.88nm through transmission into the mixed gas after, can be obtained directly by the spectrometer, A2 and A3 concentrations of A and No2 obtained by. Absorbance of A3 so you can get the NO in 225.88nm, then according to the concentration and absorbance at 225.88nm correspondence between the NO table to calculate the concentration of NO. The absorbance of A4NH3 can be obtained with the same method in 208.23nm, then according to the corresponding relationship between concentration and absorbance at 208.23nm NH3 for NH3 concentration.

CONCLUSION

Main harmful components for atmospheric environmental pollution monitoring requirements, using uv wavelength grating type continuous frequency measuring method, precision is put forward to solve the various harmful ingredient concentration of recursive iteration fast inversion algorithm, and validates the effectiveness of the algorithm. Portable flue gas analyzer based on this algorithm is based on embedded technology, sensor based on micro spectrometer data collection, using uv light source through the spectrum analysis method analysis of flue gas concentration, achieved through the use of a miniature spectrum analyzer to a variety of gas composition at the same time for the purpose of accurate measurement. The product has compact structure, high measuring accuracy, strong anti-interference, high sensitivity etc, has a broad application prospect and popularization value.

eISSN:
2470-8038
Language:
English
Publication timeframe:
4 times per year
Journal Subjects:
Computer Sciences, other