Method and apparatus, device, and medium for constructing ozone non-linear equivalent value curve diagram

Abstract

Embodiments of the present disclosure provide a method and apparatus, device, and medium for constructing an ozone non-linear contour graph. The method for constructing an ozone non-linear contour graph includes: obtaining a precursor monitoring concentration sequence monitored during the same period, and obtaining multiple sets of adjacent precursor simulated concentration sequences based on the monitoring concentration sequence; respectively obtaining multiple first ozone characteristic parameters based on the monitoring concentration sequence or a set of precursor simulated concentration sequences; performing derivative calculations based on each of the first ozone characteristic parameters to determine partial derivatives of each order; calculating second ozone characteristic parameters corresponding to all target precursor concentrations based on the partial derivatives of each order; constructing an ozone non-linear contour graph based on the first ozone characteristic parameters and the precursor concentrations at the corresponding target moments, and the second ozone characteristic parameters and the corresponding target precursor concentrations. Based on the method of this solution, the ozone non-linear contour curve can be obtained only through less data calculation, and thus the calculation cost can be reduced.

Description

Technical Field

[0001] The present disclosure relates to the technical field of environmental protection, and particularly relates to a method and device, equipment, and medium for constructing an ozone non-linear equivalent curve graph. Background Art

[0002] Ozone and pollutants in the atmospheric environment are secondary products, and there is a very complex non-linear relationship between their concentrations and the precursor volatile organic compounds (VOCs) and nitrogen oxides (NOx). Blindly reducing the concentration of ozone precursors often not only fails to reduce the ozone concentration but may instead cause a rebound in the ozone concentration. In the industry, the analysis of the impact of precursor changes on ozone is generally referred to 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 an ozone non-linear equivalent curve based on observed model simulation data.

[0003] Drawing an ozone non-linear equivalent curve based on observed model simulation data generally uses the brute-force method. This method calculates the ozone formation potential at different concentrations by continuously changing the ratio of precursor concentrations, thereby drawing an ozone non-linear equivalent curve graph for the analysis of ozone and precursor sensitivity. However, this method requires continuously changing the precursor concentration values for simulation, resulting in a large amount of calculation and long calculation time. Summary of the Invention

[0004] To solve the problem of large calculation amount and long calculation time in constructing an ozone non-linear equivalent curve graph using the brute-force method, embodiments of the present disclosure provide a new method and device, equipment, and medium for constructing an ozone non-linear equivalent curve graph.

[0005] In a first aspect, embodiments of the present disclosure provide a method for constructing an ozone non-linear equivalent curve graph, including:

[0006] Obtaining a sequence of monitored precursor concentrations, a sequence of meteorological data, and a sequence of light fluxes monitored during the same period, and obtaining multiple sets of adjacent simulated precursor concentration sequences based on the monitored concentration sequence and a preset concentration adjustment value; both the monitored concentration sequence and the simulated concentration sequence include a sequence of volatile organic compound concentrations and a sequence of nitrogen oxide concentrations;

[0007] Respectively, based on the monitored concentration sequence or a set of the simulated precursor concentration sequences, combining with the meteorological data sequence and the light flux sequence for photochemical reaction simulation calculation to obtain multiple first ozone characteristic parameters corresponding to the target time; the first ozone characteristic parameter is the ozone formation potential or the ozone concentration;

[0008] Derivative calculations are performed based on each first ozone characteristic parameter and the preset concentration adjustment value to determine the partial derivatives of each order;

[0009] Based on the partial derivatives of each order, the preset concentration adjustment value, and the precursor concentrations at the target time corresponding to each of the first ozone characteristic parameters, the second ozone characteristic parameters of the neighboring other precursor concentrations are calculated by summing according to 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, the partial derivatives, and the preset concentration adjustment parameters until the second ozone characteristic parameters corresponding to all target precursor concentrations are calculated;

[0010] Based on the first ozone characteristic parameter and the precursor concentration at the corresponding target time, and the second ozone characteristic parameter and the corresponding target precursor concentration, an ozone non-linear isoconcentration curve graph is constructed.

[0011] Optionally, the performing of the photochemical reaction simulation calculations by combining the monitoring concentration sequence or one of the precursor simulated concentration sequences with the meteorological data sequence and the light flux sequence respectively to obtain multiple first ozone characteristic parameters corresponding to the target time includes:

[0012] Performing photochemical reaction simulation calculations by combining the monitoring concentration sequence with the meteorological data sequence and the light flux sequence to obtain the candidate ozone characteristic parameters at each moment of the monitoring concentration sequence;

[0013] Selecting the maximum candidate ozone characteristic parameter as the first ozone characteristic parameter corresponding to the monitoring concentration sequence, and taking the moment corresponding to the maximum candidate ozone characteristic parameter as the target time;

[0014] Performing photochemical reaction simulation calculations by combining each of the precursor simulated concentration sequences with the meteorological data sequence and the light flux sequence respectively to obtain other first ozone characteristic parameters corresponding to the target time.

[0015] Optionally, the precursor monitoring concentration sequence includes an acetylene concentration sequence; the method further includes: obtaining an ozone concentration sequence monitored in the same period;

[0016] Before performing the photochemical reaction simulation calculations by combining the monitoring concentration sequence or one of the precursor simulated concentration sequences with the meteorological data sequence and the light flux sequence respectively, the method further includes:

[0017] Calculating 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 constructing a formaldehyde concentration sequence based on the formaldehyde concentration;

[0018] Performing photochemical reaction simulation calculations by combining the monitored concentration sequence or one of the precursor simulation concentration sequences with the meteorological data sequence and the light flux sequence includes: performing photochemical reaction simulation calculations by combining the monitored concentration sequence or one of the precursor simulation concentration sequences with the meteorological data sequence, the light flux sequence, and the formaldehyde concentration sequence.

[0019] Optionally, the precursor monitored concentration sequence includes a nitrogen dioxide concentration sequence;

[0020] Before performing photochemical reaction simulation calculations by combining the monitored concentration sequence or one of the precursor simulation concentration sequences with the meteorological data sequence and the light flux sequence, the method further includes:

[0021] 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;

[0022] Performing photochemical reaction simulation calculations by combining the monitored concentration sequence or one of the precursor simulation concentration sequences with the meteorological data sequence and the light flux sequence includes: performing photochemical reaction simulation calculations by combining the monitored concentration sequence or one of the precursor simulation concentration sequences with the meteorological data sequence, the light flux sequence, and the nitrous acid concentration sequence.

[0023] Optionally, before performing photochemical reaction simulation calculations by combining the monitored concentration sequence or one of the precursor simulation concentration sequences with the meteorological data sequence and the light flux sequence respectively, it further includes:

[0024] Performing standard mechanism substance concentration conversion on the various volatile organic compound concentration sequences in the monitored concentration sequence and the precursor simulation concentration sequence to determine the standard mechanism substance concentration sequences corresponding to the monitored concentration sequence and the precursor simulation concentration sequence;

