Method, system, device and medium for quantifying ozone source contributions within a boundary layer

Abstract

The application discloses a boundary layer ozone source contribution quantification method, system, device and medium, and relates to the field of ozone source analysis; vertical observation data of meteorological parameters and target species at different altitudes in different atmospheric layer structures at a target moment are acquired; vertical integrated column concentrations of the target species in different atmospheric layer structures at the target moment are calculated according to the vertical observation data of the meteorological parameters and the target species at different altitudes in different atmospheric layer structures at the target moment by using a vertical integration column concentration calculation formula; the entrainment contribution and the chemical production contribution of the residual layer at night are calculated according to the vertical integrated column concentrations of the target species in different atmospheric layer structures at the target moment; and the ozone source change in the boundary layer is determined according to the entrainment contribution and the chemical production contribution of the residual layer at night. The ozone source contribution quantification in the boundary layer is calculated by using the formula, the titration influence of NO on ozone is considered, and the accuracy of the ozone source contribution quantification in the boundary layer is improved.

Description

Technical Field

[0001] This invention relates to the field of ozone source analysis, and in particular to methods, systems, equipment and media for quantifying the contribution of ozone sources within the boundary layer. Background Technology

[0002] Tropospheric ozone is a significant environmental problem worldwide due to its adverse effects on human health, ecosystems, and climate change.

[0003] Tropospheric ozone is primarily a product of the photochemical reaction between nitrogen oxides (NOx = NO + NO2) and volatile organic compounds (VOCs). Within the atmospheric boundary layer, in addition to local photochemical reactions, external transport (including horizontal and vertical transport) is also a significant contributor to changes in total ozone levels. Due to the complexity of ozone formation and the wide range of its sources, the quantitative analysis of ozone source contributions is a major research focus and a pressing scientific challenge in the field of atmospheric pollution control.

[0004] To date, several methods have been proposed to quantify the contribution of ozone sources within the boundary layer, including the OSAT / APCA method in the CAMx model and the ISAM method in the CMAQ model. While these numerical model-based methods can directly provide quantitative assessment results, the lack of necessary observational data constraints and validation often results in high uncertainty, especially when used to resolve the vertical transport contribution of ozone within the boundary layer, where the uncertainty is even more difficult to control. Furthermore, some researchers have quantified the contribution of ozone sources within the boundary layer based on vertical observational data; however, this often ignores the influence of chemical reactions during the vertical mixing process of ozone (such as the titration loss of ozone by NO), leading to serious biases in the quantitative assessment results of the ozone source contribution within the boundary layer. Overall, a reliable method for resolving the contribution of ozone sources within the boundary layer based on direct vertical observational data still lacks a reliable approach.

[0005] Ozone pollution is becoming increasingly serious in many cities and regions, making it crucial to clarify the main sources of ozone within the boundary layer. Summary of the Invention

[0006] The purpose of this invention is to provide a method, system, device, and medium for quantifying the contribution of ozone sources within the boundary layer. The method calculates the contribution of ozone sources within the boundary layer using a formula and takes into account the titration effect of NO on ozone, thereby improving the accuracy of quantifying the contribution of ozone sources within the boundary layer.

[0007] To achieve the above objectives, the present invention provides the following solution:

[0008] A method for quantifying the contribution of ozone sources within the boundary layer includes:

[0009] The system acquires vertical observation data of meteorological parameters and target species at different altitudes within different atmospheric strata at a target time. The meteorological parameters include temperature and air pressure. The target species include ozone, nitric oxide, and nitrogen dioxide. The atmospheric strata include the boundary layer and the residual layer. The target time includes multiple times during the day and multiple times at night.

[0010] The vertical integral column concentration of the target species in different atmospheric strata at the target time was calculated using the vertical integral column concentration calculation formula based on meteorological parameters at different altitudes in different atmospheric strata at the target time and vertical observation data of the target species.

[0011] The nighttime residual layer entrainment contribution and chemical generation contribution were calculated based on the vertical integral column concentration of the target species in different atmospheric strata within the target time.

[0012] The sources of ozone variation within the boundary layer are determined based on the contributions from the nighttime residual layer and the contributions from chemical generation.

[0013] Optionally, vertical observation data of meteorological parameters and target species at different altitudes within different atmospheric strata at the target time can be acquired, specifically including:

[0014] Use at least one of airships, drones, and lidar to acquire vertical observation data of meteorological parameters and target species at different altitudes in different atmospheric strata at the target time.

