Air quality mode meteorological field generation method and system

Abstract

The invention belongs to the technical field of air quality numerical forecasting, and relates to an air quality mode meteorological field generation method and system, and the method comprises the steps: 1) dividing a large region into a plurality of sub-regions, and determining an optimal parameterization scheme combination; 2) setting a buffer area between adjacent sub-areas, and determining an optimal parameterization scheme combination; 3) obtaining a first meteorological field of an inner layer nested region of each sub-region and each buffer region; 4) fitting in the east-west direction to obtain a second meteorological field of the inner-layer nested region of each sub-region; 5) performing north-south fitting to obtain a final meteorological field of the inner-layer nested region of each sub-region; according to the characteristics of the meteorological field required by the air quality mode, a multi-area simulation result coupling technology is used, and the forecasting effect of the air quality mode is further improved by improving the forecasting accuracy of the meteorological mode.

Description

Technical Field

[0001] The present invention belongs to the technical field of air quality numerical forecasting, and relates to a meteorological field generation method and system, in particular to an air quality model meteorological field generation method and system. Background Art

[0002] With the development of numerical models, they have become an important approach for pollutant concentration prediction. In actual operational forecasting, third-generation air quality models are primarily used to construct air quality forecast systems. The mainstream third-generation air quality models today include the Community Multiscale Air Quality Modeling System (MODEL-3 / CAMQ), the Community Atmosphere Model (CAMx), the WRF-CHEM model, and the Nested Gridded Air Quality Prediction System (NAQPMS). These models incorporate complex and comprehensive gas-phase chemistry and photochemistry mechanisms, providing excellent simulation and prediction capabilities for the spatial and temporal distribution of pollutants.

[0003] Regardless of which third-generation air quality model is used to construct an air quality forecast system, a complete air quality forecast system primarily consists of three components: an emission source treatment system (providing emission source inputs), a meteorological model (providing meteorological fields such as temperature, pressure, humidity, and wind), and an air quality model (simulating the spatiotemporal distribution of pollutants). Therefore, meteorological fields are a crucial input to the air quality forecast system and significantly influence its accuracy. Therefore, the forecast accuracy of the meteorological model significantly influences the accuracy of air quality forecasts.

[0004] In weather forecasting, how to use weather models to obtain the most realistic meteorological conditions has always been a key research topic. Regional weather forecasting departments organize dedicated teams to adapt weather models to local conditions. Adjusting weather model parameters has always been a key aspect of local adaptation. Weather models incorporate a large number of physical parameterization schemes, which primarily describe various physical processes in the atmosphere through mathematical modeling. For example, the WRF weather model includes over ten parameterization schemes, including cumulus convection parameterization, boundary layer parameterization, microphysics parameterization, and land surface process parameterization. Each parameterization scheme offers multiple options. For example, there are over ten boundary layer parameterization schemes, including the MYJ boundary layer parameterization scheme, the MRF boundary layer parameterization scheme, the ACM2 boundary layer parameterization scheme, the QNSE boundary layer parameterization scheme, and the MYNN boundary layer parameterization scheme. This requires selecting the most suitable parameterization scheme from various physical parameterization schemes based on local actual conditions. That is, it is necessary to select one from various boundary layer schemes, one from various cumulus convection parameterization schemes, and one from various other schemes. These schemes are combined to form the optimal parameterization scheme combination, so that the meteorological model has the best forecast effect on the local area.

[0005] The local applicability of existing meteorological models is of great practical significance for weather forecasting. However, there are certain differences between the meteorological fields required by air quality models and those used by meteorological departments for forecasting. If the meteorological model for a specific region is adjusted to provide good local forecast results, it can be applied to the actual local weather forecast system. However, the meteorological field required by air quality models generally covers a relatively large area. This is mainly because the transmission of pollutants has a significant impact. For example, at the scale of a typical city, the local contribution of particulate matter is approximately 30%, while the transmission contribution is approximately 70%; the local contribution of ozone is approximately 20%, while the transmission contribution is approximately 80%. If the simulation and prediction area is small, it is difficult to accurately consider the transmission contribution of surrounding areas, and the simulation error of pollutant concentrations is generally large.

[0006] This requires selecting a combination of schemes that are applicable to a larger area during the parameter adjustment process. However, it is difficult to select a single parameterized scheme combination that is applicable to all areas, which is a significant factor in the low forecast accuracy. For example, a certain parameterized scheme combination may provide a good forecast for the meteorological field in Area A, but a poor forecast for the meteorological field in a neighboring area, Area B, which significantly contributes to the transmission of information to Area A. Using Area B's poorly forecasted meteorological field to drive the air quality model will result in large errors in the forecast of pollutant concentrations in Area B. Large errors in the forecast for Area B will lead to large errors in the transmission of information to Area A, resulting in large errors in the forecast of pollutant concentrations in Area A. Therefore, even when developing an air quality forecast system for Area A, the accuracy of the meteorological field forecasts in the surrounding areas must be considered.

[0007] Existing technologies for improving meteorological forecast accuracy primarily utilize methods such as data assimilation, ensemble forecasting, and improved parameterization schemes. Data assimilation primarily involves assimilating monitoring data into the initial field to improve the accuracy of the initial field, thereby enhancing the accuracy of meteorological forecasts. Ensemble forecasting primarily accounts for the inherent errors in both the initial field and each parameterization scheme. Therefore, during simulation and forecasting, multiple initial value ensemble forecasts are constructed based on the initial field. Physical process ensemble members are constructed based on different parameterization schemes within the same class. Each ensemble member drives the meteorological model, resulting in a series of forecast results that are averaged using the ensemble mean. Improved parameterization schemes primarily refine the physical representation of the parameterization scheme, making it more similar to the physical processes in the actual atmosphere and thus applicable over a wider area. These technologies have, to a certain extent, improved meteorological forecast accuracy, and thus the accuracy of air quality forecasts. However, these technologies still cannot address the cross-regional disparities in air quality forecasting.

[0008] The Chinese invention patent application with application number 202211577968.1 discloses a method for determining weather forecast values ​​across multiple climate zones. The method mainly divides the target simulation area into multiple sub-areas according to the projection band, and then interpolates the simulation forecast values ​​with the minimum nested resolution of each sub-area to the longitude and latitude to generate the weather forecast values ​​for each longitude and latitude of the target area. Because the two sub-areas are simulated and predicted separately, there will be differences in the simulation results of the two sub-areas for each longitude and latitude in the overlapping area of ​​the two sub-areas. For this reason, the method performs weighted averaging. However, this method is mainly used in the field of weather forecasting. As long as the relative error of the meteorological forecast value at each longitude and latitude is small, there is no need to consider the abnormal discontinuity problem of the meteorological field. For example, suppose the actual wind in the overlapping area and surrounding areas during a certain period is westerly (270 degrees 2 m / s). The wind forecast for the first sub-area (denoted as Sub-area 1) is west-northwest (280 degrees 2 m / s), and the wind forecast for the second sub-area (denoted as Sub-area 2) is west-southwest (260 degrees 2 m / s). If the forecast is performed according to the weather forecast requirements, both sub-areas have good forecast results. If the weighted average of the overlapping area is used, the forecast error for the overlapping area will be smaller. However, the wind changes from Sub-area 1 to the overlapping area and then to Sub-area 2, from northwest to westerly and then to southwest (forming a wind shear, which is an abnormal discontinuity). Although the polluted air mass still generally spreads from west to east, a weak convergence field forms in the boundary area, causing pollutants to be trapped there. Therefore, this meteorological field is not suitable for air quality model simulation.