[0025] Performing photochemical reaction simulation calculations by combining the monitored concentration sequence or one of the precursor simulation concentration sequences with the meteorological data sequence and the light flux sequence respectively includes: performing photochemical reaction simulation calculations by combining the standard mechanism substance concentration sequences and the volatile organic compound concentration sequences corresponding to the monitored concentration sequence or one of the precursor simulation concentration sequences with the meteorological data sequence and the light flux sequence.

[0026] Optionally, the method for obtaining the light flux sequence includes: obtaining the longitude and latitude information of the monitoring location, and performing orthographic projection solution of illumination based on the longitude and latitude information to determine the orthographic projection conversion coefficient;

[0027] Based on the solar light incident intensity and the positive projection conversion coefficient, the light flux sequence is obtained.

[0028] Optionally, before performing the photochemical reaction simulation calculation by combining the monitoring concentration sequence or one of the precursor simulated concentration sequences with the meteorological data sequence and the light flux sequence respectively, the method further includes: determining the mixed layer height parameter;

[0029] The performing the photochemical reaction simulation calculation by combining the monitoring concentration sequence or one of the precursor simulated concentration sequences with the meteorological data sequence and the light flux sequence respectively includes:

[0030] Using the mixed layer height parameter as the boundary layer constraint, performing the photochemical reaction simulation calculation by combining the monitoring concentration sequence or one of the precursor simulated concentration sequences with the meteorological data sequence and the light flux sequence.

[0031] In a second aspect, an ozone non - linear contour graph construction device provided by an embodiment of the present disclosure includes:

[0032] A basic data determination unit, configured to obtain a precursor monitoring concentration sequence, a meteorological data sequence, and a light flux sequence monitored at the same time period, and obtain multiple groups of adjacent precursor simulated concentration sequences based on the monitoring concentration sequence and a preset concentration adjustment value; both the monitoring concentration sequence and the simulated concentration sequence include a volatile organic compound concentration sequence and a nitrogen oxide concentration sequence;

[0033] A photochemical reaction model calculation unit, configured to perform photochemical reaction simulation calculations by combining the monitoring concentration sequence or a group of the precursor simulated concentration sequences with the meteorological data sequence and the light flux sequence respectively, to obtain multiple first ozone characteristic parameters corresponding to the target moment; the first ozone characteristic parameter is the ozone formation potential or the ozone concentration;

[0034] A derivative calculation unit, configured to perform derivative calculations based on each first ozone characteristic parameter and the preset concentration adjustment value to determine each order partial derivative;

[0035] A point value data calculation unit, configured to calculate second ozone characteristic parameters of adjacent other precursor concentrations in the manner of summing according to the Taylor expansion and ignoring the high - order expansion terms based on the each order partial derivative, the preset concentration adjustment value, and the precursor concentrations corresponding to each of the first ozone characteristic parameters at the target moment, and repeat the calculation based on the adjacent other precursor concentrations and the corresponding second ozone characteristic parameters, the partial derivative, and the preset concentration adjustment parameter until the second ozone characteristic parameters corresponding to all target precursor concentrations are calculated;

[0036] A curve graph construction unit is configured to construct a non-linear ozone 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.

[0037] In a third aspect, an embodiment of the present disclosure provides a computing device, including a processor and a memory, where the memory is configured to store a computer program; when the computer program is loaded by the processor, the processor is caused to execute the ozone non-linear contour graph construction method as described above.

[0038] In a fourth aspect, an embodiment of the present disclosure provides a computer-readable storage medium storing a computer program, which, when executed by a processor, causes the processor to implement the ozone non-linear contour graph construction method as described above.

[0039] The solution of the embodiment of the present disclosure considers that after obtaining the first ozone characteristic parameters of multiple concentration points through photochemical simulation with only a relatively small number of precursor sequences, calculating the derivatives of each order using the foregoing first ozone characteristic parameters, and then obtaining the second ozone characteristic parameters corresponding to other precursor concentrations by means of Taylor expansion using the foregoing derivatives of each order. Subsequently, a non-linear ozone contour graph can be constructed based on the first ozone characteristic parameters and the second ozone characteristic parameters. Based on the method of this solution, a non-linear ozone contour graph can be obtained only through relatively few data calculations, thereby reducing the calculation cost. Description of the Drawings

[0040] The drawings herein are incorporated into the specification and form a part of the specification, showing embodiments consistent with the present disclosure, and are used together with the specification to explain the principles of the present disclosure.

[0041] To more clearly illustrate the technical solutions in the embodiments of the present disclosure or the prior art, the following will briefly introduce the drawings required for use in the description of the embodiments or the prior art. Obviously, for those of ordinary skill in the art, other drawings can also be obtained based on these drawings without creative efforts, where

[0042] Figure 1 is a flowchart of the ozone non-linear contour graph construction method provided by an embodiment of the present disclosure;

[0043] Figure 2 is a non-linear ozone concentration contour graph constructed using the solution of the embodiment of the present disclosure;

[0044] Figure 3 is a non-linear ozone concentration contour graph constructed using the brute force method;

[0045] Figure 4is the non-linear isopleth curve diagram of ozone formation potential constructed by the solution of the present disclosure embodiment;

[0046] Figure 5 is the non-linear isopleth curve diagram of ozone formation potential constructed by the brute force method;

[0047] Figure 6 is the schematic structural diagram of the architecture based on in the specific implementation of the present disclosure of the foregoing method;

[0048] Figure 7 is the schematic diagram of the ozone non-linear isopleth curve construction device provided by the embodiment of the present disclosure;

[0049] Figure 8 is the schematic structural diagram of the computing device provided by the embodiment of the present disclosure. Detailed Description of the Invention

[0050] Embodiments of the present disclosure will be described in more detail below with reference to the accompanying drawings. Although some embodiments of the present disclosure are shown in the drawings, it should be understood that the present disclosure can be implemented in various forms and should not be construed as limited to the embodiments set forth herein. On the contrary, these embodiments are provided to more thoroughly and completely understand the present disclosure. It should be understood that the drawings and embodiments of the present disclosure are only for exemplary purposes and are not used to limit the protection scope of the present disclosure.

[0051] As used herein, the term "comprising" and its variations are open-ended, i.e., "including but not limited to". The term "based on" is "at least partially based on". The term "one embodiment" means "at least one embodiment"; the term "another embodiment" means "at least one additional embodiment"; the term "some embodiments" means "at least some embodiments". The relevant definitions of other terms will be given in the following description. In this document, relational terms such as "first" and "second" are only used to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply these embodiments.

[0052] To solve the problem that the existing brute force method for calculating and plotting the ozone non-linear isopleth curve has a large amount of calculation and is very time-consuming, the embodiment of the present disclosure provides a new method for constructing the ozone non-linear isopleth curve diagram. The method for constructing the ozone non-linear isopleth curve diagram provided by the embodiment of the present disclosure is executed by a computing device.

