Sea-land breeze ozone transport channel analysis method and system

Abstract

The present application relates to the technical field of meteorological analysis, and discloses a sea-land wind ozone transmission channel analysis method and system, the method comprising: acquiring meteorological elements and ozone concentration of a target detection area; identifying the duration of sea wind in the meteorological elements, counting the distance from the coast of the target detection area, taking the decay relationship between the duration of sea wind and the distance from the coast as the sea wind penetration intensity value of the target detection area; extracting the ozone change slope of the ozone concentration from the accumulation period to the dissipation period, dividing the target detection area into a coastal strong penetration area and an inland weak penetration area according to the sea wind penetration intensity value; determining the inversion layer height decay rate of the target detection area based on the temperature vertical gradient in the meteorological elements; superimposing the inversion layer height decay rate and the ozone change slope as the ozone subsidence contribution coefficient of the target detection area, and the present application can solve the problem that the prior art cannot support refined ozone source tracing.

Description

TECHNICAL FIELD

[0001] The present application relates to the field of meteorological analysis, and in particular to a sea-land wind ozone transport channel analysis method and system. BACKGROUND

[0002] The existing ozone transport analysis method under the influence of sea-land wind mostly relies on a single meteorological element or simple spatial distance for rough evaluation, and it is difficult to accurately depict the penetration strength of sea wind to inland areas. The traditional method often ignores the dynamic decay relationship between the duration of sea wind and the distance from the shore, and does not consider the influence of factors such as terrain undulation and atmospheric stability on the penetration ability, resulting in blurred division of the ozone transport boundary between coastal and inland areas, and inability to accurately identify strong / weak penetration areas, making it difficult to support refined ozone source tracing.

[0003] At the same time, the existing technology lacks coupling research on the decay law of inversion layer height and the change trend of ozone concentration when analyzing the vertical transport of ozone, and only analyzes the vertical gradient of temperature or the change of ozone concentration in isolation, without quantifying the contribution coefficient of the inversion layer to the sinking of ozone. This makes the ozone transport channel unable to reflect the internal mechanism of ozone accumulation and dissipation in different regions, resulting in insufficient prediction accuracy of the migration path of ozone pollution, and making it difficult to meet the decision-making needs of regional ozone pollution joint prevention and control. SUMMARY

[0004] The present application provides a sea-land wind ozone transport channel analysis method and system to solve the problems raised in the background art.

[0005] To achieve the above-mentioned purpose, the present application provides a sea-land wind ozone transport channel analysis method, which comprises:

[0006] S1: Obtain meteorological elements and ozone concentration of a target detection area;

[0007] S2: Identify the duration of sea wind in the meteorological elements, and count the distance from the shore of the target detection area, and take the decay relationship between the duration of sea wind and the distance from the shore as the sea wind penetration strength value of the target detection area;

[0008] S3: Extract the ozone change slope of the ozone concentration from the accumulation period to the dissipation period, and divide the target detection area into coastal strong penetration area and inland weak penetration area according to the sea wind penetration strength value;

[0009] S4: Determine the inversion layer height decay rate of the target detection area based on the vertical gradient of temperature in the meteorological elements;

[0010] S5: Superimpose the inversion layer height decay rate and the ozone change slope as the ozone sinking contribution coefficient of the target detection area;

[0011] S6: superimposing the coastal strong-permeation area and the inland weak-permeation area with the contribution coefficient to form a hierarchical ozone transmission channel.

[0012] Preferably, the meteorological elements and ozone concentration of the target detection area are obtained by:

[0013] Continuously collecting meteorological observation raw data of the target detection area and synchronously obtaining regional ozone information of the target detection area;

[0014] Eliminating abnormal data in the meteorological observation raw data and the regional ozone information to obtain the meteorological elements and ozone concentration of the target detection area.

[0015] Preferably, the duration of the sea wind in the meteorological elements is identified, and the off-shore distance of the target detection area is counted, including:

[0016] Identifying a wind direction change event of a wind field time domain sequence in the meteorological elements, and marking a duration of the wind direction change event as the duration of the sea wind;

[0017] Determining the off-shore distance of the target detection area based on a coordinate reference position of the target detection area.

[0018] Preferably, the attenuation relationship between the duration of the sea wind and the off-shore distance is taken as the sea wind permeation intensity value of the target detection area, including:

[0019] Time-weighted evaluation of the duration of the sea wind is performed based on a preset seasonal wind field scale parameter to obtain a time-weighted intensity of the duration of the sea wind;

[0020] Analyzing a spatial attenuation intensity value of the off-shore distance;

[0021] Synthesizing the time-weighted intensity and the spatial attenuation intensity value in terms of permeation capacity to obtain the sea wind permeation intensity value of the target detection area.

[0022] Preferably, the time-weighted intensity and the spatial attenuation intensity value are synthesized in terms of permeation capacity to obtain the sea wind permeation intensity value of the target detection area, including:

[0023]

[0024] Wherein: is the sea wind permeation intensity value, is the time-weighted intensity, is a terrain permeation transfer function, is an atmospheric stability factor, is the spatial attenuation intensity value, is an absolute value of an altitude difference, Maximum relief.

[0025] Preferably, the target detection area is divided into a coastal strong penetration area and an inland weak penetration area according to the sea wind penetration intensity value, comprising:

[0026] The sea wind penetration intensity value is taken as a penetration threshold value of the target detection area;

[0027] The geographical spatial distribution of the target detection area is marked by zones based on the penetration threshold value;

[0028] The spatial continuity of the results of the zone marking is checked, and the results that pass the check are divided into a coastal strong penetration area and an inland weak penetration area.

[0029] Preferably, the inversion layer height decay rate of the target detection area is determined based on the temperature vertical gradient in the meteorological element, comprising:

[0030] The height boundary feature of the temperature vertical gradient is extracted, and the height change difference of the height boundary feature in a continuous time sequence is analyzed synchronously;

[0031] The height change difference obtained by analysis is aggregated as the inversion layer height decay rate of the target detection area.

[0032] Preferably, the inversion layer height decay rate and the ozone change slope are superimposed as the ozone sinking contribution coefficient of the target detection area, comprising:

[0033] The ozone vertical distribution relationship between the change characteristics of the inversion layer height decay rate and the distribution characteristics of the ozone change slope is quantitatively analyzed;