[0015] Optionally, the vertical integral column concentration of the target species in different atmospheric strata within the target time is calculated using the vertical integral column concentration calculation formula based on meteorological parameters at different altitudes in different atmospheric strata within the target time and vertical observation data of the target species. Specifically, this includes:

[0016] The ideal gas law is used to calculate the gas molar volume at different altitudes in each atmospheric layer at the target time based on meteorological parameters at different altitudes in different atmospheric layers at the target time.

[0017] The vertical integral column concentration calculation formula is used to calculate the vertical integral column concentration of the target species in each atmospheric stratum at the target time based on the volume mixing ratio and gas molar volume of the target species at different altitudes in each atmospheric stratum.

[0018] Optionally, the formula for calculating the concentration of the vertical integration column is:

[0019] ;

[0020] in, This represents the vertical integral column concentration of target species i in the atmospheric stratification layer c at time t, in units of: ; This represents the volume mixing ratio of target species i at altitude h in the atmospheric stratification layer c at time t, in units of: h1 represents the lowest altitude in the c-layer atmospheric stratification, in cm; h2 represents the highest altitude in the c-layer atmospheric stratification, in cm; N A Represents Avogadro's constant; This represents the molar volume of gas at altitude h within the c-layer atmospheric stratification at time t, in units of: .

[0021] Optionally, the nighttime remnant layer entrainment contribution and chemical generation contribution are calculated based on the vertical integral column concentration of the target species in different atmospheric strata within the target time, specifically including:

[0022] The target species are classified according to preset conditions to obtain total oxidants and nitrogen oxides; the preset conditions are that ozone and nitrogen dioxide are used as total oxidants; and nitric oxide and nitrogen dioxide are used as nitrogen oxides.

[0023] Based on the vertical integral column concentrations of the target species, total oxidants, and nitrogen oxides in different atmospheric strata at the target time, the vertical integral column concentrations of nitrogen oxides and total oxidants in different atmospheric strata at the target time are obtained.

[0024] The contribution of the nighttime residual layer entrainment was calculated based on the vertical integral column concentration change of total oxidant in different atmospheric strata within the target time.

[0025] The ground emission contribution is calculated based on the vertical integral column concentration changes of nitrogen oxides in different atmospheric strata during the target time.

[0026] The chemical generation contribution was calculated based on the vertical integral column concentrations of nitrogen oxides in different atmospheric strata at the target time, the vertical integral column concentrations of total oxidants in different atmospheric strata at the target time, and the contribution from ground emissions.

[0027] Optionally, the formula for calculating the contribution of the nighttime residual layer entrainment is:

[0028] ;

[0029] The formula for calculating the contribution of surface emissions is as follows:

[0030] ;

[0031] The formula for calculating the contribution of chemical generation is as follows:

[0032] ;

[0033] in, Indicates the contribution of the nighttime residual layer entrainment, in units: ; This represents the average vertical integral column concentration of total oxidant in the residual layer during the night, in units of: ; This represents the vertical integral column concentration of total oxidant in the residual layer at time t during the day, in units of: ; Indicates the contribution of ground-level emissions, in units of: ; This represents the sum of the vertical integral column concentrations of nitrogen oxides in the residual layer and boundary layer at time t during the day, in units of: ; This represents the sum of the vertically integrated column concentrations of nitrogen oxides in the residual layer and boundary layer during the night, in units of: ER represents the percentage of nitrogen dioxide emissions in primary nitrogen oxide emissions; Indicates chemical contribution, unit: ; This represents the sum of the vertical integral column concentrations of total oxidant in the boundary layer and residual layer at time t during the day, in units of: ; This represents the sum of the average vertical integral column concentrations of total oxidants in the inner boundary layer and residual layer at night, in units of: .

[0034] A boundary layer ozone source contribution quantification system, wherein the boundary layer ozone source contribution quantification system is applied to the above-described boundary layer ozone source contribution quantification method, and the boundary layer ozone source contribution quantification system comprises:

[0035] The acquisition module is used to acquire vertical observation data of meteorological parameters and target species at different altitudes in different atmospheric strata within a target time period; the meteorological parameters include temperature and air pressure; the target species include ozone, nitric oxide, and nitrogen dioxide; the atmospheric strata include the boundary layer and the residual layer; the target time includes multiple times during the day and multiple times at night.

[0036] The first calculation module is used to calculate the vertical integral column concentration of the target species in different atmospheric strata within the target time using the vertical integral column concentration calculation formula based on meteorological parameters at different altitudes in different atmospheric strata within the target time and vertical observation data of the target species.

[0037] The second calculation module is used to calculate the nighttime residual layer entrainment contribution and chemical generation contribution based on the vertical integral column concentration of the target species in different atmospheric strata within the target time.