[0053] Figure 1 is the flowchart of the method for constructing the ozone non-linear isopleth curve provided by the embodiment of the present disclosure. As Figure 1 shown, the method for constructing the ozone non-linear isopleth curve provided by the embodiment of the present disclosure includes S110 - S150.

[0054] S110: Obtain the precursor monitoring concentration sequence, meteorological data sequence, and light flux sequence monitored during the same period, and obtain multiple groups of adjacent precursor simulation concentration sequences based on the monitoring concentration sequence and the preset concentration adjustment value.

[0055] Considering that ozone in the air is caused by photochemical reactions and light is a prerequisite for photochemical reactions, the same period mentioned in the embodiments of the present disclosure is the period with sunlight irradiation. In specific implementation, the same period can be determined according to the sunrise and sunset times of a certain area.

[0056] The precursor monitoring concentration sequence is a monitoring concentration sequence obtained by arranging the concentrations of a certain ozone-forming precursor monitored at the monitoring site during the same period. The precursor concentration monitoring sequence is a sequence obtained by sorting the precursor concentrations in chronological order.

[0057] Since the precursors for ozone formation include volatile organic compounds and nitrogen oxides, the precursor concentration sequence in the embodiments of the present disclosure includes an organic compound concentration sequence and a nitrogen oxide concentration sequence. In specific implementation, the types of volatile organic compounds corresponding to the aforementioned organic compound concentration sequence are determined, and they can be various possible organic compounds causing regional environmental pollution, which are specifically determined according to environmental pollution monitoring experience. Similarly, the types of volatile nitrogen oxides corresponding to the nitrogen oxide concentration sequence are also determined, and they can also be determined according to environmental monitoring experience. In specific implementation, volatile organic compounds can include carbon monoxide, various volatile alkane compounds, benzene compounds, and volatile oxides of the aforementioned various compounds; nitrogen oxides can be nitric oxide, nitrogen dioxide, etc.

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

[0059] Table 1 Table of Monitored Concentration Sequences of Volatile Organic Compounds

[0060]

[0061]

[0062] Table 2 Table of Nitrogen Oxide Concentration Sequences

[0063]

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

[0065] Table 3 Temperature sequence table for a certain period

[0066]

[0067] A light flux sequence is a sequence composed of light fluxes monitored or determined by data conversion during a certain period in chronological order. The aforementioned light flux is the radiant flux irradiating a certain volume of air from all directions. The light flux is related to the incident intensity of the sun. The light flux corresponding to the light flux sequence in the embodiments of the present disclosure is the sequence corresponding to the light of wavelengths that can affect the photochemical reaction process.

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

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

[0070] S112: Based on the solar incident intensity and the orthographic projection conversion coefficient, obtain the light flux sequence.

[0071] According to geographical knowledge, the longitude and latitude information and time point information of the monitoring location affect the light intensity received per unit area. When the solar incident intensity is determined, it is necessary to determine the solar incident angle in order to determine the orthographic projection light intensity at the corresponding position according to the incident angle. According to the aforementioned logic, the embodiment solution of the present disclosure performs light projection orthographic transformation 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.

[0072] After determining the orthographic projection conversion coefficient, the orthographic projection conversion coefficient can then be multiplied by the solar incident intensity to obtain the light intensity at the corresponding time point. Since the energy ratio of various wavelengths in the solar spectrum is basically determined, the light fluxes of different wavelengths at each moment can be determined according to the aforementioned method. Sorting the aforementioned light fluxes in chronological order can obtain the light flux sequence.

[0073] In the embodiments of the present disclosure, after obtaining the precursor monitoring concentration sequence, numerical addition and subtraction operations are performed based on each sequence in the foregoing set of precursor monitoring concentrations and a preset concentration adjustment value to obtain a precursor simulation concentration sequence similar to the precursor monitoring concentration sequence. As its name implies, the precursor simulation concentration sequence is not a real sequence but a hypothetical sequence. The preset concentration adjustment value is a concentration adjustment value determined according to actual experience and the required numerical accuracy. For example, under normal circumstances, the concentration change range of a certain pollutant is 0 - 300, and 16 points need to be arranged within this concentration change range for this pollutant. The corresponding preset concentration adjustment value can be 20 (this is only an example here, so the dimension problem is not considered).

[0074] It should be noted that in the embodiments of the present disclosure, when performing numerical addition and subtraction operations according to a set of precursor monitoring concentration sequences and a preset concentration adjustment value to obtain a set of precursor simulation concentration sequences, the preset concentration adjustment value is simultaneously added to or subtracted from the organic matter concentration sequences corresponding to different organic matters to obtain an organic matter concentration adjustment sequence with the same law of change, and the preset concentration adjustment value is simultaneously added to or subtracted from different nitrogen oxide concentration sequences to obtain a nitrogen oxide concentration sequence with the same law of change.

[0075] In some embodiments of the present disclosure, the number of multiple sets of adjacent precursor simulation concentration sequences determined based on the precursor monitoring concentration sequence is 8 sets. The foregoing 8 sets of adjacent precursor simulation concentration sequences can be centered on the precursor monitoring concentration sequence or located on one side of the precursor monitoring concentration sequence. Assuming that each concentration sequence is a point, the foregoing 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 embodiments of the present disclosure, the goal of determining 8 sets of adjacent precursor concentration sequences is to obtain the first-order partial derivative and the second-order partial derivative for subsequent summation calculations based on the Taylor expansion method. In other embodiments, if only the first-order partial derivative is required for the summation calculation based on the Taylor expansion, the number of adjacent precursor simulation concentration sequences only needs to be 3 sets. And if the third-order partial derivative or even higher-order partial derivatives are required for the summation calculation based on the Taylor expansion, more sets of adjacent precursor concentration sequences are needed.

[0076] In specific implementation, while monitoring the precursor monitoring concentration sequence, the concentration of pollutants such as sulfur dioxide may also be monitored. Because it has a similar form to nitrogen oxides, in practical applications, it can also be treated as nitrogen oxides and processed according to the treatment method of nitrogen oxides to obtain the corresponding simulation concentration sequence.

[0077] S120: Based on the monitored concentration sequence or a set of simulated concentration sequences of precursors respectively, combine them with the meteorological data sequence and the light flux sequence to perform photochemical reaction simulation calculations, and obtain multiple first ozone characteristic parameters corresponding to the target time.

[0078] After obtaining the aforementioned monitored concentration sequence or a simulated concentration sequence of a precursor, and after obtaining the meteorological data sequence and the light flux sequence, subsequent photochemical simulation calculations can be carried out according to the principle of photochemical reactions.

[0079] In the embodiments of the present disclosure, for the aforementioned monitored concentration sequence of precursors and each set of simulated concentration sequences of precursors, respectively combine them with the meteorological data sequence and the light flux sequence to obtain multiple first ozone characteristic parameters corresponding to the target time.