[0034] The ozone vertical distribution relationship is subjected to contribution weight generation to obtain an ozone vertical sinking value;

[0035] The ozone vertical sinking value is matched with the corresponding adaptive area to obtain the ozone sinking contribution coefficient of the target detection area.

[0036] Preferably, the coastal strong penetration area and the inland weak penetration area are spatially superimposed with the contribution coefficient to form a hierarchical ozone transport channel, comprising:

[0037] A composite penetration boundary map is constructed based on the spatial position data of the coastal strong penetration area and the inland weak penetration area;

[0038] The ozone sinking contribution coefficient and the composite penetration boundary map are subjected to contribution coupling to obtain an ozone distribution area;

[0039] dividing the ozone distribution area according to intensity to obtain intensity levels of the ozone distribution area;

[0040] dividing the ozone distribution area into different hierarchical ozone transmission channels based on the intensity levels.

[0041] A sea-land wind ozone transmission channel analysis system, the system comprising:

[0042] a data acquisition module for acquiring meteorological elements and ozone concentration of a target detection area;

[0043] a sea wind influence module for identifying a sea wind duration in the meteorological elements, counting an off-shore distance of the target detection area, and taking an attenuation relationship between the sea wind duration and the off-shore distance as a sea wind penetration intensity value of the target detection area;

[0044] a region division module for extracting an ozone change slope of the ozone concentration from an accumulation period to a dissipation period, and dividing the target detection area into a coastal strong penetration area and an inland weak penetration area according to the sea wind penetration intensity value;

[0045] a temperature analysis module for determining an inversion layer height decay rate of the target detection area based on a temperature vertical gradient in the meteorological elements;

[0046] an ozone contribution module for superimposing the inversion layer height decay rate and the ozone change slope as an ozone sinking contribution coefficient of the target detection area;

[0047] a channel generation module for spatially superimposing the coastal strong penetration area and the inland weak penetration area with the contribution coefficient to form hierarchical ozone transmission channels.

[0048] Advantages

[0049] 1. By converting the attenuation relationship between the sea wind duration and the off-shore distance into the sea wind penetration intensity value, and introducing parameters such as the terrain penetration transfer function and the atmospheric stability factor to synthesize the penetration ability, the limitation of relying on a single element in the traditional method is broken through, the coastal strong penetration area and the inland weak penetration area can be accurately divided, the ozone transmission boundary of different areas is clearly defined, and a more effective basis for identifying the ozone cross-regional transmission path is provided.

[0050] 2. By superimposing the inversion layer height decay rate and the ozone change slope as the ozone sinking contribution coefficient, and forming a hierarchical ozone transport channel through spatial superposition, the coupling analysis of ozone horizontal transport and vertical deposition process is realized, which makes up for the defects of the prior art in isolated analysis of meteorological elements or ozone concentration, not only quantifies the influence weight of the inversion layer on ozone sinking, but also distinguishes the transport channel according to the intensity level, significantly improves the analysis accuracy of the migration rule of ozone pollution, and provides more accurate technical support for regional ozone pollution joint prevention and control. BRIEF DESCRIPTION OF DRAWINGS

[0051] Figure 1 A flowchart of a sea-land wind ozone transport channel analysis method provided by an embodiment of the present application is shown in the figure.

[0052] Figure 2 A functional module diagram of a sea-land wind ozone transport channel analysis system provided by an embodiment of the present application is shown in the figure. DETAILED DESCRIPTION

[0053] It should be understood that the specific embodiments described herein are merely intended to explain the present application, and are not intended to limit the present application.

[0054] The embodiments of the present application provide a sea-land wind ozone transport channel analysis method and system. The execution subject of the sea-land wind ozone transport channel analysis method and system includes but is not limited to at least one of electronic devices such as a server, a terminal, etc., which can be configured to execute the method provided by the embodiments of the present application. In other words, the sea-land wind ozone transport channel analysis method and system can be executed by software or hardware installed in a terminal device or a server device. The server includes but is not limited to a single server, a server cluster, a cloud server, or a cloud server cluster, etc. The server can be a stand-alone server, or a cloud server providing cloud services, cloud databases, cloud computing, cloud functions, cloud storage, network services, cloud communication, middleware services, domain name services, security services, content distribution networks, and big data and artificial intelligence platforms, etc. basic cloud computing services.

[0055] Referring to Figure 1 A flowchart of a sea-land wind ozone transport channel analysis method provided by an embodiment of the present application is shown in the figure. In the embodiment, the sea-land wind ozone transport channel analysis method includes:

[0056] S1: Obtain meteorological elements and ozone concentration of a target detection area.

[0057] In the embodiment, the obtaining of the meteorological elements and ozone concentration of the target detection area includes:

[0058] continuously collect meteorological observation raw data of the target detection area, and synchronously acquire regional ozone information of the target detection area;

[0059] remove abnormal data in the meteorological observation raw data and the regional ozone information to obtain meteorological elements and ozone concentration of the target detection area.

[0060] Specifically, the meteorological elements refer to physical quantities representing the atmospheric state of the target detection area, including wind field (wind direction, wind speed), temperature (especially temperature vertical gradient), air pressure, etc. In the sea-land wind environment, the change of the wind field is the key to identifying the sea breeze. During the day, the ocean temperature is lower than the land, and the sea breeze blows from the ocean to the land. At night, the land breeze is formed in the opposite direction. The temperature vertical gradient affects the formation of the inversion layer, which hinders the vertical diffusion of ozone.

[0061] The ozone concentration is the content of ozone in the air of the target detection area. Under the action of sea-land breeze, the ozone in the coastal area may be transported inland along with the sea breeze. At the same time, during the accumulation period, such as during the day when the light is strong, more ozone is generated through photochemical reaction, and the concentration rises. During the dissipation period, such as at night when photochemical reaction is weakened, ozone gradually decomposes, and the concentration decreases.

[0062] In detail, the continuous collection of meteorological observation raw data is achieved by continuously collecting the original meteorological data of the target detection area through meteorological stations, radars and other equipment. For example, wind direction, wind speed, temperature and other data are recorded every few minutes to ensure that the entire process of the sea breeze from start to finish and the continuous change of temperature over time can be captured.

