Ozone nonlinear equivalence curve graph construction method and device, equipment and medium

Abstract

The embodiment of the invention provides an ozone nonlinear equivalence curve graph construction method and device, equipment and a medium. The ozone non-linear equivalence curve graph construction method comprises the following steps: obtaining a precursor monitoring concentration sequence obtained by monitoring in the same time period, and obtaining a plurality of groups of adjacent precursor simulation concentration sequences based on the monitoring concentration sequence; obtaining a plurality of first ozone characteristic parameters based on the monitoring concentration sequence or a group of precursor simulation concentration sequences; carrying out derivative calculation based on each first ozone characteristic parameter, and determining each order of partial derivative; second ozone characteristic parameters corresponding to the concentrations of all the target precursors are calculated based on all the orders of partial derivatives; and on the basis of the first ozone characteristic parameter, the precursor concentration at the corresponding target moment, the second ozone characteristic parameter and the corresponding target precursor concentration, constructing an ozone nonlinear equivalence curve graph. Based on the method of the scheme, the ozone nonlinear equivalence curve can be obtained only through less data calculation, so that the calculation cost can be reduced.

Description

Technical Field

[0001] The present disclosure relates to the field of environmental protection technology, and in particular to a method, device, equipment, and medium for constructing a nonlinear isovalue curve diagram of ozone. Background Art

[0002] Ozone and pollutants in the atmospheric environment are secondary products. Their concentrations have a very complex nonlinear relationship with their precursors, volatile organic compounds (VOCs) and nitrogen oxides (NOx). Blindly reducing the concentration of ozone precursors will not only fail to reduce ozone concentrations in most cases, but will cause a rebound in ozone concentrations. The industry generally refers to the analysis of the impact of precursor changes on ozone as ozone sensitivity analysis, which is the basis for formulating ozone control plans. There are many methods for ozone sensitivity analysis, and the most common method is to draw ozone nonlinear contour curves based on observation model simulation data.

[0003] The brute-force method is generally used to draw the ozone nonlinear contour curve based on the observation model simulation data. This method calculates the ozone generation potential at different concentrations by continuously changing the ratio of precursor concentrations, thereby drawing the ozone nonlinear contour curve diagram, thereby analyzing the sensitivity of ozone and precursors. However, this method requires constantly changing the precursor concentration value for simulation, which requires a large amount of calculation and causes a long calculation time. Summary of the invention

[0004] In order to solve the problem that the existing brute force method is used to construct an ozone nonlinear contour curve diagram, the calculation amount is large and the calculation is time-consuming, the embodiment of the present disclosure provides a new ozone nonlinear contour curve diagram construction method and device, equipment, and medium.

[0005] In a first aspect, an embodiment of the present disclosure provides a method for constructing an ozone nonlinear contour curve diagram, comprising: Acquire precursor monitoring concentration sequences, meteorological data sequences and light flux sequences monitored in the same period, and obtain multiple groups of adjacent precursor simulation concentration sequences based on the monitoring concentration sequences and preset concentration adjustment values; the monitoring concentration sequences and the simulation concentration sequences both include volatile organic compound concentration sequences and nitrogen oxide concentration sequences; Based on the monitoring concentration sequence or a group of the precursor simulation concentration sequences, respectively, in combination with the meteorological data sequence and the light flux sequence, a photochemical reaction simulation calculation is performed to obtain a plurality of first ozone characteristic parameters corresponding to the target time; the first ozone characteristic parameter is ozone generation potential or ozone concentration; Perform derivative calculation based on each first ozone characteristic parameter and the preset concentration adjustment value to determine partial derivatives of various orders; Based on the partial derivatives of each order, the preset concentration adjustment value, and the target moment precursor concentration corresponding to each of the first ozone characteristic parameters, the second ozone characteristic parameter of the neighboring other precursor concentrations is calculated by summing the Taylor expansion and ignoring the high-order expansion terms, and the calculation is repeated based on the neighboring other precursor concentrations and the corresponding second ozone characteristic parameter, the partial derivative, and the preset concentration adjustment parameter until the second ozone characteristic parameters corresponding to all target precursor concentrations are calculated; Based on the first ozone characteristic parameter and the precursor concentration at the corresponding target time, the second ozone characteristic parameter and the corresponding target precursor concentration, an ozone nonlinear contour curve diagram is constructed.

[0006] Optionally, the photochemical reaction simulation calculation is performed based on the monitoring concentration sequence or one of the precursor simulation concentration sequences, in combination with the meteorological data sequence and the light flux sequence, to obtain a plurality of first ozone characteristic parameters corresponding to the target time, including: Performing a photochemical reaction simulation calculation based on the combination of the monitoring concentration sequence, the meteorological data sequence and the light flux sequence to obtain the candidate ozone characteristic parameters at each moment of the monitoring concentration sequence; Selecting the largest candidate ozone characteristic parameter as the first ozone characteristic parameter corresponding to the monitoring concentration sequence, and taking the time corresponding to the largest candidate ozone characteristic parameter as the target time; Based on each of the precursor simulated concentration sequences, the simulated photochemical reaction is calculated in combination with the meteorological data sequence and the light flux sequence to obtain other first ozone characteristic parameters corresponding to the target time.

[0007] Optionally, the precursor monitoring concentration sequence includes an acetylene concentration sequence; the method further comprises: obtaining an ozone concentration sequence monitored during the same period; Before performing photochemical reaction simulation calculation based on the monitoring concentration sequence or one of the precursor simulation concentration sequences, in combination with the meteorological data sequence and the light flux sequence, the method further comprises: Calculating formaldehyde concentration based on the monitored ozone concentration at each moment of the ozone concentration sequence and the acetylene concentration at the corresponding moment of the acetylene concentration sequence, and constructing a formaldehyde concentration sequence based on the formaldehyde concentration; The photochemical reaction simulation calculation based on the monitored concentration sequence or one of the precursor simulated concentration sequences, combined with the meteorological data sequence and the light flux sequence, includes: performing photochemical reaction simulation calculation based on the monitored concentration sequence or one of the precursor simulated concentration sequence, combined with the meteorological data sequence, the light flux sequence and the formaldehyde concentration sequence.

[0008] Optionally, the precursor monitoring concentration sequence includes a nitrogen dioxide concentration sequence; Before performing photochemical reaction simulation calculation based on the monitoring concentration sequence or one of the precursor simulation concentration sequences, in combination with the meteorological data sequence and the light flux sequence, the method further comprises: Calculating nitrous acid concentration based on the monitored nitrogen dioxide concentration at each moment of the nitrogen dioxide concentration sequence, and constructing a nitrous acid concentration sequence based on the nitrous acid concentration; The photochemical reaction simulation calculation based on the monitored concentration sequence or one of the precursor simulated concentration sequences, combined with the meteorological data sequence and the light flux sequence, includes: performing photochemical reaction simulation calculation based on the monitored concentration sequence or one of the precursor simulated concentration sequence, combined with the meteorological data sequence, the light flux sequence and the nitrite concentration sequence.

[0009] Optionally, before performing the photochemical reaction simulation calculation based on the monitoring concentration sequence or one of the precursor simulation concentration sequences, in combination with the meteorological data sequence and the light flux sequence, the method further includes: Performing standard mechanism substance concentration conversion on various volatile organic compound concentration sequences in the monitoring concentration sequence and the precursor simulation concentration sequence, and determining a standard mechanism substance concentration sequence corresponding to the monitoring concentration sequence and the precursor simulation concentration sequence; The photochemical reaction simulation calculation is performed based on the monitored concentration sequence or one of the precursor simulation concentration sequences, in combination with the meteorological data sequence and the light flux sequence, including: the photochemical reaction simulation calculation is performed based on the standard mechanism substance concentration sequence and volatile organic compound concentration sequence corresponding to the monitored concentration sequence or one of the precursor simulation concentration sequence, in combination with the meteorological data sequence and the light flux sequence.

[0010] Optionally, the method for obtaining the light flux sequence includes: obtaining longitude and latitude information of the monitoring location, and performing illumination orthographic projection calculation based on the longitude and latitude information to determine an orthographic projection conversion coefficient; The luminous flux sequence is obtained based on the incident intensity of sunlight and the orthographic projection conversion coefficient.