[0080] Taking the combination of the monitored concentration sequence of precursors, the meteorological data sequence and the light flux sequence for photochemical reaction simulation calculations as an example, it performs various chemical reaction simulations based on the monitored concentration of precursors, temperature, humidity data and light flux at the current time to determine the change value of ozone concentration at the first time. And based on the historical cumulative value of ozone concentration and the change value of ozone concentration, determine a new cumulative value of ozone concentration, and use the new cumulative value of ozone concentration, the monitored concentration of precursors, temperature, humidity data and light flux at the next time as conditions to calculate the cumulative value of ozone concentration at the next time until the first ozone characteristic parameter at the target time is obtained. Similarly, by performing photochemical simulation calculations with each set of adjacent simulated concentration sequences of precursors, the meteorological data sequence and the light flux sequence, the corresponding first ozone characteristic parameters can also be obtained.

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

[0082] S01: Perform standard mechanism substance concentration conversion on the concentration sequences of various volatile organic compounds in the monitored concentration sequence and the simulated concentration sequence of precursors to determine the standard mechanism substance concentration sequences corresponding to the monitored concentration sequence and the simulated concentration sequence of precursors.

[0083] Table 4 is a conversion coefficient table of various volatile organic compounds and standard mechanism substances in carbon bond 05, where each row represents a volatile organic compound and each column represents a standard mechanism substance. The presence of a value in the corresponding row and column indicates that a unit amount of a volatile organic compound can be converted into the corresponding amount of the standard mechanism substance. It should be noted that a volatile organic compound may correspond to multiple standard mechanism substances, that is, it is converted into a combination of various standard mechanism substances to be able to characterize its actual role in photochemical reactions.

[0084] Table 4 Conversion Coefficient Table of Various Volatile Organic Compounds and Standard Mechanism Compounds in Carbon Bond 05

[0085]

[0086]

[0087] In specific implementation, the computing device can perform conversion based on the conversion coefficients in Table 4 to obtain the corresponding concentration sequence of standard mechanism compounds. Specifically, after the volatile organic compounds in the concentration sequence are converted according to the conversion coefficient table of standard mechanism compounds, the concentration data of various standard mechanism compounds are obtained. Subsequently, the concentration data of the same standard mechanism compound at the same moment are added together, and the final concentration data of the standard mechanism compound can be obtained. The foregoing conversion method is expressed by the formula as , where is the conversion concentration of the th standard mechanism compound at moment, is the conversion concentration of the th volatile organic compound at moment, is the conversion coefficient between the volatile organic compound and the th standard mechanism compound, which can be found in Table 3.

[0088] After obtaining the final concentration data of the standard mechanism compound, arranging them in chronological order can obtain the concentration sequence of the standard mechanism compound.

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

[0090] In specific implementation, the computing device can perform hardening calculations based on the concentrations of various standard mechanisms and the concentration of nitrogen oxides according to a pre-determined reaction model, determine the reaction rates of each step in the photochemical reaction simulation process, and obtain the first ozone characteristic parameter based on the foregoing reaction rates and some 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 reactant species B. The aforementioned reactant species A and reactant species B 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 this reaction, respectively, and in most cases, a and b are both 1; , is the gas constant, is the activation energy, is the reaction temperature, is the pre-exponential factor.

[0091] In the embodiments of the present disclosure, the aforementioned first ozone characteristic parameter is the ozone formation potential or the ozone concentration. In specific implementations, the ozone formation potential and the ozone concentration can be calculated respectively for subsequent construction of two different types of non-linear isometric curves.

[0092] The following analyzes how to calculate the ozone formation potential. In specific implementations, since the formation of ozone is directly related to nitric oxide, can be used to obtain the ozone formation potential, where is the ozone formation potential, is the net generation of ozone, is net consumption.

[0093] In specific implementations, ,

[0094] .

[0095] Table 5 is the chemical reaction parameter table of carbon bond 05. The parameters in the aforementioned formula can be obtained by looking up in Table 4.

[0096] Table 5 Chemical Reaction Parameter Table of Carbon Bond 05

[0097]

[0098]

[0099]

[0100]

[0101]

[0102]

[0103] In Table 5 mentioned above, 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 are respectively photolysis reactions, where the photolysis reaction rate constant , where is the light flux, is the absorption cross-sectional area for electromagnetic waves with a wavelength of , 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 .

[0104] The following is an analysis of how to calculate the ozone concentration. In specific implementation, the ozone concentration can be obtained by integrating the ozone generation rate over time, specifically , where is the ozone generation rate, is the initial concentration of ozone, which is generally set to 0. The ozone generation rate , where is the net ozone generation rate, is the net ozone consumption rate.

[0105]

[0106] .

[0107] In specific implementation, due to the monitoring mechanism, some precursors are lacking in the monitored concentration sequence of precursors, and the concentration sequence of such precursors in the photochemical reaction cannot be determined. For example, formaldehyde is not easily soluble in water, and its actual content is not easily measured by monitoring equipment. Another example is that nitrous acid is also difficult to measure because it is soluble in water.

[0108] To solve the problem mentioned above, in some embodiments of the present disclosure, before performing the aforementioned S120, the following S02 or S03 can also be performed.

[0109] 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.

[0110] To ensure the execution of S02, the previous monitored concentration sequence of precursors should include the acetylene concentration sequence, and in addition, the ozone concentration sequence monitored during the same period should also be obtained.

[0111] After obtaining the acetylene concentration sequence and the ozone concentration sequence, the acetylene concentration and the ozone concentration at the same corresponding moment can be used to calculate the corresponding formaldehyde concentration according to wherein is the formaldehyde concentration,[[]]END]] is the acetylene concentration,[[]]END]] is the ozone concentration. Subsequently, the formaldehyde concentration can be sorted according to time to obtain the formaldehyde concentration sequence.[[]]END]]

[0112] On the premise of executing the aforementioned S02, the aforementioned S120 is specifically S122 as follows.[[]]END]]

[0113] S122: Based on the monitored concentration sequence or a precursor simulated concentration sequence, combined with the meteorological data sequence, the light flux sequence and the formaldehyde concentration sequence, perform photochemical reaction simulation calculations.[[]]END]]

[0114] The specific implementation process of S122 is as described above and will not be repeated here. It should be noted that the formaldehyde concentration sequence is a sequence of the concentration of the corresponding detected substance obtained from the monitored concentration of the precursor, rather than a sequence directly for the precursor simulated concentration sequence. Considering that the formaldehyde concentration is small and has little impact on the photochemical reaction, it is directly combined with the precursor simulated concentration sequence here 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 simulated concentration sequence, and then execute S121 correspondingly.[[]]END]]

[0115] S03: Calculate the nitrous acid concentration based on the monitored nitrous oxide concentration at each moment of the nitrous oxide concentration sequence, and construct a nitrous acid concentration sequence based on the nitrous acid concentration.[[]]END]]

