Sea and land wind ozone transmission channel analysis method and system

Abstract

The invention relates to the technical field of meteorological analysis, and discloses a sea-land wind ozone transmission channel analysis method and system, and the method comprises the steps: obtaining meteorological elements and ozone concentration of a target detection region; the sea wind duration in the meteorological elements is identified, the offshore distance of the target detection area is counted, and the attenuation relation between the sea wind duration and the offshore distance serves as the sea wind permeability strength value of the target detection area; 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 permeation area and an inland weak permeation area according to the sea wind permeation intensity value; determining the thermal inversion layer height attenuation rate of the target detection area based on the temperature vertical gradient in the meteorological elements; and superposing the thermal inversion layer height attenuation rate and the ozone change slope into an ozone sinking contribution coefficient of the target detection area. The method can solve the problem that refined ozone source tracing is difficult to support in the prior art.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of meteorological analysis, and particularly relates 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 a 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 attenuation relationship between the sea wind duration and the off-shore distance, and does not consider the influence of factors such as terrain undulation and atmospheric stability on the penetration ability, resulting in a fuzzy division of the ozone transport boundary between coastal and inland areas, and an inability to accurately identify strong / weak penetration areas, which makes it difficult to support refined ozone source tracing.

[0003] At the same time, the existing technology lacks coupling research on the inversion layer height attenuation law and the ozone concentration change trend when analyzing the vertical transport of ozone, and only analyzes the temperature vertical gradient or the change of ozone concentration in isolation, without quantifying the contribution coefficient of the inversion layer to the ozone sinking. 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 ozone pollution migration path, 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.

[0005] To achieve the above-mentioned purpose, the present application provides a sea-land wind ozone transport channel analysis method, which comprises: S1: obtaining meteorological elements and ozone concentration of a target detection area; S2: identifying the sea wind duration in the meteorological elements, and counting the off-shore distance of the target detection area, taking the attenuation relationship between the sea wind duration and the off-shore distance as the sea wind penetration strength value of the target detection area; S3: extracting the ozone change slope of the ozone concentration from the accumulation period to the dissipation period, and dividing the target detection area into coastal strong penetration area and inland weak penetration area according to the sea wind penetration strength value; S4: determining the inversion layer height attenuation rate of the target detection area based on the temperature vertical gradient in the meteorological elements; S5: superimposing the inversion layer height attenuation rate and the ozone change slope as the 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.

[0006] Preferably, the meteorological elements and ozone concentration of the target detection area are obtained, comprising: Continuously collecting meteorological observation raw data of the target detection area, and synchronously obtaining regional ozone information of the target detection area; Eliminating 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.

[0007] Preferably, the duration of the sea breeze in the meteorological elements is identified, and the off-shore distance of the target detection area is counted, comprising: 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 a sea breeze duration; Determining the off-shore distance of the target detection area based on a coordinate reference position of the target detection area.

[0008] Preferably, the sea breeze penetration intensity value of the target detection area is obtained by taking the attenuation relationship between the sea breeze duration and the off-shore distance as the sea breeze penetration intensity value, comprising: Time-weighted evaluation of the sea breeze duration based on a preset seasonal wind field scale parameter to obtain a time-weighted intensity of the sea breeze 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 breeze penetration intensity value of the target detection area.

[0009] Preferably, the sea breeze penetration intensity value of the target detection area is obtained by synthesizing the time-weighted intensity and the spatial attenuation intensity value, comprising:

[0010] Wherein: is a sea breeze penetration intensity value, is the time-weighted intensity, is a terrain penetration transfer function, is an atmospheric stability factor, is the spatial attenuation intensity value, is an absolute value of an altitude difference, is a maximum terrain undulation.

[0011] Preferably, the target detection area is divided into a coastal strong penetration area and an inland weak penetration area according to the sea breeze penetration intensity value, comprising: Taking the sea breeze penetration intensity value as a penetration threshold value of the target detection area; partition marking a geospatial distribution of the target detection region based on the penetration threshold value; checking spatial continuity of a result of the partition marking, and partitioning a result passing the checking into a coastal strong penetration area and an inland weak penetration area.

[0012] Preferably, the determining the inversion layer height decay rate of the target detection region based on the temperature vertical gradient in the meteorological element comprises: extracting a height boundary feature of the temperature vertical gradient, and synchronously analyzing a height change difference of the height boundary feature within a continuous time sequence; aggregating the height change difference obtained through the analysis into the inversion layer height decay rate of the target detection region.