[0038] The third calculation module is used to determine the source of ozone variation in the boundary layer based on the contribution of the nighttime residual layer and the contribution of chemical generation.

[0039] Optionally, the acquisition module specifically includes:

[0040] Use at least one of airships, drones, and lidar to acquire vertical observation data of meteorological parameters and target species at different altitudes in different atmospheric strata at the target time.

[0041] An electronic device includes a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor, when executing the computer program, implements the aforementioned method for quantifying the contribution of ozone sources within the boundary layer.

[0042] A computer-readable storage medium storing a computer program that, when executed, implements the method for quantifying the contribution of ozone sources within the boundary layer as described above.

[0043] According to specific embodiments provided by the present invention, the present invention discloses the following technical effects:

[0044] This invention discloses a method, system, device, and medium for quantifying the contribution of ozone sources within the boundary layer. Based on vertical observation data of meteorological parameters at different altitudes and target species in different atmospheric strata at a target time, the nighttime remnant layer entrainment contribution and chemical formation contribution are calculated using formulas. These two contributions are then added together to determine the contribution of ozone source variations within the boundary layer. NO is included as a target species in this invention, and its titration effect on ozone is taken into account, improving the accuracy of quantifying the ozone source contribution within the boundary layer. Attached Figure Description

[0045] To more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the drawings used in the embodiments will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0046] Figure 1 This is a flowchart of a method for quantifying the contribution of ozone sources within the boundary layer in an embodiment of the present invention;

[0047] Figure 2 This is a schematic diagram showing the boundary layer and residual layer heights determined using simulation data from the WRF v4.1.2 model in this invention.

[0048] Figure 3 This is a schematic diagram of the time series of Ox vertical integral column concentration within the case period calculated in this embodiment of the invention;

[0049] Figure 4This is a scatter plot and a linear fitting result plot of the concentrations of nitric oxide (NO) and nitrogen oxides (NOx) during the time period of the present invention embodiment;

[0050] Figure 5 This is a first schematic diagram illustrating the contribution and percentage of residual layer entrainment, chemical generation, and ground emissions to the increase in Ox vertical integral column concentration in the daytime boundary layer, calculated in this embodiment of the invention.

[0051] Figure 6 This is a second schematic diagram illustrating the contribution and percentage of residual layer entrainment, chemical generation, and ground emissions to the increase in Ox vertical integral column concentration in the daytime boundary layer, calculated in this embodiment of the invention.

[0052] Figure 7 This is a third schematic diagram illustrating the contribution and percentage of residual layer entrainment, chemical generation, and ground emissions to the increase in Ox vertical integral column concentration in the daytime boundary layer, calculated in this embodiment of the invention.

[0053] Figure 8 The fourth schematic diagram shows the contribution and percentage of residual layer entrainment, chemical generation, and ground emissions to the increase of Ox vertical integral column concentration in the daytime boundary layer, calculated in this embodiment of the invention. Detailed Implementation

[0054] The technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0055] The purpose of this invention is to provide a method, system, device, and medium for quantifying the contribution of ozone sources within the boundary layer. The method calculates the contribution of ozone sources within the boundary layer using a formula and takes into account the titration effect of NO on ozone, thereby improving the accuracy of quantifying the contribution of ozone sources within the boundary layer.

[0056] This invention provides a novel method to quantify the main source contribution of ozone within the boundary layer during diurnal atmospheric boundary layer uplift, specifically a method for quantifying the source contribution of ozone within the boundary layer. This invention is based on the changes in the vertical integrated concentration (CIC) of total oxidants (Ox, Ox = O3 + NO2) and nitrogen oxides (NOx, NOx = NO + NO2) in the boundary layer and residual layer.

[0057] The period of daytime atmospheric boundary layer uplift refers to the time range during which the height of the daytime boundary layer continues to increase over time.

[0058] To make the above-mentioned objects, features and advantages of the present invention more apparent and understandable, the present invention will be further described in detail below with reference to the accompanying drawings and specific embodiments.

[0059] Example 1

[0060] like Figure 1 As shown, the present invention provides a method for quantifying the contribution of ozone sources within the boundary layer, comprising:

[0061] Step 101: Obtain vertical observation data of meteorological parameters and target species at different altitudes within different atmospheric strata at the target time. The meteorological parameters include temperature and air pressure; the target species include ozone (O3), nitric oxide (NO), and nitrogen dioxide (NO2); the atmospheric strata include the boundary layer and the residual layer; the target time includes multiple times during the day and multiple times at night.

