Method for estimating concentration of nitrogen trioxide in atmosphere

Abstract

The application provides a method for estimating the concentration of nitrogen trioxide in the atmosphere, which comprises the following steps: S10, measuring irradiance D1, D2, D3,..., Dn; S20, measuring nitrogen trioxide photolysis rate values K1, K2, K3,..., Kn; S30, establishing a first relationship Y1 between the D1-Dx and the corresponding nitrogen trioxide photolysis rate K1-Kx; establishing a second relationship Y2 between the Dx-Dn and the corresponding nitrogen trioxide photolysis rate Kx-Kn; S40, establishing a third relationship Y3 and a fourth relationship Y4; S50, final relationship K=(a1*D+b1*D2)+(a2*D+b2*D2); S60, obtaining the irradiance value Dm at a specific time, and bringing the Dm into the K to calculate the nitrogen trioxide concentration at the specific time. The method for estimating the concentration of nitrogen trioxide in the atmosphere provided by the application is accurate in measurement.

Description

TECHNICAL FIELD

[0001] The present application belongs to the field of environmental protection, and particularly relates to a method for estimating the concentration of nitrogen trioxide in the atmosphere. BACKGROUND

[0002] Photochemical reactions play an important role in atmospheric circulation. Free radicals generated by the photolysis of chemical substances in the atmosphere can induce or participate in the reactions of other chemical substances, promoting the conversion of primary pollutants into secondary pollutants.

[0003] Nitrogen trioxide (NO3) is also an important intermediate in the photochemical reaction cycle in the atmosphere. Nitrogen dioxide (NO2) can react with ozone (O3) to generate NO3 and oxygen (O2) and further generate nitrogen monoxide (NO) and O2. This reaction is an important reaction in the troposphere, and when the concentrations of NO2 and O3 are high, it is the main source of NO3 in the atmosphere. This reaction does not require light and can occur at night. When the concentrations of OH and NO are not high at night, NO2 and NO3 can further react to generate N2O5. NO3 radicals and N2O5 are the main species in the nighttime atmospheric environmental chemical cycle, controlling the removal and conversion of nitrogen oxides and volatile organic compounds (VOCs), and making important contributions to the generation of nitrate and secondary organic aerosols at night. It is one of the important contents of current atmospheric complex pollution research in China. Due to the large concentration level span of NO3, the short atmospheric lifetime, the rapid change with time, and the interference of other nitrogen oxides during measurement, the measurement of NO3 is difficult. Therefore, there are still great challenges in accurate measurement. The existing online monitoring technologies for nitrogen trioxide in the environment mainly include wet chemical method and spectroscopic method. The wet chemical method consumes chemical reagents and requires frequent reagent preparation. The generated waste liquid is acidic, polluting the environment, and the instrument maintenance is not convenient. The spectroscopic method does not require reagents, but the optical principle instrument is easily affected by particulate matter, humidity, and interfering substances. SUMMARY

[0004] The present application provides a method for estimating the concentration of nitrogen trioxide in the atmosphere, which aims to estimate the concentration of nitrogen trioxide in the atmosphere, improve the accuracy of estimation, and reduce environmental pollution.

[0005] To achieve the above purpose, the technical solution adopted by the present application is to provide a method for estimating the concentration of nitrogen trioxide in the atmosphere, comprising:

[0006] S10: measuring the irradiance D1, D2, D3,..., Dn corresponding to the wavelength bands λ1, λ2,..., λn respectively by using a radiance monitoring device;

[0007] S20: Measure the nitric oxide photolysis rate values K1, K2, K3, …, Kn corresponding to the irradiance D1, D2, D3, …, Dn respectively by using a chemical luminometer;

[0008] S30: Select the irradiance D1~Dx corresponding to the wave band λ1~λx, and establish a first relationship Y1 between the D1~Dx and the corresponding nitric oxide photolysis rate K1~Kx;

[0009] Select the irradiance Dx~Dn corresponding to the wave band λx~λn, and establish a second relationship Y2 between the Dx~Dn and the corresponding nitric oxide photolysis rate Kx~Kn, wherein 1

[0010] S40: Measure the nitric oxide concentration C1, C2, C3, …, Cn corresponding to the nitric oxide photolysis rate K1, K2, K3, …, Kn respectively by using a nitric oxide analyzer, establish a third relationship Y3 between the nitric oxide concentration C1~Cx and the corresponding nitric oxide photolysis rate value K1~Kx, and establish a fourth relationship Y4 between the nitric oxide concentration Cx~Cn and the corresponding nitric oxide photolysis rate value Kx~Kn;