[0013] Preferably, the superimposing the inversion layer height decay rate and the ozone change slope into the ozone sinking contribution coefficient of the target detection region comprises: quantitative analysis of an ozone vertical distribution relationship between a change feature of the inversion layer height decay rate and a distribution feature of the ozone change slope; generating a contribution weight for the ozone vertical distribution relationship to obtain an ozone vertical sinking value; matching an adaptive region corresponding to the ozone vertical sinking value to obtain the ozone sinking contribution coefficient of the target detection region.

[0014] Preferably, the spatial superimposition of the coastal strong penetration area and the inland weak penetration area and the contribution coefficient to form a hierarchical ozone transport channel comprises: constructing a composite penetration boundary map based on spatial position data of the coastal strong penetration area and the inland weak penetration area; coupling a contribution degree of the ozone sinking contribution coefficient and the composite penetration boundary mapping map to obtain an ozone distribution region; dividing the ozone distribution region according to intensity to obtain an intensity level of the ozone distribution region; dividing the ozone distribution region into different hierarchical ozone transport channels based on the intensity level.

[0015] A sea-land wind ozone transport channel analysis system, the system comprising: a data acquisition module for acquiring meteorological elements and ozone concentration of a target detection region; a sea wind influence module for identifying a sea wind duration in the meteorological elements, counting an off-shore distance of the target detection region, 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 region; The regional division module is configured to extract an ozone change slope from the ozone concentration from the accumulation period to the dissipation period, and divide the target detection region into a coastal strong penetration area and an inland weak penetration area according to the sea wind penetration intensity value. The temperature analysis module is configured to determine an inversion layer height decay rate of the target detection region based on a temperature vertical gradient in the meteorological element. The ozone contribution module is configured to superimpose the inversion layer height decay rate and the ozone change slope as an ozone sinking contribution coefficient of the target detection region. The channel generation module 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 transport channel.

[0016] Advantages 1. The sea wind penetration intensity value is converted from the attenuation relationship between the sea wind duration and the offshore distance, and the penetration ability is synthesized by introducing parameters such as the terrain penetration transfer function and the atmospheric stability factor, which breaks through the limitation of the traditional method relying on a single element, can accurately divide the coastal strong penetration area and the inland weak penetration area, clearly define the ozone transport boundary of different regions, and provide a more effective basis for identifying the ozone cross-regional transport path.

[0017] 2. The inversion layer height decay rate and the ozone change slope are superimposed as the ozone sinking contribution coefficient, and the hierarchical ozone transport channel is formed by spatial superposition, which realizes the coupling analysis of the ozone horizontal transport and the vertical deposition process, makes up for the defects of the prior art in isolated analysis of meteorological elements or ozone concentration, can not only quantify the influence weight of the inversion layer on the ozone sinking, but also can distinguish the transport channel according to the intensity level, significantly improves the analysis accuracy of the ozone pollution migration rule, and provides more accurate technical support for regional ozone pollution joint prevention and control. BRIEF DESCRIPTION OF DRAWINGS

[0018] Figure 1 A flowchart of a sea-land wind ozone transport channel analysis method provided by an embodiment of the present application is shown. 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. DETAILED DESCRIPTION

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

[0020] Embodiments of the present application provide a sea-land breeze ozone transport channel analysis method and system. The execution subject of the sea-land breeze 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. that can be configured to execute the method provided by the embodiments of the present application. In other words, the sea-land breeze 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 that provides 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 basic cloud computing services such as big data and artificial intelligence platforms, etc.

[0021] Referring to Figure 1 FIG. 1 is a flowchart of a sea-land breeze ozone transport channel analysis method provided by an embodiment of the present application. In this embodiment, the sea-land breeze ozone transport channel analysis method includes the following steps. S1: Obtain meteorological elements and ozone concentration of a target detection area.

[0022] In this embodiment, the obtaining of the meteorological elements and ozone concentration of the target detection area includes the following steps. Continuously collect meteorological observation raw data of the target detection area, and synchronously obtain regional ozone information of the target detection area; 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.

[0023] 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 breeze environment, the change of the wind field is the key to identifying the sea breeze. During the day, the sea temperature is lower than the land temperature, and the sea breeze blows from the sea to the land. At night, the situation is reversed to form a land breeze. The temperature vertical gradient affects the formation of the inversion layer, which hinders the vertical diffusion of ozone.

[0024] The ozone concentration refers to the content of ozone in the air of the target detection area. Under the action of the 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 the photochemical reaction is weakened, the ozone gradually decomposes, and the concentration decreases.

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

[0026] The synchronous acquisition of regional ozone information is to collect ozone concentration data in the same time period and the same area as the meteorological data collection by using ozone monitors and other equipment. The matching of meteorological elements and ozone concentration data in time and space is ensured to analyze the correlation between the two, such as the relationship between sea wind intensity and ozone transport volume.

[0027] Further, the abnormal data is removed. The abnormal data of meteorological observation raw data may include sudden change of wind speed to 0 or far beyond the normal range due to equipment failure, and irregular and sharp change of wind direction.

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