[0009] At the same time, the meteorological field used for weather forecasting only needs to splice the results of the inner nested areas of the meteorological model, and the simulation results of the outer nested areas only provide initial values ​​and boundary values ​​for the inner nested areas, and are not used for specific analysis. Therefore, the method in the Chinese invention patent application with application number 202211577968.1 does not need to process the outer nested areas. When driving the air quality model, not only the meteorological background field of the inner nested area is required, but also the meteorological field of the outer nested area. The meteorological fields of each nested area drive the air quality model of each nested area respectively. The simulation results of the air quality model of the outer nested area provide initial value and boundary value conditions for the inner nested area. If the prediction error of the outer nested area for a certain area in the air quality model is large, it will affect the prediction results of the inner nested area. If the method in the above-mentioned invention patent application is used to splice the meteorological field of the outer nested area, not only will the meteorological field be discontinuous at the boundary area, but the inner nested area and the outer nested area will be spliced ​​separately, which will cause the meteorological field mismatch between the outer nested area and the inner nested area (in the original simulation results of the meteorological model, the outer nested and the inner nested directly feedback each other and there will be no mismatch).

[0010] Therefore, the technology disclosed in Chinese invention patent application number 202211577968.1 struggles to meet the requirement for simultaneously achieving good simulation and prediction results for each sub-region within a large region when forecasting regional meteorological fields, and thus struggles to meet the needs of air quality model forecasting. Furthermore, simple splicing not only creates discontinuities in the meteorological field at the splicing points, but also leads to mismatches between the meteorological fields of the nested regions.

[0011] Therefore, in view of the defects existing in the above-mentioned existing technologies, there is an urgent need for a new air quality model meteorological field generation method and system to improve the accuracy of meteorological field prediction in each region, thereby improving the air quality model prediction effect. Summary of the Invention

[0012] In response to the shortcomings of the existing technology, the present invention proposes a method and system for generating a meteorological field for an air quality model. On the basis of fully considering the characteristics of the meteorological field required by the air quality model, a buffer area is set to perform simulation and prediction in different regions, and the meteorological field generated in each region is fitted to obtain the meteorological field of a large region to drive the air quality model, thereby improving the accuracy of the air quality model forecast.

[0013] In order to achieve the above object, the present invention provides the following technical solutions: A method for generating an air quality model meteorological field, characterized by comprising the following steps: 1) Divide a large area into multiple sub-areas and determine the optimal parameterization scheme combination for each sub-area based on the simulation results of historical meteorological field data over many years; 2) Set up buffer zones between two adjacent sub-regions and determine the optimal parameterization scheme combination for each buffer zone based on simulation results of historical meteorological data over many years; 3) Use the optimal parameterization scheme combination of each sub-region and each buffer area to predict future weather and obtain the first meteorological field of the inner nested area of ​​each sub-region and each buffer area; 4) Combine the first meteorological field of the inner nested area of ​​the buffer area between each sub-area divided in the east-west direction and the first meteorological field of the inner nested area of ​​each sub-area to perform east-west fitting to obtain the second meteorological field of the inner nested area of ​​each sub-area; 5) Combine the first meteorological field of the inner nested area of ​​the buffer area between each sub-area divided in the north-south direction and the second meteorological field of the inner nested area of ​​each sub-area to perform north-south fitting to obtain the final meteorological field of the inner nested area of ​​each sub-area; 6) splicing the final meteorological fields of the inner nested areas of each sub-area to obtain the meteorological field of the large area; 7) intercepting the large area according to the geographical range of each nested area of ​​the target area under the air quality model, and obtaining the meteorological field of each nested area of ​​the target area under the air quality model based on the meteorological field of the intercepted large area.

[0014] Preferably, in step 1), determining the optimal parameterization scheme combination for each sub-region based on the simulation effect of historical meteorological field data over many years specifically includes: 11) Based on a sub-region, different parameterization scheme combinations are used to simulate the historical meteorological field of the sub-region over many years. The simulation results are tested using monitoring data to select the optimal parameterization scheme combination for the sub-region; 12) Using different parameterization scheme combinations to simulate the historical meteorological fields of the adjacent sub-regions of the sub-region over many years, using monitoring data to verify the simulation results, screening out the better parameterization scheme combinations of the adjacent sub-regions of the sub-region, and referring to the optimal parameterization scheme combination of the sub-region, determining the optimal parameterization scheme combination of the adjacent sub-regions of the sub-region based on the principle of minimizing the adjustment relative to the optimal parameterization scheme combination of the sub-region; 13) Use different parameterization scheme combinations to simulate the historical meteorological fields of other adjacent sub-regions of each adjacent sub-region of the sub-region, use monitoring data to test the simulation effect, screen out the better parameterization scheme combinations of other adjacent sub-regions of each adjacent sub-region of the sub-region, and refer to the optimal parameterization scheme combination of each adjacent sub-region of the sub-region, so as to minimize the adjustment relative to the optimal parameterization scheme combination of each adjacent sub-region of the sub-region, and determine the optimal parameterization scheme combination of other adjacent sub-regions of each adjacent sub-region of the sub-region.

[0015] Preferably, in step 2), determining the optimal parameterization scheme combination for each buffer area based on the simulation effect of the historical meteorological field data for many years specifically includes: simulating the historical meteorological field of the buffer area for many years using different parameterization scheme combinations, using monitoring data to test the simulation effect, screening out the better parameterization scheme combination for the buffer area, and referring to the optimal parameterization scheme combination of two adjacent sub-areas of the buffer area, so as to determine the optimal parameterization scheme combination for the buffer area based on the principle of minimizing the adjustment relative to the optimal parameterization scheme combination of the two adjacent sub-areas of the buffer area.

[0016] Preferably, the step 4) specifically includes: 41) Divide the buffer area between each sub-area divided in the east-west direction into a central area and two edge areas, wherein the first meteorological field of the inner nested area of ​​the buffer area is used as the second meteorological field of the inner nested area of ​​the central area; 42) using the fitting results of the first meteorological field of the inner nested area of ​​the buffer area and the first meteorological field of the inner nested area of ​​the sub-area that overlaps with the edge area as the second meteorological field of the inner nested area of ​​the edge area; 43) using the first meteorological field of the inner nested area of ​​each sub-area as the second meteorological field of the inner nested area of ​​the portion of each sub-area that does not overlap with the buffer area; 44) Integrate the second meteorological field of the inner nested area of ​​the central area, the second meteorological field of the inner nested area of ​​the edge area, and the second meteorological field of the inner nested area of ​​each sub-area that does not overlap with the buffer area to obtain the second meteorological field of the inner nested area of ​​each sub-area.