[0116] In specific implementation, the computing device calculates according to to obtain the nitrous acid concentration wherein is the nitrous oxide concentration, and MR is the relative molecular mass of nitrous oxide. After obtaining the nitrous acid concentration at each moment, sort the nitrous acid concentrations in chronological order to obtain the nitrous acid concentration sequence.[[]]END]]

[0117] On the premise of executing the aforementioned S03, the aforementioned S120 is specifically S123 as follows.[[]]END]]

[0118] S123: Based on the monitored concentration sequence or a precursor simulated concentration sequence, combined with the meteorological data sequence, the light flux sequence and the nitrous acid concentration sequence, perform photochemical reaction simulation calculations.[[]]END]]

[0119] The specific implementation process of S123 is as described above, and will not be repeated here. It should be noted that the aforementioned nitrous acid concentration sequence is a sequence corresponding to the concentration of the analyte obtained from the monitored concentration of the precursor, rather than a sequence directly targeting the simulated concentration sequence of the precursor. Considering that the nitrous acid concentration is relatively small and has a relatively small impact on the photochemical reaction, it is directly combined with the simulated concentration sequence of the precursor here to perform photochemical reaction simulation calculations to obtain the corresponding first ozone characteristic parameter. In other embodiments, the nitrous acid concentration sequence can also be numerically adjusted according to the preset concentration adjustment value to obtain the nitrous acid concentration sequence corresponding to each simulated concentration sequence of the precursor, and then execute S122 correspondingly.

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

[0121] After completing the execution of S120, S130 can be executed subsequently.

[0122] S130: Perform derivative calculations based on each first ozone characteristic parameter and the preset concentration adjustment value to determine each order of partial derivatives.

[0123] In the embodiments of the present disclosure, each order of partial derivatives is the partial derivative of the first ozone characteristic parameter with respect to the concentration of volatile organic compounds and / or volatile nitrogen oxides.

[0124] 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 specific implementations, the difference between the first ozone characteristic parameters corresponding to the adjacent concentrations of volatile organic compounds and / or volatile nitrogen oxides can be calculated, and then the first-order partial derivative and the second-order partial derivative are calculated using the aforementioned difference.

[0125] According to the aforementioned analysis, taking the first ozone characteristic parameter in the middle of the assumed square as the reference, the differences calculated based on the first ozone characteristic parameters on both sides are respectively and , then it can be determined that the first-order partial derivative adopts , and the second-order partial derivative is . The first-order partial derivative and the second-order partial derivative for the change in the concentration of volatile organic compounds and the change in the concentration of nitrogen oxides can be calculated according to the aforementioned method.

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

[0127] S140: Based on the partial derivatives of each order, the preset concentration adjustment values, and the precursor concentrations at the target time corresponding to each first ozone characteristic parameter, calculate the second ozone characteristic parameter for the precursor concentration near other precursors by summing according to the Taylor expansion and ignoring the high-order expansion terms, and repeat the calculation based on the precursor concentration near other precursors and the corresponding second ozone characteristic parameter, partial derivative, and preset concentration adjustment parameter until the second ozone characteristic parameters corresponding to all target precursor concentrations are calculated.

[0128] According to the principle of calculus, based on the Taylor expansion formula, near the second ozone characteristic parameter at the position can be obtained through the Taylor expansion formula, and the corresponding terms of the high-order derivatives in the Taylor expansion formula are relative to the preceding constant terms and low-order derivative terms. Therefore, the high-order derivative terms can be ignored and will not have a great impact on the result. For example, using its difference from the true is relatively small.

[0129] In addition, through the analysis of a large number of ozone non-linear contour plots, it is found that except for the regions with relatively low precursor concentrations, the shapes of the ozone non-linear contour curves at other regions are basically similar, and the numerical differences of the low-order derivatives of each order are not significant.

[0130] Based on the foregoing analysis, the second ozone characteristic parameter of the adjacent region corresponding to the region of the first ozone characteristic parameter can be calculated by means of Taylor expansion and ignoring the higher-order infinitesimals. For example, if the foregoing second-order derivative is used, the second ozone characteristic parameter at the adjacent position can be calculated using and obtained.

[0131] Similarly, according to the foregoing principle, after obtaining the second ozone characteristic parameter at the adjacent position, the second ozone characteristic parameter at the position adjacent to the adjacent position can be obtained, and by executing the loop, the second ozone characteristic parameters corresponding to all target precursor concentrations can be obtained.

[0132] S150: Based on the first ozone characteristic parameter and the precursor concentration at the corresponding target time, and the second ozone characteristic parameter and the corresponding target precursor concentration, construct an ozone non-linear contour plot.

[0133] After obtaining the foregoing first ozone characteristic parameters and second ozone characteristic parameters, and the precursor concentrations at the target time corresponding to the foregoing ozone characteristic parameters, subsequently, the foregoing precursor concentrations can be standardized (specifically, converting the concentrations of various precursors into a certain specific type of precursor or a specific functional group), and then the ozone non-linear contour curve can be constructed. In the embodiments of the present disclosure, the two coordinate axes of the ozone non-linear contour curve are respectively the normalized values of the volatile organic compound concentration and the normalized value of the volatile nitrogen oxide concentration.

[0134] As analyzed above, by analyzing a large number of ozone non-linear isoconcentration curves obtained by the brute-force method, it is determined that the curve characteristics at various positions of the ozone non-linear isoconcentration curve have the same gradient change law. At the same time, the goal of plotting the ozone non-linear isoconcentration curve is not to obtain absolutely accurate results, but rather to obtain a curve trend chart that can guide practical applications. Based on the foregoing premise, in the solution of the embodiment of the present disclosure, it is considered that after obtaining the first ozone characteristic parameters at multiple concentration points through photochemical simulation with only a relatively small number of precursor sequences, after calculating the derivatives of each order using the foregoing first ozone characteristic parameters, the foregoing derivatives of each order can be used to obtain the second ozone characteristic parameters corresponding to other precursor concentrations according to the Taylor expansion method. Subsequently, the ozone non-linear isoconcentration curve can be constructed based on the first ozone characteristic parameters and the second ozone characteristic parameters. Based on the method of this solution, the ozone non-linear isoconcentration curve can be obtained only by calculating with relatively few data, thus reducing the calculation cost.

[0135] As mentioned above, it is necessary to obtain multiple first ozone characteristic parameters corresponding to the target time. In specific implementation, the target time is a time determined by artificial selection or a time determined according to a pre-set rule. In some embodiments, the computing device may determine the target time by the following S124-S125.

[0136] S124: Based on the combination of the monitored concentration sequence, the meteorological data sequence, and the light flux sequence, perform photochemical reaction simulation calculations to obtain the candidate ozone characteristic parameters at each moment of the monitored concentration sequence.

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

[0138] S125: Select the largest candidate ozone characteristic parameter as the first ozone characteristic parameter corresponding to the said monitored concentration sequence, and use the moment corresponding to the largest candidate ozone characteristic parameter as the target time.