[0011] Optionally, before performing photochemical reaction simulation calculation based on the monitoring concentration sequence or one of the precursor simulation concentration sequences, in combination with the meteorological data sequence and the light flux sequence, the method further comprises: determining a mixing layer height parameter; The performing photochemical reaction simulation calculation based on the monitoring concentration sequence or one of the precursor simulation concentration sequences, in combination with the meteorological data sequence and the light flux sequence, comprises: The mixing layer height parameter is used as a boundary layer constraint, and based on the monitored concentration sequence or a precursor simulated concentration sequence, a photochemical reaction simulation calculation is performed in combination with the meteorological data sequence and the light flux sequence.

[0012] In a second aspect, an embodiment of the present disclosure provides a device for constructing a nonlinear contour graph of ozone, comprising: A basic data determination unit, used to obtain a precursor monitoring concentration sequence, a meteorological data sequence and a light flux sequence monitored in the same period, and obtain multiple groups of adjacent precursor simulation concentration sequences based on the monitoring concentration sequence and a preset concentration adjustment value; the monitoring concentration sequence and the simulation concentration sequence both include a volatile organic compound concentration sequence and a nitrogen oxide concentration sequence; a photochemical reaction model calculation unit, configured to perform photochemical reaction simulation calculations based on the monitoring concentration sequence or a group of precursor simulation concentration sequences, in combination with the meteorological data sequence and the light flux sequence, to obtain a plurality of first ozone characteristic parameters corresponding to a target time; the first ozone characteristic parameter is ozone generation potential or ozone concentration; A derivative calculation unit, used to perform derivative calculation based on each first ozone characteristic parameter and the preset concentration adjustment value to determine partial derivatives of various orders; a point value data calculation unit, for calculating the second ozone characteristic parameter of the neighboring other precursor concentrations in a manner of Taylor expansion summation and ignoring high-order expansion terms based on the partial derivatives of each order, the preset concentration adjustment value, and the target moment precursor concentration corresponding to each of the first ozone characteristic parameters, and repeating the calculation based on the neighboring other precursor concentrations and the corresponding second ozone characteristic parameter, the partial derivative, and the preset concentration adjustment parameter until the second ozone characteristic parameters corresponding to all target precursor concentrations are calculated; The curve graph construction unit is used to construct an ozone nonlinear contour curve graph based on the first ozone characteristic parameter and the precursor concentration at the corresponding target time, the second ozone characteristic parameter and the corresponding target precursor concentration.

[0013] In a third aspect, an embodiment of the present disclosure provides a computing device, including a processor and a memory, wherein the memory is used to store a computer program; when the computer program is loaded by the processor, the processor executes the method for constructing an ozone nonlinear contour curve diagram as described above.

[0014] In a fourth aspect, an embodiment of the present disclosure provides a computer-readable storage medium, wherein the storage medium stores a computer program. When the computer program is executed by a processor, the processor implements the method for constructing a nonlinear contour curve diagram of ozone as described above.

[0015] The disclosed embodiment scheme considers that only a small number of precursor sequences are needed to perform photochemical simulation to obtain the first ozone characteristic parameters of multiple concentration points, and after each order derivative is calculated using the aforementioned first ozone characteristic parameters, the second ozone characteristic parameters corresponding to other precursor concentrations are obtained using the aforementioned each order derivative according to the Taylor expansion method, and then the ozone nonlinear isovalue curve can be constructed based on the first ozone characteristic parameter and the second ozone characteristic parameter. Based on the method of this scheme, the ozone nonlinear isovalue curve can be obtained by calculating only a small amount of data, which can also reduce the calculation cost. BRIEF DESCRIPTION OF THE DRAWINGS

[0016] The accompanying drawings, which are incorporated in and constitute a part of this specification, illustrate embodiments consistent with the present disclosure and, together with the description, serve to explain the principles of the present disclosure.

[0017] In order to more clearly illustrate the technical solutions in the embodiments of the present disclosure or the prior art, the following briefly introduces the drawings required for use in the embodiments or the prior art description. Obviously, for ordinary technicians in this field, other drawings can be obtained based on these drawings without creative labor, including Figure 1 is a flow chart of a method for constructing a nonlinear isovalue curve of ozone provided in an embodiment of the present disclosure; Figure 2 is a nonlinear contour graph of ozone concentration constructed using the embodiment of the present disclosure; Figure 3 It is a nonlinear contour plot of ozone concentration constructed using the brute force method; Figure 4 is a nonlinear contour graph of ozone generation potential constructed using the embodiment scheme of the present disclosure; Figure 5 It is a nonlinear contour plot of ozone generation potential constructed using the brute force method; Figure 6 It is a schematic diagram of the architectural structure based on which the present disclosure is based in the specific implementation of the aforementioned method; Figure 7It is a schematic diagram of a device for constructing a nonlinear isovalue curve of ozone provided in an embodiment of the present disclosure; Figure 8 It is a schematic diagram of the structure of a computing device provided in an embodiment of the present disclosure. DETAILED DESCRIPTION

[0018] Embodiments of the present disclosure will be described in more detail below with reference to the accompanying drawings. Although certain embodiments of the present disclosure are shown in the accompanying drawings, it should be understood that the present disclosure can be implemented in various forms and should not be construed as being limited to the embodiments described herein, which are instead provided to provide a more thorough and complete understanding of the present disclosure. It should be understood that the drawings and embodiments of the present disclosure are only for exemplary purposes and are not intended to limit the scope of protection of the present disclosure.

[0019] The term "including" and its variations used in this document are open inclusions, i.e., "including but not limited to". The term "based on" means "based at least in part on". The term "one embodiment" means "at least one embodiment"; the term "another embodiment" means "at least one other embodiment"; the term "some embodiments" means "at least some embodiments". Relevant definitions of other terms will be given in the following description. In this document, relational terms such as "first" and "second" are merely used to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply these embodiments.

[0020] In order to solve the problem that the existing brute force method is very time-consuming and requires a lot of calculation to calculate and draw ozone nonlinear contour curves, the present disclosure provides a new method for constructing ozone nonlinear contour curves. The method for constructing ozone nonlinear contour curves provided by the present disclosure is executed by a computing device.

[0021] Figure 1 : is a flow chart of the method for constructing the ozone nonlinear contour curve provided by the embodiment of the present disclosure. Figure 1 As shown, the method for constructing the ozone nonlinear contour curve provided by the embodiment of the present disclosure includes S110-S150.

[0022] S110: Acquire precursor monitoring concentration sequences, meteorological data sequences and light flux sequences monitored in the same period, and obtain multiple groups of adjacent precursor simulation concentration sequences based on the monitoring concentration sequences and preset concentration adjustment values.

[0023] Considering that ozone in the air is caused by photochemical reactions, light is a prerequisite for achieving photochemical reactions, the same period of time in the embodiments of the present disclosure refers to a period of sunlight. In specific implementations, the same period of time can be determined based on the sunrise and sunset times of a certain area.

[0024] The precursor monitoring concentration sequence is a monitoring concentration sequence obtained by arranging the concentrations of a certain ozone-generating precursor monitored by the monitoring station within the same time period. The precursor concentration monitoring sequence is a sequence in which the precursor concentrations are sorted in chronological order.

[0025] Because the precursors for ozone formation include volatile organic compounds and nitrogen oxides, the precursor concentration sequence in the embodiment of the present disclosure includes an organic concentration sequence and a nitrogen oxide concentration sequence. In a specific implementation, the type of volatile organic compounds corresponding to the aforementioned organic concentration sequence is determined, which can be various possible organic compounds that cause regional environmental pollution, specifically determined based on environmental pollution monitoring experience. Similarly, the type of volatile nitrogen oxides corresponding to the nitrogen oxide concentration sequence is also determined, which can also be determined based on environmental monitoring experience. In a specific implementation, volatile organic compounds may include carbon monoxide, various volatile alkane compounds, benzene compounds, and volatile oxides of the aforementioned various compounds; nitrogen oxides may be nitric oxide, nitrogen dioxide, etc.