[0062] Step 102: Using the vertical integral column concentration calculation formula, the vertical integral column concentration of the target species in different atmospheric strata at different altitudes within the target time is calculated based on the meteorological parameters at different altitudes and the vertical observation data of the target species within the target time.

[0063] The vertical integral column concentration (CIC) is calculated using formula (1) based on meteorological parameters at different altitudes and vertical observation data of the target species. It is used to characterize the total amount of species molecules in a certain height space per unit area and serves as a basic parameter for constructing the ozone source contribution quantification method within the boundary layer.

[0064] Step 103: Calculate the nighttime remnant layer entrainment contribution and chemical generation contribution based on the vertical integral column concentration of the target species in different atmospheric strata within the target time. The contribution represents the proportion of ozone concentration change within the boundary layer.

[0065] Step 104: Determine the source of ozone variation within the boundary layer based on the contribution of the nighttime residual layer and the contribution of chemical generation.

[0066] As one embodiment, acquiring vertical observation data of meteorological parameters and target species at different altitudes within different atmospheric strata at a target time specifically includes:

[0067] Use at least one of airships, drones, and lidar to acquire vertical observation data of meteorological parameters and target species at different altitudes in different atmospheric strata at the target time.

[0068] In practice, data acquisition devices can be installed on high-rise buildings or towers to obtain vertical observation data of meteorological parameters and target species at different altitudes in different atmospheric strata at the target time.

[0069] In practice, the vertical distribution data of O3, NO, NO2, temperature and pressure, as well as parameters such as boundary layer height, obtained by using various observation platforms or technologies such as airships, drones, lidar, and numerical models can provide data support for the calculation of Ox (O3+NO2) and NOx (NO+NO2) vertical integral column concentrations.

[0070] As one embodiment, the vertical integral column concentration calculation formula is used to calculate the vertical integral column concentration of the target species in different atmospheric strata at the target time based on meteorological parameters at different altitudes in different atmospheric strata and vertical observation data of the target species at the target time. Specifically, it includes:

[0071] The ideal gas equation of state is used to calculate the gas molar volume at different altitudes in each atmospheric stratum at the target time based on meteorological parameters at different altitudes in different atmospheric strata at the target time.

[0072] The vertical integral column concentration calculation formula is used to calculate the vertical integral column concentration of the target species in each atmospheric stratum at the target time based on the volume mixing ratio and gas molar volume of the target species at different altitudes in each atmospheric stratum.

[0073] As one embodiment, the formula for calculating the concentration of the vertical integration column is:

[0074] (1)

[0075] in, This represents the vertical integral column concentration of target species i in the atmospheric stratification layer c at time t, in units of: ; This represents the volume mixing ratio of target species i at altitude h in the atmospheric stratification layer c at time t, in units of: h1 represents the lowest altitude in the c-layer atmospheric stratification, in cm; h2 represents the highest altitude in the c-layer atmospheric stratification, in cm; N A Represents Avogadro's constant; This represents the molar volume of gas at altitude h within the c-layer atmospheric stratification at time t, in units of: .

[0076] The methods for obtaining height data of different atmospheric strata include, but are not limited to, atmospheric boundary layer height data obtained by simulation using the WRF v4.1.2 model, which can provide boundary layer and residual layer height information for the calculation of vertical integral column concentration of relevant substances.

[0077] In the formula (1) The variable is one that varies with altitude h and can be calculated using the ideal gas law based on the atmospheric temperature (in K) and atmospheric pressure (in Pa) measured at altitude h.

[0078] The species i in the formula (1) includes nitric oxide (NO) and nitrogen dioxide (NO2), etc.

[0079] When the vertical observation data height of species i in formula (1) is less than the maximum height h2 of the atmospheric stratification, in formula (1) The calculation follows these four principles: ① The Ox concentration between the maximum measurement altitude at time t and the highest atmospheric layer height h2 is considered to be the same as the concentration at the maximum measurement altitude at time t; ② The NOx concentration between the maximum measurement altitude at time t and the highest atmospheric layer height h2 can be estimated using a linear fitting method based on observation data within the observation altitude at time t; ③ The concentration of species i within the remnant layer at time t does not change significantly relative to its concentration within the nighttime remnant layer; ④ The concentration distribution of species i between two adjacent observation altitudes is estimated using a linear interpolation method.