[0139] After obtaining the candidate ozone characteristic parameters at each moment, then compare the sizes of the candidate ozone characteristic parameters at each moment, determine the maximum value among them, and use the maximum value as the first ozone characteristic parameter. Correspondingly, the moment corresponding to the first ozone characteristic parameter is used as the target time.

[0140] By analyzing the normal line of the ozone non-linear characteristic curve determined by the brute-force method, in the region where the precursor concentration is relatively high (correspondingly, the region where the ozone characteristic parameter is relatively large), the curvature of the ozone non-linear characteristic curve is greater, and the gradient certification is more obvious. Correspondingly, by adopting the aforementioned S124-S125, the first ozone characteristic parameter can be positioned in the region where the aforementioned gradient characteristic is more obvious, making the values of the subsequent calculated derivatives more typical and representative.

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

[0142] S126: Based on each precursor simulated concentration sequence respectively, combine it with the meteorological data sequence and the light flux sequence to perform photochemical reaction simulation calculations to obtain other first ozone characteristic parameters corresponding to the target moment.

[0143] The execution method of S126 is the same as the steps explained above. However, when performing the photochemical reaction simulation, the simulation calculation is stopped when the target moment is reached, and only the first ozone characteristic parameters until the target moment are obtained.

[0144] It should be noted here that the implementation processes of S124-S126 do not conflict with S01, S02, S03, S121, S122, and S123 in the previous text. In specific implementation, they can be logically combined in the execution process of each detailed step to form a new scheme integrating each execution step. That is to say, the formaldehyde concentration sequence and the nitrite concentration sequence can be determined first, and then each precursor concentration sequence and the aforementioned formaldehyde concentration sequence are converted into a standard mechanism substance concentration sequence, and the photochemical reaction simulation calculation is performed using the standard mechanism substance concentration sequence to obtain the corresponding first ozone characteristic parameters.

[0145] When performing the photochemical reaction simulation, in addition to the aforementioned precursor concentration sequence, meteorological data sequence, and light flux sequence, in actual situations, due to reasons such as meteorological conditions, the height of the pollutant gas boundary layer will also change, and the change 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.

[0146] In some embodiments of the present disclosure, before performing the foregoing S120, the computing device further 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 ground surface. On the premise of determining the height adjustment characteristics of the mixing layer, when the computing device performs the foregoing S120, specifically: using the mixing layer height parameter as a boundary layer constraint, based on the monitored concentration sequence or a precursor simulated concentration sequence, combined with the meteorological data sequence and the light flux sequence, performing a photochemical reaction simulation calculation to obtain a first ozone characteristic parameter.

[0147] In a specific implementation, in order to verify the rationality of the foregoing method, in a specific application, the computing device also uses the brute force method to draw an ozone non-linear isoconcentration curve. Specifically, based on the maximum allowable limits of volatile organic compound concentration and nitrogen oxide concentration determined by the existing environmental control and the set percentage change, and the various precursor types obtained in the foregoing scheme, various measurement bases are obtained by adjusting according to the set percentage, 441 data points are obtained, and photochemical simulations are performed on the 441 data points according to the foregoing environmental conditions to obtain the first ozone characteristic parameter of each point, and the ozone non-linear isoconcentration curve is directly drawn based on all the first ozone characteristic parameters. Through the foregoing method, the constructed ozone non-linear isoconcentration curve is used to verify the rationality of the change characteristics of each position of the ozone non-linear isoconcentration curve constructed in S110 - S150 above.

[0148] Figure 2 is the non-linear isoconcentration curve graph of ozone concentration constructed by using the solution of the embodiment of the present disclosure. Figure 3 is the non-linear isoconcentration curve graph of ozone concentration constructed by using the brute force method. Figure 4 is the non-linear isoconcentration curve graph of ozone formation potential constructed by using the solution of the embodiment of the present disclosure. Figure 5 is the non-linear isoconcentration curve graph of ozone formation potential constructed by using the brute force method. Comparison Figure 2 and Figure 3 and Figure 4 and Figure 5 It can be seen that the ozone non-linear isoconcentration curve constructed by using this solution has the same trend as the corresponding area in the curve graph drawn by the brute force method in the case of higher precursor concentration, and the same trend as the corresponding area in the curve graph drawn by the brute force method in the case of lower precursor concentration (of course, this has something to do with the selection of the precursor concentration used to calculate the first ozone characteristic curve. It should be noted here that because the new precursor concentration points are determined by adjusting the precursor concentration percentage by 5% in the calculation, Figure 2 and Figure 4The axis data in [it] are not the same, but the precursor concentrations reflected by the binary essence are corresponding). Considering that the ozone generation rate is very low when the precursor concentration is low and meets the conventional environmental protection requirements, this is not the focus of actual environmental monitoring and pollution source control, and the aforementioned differences can be ignored in practical applications.

[0149] Figure 6 It is a schematic structural diagram of the architecture based on which the present disclosure implements the foregoing method in a specific embodiment. As Figure 6 shown, the present disclosure first needs to obtain environmental observation data, meteorological observation data (including the mixing layer height coefficient, that is, boundary layer data), photolysis parameters, and geographic information data. Subsequently, after data collation, model calculation methods are used to calculate the first ozone characteristic parameter and the second ozone characteristic parameter. Subsequently, an ozone non-linear isoconcentration curve graph is constructed according to the method of this embodiment and the strong force method respectively, and a sensitivity analysis is performed to determine the feasibility of implementing the method of the solution of this embodiment.

[0150] In addition to providing the foregoing method for constructing an ozone non-linear isoconcentration curve graph, the embodiment of the present disclosure also provides an ozone non-linear isoconcentration curve construction device. Figure 7 It is a schematic diagram of the ozone non-linear isoconcentration curve construction device provided by the embodiment of the present disclosure. As Figure 7 shown, the ozone non-linear isoconcentration 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.

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

[0152] The photochemical reaction model calculation unit 702 is used to perform photochemical reaction simulation calculations by combining the monitoring concentration sequence or a group of precursor simulation concentration sequences with the meteorological data sequence and the light flux sequence respectively, to obtain multiple first ozone characteristic parameters corresponding to the target time; the first ozone characteristic parameter is the ozone generation potential or the ozone concentration.

[0153] The derivative calculation unit 703 is used to perform derivative calculations based on each first ozone characteristic parameter and a preset concentration adjustment value to determine each order partial derivative.

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

[0155] The curve graph construction unit 705 is configured to construct a non-linear ozone isoconcentration graph based on the first ozone characteristic parameter and the precursor concentration at the corresponding target moment, and the second ozone characteristic parameter and the corresponding target precursor concentration.

[0156] In some embodiments, the photochemical reaction model calculation unit 702 first performs photochemical reaction simulation calculations based on the combination of the monitored concentration sequence, the meteorological data sequence, and the light flux sequence to obtain the 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 uses the moment corresponding to the largest candidate ozone characteristic parameter as the target moment; and then performs photochemical reaction simulation calculations respectively based on each precursor simulated concentration sequence, in combination with the meteorological data sequence and the light flux sequence, to obtain the other first ozone characteristic parameters corresponding to the target moment.