[0026] Table 1 is a table of volatile organic compound monitoring concentration sequence data provided by some embodiments. It can be seen from Tables 1-2 that the aforementioned same time period is set to 7:00-12:00. In order to avoid leakage of actual environmental monitoring data, the volatile organic compounds in Table 1 are represented by codes, where VOC1, VOC4, VOC14, ..., VOC49 represent a volatile organic compound respectively. Table 2 is a table of nitrogen oxide monitoring concentration sequence data provided by some embodiments, in which the nitrogen oxides included are nitric oxide and nitrogen dioxide, etc.

[0027] Table 1 Volatile organic compound monitoring concentration sequence table

[0028]

[0029] Table 2 Nitrogen oxide concentration sequence table

[0030] The meteorological data sequence is a sequence composed of the meteorological data of a certain period of time monitored in chronological order. The meteorological data sequence may include parameters such as temperature, air humidity, relative air humidity, wind speed and wind direction. Considering the influence of ambient temperature and air humidity on the photochemical reaction rate, the meteorological data at least includes temperature data, air humidity or relative air humidity data. Table 3 is a table of temperature data sequences provided by some embodiments. It reflects the temperature changes in a specific area within a certain period of time.

[0031] Table 3 Temperature sequence table for a certain period of time

[0032] The luminous flux sequence is a sequence of luminous fluxes in a certain period of time determined by monitoring or data conversion in a time sequence. The aforementioned luminous flux is the radiation flux irradiated to a certain volume of air from all directions. The luminous flux is related to the incident intensity of the sun. The luminous flux corresponding to the luminous flux sequence in the embodiment of the present disclosure is a sequence corresponding to light of a wavelength that can affect the progress of the photochemical reaction.

[0033] In some more embodiments of the present disclosure, the light flux sequence can be determined by data conversion, specifically including the following S111-S112.

[0034] S111: Obtain the longitude and latitude information of the monitoring location, and perform illumination orthographic projection calculation based on the longitude and latitude information to determine the orthographic projection conversion coefficient.

[0035] S112: Obtain a luminous flux sequence based on the incident intensity of sunlight and the orthographic projection conversion coefficient.

[0036] According to geographical knowledge, the latitude and longitude information and time point information of the monitoring location have an impact on the light intensity received by the unit area. When the incident intensity of sunlight is determined, it is necessary to determine the incident angle of sunlight in order to determine the orthographic projection light intensity at the corresponding position based on the incident angle. According to the aforementioned logic, the embodiment of the present disclosure performs an orthographic projection transformation of the illumination based on the longitude and latitude information of the monitoring point to determine the orthographic projection conversion coefficient. The orthographic projection conversion coefficient is a coefficient that reflects the conversion of the corresponding light intensity to the orthographic projection light intensity when the unit light intensity is incident at a specific incident angle.

[0037] After determining the orthographic projection conversion coefficient, the orthographic projection conversion coefficient and the incident intensity of sunlight can be multiplied to obtain the illumination intensity at the corresponding time point. Because the energy ratio of various wavelengths in the sunlight spectrum is basically determined, the luminous flux of different wavelengths at each time can be determined according to the above method. By sorting the above luminous flux in chronological order, a luminous flux sequence can be obtained.

[0038] In the disclosed embodiment, after obtaining the precursor monitoring concentration sequence, numerical addition and subtraction processing will be performed based on each sequence in the aforementioned set of precursor monitoring concentrations and the preset concentration adjustment value to obtain a precursor simulation concentration sequence similar to the precursor monitoring concentration sequence. As the name suggests, the precursor simulation concentration sequence is not a real sequence, but an assumed sequence. The preset concentration adjustment value is a concentration adjustment value determined based on actual experience and the required numerical accuracy. For example, under normal circumstances, the concentration variation range of a certain pollutant is 0-300. For this pollutant, 16 points need to be arranged within its concentration variation range. The corresponding preset concentration adjustment value can be 20 (this is just an example here, so the dimension problem is not considered).

[0039] It should be noted that in the embodiments of the present disclosure, when a set of precursor simulation concentration sequences are obtained by performing numerical addition and subtraction processing according to a set of precursor monitoring concentration sequences and preset concentration adjustment values, the preset concentration adjustment values ​​are simultaneously increased or decreased for the organic concentration sequences corresponding to different organic substances therein to obtain the organic concentration adjustment sequences that change in the same regularity, and the preset concentration adjustment values ​​are simultaneously increased or decreased for different nitrogen oxide concentration sequences to obtain the nitrogen oxide concentration sequences that change in the same regularity.

[0040] In some embodiments of the present disclosure, the number of multiple groups of adjacent precursor simulation concentration sequences determined based on the precursor monitoring concentration sequence is 8 groups. The aforementioned 8 groups of adjacent precursor simulation concentration sequences can be centered on the precursor monitoring concentration sequence, or can be located on one side of the precursor monitoring concentration sequence. Assuming that each concentration sequence is a point, the aforementioned 9 concentration sequences form a square with a side length of 2 times the preset concentration adjustment value (at this time, one side of the square represents the change in the concentration of volatile organic compounds, and the other side represents the change in the concentration of volatile nitrogen oxides). In the embodiment of the present disclosure, the goal of determining 8 groups of adjacent precursor concentration sequences is to obtain first-order partial derivatives and second-order partial derivatives for subsequent summation calculation based on Taylor expansion. In other embodiments, if only first-order partial derivatives are required when performing summation calculation based on Taylor expansion, the number of adjacent precursor simulation concentration sequences only needs 3 groups. If third-order partial derivatives or even higher-order partial derivatives are required when performing summation calculation based on Taylor expansion, more groups of adjacent precursor simulation concentration sequences are required.

[0041] In specific implementations, while monitoring the precursor concentration sequence, the concentration of pollutants such as sulfur dioxide may also be monitored. Because they have a similar form to nitrogen oxides, they can also be treated as nitrogen oxides in actual applications and processed in the same way as nitrogen oxides to obtain the corresponding simulated concentration sequence.

[0042] S120: performing photochemical reaction simulation calculations based on the monitoring concentration sequence or a set of precursor simulation concentration sequences, in combination with the meteorological data sequence and the light flux sequence, to obtain a plurality of first ozone characteristic parameters corresponding to the target time.

[0043] After obtaining the aforementioned monitoring concentration sequence or a precursor simulation concentration sequence, and after obtaining the meteorological data sequence and the light flux sequence, photochemical simulation calculations can then be performed according to the principle of photochemical reaction.

[0044] In the disclosed embodiment, the aforementioned precursor monitoring concentration sequence and each group of precursor simulation concentration sequences are respectively combined with a meteorological data sequence and a light flux sequence to obtain a plurality of first ozone characteristic parameters corresponding to the target time.

[0045] Taking the combination of precursor monitoring concentration sequence, meteorological data sequence and light flux sequence for photochemical reaction simulation calculation as an example, various chemical reaction simulations are performed based on the current precursor monitoring concentration, temperature and humidity data and light flux as conditions to determine the ozone concentration change value at the first moment. A new ozone concentration cumulative value is determined based on the historical ozone concentration cumulative value and the ozone concentration change value, and the new ozone concentration cumulative value and the precursor monitoring concentration, temperature and humidity data and light flux at the next moment are used as conditions to calculate the ozone concentration cumulative value at the next moment until the first ozone characteristic parameter at the target moment is obtained. Similarly, by combining each group of adjacent precursor simulation concentration sequences, meteorological data sequences and light flux sequences for photochemical simulation calculations, the corresponding first ozone characteristic parameters can be obtained.

[0046] In some embodiments of the present disclosure, in order to obtain the first ozone characteristic parameter based on a pre-constructed photochemical reaction model without having to reconstruct the photochemical reaction model due to changes in the types of monitored precursors (especially volatile organic compounds), in some specific embodiments, before executing the aforementioned S120, the following S01 step may be executed.

[0047] S01: convert various volatile organic compound concentration sequences in the monitoring concentration sequence and the precursor simulation concentration sequence into standard mechanism substance concentration sequences, and determine the standard mechanism substance concentration sequences corresponding to the monitoring concentration sequence and the precursor simulation concentration sequence.