[0080] The four principles can be specifically defined as follows: ① The Ox concentration between the maximum measurement height (e.g., 335 m) and the top of the daytime boundary layer is equal to the concentration at the maximum measurement height; ② The NOx concentration between the maximum measurement height and the maximum height of the daytime boundary layer can be estimated using a linear fitting method based on the observation results within the observation height (e.g., 5~335 m); ③ The NOx and Ox concentrations within the daytime residual layer are equal to their concentrations within the nighttime residual layer; ④ The concentration distribution of Ox and NOx between two adjacent observation heights is estimated using a linear interpolation method.

[0081] like Figure 2 As shown, the height of the residual layer is equal to the smaller of the maximum heights of the boundary layer over two consecutive days.

[0082] As one example, the nighttime remnant layer entrainment contribution and chemical generation contribution are calculated based on the vertical integral column concentration of the target species in different atmospheric strata within the target time, specifically including:

[0083] The target species are classified according to preset conditions to obtain total oxidants (Ox) and nitrogen oxides (NOx); the preset conditions are that O3 and NO2 are regarded as Ox; and NO and NO2 are regarded as NOx.

[0084] Based on the vertical integral column concentrations of the target species, Ox, and NOx in different atmospheric strata at the target time, the vertical integral column concentrations of NOx and Ox in different atmospheric strata at the target time are obtained.

[0085] The nighttime residual layer entrainment contribution was calculated based on the vertical integral column concentration changes of Ox in different atmospheric strata during the target time.

[0086] The ground emission contribution was calculated based on the vertical integral column concentration changes of NOx in different atmospheric strata during the target time.

[0087] The chemical formation contribution was calculated based on the vertical integral column concentrations of NOx in different atmospheric strata at the target time, the vertical integral column concentrations of Ox in different atmospheric strata at the target time, and the contribution of ground emissions.

[0088] In practical implementation, formula (2) can be used to calculate the change in the vertical integral column concentration of Ox in the boundary layer at time t during the day relative to its vertical integral column concentration in the boundary layer at night, based on the vertical integral column concentration of Ox in different atmospheric strata at the target time. The unit is: .

[0089] (2)

[0090] in, This represents the change in the vertical integral column concentration of Ox in the boundary layer at time t during the day relative to its vertical integral column concentration in the boundary layer at night; the first term on the right. This represents the vertical integral column concentration of Ox in the boundary layer at time t during the day, calculated using the formula (1) (unit: ), the second item on the right This represents the vertical integral column concentration of Ox in the nighttime boundary layer calculated using the formula (1) (unit: ).

[0091] The left-hand term of formula (2) is the concentration change of the Ox vertical integral column. The main contributing factors are: residual layer entrainment contribution, surface emission contribution, and chemical formation contribution (unit: ).

[0092] The average vertical integral column concentration of the nighttime boundary layer Ox in formula (2) The concentration of Ox in the boundary layer vertical integral column remains relatively stable over a period of time at night (e.g.) Figure 3 The average value (shown between 4:00 and 6:00) is used to represent the average total amount of Ox in the nighttime boundary layer.

[0093] Then, based on formula (3), an accuracy assessment is made of the calculated contributions of ground emissions, nighttime residual layer entrainment, and chemical generation.

[0094] (3)

[0095] As one embodiment, the formula for calculating the contribution of the nighttime residual layer entrainment is:

[0096] (4)

[0097] The formula for calculating the contribution of surface emissions is as follows:

[0098] (5)

[0099] Specifically, the effect of a single ground-based NO2 emission on the Ox vertical integral column concentration change within the boundary layer. . contributions.

[0100] In formula (5), the contribution ratio ER of primary NO2 emissions represents the contribution ratio of primary NO2 emissions from ground-based NOx emission sources to total NOx emissions. This can be used to linearly fit the NO and NOx concentrations observed at ground level during the nighttime period (22:00-6:00), such as... Figure 4 As shown, the emission contribution ratio of NO in the ground emission source (1-ER) is obtained, and then the primary NO2 emission contribution ratio ER in the ground emission source is calculated.

[0101] Chemical contributions include ozone photochemical reaction formation, chemical removal of ozone and NO2, and the effect of dry ground deposition on the Ox vertical integral column concentration change in the boundary layer of equation (2). The net contribution is specifically the difference between the boundary layer and residual layer Ox vertical integral column concentrations during the day and night periods, calculated using formula (1), minus the surface emission contribution in formula (5). The chemical contribution can be calculated using formula (6), and the formula for calculating the chemical contribution is:

[0102] (6)

[0103] The chemical generation contribution represents the net contribution of diurnal boundary layer uplift, diurnal photochemical reaction generation, chemical removal, and ground subsidence to the change in the vertical integral column concentration of Ox in the boundary layer.