[0029] Through data verification algorithms such as threshold method, mean standard deviation method and the like, these abnormal data are identified and removed to ensure that the obtained meteorological elements and ozone concentration data are true and reliable, and to provide accurate basic data for subsequent analysis.

[0030] S2: identifying the duration of sea wind in the meteorological elements, counting the offshore distance of the target detection area, and taking the attenuation relationship between the duration of sea wind and the offshore distance as the sea wind penetration intensity value of the target detection area.

[0031] In this embodiment, the identification of the duration of sea wind in the meteorological elements and the counting of the offshore distance of the target detection area include: identifying the wind direction change event of the wind field time sequence in the meteorological elements, and marking the duration of the wind direction change event as the duration of sea wind; determining the offshore distance of the target detection area based on the coordinate reference positioning of the target detection area.

[0032] Specifically, the wind field time sequence is the wind direction vector data of the target detection area recorded in time sequence, which 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 change of wind field with day and night alternation, and the wind field time sequence can capture this rule.

[0033] The wind direction change event refers to the critical process of the wind direction changing 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, so sea wind begins to dominate, and the wind direction change event is the identifier of the sea wind start-continue-end cycle.

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

[0035] The wind direction change recognition is to use an algorithm such as a sliding window + threshold judgment to detect the time point at which the wind direction changes from a land wind feature such as a wind direction pointing to the ocean direction to a sea wind feature in combination with the definition of the regional sea-land geographical orientation, and to mark the start and end time of the change.

[0036] The duration calibration is to calculate the time span from the start to the end of the wind direction change event, that is, the duration of the sea wind. In the physical environment, the duration of the sea wind is affected by the sea-land temperature difference intensity and the terrain (such as the blocking of coastal mountains). When the sea-land temperature difference is large on a sunny day, the duration of the sea wind is longer.

[0037] Specifically, the coordinate reference positioning is based on the geographical coordinates of the target detection area, and is associated with the coordinate reference of the coastline, such as marking the coastline as a reference line through a GIS map. The coastline is the geographical boundary of the sea-land wind, and the distance from the coast directly affects the path length and attenuation degree of the sea wind transporting ozone.

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

[0039] The off-shore distance calculation is to use the GIS spatial analysis function to calculate the perpendicular distance (i.e. the off-shore distance) between the target detection area and the coastline. The farther the off-shore distance is, the stronger the influence of the land underlying surface such as the urban heat island and vegetation cover on the sea wind, and the wind field intensity and ozone transmission capacity attenuate with the off-shore distance.

[0040] In this embodiment, the sea wind duration and the attenuation relationship of the off-shore distance are taken as the sea wind penetration intensity value of the target detection area, which comprises: 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 the 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.

[0041] Specifically, the seasonal wind field scaling parameter is the difference in sea-land thermal difference in different seasons, the background wind field, such as the intensity of monsoon, and the preset parameter is obtained by long-term meteorological observation statistics, such as the sea-land temperature difference is large in summer, the basic weight coefficient of sea wind is higher, which is used to quantify the influence of season on the duration of sea wind. For example, the scaling parameter is set to 1.2 in summer and 0.8 in spring, which reflects that the sea wind in summer is more significantly driven by thermal force.

[0042] The time-weighted intensity is the quantitative value of the duration of sea wind after the seasonal correction, reflecting the effective intensity of the duration of sea wind under the background of season.

[0043] In detail, the weighted calculation is calculated by the formula, for example: time-weighted intensity = sea wind duration × seasonal wind field scaling parameter.

[0044] In summer, the land warms up quickly, the sea-land temperature difference is large, and the duration of sea wind is longer, and the scaling parameter further strengthens this seasonal characteristic, so that the evaluation is more in line with the actual sea-land wind dynamic mechanism, such as the duration of sea wind in summer is 20 hours, the scaling parameter is 1.2, and then the weighted value is 24; the scaling parameter is 0.7 in winter, and the weighted value of 15 hours of sea wind is 10.5.

[0045] Specifically, the spatial decay intensity value is a quantitative value of the decay of sea wind intensity with the increase of offshore distance, reflecting the weakening degree of geography on sea wind.

[0046] Based on the physical law, the deeper the sea wind penetrates into the land, the faster it decays due to the influence of the friction of the underlying surface of the land and the reduction of the sea-land temperature difference caused by the warming of the land, and a spatial decay model is established, such as exponential decay: spatial decay intensity value = initial intensity × e^(-k × offshore distance), k is the decay coefficient, which is fitted by historical data.