[0063] Synchronously acquiring regional ozone information is achieved by collecting ozone concentration data in the same time period and the same area as the meteorological data collection using ozone monitors and other equipment. The matching of meteorological elements and ozone concentration data in time and space is ensured to facilitate subsequent analysis of the relationship between the two, such as the relationship between sea breeze intensity and ozone transport amount.

[0064] Further, the removal of abnormal data is for meteorological observation raw data. Abnormal data may include sudden changes in wind speed to 0 or far beyond the normal range due to equipment failure, and irregular and sharp changes in wind direction.

[0065] For regional ozone information, abnormal data may be a sudden abnormal increase or decrease in ozone concentration due to improper calibration of the monitoring equipment.

[0066] These abnormal data are identified and removed through data verification algorithms such as threshold method, mean standard deviation method, etc. to ensure that the obtained meteorological elements and ozone concentration data are true and reliable, providing accurate basic data for subsequent analysis.

[0067] S2: identify the sea wind duration in the meteorological element, count the off-shore distance of the target detection area, and take the attenuation relationship between the sea wind duration and the off-shore distance as the sea wind penetration intensity value of the target detection area.

[0068] In the embodiment, the identification of the sea wind duration in the meteorological element and the counting of the off-shore distance of the target detection area include:

[0069] Identify the wind direction transition event of the wind field time sequence in the meteorological element, and mark the duration of the wind direction transition event as the sea wind duration.

[0070] Determine the off-shore distance of the target detection area based on the coordinate reference positioning of the target detection area.

[0071] Specifically, the wind field time sequence is the target detection area wind direction vector data recorded in time sequence, and is the basic data source for identifying the dynamic change of sea wind. In the physical environment, the thermal difference between sea and land drives the periodic transition of wind field with day and night alternation, and the wind field time sequence can capture this rule.

[0072] The wind direction transition event refers to the critical process that the wind direction changes from land wind to sea wind. Driven by the sea-land temperature gradient, the sea is slow to warm up and the air pressure is high in the daytime, and the sea wind begins to dominate. The wind direction transition event is the identification of the sea wind start-continue-end cycle.

[0073] In detail, data collection and preprocessing is to continuously collect wind field data of the target area by means of meteorological stations, wind profile radars and other equipment, form a continuous time sequence, and filter abnormal values such as wind direction jump caused by equipment failure.

[0074] Wind direction transition identification is to use algorithms such as sliding window + threshold judgment to detect the time point when the wind direction changes from land wind characteristics such as wind direction pointing to the ocean direction to sea wind characteristics in the wind field time sequence, and mark the start and end time of the transition.

[0075] Duration marking is to calculate the time span from the start to the end of the wind direction transition event, that is, the sea wind duration. In the physical environment, the sea wind duration is affected by the sea-land temperature difference intensity and terrain (such as coastal mountain blockage). When the sea-land temperature difference is large on a sunny day, the sea wind duration is longer.

[0076] Specifically, the coordinate reference positioning is based on the geographical coordinates of the target detection area, and is related to the coordinate reference of the coastline, such as marking the coastline as a reference line through GIS map. The coastline is the geographical boundary of sea and land wind, and the off-shore distance directly affects the path length and attenuation degree of sea wind transmission ozone.

[0077] In detail, the geographic information matching is to import the latitude and longitude coordinates of the target detection region into a GIS geographic information system, and superimpose the coastline vector data, which can be obtained through an official geographic database.

[0078] The off-shore distance calculation is to calculate the vertical distance (i.e., off-shore distance) between the target detection region and the coastline by using the GIS spatial analysis function. The farther the off-shore distance, the stronger the influence of the land underlying surface, such as the urban heat island and vegetation coverage, on the sea breeze, and the wind field intensity and ozone transport capacity decrease with the off-shore distance.

[0079] In the embodiment, the attenuation relationship between the sea breeze duration and the off-shore distance is taken as the sea breeze penetration intensity value of the target detection region, which includes:

[0080] The sea breeze duration is time-weighted evaluated based on a preset seasonal wind field scale parameter to obtain a time-weighted intensity of the sea breeze duration;

[0081] The spatial attenuation intensity value of the off-shore distance is analyzed;

[0082] The time-weighted intensity and the spatial attenuation intensity value are synthesized in penetration ability to obtain the sea breeze penetration intensity value of the target detection region.

[0083] Specifically, the seasonal wind field scale parameter is the difference in heat between sea and land in different seasons, and the background wind field, such as the intensity of monsoon, needs to be preset. The parameter is obtained by long-term meteorological observation and statistics, such as a higher basic weight coefficient of sea breeze in summer when the sea-land temperature difference is large, which is used to quantify the influence of seasons on the sea breeze duration. For example, the scale parameter is set to 1.2 in summer and 0.8 in spring, which reflects that the sea breeze is more significantly driven by heat in summer.

[0084] The time-weighted intensity is a quantitative value of the sea breeze duration after the season correction, which reflects the effective intensity of the sea breeze duration under the seasonal background.

[0085] In detail, the weighted calculation is performed by using a formula, for example: time-weighted intensity = sea breeze duration × seasonal wind field scale parameter.

[0086] In summer, the land warms up quickly, and the sea-land temperature difference is large, so the sea breeze duration is longer, and the scale parameter further strengthens this seasonal characteristic, so that the evaluation is more consistent with the actual sea-land wind dynamic mechanism. For example, if the sea breeze lasts for 20 hours in summer, and the scale parameter is 1.2, then the weighted value is 24; if the scale parameter is 0.7 in winter, and the sea breeze lasts for 15 hours, then the weighted value is 10.5.

[0087] Specifically, the spatial attenuation intensity value is a quantitative value of the attenuation of the sea breeze intensity with the increase of the off-shore distance, which reflects the weakening degree of the sea breeze in the geographical space.

[0088] Based on physical laws, the deeper the sea breeze penetrates inland, the faster it decays due to the influence of land surface friction and the reduction of sea-land temperature difference caused by land warming. A spatial decay model is established, such as exponential decay: spatial decay intensity value = initial intensity x e^(-k x distance from shore), k is the decay coefficient, which is fitted through historical data.

[0089] The distance from shore of the target area is substituted into the model to calculate the spatial decay intensity value. For example, 10 km from the coast, k = 0.05, initial intensity 100, after decay, 100 x e^(-0.5) ≈ 60.7; 50 km inland, 100 x e^(-2.5) ≈ 8.2, which shows that the farther from the shore, the weaker the sea breeze penetration, and the lower the ozone transport power.