[0017] Preferably, in step 42), the fitting results of the first meteorological field of the inner nested area of ​​the buffer area and the first meteorological field of the inner nested area of ​​the sub-area overlapping with the edge area are specifically: Where, is the meteorological field of the i-th grid in a row of grids in the east-west direction of the inner nested area of ​​the edge area, starting from the central area, and n is the total number of grids in the row of grids in the east-west direction of the inner nested area of ​​the edge area. is the meteorological field of the corresponding grid of the inner nested area of ​​the buffer area, It is the meteorological field of the corresponding grid of the inner nested area of ​​the sub-area that overlaps with the edge area.

[0018] Preferably, the step 7) specifically includes: 71) intercepting the large area according to the geographical scope of the inner nested area of ​​the target area under the air quality model, and using the meteorological field of the intercepted large area as the meteorological field of the inner nested area of ​​the target area under the air quality model; 72) The large area is intercepted based on the geographical scope of other nested areas of the target area under the air quality model, and the meteorological fields of other nested areas of the target area under the air quality model are processed based on the meteorological field of the intercepted large area.

[0019] Preferably, in step 72), based on the intercepted meteorological field of the large area, obtaining meteorological fields of other nested areas of the target area in the air quality mode specifically includes: 721) The meridional wind in the meteorological field of each grid in the other nested regions is obtained by taking a weighted average of the meridional winds in the meteorological field of the westernmost column of grids in the large region contained in each grid in the other nested regions; 722) The latitudinal wind in the meteorological field of each grid in the other nested region is obtained by taking a weighted average of the latitudinal winds in the meteorological field of the northernmost column of grids in the large region contained in each grid in the other nested region; 723) The prediction results of the remaining meteorological elements in the meteorological fields of the grids in the other nested regions are obtained by taking the weighted average of the prediction results of the corresponding meteorological elements in the meteorological fields of all the grids in the large region contained in the grids in the other nested regions.

[0020] In addition, the present invention also provides an air quality model meteorological field generation system, which is characterized by comprising: The module for determining the optimal parameterization scheme combination for a sub-region is used to divide a large region into multiple sub-regions and determine the optimal parameterization scheme combination for each sub-region based on the simulation results of meteorological field data over many years of history; A module for determining the optimal parameterization scheme combination for buffer areas, which is used to set up buffer areas between two adjacent sub-areas and determine the optimal parameterization scheme combination for each buffer area based on the simulation results of historical meteorological field data over many years; A first meteorological field determination module is used to predict future weather using an optimal parameterization scheme combination of each sub-region and each buffer region to obtain a first meteorological field for an inner nested region of each sub-region and each buffer region; An east-west fitting module is used to perform east-west fitting on the first meteorological field of the inner nested area of ​​the buffer area between each sub-area divided in the east-west direction and the first meteorological field of the inner nested area of ​​each sub-area to obtain the second meteorological field of the inner nested area of ​​each sub-area; A north-south fitting module is used to combine the first meteorological field of the inner nested area of ​​the buffer area between each sub-area divided in the north-south direction and the second meteorological field of the inner nested area of ​​each sub-area to perform north-south fitting to obtain the final meteorological field of the inner nested area of ​​each sub-area; A large-area meteorological field determination module, configured to combine the obtained meteorological fields of the inner nested areas of each sub-area to obtain the meteorological field of the large area; The meteorological field determination module for each nested area is used to intercept the large area according to the geographical scope of each nested area of ​​the target area under the air quality mode, and obtain the meteorological field of each nested area of ​​the target area under the air quality mode based on the meteorological field of the intercepted large area.

[0021] Furthermore, the present invention also provides an air quality model meteorological field generating device, characterized in that it includes: one or more processors; a memory for storing one or more programs; When the one or more programs are executed by the one or more processors, the one or more processors are enabled to implement the air quality model meteorological field generation method as described above. Finally, the present invention also provides a computer-readable storage medium having a computer program stored thereon, characterized in that when the program is executed by a processor, the steps of the air quality model meteorological field generation method as described above are implemented.

[0022] Compared with the prior art, the air quality model meteorological field generation method and system of the present invention has one or more of the following beneficial technical effects: 1. The present invention improves the accuracy of the simulation prediction effect of each sub-region by combining different parameterization schemes in different sub-regions.

[0023] 2. The present invention reduces the impact of abnormal discontinuity of meteorological fields in adjacent sub-regions caused by the use of different parameterization scheme combinations on air quality model simulation by setting buffer areas and performing east-west and north-south fitting.

[0024] 3. The present invention obtains the meteorological field of the outer nested area by processing the meteorological field of the inner nested area, thereby solving the problem of mismatch between the meteorological fields of the inner and outer nested areas.

[0025] 4. The present invention ultimately obtains a meteorological field for driving the air quality model, thereby improving the accuracy of air quality model prediction. BRIEF DESCRIPTION OF THE DRAWINGS

[0026] Figure 1 It is a flow chart of the method for generating an air quality model meteorological field of the present invention.

[0027] Figure 2 A schematic diagram showing the large-area division of the present invention.

[0028] Figure 3 A nested schematic diagram of sub-regions of the present invention is shown.

[0029] Figure 4 A schematic diagram showing the buffer area arrangement of the present invention is shown.

[0030] Figure 5 A schematic diagram showing the east-west fitting of the present invention is shown.

[0031] Figure 6 A schematic diagram of determining the meridional wind of each grid in the outer nested area according to the present invention is shown.

[0032] Figure 7 A schematic diagram of determining the zonal wind of each grid in the outer nested area according to the present invention is shown.

[0033] Figure 8A schematic diagram showing the prediction results of other meteorological elements of each grid in the outer nested area according to the present invention is shown.

[0034] Figure 9 It is a schematic diagram of the composition of the air quality model meteorological field generation system of the present invention.

[0035] in, Figure 6-8 In the figure, the left side is a schematic diagram of the grid in the large area, and the right side is a schematic diagram of the grid of the outer nested area of ​​the air quality model. One grid in the outer nested area includes 9 grids in the large area. DETAILED DESCRIPTION

[0036] Before describing in detail any embodiment of the present invention, it should be understood that the present invention is not limited in its application to the construction and arrangement details of the components set forth in the following description or illustrated in the following figures. The present invention is capable of other embodiments and can be practiced or carried out in various ways. In addition, it should be understood that the words and terms used herein are for descriptive purposes and should not be considered restrictive. The use of "including" or "having" and their variations herein is intended to cover the items and their equivalents set forth below and additional items. Unless otherwise specified or limited, the terms "mounted", "connected", "supported" and "coupled" and their variations are used broadly and cover direct mounting and indirect mounting, connection, support and coupling. In addition, "connected" and "coupled" are not limited to physical or mechanical connections or couplings. Furthermore, on the first hand, in the disclosure of the present invention, the terms "longitudinal", "transverse", "up", "down", "front", "back", "left", "right", "vertical", "horizontal", "top", "bottom", "inside", "outside" and the like to indicate orientations or positional relationships are based on the orientations or positional relationships shown in the accompanying drawings, which are only for the convenience of describing the present invention and simplifying the description, rather than indicating or implying that the device or element referred to must have a specific orientation, be constructed and operated in a specific orientation, and therefore the above terms cannot be understood as limitations on the present invention; on the second hand, the term "one" should be understood as "at least one" or "one or more", that is, in one embodiment, the number of an element may be one, while in another embodiment, the number of the element may be multiple, and the term "one" cannot be understood as a limitation on the quantity.