[0047] The offshore distance of the target area is substituted into the model to calculate the spatial decay intensity value. For example, k = 0.05 at 10 km offshore, the initial intensity is 100, and the decay is 100 × e^(-0.5) ≈ 60.7; at 50 km inland, it is 100 × e^(-2.5) ≈ 8.2, which reflects that the farther the offshore distance, the weaker the spatial penetration of sea wind, and the lower the ozone transport power.

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

[0049] Further, the time-weighted intensity reflects how long the sea breeze can drive the transmission, and the spatial decay intensity reflects how far the transmission becomes weak. The synthesized value 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, the spatial decay intensity is 60.7, and the synthesis is 1456.8. The higher the value, the more effectively the sea breeze can transmit ozone in the area, which is the core indicator for determining whether the ozone transmission channel is active.

[0050] More specifically, the sea-land breeze is a local circulation driven by the sea-land thermal difference. The seasonal effect is the sea-land temperature difference, and the distance from the coast affects the degree of interference of the sea breeze by the land. Through time weighting + spatial decay + synthesis, the ability of the sea breeze to carry ozone transmission in the real sea-land environment is accurately quantified, providing key parameters for subsequent identification of ozone transmission channels and analysis of transmission rules.

[0051] 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, including:

[0052] 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 elevation difference, is the maximum terrain relief.

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

[0054] Specifically, the time-weighted intensity reflects the effective dynamic intensity of the duration of the sea breeze in the time dimension of seasons, day and night, etc. It is affected by the seasonal variation of the sea-land thermal difference.

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

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

[0057] The spatial decay intensity value reflects the degree of sea wind decay with increasing distance from the shore. The farther the distance from the shore, the greater the sea wind power, which is calculated by the spatial decay step of the distance from the shore, reflecting the weakening of sea wind by geographical space.

[0058] In detail, the exponential function simulates the decay law of sea wind power with atmospheric stability + distance from the shore, or the greater, the closer to 0, the closer to 1.

[0059] For example: coastal areas small, if the atmosphere is unstable, ≈1, the correction term ≈0, but the actual coastal sea wind power is strong, indicating that the combined action is needed; inland areas large, if the atmosphere is stable, ≈0, the correction term ≈1, reflecting the decay of sea wind power to the limit, approaching pure land wind, and the weak influence of sea wind on ozone transport.

[0060] Specifically, the formula for calculating the terrain penetration transfer function is:

[0061] where: is the terrain penetration transfer function, is the coastal curvature coefficient.

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

[0063] This design makes the formula have double precise matching ability: when the coastal terrain is flat ( ≈1), it degenerates into an exponential decay function, which conforms to the energy dissipation characteristics of gentle climbing; when facing curved coastlines ( >1), the increase of the value significantly improves the decay rate of the function through the multiplication factor , which accurately quantifies the additional kinetic energy loss of sea wind flow lines due to the tortuosity of the coastline.

[0064] Especially for the problem of back slope vortex mutation that the traditional model completely fails, when the target point is close to the terrain peak, the function value sharply rises to positive infinity at the critical point (x=0.95), causing The cliff-like drop is realized.

[0065] This mathematical mutation behavior corresponds exactly to the phenomenon of flow separation caused by sudden change of air pressure after the airflow crosses the ridge in reality.

[0066] In more detail, by normalizing the parameter Eliminate the area dependence of the absolute value of the altitude, with global application range.

[0067] By the coefficient Generated by the fractal dimension of the coastline, when the tortuosity of the coastline rises, such as the fractal dimension of the bedrock coast 1.6→1.8, The value increases from 0.4 to 0.2, and under the same topographic conditions, the penetration distance of sea breeze is shortened, accurately reflecting the blocking effect of the spiral-shaped bay coast on the airflow.

[0068] When <0.3, ≈1 (airflow unobstructed penetration), >0.6, there is only one inflection point (corresponding to the triggering line of the canyon acceleration effect), >0.95, the exponential decay mechanism is automatically activated, which enables the early warning system to respond in advance to mountain vortex disasters.

[0069] 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 breeze penetration intensity value.

[0070] In this embodiment, the target detection area is divided into coastal strong penetration area and inland weak penetration area according to the sea breeze penetration intensity value, comprising: The sea breeze penetration intensity value is used as the penetration threshold of the target detection area; Based on the penetration threshold, the geographical spatial distribution of the target detection area is marked by zoning; Verify the spatial continuity of the zoning marking result, and divide the result that passes the verification into coastal strong penetration area and inland weak penetration area.

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

[0072] The generation and decomposition of ozone are affected by light, temperature, concentrations of volatile organic compounds (VOCs) and nitrogen oxides (NOx) and other precursor substances. 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, the ozone gradually decomposes, and it enters the dissipation period. The greater the change slope of ozone, the more intense the change of ozone concentration from the accumulation period to the dissipation period.

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

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