[0090] In detail, the sea breeze penetration intensity value is the final synthesis result, which integrates time and space dimensions and quantifies the ability of sea breeze to carry ozone inland.

[0091] Further, the time-weighted intensity reflects how long the sea breeze can last to drive the transmission, and the spatial decay intensity reflects how far the transmission can go before the power weakens. The synthesized value comprehensively reflects the effective ability of the sea breeze to carry ozone from the coast to the target area under seasonal background. For example, in summer, the time-weighted intensity of a certain area is 24, and the spatial decay intensity is 60.7, the synthesis is 1456.8, the higher the value, the more effective the sea breeze in transporting ozone, which is the core indicator to judge whether the ozone transport channel is active.

[0092] More specifically, sea-land breeze is a local circulation driven by sea-land thermal difference. Seasonal effects on sea-land temperature difference and distance from shore affect the degree of sea breeze interference from land. Through time weighting + spatial decay + synthesis, the ability of sea breeze to carry ozone inland in real sea-land environment is accurately quantified, providing key parameters for subsequent identification of ozone transport channels and analysis of transport rules.

[0093] In the embodiment, the time-weighted intensity and the spatial decay intensity value are synthesized to obtain the sea breeze penetration intensity value of the target detection area, comprising:

[0094]

[0095] Wherein: is the sea breeze penetration intensity value, is the time-weighted intensity, is the terrain penetration transfer function, is the atmospheric stability factor, is the spatial decay intensity value, is the absolute value of the altitude difference, is the maximum terrain relief.

[0096] Specifically, the sea wind penetration intensity value is a model that quantifies the penetration ability of sea wind carrying ozone to the target area in the three dimensions of time, space, and terrain. The core is to use the influence of various physical environmental factors on sea wind transmission to identify the ozone transmission channel.

[0097] Specifically, the time-weighted intensity reflects the effective dynamic intensity of the duration of sea wind in the time dimension of season, day and night, etc. Influenced by the seasonal changes in the thermal difference between sea and land.

[0098] The time-weighted intensity is obtained by weighting the duration of sea wind based on the seasonal wind field scale parameter, that is, the original duration of sea wind is converted into time-weighted intensity by combining the seasonal correction coefficient, which reflects the difference in the ability of sea wind to drive ozone transmission in different seasons.

[0099] The atmospheric stability factor reflects the influence of the vertical stratification of the atmosphere on sea wind transmission. If the atmosphere is stable, such as the existence of an inversion layer, the sea wind is easily inhibited, the value is large, and when it is unstable, such as strong convection, the value is small.

[0100] The spatial decay intensity value reflects the degree of attenuation of sea wind with increasing distance from the shore. The farther the distance from the shore, the larger the value, the weaker the sea wind power, and the spatial decay step is calculated by analyzing the distance from the shore, which reflects the weakening of the sea wind in geographical space.

[0101] In detail, the exponential function simulates the attenuation law of sea wind power with atmospheric stability + distance from the shore, or the larger the value, the closer to 0, the closer to 1.

[0102] For example: in coastal areas the value is small, and if the atmosphere is unstable, ≈1, the correction term ≈0, but the actual coastal sea wind power is strong, indicating that the should be combined together; in inland areas the value is large, and if the atmosphere is stable, ≈0, the correction term ≈1, which reflects that the sea wind power has decayed to the limit and is close to pure land wind, and the ozone transmission is weakly affected by the sea wind.

[0103] Specifically, the calculation formula of the terrain penetration transfer function is:

[0104]

[0105] Among them: is the terrain penetration transfer function, is the coast curvature coefficient.

[0106] Furthermore, through the composite structure of the hyperbolic cosine function (cosh) and its inverse function (arccosh), and after conformal mapping of complex functions, the original three-dimensional terrain curvature is transformed into a particular solution of the standard hyperbolic differential equation.

[0107] This design gives the formula a dual precise matching capability: when the coastal terrain is flat ( When ≈1), It degenerates into an exponentially decaying function, which is consistent with the energy dissipation characteristics of a gentle uphill slope; when facing a curved coastline ( When >1), The value increases through the multiplication factor. It significantly improves the function decay rate and accurately quantifies the additional kinetic energy loss of the sea breeze streamline caused by the tortuous coastline.

[0108] Especially for the problem of abrupt vortex change on the leeward slope where traditional models completely fail, when the target point is close to the topographic peak, The function value rises sharply to positive infinity at the critical point (x=0.95), leading to... Achieve a precipitous drop.

[0109] This mathematical mutation behavior corresponds precisely to the flow separation phenomenon that occurs in reality when airflow crosses a mountain ridge due to a sudden change in air pressure.

[0110] More specifically, through normalized parameters It eliminates the regional dependence of absolute altitude values ​​and has global applicability.

[0111] Through coefficients The fractal dimension is generated from the coastline. As the coastline's tortuosity increases, such as the fractal dimension of a bedrock coastline changing from 1.6 to 1.8, The value increased from 0.4 to 0.2, which shortened the sea breeze penetration distance under the same terrain conditions, accurately reflecting the obstruction effect of the spiral bay shore on airflow.

[0112] when <0.3 ≈1 (unobstructed airflow penetration), There is a unique inflection point when the value is >0.6 (corresponding to the canyon acceleration effect trigger line). When the value is greater than 0.95, the exponential decay mechanism is automatically activated, which enables the early warning system to respond to mountain vortex disasters in advance.

[0113] S3: Extract the ozone concentration change slope from the accumulation period to the dissipation period, and divide the target detection area into a coastal high-permeability zone and an inland low-permeability zone according to the sea breeze infiltration intensity value.

[0114] In the embodiment, the target detection area is divided into a coastal strong penetration area and an inland weak penetration area according to the sea wind penetration intensity value, including:

[0115] The sea wind penetration intensity value is taken as a penetration threshold value of the target detection area;

[0116] The geographical spatial distribution of the target detection area is marked by zones based on the penetration threshold value;

[0117] The spatial continuity of the marking result is checked, and the result passing the check is divided into a coastal strong penetration area and an inland weak penetration area.

[0118] Specifically, the ozone change slope is used to measure the change rate of the ozone concentration from the accumulation period to the dissipation period, reflecting the dynamic change trend of the ozone concentration over time.