[0037] In response to the defects of the prior art, the present invention provides a method and system for generating a meteorological field for an air quality model. Sub-regions are set and buffer zones are set between the sub-regions based on full consideration of the characteristics of the meteorological field required for the air quality model. Simulation and prediction are performed in different regions based on a combination of parameterized schemes, that is, each sub-region adopts a different combination of parameterized schemes for prediction. On the basis of considering that the atmosphere is a fluid, based on the principle of continuity of the meteorological field, the meteorological field generated by each sub-region is quadratically fitted, and a high-resolution meteorological field of a large area is obtained by splicing. This high-resolution meteorological background field is then combined with the area simulated by the air quality model to obtain the meteorological background fields of each nested layer for driving the air quality model, thereby driving the air quality model to improve the accuracy of the air quality model forecast. Figure 1 FIG. 1 is a flow chart showing the method for generating an air quality model meteorological field according to the present invention. Figure 1 As shown, the air quality model meteorological field generation method of the present invention includes the following steps: 1. Determine the optimal parameterization scheme combination for the sub-region.

[0038] The large area is divided into multiple sub-areas, and the optimal parameterization scheme combination for each sub-area is determined based on the simulation results of historical meteorological field data over many years.

[0039] In this invention, a specific implementation method is described by dividing a large area (e.g., a large district) into four sub-areas, with each sub-area implementing double nesting (the outer nested area is gridded with a grid of 9 km, i.e., a square grid with a side length of 9 km, also known as a grid with a 9 km resolution; the inner nested area is also gridded with a grid of 3 km, i.e., a square grid with a side length of 3 km, also known as a grid with a 3 km resolution; thus, one grid in the outer nested area includes nine grids in the inner nested area). In practical applications, multiple sub-areas (e.g., 8 or 16) and multiple nesting (e.g., triple nesting, with a grid of 27 km in the outer nested area, a grid of 9 km in the middle nested area, and a grid of 3 km in the inner nested area) can be set according to actual needs. The specific processing method is the same as that for four sub-areas and double nesting.

[0040] For ease of expression, Figure 2 As shown in the figure, the present invention divides the large area into four sub-areas (northeast, northwest, southeast and southwest, respectively denoted as NE, NW, SE and SW). Figure 3 As shown, each sub-area includes an inner nested area and an outer nested area. In addition, when setting, ensure that the innermost nested areas of each sub-area are closely connected and do not overlap (such as Figure 2), the outer nested areas can overlap. Because the outer nested area of ​​the present invention only provides initial values ​​and boundary values ​​to the inner nested area, the data of the outer nested area is not processed subsequently. For the convenience of display, Figure 2 The outer nested area is hidden and only the inner nested area is displayed.

[0041] In the present invention, the optimal parameterization scheme combination for each sub-region is determined based on the simulation effect of historical meteorological field data over many years, specifically including: 1. Based on a certain sub-region, different parameterization scheme combinations are used to simulate the historical meteorological field of the sub-region for many years. The simulation effect is tested using monitoring data to screen out the optimal parameterization scheme combination for the sub-region.

[0042] For example, based on the Northeast sub-region, the historical three-year meteorological field of the Northeast sub-region is simulated using WRF or other meteorological models with different parameterization scheme combinations. The simulation results are tested using monitoring data to screen out the optimal parameterization scheme combination for the Northeast sub-region.

[0043] Among them, when using the WRF meteorological model, the optimal parameterization scheme combination includes more than 10 types of parameterization schemes, such as cumulus convection parameterization scheme, boundary layer parameterization scheme, microphysical parameterization scheme, and land surface process parameterization scheme. Each parameterization scheme includes the optimal parameterization scheme adopted.

[0044] 2. Use different parameterization scheme combinations to simulate the historical meteorological fields of adjacent sub-regions within the sub-region. Use monitoring data to verify the simulation results, screen out the optimal parameterization scheme combinations for adjacent sub-regions within the sub-region, and refer to the optimal parameterization scheme combination of the sub-region. Based on the principle of minimizing the adjustment relative to the optimal parameterization scheme combination of the sub-region, determine the optimal parameterization scheme combination for adjacent sub-regions within the sub-region. This process ensures that the optimal parameterization scheme combinations of each sub-region are relatively small, thereby reducing abnormal discontinuities in the meteorological fields of adjacent sub-regions.

[0045] For example, use WRF or other meteorological models with different parameterization scheme combinations to simulate the historical three-year meteorological fields of the adjacent sub-regions of the northeast sub-region, that is, the northwest and southeast sub-regions, and use monitoring data to test the simulation results to screen out the better parameterization scheme combinations for the northwest and southeast sub-regions.

[0046] At the same time, referring to the optimal parameterization scheme combination in the Northeast subregion, the optimal parameterization scheme combination for the Northwest and Southeast subregions was determined, with the principle of minimizing adjustments relative to the optimal parameterization scheme combination in the Northeast subregion. For example, if the boundary layer parameterization scheme in the Northeast subregion is the MYJ boundary layer parameterization scheme, and the preferred boundary layer parameterization schemes selected in the Northwest subregion are comparable to the MYJ boundary layer parameterization scheme and the MRF boundary layer parameterization scheme, the final boundary layer parameterization scheme for the Northwest subregion should be the MYJ boundary layer parameterization scheme to maintain consistency with the boundary layer parameterization scheme in the Northeast subregion as much as possible.

[0047] 3. Use different parameterization scheme combinations to simulate the historical meteorological fields of other adjacent sub-regions of each adjacent sub-region of the sub-region for many years, use monitoring data to test the simulation effect, screen out the better parameterization scheme combinations of other adjacent sub-regions of each adjacent sub-region of the sub-region, and refer to the optimal parameterization scheme combination of each adjacent sub-region of the sub-region to minimize the adjustment relative to the optimal parameterization scheme combination of each adjacent sub-region of the sub-region, and determine the optimal parameterization scheme combination of other adjacent sub-regions of each adjacent sub-region of the sub-region.

[0048] For example, use WRF or other meteorological models with different parameterization scheme combinations to simulate the historical three-year meteorological fields of the adjacent sub-regions of the northwest and southeast sub-regions except the northeast sub-region, that is, the southwest sub-region. Use monitoring data to test the simulation results and screen out the better parameterization scheme combination for the southwest sub-region.

[0049] At the same time, the optimal parameterization scheme combination for the southwest subregion is determined with reference to the optimal parameterization scheme combination for the northwest and southeast subregions, with the principle of minimizing adjustments relative to the optimal parameterization scheme combination for the northwest and southeast subregions. For example, if the boundary layer parameterization scheme for the northwest and southeast subregions is the MYJ boundary layer parameterization scheme, and the boundary layer parameterization scheme selected for the southwest subregion has comparable performance to the ACM2 boundary layer parameterization scheme, the MYJ boundary layer parameterization scheme should be selected as the final boundary layer parameterization scheme for the southwest subregion to maintain consistency with the boundary layer parameterization schemes for the northwest and southeast subregions as much as possible.