[0075] The change slope of ozone concentration with time in these two stages is calculated. Linear regression algorithm can be used to fit the ozone concentration data at multiple time points in the accumulation period and the dissipation period, and the slope of the fitted straight line is the change slope of ozone.

[0076] Specifically, the sea wind penetration intensity value is a comprehensive consideration of the duration of sea wind, the distance from the shore, the terrain and the 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 is the key basis for dividing the coastal strong penetration area and the inland weak penetration area.

[0077] 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 effective penetration of sea wind, and is used to distinguish the degree of influence of sea wind on different areas.

[0078] 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. For each geographic unit, compare its sea wind penetration intensity value with the penetration threshold.

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

[0080] Finally, check whether there are fragmented or isolated parts in the marked area to ensure that the partitioning results are spatially continuous and reasonable. If there are discontinuous areas, it may be due to factors such as data anomalies or special terrain, and further analysis and processing are required.

[0081] For example, discontinuous boundaries can be adjusted by checking the marking status of adjacent areas. Ultimately, the results that pass the verification will be formally divided into coastal strong infiltration areas and inland weak infiltration areas. This division can more accurately reflect the different effects of sea breeze on ozone transmission in different regions.

[0082] S4: Determine the inversion layer height attenuation rate of the target detection area based on the vertical temperature gradient in the meteorological element.

[0083] In this embodiment, determining the inversion layer height attenuation rate of the target detection area based on the vertical temperature gradient in the meteorological element includes: Extracting the height boundary feature of the vertical temperature gradient, and simultaneously analyzing the height change difference of the height boundary feature in a continuous time series; The height change differences obtained through analysis are aggregated into the inversion layer height attenuation rate of the target detection area.

[0084] Specifically, the team focused on calculating the attenuation rate at the inversion layer's height. The inversion layer (a layer where the temperature rises with altitude) inhibits the vertical diffusion of ozone, while the thermal differences between sea and land breezes dynamically alter the inversion layer's height. By analyzing the relationship between vertical temperature gradient, changes in inversion layer height, and the attenuation rate, the team provided a basis for determining the extent of vertical ozone transmission obstruction. The higher the attenuation rate, the more unstable the inversion layer, and the easier it is for ozone to break through and diffuse.

[0085] Specifically, the vertical temperature gradient is the rate of change of atmospheric temperature with altitude. In a sea-land breeze environment, land warms rapidly during the day, resulting in a large vertical temperature gradient near the surface (prone to unstable stratification). At night, the temperature difference between land and sea is small, making inversions (vertical temperature gradient ≤ 0) more likely to form near the surface.

[0086] The height boundary characteristics are the top and bottom heights of the inversion layer. They are the critical heights where the vertical temperature gradient changes from an inversion characteristic to a normal stratification, reflecting the distribution range of the inversion layer in the vertical space.

[0087] Furthermore, vertical sounding data collection is to obtain temperature data of continuous altitude layers in the target detection area with the help of sounding balloons, weather radars or vertical observation stations.

[0088] Driven by sea and land breezes, the vertical distribution of near-surface temperature changes dynamically with day and night and the position of land and sea, requiring high-density vertical sampling to capture the inversion layer.

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

[0090] In the embodiment, the superimposing the inversion layer height decay rate and the ozone change slope as the ozone sinking contribution coefficient of the target detection area comprises: 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; generating a contribution weight of 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.

[0091] 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, thereby providing a basis for identifying the vertical sinking effect of the ozone transport channel.

[0092] Specifically, the ozone vertical distribution relationship is the correlation law 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 likely to sink and accumulate when the inversion layer is pressed down.

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

[0094] The sea-land wind is a local circulation, and the inversion layer and the ozone change characteristics of different grids are significantly different, so that fine gridding is required.

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

[0096] Statistical methods such as Pearson correlation analysis and multiple linear regression are used to analyze the correlation between the inversion layer height decay rate change and the ozone concentration change slope in the vertical profile, and an ozone vertical distribution relationship model is established, such as: the ozone change slope = a x the inversion layer height decay rate + b, a and b are fitting coefficients.

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

[0098] The stronger the inversion layer, the higher the weight of the capping effect on ozone sinking; the weaker the inversion layer, the higher the weight of the releasing effect on ozone diffusion.

[0099] The ozone vertical sinking value is a quantitative value of the net movement intensity of ozone in the vertical direction, sinking to the near surface or diffusing to the high altitude. A positive sinking value indicates that ozone converges to the near surface, and a negative sinking value indicates that ozone diffuses to the high altitude.

[0100] Further, based on the ozone vertical distribution relationship, the response degree of the ozone change slope is observed by sensitivity analysis, such as changing the inversion layer height decay rate, to calibrate the contribution weight of the inversion layer height decay rate.