[0119] The generation and decomposition of ozone are affected by light, temperature, volatile organic compounds (VOCs), and nitrogen oxide (NOx) concentrations. During the day, the light is strong, the photochemical reaction is active, the ozone generation rate is greater than the decomposition rate, and it is in the accumulation period. At night, the light weakens, the photochemical reaction slows down, and the ozone gradually decomposes, entering the dissipation period. The larger the ozone change slope, the more dramatic the change in ozone concentration from the accumulation period to the dissipation period.

[0120] In detail, the accumulation period and the dissipation period of the ozone concentration are accurately determined according to the time sequence. Usually, based on the analysis of local meteorological conditions and historical data of ozone concentration, combined with the change law of light intensity to determine.

[0121] For example, after sunrise, with the increase of light, the ozone concentration begins to rise, and it enters the accumulation period; after sunset, the light weakens, and the ozone concentration gradually decreases, which is the dissipation period.

[0122] The change slope of the ozone concentration over time in these two stages is calculated. Linear regression and other algorithms can be used to fit the ozone concentration data at multiple time points in the accumulation period and the dissipation period to obtain the slope of the fitting straight line, which is the ozone change slope.

[0123] Specifically, the sea wind penetration intensity value is a comprehensive consideration of the duration of sea wind, distance from the shore, terrain, and atmospheric stability, and quantifies the ability of sea wind penetration to inland. In the physical environment, the longer the duration of sea wind, the closer the distance from the shore, the flatter the terrain, and the more unstable the atmosphere, the higher the sea wind penetration intensity value. The sea wind penetration intensity value as a penetration threshold value is the key basis for dividing the coastal strong penetration area and the inland weak penetration area.

[0124] In detail, the calculated sea wind penetration intensity value is set as the penetration threshold of the target detection area. This threshold represents the ability limit of the sea wind that can effectively penetrate, and is used to distinguish the degree of influence of different areas by sea wind.

[0125] At the same time, with the help of geographic information system technology, based on the set penetration threshold, the geographic spatial distribution of the target detection area is marked by partition. For each geographic unit, compare its sea wind penetration intensity value with the penetration threshold.

[0126] If the sea wind penetration intensity value of a certain area is greater than or equal to the penetration threshold, it is marked as possibly belonging to the coastal strong penetration area; if it is less than the penetration threshold, it is marked as possibly belonging to the inland weak penetration area. In the marking process, the continuity of geographic space will be considered to avoid isolated marking points.

[0127] Finally, check whether the marked area has broken or isolated parts to ensure that the partition result has continuity and reasonableness in space. If there are discontinuous areas, it may be due to data anomalies or special topography, etc. Further analysis and processing are needed.

[0128] For example, the discontinuous boundary can be adjusted by checking the marking situation of adjacent areas. Finally, the result that passes the verification is formally divided into coastal strong penetration area and inland weak penetration area, which can more accurately reflect the difference in the influence of sea wind on ozone transport in different areas.

[0129] S4: determining the inversion layer height decay rate of the target detection area based on the temperature vertical gradient in the meteorological elements.

[0130] In this embodiment, the determination of the inversion layer height decay rate of the target detection area based on the temperature vertical gradient in the meteorological elements includes:

[0131] extracting the height boundary feature of the temperature vertical gradient, and synchronously analyzing the height change difference of the height boundary feature in the continuous time sequence;

[0132] aggregating the height change difference obtained by analysis as the inversion layer height decay rate of the target detection area.

[0133] Specifically, focusing on the calculation of the inversion layer height decay rate, the inversion layer (the air layer whose temperature increases with height) will inhibit the vertical diffusion of ozone, and the thermal difference of sea-land wind will dynamically change the inversion layer height. By analyzing the temperature vertical gradient→ inversion layer height change→ decay rate, a basis is provided for subsequent judgment of the degree of obstruction of ozone transmission in the vertical direction. The higher the decay rate, the more unstable the inversion layer, and the easier it is for ozone to break through the inversion layer and diffuse.

[0134] In detail, the temperature vertical gradient is the rate of change of temperature in the atmosphere with height. In the sea-land wind environment, the land warms up quickly during the day, and the temperature vertical gradient near the ground is large (unstable stratification is prone to occur); at night, the sea-land temperature difference is small, and the near-surface is prone to form an inversion (temperature vertical gradient ≤ 0).

[0135] The height boundary feature is the top and bottom height of the inversion layer, which is the critical height at which the temperature vertical gradient changes from the inversion feature to the normal stratification, and reflects the distribution range of the inversion layer in the vertical space.

[0136] Further, the vertical sounding data acquisition is to obtain temperature data of the target detection area at continuous height layers by means of a sounding balloon, a weather radar or a vertical observation station.

[0137] Under the driving of the sea-land wind, the vertical distribution of the near-surface temperature changes dynamically with the day and night and the sea-land position, and a high-density vertical sampling is needed to capture the inversion layer.

[0138] S5: superimposing the inversion layer height decay rate and the ozone change slope as an ozone sinking contribution coefficient of the target detection area.

[0139] In this embodiment, the superimposing the inversion layer height decay rate and the ozone change slope as an ozone sinking contribution coefficient of the target detection area comprises:

[0140] quantitative analysis of the ozone vertical distribution relationship of the change characteristics of the inversion layer height decay rate and the distribution characteristics of the ozone change slope;

[0141] generating a contribution weight of the ozone vertical distribution relationship to obtain an ozone vertical sinking value;

[0142] matching an adaptive area corresponding to the ozone vertical sinking value to obtain the ozone sinking contribution coefficient of the target detection area.

[0143] Specifically, through the ozone sinking contribution coefficient calculation, the inversion layer height decay rate reflects the inhibition / promotion dynamics of the ozone vertical diffusion, and the ozone change slope reflects the accumulation / dissipation rate of the ozone concentration. By superimposing the two, the contribution degree of the ozone sinking to the near-surface or the diffusion to the high altitude under the synergistic action of the inversion layer and the ozone concentration change in the sea-land wind environment is quantified, which provides a basis for identifying the vertical sinking effect of the ozone transport channel.

[0144] Specifically, the ozone vertical distribution relationship is the correlation rule of the inversion layer height decay rate and the ozone change slope in the vertical space, such as whether the ozone change slope increases synchronously when the inversion layer rises, and whether the ozone is more prone to sink and accumulate when the inversion layer is pressed down.