[0157] In some embodiments, the basic data determination unit 701 obtains the ozone concentration sequence monitored in the same period, and 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 a formaldehyde concentration sequence based on the formaldehyde concentration. The photochemical reaction model calculation unit 702 performs photochemical reaction simulation calculations based on the combination of the monitored concentration sequence or a precursor simulated concentration sequence, the meteorological data sequence, the light flux sequence, and the formaldehyde concentration sequence.

[0158] In some embodiments, the precursor monitored 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. Correspondingly, the photochemical reaction model calculation unit 702 performs photochemical reaction simulation calculations based on the combination of the monitored concentration sequence or a precursor simulated concentration sequence, the meteorological data sequence, the light flux sequence, and the nitrous acid concentration sequence.

[0159] In some embodiments, the photochemical reaction model calculation unit 702 first performs a standard mechanism substance concentration conversion on the concentration sequences of various volatile organic compounds in the monitored concentration sequence and the precursor simulated concentration sequence to determine the standard mechanism substance concentration sequences corresponding to the monitored concentration sequence and the precursor simulated concentration sequence. Subsequently, based on the standard mechanism substance concentration sequence and the volatile organic compound concentration sequence corresponding to 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.

[0160] In some embodiments, the basic data determination unit 701 obtains the longitude and latitude information of the monitoring location, and performs a light projection orthographic calculation based on the longitude and latitude information to determine the orthographic projection conversion coefficient. Subsequently, based on the solar incident intensity and the orthographic projection conversion coefficient, a light flux sequence is obtained.

[0161] 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 a photochemical reaction simulation calculation based on the monitored concentration sequence or a precursor simulated concentration sequence in combination with the meteorological data sequence and the light flux sequence.

[0162] The embodiments of the present disclosure also provide a computing device for implementing the foregoing method. Figure 8 It is a schematic structural diagram of the computing device provided by the embodiments of the present disclosure. Specifically refer to the following Figure 8 which shows a schematic structural diagram of a computing device 800 suitable for implementing the present disclosure. Figure 8 The shown computing device is only an example and should not impose any limitation on the functions and usage scope of the embodiments of the present disclosure.

[0163] As Figure 8 shown, the computing device 800 may include a processing device (such as a central processing unit, a graphics processing unit, etc.) 801, which can perform various appropriate actions and processes according to the program stored in the read-only memory ROM 802 or the program loaded from the storage device 808 into the 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 through a bus 804. The input / output I / O interface 805 is also connected to the bus 804.

[0164] Typically, the following devices can be connected to the I / O interface 805: input devices 805 including, for example, a touch screen, a touchpad, 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, magnetic tapes, hard disks, etc.; and communication devices 809. The communication device 809 can allow the computing device 800 to communicate with other devices wirelessly or wiredly to exchange data. Although Figure 8 the computing device 800 with various devices is shown, it should be understood that it is not required to implement or have all the shown devices. Instead, more or fewer devices can be implemented or had.

[0165] Specifically, according to an embodiment of the present disclosure, the processes described above with reference to the flowcharts can be implemented as computer software programs. For example, an embodiment of the present disclosure includes a computer program product that includes a computer program carried on a non-transitory computer-readable medium, and the computer program includes program codes for performing the methods shown in the flowcharts. In such an embodiment, the computer program can be downloaded and installed from a network through the communication device 809, or installed from the storage device 808, or installed from the ROM 802. When the computer program is executed by the processing device 801, the above-mentioned functions defined in the methods of the embodiments of the present disclosure are executed.

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

[0167] The computer-readable storage medium can be, for example, but not limited to, an electrical, magnetic, optical, electromagnetic, infrared, or semiconductor system, apparatus, or device, or any combination of the above. More specific examples of the computer-readable storage medium can 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 of the above. In the present disclosure, the computer-readable storage medium can be any tangible medium that contains or stores a program that can be used by or in combination with an instruction execution system, apparatus, or device.

[0168] A computer-readable signal medium may include a data signal propagated in a baseband or as part of a carrier wave, in which computer-readable program code is carried. Such a propagated data signal may take many forms, including but not limited to electromagnetic signals, optical signals, or any suitable combination of the foregoing. The computer-readable signal medium may also be any computer-readable medium other than a computer-readable storage medium, which can send, propagate, or transmit a program for use by or in connection with an instruction execution system, apparatus, or device. The program code contained on the computer-readable medium may be transmitted using any appropriate medium, including but not limited to wires, optical cables, RF (radio frequency), etc., or any suitable combination of the foregoing.

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

[0170] The above-mentioned computer-readable medium may be included in the above-mentioned computing device; or may exist separately without being assembled into the computing device.

[0171] Computer program code for performing the operations of the present disclosure may be written in one or more programming languages or combinations thereof, including but not limited to object-oriented programming languages such as Java, Smalltalk, C++, and also including conventional procedural programming languages such as the "C" language or similar programming languages. The program code may be executed entirely on the tester's computer, partially on the tester's computer, executed as a stand-alone software package, partially on the tester's computer and partially on a remote computer, or entirely on a remote computer or server. In the case of 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., by using an Internet service provider to connect through the Internet).

[0172] The flowcharts and block diagrams in the accompanying drawings illustrate the possible architectures, functions, and operations of systems, methods, and computer program products according to various embodiments of the present disclosure. In this regard, each block in the flowchart or block diagram may represent a module, a segment of a program, or a portion of code that contains one or more executable instructions for implementing the specified logical function. It should also be noted that, in some alternative implementations, the functions noted in the blocks may occur in an order different from that noted in the accompanying drawings. For example, two consecutive blocks shown may actually be executed substantially in parallel, or they may sometimes be executed in the reverse order, depending on the functions involved. It should also be noted that each block in the block diagram and / or flowchart, and combinations of blocks in the block diagram and / or flowchart, may be implemented by a dedicated hardware-based system that performs the specified functions or operations, or may be implemented by a combination of dedicated hardware and computer instructions.

[0173] The units described in the embodiments of the present disclosure can be implemented in software or in hardware. Among them, the name of the unit does not constitute a limitation on the unit itself in some cases. The functions described above herein can 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 arrays (FPGAs), application specific integrated circuits (ASICs), application specific standard products (ASSPs), systems on a chip (SOCs), complex programmable logic devices (CPLDs), and so on.

[0174] The above are only specific implementation manners of the present disclosure, enabling those skilled in the art to understand or implement the present disclosure. Various modifications to these embodiments will be obvious 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 rather to the widest scope consistent with the principles and novel features disclosed herein.

Claims