[0104] The sum of the residual layer entrainment contribution, the ground emission contribution, and the chemical generation contribution at a certain moment during the day is equal to the total change in the Ox vertical integral column concentration in the daytime boundary layer at that moment. That is, the sum of the right-hand side terms of the formulas (4), (5), and (6) is equal to the calculated value of the right-hand side term of the formula (2).

[0105] Using the formulas (4), (5), and (6), the contributions and percentages of the residual layer entrainment contribution, surface emission contribution, and chemical generation contribution to the change in Ox vertical integral column concentration within the boundary layer during diurnal boundary layer uplift can be calculated, such as... Figure 5 , Figure 6 , Figure 7 and Figure 8 As shown.

[0106] in, This represents the contribution of the nighttime residual layer entrainment (the vertical integral column concentration of total oxidant in the residual layer entrainment at time t during the night), in units of: ; This represents the average vertical integral column concentration of total oxidant in the residual layer during the night; This represents the vertical integral column concentration of total oxidant in the residual layer at time t during the day, in units of: ; This represents the contribution of ground emissions (the vertical integral column concentration of total oxidant in ground emissions at time t), in units of: ; This represents the sum of the vertical integral column concentrations of nitrogen oxides in the residual layer and boundary layer at time t during the day, in units of: ; This represents the sum of the vertically integrated column concentrations of nitrogen oxides in the residual layer and boundary layer during the night, in units of: ER represents the percentage of nitrogen dioxide emissions in primary nitrogen oxide emissions; The contribution of chemical generation (concentration of total oxidant in the vertical integral column during time t) is expressed in units of: ; This represents the sum of the vertical integral column concentrations of total oxidant in the boundary layer and residual layer at time t during the day, in units of: ; This represents the sum of the average vertical integral column concentrations of total oxidants in the inner boundary layer and residual layer at night, in units of: .

[0107] The sum of the calculation results of the right-hand side terms in formulas (4), (5), and (6) is equal to the calculation result of the right-hand side of formula (3), thus completing the calculation loop. According to formulas (1) to (6), the change in the Ox vertical integral column concentration in the boundary layer at time t during the day relative to its vertical integral column concentration in the boundary layer at night can be obtained. ) and the contribution of the three factors ( , and ).

[0108] The Ox vertical integral column concentration variation method takes into account the titration effect of NO on ozone, and after eliminating the contribution of ground emissions, the chemical formation contribution of the Ox column concentration variation in the boundary layer is... and residual layer entrainment contribution It can also be regarded as its contribution to the change of ozone column concentration in the boundary layer. It can be calculated using the formulas (4) and (6) respectively, thus achieving the purpose of quantifying the contribution of ozone source in the boundary layer.

[0109] Example 2

[0110] A boundary layer ozone source contribution quantification system is provided, wherein the boundary layer ozone source contribution quantification system is applied to the boundary layer ozone source contribution quantification method described in Example 1, and the boundary layer ozone source contribution quantification system includes:

[0111] The acquisition module is used to acquire vertical observation data of meteorological parameters and target species at different altitudes in different atmospheric strata within a target time period; the meteorological parameters include temperature and air pressure; the target species include ozone, nitric oxide, and nitrogen dioxide; the atmospheric strata include the boundary layer and the residual layer; and the target time includes multiple times during the day and multiple times at night.

[0112] The first calculation module is used to calculate the vertical integral column concentration of the target species in different atmospheric strata within the target time using the vertical integral column concentration calculation formula based on meteorological parameters at different altitudes in different atmospheric strata within the target time and vertical observation data of the target species.

[0113] The second calculation module is used to calculate the nighttime residual layer entrainment contribution and chemical generation contribution based on the vertical integral column concentration of the target species in different atmospheric strata at the target time.

[0114] The third calculation module is used to determine the source of ozone variation in the boundary layer based on the contribution of the nighttime residual layer and the contribution of chemical generation.

[0115] As one embodiment, the acquisition module specifically includes:

[0116] Use at least one of airships, drones, and lidar to acquire vertical observation data of meteorological parameters and target species at different altitudes in different atmospheric strata at the target time.

[0117] An electronic device includes a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor, when executing the computer program, implements the boundary layer ozone source contribution quantification method described in Embodiment 1.

[0118] A computer-readable storage medium storing a computer program that, when executed, implements the method for quantifying the contribution of ozone sources within the boundary layer as described in Embodiment 1.

[0119] The method of this invention can be applied to vertical profile data of Ox and NOx obtained by various vertical observation techniques or model simulations, so as to quantify the contribution of residual layer entrainment and chemical generation to ozone growth in the boundary layer. It is applicable to the need for source apportionment of ozone in the daytime boundary layer under different scenarios in the fields of atmospheric chemistry and air pollution control.