[0145] In detail, the target detection area is divided into geographical grids, and in each grid, the inversion layer height decay rate time series and the ozone concentration time series are aligned.

[0146] Sea-land breeze is a local circulation, and the inversion layer and ozone change characteristics of different grids are significantly different, which requires fine gridding.

[0147] For each grid, the vertical direction inversion layer height decay rate is extracted, such as the inversion layer change and ozone concentration data of different height layers, such as near the ground, the top / bottom of the inversion layer.

[0148] Using statistical methods, such as Pearson correlation analysis and multiple linear regression, the correlation between the change of the inversion layer height decay rate and the slope of the ozone concentration change in the vertical profile is analyzed, and a model of the vertical distribution of ozone is established, such as: the slope of the change of ozone=a×inversion layer height decay rate+b, a and b are fitting coefficients.

[0149] Specifically, the contribution weight is the weight coefficient of the influence of the inversion layer height decay rate on the vertical distribution of ozone, reflecting the proportion of the role of the dynamic of the inversion layer in the sinking / diffusion of ozone.

[0150] The stronger the inversion layer, the higher the weight of the sealing effect of ozone sinking; the weaker the inversion layer, the higher the weight of the release effect of ozone diffusion.

[0151] The vertical sinking value of ozone is the quantitative net movement intensity of ozone in the vertical direction, sinking towards the ground or diffusing towards the high altitude, and the sinking value is positive, indicating that ozone converges towards the ground, and negative, indicating that it diffuses towards the high altitude.

[0152] Further, based on the vertical distribution of ozone, through sensitivity analysis, such as changing the inversion layer height decay rate, the response degree of the slope of the change of ozone is observed, and the contribution weight of the inversion layer height decay rate is calibrated.

[0153] By iterating each grid, the weight and parameters are substituted to calculate the vertical sinking value of ozone, reflecting the net movement trend of ozone in the vertical direction in the grid.

[0154] Specifically, the adaptive region is a continuous region in geographical space with similar characteristics of the vertical sinking value of ozone, such as strong sinking area and weak diffusion area, which is affected by sea-land breeze, terrain, and inversion layer, and the characteristics of the ozone transport channel of the adaptive region are consistent.

[0155] The ozone sinking contribution coefficient is the final quantitative result, reflecting the comprehensive contribution degree of ozone sinking towards the ground in the target detection area under the synergistic action of the inversion layer and the change of ozone concentration, and the higher the coefficient, the more significant the effect of ozone sinking.

[0156] In detail, using spatial clustering algorithms such as K-Means clustering and DBSCAN, the vertical sinking value of ozone of all grids is clustered, and adaptive regions such as strong sinking area, weak sinking area, and diffusion area are divided.

[0157] Coastal areas may form a continuous strong sinking area due to the superposition of sea breeze and inversion temperature; the complex terrain of inland mountainous areas may form a broken adaptive region. For each adaptive region, the average value / median value of the vertical ozone sinking value in the region is calculated, and the common characteristics of the strong sinking area, such as the stability of the inversion layer and the large accumulation slope of ozone, are combined with the common physical mechanisms of ozone change, such as the stability of the inversion layer and the large accumulation slope of ozone, to calibrate the ozone sinking contribution coefficient.

[0158] The ozone sinking contribution coefficient of the output target detection area is used to judge the influence strength of the vertical sinking process in the ozone transport channel in the subsequent process. The higher the coefficient, the easier the ozone is to sink and accumulate near the ground, and the greater the influence on the ground ozone concentration.

[0159] S6: Spatially superimpose the coastal strong penetration area and the inland weak penetration area with the contribution coefficient to form a hierarchical ozone transport channel.

[0160] In this embodiment, the spatial superposition of the coastal strong penetration area and the inland weak penetration area with the contribution coefficient to form a hierarchical ozone transport channel comprises:

[0161] Construct a composite penetration boundary map based on the spatial position data of the coastal strong penetration area and the inland weak penetration area;

[0162] Contribute to the coupling of the ozone sinking contribution coefficient and the composite penetration boundary map to obtain an ozone distribution area;

[0163] Divide the ozone distribution area according to the intensity to obtain the intensity level of the ozone distribution area;

[0164] Divide the ozone distribution area into different hierarchical ozone transport channels based on the intensity level.

[0165] Specifically, the above steps are the final synthesis link of the sea-land wind ozone transport channel analysis. Through spatial superposition, the coastal strong / weak penetration area divided in the early stage is fused with the ozone sinking contribution coefficient to identify ozone transport channels of different intensities and different spatial distributions, providing accurate basis for ozone pollution prevention and control and transmission path tracing.

[0166] Specifically, the coastal strong penetration area / inland weak penetration area is a geographical division in step S3. The strong penetration area has strong sea wind power and is the main channel for horizontal input of ozone; the weak penetration area has weak sea wind influence, and the horizontal transmission of ozone is mainly local generation.

[0167] The composite penetration boundary map is a GIS layer that fuses the spatial boundaries of strong / weak penetration areas, clearly marking the penetration influence range of sea breeze on different areas, and is the horizontal spatial basis for subsequent analysis of ozone transport paths.

[0168] Further, by extracting the spatial vector boundary of the coastal strong permeation area and the inland weak permeation area in step S3, such as polygon coordinates and geographical range;

[0169] In GIS software such as ArcGIS and QGIS, the boundary data of the strong and weak permeation areas are superimposed on the same base map containing the coastline, administrative division and other geographical references to construct a composite permeation boundary map.

[0170] The coastal strong permeation area is adjacent to the coastline and is distributed in a strip shape; the inland weak permeation area extends inland, and the boundary is affected by the terrain and city distribution and needs to be accurately drawn.

[0171] Combined with the actual terrain such as the trend of coastal mountains and the location of river mouths, the boundary rationality is verified. If the strong permeation area is cut off by mountains, the boundary needs to be corrected because the terrain will hinder the penetration of sea breeze, and the actual strong permeation area cannot extend across the mountains.

[0172] The corrected composite permeation boundary map is output to ensure consistency with the spatial characteristics of the sea-land wind physical action.

[0173] Specifically, the ozone distribution area is a geographical area that combines the horizontal penetration potential of sea breeze and the vertical deposition effect of ozone, reflecting the spatial distribution of ozone under the combined action of horizontal transmission and vertical deposition driven by sea-land breeze.