[0050] 2. Determine the optimal parameterization scheme combination for the buffer area.

[0051] A buffer area is set between two adjacent sub-areas, and the optimal parameterization scheme combination for each buffer area is determined based on the simulation results of historical meteorological field data over many years. Specifically, different parameterization scheme combinations are used to simulate the historical meteorological field of the buffer area over many years. The simulation results are verified using monitoring data, and the optimal parameterization scheme combination for the buffer area is screened out. The optimal parameterization scheme combination of the two adjacent sub-areas where the buffer area is set is then referenced, and the optimal parameterization scheme combination of the buffer area is determined based on the principle of minimizing the adjustment relative to the optimal parameterization scheme combination of the two adjacent sub-areas where the buffer area is set.

[0052] In an example of the present invention, it specifically includes: 1. If Figure 4 As shown, a buffer area 1 is set between the northeast and northwest sub-areas (buffer area 1 is also set to double nesting), which is recorded as HC1. Among them, the inner nested area of ​​buffer area 1 needs to completely overlap with the grid of the inner nested area of ​​the northeast and northwest sub-areas (such as Figure 4 , because the outer nested area of ​​the present invention only provides initial values ​​and boundary values ​​to the inner nested area, the data of the outer nested area will not be processed subsequently. For the convenience of display, Figure 4 The outer nested regions are hidden, showing only the inner nested regions. Three years of historical meteorological fields in buffer region 1 were simulated using WRF or other meteorological models with different parameterization scheme combinations. Monitoring data were used to verify the simulation results and identify the optimal parameterization scheme combination for buffer region 1. Furthermore, the optimal parameterization scheme combination for buffer region 1 was determined by minimizing adjustments to the optimal parameterization scheme combinations for the northeast and northwest subregions.

[0053] 2. Similarly, set up Buffer Area 2 between the southeast and southwest sub-areas, Buffer Area 3 between the northeast and southeast sub-areas, and Buffer Area 4 between the southwest and northwest sub-areas (Buffer Areas 2, 3, and 4 are also doubly nested). Buffer Areas 1 and 2 are buffers between sub-areas divided in the east-west direction, while Buffer Areas 3 and 4 are buffers between sub-areas divided in the north-south direction. Similarly to step 1 above, select the optimal parameterization scheme combination for Buffer Areas 2, 3, and 4.

[0054] In this paper, the purpose of setting up a buffer zone is to establish a transition zone between two connected sub-regions, allowing meteorological elements to evolve slowly and avoiding abnormal discontinuities in the meteorological field at the junction of the two adjacent sub-regions. For example, if the northeastern and northwest sub-regions are directly connected, and the two sub-regions use different parameterization schemes, abnormal discontinuities in the meteorological field at the junction will exist. This abnormal discontinuity may form a false convergence zone at the junction. Using this meteorological field to drive the air quality model will increase simulation errors.

[0055] 3. Determine the first meteorological field.

[0056] With the optimal parameterization scheme combination of each sub-region and each buffer area, the optimal parameterization scheme combination of each sub-region and each buffer area can be used to predict future weather and obtain the first meteorological field of each sub-region and each buffer area.

[0057] Because each sub-region and each buffer region is multiply nested (double nested in the present invention), future weather forecasting can yield both the first meteorological field for the innermost nested region of each sub-region and each buffer region, and the first meteorological field for the outermost nested region of each sub-region and each buffer region. However, the present invention only uses the first meteorological field for the innermost nested region of each sub-region and each buffer region in subsequent processing.

[0058] 4. East-west fitting.

[0059] Combine the first meteorological field of the inner nested area of ​​the buffer area between each sub-area divided in the east-west direction and the first meteorological field of the inner nested area of ​​each sub-area for east-west fitting to obtain the second meteorological field of the inner nested area of ​​each sub-area. The specific implementation steps are as follows: 1. Divide the buffer area between each sub-area divided east-west into a central area and two edge areas. The central area and edge areas are also doubly nested. The first meteorological field of the inner nested area of ​​the buffer area is used as the second meteorological field of the inner nested area of ​​the central area.

[0060] For example, Figure 5 As shown, the buffer zone 1 is divided into a central zone C and two edge zones D and E. The central zone C and the edge zones D and E are also doubly nested. The widths of the central zone D and the two edge zones D and E can each account for one-third of the width of the buffer zone 1. The first meteorological field of the inner nested zone of the buffer zone 1 is used as the second meteorological field of the inner nested zone of the central zone C of the buffer zone 1.

[0061] 2. The fitting results of the first meteorological field of the inner nested area of ​​the buffer area and the first meteorological field of the inner nested area of ​​the sub-area overlapping with the edge area are used as the second meteorological field of the inner nested area of ​​the edge area.

[0062] For example, Figure 5 As shown, the area outside the middle area C of buffer area 1, that is, area D overlapping with the northeastern sub-area and area E overlapping with the northwest sub-area need to be fitted to solve the abnormal discontinuity problem of the meteorological field in the east-west direction in adjacent areas and improve the accuracy of the meteorological field.

[0063] For each grid within the inner nested area of ​​region D, there are two meteorological element prediction results (i.e., meteorological fields): the meteorological element prediction results for the inner nested area of ​​the northeastern sub-region and the meteorological element prediction results for the inner nested area of ​​buffer region 1. When fitting region D, to ensure the continuity of the meteorological field, the prediction results of the meteorological elements (temperature, air pressure, humidity, radial wind, zonal wind, and more than 100 meteorological elements) in the inner nested area of ​​buffer region 1 should be given greater weights for grids in the inner nested area of ​​region D close to region C. For grids in the inner nested area of ​​region B, the prediction results of the meteorological elements in the inner nested area of ​​the northeastern sub-region should be given greater weights. From C to B, the weights of the meteorological element prediction results in the inner nested area of ​​buffer region 1 gradually decrease, while the weights of the meteorological element prediction results in the inner nested area of ​​the northeastern sub-region gradually increase. The fitting formula is as follows: Where, is the forecast result of the meteorological elements of the i-th grid from the central area in a certain row of grids in the east-west direction of the inner nested area of ​​the edge area D, n is the total number of grids in the row of grids in the east-west direction of the inner nested area of ​​the edge area D, is the forecast result of the meteorological elements of the corresponding grid of the inner nested area of ​​the buffer area, It is the forecast result of the meteorological elements of the corresponding grid of the inner nested area of ​​the sub-area that overlaps with the edge area, and i is 1, 2, 3...n from C to B.

[0064] By using the above fitting formula to fit all grids in a certain row of grids in the east-west direction of the inner nested area of ​​the edge zone, the forecast results of the meteorological elements of all grids in the row of grids can be obtained.

[0065] It can be seen from the formula that the first grid in the inner nested area near area C on the left side of area D uses the prediction results of the meteorological elements of the inner nested area of ​​buffer area 1, and the first grid in the inner nested area near area B on the right side uses the prediction results of the meteorological elements of the inner nested area of ​​the northeast sub-area. The weights of the prediction results of the inner nested areas of the buffer areas from C to B gradually decrease, and the weights of the prediction results of the meteorological elements of the inner nested areas of the northeast sub-area gradually increase.

[0066] The above fitting formula is used to fit all grids in each row of grids in the east-west direction of the inner nested area of ​​the edge area, so as to obtain the forecast results of the meteorological elements of the inner nested area of ​​the edge area after fitting.

[0067] In the same way, the forecast results of meteorological elements in the inner nested areas of area E can be obtained.

[0068] 3. Using the first meteorological field of the inner nested area of ​​each sub-area as the second meteorological field of the inner nested area of ​​the portion of each sub-area that does not overlap with the buffer area.

[0069] For example, outside buffer area 1, area B of the northeastern sub-area and area A of the northwest sub-area use the forecast results of meteorological elements of the inner nested areas of their respective sub-areas. 4. Integrate the second meteorological field of the inner nested area of ​​the central area, the second meteorological field of the inner nested area of ​​the edge area, and the second meteorological field of the inner nested area of ​​each sub-area that does not overlap with the buffer area to obtain the second meteorological field of the inner nested area of ​​each sub-area.

[0070] That is, the prediction results of the meteorological elements of the inner nested areas of area C obtained in steps 2-3, the prediction results of the meteorological elements of the inner nested areas of area D, the prediction results of the meteorological elements of the inner nested areas of area E, the prediction results of the meteorological elements of the inner nested areas of area A, and the prediction results of the meteorological elements of the inner nested areas of area B are integrated to obtain the prediction results of the meteorological elements of the inner nested areas of the new northeast and northwest sub-areas, that is, the second meteorological field.

[0071] Similarly, by fitting the buffer area 2 between the southwest sub-region and the southeast sub-region using the above steps 1-4, the forecast results of the meteorological elements of the inner nested areas of the new southwest and southeast sub-regions, that is, the second meteorological field, can be obtained.

[0072] However, the second meteorological field is only fitted in the east-west direction for the inner nested areas of the buffer areas of the northeastern sub-region and the northwest sub-region and the buffer areas of the southeastern sub-region and the southwest sub-region. The prediction results of meteorological elements in the north-south direction of the inner nested areas of the northeastern sub-region and the southeastern sub-region, and the northwest sub-region and the southwest sub-region may still be quite different, which will cause discontinuity of the meteorological field and require north-south fitting.

[0073] 5. North-South fitting.

[0074] The first meteorological field of the inner nested area of ​​the buffer area between the sub-areas divided in the north-south direction and the second meteorological field of the inner nested area of ​​each sub-area are combined to perform north-south fitting to obtain the final meteorological field of the inner nested area of ​​each sub-area, that is, the prediction results of the meteorological elements of the inner nested area of ​​each sub-area.

[0075] In the present invention, similar to the east-west fitting method, north-south fitting is performed based on the second meteorological field of the inner nested area of ​​the buffer area 4 and the northwest sub-area and the southwest sub-area, and north-south fitting is performed based on the second meteorological field of the inner nested area of ​​the buffer area 3 and the northeast sub-area and the southeast sub-area, so as to obtain the final meteorological field of the inner nested area of ​​each sub-area, that is, the prediction results of the meteorological elements in the inner nested area of ​​each sub-area.

[0076] 6. Determine the meteorological field of a large area.

[0077] After obtaining the final meteorological fields of the inner nested regions of each sub-region, the meteorological fields of the inner nested regions of each sub-region are spliced ​​to obtain the meteorological field of the large region.

[0078] In the present invention, since the previous processing is based on the inner nested area, it can be seen that the meteorological field of the large area is only the meteorological field of the inner nested area, that is, the meteorological field of the 3 km resolution grid.

[0079] 7. Determine the meteorological field of each nested area.

[0080] When performing an air quality model simulation, you first need to determine the area for the air quality model simulation, that is, the target area in the air quality model. For example, if you are simulating an air quality model for a certain location, the target area is the geographical area within a certain distance of the location and its surroundings.

[0081] At the same time, air quality model simulations are generally multi-nested, for example, double nesting. In this case, the outer nested area (e.g., a certain place) is larger and the grid of the outer nested area is larger, for example, the grid of the outer nested area is 9 km, but the resolution is lower. The inner nested area (e.g., a certain area within the said place) is smaller than the outer nested area and the grid of the inner nested area is smaller than the grid of the outer nested area, for example, the grid of the inner nested area is 3 km, but the resolution is higher. Of course, it can also be triple nested or other heavily nested areas. In triple nesting, the grid of the outer nested area is 27 km, the grid of the middle nested area is 9 km, and the grid of the inner nested area is 3 km.

[0082] The present invention, through the processing of steps 1 to 6, can obtain a 3km-resolution grid meteorological field for a large area (e.g., a certain region). Thus, when performing an air quality model simulation, a target area in the air quality model (e.g., a certain location in Hebei Province or a certain city within the certain region) can directly obtain a 3km-resolution grid meteorological field within its own geographical scope based on the meteorological field of the large region. This eliminates the need for each target area to separately obtain a 3km-resolution grid meteorological field within its own scope through steps 1 to 6, thus saving significant costs.

[0083] In the present invention, when performing air quality model simulation, the large area can be intercepted according to the geographical range of each nested area of ​​the target area under the air quality model, and the meteorological field of each nested area of ​​the target area under the air quality model is obtained based on the meteorological field of the intercepted large area, which specifically includes: 1. Directly intercept the large area based on the geographical scope of the inner nested area of ​​the target area under the air quality mode, and use the meteorological field of the intercepted large area as the meteorological field of the inner nested area of ​​the target area under the air quality mode.

[0084] For example, the large area is directly intercepted according to the size of the inner nested area of ​​the target area under the air quality mode, and the meteorological field of each grid in the intercepted large area (that is, each 3km resolution grid of the large area) is used as the meteorological field of each grid in the inner nested area of ​​the target area under the air quality mode (which is also a 3km resolution grid).

[0085] 2. Based on other nested areas of the air quality model, for example, the geographical scope of the outer nested area intercepts the large area, and based on the meteorological fields of each grid in the large area contained in each grid of the other nested areas, the meteorological fields of each grid in other layers of nested areas of the air quality model are processed.

[0086] For example, in the present invention, the grid of the outer nested area of ​​the target area under the air quality mode is a grid with a resolution of 9 km, and the grid in the large area is a grid with a resolution of 3 km. Therefore, one grid of the outer nested area of ​​the target area under the air quality mode includes 9 grids in the large area. By processing the meteorological fields of the 9 grids in the large area, the meteorological field of the grid of the outer nested area of ​​the target area under the air quality mode can be obtained.

[0087] In the present invention, based on the meteorological fields of the grids in the large area contained in the grids of the other nested areas, processing to obtain the meteorological fields of the grids in the other nested areas of the target area in the air quality mode specifically includes: 1) The meridional wind in the meteorological field of each grid in the other nested regions is obtained by taking a weighted average of the meridional winds in the meteorological field of the westernmost column of grids in the large region contained in each grid in the other nested regions.