1. A method for constructing an ozone non-linear equivalent curve graph, characterized in that Including: Obtaining a precursor monitoring concentration sequence, a meteorological data sequence, and a light flux sequence monitored in the same period, and obtaining multiple sets of adjacent precursor simulation concentration sequences based on the monitoring concentration sequence and a preset concentration adjustment value; both the monitoring concentration sequence and the simulation concentration sequence include a volatile organic compound concentration sequence and a nitrogen oxide concentration sequence; Respectively performing photochemical reaction simulation calculations by combining the monitoring concentration sequence or a set of the precursor simulation concentration sequences with the meteorological data sequence and the light flux sequence to obtain multiple first ozone characteristic parameters corresponding to a target time; the first ozone characteristic parameters are ozone formation potential or ozone concentration; Performing derivative calculations based on each first ozone characteristic parameter and the preset concentration adjustment value to determine each order partial derivative; Calculating a second ozone characteristic parameter of adjacent other precursor concentrations in a manner of summing according to Taylor expansion and ignoring high-order expansion terms based on the each order partial derivative, the preset concentration adjustment value, and the precursor concentrations corresponding to each of the first ozone characteristic parameters at the target time, and repeating the calculation based on the adjacent other precursor concentrations and the corresponding second ozone characteristic parameters, the partial derivative, and the preset concentration adjustment parameter until calculating the second ozone characteristic parameters corresponding to all target precursor concentrations; Constructing an ozone non-linear contour graph based on the first ozone characteristic parameters and the precursor concentrations at the corresponding target time, and the second ozone characteristic parameters and the corresponding target precursor concentrations.

2. The method according to claim 1, characterized in that, The respectively performing photochemical reaction simulation calculations by combining the monitoring concentration sequence or a precursor simulation concentration sequence with the meteorological data sequence and the light flux sequence to obtain multiple first ozone characteristic parameters corresponding to a target time includes: Performing photochemical reaction simulation calculations by combining the monitoring concentration sequence with the meteorological data sequence and the light flux sequence to obtain 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 moment corresponding to the largest candidate ozone characteristic parameter as the target time; Respectively performing photochemical reaction simulation calculations by combining each of the precursor simulation concentration sequences 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, characterized in that, The precursor monitoring concentration sequence includes an acetylene concentration sequence; the method further includes: obtaining an ozone concentration sequence monitored in the same period; Before respectively performing photochemical reaction simulation calculations by combining the monitoring concentration sequence or a precursor simulation concentration sequence with the meteorological data sequence and the light flux sequence, the method further includes: 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; Performing photochemical reaction simulation calculations by combining the monitored concentration sequence or one of the precursor simulated concentration sequences with the meteorological data sequence and the light flux sequence includes: performing photochemical reaction simulation calculations by combining the monitored concentration sequence or one of the precursor simulated concentration sequences with the meteorological data sequence, the light flux sequence, and the formaldehyde concentration sequence.

4. The method according to claim 1 or 2, characterized in that, The precursor monitored concentration sequence includes a nitrogen dioxide concentration sequence; Before performing photochemical reaction simulation calculations by combining the monitored concentration sequence or one of the precursor simulated concentration sequences with the meteorological data sequence and the light flux sequence respectively, the method further includes: 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; Performing photochemical reaction simulation calculations by combining the monitored concentration sequence or one of the precursor simulated concentration sequences with the meteorological data sequence and the light flux sequence includes: performing photochemical reaction simulation calculations by combining the monitored concentration sequence or one of the precursor simulated concentration sequences with the meteorological data sequence, the light flux sequence, and the nitrous acid concentration sequence.

5. The method according to claim 1 or 2, characterized in that, Before performing photochemical reaction simulation calculations by combining the monitored concentration sequence or one of the precursor simulated concentration sequences with the meteorological data sequence and the light flux sequence respectively, it further includes: Performing standard mechanism substance concentration conversion on various volatile organic compound concentration sequences in the monitored concentration sequence and the precursor simulated concentration sequence to determine the standard mechanism substance concentration sequences corresponding to the monitored concentration sequence and the precursor simulated concentration sequence; Performing photochemical reaction simulation calculations by combining the monitored concentration sequence or one of the precursor simulated concentration sequences with the meteorological data sequence and the light flux sequence respectively includes: performing photochemical reaction simulation calculations by combining the standard mechanism substance concentration sequences and volatile organic compound concentration sequences corresponding to the monitored concentration sequence or one of the precursor simulated concentration sequences with the meteorological data sequence and the light flux sequence.

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

7. The method according to claim 1 or 2, characterized in that, Before performing photochemical reaction simulation calculations by combining the monitored concentration sequence or one of the precursor simulated concentration sequences with the meteorological data sequence and the light flux sequence respectively, the method further includes: determining the mixed layer height parameter; Performing photochemical reaction simulation calculations by combining the monitored concentration sequence or one of the precursor simulated concentration sequences with the meteorological data sequence and the light flux sequence respectively includes: Using the mixed layer height parameter as a boundary layer constraint, based on the monitored concentration sequence or one of the precursor simulated concentration sequences, perform photochemical reaction simulation calculations in combination with the meteorological data sequence and the light flux sequence.

8. An ozone non-linear equivalent curve graph construction device, characterized in that, Including: A basic data determination unit, configured to obtain the precursor monitored concentration sequence, the meteorological data sequence, and the light flux sequence monitored at the same time period, and obtain multiple sets of adjacent precursor simulated concentration sequences based on the monitored concentration sequence and a preset concentration adjustment value; both the monitored concentration sequence and the simulated concentration sequence 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 respectively based on the monitored concentration sequence or one of the precursor simulated concentration sequences, in combination with the meteorological data sequence and the light flux sequence, to obtain multiple first ozone characteristic parameters corresponding to the target time; the first ozone characteristic parameter is the ozone formation potential or the ozone concentration; A derivative calculation unit, configured to perform derivative calculations based on each first ozone characteristic parameter and the preset concentration adjustment value to determine the partial derivatives of each order; A point value data calculation unit, configured to calculate the second ozone characteristic parameter of the adjacent other precursor concentrations by summing according to Taylor expansion and ignoring the high-order expansion terms based on the partial derivatives of each order, the preset concentration adjustment value, and the precursor concentrations corresponding to each first ozone characteristic parameter at the target time, and repeat the calculation based on the adjacent other precursor concentrations and the corresponding second ozone characteristic parameters, the partial derivatives, and the preset concentration adjustment parameters until the second ozone characteristic parameters corresponding to all target precursor concentrations are calculated; A curve graph construction unit, configured to construct an ozone non-linear contour graph based on the first ozone characteristic parameter and the precursor concentration at the corresponding target time, and the second ozone characteristic parameter and the corresponding target precursor concentration.

9. A computing device, characterized in that, Including a processor and a memory, the memory is used to store a computer program; when the computer program is loaded by the processor, the processor executes the ozone non-linear contour graph construction method according to any one of claims 1-7.

10. A computer-readable storage medium, characterized in that, The storage medium stores a computer program, and when the computer program is executed by a processor, the processor realizes the ozone non-linear contour graph construction method according to any one of claims 1-7.

Images