CN108205164A - A kind of atmospheric visibility parametrization Forecasting Methodology based on WRF-Chem - Google Patents
A kind of atmospheric visibility parametrization Forecasting Methodology based on WRF-Chem Download PDFInfo
- Publication number
- CN108205164A CN108205164A CN201711259479.0A CN201711259479A CN108205164A CN 108205164 A CN108205164 A CN 108205164A CN 201711259479 A CN201711259479 A CN 201711259479A CN 108205164 A CN108205164 A CN 108205164A
- Authority
- CN
- China
- Prior art keywords
- wrf
- atmospheric
- chem
- scheme
- visibility
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01W—METEOROLOGY
- G01W1/00—Meteorology
- G01W1/10—Devices for predicting weather conditions
Landscapes
- Environmental & Geological Engineering (AREA)
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Atmospheric Sciences (AREA)
- Biodiversity & Conservation Biology (AREA)
- Ecology (AREA)
- Environmental Sciences (AREA)
- Organic Low-Molecular-Weight Compounds And Preparation Thereof (AREA)
Abstract
The invention discloses a kind of atmospheric visibilities based on WRF Chem to parameterize Forecasting Methodology, obtains meteorological data and topographic(al) data in simulated time, sets assessment area and meteorological data and topographic(al) data are pre-processed；It selects to combine to simulating the different parameters scheme with larger impact, suitable Parameterization Scheme is selected according to simulation area；WRF Chem are run, using WRF Chem pattern simulation particle concentrations, according to Mie theoretical calculations particulate matter to the extinction coefficient of 550 nm, and consider the delustring of atmospheric molecule, atmospheric visibility is calculated with improved atmospheric visibility Parameterization Scheme.Due to combining at present research pair both at home and abroadKValue is corrected, and atmospheric extinction coefficient considers particulate matter extinction coefficient and gas molecular extinction coefficient simultaneously, and this method more can be derived that more accurate atmospheric visibility.
Description
Technical field
The invention belongs to atmospheric environment pattern technology fields, and in particular to a kind of atmospheric visibility ginseng based on WRF-Chem
Numberization Forecasting Methodology.
Background technology
Atmospheric visibility characterizes atmospheric transparency, can reflect the clean level of Air Close To The Earth Surface, with social life and
Traffic safety is closely related.Since existing atmosphere pollution observation time sequence is shorter, the observation and analysis of atmospheric visibility
It is still the main path of understanding atmospheric environment longer term climatic variation characteristic at present.The forecasting procedure of atmospheric visibility mainly has system
Meter forecast and numerical forecast.Statistical fluctuation changes generally by analysis mist, atmospheric visibility of the haze when weather phenomena occur
Rule establishes the statistical relationship of the meteorological elements such as temperature, humidity, wind speed, pressure and atmospheric visibility, establishes atmospheric visibility
Prognostic equation numerical forecast mainly utilizes the elements such as pollutant, humidity, the Liquid water content in numerical model simulated atmosphere, according to
According to atmospheric optics theory, its contribution to atmospheric extinction, diagnosis forecast atmospheric visibility are calculated.
At present, domestic and international main atmospheric visibility computational methods are the relation formulas of atmospheric visibility and extinction coefficient,
I.e. object visual range theory atmospheric extinction coefficient is mainly proposed by IMPROVE projects (U.S.'s large size visibility surveillance program)
IMPROVE empirical equations calculate.CHEN etc. by Tianjin Wuqing area atmospheric visibility to aerosol fraction and suction
A kind of sensitivity analysis of wet growth factor, it is proposed that the Parameterization Scheme that low visibility delustring calculates under haze weather.Kunkel
Deng passing through mist observation experiment in 1984, it is proposed that the relational expression of atmospheric visibility and Liquid water content.Gultepe etc. gives
Atmospheric visibility formula based on fog content and droplet concentration.ZHOU etc. gives liquid water on the basis of radiation fog is studied
The diagnostic method of content (LWC), and then propose the atmospheric visibility DATA PROCESSING IN ENSEMBLE PREDICTION SYSTEM method of mist.
For being now widely used in the WRF-Chem isotypes of air quality model forecast, existing calculating atmospheric visibility
Parameterization Scheme availability it is not high, cause pattern poor to the computational accuracy of atmospheric visibility.
Invention content
Purpose：In order to overcome the deficiencies in the prior art, the present invention provides a kind of air energy based on WRF-Chem
Degree of opinion parameterizes Forecasting Methodology, realizes the coupling of atmospheric visibility Parameterization Scheme and WRF-Chem patterns, this method can be more
Add the prediction for accurately carrying out atmospheric visibility.
Technical solution：In order to solve the above technical problems, the technical solution adopted by the present invention is：
A kind of atmospheric visibility parametrization Forecasting Methodology based on WRF-Chem, which is characterized in that include the following steps：
(1) meteorological data and topographic(al) data in simulated time are obtained, assessment area is set and to meteorological data and landform
Data is pre-processed；
(2) selection is combined to simulating the different parameters scheme with larger impact, suitable according to simulation area selection
Parameterization Scheme；
Specifically, the Parameterization Scheme includes Microphysical scheme, long-wave radiation Parameterization Scheme, Shortwave radiative parameter
Scheme, land surface scheme, Different Boundary Layer Parameterization Schemes, Convective Parameterization Schemes, aerosol scheme gas-phase chemical reaction
Mechanism.
(3) using WRF-Chem pattern simulation particle concentrations, according to delustring of the Mie theoretical calculations particulate matter to 550nm
Coefficient, and consider the delustring of atmospheric molecule, calculate atmospheric visibility with improved atmospheric visibility Parameterization Scheme.
Atmospheric visibility VRCalculation formula be：
In formula：bpFor particulate matter extinction coefficient, unit Mm-1；bsgFor gas molecule scattering coefficient, unit Mm-1；bagFor gas
Body molecular absorption coefficient, unit Mm-1。
Particulate matter extinction coefficient is extinction coefficient of the particulate matter to 550nm wavelength, is obtained by Mie theoretical calculations, formula
It is as follows：
In formula：I is grain size section, value 1,2,3,4, corresponding grain size is respectively 0.039~0.156,0.156~
0.625th, 0.625~2.5,2.5~10 μm；QextFor extinction efficiency factor；DpFor the average diameter of particulate matter, m；N is particulate matter
Particle density, a/m3；voldryFor the volume fraction of dry particl object, %；volwaterFor liquid water volume fraction, %.
bsgFor the Rayleigh scattering of gas molecule scattering coefficient, mainly atmospheric molecule, constant is usually regarded as, it is preferred that bsg
Value is 13Mm-1。
bagFor gas molecules sorb coefficient, mainly by NO2Pollution contribution, absorption coefficient are about NO2Mass concentration
0.33 times.
Advantageous effect：Atmospheric visibility parametrization Forecasting Methodology provided by the invention based on WRF-Chem, based on existing
The total extinction coefficient of air and atmospheric visibility relationship, with reference to east China air quality and atmospheric visibility observational study,
According to Mie theories and the Extinction Characteristic of atmospheric molecule, particulate matter and NO are calculated2Extinction coefficient, propose a kind of improved air
Visibility Parameterization Scheme, and atmospheric extinction coefficient considers particulate matter extinction coefficient and gas molecular extinction coefficient simultaneously,
Gas molecule extinction coefficient is the sum of gas molecule scattering coefficient and gas molecular absorption coefficient.It is said from overall numerical procedure, it should
Pattern is to calculate the preferable selection of atmospheric visibility.
With Nanjing of China 25 days-December 10 November in 2013 and two haze contamination accidents on December 20-29th, 2013
Example shows the Average normalized deviation and average deviation of two hazes pollution example visibility of modified scheme simulation
In respectively 17.19%, 3.18% and 517m, 173m, and and observation visibility related coefficient be respectively increased to 0.76,
0.87, the accuracy of visibility simulation better than other two kinds of existing Parameterization Schemes, and in different relative humidity (RH) and
The standardization mean error that modified scheme is simulated in visibility range is respectively less than other two kinds of Parameterization Schemes, wherein
In RH ＜ 80% and visibility >=1km range internal standardization mean errors less than 50%.The improved air of the present invention can be shown in
Degree Parameterization Scheme can effectively improve the forecast accuracy of the haze pollution ambient air visibility of air quality model.
Description of the drawings
Fig. 1 is the method for the present invention flow chart；
Fig. 2 is the comparison of three kinds of parametrization program simulation visibility and observation in the different haze pollution periods；
Fig. 3 is three kinds of parametrization program simulation results in different visibility scales and the average deviation of humidity range；
Fig. 4 is that three kinds of parametrization program simulation results are averagely missed in the standardization of different visibility scales and humidity range
Difference.
Specific embodiment
With reference to embodiment, the invention will be further described.Following embodiment is only used for clearly illustrating this hair
Bright technical solution, and be not intended to limit the protection scope of the present invention and limit the scope of the invention.
As shown in Figure 1, a kind of atmospheric visibility parametrization Forecasting Methodology based on WRF-Chem, includes the following steps：It obtains
Meteorological data and topographic(al) data in the modulus pseudotime set assessment area and meteorological data and topographic(al) data are located in advance
Reason.
It selects to combine to simulating the different parameters scheme with larger impact, suitable parameter is selected according to simulation area
Change scheme.
It is as follows using the grid setting of WRF-Chem3.7 version analog studies：Vertical direction is divided into 28 layers, and top of model is
100hPa；Horizontal direction sets three layers of nested grid, resolution ratio ecto-entad is respectively 45,15,5km, covering China is absolutely successively
Most area, East China and Central China some areas, Jiangsu Province and Anhui Zhejiang some areas；Regional center be located at Nanjing (32 ° of N,
118°E).Table 1 gives the physical and chemical process Parameterization Scheme in pattern simulation setting.Emission inventory is used by population
INTEX-B data in 2006 after distribution, chemical initial fields use MOZART, simulate 25 days-December of November in 2013 10 respectively
Day and two haze contamination accidents on December 20-29th, 2013.
WRF-Chem mode parameters scheme setting such as following table：
| Scheme type | Pattern options |
| Microphysical scheme | Morrison schemes |
| Long-wave radiation Parameterization Scheme | RRTM schemes |
| Shortwave radiative parameter scheme | RRTM schemes |
| Land surface scheme | Noah schemes |
| Different Boundary Layer Parameterization Schemes | YSU schemes |
| Convective Parameterization Schemes | Grell 3D schemes |
| Aerosol scheme | MOSAIC schemes |
| Gas-phase chemical reaction mechanism | CBMZ mechanism |
WRF-Chem is run, using WRF-Chem pattern simulation particle concentrations, according to Mie theoretical calculation particulate matters pair
The extinction coefficient of 550nm, and consider the delustring of atmospheric molecule, it is counted with improved atmospheric visibility Parameterization Scheme (scheme C)
Calculate atmospheric visibility.And compared with the visibility Parameterization Scheme (option b) of existing IMPROVE schemes (option A) and CHEN,
As a result such as following table：
The atmospheric visibility that scheme C (improved atmospheric visibility Parameterization Scheme) is simulated it can be seen from analog result
With observe it is closest, average deviation and standardization average deviation be minimum, respectively 517m, 17.19%, for it is entire when
Between section simulation effect it is best.Atmospheric visibility moulds of the scheme C (improved atmospheric visibility Parameterization Scheme) within the entire period
Intend deviation to be substantially reduced, simulation is the most accurate, and related coefficient is up to 0.87, and average deviation is only 173m, standardizes average deviation
It is respectively 3.18%, 27.16% with standardization mean error.
Fig. 2 is the comparison of three kinds of parametrization program simulation visibility and observation in the different haze pollution periods；As seen from the figure,
Scheme C (improved atmospheric visibility Parameterization Scheme) and measured value are integrally closest, are better than and option A (IMPROVE side
Case) and option b (the visibility Parameterization Scheme of CHEN) result of calculation, scheme C (improved atmospheric visibility parametrization sides
Case) it can effectively improve the accuracy in computation of atmospheric visibility.
Fig. 3 is three kinds and parameterizes program simulation results in different visibility scales and the average deviation of humidity range, and Fig. 4 is
Three kinds of parametrization program simulation results are in different visibility scales and the standardization mean error of humidity range.
Scheme C is under the conditions of each visibility range and different RH to the possesses good fitting of atmospheric visibility, average deviation and mark
Standardization mean error is minimum, and optical principle Mie theoretical calculation particulate matter extinction coefficients are based on better than other two schemes schemes C,
Consider the Extinction Characteristic of atmospheric molecule, diagnosis obtains visibility, and result is more accurate, reasonable on the whole.
The above is only the preferred embodiment of the present invention, it should be pointed out that：For the ordinary skill people of the art
For member, various improvements and modifications may be made without departing from the principle of the present invention, these improvements and modifications also should
It is considered as protection scope of the present invention.
Claims (6)
1. a kind of atmospheric visibility parametrization Forecasting Methodology based on WRF-Chem, which is characterized in that include the following steps：
(1) meteorological data and topographic(al) data in simulated time are obtained, assessment area is set and to meteorological data and topographic(al) data
It is pre-processed；
(2) selection is combined to simulating the different parameters scheme with larger impact, and suitable parameter is selected according to simulation area
Change scheme；
(3) using WRF-Chem pattern simulation particle concentrations, according to Mie theoretical calculations particulate matter to the extinction coefficient of 550nm,
And consider the delustring of atmospheric molecule, calculate atmospheric visibility with improved atmospheric visibility Parameterization Scheme.
2. the atmospheric visibility parametrization Forecasting Methodology according to claim 1 based on WRF-Chem, it is characterised in that：Institute
It states Parameterization Scheme and includes Microphysical scheme, long-wave radiation Parameterization Scheme, Shortwave radiative parameter scheme, land surface emissivity side
Case, Different Boundary Layer Parameterization Schemes, Convective Parameterization Schemes, aerosol scheme gas-phase chemical reaction mechanism.
3. the atmospheric visibility parametrization Forecasting Methodology according to claim 1 based on WRF-Chem, it is characterised in that：Greatly
Gas visibility VRCalculation formula be：
In formula：bpFor particulate matter extinction coefficient, unit Mm-1；bsgFor gas molecule scattering coefficient, unit Mm-1；bagFor gas point
Sub- absorption coefficient, unit Mm-1。
4. the atmospheric visibility parametrization Forecasting Methodology according to claim 3 based on WRF-Chem, it is characterised in that：
Grain object extinction coefficient is extinction coefficient of the particulate matter to 550nm wavelength, is obtained by Mie theoretical calculations, formula is as follows：
In formula：I is grain size section, value 1,2,3,4, corresponding grain size is respectively 0.039~0.156,0.156~0.625,
0.625~2.5,2.5~10 μm；QextFor extinction efficiency factor；DpFor the average diameter of particulate matter, m；N is that the number of particulate matter is dense
Degree, a/m3；voldryFor the volume fraction of dry particl object, %；volwaterFor liquid water volume fraction, %.
5. the atmospheric visibility parametrization Forecasting Methodology according to claim 3 based on WRF-Chem, it is characterised in that：bsg
Value is 13Mm-1。
6. the atmospheric visibility parametrization Forecasting Methodology according to claim 3 based on WRF-Chem, it is characterised in that：bag
For gas molecules sorb coefficient, mainly by NO2Pollution contribution, absorption coefficient NO20.33 times of mass concentration.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711259479.0A CN108205164A (en) | 2017-12-04 | 2017-12-04 | A kind of atmospheric visibility parametrization Forecasting Methodology based on WRF-Chem |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711259479.0A CN108205164A (en) | 2017-12-04 | 2017-12-04 | A kind of atmospheric visibility parametrization Forecasting Methodology based on WRF-Chem |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN108205164A true CN108205164A (en) | 2018-06-26 |
Family
ID=62604664
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711259479.0A Pending CN108205164A (en) | 2017-12-04 | 2017-12-04 | A kind of atmospheric visibility parametrization Forecasting Methodology based on WRF-Chem |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108205164A (en) |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109444986A (en) * | 2018-09-30 | 2019-03-08 | 国网天津市电力公司电力科学研究院 | A kind of photovoltaic power station foggy weather visibility prediction technique |
| CN109543906A (en) * | 2018-11-23 | 2019-03-29 | 长三角环境气象预报预警中心（上海市环境气象中心） | A kind of method and apparatus of atmospheric visibility prediction |
| CN109597969A (en) * | 2019-01-25 | 2019-04-09 | 南京大学 | A kind of surface ozone Concentration Estimation Method |
| CN109932988A (en) * | 2019-03-27 | 2019-06-25 | 四川瞭望工业自动化控制技术有限公司 | A kind of city raised dust contamination forecasting system and method |
| CN111382506A (en) * | 2020-03-02 | 2020-07-07 | 苏州工业园区洛加大先进技术研究院 | Method for evaluating influence of aerosol and radiation interaction on atomization effect |
| CN111398109A (en) * | 2020-03-10 | 2020-07-10 | 上海眼控科技股份有限公司 | Atmospheric visibility measuring method, sensor module, system and storage medium |
| CN115204507A (en) * | 2022-07-26 | 2022-10-18 | 北京中科三清环境技术有限公司 | Atmospheric visibility prediction method, device, equipment and storage medium |
| CN117554992A (en) * | 2024-01-10 | 2024-02-13 | 青岛镭测创芯科技有限公司 | Extinction coefficient acquisition method and system based on laser radar |
-
2017
- 2017-12-04 CN CN201711259479.0A patent/CN108205164A/en active Pending
Non-Patent Citations (2)
| Title |
|---|
| 胡俊，赵天良，张泽锋，邱玉珺，谭成好，吴峡生等: "霾污染环境大气能见度参数化方案的改进", 《环境科学研究》 * |
| 荣昕: "基于WRF-CHEM模式的连续雾霾过程数值模拟及其能见度参数化", 《基础科学辑》 * |
Cited By (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109444986A (en) * | 2018-09-30 | 2019-03-08 | 国网天津市电力公司电力科学研究院 | A kind of photovoltaic power station foggy weather visibility prediction technique |
| CN109444986B (en) * | 2018-09-30 | 2021-01-19 | 国网天津市电力公司电力科学研究院 | Visibility prediction method for foggy weather of distributed photovoltaic power station |
| CN109543906A (en) * | 2018-11-23 | 2019-03-29 | 长三角环境气象预报预警中心（上海市环境气象中心） | A kind of method and apparatus of atmospheric visibility prediction |
| CN109543906B (en) * | 2018-11-23 | 2024-04-16 | 长三角环境气象预报预警中心（上海市环境气象中心） | Atmospheric visibility prediction method and equipment |
| CN109597969A (en) * | 2019-01-25 | 2019-04-09 | 南京大学 | A kind of surface ozone Concentration Estimation Method |
| CN109932988A (en) * | 2019-03-27 | 2019-06-25 | 四川瞭望工业自动化控制技术有限公司 | A kind of city raised dust contamination forecasting system and method |
| CN109932988B (en) * | 2019-03-27 | 2020-11-24 | 四川瞭望工业自动化控制技术有限公司 | Urban raise dust pollution diffusion prediction system and method |
| CN111382506A (en) * | 2020-03-02 | 2020-07-07 | 苏州工业园区洛加大先进技术研究院 | Method for evaluating influence of aerosol and radiation interaction on atomization effect |
| CN111398109A (en) * | 2020-03-10 | 2020-07-10 | 上海眼控科技股份有限公司 | Atmospheric visibility measuring method, sensor module, system and storage medium |
| CN115204507A (en) * | 2022-07-26 | 2022-10-18 | 北京中科三清环境技术有限公司 | Atmospheric visibility prediction method, device, equipment and storage medium |
| CN117554992A (en) * | 2024-01-10 | 2024-02-13 | 青岛镭测创芯科技有限公司 | Extinction coefficient acquisition method and system based on laser radar |
| CN117554992B (en) * | 2024-01-10 | 2024-04-09 | 青岛镭测创芯科技有限公司 | Extinction coefficient acquisition method and system based on laser radar |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108205164A (en) | A kind of atmospheric visibility parametrization Forecasting Methodology based on WRF-Chem | |
| Tang et al. | Impact of emission controls on air quality in Beijing during APEC 2014: lidar ceilometer observations | |
| Soegaard et al. | Towards a spatial CO2 budget of a metropolitan region based on textural image classification and flux measurements | |
| Yan et al. | The measurement of aerosol optical properties at a rural site in Northern China | |
| Xiang-Ao et al. | Aerosol properties and their spatial and temporal variations over North China in spring 2001 | |
| Hoffmann et al. | Effects of grazing and topography on dust flux and deposition in the Xilingele grassland, Inner Mongolia | |
| Quinn et al. | Comparison of measured and calculated aerosol properties relevant to the direct radiative forcing of tropospheric sulfate aerosol on climate | |
| CN104897853A (en) | Thermal power plant pollutant discharging monitoring display method based on tower type diffusion model | |
| CN110929228B (en) | Inversion algorithm for moisture absorption growth factor of uniformly mixed aerosol | |
| Goossens et al. | Techniques to measure the dry aeolian deposition of dust in arid and semi‐arid landscapes: a comparative study in West Niger | |
| CN107561554A (en) | Inversion method with multi-wavelength laser radar data is counted based on solar luminosity | |
| CN112180472A (en) | Atmospheric visibility integrated forecasting method based on deep learning | |
| CN103234882B (en) | A kind of Atmospheric particulates mass concentration inversion method based on the particulate matter flight time | |
| CN103336995A (en) | Method for constructing real-time light metering network of million kilowatt level photovoltaic power generation base | |
| Zhou et al. | Improving aerosol interaction with clouds and precipitation in a regional chemical weather modeling system | |
| Patterson et al. | Monte Carlo simulation of daily regional sulfur distribution: comparison with SURE sulfate data and visual range observations during August 1977 | |
| CN108663727A (en) | The method for estimating height of evaporation duct within the scope of world marine site using evaporation rate | |
| Zou et al. | Impact of eddy characteristics on turbulent heat and momentum fluxes in the urban roughness sublayer | |
| Cleugh | Field measurements of windbreak effects on airflow, turbulent exchanges and microclimates | |
| CN110907319B (en) | Attribution analysis method for near-surface fine particulate matters | |
| Virkkula et al. | The influence of Kola Peninsula, continental European and marine sources on the number concentrations and scattering coefficients of the atmospheric aerosol in Finnish Lapland | |
| CN110907318B (en) | Near-ground atmospheric total suspended particulate matter mass concentration remote sensing physical estimation method | |
| CN109948175B (en) | Satellite remote sensing albedo missing value inversion method based on meteorological data | |
| CN109918770B (en) | Artificial precipitation removal fine particle prediction model | |
| KR101194677B1 (en) | Estimation method of source of high concentration fine dust in urban area |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180626 |