[0088] For example, Figure 6 As shown, the meridional wind u in a grid 1 of the outer nested area is generated by the westernmost grid of the nine grids in the large area contained by the grid 1 of the outer nested area, that is, Figure 6 The weighted average of the meridional winds of grid 1, grid 2, and grid 3 in the large area is obtained.

[0089] 2) The latitudinal wind in the meteorological field of each grid in the other layer nested area is obtained by taking the weighted average of the latitudinal winds in the meteorological field of the northernmost column of grids in the large area contained in the grids of the gas layer nested area.

[0090] For example, Figure 7 As shown, the zonal wind v in the grid 1 of the outer nested area is the westernmost grid of the 9 grids in the large area contained by the grid 1 of the outer nested area, that is, Figure 7 The weighted average of the zonal winds of grid 1, grid 2, and grid 3 in the large area is obtained.

[0091] 3) The prediction results of the remaining meteorological elements in the meteorological fields of the grids in the other nested regions are obtained by taking the weighted average of the prediction results of the corresponding meteorological elements in the meteorological fields of all grids in the large region contained in the grids in the outer nested region.

[0092] For example, Figure 8 As shown, the forecast results of the remaining meteorological elements except the meridional wind u and the zonal wind v in the grid 1 of the outer nested area are obtained by the 9 grids in the large area contained in the grid 1 of the outer nested area, that is, Figure 8The prediction results of the corresponding meteorological elements of grid 1 in the large area, grid 2 in the large area, grid 3 in the large area, ..., grid 9 in the large area are obtained by taking the weighted average.