[0174] The ozone sinking contribution coefficient is mapped onto the composite permeation boundary map according to the grid, and the spatial correlation between the contribution coefficient and the permeation partition is established;

[0175] A coupling algorithm is designed, such as: ozone distribution intensity = sea breeze penetration intensity level x contribution coefficient weight, which combines the effects of horizontal penetration and vertical deposition.

[0176] For example, if a high contribution coefficient is superimposed on the coastal strong permeation area, then the ozone distribution intensity = strong penetration weight, such as 0.7 x high contribution coefficient, such as 0.8 = 0.56, representing an active ozone transport area.

[0177] Each grid of the composite permeation boundary map is traversed to calculate the coupled ozone distribution intensity, and an ozone distribution area layer such as a heat map is generated, with darker colors indicating stronger ozone distribution;

[0178] If the contribution coefficient is high in the coastal strong permeation area, the ozone horizontal input and vertical sinking double action will easily form a high value area; if the contribution coefficient is low in the inland weak permeation area, the ozone will mainly diffuse locally, and the distribution intensity will be weak.

[0179] Specifically, the intensity level is a level divided according to the coupling intensity value of the ozone distribution area, reflecting the activity degree of the ozone transport channel.

[0180] In detail, based on statistical methods such as natural breakpoint method and equal interval method, the coupling strength value distribution of the ozone distribution area is analyzed, the strength level threshold is set, under the influence of sea-land breeze, the strong penetration area along the coast is mostly strong level candidate, the weak penetration area in the inland is mostly weak level, and the threshold rationality needs to be verified in combination with historical ozone concentration data.

[0181] For each grid of the ozone distribution area, the strength level is divided according to the threshold value, and is marked as a strong transmission area, a medium transmission area and a weak transmission area, and an ozone distribution area graph with a strength level is output, and the ozone transmission potential of different levels is directly presented.

[0182] Specifically, the hierarchical ozone transmission channel is the final result, which is a geographical channel with different ozone transmission characteristics divided according to the strength level, and is a concrete presentation of the horizontal-vertical coordinated transmission of ozone under the driving of sea-land breeze.

[0183] In detail, by defining hierarchical rules, for example: first-level channel→strong penetration area+strong contribution coefficient+strong strength level (ozone horizontal input is more, vertical sinking is significant, and transmission is most active);

[0184] Second-level channel→strong penetration area edge+medium contribution coefficient+medium strength level (transmission capacity is second);

[0185] Third-level channel→weak penetration area+weak contribution coefficient+weak strength level (transmission is mainly local);

[0186] Further, the first-level channel corresponds to an ozone efficient transmission path dominated by sea breeze, such as a coastal city belt, sea breeze carries ozone+inversion sinking, and high concentration is easy to form; the third-level channel corresponds to an inland ozone weak transmission area, such as a mountainous area, sea breeze has weak influence, and ozone is difficult to input on a large scale.

[0187] As shown in FIG. Figure 2 It is a functional module diagram of an analysis system for sea-land breeze ozone transmission channel provided by an embodiment of the present application.

[0188] The analysis system 100 for sea-land breeze ozone transmission channel can be installed in an electronic device. According to the functions to be realized, the analysis system 100 for sea-land breeze ozone transmission channel can include a data acquisition module 101, a sea breeze influence module 102, an area division module 103, a temperature analysis module 104, an ozone contribution module 105 and a channel generation module 106. The modules of the present application can also be referred to as units, which refer to a series of computer program segments that can be executed by an electronic device processor and can complete a fixed function, and are stored in the memory of the electronic device.

[0189] In the embodiment, the functions of the modules / units are as follows:

[0190] The data acquisition module 101 is configured to acquire meteorological elements and ozone concentration of a target detection area;

[0191] The sea wind influence module 102 is configured to identify a sea wind duration in the meteorological elements, count an off-shore distance of the target detection area, and take an attenuation relationship between the sea wind duration and the off-shore distance as a sea wind penetration intensity value of the target detection area;

[0192] The area division module 103 is configured to extract an ozone change slope of the ozone concentration from an accumulation period to a dissipation period, and divide the target detection area into a coastal strong penetration area and an inland weak penetration area according to the sea wind penetration intensity value;

[0193] The temperature analysis module 104 is configured to determine an inversion layer height attenuation rate of the target detection area based on a temperature vertical gradient in the meteorological elements;

[0194] The ozone contribution module 105 is configured to superimpose the inversion layer height attenuation rate and the ozone change slope as an ozone sinking contribution coefficient of the target detection area;

[0195] The channel generation module 106 is configured to spatially superimpose the coastal strong penetration area and the inland weak penetration area with the contribution coefficient to form a hierarchical ozone transmission channel.

[0196] In several embodiments provided in the present application, it should be understood that the disclosed method and system can be implemented in other ways. For example, the system embodiments described above are only illustrative, for example, the division of the modules is only a logical function division, and another division mode can be used in actual implementation.

[0197] The modules described as separate components can or can not be physically separated, and the components displayed as modules can or can not be physical units, that is, they can be located in one place or distributed on multiple network units. Part or all of the modules can be selected to achieve the purpose of the embodiment according to actual needs.

[0198] In addition, the function modules in each embodiment of the present application can be integrated in one processing unit, or each unit can exist physically, or two or more units can be integrated in one unit. The integrated unit can be realized in the form of hardware or in the form of hardware plus software function modules.

[0199] It is obvious for those skilled in the art that the present application is not limited to the details of the above exemplary embodiments, and the present application can be implemented in other specific forms without departing from the spirit or essential characteristics of the present application.

[0200] The embodiments of the present application can acquire and process related data based on artificial intelligence technology. The artificial intelligence is a theory, method, technology and application system for simulating, extending and expanding human intelligence by using a digital computer or a machine controlled by a digital computer, perceiving an environment, acquiring knowledge and using the knowledge to obtain optimal results.

[0201] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present application and are not limiting. Although the present application has been described in detail with reference to the preferred embodiments, it should be understood by those skilled in the art that the technical solutions of the present application can be modified or replaced equivalently without departing from the spirit and scope of the technical solutions of the present application.

Claims