[0048] Table 4 is a table of conversion coefficients of various volatile organic compounds and standard mechanism substances in carbon bond 05, in which each row represents a volatile organic compound, each column represents a standard mechanism substance, and the corresponding rows and columns have values ​​indicating that a unit amount of volatile organic compound can be converted into a corresponding amount of standard mechanism substances. It should be noted that a volatile organic compound may correspond to multiple standard mechanism substances, that is, it can be converted into a combination of multiple standard mechanism substances, which can characterize its actual role in the photochemical reaction.

[0049] Table 4 Conversion coefficients of various volatile organic compounds and standard mechanism substances in carbon bond 05

[0050]

[0051] In a specific implementation, the computing device can convert based on the conversion coefficients in Table 4 to obtain the corresponding standard mechanism substance concentration sequence. Specifically, the computing device converts the volatile organic compounds in the concentration sequence according to the standard mechanism substance conversion coefficient table to obtain the concentration data corresponding to various standard mechanism substances. Then, the concentration data of the same standard mechanism substance at the same time are added together to obtain the final standard mechanism substance concentration data. The above conversion method is expressed by the formula: ,in For the Standard mechanism The conversion concentration at the moment, For the Volatile organic compounds in The conversion concentration at the moment, for Volatile organic compounds and The conversion coefficient of the standard mechanism substance, It can be found in Table 3.

[0052] After obtaining the final standard mechanism substance concentration data, they are arranged in chronological order to obtain the standard mechanism substance concentration sequence.

[0053] On the premise of executing the aforementioned S01, the aforementioned S120 specifically includes the following S121: performing photochemical reaction simulation calculation based on the standard mechanism substance concentration sequence and the volatile organic compound concentration sequence corresponding to the monitoring concentration sequence or a precursor simulation concentration sequence, combined with the meteorological data sequence and the light flux sequence. Of course, in other implementations, the aforementioned S01 may not be executed, but the photochemical reaction simulation may be performed directly based on the chemical reaction principles of various volatile organic compounds.

[0054] In a specific implementation, the computing device can perform hard calculations based on the concentrations of various standard mechanisms and the concentrations of nitrogen oxides according to a predetermined reaction model to determine the reaction rates of each step in the photochemical reaction simulation process, and obtain the first ozone characteristic parameter based on the aforementioned reaction rates and other necessary data. The basic chemical reaction rate is , is the reaction rate constant, is the concentration of reactant species A, is the concentration of the reaction species B. The reaction species A and reaction species B mentioned above may be standard mechanism substances or nitrogen oxides according to the specific reaction mechanism. a and b are the orders of substance A and substance B in the reaction, respectively. In most cases, a and b are 1; , is the gas constant, is the activation energy, is the reaction temperature, is the pre-exponential factor.

[0055] In the embodiment of the present disclosure, the first ozone characteristic parameter is ozone generation potential or ozone concentration. In specific implementation, ozone generation potential and ozone concentration can be calculated respectively, so as to be used to subsequently construct two different types of nonlinear contour curves respectively.

[0056] The following is an analysis of how to calculate the potential before ozone generation. In specific implementation, because the generation of ozone is directly related to nitric oxide, it can be used The ozone generation potential is obtained, where is the ozone generation potential, is the net generation of ozone, yes net consumption.

[0057] In specific implementation, , .

[0058] Table 5 is the chemical reaction parameter table of carbon bond 05. The parameters in the above formula can be found through Table 4.

[0059] Table 5 Carbon bond 05 chemical reaction parameters

[0060]

[0061]

[0062]

[0063]

[0064]

[0065] In the above Table 5, k@289K represents the reaction rate constant at 25°C, R1, R8, R9, R14, R15, R25, R36, R51, R52, R53, R62, R65, R72, R75, R76, R87, R91, R97, R102, R106, R138, R143, and R152 represent the photolysis reaction, and the photolysis reaction rate constant is ,in is the luminous flux, For wavelength The electromagnetic wave absorption cross-sectional area, is the quantum yield of the photolysis reaction. According to the formula, the absorption cross section and the quantum yield of the photolysis reaction are related to the wavelength. related.

[0066] The following is an analysis of how to calculate ozone concentration. In specific implementation, ozone concentration can be obtained by integrating the ozone generation rate and time, specifically: ,in is the ozone generation rate, is the initial concentration of ozone, usually set to 0. Ozone generation rate ,in is the net ozone generation rate, is the net consumption rate of ozone.

[0067]

[0068] .

[0069] In specific implementation, due to the monitoring mechanism, the precursor monitoring concentration sequence lacks monitoring of certain precursors, and the concentration sequence of such precursors in the photochemical reaction cannot be determined. For example, formaldehyde is easily soluble in water, so its actual content is not easy to be measured by the monitoring equipment. For example, nitrous acid is also difficult to measure because it is soluble in water.

[0070] To solve the problem mentioned in the previous paragraph, in some embodiments of the present disclosure, before executing the aforementioned S120, the following S02 or S03 may also be executed.

[0071] S02: Calculate the formaldehyde concentration based on the monitored ozone concentration at each moment of the ozone concentration sequence and the acetylene concentration at the corresponding moment of the acetylene concentration sequence, and construct a formaldehyde concentration sequence based on the formaldehyde concentration.

[0072] To ensure the implementation of S02, the previous precursor monitoring concentration sequence should include the acetylene concentration sequence, and in addition, the ozone concentration sequence monitored during the same period should also be obtained.

[0073] After obtaining the acetylene concentration sequence and the ozone concentration sequence, the acetylene concentration and the ozone concentration corresponding to the same time can be calculated according to The corresponding formaldehyde concentration is calculated, where is the formaldehyde concentration, is the acetylene concentration, is the ozone concentration. Then the formaldehyde concentration is sorted by time to obtain the formaldehyde concentration sequence.

[0074] On the premise of executing the aforementioned S02, the aforementioned S120 is specifically the following S122.

[0075] S122: Based on the monitored concentration sequence or a precursor simulated concentration sequence, a photochemical reaction simulation calculation is performed in combination with a meteorological data sequence, a light flux sequence and a formaldehyde concentration sequence.

[0076] The specific implementation process of S122 is as mentioned above and will not be repeated here. It should be noted that the formaldehyde concentration sequence is a sequence of corresponding detection substance concentrations obtained for the precursor monitoring concentration, rather than a sequence directly for the precursor simulation concentration sequence. Considering that the formaldehyde concentration is small and has little effect on the photochemical reaction, it is directly combined with the precursor simulation concentration sequence to perform photochemical reaction simulation calculations to obtain the corresponding first ozone characteristic parameter. In other embodiments, the formaldehyde concentration sequence can also be numerically adjusted according to the preset concentration adjustment value to obtain the formaldehyde concentration sequence corresponding to each precursor simulation concentration sequence, and S121 is executed accordingly.

[0077] S03: Calculating nitrite concentration based on the monitored nitrogen dioxide concentration at each moment of the nitrogen dioxide concentration sequence, and constructing a nitrite concentration sequence based on the nitrite concentration.

[0078] In a specific implementation, the computing device is Calculate the nitrite concentration ,in is the concentration of nitrogen dioxide, and MR is the relative molecular mass of nitrogen dioxide. After obtaining the nitrite concentration at each moment, the nitrite concentrations are sorted in chronological order to obtain a nitrite concentration sequence.

[0079] On the premise of executing the aforementioned S03, the aforementioned S120 is specifically implemented as follows in S123.

[0080] S123: Based on the monitored concentration sequence or a precursor simulated concentration sequence, a meteorological data sequence, a light flux sequence and a nitrite concentration sequence are combined to perform a photochemical reaction simulation calculation.