[0011] S50: Convert the first relationship Y1 and the third relationship Y3 to each other to obtain the relationship K=a1*D+b1*D corresponding to the wave band λ1~λx 2 , and convert the second relationship Y2 and the fourth relationship Y4 to each other to obtain the relationship K=a2*D+b2*D corresponding to the wave band λx~λn 2 , wherein a1, a2, b1, b2 are constants, and the final relationship K=(a1*D+b1*D 2 )+(a2*D+b2*D 2 ) is obtained;

[0012] S60: Obtain the irradiance value Dm at a specific time, and bring the irradiance value Dm into the final relationship to obtain the nitric oxide concentration Km at the specific time.

[0013] In a possible implementation, the first relationship Y1, the second relationship Y2, the third relationship Y3, and the fourth relationship Y4 are modeled and fitted by using MATLAB software or SCILAB software.

[0014] In a possible implementation, the nitric oxide analyzer measures data C11, C12, C13, …, C1n, and the C1=(C11+C12+C13+…+C1n) / n.

[0015] …;

[0016] The nitric oxide analyzer measurement data Cn1, Cn2, Cn3, …, Cnn, Cn = (Cn1+Cn2+Cn3+…+Cnn) / n.

[0017] In a possible implementation, the data Cn1, Cn2, Cn3, …, Cnn are all data within a preset measurement duration.

[0018] In a possible implementation, the chemiluminescence meter measurement data K11, K12, K13, …, K1n, the nitric oxide photolysis rate value K1 is an integral value of the data K11~K1n;

[0019] …;

[0020] The chemiluminescence meter measurement data Kn1, Kn2, Kn3, …, Knn, the nitric oxide photolysis rate value Kn is an integral value of the data Kn1~Knn.

[0021] In a possible implementation, the chemiluminescence meter measurement data K11, K12, K13, …, K1n, the integral value of the data K11~K1n obtains data K1';

[0022] …;

[0023] The chemiluminescence meter measurement data Kn1, Kn2, Kn3, …, Knn, the integral value of the data Kn1~Knn obtains data Kn';

[0024] The nitric oxide photolysis rate value K1 = (K1'+K2'+K3'+K4'+K5') / 5, …, the nitric oxide photolysis rate value Kn = (K(n-4)' + K(n-3)' + K(n-2)' + K(n-1)' + Kn') / 5.

[0025] In a possible implementation, the data Kn1, Kn2, Kn3, …, Knn are all data within a preset measurement duration.

[0026] In a possible implementation, the irradiance data, the nitric oxide photolysis rate value data, and the nitric oxide concentration data are classified according to weather 1, weather 2, …, weather n, and the steps S10-S50 are repeated to obtain the final relationship K corresponding to the weather 1, the weather 2, …, the weather n, respectively.

[0027] In a possible implementation, the irradiance monitoring device, the chemical luminometer and the nitrogen trioxide analyzer are powered on to repeatedly perform the steps S10-S50 to obtain a relationship K', and a floating value W is obtained according to a formula W=K' / K, and finally the nitrogen trioxide concentration = the nitrogen trioxide concentration Km*W at a specific time measured in the step S60.

[0028] In a possible implementation, the step S40 further includes:

[0029] According to a formula

[0030]

[0031]

[0032] to obtain a weight M1 and a weight M2, and the final relationship K= (a1*D+b1*D )*M1+ *M2.