[0101] By substituting the weight and the parameter into each grid, the ozone vertical sinking value is calculated to reflect the net movement trend of ozone in the vertical direction in the grid.

[0102] Specifically, the adaptive region is a continuous region in the geographical space with similar characteristics of ozone vertical sinking value, such as a strong sinking region and a weak diffusion region, which is jointly affected by sea-land wind, terrain, and inversion layer. The ozone transport channel characteristics of the adaptive region have consistency.

[0103] The ozone sinking contribution coefficient is the final quantitative result, which reflects the comprehensive contribution degree of ozone sinking to the near surface in the target detection region under the synergistic action of the inversion layer and the ozone concentration change. The higher the coefficient, the more significant the ozone sinking effect.

[0104] In detail, the spatial clustering algorithm, such as K-Means clustering and DBSCAN, is used to cluster the ozone vertical sinking values of all grids to divide the adaptive regions, such as strong sinking regions, weak sinking regions, and diffusion regions.

[0105] The coastal area may form a continuous strong sinking region due to the superposition of sea wind and inversion, and the inland mountainous area may form a broken adaptive region due to complex terrain. For each adaptive region, the average value / median of the ozone vertical sinking value in the region is calculated, and the common characteristics of the physical mechanism of the inversion layer and ozone change in the region are combined, such as the common characteristics of the strong sinking region: stable inversion layer and large ozone accumulation slope, to calibrate the ozone sinking contribution coefficient.

[0106] The ozone sinking contribution coefficient of the target detection region is output, which is used to judge the influence intensity of the vertical sinking process in the ozone transport channel in the subsequent process. The higher the coefficient, the easier the ozone accumulates in the near surface, which affects the ground ozone concentration.

[0107] S6: Spatially overlaying the coastal strong penetration region and the inland weak penetration region with the contribution coefficient to form a hierarchical ozone transport channel.

[0108] In this embodiment, the spatial superposition of the coastal strong penetration area and the inland weak penetration area with the contribution coefficient forms a hierarchical ozone transport channel, which includes: Constructing a composite penetration boundary map based on the spatial position data of the coastal strong penetration area and the inland weak penetration area; Coupling the ozone sinking contribution coefficient and the composite penetration boundary map 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; Based on the intensity level, the ozone distribution area is divided into different hierarchical ozone transport channels.

[0109] 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 with different intensities and different spatial distributions, providing accurate basis for ozone pollution prevention and control and transport path tracing.

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

[0111] 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 wind on different areas, and is the horizontal spatial basis for subsequent analysis of ozone transport paths.

[0112] Further, the spatial vector boundaries of the coastal strong penetration area and the inland weak penetration area in step S3 are extracted, such as polygon coordinates and geographical range. In GIS software such as ArcGIS and QGIS, the boundary data of the strong and weak penetration areas are superimposed on the same base map containing coastlines, administrative divisions, and other geographical references to construct a composite penetration boundary map.

[0113] The coastal strong penetration area is adjacent to the coastline and is distributed in a strip shape. The inland weak penetration area extends inland, and the boundary is affected by topography and urban distribution, which needs to be accurately drawn.

[0114] 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 penetration area is truncated by mountains, the boundary needs to be corrected because the terrain will hinder the penetration of sea wind and the actual strong penetration area cannot extend across the mountains.

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

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

[0117] The ozone sinking contribution coefficient is mapped onto the composite penetration boundary map by aligning the grid with the grid of the composite penetration boundary map, establishing the spatial correlation of the contribution coefficient-penetration partition; The coupling algorithm is designed, such as: ozone distribution intensity = sea breeze penetration intensity level x contribution coefficient weight, combining the effects of horizontal penetration and vertical deposition.

[0118] For example, in the coastal strong penetration area, the high contribution coefficient is superimposed, and the ozone distribution intensity = strong penetration weight, such as 0.7 x high contribution coefficient, such as 0.8 = 0.56, representing the ozone transmission active area.

[0119] Traverse each grid of the composite penetration boundary map, calculate the ozone distribution intensity after coupling, and generate the ozone distribution area layer, such as the heat map, the deeper the color, the stronger the ozone distribution; If the contribution coefficient is high in the coastal strong penetration area, the ozone horizontal input + vertical sinking double action, it is easy to form a high value area; if the contribution coefficient is low in the inland weak penetration area, the ozone is mainly local diffusion, and the distribution intensity is weak.

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

[0121] In detail, based on statistical methods such as natural breakpoints method and equal interval method, analyze the distribution of coupling intensity value of ozone distribution area, set the intensity level threshold, under the influence of sea-land breeze, the coastal strong penetration area is mostly strong level candidate, the inland weak penetration area is mostly weak level, the threshold rationality needs to be verified combined with historical ozone concentration data.