[0081] The specific implementation process of S123 is as described above and will not be repeated here. It should be noted that the aforementioned nitrite concentration sequence is a sequence of corresponding detection substance concentrations obtained for the precursor monitoring concentration, rather than a sequence directly for the precursor simulation concentration sequence. Considering that the nitrite concentration is small and has little effect on the photochemical reaction, it is directly combined with the precursor simulation concentration sequence to perform photochemical reaction simulation calculations to obtain the corresponding first ozone characteristic parameter. In other embodiments, the nitrite concentration sequence can also be numerically adjusted according to the preset concentration adjustment value to obtain the nitrite concentration sequence corresponding to each precursor simulation concentration sequence, and S122 is executed accordingly.

[0082] It should also be noted that the aforementioned S02 and S03 do not conflict, and can be executed simultaneously in specific implementations, and the corresponding S122 and S123 can be combined and executed to calculate the corresponding first ozone characteristic parameter. In addition, in some embodiments, after executing the aforementioned S02 and S03, the aforementioned S01 is executed for standard mechanism conversion, and then S122 or S123 is executed in the same manner as the aforementioned S121.

[0083] After completing the execution of S120 , the following S130 may be subsequently executed.

[0084] S130: Perform derivative calculation based on each first ozone characteristic parameter and a preset concentration adjustment value to determine partial derivatives of various orders.

[0085] The partial derivatives of each bound in the embodiment of the present disclosure are partial derivatives of the first ozone characteristic parameter with respect to the concentration of volatile organic compounds and / or volatile nitrogen oxides.

[0086] As analyzed above, when the first ozone characteristic parameter is 9 data, the first-order partial derivative and the second-order partial derivative can be calculated. In a specific implementation, the difference of the first ozone characteristic parameter corresponding to the adjacent volatile organic compound concentration and / or volatile nitrogen oxide concentration can be calculated, and the first-order partial derivative and the second-order partial derivative can be calculated using the above difference.

[0087] According to the above analysis, taking the first ozone characteristic parameter in the middle of the assumed square as the benchmark, the differences calculated according to the first ozone characteristic parameters on both sides are and , then we can determine the first-order partial derivative using , the second-order partial derivative is According to the above method, the first-order partial derivative and the second-order partial derivative for the change of volatile organic compound concentration and the change of nitrogen oxide concentration can be calculated.

[0088] As analyzed above, if the first ozone characteristic parameter is other quantities, it may be necessary to calculate the partial derivative of the corresponding order.

[0089] S140: Based on the partial derivatives of each order, the preset concentration adjustment value, and the target moment precursor concentration corresponding to each first ozone characteristic parameter, the second ozone characteristic parameter of other neighboring precursor concentrations is calculated by summing the Taylor expansion and ignoring the high-order expansion terms, and the calculation is repeated based on the neighboring other precursor concentrations and the corresponding second ozone characteristic parameters, partial derivatives, and preset concentration adjustment parameters until the second ozone characteristic parameters corresponding to all target precursor concentrations are calculated.

[0090] According to the principle of calculus, based on Taylor expansion formula, near The second ozone characteristic parameter at the position can be obtained by the Taylor expansion formula, and the corresponding terms of the high-order derivatives in the Taylor expansion formula are relative to the previous constant terms and low-order derivative terms, so the high-order derivative terms can be ignored and will not have a significant impact on the results. With the truth The difference is smaller.

[0091] In addition, through analysis based on a large number of ozone nonlinear contour curves, it was found that, except for the areas with low precursor concentrations, the shapes of the ozone nonlinear contour curves in other areas are basically similar, and the corresponding values ​​of the low-order derivatives are not much different.

[0092] Based on the above analysis, the second ozone characteristic parameter of the adjacent area corresponding to the first ozone characteristic parameter area can be calculated by Taylor expansion and ignoring high-order infinitesimals. For example, if the second-order derivative is used, the second ozone characteristic parameter of the adjacent position can be calculated using Calculated.

[0093] Similarly, according to the aforementioned principle, after obtaining the second ozone characteristic parameter of the adjacent position, the second ozone characteristic parameter of the position adjacent to the adjacent position can be obtained, and the second ozone characteristic parameters corresponding to all target precursor concentrations can be obtained by executing the cycle.

[0094] S150: constructing an ozone nonlinear contour curve diagram based on the first ozone characteristic parameter and the precursor concentration at the corresponding target time, the second ozone characteristic parameter and the corresponding target precursor concentration.

[0095] After obtaining the aforementioned first ozone characteristic parameters and the second ozone characteristic parameters, as well as the precursor concentrations at the target time corresponding to the aforementioned ozone characteristic parameters, the aforementioned precursor concentrations can then be standardized (specifically, the various precursor concentrations are converted into a certain type of precursor or a specific functional group), and the ozone nonlinear isovalue curve can be constructed. In the disclosed embodiment, the two coordinate axes of the ozone nonlinear isovalue curve are the normalized values ​​of the volatile organic compound concentration and the normalized values ​​of the volatile nitrogen oxide concentration.

[0096] As analyzed above, by analyzing a large number of ozone nonlinear contour curves obtained by the brute force method, it is determined that the curve characteristics at each position of the ozone nonlinear contour curve have the same gradient change law. At the same time, the goal of drawing the ozone nonlinear contour curve is not to obtain absolutely accurate results, but more to obtain a curve trend diagram that can guide practical applications. Based on the above premise, the embodiment scheme of the present disclosure considers that only a small number of precursor sequences are needed to perform photochemical simulation to obtain the first ozone characteristic parameters of multiple concentration points, and then the first ozone characteristic parameters are used to calculate the derivatives of each order. The derivatives of each order can be used to obtain the second ozone characteristic parameters corresponding to the concentration of other precursors according to the Taylor expansion method, and then the ozone nonlinear contour curve can be constructed based on the first ozone characteristic parameter and the second ozone characteristic parameter. Based on the method of this scheme, the ozone nonlinear contour curve can be obtained by calculating only a small number of data, which can also reduce the calculation cost.

[0097] As mentioned above, it is necessary to obtain a plurality of first ozone characteristic parameters corresponding to the target time. In a specific implementation, the target time is a time selected and determined manually, or a time determined according to a predetermined rule. In some embodiments, the computing device may determine the target time using the following S124-S125.

[0098] S124: Perform photochemical reaction simulation calculation based on the combination of the monitoring concentration sequence, the meteorological data sequence and the light flux sequence to obtain candidate ozone characteristic parameters at each moment of the monitoring concentration sequence.

[0099] The execution process of S124 is as analyzed above, except that the computing device will perform photochemical reaction simulation calculations on all the time data within the same period mentioned above to determine the candidate ozone characteristic parameters at each time.

[0100] S125: Selecting the largest candidate ozone characteristic parameter as the first ozone characteristic parameter corresponding to the monitoring concentration sequence, and taking the time corresponding to the largest candidate ozone characteristic parameter as the target time.

[0101] After obtaining the candidate ozone characteristic parameters at each time, the candidate ozone characteristic parameters are then compared to determine the maximum value, and the maximum value is used as the first ozone characteristic parameter. Accordingly, the time corresponding to the first ozone characteristic parameter is used as the target time.

[0102] By analyzing the normal line of the ozone nonlinear characteristic curve determined by the brute force method, in the area with higher precursor concentration (corresponding to the area with larger ozone characteristic parameters), the curvature of the ozone nonlinear characteristic curve is larger and the gradient authentication is more obvious. Accordingly, by adopting the aforementioned S124-S125, the first ozone characteristic parameter can be positioned in the area corresponding to the aforementioned more obvious gradient characteristics, so that the values ​​of the derivatives of each order calculated subsequently are more typical and representative.

[0103] On the premise of executing the aforementioned S124-S125, when processing the simulated concentration sequence of each precursor, the following S126 can also be executed.

[0104] S126: Perform photochemical reaction simulation calculation based on the simulated concentration sequence of each precursor, combined with the meteorological data sequence and the light flux sequence, to obtain other first ozone characteristic parameters corresponding to the target time.

[0105] The execution method of S126 is the same as the steps explained above, except that when performing photochemical reaction simulation, the simulation calculation is stopped when the target time is reached, and only the first ozone characteristic parameter up to the target time is obtained.