[0120] The various embodiments in this specification are described in a progressive manner, with each embodiment focusing on its differences from other embodiments. Similar or identical parts between embodiments can be referred to interchangeably. For the systems disclosed in the embodiments, since they correspond to the methods disclosed in the embodiments, the descriptions are relatively simple; relevant parts can be referred to the method section.

[0121] This document uses specific examples to illustrate the principles and implementation methods of the present invention. The descriptions of the above embodiments are only for the purpose of helping to understand the method and core ideas of the present invention. Furthermore, those skilled in the art will recognize that, based on the ideas of the present invention, there will be changes in the specific implementation methods and application scope. Therefore, the content of this specification should not be construed as a limitation of the present invention.

Claims

1. A method for quantifying the contribution of ozone sources within the boundary layer, characterized in that, The method for quantifying the contribution of ozone sources within the boundary layer includes: The system acquires vertical observation data of meteorological parameters and target species at different altitudes within different atmospheric strata at a target time. The meteorological parameters include temperature and air pressure. The target species include ozone, nitric oxide, and nitrogen dioxide. The atmospheric strata include the boundary layer and the residual layer. The target time includes multiple times during the day and multiple times at night. The vertical integral column concentration of the target species in different atmospheric strata at the target time is calculated using the vertical integral column concentration calculation formula based on meteorological parameters at different altitudes in different atmospheric strata within the target time and vertical observation data of the target species. The vertical integral column concentration calculation formula is as follows: ； in, This represents the vertical integral column concentration of target species i in the atmospheric stratification layer c at time t, in molecules cm⁻¹. -2 ; The volume mixing ratio of target species i at altitude h in the atmospheric stratification layer c at time t is dimensionless; h1 represents the lowest altitude in the atmospheric stratification layer c, in cm; h2 represents the highest altitude in the atmospheric stratification layer c, in cm; N A Represents Avogadro's constant; This represents the molar volume of gas at altitude h within the c-layer atmospheric stratification at time t, in cm³. 3 mol -1 ; The nighttime remnant layer entrainment contribution and chemical generation contribution are calculated based on the vertical integral column concentrations of target species in different atmospheric strata within the target time period. Specifically, this includes: classifying target species according to preset conditions to obtain total oxidants and nitrogen oxides; the preset conditions are ozone and nitrogen dioxide as total oxidants; and nitric oxide and nitrogen dioxide as nitrogen oxides; obtaining the vertical integral column concentrations of nitrogen oxides and total oxidants in different atmospheric strata within the target time period based on the vertical integral column concentrations of target species, total oxidants, and nitrogen oxides in different atmospheric strata within the target time period; calculating the nighttime remnant layer entrainment contribution based on the changes in the vertical integral column concentrations of total oxidants in different atmospheric strata within the target time period; calculating the ground emission contribution based on the changes in the vertical integral column concentrations of nitrogen oxides in different atmospheric strata within the target time period; and calculating the chemical generation contribution based on the vertical integral column concentrations of nitrogen oxides, total oxidants, and ground emission contributions in different atmospheric strata within the target time period. The sources of ozone variation within the boundary layer are determined based on the contributions from the nighttime residual layer and the contributions from chemical generation.

2. The method for quantifying the contribution of ozone sources within the boundary layer according to claim 1, characterized in that, Acquire vertical observation data of meteorological parameters and target species at different altitudes within different atmospheric strata at the target time, specifically including: Use at least one of airships, drones, and lidar to acquire vertical observation data of meteorological parameters and target species at different altitudes in different atmospheric strata at the target time.

3. The method for quantifying the contribution of ozone sources within the boundary layer according to claim 1, characterized in that, The vertical integral column concentration calculation formula is used to calculate the vertical integral column concentration of the target species in different atmospheric strata at the target time, based on meteorological parameters at different altitudes in different atmospheric strata and vertical observation data of the target species. Specifically, it includes: The ideal gas law is used to calculate the gas molar volume at different altitudes in each atmospheric layer at the target time based on meteorological parameters at different altitudes in different atmospheric layers at the target time. The vertical integral column concentration calculation formula is used to calculate the vertical integral column concentration of the target species in each atmospheric stratum at the target time based on the volume mixing ratio and gas molar volume of the target species at different altitudes in each atmospheric stratum.