[0122] For each grid of the ozone distribution area, divide the intensity level according to the threshold, and label it as strong transmission area, medium transmission area, and weak transmission area, output the ozone distribution area map with intensity level, and intuitively present the ozone transmission potential of different levels.

[0123] Specifically, the graded ozone transmission channel is the final result, which is a geographical channel with different ozone transmission characteristics divided by intensity level, and is a concrete representation of ozone horizontal-vertical coordinated transmission driven by sea-land breeze.

[0124] In detail, by defining the grading rules, for example: first level channel → strong penetration area + strong contribution coefficient + strong intensity level (ozone horizontal input is more, vertical sinking is significant, transmission is most active); Second level channel → strong penetration area edge + medium contribution coefficient + medium intensity level (transmission capacity is second); The third channel is a weak penetration zone, a weak contribution coefficient and a weak intensity level (transmission is mainly local). Further, the first channel corresponds to an ozone efficient transmission path dominated by sea breeze, such as a coastal city belt, and sea breeze carries ozone + inversion sinking, which is easy to form high concentration.

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

[0126] The sea-land wind ozone transmission channel analysis system 100 can be installed in an electronic device. According to the functions implemented, the sea-land wind ozone transmission channel analysis system 100 can include a data acquisition module 101, a sea wind influence module 102, a region 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 fixed functions, which are stored in the memory of the electronic device.

[0127] In the present embodiment, the functions of each module / unit are as follows: The data acquisition module 101 is used to acquire meteorological elements and ozone concentration of a target detection area. The sea wind influence module 102 is used to identify the duration of sea wind in the meteorological elements, to count the distance from the shore of the target detection area, and to take the decay relationship between the duration of sea wind and the distance from the shore as the sea wind penetration intensity value of the target detection area. The region division module 103 is used to extract the ozone change slope of the target detection area from the accumulation period to the dissipation period, and to divide the target detection area into a coastal strong penetration zone and an inland weak penetration zone according to the sea wind penetration intensity value. The temperature analysis module 104 is used to determine the inversion layer height decay rate of the target detection area based on the temperature vertical gradient in the meteorological elements. The ozone contribution module 105 is used to superimpose the inversion layer height decay rate and the ozone change slope as the ozone sinking contribution coefficient of the target detection area. The channel generation module 106 is used to spatially superimpose the coastal strong penetration zone and the inland weak penetration zone with the contribution coefficient to form a hierarchical ozone transmission channel.

[0128] In several embodiments provided by the present application, it should be understood that the disclosed method and system can be implemented in other manners. For example, the division of the system embodiments described above is merely illustrative, and the division of the modules can be other division manners in actual implementation.

[0129] The modules described as separated components can or can not be physically separated, and the components displayed as modules can or can not be physical units, i.e., can be located in one place, or can be distributed on a plurality of network units. Part or all of the modules can be selected according to actual needs to achieve the purpose of the embodiments.

[0130] In addition, each functional module in each embodiment of the present application can be integrated in one processing unit, or each unit can be physically present separately, 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 functional modules.

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

[0132] The embodiments of the present application can acquire and process related data based on artificial intelligence technology. Among them, artificial intelligence is to use digital computers or machine controlled by digital computers to simulate, extend and expand human intelligence, perceive environment, acquire knowledge and use knowledge to obtain the best results.

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

Claims

1. A method for analyzing sea-land breeze ozone transmission channels, characterized in that: The method comprises: S1: Obtain meteorological elements and ozone concentration in the target detection area; S2: Identify the duration of the sea breeze in the meteorological element, calculate the offshore distance of the target detection area, and use the attenuation relationship between the duration of the sea breeze and the offshore distance as the sea breeze penetration intensity value of the target detection area; S3: extracting the ozone change slope from the accumulation period to the dissipation period of the ozone concentration, and dividing the target detection area into a coastal strong permeability area and an inland weak permeability area according to the sea breeze penetration intensity value; S4: determining the inversion layer height attenuation rate of the target detection area based on the vertical temperature gradient in the meteorological element; S5: superimposing the inversion layer height attenuation rate and the ozone change slope to obtain the ozone sinking contribution coefficient of the target detection area; S6: Spatially superimpose the coastal strong permeability area and the inland weak permeability area with the contribution coefficient to form a graded ozone transmission channel.

2. The method for analyzing the sea-land breeze ozone transmission channel according to claim 1, characterized in that: The obtaining of meteorological elements and ozone concentration in the target detection area includes: Continuously collecting meteorological observation raw data of the target detection area and synchronously obtaining regional ozone information of the target detection area; Abnormal data in the meteorological observation raw data and the regional ozone information are eliminated to obtain the meteorological elements and ozone concentration of the target detection area.