[0106] It should be noted here that the implementation process of S124-S126 does not conflict with S01, S02, S03, S121, S122 and S123 in the previous text. In the specific implementation, the execution process of each detailed step can be logically combined to form a new solution that integrates each execution step. That is to say, the formaldehyde concentration sequence and nitrite concentration sequence can be determined first, and then the concentration sequences of each precursor and the aforementioned formaldehyde concentration sequence can be converted into a standard mechanism concentration sequence. The standard mechanism concentration sequence is used to perform photochemical reaction simulation calculations to obtain the corresponding first ozone characteristic parameters.

[0107] When simulating photochemical reactions, in addition to the aforementioned precursor concentration sequence, meteorological data sequence, and light flux sequence, in actual situations, meteorological conditions and other reasons may also cause changes in the height of the pollutant gas boundary layer. Changes in the pollutant gas boundary layer will affect the pollutant gas concentration and the projection and scattering of sunlight, thereby causing changes in the actual photochemical reaction process.

[0108] In some embodiments of the present disclosure, before executing the aforementioned S120, the computing device also determines a mixing layer height coefficient. The mixing layer height coefficient is used to characterize the height adjustment characteristics of the precursor mixing layer near the surface. On the premise of determining the height adjustment characteristics of the mixing layer, the computing device performs the aforementioned S120 as follows: using the mixing layer height parameter as a boundary layer constraint, based on a monitored concentration sequence or a precursor simulated concentration sequence, combined with a meteorological data sequence and a light flux sequence, a photochemical reaction simulation calculation is performed to obtain a first ozone characteristic parameter.

[0109] In the specific implementation, in order to verify the rationality of the aforementioned method, in a specific application, the computing device also uses the brute force method to draw the ozone nonlinear contour curve. Specifically, it is based on the maximum allowable limit of volatile organic compound concentration and the maximum allowable limit of nitrogen oxide concentration determined by the existing environmental control and the set percentage change, as well as the various precursor types obtained in the aforementioned scheme. Various measurement bases are adjusted according to the set percentage to obtain 441 data points, and the 441 data points are photochemically simulated according to the aforementioned environmental conditions to obtain the first ozone characteristic parameters of each point, and the ozone nonlinear contour curve is directly drawn according to all the first ozone characteristic parameters. The ozone nonlinear contour curve constructed by the aforementioned method is used to verify the rationality of the change characteristics of each position of the ozone nonlinear contour curve constructed in the previous S110-S150.

[0110] Figure 2 It is a nonlinear contour curve diagram of ozone concentration constructed using the embodiment scheme of the present disclosure. Figure 3 It is a nonlinear contour plot of ozone concentration constructed using the brute force method. Figure 4 It is a nonlinear contour curve diagram of ozone generation potential constructed using the embodiment scheme of the present disclosure. Figure 5 It is a nonlinear contour curve of ozone generation potential constructed by the brute force method. Figure 2 and Figure 3 ,as well as Figure 4 and Figure 5 It can be seen that the ozone nonlinear isovalue curve constructed by this scheme has the same trend as the corresponding area in the curve diagram drawn by the strong force method when the precursor concentration is high, and the same trend as the corresponding area in the curve diagram drawn by the strong force method when the precursor concentration is low (of course, this is somewhat related to the choice of precursor concentration used to calculate the first ozone characteristic curve. It should be noted here that the new precursor concentration point is determined by adjusting the precursor concentration percentage by 5% in the calculation, Figure 2 and Figure 4The coordinate axis data in the figure are not the same, but the precursor concentrations reflected by the binary values ​​are corresponding). Considering that the ozone generation rate is very low when the precursor concentration is low, it meets the conventional environmental protection requirements. This is not the focus of actual environmental monitoring and pollution source control. The above differences can be ignored in practical applications.

[0111] Figure 6 It is a schematic diagram of the architectural structure based on which the present disclosure is based in the specific implementation of the aforementioned method. Figure 6 As shown, the present disclosure first needs to obtain environmental observation data, meteorological observation data (including mixing layer height coefficient, that is, boundary layer data), photolysis parameters and geographic information data, and then after data sorting, a model calculation method is used to calculate the first ozone characteristic parameter and the second ozone characteristic parameter, and then the ozone nonlinear contour curve diagram is constructed according to the method of this embodiment and the brute force method, and a sensitivity analysis is performed to determine the feasibility of implementing the method of this embodiment.

[0112] In addition to providing the aforementioned method for constructing an ozone nonlinear contour curve diagram, the embodiment of the present disclosure also provides an ozone nonlinear contour curve construction device. Figure 7 Schematic diagram of a device for constructing an ozone nonlinear contour curve provided by an embodiment of the present disclosure. Figure 7 As shown, the ozone nonlinear isovalue curve construction device 700 includes a basic data determination unit 701, a photochemical reaction model calculation unit 702, a derivative calculation unit 703, a point value data calculation unit 704 and a curve graph construction unit 705.

[0113] The basic data determination unit 701 is used to obtain the precursor monitoring concentration sequence, meteorological data sequence and light flux sequence monitored in the same time period, and obtain multiple groups of adjacent precursor simulation concentration sequences based on the monitoring concentration sequence and the preset concentration adjustment value; the monitoring concentration sequence and the simulation concentration sequence both include volatile organic compound concentration sequence and nitrogen oxide concentration sequence.

[0114] The photochemical reaction model calculation unit 702 is used to perform photochemical reaction simulation calculations based on a monitored concentration sequence or a set of precursor simulated concentration sequences, in combination with a meteorological data sequence and a light flux sequence, to obtain a plurality of first ozone characteristic parameters corresponding to the target time; the first ozone characteristic parameter is the ozone generation potential or the ozone concentration.

[0115] The derivative calculation unit 703 is used to perform derivative calculation based on each first ozone characteristic parameter and a preset concentration adjustment value to determine partial derivatives of various orders.

[0116] The point value data calculation unit 704 is used to calculate the second ozone characteristic parameter of other neighboring precursor concentrations based on the partial derivatives of each order, the preset concentration adjustment value, and the target moment precursor concentration corresponding to each first ozone characteristic parameter, by summing up the Taylor expansion and ignoring the higher-order expansion terms, and repeat the calculation based on the neighboring other precursor concentrations and the corresponding second ozone characteristic parameters, partial derivatives, and preset concentration adjustment parameters until the second ozone characteristic parameters corresponding to all target precursor concentrations are calculated.

[0117] The graph construction unit 705 is used to construct an ozone nonlinear contour graph based on the first ozone characteristic parameter and the precursor concentration at the corresponding target time, the second ozone characteristic parameter and the corresponding target precursor concentration.

[0118] In some embodiments, the photochemical reaction model calculation unit 702 first performs a photochemical reaction simulation calculation based on a combination of a monitored concentration sequence with a meteorological data sequence and a light flux sequence to obtain candidate ozone characteristic parameters at each moment of the monitored concentration sequence; then selects the largest candidate ozone characteristic parameter as the first ozone characteristic parameter corresponding to the monitored concentration sequence, and takes the moment corresponding to the largest candidate ozone characteristic parameter as the target moment; then performs a photochemical reaction simulation calculation based on each precursor simulated concentration sequence, in combination with a meteorological data sequence and a light flux sequence, to obtain other first ozone characteristic parameters corresponding to the target moment.

[0119] In some embodiments, the basic data determination unit 701 obtains the ozone concentration sequence monitored in the same period, then calculates the formaldehyde concentration based on the monitored ozone concentration at each moment of the ozone concentration sequence and the acetylene concentration at the corresponding moment of the acetylene concentration sequence, and constructs the formaldehyde concentration sequence based on the formaldehyde concentration. The photochemical reaction model calculation unit 702 performs photochemical reaction simulation calculation based on the monitored concentration sequence or a precursor simulation concentration sequence, in combination with the meteorological data sequence, the light flux sequence and the formaldehyde concentration sequence.

[0120] In some embodiments, the precursor monitoring concentration sequence includes a nitrogen dioxide concentration sequence. The basic data determination unit 701 calculates the nitrous acid concentration based on the monitored nitrogen dioxide concentration at each moment of the nitrogen dioxide concentration sequence, and constructs a nitrous acid concentration sequence based on the nitrous acid concentration. Accordingly, the photochemical reaction model calculation unit 702 performs a photochemical reaction simulation calculation based on the monitoring concentration sequence or a precursor simulation concentration sequence, in combination with the meteorological data sequence, the light flux sequence and the nitrous acid concentration sequence.