[0093] In general meteorological model simulations and forecasts, the outer nested region (using a 9km resolution grid) is larger than the inner nested region (using a 3km resolution grid). Initial and boundary values ​​are provided to the inner nested region, and the inner nested region provides feedback to the outer nested region to ensure that the meteorological elements of the grids in the overlapping portion of the two nested regions match. However, when obtaining the forecast results for the meteorological elements of the grids in the inner nested region, the present invention reduces the impact of the discontinuity of the meteorological field between adjacent sub-regions by setting a buffer and performing quadratic fitting. The resulting forecast results for the meteorological elements of the grids in the overlapping region deviate from the original forecast results. If the outer nested region is processed and spliced ​​in a similar manner to the inner nested region to obtain the forecast results for the meteorological elements of the outer nested region, the forecast results for the meteorological elements of the grids in the outer nested region where the two sub-regions overlap or are adjacent will also deviate from the original forecast results. This will cause the forecast results for the meteorological elements of the 9km resolution grid and the 3km resolution grid in the overlapping portion of the inner and outer nested regions to not match. Therefore, the outer nested region cannot be directly processed. Therefore, the prediction results of meteorological elements in the inner and outer nested areas of the present invention are obtained based on the prediction results of meteorological elements of the 3km resolution grid (i.e., the 3km resolution grid within the large area of ​​the meteorological model obtained in step six), and can be arbitrarily intercepted based on the needs of air quality model simulation.

[0094] Figure 9 FIG. 1 shows a schematic diagram of the air quality model meteorological field generation system of the present invention. Figure 9 As shown, the air quality model meteorological field generation system of the present invention includes: 1. Module for determining the optimal parameterized scheme combination for sub-regions.

[0095] The sub-region optimal parameterization scheme combination determination module is used to divide a large area into multiple sub-regions, and determine the optimal parameterization scheme combination for each sub-region based on the simulation effect of meteorological field data over many years of history.

[0096] 2. Module for determining the optimal parameterized scheme combination for the buffer area.

[0097] The module for determining the optimal parameterized scheme combination for the buffer area is used to set a buffer area between two adjacent sub-areas and determine the optimal parameterized scheme combination for each buffer area based on the simulation effect of meteorological field data over many years of history.

[0098] 3. The first meteorological field determination module.

[0099] The first meteorological field determination module is used to predict future weather using the optimal parameterization scheme combination of each sub-region and each buffer area to obtain the first meteorological field of the inner nested area of ​​each sub-region and each buffer area.

[0100] 4. East-west fitting module.

[0101] The east-west fitting module is used to combine the first meteorological field of the inner nested area of ​​the buffer area between the sub-areas divided in the east-west direction and the first meteorological field of the inner nested area of ​​each sub-area to perform east-west fitting to obtain the second meteorological field of the inner nested area of ​​each sub-area.

[0102] 5. North-South fitting module.

[0103] The north-south fitting module is used to combine the first meteorological field of the inner nested area of ​​the buffer area between the sub-areas divided in the north-south direction and the second meteorological field of the inner nested area of ​​each sub-area to perform north-south fitting to obtain the final meteorological field of the inner nested area of ​​each sub-area.

[0104] 6. Meteorological field determination module for large areas.

[0105] The large-area meteorological field determination module is used to splice the obtained meteorological fields of the inner nested areas of each sub-area to obtain the meteorological field of the large area.

[0106] 7. Meteorological field determination module for each nested area.

[0107] The module for determining the meteorological fields of each nested area is used to intercept the large area according to the geographical scope of each nested area of ​​the target area under the air quality mode, and obtain the meteorological fields of each nested area of ​​the target area under the air quality mode based on the meteorological fields of the intercepted large area.

[0108] In addition, the present invention also provides an air quality pattern meteorological field generation device, which includes: one or more processors; a memory for storing one or more programs; when the one or more programs are executed by the one or more processors, the one or more processors implement the air quality pattern meteorological field generation method as described above. Finally, the present invention also provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the steps of the method for generating an air quality model meteorological field as described above.

[0109] Finally, it should be noted that the above embodiments are intended only to illustrate the technical solutions of the present invention and are not intended to limit the scope of protection of the present invention. Those skilled in the art may, based on the principles of the present invention, modify or replace the technical solutions of the present invention with equivalents without departing from the essence and scope of the technical solutions of the present invention.

Claims

1. A method for generating an air quality model meteorological field, characterized in that: The following steps are involved: 1) Divide a large area into multiple sub-areas and determine the optimal parameterization scheme combination for each sub-area based on the simulation results of historical meteorological field data over many years; 2) Set up buffer zones between two adjacent sub-regions and determine the optimal parameterization scheme combination for each buffer zone based on simulation results of historical meteorological data over many years; 3) Use the optimal parameterization scheme combination of each sub-region and each buffer area to predict future weather and obtain the first meteorological field of the inner nested area of ​​each sub-region and each buffer area; 4) Combine the first meteorological field of the inner nested area of ​​the buffer area between each sub-area divided in the east-west direction and the first meteorological field of the inner nested area of ​​each sub-area to perform east-west fitting to obtain the second meteorological field of the inner nested area of ​​each sub-area; 5) Combine the first meteorological field of the inner nested area of ​​the buffer area between each sub-area divided in the north-south direction and the second meteorological field of the inner nested area of ​​each sub-area to perform north-south fitting to obtain the final meteorological field of the inner nested area of ​​each sub-area; 6) splicing the final meteorological fields of the inner nested areas of each sub-area to obtain the meteorological field of the large area; 7) intercepting the large area according to the geographical range of each nested area of ​​the target area under the air quality model, and obtaining the meteorological field of each nested area of ​​the target area under the air quality model based on the meteorological field of the intercepted large area.