[0033] In the embodiments of the present application, compared with the prior art, the measurement is performed by an online monitoring technology in the prior art, and the existing online monitoring technology includes a cadmium column reduction-naphthalene ethylenediamine hydrochloride spectrophotometry and a spectrum method (such as a tunable diode laser absorption spectrum (TDLAS) and a light cavity ring-down spectrum technology, that is, two kinds of wet chemical method and spectrum method, wherein the wet chemical method needs to consume chemical reagents, frequently configure reagents, produce acidic waste liquid to pollute the environment, and is inconvenient to maintain the equipment; although the spectrum method does not need reagents, the optical remote instrument is easily affected by particulate matter, humidity and interference substances and the like. The estimation method for the nitrogen trioxide concentration in the atmosphere provided in the present application establishes a final relationship K= (a1*D+b1*D 2 ) + (a2*D+b2*D 2 ) through the steps S10-S50 to obtain the relationship between the irradiance and the nitrogen trioxide concentration, and then estimates the nitrogen trioxide concentration by directly measuring the light irradiance. Although the instruments related to the wet chemical method and the spectrum method are used in the process of establishing the final relationship K, the instruments are only used for a short time in the process of establishing the final relationship K, and the estimation of the nitrogen trioxide concentration in the subsequent sequence still relies on the acquisition of the irradiance value, so compared with the existing spectrum method or wet chemical method, the measurement data is stable, and the pollution to the environment is reduced. In the process of establishing the final relationship K, two photolysis modes of the nitrogen trioxide in the atmosphere are considered, the irradiance corresponding to different wavelength ranges is brought into the relationship K=a1*D+b1*D 2 and the relationship K=a2*D+b2*D 2, respectively, and then the concentration of nitrogen trioxide is obtained by taking the sum of the concentrations of nitrogen dioxide and nitrogen monoxide, so that the accuracy of the estimated concentration of nitrogen trioxide is improved. DETAILED DESCRIPTION

[0034] In order to make the technical problems to be solved by the present application, technical solutions and beneficial effects clearer, the following embodiments are used to further illustrate the present application. It should be understood that the specific embodiments described herein are only used to explain the present application, and are not used to limit the present application.

[0035] The method for estimating the concentration of nitrogen trioxide in the atmosphere provided by the present application is described. The method for estimating the concentration of nitrogen trioxide in the atmosphere comprises the following steps:

[0036] S10: The irradiance D1, D2, D3, …, Dn corresponding to the wavelength bands λ1, λ2, …, λn is measured by using a radiance monitoring device;

[0037] S20: The nitrogen trioxide photolysis rate values K1, K2, K3, …, Kn corresponding to the irradiance D1, D2, D3, …, Dn are measured by using a chemical luminometer;

[0038] S30: The irradiance D1~Dx corresponding to the wavelength bands λ1~λx is selected, and the first relationship Y1 between D1~Dx and the corresponding nitrogen trioxide photolysis rate K1~Kx is established;

[0039] The irradiance Dx~Dn corresponding to the wavelength bands λx~λn is selected, and the second relationship Y2 between Dx~Dn and the corresponding nitrogen trioxide photolysis rate Kx~Kn is established, wherein 1

[0040] S40: The nitrogen trioxide concentrations C1, C2, C3, …, Cn corresponding to the nitrogen trioxide photolysis rates K1, K2, K3, …, Kn are measured by using a nitrogen trioxide analyzer, the third relationship Y3 between the nitrogen trioxide concentrations C1~Cx and the corresponding nitrogen trioxide photolysis rate values K1~Kx is established, and the fourth relationship Y4 between the nitrogen trioxide concentrations Cx~Cn and the corresponding nitrogen trioxide photolysis rate values Kx~Kn is established;

[0041] S50: The first relationship Y1 and the third relationship Y3 are converted with each other to obtain the relationship K=a1*D+b1*D corresponding to the wavelength bands λ1~λx 2 The second relationship Y2 and the fourth relationship Y4 are converted with each other to obtain the relationship K=a2*D+b2*D corresponding to the wavelength bands λx~λn 2 , wherein a1, a2, b1, b2 are constants, and the final relationship K=(a1*D+b1*D 2)+ (a2*D+b2*D 2 );

[0042] S60: Obtain the irradiance value Dm at a specific time, and bring the irradiance value Dm into the final relationship formula to obtain the concentration Km of nitrogen trioxide at the specific time.

[0043] It should be noted that, since the concentration of nitrogen trioxide and irradiance are not directly related, direct establishment will result in a large error of the relationship formula; the embodiment first establishes the relationship between irradiance and the photolysis rate of nitrogen trioxide, and then establishes the relationship between the photolysis rate of nitrogen trioxide and the concentration of nitrogen trioxide, thereby improving the accuracy of the relationship formula between irradiance and the concentration of nitrogen trioxide, and further improving the estimation accuracy of the concentration of nitrogen trioxide.

[0044] Compared with the prior art, the estimation method of the concentration of nitrogen trioxide in the atmosphere provided in the embodiment is measured by online monitoring technology. The existing online monitoring technology includes cadmium column reduction-naphthalene ethylenediamine hydrochloride spectrophotometry and spectroscopy (such as tunable diode laser absorption spectroscopy (TDLAS) and optical cavity ring-down spectroscopy technology, i.e., two kinds of wet chemical method and spectroscopy method, wherein the wet chemical method needs to consume chemical reagents, frequently configure reagents, produce acidic waste liquid to pollute the environment, and is inconvenient to maintain the equipment; although the spectroscopy method does not need reagents, the optical remote instrument is easily affected by particulate matter, humidity and interference substances. The estimation method of the concentration of nitrogen trioxide in the atmosphere provided in the application establishes the final relationship formula K=(a1*D+b1*D 2 )+ (a2*D+b2*D 2 ) through steps S10-S50 to obtain the relationship between irradiance and the concentration of nitrogen trioxide, and further estimates the concentration of nitrogen trioxide by directly measuring the light irradiance. Although the instruments related to the wet chemical method and the spectroscopy method are used in the process of establishing the final relationship formula K, such instruments are only used for a short period in the process of establishing the final relationship formula K, and the subsequent estimation of the concentration of nitrogen trioxide still relies on the acquisition of the irradiance value, so compared with the existing spectroscopy method or wet chemical method, the measurement data is stable, and the pollution to the environment is reduced; in the process of establishing the final relationship formula K, two photolysis modes of nitrogen trioxide in the atmosphere are considered, the irradiance corresponding to different wavelength ranges is brought into the relationship formula K=a1*D+b1*D 2 and the relationship formula K=a2*D+b2*D 2 , respectively, to obtain the concentration values of nitrogen dioxide and nitrogen monoxide, and then the sum of the concentration values of nitrogen monoxide and nitrogen dioxide is obtained to obtain the concentration of nitrogen trioxide, thereby improving the accuracy of the estimated concentration of nitrogen trioxide.

[0045] In some embodiments, the first relationship Y1, the second relationship Y2, the third relationship Y3 and the fourth relationship Y4 are modeled and fitted by using MATLAB software or SCILAB software. By inputting the dot array information into the modeling software, the modeling software can automatically generate a linear relationship curve. For example, the radiometric monitoring device and the chemical luminometer are both connected to a computer, and the measurement data of the radiometric monitoring device and the chemical luminometer form dot array information (D1, K1), (D2, K2), …, (Dn, Kn) in the computer, and the first relationship Y1 and the second relationship Y2 are automatically generated. The third relationship Y3 and the fourth relationship Y4 are the same.

[0046] It should be noted that the first relationship Y1, the second relationship Y2, the third relationship Y3 and the fourth relationship Y4 are all monomial quadratic equations.

[0047] For example, when using MATLAB, a function curve relationship graph between irradiance and nitric oxide photolysis rate is mainly drawn, and when actually measuring, the nitric oxide photolysis rate value can be directly obtained by bringing the measured irradiance value into the function curve relationship graph, thereby reducing the labor intensity, and the curve relationship graph drawn by the MATLAB software or the SCILAB software is more accurate, and the measurement accuracy of the nitric oxide concentration is improved.

[0048] In some embodiments, the nitric oxide analyzer measurement data C11, C12, C13, …, C1n, C1= (C11+C12+C13+…+C1n) / n.

[0049] …

[0050] The nitric oxide analyzer measurement data Cn1, Cn2, Cn3, …, Cnn, Cn= (Cn1+Cn2+Cn3+…+Cnn) / n.

[0051] That is, each value of the nitric oxide concentration C1-Cn is an average value calculated by multiple data, and each nitric oxide concentration value is calculated by the average value of multiple data, which not only reduces the calculation workload in the early stage, but also ensures the accuracy of the data.

[0052] It should be noted that the number of data included in C1-Cn can be the same or different; however, when the third relationship Y3 and the fourth relationship Y4 are established between the nitric oxide concentration C and the nitric oxide photolysis rate K, the number of data corresponding to C1 and K1 is the same, and K1 is also obtained by using the same average value as C1, and so on, and Cn and Kn are the same.

[0053] In some embodiments, the data Cn1, Cn2, Cn3, …, Cnn is all data within a preset measurement time length. Optionally, the preset time length is 5 minutes, and each of C1-Cn is an average value of all data obtained within 5 minutes. By dividing the data by time period and then calculating the average value, the data can be conveniently sorted by time period, which is conducive to the orderliness of calculation and facilitates review.

[0054] In some embodiments, the chemiluminometer measures data K11, K12, K13, …, K1n, and the nitrogen trioxide photolysis rate value K1 is an integral value of the data K11-K1n.

[0055] …;

[0056] The chemiluminometer measures data Kn1, Kn2, Kn3, …, Knn, and the nitrogen trioxide photolysis rate value Kn is an integral value of the data Kn1-Knn.

[0057] That is, each value of the nitrogen trioxide photolysis rate K1-Kn is an integral value calculated from multiple data. Since the nitrogen trioxide photolysis rate gradually accumulates over a period of time, calculating each nitrogen trioxide photolysis rate from the integral value of multiple data not only reduces the calculation workload in the early stage, but also ensures the accuracy of the data.

[0058] It should be noted that the number of data included in K1-Kn can be the same or different.

[0059] In some embodiments, the chemiluminometer measures data K11, K12, K13, …, K1n, and the integral value of the data K11-K1n obtains data K1’;

[0060] …;

[0061] The chemiluminometer measures data Kn1, Kn2, Kn3, …, Knn, and the integral value of the data Kn1-Knn obtains data Kn’;

[0062] The nitrogen trioxide photolysis rate value K1=(K1’+K2’+K3’+K4’+K5’) / 5, …, the nitrogen trioxide photolysis rate value Kn=(K(n-4)’+K(n-3)’+K(n-2)’+K(n-1)’+Kn’) / 5.

[0063] Optionally, the integral value is taken for each 60 seconds of measured data, the average value is taken for every 5 integral values, and finally the nitrogen trioxide photolysis rate value used for calculation is obtained, which reduces the calculation of data, further reduces the calculation workload, reduces the dot matrix information of the dot matrix data, ensures the calculation accuracy while reducing the labor intensity.

[0064] It should be noted that the embodiment calculates the average value for every 5 integral values, and in the actual calculation process, it can also be 4 times, 6 times, etc.

[0065] In some embodiments, the data Kn1, Kn2, Kn3, …, Knn is all data within a preset measurement time length. Optionally, the preset time length is 60 seconds, and each of K1~Kn is an integral value of all data within 60 seconds, or an average value of five 60-second integral values. By dividing the data by time period to calculate the integral value or the average value of five integral values, the data can be conveniently arranged in sequence for subsequent time period, which is conducive to the orderliness and convenience of calculation.

[0066] In some embodiments, the irradiance data, the nitric oxide trioxide photolysis rate value data, and the nitric oxide trioxide concentration data are classified according to weather 1, weather 2, …, weather n, and the steps S10~S50 are repeated to obtain the final relationship K corresponding to weather 1, weather 2, …, weather n, respectively. According to the specific weather when the final relationship K is constructed, the final relationship K under different weather conditions can be constructed, such as sunny, cloudy, rainy, etc. According to different weather conditions, the corresponding final relationship K is selected according to the weather at the time of actual measurement, which is conducive to the accuracy of nitric oxide trioxide concentration measurement under different weather conditions.

[0067] In some embodiments, the irradiance monitoring device, the chemical luminometer, and the nitric oxide trioxide analyzer are turned on, and the steps S10~S50 are repeated to obtain the relationship K’, the floating value W is obtained according to the formula W=K’ / K, and the final nitric oxide trioxide concentration=S60 step measured nitric oxide trioxide concentration Km*W at a specific time. Here, K’ and K both refer to the final relationship K=(a1*D+b1*D 2 ) + (a2*D+b2*D 2 ), but K’ and K are not completely the same in , , , The establishment of the final relationship K’ is repeated in the initial measurement of the instrument (here, the instrument includes the irradiance monitoring device, the chemical luminometer, and the nitric oxide trioxide analyzer, and the instrument hereinafter is the same as here). Since the instrument is just turned on at this time, the data after multiple subsequent measurements are compared to obtain the floating value W, and then the nitric oxide trioxide concentration value calculated by the final relationship K is multiplied by the floating value W to obtain the final nitric oxide trioxide concentration, which effectively improves the accuracy of the nitric oxide trioxide concentration estimation.

[0068] In some embodiments,

[0069] The step S40 further includes:

[0070] According to the formula

[0071]

[0072]

[0073] The weight M1 and the weight M2 are obtained, and the final relationship K= (D1+D2+D3+…+D311) *M1+ (D231+D232+D233+…+D311) *M2.

[0074] Since the decomposition of nitrogen trioxide in the air is mainly in the form of nitrogen monoxide and nitrogen dioxide, by calculating the weight, the respective weights of nitrogen monoxide and nitrogen dioxide in the air can be obtained, and then the final relationship obtained can be used to obtain the concentration of nitrogen trioxide according to the actual weight. According to different regions, the weight is different, and the weight can be recalculated in the process of constructing the final relationship, so as to ensure the accuracy of the final relationship K obtained by calculation and improve the estimation accuracy of the concentration of nitrogen trioxide.

[0075] As an example of a specific embodiment in the present application:

[0076] (1) The irradiance D1, D2, D3, …, D311 corresponding to each wave band in the wave band range 470~780nm is measured by using a radiance monitoring device.

[0077] (2) The nitrogen trioxide photolysis rate values K1, K2, K3, …, K311 corresponding to the irradiance D1, D2, D3, …, D311 are measured by using a chemical luminometer, and K1, K2, K3, …, K216 corresponding to the wave band range 470~685nm are selected, and K231, K232, K233, …, K311 corresponding to the wave band range 700~780nm are selected.

[0078] (3) The K1, K2, K3, …, K216 corresponding to the wave band range 470~685nm and the D1, D2, D3, …, D216 are input into the MATLAB software to form a first relationship Y1 (a quadratic equation of one variable); and the K231, K232, K233, …, K311 corresponding to the wave band range 700~780nm and the D231, D232, D233, …, D311 are input into the MATLAB software to form a second relationship Y2 (a quadratic equation of one variable).

[0079] ​​(4) using the nitric oxide analyzer to measure the nitric oxide concentration C1, C2, C3, …, C311 corresponding to K1, K2, K3, …, K311 respectively, and establishing a third relationship Y3 between the nitric oxide concentration C1-C216 and the nitric oxide photolysis rate K1-K216, and establishing a fourth relationship Y4 between the nitric oxide concentration C231-C311 and the nitric oxide photolysis rate K231-K311;

[0080] (5) converting the first relationship and the third relationship to each other to obtain the relationship K=a1*D+b1*D 2 ;

[0081] Converting the second relationship and the fourth relationship to each other to obtain the relationship K=a2*D+b2*D 2 ;

[0082] According to the formula:

[0083]

[0084]

[0085] The weights M1 and M2 are calculated;

[0086] The final relationship K=(a1*D+b1*D 2 )*M1+(a2*D+b2*D 2 )*M2;

[0087] (6) The steps (1)-(5) are repeated to obtain the relationship Y3' when the irradiance monitoring device, the chemical luminometer, and the nitric oxide analyzer are turned on, and the floating value W is obtained according to the formula W=Y3' / Y3;

[0088] (7) The irradiance value Dm at a specific time is obtained, and the irradiance value Dm is brought into the final relationship K=(a1*D+b1*D 2 )*M1+(a2*D+b2*D 2 )*M2 to calculate the nitric oxide concentration at the specific time, and the final nitric oxide concentration is obtained by multiplying the nitric oxide concentration at the specific time by the floating value W.

[0089] It should be noted that all the data involving n in the embodiments of the present application, wherein n represents an unknown integer greater than 1, the n value corresponding to each step may be the same or different.

[0090] The above only describes the preferred embodiments of the present application and is not intended to limit the present application, and any modifications, equivalent replacements and improvements made within the spirit and principles of the present application shall be included in the protection scope of the present application.

Claims

1. A method of estimating the concentration of nitrogen trioxide in the atmosphere, characterized by, Comprising the following steps: S10: Measure the irradiance D1, D2, D3, …, Dn corresponding to the wave band λ1, λ2, …, λn respectively by using the radiance monitoring device; S20: Measure the nitric oxide photolysis rate value K1, K2, K3, …, Kn corresponding to the irradiance D1, D2, D3, …, Dn respectively by using the chemical photometer; S30: Select the irradiance D1~Dx corresponding to the wave band λ1~λx, and establish a first relationship Y1 between the D1~Dx and the corresponding nitric oxide photolysis rate K1~Kx; Select the irradiance Dx~Dn corresponding to the wave band λx~λn, and establish a second relationship Y2 between the Dx~Dn and the corresponding nitric oxide photolysis rate Kx~Kn, wherein 1 S40: Measure the nitric oxide concentration C1, C2, C3, …, Cn corresponding to the nitric oxide photolysis rate K1, K2, K3, …, Kn respectively by using the nitric oxide analyzer, establish a third relationship Y3 between the nitric oxide concentration C1~Cx and the corresponding nitric oxide photolysis rate value K1~Kx, and establish a fourth relationship Y4 between the nitric oxide concentration Cx~Cn and the corresponding nitric oxide photolysis rate value Kx~Kn; S50: converting the first relationship Y1 and the third relationship Y3 to each other to obtain a relationship K=a1*D+b1*D corresponding to the wave band λ1~λx 2 , converting the second relationship Y2 and the fourth relationship Y4 to each other to obtain a relationship K=a2*D+b2*D corresponding to the wave band λx~λn 2 , wherein a1, a2, b1, b2 are constants, and the final relationship K=(a1*D+b1*D 2 )+(a2*D+b2*D 2 ) S60: Obtain the irradiance value Dm at a specific time, and calculate the nitric oxide concentration Km at the specific time by bringing the irradiance value Dm into the final relationship.

2. The method of estimating the concentration of nitrogen trioxide in the atmosphere according to claim 1, characterized by, The first relationship Y1, the second relationship Y2, the third relationship Y3, and the fourth relationship Y4 are modeled and fitted by using MATLAB software or SCILAB software.

3. The method of estimating concentration of nitrogen trioxide in the atmosphere according to claim 1, wherein The nitric oxide analyzer measurement data C11, C12, C13, …, C1n, and C1=(C11+C12+C13+…+C1n) / n; ……； The nitric oxide analyzer measurement data Cn1, Cn2, Cn3, …, Cnn, and Cn=(Cn1+Cn2+Cn3+…+Cnn) / n.

4. The method of estimating concentration of nitrogen trioxide in the atmosphere according to claim 3, wherein The data Cn1, Cn2, Cn3, …, Cnn are all data within a preset measurement duration.

5. The method of estimating concentration of nitrogen trioxide in the atmosphere according to claim 3, wherein The chemical photometer measurement data K11, K12, K13, …, K1n, and the nitric oxide photolysis rate value K1 is the integral value of the data K11~K1n; ……； The chemical photometer measurement data Kn1, Kn2, Kn3, …, Knn, and the nitric oxide photolysis rate value Kn is the integral value of the data Kn1~Knn.

6. The method of estimating concentration of nitrogen trioxide in the atmosphere according to claim 3, wherein The chemical photometer measurement data K11, K12, K13, …, K1n, and the integral value of the data K11~K1n is data K1'; ……； The chemical photometer measurement data Kn1, Kn2, Kn3, …, Knn, and the integral value of the data Kn1~Knn is data Kn'; The nitric oxide photolysis rate value K1=(K1'+K2'+K3'+K4'+K5') / 5, and the nitric oxide photolysis rate value Kn=(K(n-4')+K(n-3')+K(n-2')+K(n-1')+Kn') / 5.

7. The method of estimating the concentration of nitrogen trioxide in the atmosphere according to claim 5 or 6, characterized by, The data Kn1, Kn2, Kn3, …, Knn are all data within a preset measurement duration.

8. The method of estimating concentration of nitrogen trioxide in the atmosphere according to claim 1, wherein The irradiance data, the nitric oxide photolysis rate value data, and the nitric oxide concentration data are classified according to weather 1, weather 2, …, weather n, and the steps S10-S50 are repeated to obtain the final relationship K corresponding to the weather 1, the weather 2, …, and the weather n, respectively.

9. The method of estimating concentration of nitrogen trioxide in the atmosphere according to claim 1, wherein The irradiance monitoring device, the chemical luminometer, and the nitric oxide analyzer are turned on, the steps S10-S50 are repeated to obtain the relationship K', the floating value W is obtained according to the formula W=K' / K, and the final nitric oxide concentration=the nitric oxide concentration Km measured at a specific time in the step S60*W.

10. The method of estimating concentration of nitrogen trioxide in the atmosphere according to claim 1, wherein The step S40 further comprises: According to the formula obtaining a weight M1 and a weight M2, the final relationship K = ( ) * M1 + * M2.​

Images