[0121] In some embodiments, the photochemical reaction model calculation unit 702 first converts various volatile organic compound concentration sequences in the monitoring concentration sequence and the precursor simulation concentration sequence into standard mechanism substance concentrations, and determines the standard mechanism substance concentration sequence corresponding to the monitoring concentration sequence and the precursor simulation concentration sequence; then, based on the monitoring concentration sequence or the standard mechanism substance concentration sequence and the volatile organic compound concentration sequence corresponding to a precursor simulation concentration sequence, the photochemical reaction simulation calculation is performed in combination with the meteorological data sequence and the light flux sequence.

[0122] In some embodiments, the basic data determination unit 701 obtains the latitude and longitude information of the monitoring location, and performs illumination orthographic projection solution based on the latitude and longitude information to determine the orthographic projection conversion coefficient; then, based on the sunlight incident intensity and the orthographic projection conversion coefficient, a luminous flux sequence is obtained.

[0123] In some embodiments, the basic data determination unit 701 determines the mixing layer height parameter. The photochemical reaction model calculation unit 702 uses the mixing layer height parameter as a boundary layer constraint, and performs photochemical reaction simulation calculation based on a monitored concentration sequence or a precursor simulated concentration sequence, in combination with a meteorological data sequence and a light flux sequence.

[0124] The embodiment of the present disclosure also provides a computing device for implementing the aforementioned method. Figure 8 Schematic diagram of the structure of the computing device provided by the embodiment of the present disclosure. Figure 8 , which shows a schematic diagram of the structure of a computing device 800 suitable for implementing the embodiments of the present disclosure. Figure 8 The computing device shown is only an example and should not bring any limitation to the functionality and scope of use of the embodiments of the present disclosure.

[0125] like Figure 8 As shown, the computing device 800 may include a processing device (e.g., a central processing unit, a graphics processing unit, etc.) 801, which can perform various appropriate actions and processes according to a program stored in a read-only memory ROM 802 or a program loaded from a storage device 808 to a random access memory RAM 803. In the RAM 803, various programs and data required for the operation of the computing device 800 are also stored. The processing device 801, the ROM 802, and the RAM 803 are connected to each other via a bus 804. An input / output I / O interface 805 is also connected to the bus 804.

[0126] Typically, the following devices may be connected to the I / O interface 805: input devices 805 including, for example, a touch screen, a touch pad, a camera, a microphone, etc.; output devices 807 including, for example, a liquid crystal display (LCD), a speaker, a vibrator, etc.; storage devices 808 including, for example, a magnetic tape, a hard disk, etc.; and communication devices 809. The communication devices 809 may allow the computing device 800 to communicate with other devices wirelessly or by wire to exchange data. Although Figure 8 The computing device 800 having various devices is shown, but it should be understood that it is not required to implement or possess all the devices shown. More or fewer devices may be implemented or possessed instead.

[0127] In particular, according to an embodiment of the present disclosure, the process described above with reference to the flowchart can be implemented as a computer software program. For example, an embodiment of the present disclosure includes a computer program product, which includes a computer program carried on a non-transitory computer-readable medium, and the computer program contains program code for executing the method shown in the flowchart. In such an embodiment, the computer program can be downloaded and installed from a network through a communication device 809, or installed from a storage device 808, or installed from a ROM 802. When the computer program is executed by the processing device 801, the above-mentioned functions defined in the method of the embodiment of the present disclosure are executed.

[0128] It should be noted that the computer-readable medium mentioned above in the present disclosure may be a computer-readable storage medium, a computer-readable signal medium, or any combination of the above two.

[0129] The computer-readable storage medium may be, for example, but not limited to, an electrical, magnetic, optical, electromagnetic, infrared, or semiconductor system, device or device, or any combination thereof. More specific examples of computer-readable storage media may include, but are not limited to, an electrical connection with one or more wires, a portable computer disk, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber, a portable compact disk read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination thereof. In the present disclosure, a computer-readable storage medium may be any tangible medium containing or storing a program that may be used by or in conjunction with an instruction execution system, device, or device.

[0130] Computer readable signal media may include data signals propagated in baseband or as part of a carrier wave, which carry computer readable program code. Such propagated data signals may take a variety of forms, including but not limited to electromagnetic signals, optical signals, or any suitable combination of the above. Computer readable signal media may also be any computer readable medium other than a computer readable storage medium, which may send, propagate, or transmit a program for use by or in conjunction with an instruction execution system, apparatus, or device. The program code contained on the computer readable medium may be transmitted using any suitable medium, including but not limited to: wires, optical cables, RF (radio frequency), etc., or any suitable combination of the above.

[0131] In some embodiments, the client and the server may communicate using any currently known or future developed network protocol such as HTTP (HyperText Transfer Protocol), and may be interconnected with any form or medium of digital data communication (e.g., a communication network). Examples of communication networks include a local area network ("LAN"), a wide area network ("WAN"), an internet (e.g., the Internet), and a peer-to-peer network (e.g., an ad hoc peer-to-peer network), as well as any currently known or future developed network.

[0132] The computer-readable medium may be included in the computing device, or may exist independently without being incorporated into the computing device.

[0133] Computer program code for performing the operations of the present disclosure may be written in one or more programming languages ​​or a combination thereof, including, but not limited to, object-oriented programming languages, such as Java, Smalltalk, C++, and conventional procedural programming languages, such as "C" or similar programming languages. The program code may be executed entirely on the tester's computer, partially on the tester's computer, as a separate software package, partially on the tester's computer and partially on a remote computer, or entirely on a remote computer or server. In cases involving a remote computer, the remote computer may be connected to the tester's computer through any type of network, including a local area network (LAN) or a wide area network (WAN), or may be connected to an external computer (e.g., through the Internet using an Internet service provider).

[0134] The flow chart and block diagram in the accompanying drawings illustrate the possible architecture, function and operation of the system, method and computer program product according to each embodiment of the present disclosure. In this regard, each box in the flow chart or block diagram can represent a module, a program segment or a part of a code, and the module, the program segment or a part of the code contains one or more executable instructions for realizing the specified logical function. It should also be noted that in some implementations as replacements, the functions marked in the box can also occur in a sequence different from that marked in the accompanying drawings. For example, two boxes represented in succession can actually be executed substantially in parallel, and they can sometimes be executed in the opposite order, depending on the functions involved. It should also be noted that each box in the block diagram and / or flow chart, and the combination of the boxes in the block diagram and / or flow chart can be implemented with a dedicated system based on hardware that performs the specified function or operation, or can be implemented with a combination of dedicated hardware and computer instructions.

[0135] The units involved in the embodiments described in the present disclosure may be implemented by software or by hardware. The name of the unit does not constitute a limitation on the unit itself in some cases. The functions described above in this article may be performed at least in part by one or more hardware logic components. For example, without limitation, exemplary types of hardware logic components that can be used include: field programmable gate array (FPGA), application specific integrated circuit (ASIC), application specific standard product (ASSP), system on chip (SOC), complex programmable logic device (CPLD), etc.

[0136] The above are only specific embodiments of the present disclosure, so that those skilled in the art can understand or implement the present disclosure. Various modifications to these embodiments will be apparent to those skilled in the art, and the general principles defined herein can be implemented in other embodiments without departing from the spirit or scope of the present disclosure. Therefore, the present disclosure will not be limited to these embodiments herein, but will conform to the widest scope consistent with the principles and novel features disclosed herein.

Claims