2. The method for generating an air quality model meteorological field according to claim 1, wherein: Determining the optimal parameterization scheme combination for each sub-region based on the simulation results of historical meteorological field data over many years in step 1) specifically includes: 11) Based on a sub-region, different parameterization scheme combinations are used to simulate the historical meteorological field of the sub-region over many years. The simulation results are tested using monitoring data to select the optimal parameterization scheme combination for the sub-region; 12) Using different parameterization scheme combinations to simulate the historical meteorological fields of the adjacent sub-regions of the sub-region over many years, using monitoring data to verify the simulation results, screening out the better parameterization scheme combinations of the adjacent sub-regions of the sub-region, and referring to the optimal parameterization scheme combination of the sub-region, determining the optimal parameterization scheme combination of the adjacent sub-regions of the sub-region based on the principle of minimizing the adjustment relative to the optimal parameterization scheme combination of the sub-region; 13) Use different parameterization scheme combinations to simulate the historical meteorological fields of other adjacent sub-regions of each adjacent sub-region of the sub-region, use monitoring data to test the simulation effect, screen out the better parameterization scheme combinations of other adjacent sub-regions of each adjacent sub-region of the sub-region, and refer to the optimal parameterization scheme combination of each adjacent sub-region of the sub-region, so as to minimize the adjustment relative to the optimal parameterization scheme combination of each adjacent sub-region of the sub-region, and determine the optimal parameterization scheme combination of other adjacent sub-regions of each adjacent sub-region of the sub-region.

3. The method for generating an air quality model meteorological field according to claim 1, wherein: Determining the optimal parameterization scheme combination for each buffer area based on the simulation effect of the historical meteorological field data over many years in the step 2) specifically includes: simulating the historical meteorological field of the buffer area over many years using different parameterization scheme combinations, verifying the simulation effect using monitoring data, screening out a better parameterization scheme combination for the buffer area, and referring to the optimal parameterization scheme combination of two adjacent sub-areas of the buffer area, so as to determine the optimal parameterization scheme combination for the buffer area based on the principle of minimizing the adjustment relative to the optimal parameterization scheme combination of the two adjacent sub-areas of the buffer area.

4. The method for generating an air quality model meteorological field according to claim 1, wherein: The step 4) specifically includes: 41) Divide the buffer area between each sub-area divided in the east-west direction into a central area and two edge areas, wherein the first meteorological field of the inner nested area of ​​the buffer area is used as the second meteorological field of the inner nested area of ​​the central area; 42) using the fitting results of the first meteorological field of the inner nested area of ​​the buffer area and the first meteorological field of the inner nested area of ​​the sub-area that overlaps with the edge area as the second meteorological field of the inner nested area of ​​the edge area; 43) using the first meteorological field of the inner nested area of ​​each sub-area as the second meteorological field of the inner nested area of ​​the portion of each sub-area that does not overlap with the buffer area; 44) Integrate the second meteorological field of the inner nested area of ​​the central area, the second meteorological field of the inner nested area of ​​the edge area, and the second meteorological field of the inner nested area of ​​each sub-area that does not overlap with the buffer area to obtain the second meteorological field of the inner nested area of ​​each sub-area.

5. The method for generating an air quality model meteorological field according to claim 4, wherein: In step 42), the fitting results of the first meteorological field of the inner nested area of ​​the buffer area and the first meteorological field of the inner nested area of ​​the sub-area overlapping with the edge area are specifically: Where, is the meteorological field of the i-th grid in a row of grids in the east-west direction of the inner nested area of ​​the edge area, starting from the central area, and n is the total number of grids in the row of grids in the east-west direction of the inner nested area of ​​the edge area. is the meteorological field of the corresponding grid of the inner nested area of ​​the buffer area, It is the meteorological field of the corresponding grid of the inner nested area of ​​the sub-area that overlaps with the edge area.

6. The method for generating an air quality model meteorological field according to claim 1, wherein: The step 7) specifically includes: 71) intercepting the large area according to the geographical scope of the inner nested area of ​​the target area under the air quality model, and using the meteorological field of the intercepted large area as the meteorological field of the inner nested area of ​​the target area under the air quality model; 72) The large area is intercepted based on the geographical scope of other nested areas of the target area under the air quality model, and the meteorological fields of other nested areas of the target area under the air quality model are processed based on the meteorological field of the intercepted large area.

7. The method for generating an air quality model meteorological field according to claim 6, characterized in that: The step 72) of processing the intercepted meteorological field of the large area to obtain meteorological fields of other nested areas of the target area under the air quality mode specifically includes: 721) The meridional wind in the meteorological field of each grid in the other nested regions is obtained by taking a weighted average of the meridional winds in the meteorological field of the westernmost column of grids in the large region contained in each grid in the other nested regions; 722) The latitudinal wind in the meteorological field of each grid in the other nested region is obtained by taking a weighted average of the latitudinal winds in the meteorological field of the northernmost column of grids in the large region contained in each grid in the other nested region; 723) The prediction results of the remaining meteorological elements in the meteorological fields of the grids in the other nested regions are obtained by taking the weighted average of the prediction results of the corresponding meteorological elements in the meteorological fields of all the grids in the large region contained in the grids in the other nested regions.

8. An air quality model meteorological field generation system, characterized in that: include: The module for determining the optimal parameterization scheme combination for a sub-region is used to divide a large region into multiple sub-regions and determine the optimal parameterization scheme combination for each sub-region based on the simulation results of meteorological field data over many years of history; A module for determining the optimal parameterization scheme combination for buffer areas, which is used to set up buffer areas between two adjacent sub-areas and determine the optimal parameterization scheme combination for each buffer area based on the simulation results of historical meteorological field data over many years; A first meteorological field determination module is used to predict future weather using an optimal parameterization scheme combination of each sub-region and each buffer region to obtain a first meteorological field for an inner nested region of each sub-region and each buffer region; An east-west fitting module is used to perform east-west fitting on the first meteorological field of the inner nested area of ​​the buffer area between each sub-area divided in the east-west direction and the first meteorological field of the inner nested area of ​​each sub-area to obtain the second meteorological field of the inner nested area of ​​each sub-area; A north-south fitting module is used to combine the first meteorological field of the inner nested area of ​​the buffer area between each sub-area divided in the north-south direction and the second meteorological field of the inner nested area of ​​each sub-area to perform north-south fitting to obtain the final meteorological field of the inner nested area of ​​each sub-area; A large-area meteorological field determination module, which is used to splice the meteorological fields of the inner nested areas of each sub-area to obtain the meteorological field of the large area; The meteorological field determination module for each nested area is used to intercept the large area according to the geographical scope of each nested area of ​​the target area under the air quality mode, and obtain the meteorological field of each nested area of ​​the target area under the air quality mode based on the meteorological field of the intercepted large area.

9. An air quality model meteorological field generating device, characterized in that: include: one or more processors; a memory for storing one or more programs; When the one or more programs are executed by the one or more processors, the one or more processors are enabled to implement the air quality model meteorological field generating method according to any one of claims 1 to 7.

10. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the program is executed by a processor, the steps of the air quality model meteorological field generation method according to any one of claims 1 to 7 are implemented.

Images