3. The method for analyzing the sea-land breeze ozone transmission channel according to claim 1, wherein: The identifying the duration of the sea breeze in the meteorological element and calculating the offshore distance of the target detection area include: Identifying a wind direction change event in a wind field time domain sequence in the meteorological element, and calibrating the duration of the wind direction change event as the duration of the sea breeze; The offshore distance of the target detection area is determined based on the coordinate reference positioning of the target detection area.

4. The method for analyzing the sea-land breeze ozone transmission channel according to claim 1, wherein: The method of using the attenuation relationship between the duration of the sea breeze and the offshore distance as the sea breeze penetration intensity value of the target detection area includes: Performing a time-weighted evaluation on the duration of the sea breeze based on a preset seasonal wind field scaling parameter to obtain a time-weighted intensity of the duration of the sea breeze; Analyzing the spatial attenuation intensity value of the offshore distance; The time-weighted intensity and the spatial attenuation intensity value are synthesized by penetration capacity to obtain the sea breeze penetration intensity value of the target detection area.

5. The method for analyzing the sea-land breeze ozone transmission channel according to claim 4, characterized in that: The performing penetration capacity synthesis on the time-weighted intensity and the spatial attenuation intensity value to obtain the sea breeze penetration intensity value of the target detection area includes: ， in: 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 attenuation intensity value, is the absolute value of the altitude difference, is the maximum terrain relief.

6. The method for analyzing the sea-land breeze ozone transmission channel according to claim 5, characterized in that: The method of dividing the target detection area into a coastal strong permeability area and an inland weak permeability area according to the sea breeze permeability intensity value includes: Using the sea breeze penetration intensity value as the penetration threshold of the target detection area; Performing partition marking on the geographic spatial distribution of the target detection area based on the penetration threshold; The spatial continuity of the results of the partition marking is verified, and the results that pass the verification are divided into coastal strong permeability areas and inland weak permeability areas.

7. The method for analyzing sea-land breeze ozone transmission channels according to claim 1, characterized in that: The determining of the inversion layer height attenuation rate of the target detection area based on the vertical temperature gradient in the meteorological element includes: Extracting the height boundary feature of the vertical temperature gradient, and simultaneously analyzing the height change difference of the height boundary feature in a continuous time series; The height change differences obtained through analysis are aggregated into the inversion layer height attenuation rate of the target detection area.

8. The method for analyzing sea-land breeze ozone transmission channels according to claim 6, characterized in that: The step of superimposing the inversion layer height attenuation rate and the ozone change slope as the ozone sinking contribution coefficient of the target detection area includes: Quantitatively analyzing the relationship between the vertical distribution of ozone and the distribution characteristics of the inversion layer height attenuation rate and the ozone change slope; Generating contribution weights for the ozone vertical distribution relationship to obtain an ozone vertical sinking value; The adaptation area corresponding to the ozone vertical sink value is matched to obtain the ozone sink contribution coefficient of the target detection area.

9. The method for analyzing sea-land breeze ozone transmission channels according to claim 8, characterized in that: The step of spatially superimposing the coastal strong permeability zone and the inland weak permeability zone with the contribution coefficient to form a hierarchical ozone transmission channel comprises: Constructing a composite permeability boundary map based on the spatial location data of the coastal high permeability area and the inland low permeability area; performing contribution coupling on the ozone sinking contribution coefficient and the composite infiltration boundary map to obtain an ozone distribution area; Dividing the ozone distribution area according to intensity to obtain intensity levels of the ozone distribution area; The ozone distribution area is divided into different graded ozone transmission channels based on the intensity levels.

10. A sea-land breeze ozone transmission channel analysis system, characterized in that: The system comprises: Data acquisition module: used to obtain meteorological elements and ozone concentration in the target detection area; Sea breeze impact module: used to identify the duration of the sea breeze in the meteorological element, calculate the offshore distance of the target detection area, and use the attenuation relationship between the duration of the sea breeze and the offshore distance as the sea breeze penetration intensity value of the target detection area; A region division module is used to extract the ozone change slope from the accumulation period to the dissipation period of the ozone concentration, and divide the target detection area into a coastal strong permeability area and an inland weak permeability area according to the sea breeze penetration intensity value; Temperature analysis module: used to determine the inversion layer height attenuation rate of the target detection area based on the vertical temperature gradient in the meteorological element; Ozone contribution module: used to superimpose the inversion layer height attenuation rate and the ozone change slope to obtain the ozone sinking contribution coefficient of the target detection area; Channel generation module: used for spatially superimposing the coastal strong permeability area and the inland weak permeability area with the contribution coefficient to form a graded ozone transmission channel.

Images