1. A method for constructing nonlinear contour curves of ozone. It is characterized in that include: Acquire precursor monitoring concentration sequences, meteorological data sequences and light flux sequences monitored in the same period, and obtain multiple groups of adjacent precursor simulation concentration sequences based on the monitoring concentration sequences and preset concentration adjustment values; the monitoring concentration sequences and the simulation concentration sequences both include volatile organic compound concentration sequences and nitrogen oxide concentration sequences; Based on the monitoring concentration sequence or a group of the precursor simulation concentration sequences, respectively, in combination with the meteorological data sequence and the light flux sequence, a photochemical reaction simulation calculation is performed to obtain a plurality of first ozone characteristic parameters corresponding to the target time; the first ozone characteristic parameter is ozone generation potential or ozone concentration; Perform derivative calculation based on each first ozone characteristic parameter and the preset concentration adjustment value to determine partial derivatives of various orders; Based on the partial derivatives of each order, the preset concentration adjustment value, and the target moment precursor concentration corresponding to each of the first ozone characteristic parameters, the second ozone characteristic parameter of the neighboring other precursor concentrations is calculated by summing the Taylor expansion and ignoring the high-order expansion terms, and the calculation is repeated based on the neighboring other precursor concentrations and the corresponding second ozone characteristic parameter, the partial derivative, and the preset concentration adjustment parameter until the second ozone characteristic parameters corresponding to all target precursor concentrations are calculated; Based on the first ozone characteristic parameter and the precursor concentration at the corresponding target time, the second ozone characteristic parameter and the corresponding target precursor concentration, an ozone nonlinear contour curve diagram is constructed.

2. The method according to claim 1, It is characterized in that The photochemical reaction simulation calculation is performed based on the monitoring concentration sequence or one of the precursor simulation concentration sequences, in combination with the meteorological data sequence and the light flux sequence, to obtain a plurality of first ozone characteristic parameters corresponding to the target time, including: Performing a photochemical reaction simulation calculation based on the combination of the monitoring concentration sequence, the meteorological data sequence and the light flux sequence to obtain the candidate ozone characteristic parameters at each moment of the monitoring concentration sequence; Selecting the largest candidate ozone characteristic parameter as the first ozone characteristic parameter corresponding to the monitoring concentration sequence, and taking the time corresponding to the largest candidate ozone characteristic parameter as the target time; Based on each of the precursor simulated concentration sequences, the simulated photochemical reaction is calculated in combination with the meteorological data sequence and the light flux sequence to obtain other first ozone characteristic parameters corresponding to the target time.

3. The method according to claim 1 or 2, It is characterized in that The precursor monitoring concentration sequence includes an acetylene concentration sequence; the method further includes: obtaining an ozone concentration sequence monitored during the same period; Before performing photochemical reaction simulation calculation based on the monitoring concentration sequence or one of the precursor simulation concentration sequences, in combination with the meteorological data sequence and the light flux sequence, the method further comprises: Calculating formaldehyde concentration based on the monitored ozone concentration at each moment of the ozone concentration sequence and the acetylene concentration at the corresponding moment of the acetylene concentration sequence, and constructing a formaldehyde concentration sequence based on the formaldehyde concentration; The photochemical reaction simulation calculation based on the monitored concentration sequence or one of the precursor simulated concentration sequences, combined with the meteorological data sequence and the light flux sequence, includes: performing photochemical reaction simulation calculation based on the monitored concentration sequence or one of the precursor simulated concentration sequence, combined with the meteorological data sequence, the light flux sequence and the formaldehyde concentration sequence.

4. The method according to claim 1 or 2, It is characterized in that The precursor monitoring concentration sequence includes a nitrogen dioxide concentration sequence; Before performing photochemical reaction simulation calculation based on the monitoring concentration sequence or one of the precursor simulation concentration sequences, in combination with the meteorological data sequence and the light flux sequence, the method further comprises: Calculating nitrous acid concentration based on the monitored nitrogen dioxide concentration at each moment of the nitrogen dioxide concentration sequence, and constructing a nitrous acid concentration sequence based on the nitrous acid concentration; The photochemical reaction simulation calculation based on the monitored concentration sequence or one of the precursor simulated concentration sequences, combined with the meteorological data sequence and the light flux sequence, includes: performing photochemical reaction simulation calculation based on the monitored concentration sequence or one of the precursor simulated concentration sequence, combined with the meteorological data sequence, the light flux sequence and the nitrite concentration sequence.

5. The method according to claim 1 or 2, It is characterized in that Before performing photochemical reaction simulation calculation based on the monitoring concentration sequence or one of the precursor simulation concentration sequences, in combination with the meteorological data sequence and the light flux sequence, the method further includes: Performing standard mechanism substance concentration conversion on various volatile organic compound concentration sequences in the monitoring concentration sequence and the precursor simulation concentration sequence, and determining a standard mechanism substance concentration sequence corresponding to the monitoring concentration sequence and the precursor simulation concentration sequence; The photochemical reaction simulation calculation is performed based on the monitored concentration sequence or one of the precursor simulation concentration sequences, in combination with the meteorological data sequence and the light flux sequence, including: the photochemical reaction simulation calculation is performed based on the standard mechanism substance concentration sequence and volatile organic compound concentration sequence corresponding to the monitored concentration sequence or one of the precursor simulation concentration sequence, in combination with the meteorological data sequence and the light flux sequence.

6. The method according to claim 1 or 2, It is characterized in that The method for obtaining the light flux sequence comprises: Obtaining the latitude and longitude information of the monitoring location, and performing illumination orthographic projection calculation based on the latitude and longitude information to determine the orthographic projection conversion coefficient; The luminous flux sequence is obtained based on the incident intensity of sunlight and the orthographic projection conversion coefficient.

7. The method according to claim 1 or 2, It is characterized in that Before performing photochemical reaction simulation calculation based on the monitoring concentration sequence or one of the precursor simulation concentration sequences, respectively, in combination with the meteorological data sequence and the light flux sequence, the method further comprises: determining a mixing layer height parameter; The performing photochemical reaction simulation calculation based on the monitoring concentration sequence or one of the precursor simulation concentration sequences, in combination with the meteorological data sequence and the light flux sequence, comprises: The mixing layer height parameter is used as a boundary layer constraint, and based on the monitored concentration sequence or a precursor simulated concentration sequence, a photochemical reaction simulation calculation is performed in combination with the meteorological data sequence and the light flux sequence.

8. A device for constructing nonlinear contour graphs of ozone, It is characterized in that include: A basic data determination unit, used to obtain a precursor monitoring concentration sequence, a meteorological data sequence and a light flux sequence monitored in the same period, and obtain multiple groups of adjacent precursor simulation concentration sequences based on the monitoring concentration sequence and a preset concentration adjustment value; the monitoring concentration sequence and the simulation concentration sequence both include a volatile organic compound concentration sequence and a nitrogen oxide concentration sequence; a photochemical reaction model calculation unit, configured to perform photochemical reaction simulation calculations based on the monitoring concentration sequence or a group of precursor simulation concentration sequences, in combination with the meteorological data sequence and the light flux sequence, to obtain a plurality of first ozone characteristic parameters corresponding to a target time; the first ozone characteristic parameter is ozone generation potential or ozone concentration; A derivative calculation unit, used to perform derivative calculation based on each first ozone characteristic parameter and the preset concentration adjustment value to determine partial derivatives of various orders; a point value data calculation unit, for calculating the second ozone characteristic parameter of the neighboring other precursor concentrations in a manner of Taylor expansion summation and ignoring high-order expansion terms based on the partial derivatives of each order, the preset concentration adjustment value, and the target moment precursor concentration corresponding to each of the first ozone characteristic parameters, and repeating the calculation based on the neighboring other precursor concentrations and the corresponding second ozone characteristic parameter, the partial derivative, and the preset concentration adjustment parameter until the second ozone characteristic parameters corresponding to all target precursor concentrations are calculated; The curve graph construction unit is used to construct an ozone nonlinear contour curve graph based on the first ozone characteristic parameter and the precursor concentration at the corresponding target time, the second ozone characteristic parameter and the corresponding target precursor concentration.

9. A computing device, It is characterized in that The method comprises a processor and a memory, wherein the memory is used to store a computer program; when the computer program is loaded by the processor, the processor executes the method for constructing a nonlinear isovalue curve diagram of ozone according to any one of claims 1 to 7.

10. A computer-readable storage medium, It is characterized in that The storage medium stores a computer program, and when the computer program is executed by the processor, the processor implements the method for constructing an ozone nonlinear contour curve map as described in any one of claims 1 to 7.

Images