4. The method for quantifying the contribution of ozone sources within the boundary layer according to claim 1, characterized in that, The formula for calculating the contribution of the nighttime residual layer entrainment is as follows: ； The formula for calculating the contribution of surface emissions is as follows: ； The formula for calculating the contribution of chemical generation is as follows: ； in, Indicates the contribution of the nighttime residual layer entrainment, unit: molecule cm -2 ; This represents the average vertical integral column concentration of total oxidant in the residual layer during the night, in molecules cm⁻¹. -2 ; This represents the vertical integral column concentration of total oxidant in the residual layer at time t during the day, in molecules cm⁻¹. -2 ; Indicates the contribution of ground-level emissions, unit: molecule cm -2 ; This represents the sum of the vertical integral column concentrations of nitrogen oxides in the residual layer and boundary layer at time t during the day, in molecules cm⁻¹. -2 ; This represents the sum of the vertical integral column concentrations of nitrogen oxides in the residual layer and boundary layer during the night, in molecules cm⁻¹. -2 ER represents the percentage of nitrogen dioxide emissions in primary nitrogen oxide emissions; Indicates chemical contribution, unit: molecule cm -2 ; This represents the sum of the vertical integral column concentrations of total oxidant in the boundary layer and residual layer at time t during the day, in molecules cm⁻¹. -2 ; This represents the sum of the average vertical integral column concentrations of total oxidants in the inner boundary layer and residual layer at night, in molecules cm⁻¹. -2 .

5. A system for quantifying the contribution of ozone sources within the boundary layer, characterized in that, The boundary layer ozone source contribution quantification system is applied to the boundary layer ozone source contribution quantification method according to any one of claims 1-4, and the boundary layer ozone source contribution quantification system includes: The acquisition module is used to acquire vertical observation data of meteorological parameters and target species at different altitudes in different atmospheric strata within a target time period; the meteorological parameters include temperature and air pressure; the target species include ozone, nitric oxide, and nitrogen dioxide; the atmospheric strata include the boundary layer and the residual layer; the target time includes multiple times during the day and multiple times at night. The first calculation module is used to calculate the vertical integral column concentration of the target species in different atmospheric strata within the target time using the vertical integral column concentration calculation formula based on meteorological parameters at different altitudes in different atmospheric strata within the target time and vertical observation data of the target species; the vertical integral column concentration calculation formula is as follows: ； in, This represents the vertical integral column concentration of target species i in the atmospheric stratification layer c at time t, in molecules cm⁻¹. -2 ; The volume mixing ratio of target species i at altitude h in the atmospheric stratification layer c at time t is dimensionless; h1 represents the lowest altitude in the atmospheric stratification layer c, in cm; h2 represents the highest altitude in the atmospheric stratification layer c, in cm; N A Represents Avogadro's constant; This represents the molar volume of gas at altitude h within the c-layer atmospheric stratification at time t, in cm³. 3 mol -1 ; The second calculation module is used to calculate the nighttime remnant layer entrainment contribution and chemical generation contribution based on the vertical integral column concentration of target species in different atmospheric strata within the target time. Specifically, this includes: classifying target species according to preset conditions to obtain total oxidants and nitrogen oxides; the preset conditions are: ozone and nitrogen dioxide as total oxidants; nitric oxide and nitrogen dioxide as nitrogen oxides; obtaining the vertical integral column concentrations of nitrogen oxides and total oxidants in different atmospheric strata within the target time based on the vertical integral column concentrations of target species, total oxidants, and nitrogen oxides in different atmospheric strata within the target time; calculating the nighttime remnant layer entrainment contribution based on the change in the vertical integral column concentration of total oxidants in different atmospheric strata within the target time; calculating the ground emission contribution based on the change in the vertical integral column concentration of nitrogen oxides in different atmospheric strata within the target time; and calculating the chemical generation contribution based on the vertical integral column concentrations of nitrogen oxides, total oxidants, and ground emission contributions in different atmospheric strata within the target time. The third calculation module is used to determine the source of ozone variation in the boundary layer based on the contribution of the nighttime residual layer and the contribution of chemical generation.

6. The boundary layer ozone source contribution quantification system according to claim 5, characterized in that, The acquisition module specifically includes: Use at least one of airships, drones, and lidar to acquire vertical observation data of meteorological parameters and target species at different altitudes in different atmospheric strata at the target time.

7. An electronic device, characterized in that, It includes a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor, when executing the computer program, implements the method for quantifying the contribution of ozone sources within the boundary layer as claimed in any one of claims 1 to 4.

8. A computer-readable storage medium, characterized in that, The storage medium stores a computer program that, when executed, implements the method for quantifying the contribution of ozone sources within the boundary layer as described in any one of claims 1 to 4.

Images