1. A sea-land breeze ozone transport channel analysis method, characterized by, The method comprises: S1: obtaining meteorological elements and ozone concentration of a target detection area; S2: identifying a sea wind duration in the meteorological elements, counting an off-shore distance of the target detection area, and taking an attenuation relationship between the sea wind duration and the off-shore distance as a sea wind penetration intensity value of the target detection area; The sea wind duration and the off-shore distance are taken as the sea wind penetration intensity value of the target detection area, comprising: time-weighted evaluation of the sea wind duration based on a preset seasonal wind field scale parameter to obtain a time-weighted intensity of the sea wind duration; analyzing a spatial attenuation intensity value of the off-shore distance; synthesizing the time-weighted intensity and the spatial attenuation intensity value to obtain the sea wind penetration intensity value of the target detection area; The time-weighted intensity and the spatial attenuation intensity value are synthesized to obtain the sea wind penetration intensity value of the target detection area, comprising: ， wherein: is a sea wind penetration intensity value, is the time weighted intensity, is a terrain penetration transfer function, is an atmospheric stability factor, is the spatial decay intensity value, is an absolute value of an elevation difference, is a maximum terrain roughness; S3: extracting an ozone change slope of the ozone concentration from an accumulation period to a dissipation period, and dividing the target detection area into a coastal strong penetration area and an inland weak penetration area according to the sea wind penetration intensity value; S4: determining an inversion layer height attenuation rate of the target detection area based on a temperature vertical gradient in the meteorological elements; S5: superimposing the inversion layer height attenuation rate and the ozone change slope as an ozone sinking contribution coefficient of the target detection area; S6: spatially superimposing the coastal strong penetration area and the inland weak penetration area with the contribution coefficient to form a hierarchical ozone transport channel.

2. A sea-breeze ozone transport channel analysis method as claimed in claim 1, characterized in that, The meteorological elements and the ozone concentration of the target detection area are obtained, comprising: continuously collecting meteorological observation original data of the target detection area and synchronously obtaining regional ozone information of the target detection area; eliminating abnormal data in the meteorological observation original data and the regional ozone information to obtain the meteorological elements and the ozone concentration of the target detection area.

3. A method of sea-breeze ozone transport channel analysis as claimed in claim 1, wherein, The sea wind duration in the meteorological elements is identified, and the off-shore distance of the target detection area is counted, comprising: identifying a wind direction transition event of a wind field time domain sequence in the meteorological elements, and marking a duration of the wind direction transition event as a sea wind duration; determining an off-shore distance of the target detection area based on a coordinate reference positioning of the target detection area.

4. A method of sea-breeze ozone transport channel analysis as claimed in claim 1, wherein, The target detection area is divided into a coastal strong penetration area and an inland weak penetration area according to the sea wind penetration intensity value, comprising: taking the sea wind penetration intensity value as a penetration threshold value of the target detection area; based on the penetration threshold value, marking a geographical spatial distribution of the target detection area; verifying the spatial continuity of the result of the marking, and dividing the result that passes the verification into the coastal strong penetration area and the inland weak penetration area.

5. A method of sea-breeze ozone transport channel analysis as claimed in claim 1, wherein, The inversion layer height attenuation rate of the target detection area is determined based on the temperature vertical gradient in the meteorological elements, comprising: extracting a height boundary feature of the temperature vertical gradient, and synchronously analyzing a height change difference of the height boundary feature in a continuous time sequence; The height change difference obtained by analysis is aggregated as an inversion layer height decay rate of the target detection area.

6. A method of sea-breeze ozone transport channel analysis as claimed in claim 1, wherein, The inversion layer height decay rate and the ozone change slope are superimposed as an ozone sinking contribution coefficient of the target detection area, including: Quantitative analysis of the ozone vertical distribution relationship of the change characteristics of the inversion layer height decay rate and the distribution characteristics of the ozone change slope; Performing contribution weight generation on the ozone vertical distribution relationship to obtain an ozone vertical sinking value; Matching the ozone vertical sinking value corresponding to the adaptive area to obtain the ozone sinking contribution coefficient of the target detection area.

7. A sea-breeze ozone transport channel analysis method as claimed in claim 6, characterized in that, The coastal strong permeation area and the inland weak permeation area are spatially superimposed with the contribution coefficient to form a hierarchical ozone transmission channel, including: Constructing a composite permeation boundary map based on the spatial position data of the coastal strong permeation area and the inland weak permeation area; Coupling the contribution degree of the ozone sinking contribution coefficient and the composite permeation boundary mapping to obtain an ozone distribution area; Dividing the ozone distribution area according to the intensity to obtain the intensity level of the ozone distribution area; Dividing the ozone distribution area into different hierarchical ozone transmission channels based on the intensity level.

8. A sea-land breeze ozone transport channel analysis system characterized by, The system includes: A data acquisition module for acquiring meteorological elements and ozone concentration of a target detection area; A sea wind influence module for identifying the duration of sea wind in the meteorological elements, counting the distance from shore of the target detection area, and taking the decay relationship between the duration of sea wind and the distance from shore as the sea wind permeation intensity value of the target detection area; The sea wind duration and the distance from shore are taken as the sea wind permeation intensity value of the target detection area, including: Time-weighted evaluation of the duration of sea wind based on a preset seasonal wind field scale parameter to obtain a time-weighted intensity of the duration of sea wind; Analyzing the spatial decay intensity value of the distance from shore; Synthesizing the time-weighted intensity and the spatial decay intensity value to obtain the sea wind permeation intensity value of the target detection area; The time-weighted intensity and the spatial decay intensity value are synthesized to obtain the sea wind permeation intensity value of the target detection area, including: ， wherein: is a sea wind penetration intensity value, is the time weighted intensity, is a terrain penetration transfer function, is an atmospheric stability factor, is the spatial decay intensity value, is an absolute value of an elevation difference, is a maximum terrain roughness; A region division module for extracting the ozone change slope of the ozone concentration from the accumulation period to the dissipation period, and dividing the target detection area into a coastal strong permeation area and an inland weak permeation area according to the sea wind permeation intensity value; A temperature analysis module for determining the inversion layer height decay rate of the target detection area based on the temperature vertical gradient in the meteorological elements; An ozone contribution module for superimposing the inversion layer height decay rate and the ozone change slope as an ozone sinking contribution coefficient of the target detection area; A channel generation module for spatially superimposing the coastal strong permeation area and the inland weak permeation area with the contribution coefficient to form a hierarchical ozone transmission channel.

Images