CN109858686A - A kind of ground emission inventories inverting optimization method based on EnKF - Google Patents

A kind of ground emission inventories inverting optimization method based on EnKF Download PDF

Info

Publication number
CN109858686A
CN109858686A CN201910030086.5A CN201910030086A CN109858686A CN 109858686 A CN109858686 A CN 109858686A CN 201910030086 A CN201910030086 A CN 201910030086A CN 109858686 A CN109858686 A CN 109858686A
Authority
CN
China
Prior art keywords
inverting
inventory
observation
enkf
error
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
Application number
CN201910030086.5A
Other languages
Chinese (zh)
Inventor
冯树状
江飞
王恒茂
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing University
Original Assignee
Nanjing University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Nanjing University filed Critical Nanjing University
Priority to CN201910030086.5A priority Critical patent/CN109858686A/en
Publication of CN109858686A publication Critical patent/CN109858686A/en
Pending legal-status Critical Current

Links

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention discloses a kind of ground emission inventories inverting optimization method based on EnKF, comprising: building EnKF inventory Inversion System is simultaneously coupled to WRF-CMAQ modular system;Inventory after daily inverting is re-applied to same day model predictions and provides initial fields, the average inventory as final optimization pass of more days inverting inventories for the DATA PROCESSING IN ENSEMBLE PREDICTION SYSTEM in inverting in next day;Assess inventory regional change and uncertainty;Forecast when inventory after optimization is carried out long within the assimilation period is forecast to compare with source list, improvement effect of the inventory after assessment optimization to forecast.The present invention constructs set inverting assimilation system by EnKF methodology to optimize SO simultaneously2、NOX、PM2.5Emission inventories effectively correct for discharge error and improve model predictions.

Description

A kind of ground emission inventories inverting optimization method based on EnKF
Technical field
The invention belongs to field of environment protection, especially a kind of emission inventories inverting optimization method.
Background technique
Conventional exhaust inventory usually passes through on-site inspection method and obtains, since the distribution and variation of pollution sources are sufficiently complex, Establishment inventory usually requires longer time and a large amount of manpower and is estimated to check each pollutant and establish model, Emission inventories data are fixed and invariable in a long time, so can be influenced by industrial expansion or policy, actual discharge Fluctuation is big, causes the biggish error of mode.
Summary of the invention
Goal of the invention: a kind of ground emission inventories inverting optimization method based on EnKF is provided, is deposited with solving the prior art The above problem.
A kind of technical solution: ground emission inventories inverting optimization method based on EnKF, comprising:
Step 1, building EnKF inventory Inversion System are simultaneously coupled to WRF-CMAQ modular system;
Step 2, by the inventory after daily inverting be re-applied to the same day model predictions be inverting in next day in set it is pre- Report provides initial fields, the average inventory as final optimization pass of more days inverting inventories;
Step 3, assessment inventory regional change and uncertainty;It is pre- when inventory after optimization is carried out long within the assimilation period Report forecasts to compare with source list, improvement effect of the inventory after assessment optimization to forecast.
According to an aspect of the present invention, the step of building EnKF inventory Inversion System includes:
Step 11 generates source emission set sample
For source list, δ x is random perturbation field, and i is set sample number;
Step 12 considers PM2.5The influence of precursor, while to SO2、NOxDischarge carry out inverting correct, calculate source row The background error covariance matrix put For source emission ensemble average before inverting, N is natural number;
Step 13 is updated based on aggregate error statistics and measurement vector, the inventory after inverting by following equation:
Wherein, H Observation Operators, R are observation error covariance matrix, and K is gain matrix, determine ambient field and observation Weight,Ensemble average for the analysis field of ith member, i.e. emission inventories after optimization, all members is estimated to be optimal Meter, provides precision higher emission inventories for inverting next time.
According to an aspect of the present invention, further include step 14, carry out inverting using average daily concentration, assimilation window is set as 1 It, the influence to avoid model predictions diurnal variation error to inverting;Super observation is carried out to observation based on optimal estimation theory Processing,
M observation y of same mesh point will be located atiPredicted value corresponding with n memberConstitute a super observation Value and corresponding predicted value, i=1,2 ..., m;K=1,2 ..., n, it is assumed that the observation error of different website different moments is mutual It is independent, yiCorresponding observation standard deviation is ri, then new to constitute observation ynew, observe standard deviation rnewAnd corresponding forecast Meet following relationship:
According to an aspect of the present invention, further include step 15, single website successively assimilated using sequence assimilation mode, i.e., Assimilate the analysis field updated after an observation to continue to assimilate as the ambient field assimilated next time.
According to an aspect of the present invention, observation error covariance R includes measurement error and representive error,
Representive error
ε0For measurement error, ε0=ermax+ermin* Π0, wherein Π0For observation, ermax is elementary error;
γ is tuning scale factor, and Δ x is sizing grid, and L indicates the radius of influence of observation website, is equipped with observation bit It closes, total observation error is defined as
According to an aspect of the present invention, CMAQ model system includes:
Chemical Transport module, transmission process, chemical process and infall process for simulating pollution object;
Primary condition module and boundary condition module, for providing pollutant initial fields and boundary field for CCTM;
Photochemical breakdown rate module, for calculating photochemical breakdown rate;
Meteorological-chemical interface module, is meteorologic model and the interface of CCTM, can for meteorological data to be converted into CCTM The data format of reading.
The utility model has the advantages that constructing set inverting assimilation system by EnKF methodology to optimize SO simultaneously2、NOX、PM2.5Discharge is clear It is single, it effectively corrects for discharge error and improves model predictions.
Detailed description of the invention
Fig. 1 is the principle of the present invention schematic diagram.
Fig. 2 is 10-m wind speed (WS10), the 2-m positioned at the Nanjing (a, c, e) and the observation of (b, d, f) Shanghai website with simulation Temperature (T2) and 2-m relative humidity (RH2) time series comparison diagram.
Fig. 3 be China's Mainland, the region NCP, YRD and SCB be averaged SO2, NOX, a PM2.5 inverting discharge ratio Factor time variation diagram.
Fig. 4 is inventory and source list (a, c, e) discharge difference (106g/d) and (b, d, f) discharge ratio after inverting Factor space distribution map.
Fig. 5 is SO of 25 to 30 December in 2015 using inventory after inverting and source list simulation2Concentration and RMSE comparison diagram.
Fig. 6 is NO of 25 to 30 December in 2015 using inventory after inverting and source list simulation2Concentration and RMSE comparison diagram.
Fig. 7 is PM of 25 to 30 December in 2015 using inventory after inverting and source list simulation2.5Concentration and RMSE comparison diagram.
Fig. 8 is inventory simulation (a, b) SO after source list and inverting2, (c, d) NO2, (e, f) PM2.5Average deviation pair Than figure.
Fig. 9 is assimilation and does not assimilate SO2、NO2Observe the PM of inverting2.5(a) inventory difference (106G/d it) and (b) discharges Scale factor spatial distribution map.
Figure 10 is SO2、NO2And PM2.5Inverting after localization scale becomes 180,108 and 180km from 108,72,108km PM2.5(a) inventory difference (106G/d scale factor spatial distribution map) and (b) is discharged.
Specific embodiment
Details and principle of the invention is described in detail below.
As shown in Figure 1, constructing EnKF inventory Inversion System first and being coupled to WRF-CMAQ mode.Due to WRF-CMAQ The offline coupling of model, ignoring the meteorological mutual feedback with chemistry influences, for the simulation error for avoiding Extreme Weather Events, this reality It applies example and chooses the research period as 25 to 30 December in 2015.Secondly to a PM in emission inventories2.5And precursor SO2、 NOXIt carries out random perturbation 30 collection and is merged into row set forecast, be set as one to avoid the influence of discharge diurnal variation from assimilating window It.Inventory after daily inverting is re-applied to same day model predictions and provides initial fields for the DATA PROCESSING IN ENSEMBLE PREDICTION SYSTEM in inverting in next day, The average inventory as final optimization pass of more days inverting inventories.Assess inventory regional change and uncertainty;It will be clear after optimization Forecast when singly carrying out long within the assimilation period is forecast to compare with source list, raising effect of the inventory after assessment optimization to forecast Fruit.Aerosol concentration and component are forecast using WRF-CMAQ offline area air quality model.WRF modular system is state, the U.S. Family's Center for Atmospheric Research (NCAR), U.S.National Oceanic and atmosphere office (NOAA) and NCEP combine multiple scientific research institution's exploitations Mesoscale high-resolution weather research of new generation and Forecast Mode, can accurately reproduce the various synoptic process in real atmosphere, extensively It is general to be applied to atmospheric science research and operational forecast demand.Research of the WRF model integrated so far in mesoscale at Fruit, including two dynamic cores: the system knot of a data assimilation system and support a parallel computation and the system expandability Structure.WRF horizontal direction uses Arakawa C grid, and vertical direction then uses terrain following mass coordinate by spatial spreading Change, is complete compressible non-hydrostatic atmospheric model, can preferably improve the simulation and forecast to Study of Meso Scale Weather.WRF model Inner parameter scheme is more richer than other mesoscale models, and the physical process considered is also in further detail.It is introduced in mode Physical parameter scheme mainly includes microphysical processes scheme, cumulus convection scheme, PBL scheme, long wave shortwave radiation mould Quasi- scheme and land face scheme etc..Simulation and real-time prediction experiments have shown that, WRF modular system is in terms of forecasting various weather conditions Do well.
The third generation prediction of air quality model that CMAQ is developed by US Gov Env Protection Agency uses " atmosphere " Design concept, physical and chemical process complicated between atmosphere difference pollutant is combined together consideration, meteorologic model is not only simulated Driving under the conveying of pollutant, diffusion, conversion and transition process influence, combined dirty between region and City-scale Contaminate the various chemical processes of object in an atmosphere, including liquid-phase chemistry, heterogeneous phase chemistry process, aerosol processes and dry and wet Influence of the deposition process to concentration distribution, to substantially increase the precision of CMAQ simulation.CMAQ can be in different spaces scale model In enclosing while the various problem of environmental pollutions such as the various air pollutants such as PM, ozone and acid deposition, visibility are simulated, may be used also To be used to study pollutant sources, mechanism of production, diffusion, the area transmissions and assessment pollution reduction effect of pollutant are disclosed The influence of fruit, Environment Controlling Strategy to air quality, to work out best decision.CMAQ model is designed using high modularization, It is mainly made of 5 parts, core is Chemical Transport module CCTM (CMAQ Chemical-Transport Model It Processor), can be with the transmission process, chemical process and infall process etc. of simulating pollution object;Primary condition module I CON (Initial Conditions Processor) and boundary condition module BCON (Boundary Conditions Processor pollutant initial fields and boundary field) are provided for CCTM;Photochemical breakdown rate module J PROC (Photolysis Rate Processor) calculate photochemical breakdown rate;Meteorology-chemistry interface module MCIP (Meteorology-Chemistrry Interface Processor) it is meteorologic model and the interface of CCTM, meteorological data is converted into the data that CCTM can be read Format.For CMAQ relative to online coupling model, computational efficiency is higher, has been widely used in the research and business of air pollution Forecast work.
WRF pattern simulation region, East and West direction include 169 mesh points, and north-south includes 129 mesh points, which covers East Asia most area is covered, resolution ratio is 36km × 36km;Vertical direction is divided into 51 layers, a height of 50hPa of top of model.CMAQ Each few three grids of the every side ratio WRF model of the simulated domain of model, it is vertical on be compressed to 15 layers from 51 WRF layers, wherein Seven layers are located at boundary layer, and remaining in free convection layer.Temporal resolution is 6 hours, 1 ° × 1 ° of spatial resolution.Gas is molten Glue initial fields are from the forecast for recycling 6 hours before.It is further noted that artificial source emission was worked out from Tsinghua University The artificial source emission inventory (MIX) in Asia in 2010 (Li et al., 2017), provide five industries (electric power, industry, it is civilian, hand over Logical, agricultural) 0.25 ° of spatial resolution monthly gridding emissions data.Biological source emission passes through MEGAN model (Guenther Et al., 2006) off-line calculation.AERO4 aerosol module is used in the present embodiment, includes PM2.5Seven aerosol species, Including organic carbon, black wood charcoal, sulfate, nitrate, ammonium salt, sea salt and also unformed aerosol species.
The random sample that the error for the ambient field that EnKF is inputted according to mode generates one group of finite aggregate is not true to represent it It is qualitative, B matrix is calculated by model predictions result, not only solves what Kalman filtering faced in high dimensional pattern system B matrix stores and calculates the excessive problem of cost, so that B matrix is developed with the update of dynamic mode, both can be applied to height Nonlinear modular system is spent, adjoint mode need not be also write, be easier implementation and application.A committed step of EnKF It is to generate discharge set sample to provide Ensemble simulation, different from Ensemble Numerical Weather Prediction generating mode, emission source does not need full The foot certain condition of continuity and equilibrium condition, therefore the present embodiment generates set disturbance using Monte Carlo method.Source emission collection Sample is closed to obtain by disturbance source list:
Wherein,For source list, δ X is random perturbation field, and in the present embodiment, perturbation amplitude is set as source list 15%, value is set as 0, i for set sample number, for computational efficiency and assimilation effect, this reality is better balanced lower than 0 after disturbance It tests set number and is set as 30.In addition, the present embodiment considers PM2.5The influence of precursor, while to SO2、 NOXDischarge carry out it is anti- It drills and corrects, since the correlation for directly acquiring each state variable is more difficult, so state variable is mutual between assuming It is independent.The background error covariance matrix of source emission is obtained according to following equation:
Wherein,For source emission ensemble average before inverting.Inventory based on aggregate error statistics and measurement vector, after inverting It is updated by following equation:
Wherein, H Observation Operators, R are observation error covariance matrix, and K is gain matrix, determine ambient field and observation Weight,Ensemble average for the analysis field of ith member, i.e. emission inventories after optimization, all members is estimated to be optimal Meter, provides precision higher emission inventories, y is observation matrix for inverting next time.
To avoid influence of the model predictions diurnal variation error to inverting, inverting is carried out using average daily concentration, assimilation window is set It is 1 day.To reduce the influence for calculating cost and representive error, observation is carried out based on optimal estimation theory " super Observation " processing, i.e., will be located at m observation y of same mesh pointiPredicted value corresponding with n member(i=1,2 ..., m;K=1,2 ..., n) constitute a super observation and corresponding predicted value.Assuming that the observation of different website different moments misses It is poor mutually indepedent, yiCorresponding observation standard deviation is ri, then new to constitute observation ynew, observe standard deviation rnewAnd it is corresponding pre- ReportMeet following relationship:
Biggish matrix is stored and converted when to avoid analysis, it is contemplated that mutually indepedent between observation data, this implementation Example successively assimilates single website using sequence assimilation mode, that is, assimilates the analysis field updated after an observation as next time The ambient field of assimilation continues to assimilate.A major defect of EnKF is influenced by sampling error and mode error, and meeting exists Introduced between two independent variables false correlation cause analytical error covariance underestimate and to cause state variable to generate false Analysis increment.In order to reduce the influence for the false increment that finite aggregate generates, analysis is confined to see using localization scheme It surveys in a certain range of grid, advantage is to facilitate the weak influence for eliminating that observation generates at a distance.Optimal localization scale and same It is related to change the retention time of selection, dynamical system and the chemical species of window in an atmosphere, it is contemplated that NOxIt survives in an atmosphere Period, localization scale was set as apart from the distance for observing 2 grids of website, and PM at one day or so2.5And SO2It is set as 3 The distance of grid.EnKF another serious problem in assimilation process is filtering divergence, and it is too small and neglect to show as set deviation Depending on observation information.One method is expansion background error covariance, but due to lacking physical basis and in the area far from observation Domain introduces the linearly increasing of background error falseness, keeps the effect of this method limited.In the present embodiment, we are using every Identical disturbance all is carried out to source emission in the inverting of one step.In addition the covariance caused by sampling error in refutation process is Negative value (indicating that forecast concentration is negative with source emission amount correlation) is set as 0.Unreasonable correction value in order to prevent, it is anti-every time The discharge Dynamic gene control drilled is between 0.2~5 times of original discharge.
Observe data and error: firstly, for PM2.5, only measured value is located at 6~425 μ g/m3Within range Data just assimilate, and for SO2、NO2, the threshold range of setting is respectively 1~400 μ g/m3, 1~185 μ g/m3, it is more than this A range is considered as unreasonable is removed.Then pass through | O (t)-O (t ± 1) |≤m (t) carries out time continuity inspection Excluding outlier, wherein O (t) and O (t+1) are value of the observation at t the and t+1 moment, and m (t) is according to Jiang et a1. (2013) Experience is set as 50+0.15O (t).In addition, PM2.5、SO2、NO2Observation and predicted value deviation be greater than 150,150,300 μ g/ m3When also do not assimilate.
Observation error covariance R comes from measurement error and representive error, and usual biggish measured value has biggish measurement Error, according to the research of Schwartz et al. (2012), measurement error is defined as ε0=ermax+ermin* Π0, wherein Π0For observation, ermax is elementary error, PM2.5、SO2、NO2It is respectively set to 1.5,1.0,1.0 μ g/m3, ermin setting It is 0.0075,0.005,0.005.Representive error εrAfter parametrization, obtain,
Wherein, γ is tuning scale factor, is set as 0.5 according to sensitivity experiments, Δ x is sizing grid (36km), L The radius of influence for indicating observation website, related (being reduced to 3km) with observation position, total observation error is defined as
In short, WRF/CMAQ/EnKF assimilation Inversion System target is exactly assimilation observation data, Lai Youhua emission inventories make Discharge of the inventory of optimization closer to the current generation.EnKF emission inventories inverting is to PM2.5The influence of forecast is studied.Based on EnKF Method and WRF-CMAQ mode construct area alignment inventory inverting optimization system.Assimilate 25 days to 2015 December in 2015 first On December 30, ground SO2, NO2, PM2.5Data are observed, inverting optimizes SO2、 NO2With a PM2.5The emission inventories of discharge, grind Study carefully the regional change of inventory after inverting, and carries out pollutant concentration forecast in the inventory of optimization input WRF-CMAQ, assessment is instead Inventory after drilling is to raising PM2.5The influence of forecast.
Embodiment
Have chosen the average 2-m temperature of the meteorological observation website of CHINESE REGION 36,2-m relative humidity and 10-m wind speed pair The analogue value is assessed.Table 1, which gives, compares statistical result.All in all, the 10-m wind speed of simulation and 2-m temperature are than seeing Measured value is slightly lower, and BIAS is respectively -1m/s and 4.3%, and 2-m temperature is slightly higher, and BIAS is 0.4 DEG C.10-m wind speed 2-m temperature and The CORR of 2-m relative humidity is respectively 0.62,0.96 and 0.82, and it is good that this illustrates that the meteorological field of simulation and observation field have Consistency.
Table 1
As shown in Fig. 2, giving the time sequence of the meteorological field of the simulation and observation in two website Nanjing, East China and Shanghai The comparison of column.The simulation of temperature and relative humidity and observation more coincide on the whole, and wind speed is slightly lower, and the wind speed underestimated may Lead to positive deviation of simulation.It should be noted that since observation data precision is 1m/s, and air speed value itself is smaller, therefore can Very big statistic bias can be will cause.Research in, the wind speed of simulation is respectively in NCP and entire regional error 27.4% and 22.6%, using 36km resolution analog regional, mean wind speed has been over-evaluated more than 30%, and uses nested The wind speed of 12km resolution analog NCP is over-evaluated and is reduced to 10% hereinafter, therefore may be helped using the simulation of high-resolution nesting In the raising of wind speed simulation.On the whole, the meteorological simulation of the present embodiment does not have apparent systematic error, can reproduce meteorological field Time-varying process.
As shown in figure 3, SO2、NOXWith a PM2.5It discharges and carries out inverting by adjusting daily priori emission inventories.Such as Shown in figure, three species are in discharge scale factor (inventory/source list after the inverting) estimation analyzed every time, for SO2, four A region is respectively less than 1 in the scale factor of all window invertings, shows SO2Exist in all areas and significantly over-evaluate, In, the SCB range of decrease is the most obvious, after inverting twice, stablizes 10% or so of source list.And NOXIt is emitted on entire big Land region, the discharge of the area NCP and YRD are basically unchanged, and scale factor is respectively less than 1 after the inverting of the area YRD, and discharge decreases. For a PM2.5Inverting is discharged, although the whole continent average emission is basically unchanged, three area alignment variations are obvious.Instead After drilling, NCP discharge is overall in rising trend, and YRD and SCB discharge are on a declining curve, wherein the PM in the area YRD2.5Discharge drop It is low more obvious.There is apparent fluctuation in three species different windows in refutation process, however, fluctuation range is relatively existing There is technology smaller using assimilation window inverting CO discharge in one hour, these are the result shows that inverting discharge estimation can in different windows There can be higher uncertainty, rapidly fluctuation may be intrinsic with random model error (for example, meteorological field) and reality Discharge variation it is related.From the point of view of three regions, NCP fluctuation is bigger, may be related with complicated regional pollution.
Under limited observation, the time change for studying discharge is relatively difficult, and the purpose of the application is to determine and drop Low SO2、NOX、PM2.5The error of priori emission inventories.Therefore, after to continuous assessment 25 to 30 December in 2015, Time is carried out to all analyses in all assimilation windows and space average obtains the estimation of the discharge after inverting, is post-processed in this way Be conducive to filter out the time change of influence and the discharge estimation of Model by Stochastic Model Error.Average three species row after final inverting Estimation and scale factor are put as shown in figure 3, table 1 while giving and discharge Variant statistical before and after inverting.
For SO2Discharge, in addition to northeast and the Northwest, national most area is in different degrees of downward trend, Generally, China's Mainland averagely reduces 32.2%.NCP, YRD, SCB, Central China and the discharge of Pearl River Delta area reduce more bright Aobvious, some areas discharge is reduced to former 50% or less discharge.Wherein, it is dropped respectively after the inverting of the area NCP, YRD and SCB 40.3%, 71.6%, 84%.The prior art is it has also been found that Chongqing region SO2In the presence of significantly over-evaluating, and it is anti-by EnKF methodology Drill Chongqing region SO2Discharge, result are similar with the application.In addition, also there is a lot of other evidences to show recent years (especially After being 2010-2011) since flue gas desulfurization, SO is widely used in Chinese power plant2Discharge reduces rapidly, such as from bottom to top Emission inventories estimation and based on moonscope etc..From 2010 to 2014 year, East China SO2Vertical column density 30% or so are reduced, and a year discharge is reduced up to 50~60% from 2007-2008 to 2014.NOXIt is relatively small to discharge difference, The faint increase in central and west regions, and there is smaller reduction in eastern region, wherein the YRD range of decrease is obvious, reduces than source emission rate 30.1%.NOXDischarge is mainly derived from motor vehicle and plant emissions, based on emission inventories and satellite methods table from bottom to top Bright interior discharge nationwide before 2010 generally increases.And in western part of China, until discharge in 2013 still increases.In State's government goal is in 12 period NOXDischarge reduces by 10%, and moonscope is shown in East China NO in 2011XDischarge Reach peak, and is begun to decline during 2012.Relative to 2010, east China area NOX, VCDs is under 2014 Drop 20% or so, is consistent with the application substantially.For PM2.5Discharge, spatial distribution and PM2.5Analyze field spatial distribution more Unanimously, the reliability that the application uses inversion method is further demonstrated.On the whole, discharge increase be mainly distributed on North China and Northeast and the Northwest, and discharge reduce be mainly distributed on Henan and on the south southern area, the whole country averagely increase 12.5% PM2.5Discharge, wherein the area NCP increases 33.4%, may partially influence since winter heating etc. is uncertain. And in the area YRD and SCB, PM2.5Discharge reduces 88.3%, 41.5% respectively.Table 2
Emission inventories inverting optimizes the influence to forecast: in order to probe into the influence that inverting emission inventories simulate pollutant, Mode is re-entered into the inventory after optimization, and is compared with source list simulation.Using original during respectively studying The SO that inventory is simulated after inventory and optimization2、NO2、PM2.5Concentration and RMSE.It can be seen that being used in different research areas The SO that inventory is simulated after inverting2、NO2And PM2.5It is substantially better than source list simulation, in addition, it is different from initial fields assimilation, it uses Inventory after inverting can constantly improve model simulation results in the case where outside drain is stablized.For SO2Simulation, it is entire in State continent BIAS is by 20.2 μ g/m3It is reduced to -6.7 μ g/m3, that is to say, that error reduces 67%, in NCP, YRD and SCB, BIAS is by 36.7,46.2,96 μ g/m3It is reduced to -10.5, -1.7,1.3 μ g/m3, reduce 71.4% respectively, 96.3%, 98.6% error, hence it is evident that the system of correcting for is over-evaluated.For NO2Simulation, although continent BIAS statistics it is constant, RMSE and CORR increases, and three sub-regions source lists simulate BIAS very little, but are still effectively increased using the inventory after inverting Simulation, BIAS is respectively by 2.5,17.0, -7.7 μ g/m3Drop to -1.2,2.7, -3.1 μ g/m3.For PM2.5Simulation, four areas The BIAS in domain is respectively by -9.6, -34.1,20.9,36.5 μ g/m3Drop to -4.9, -17.6,9.7 μ g/m3、27.2μg/m3, remove Outside SCB, other three region differences all have dropped 50% or so substantially.For SCB, although inventory compares source list after inverting 41.5% is had dropped, but BIAS only reduces 25.5%, however it remains more serious over-evaluates, this may be with the meteorology of simulation And chemical deviation is in relation to (4.4 discuss).In addition, all species RMSE in four regions are declined, and CORR removes YRD NO2It has dropped other than 0.01, also there is different degrees of raising.On the whole, the application is clear by the discharge of EnKF methodology inverting List can effectively improve forecast result.
Table 3
In order to further verify influence of the source emission of assimilation to pollutant simulation is improved, three kinds of pollutants are in each city Forecast result, for SO2, existed using the analog result of source list in East China, Central China and SCB, Pearl River Delta area Significantly over-evaluate, and exist in northeast and the Northwest and significantly underestimate, using the emission inventories after assimilation, removes Gansu, mountain Outside the province, the BIAS of the regional analog result of other overwhelming majority is substantially in 10 μ g/m in west3Within.Although the SO of latter two province of inverting2 Discharge has certain increase, and there are still bigger to underestimate for the analog result after inverting, it is thus possible to also need the row of further increasing It puts.For NO2Simulation effectively corrects for central and east simulation using the emission inventories after inverting and over-evaluates, and analog result is also distributed In 10 μ g/m3Within, and for west area, however it remains certain minus deviation.PM2.5Although the Emission Optimization effectively reduces Over-evaluating for southern most area and underestimating for northern area, but big deviation is still present in Jing-jin-ji region, SCB and Xinjiang Etc. ground.Due to the complexity of Beijing-tianjin-hebei Region pollution, meteorological simulation error, chemical reaction missing and PM2.5Non-linear two It is secondary generate etc. factors all limit inverting as a result, and for Xinjiang region, anthropogenic discharge is relatively fewer, although new after inverting Boundary most area improves discharge, but is affected by the non-artificial discharge such as dust storm, therefore still deposit in Xinjiang most area Underestimate serious.On the whole, using the inventory of inverting after assimilation, it is nationwide in have 85.8% respectively, 82.5%, 85.4% city improves SO2、NO2And PM2.5Simulation.
In short, constructing set inverting assimilation system by EnKF methodology to optimize SO simultaneously2、NOX、PM2.5Emission inventories, It effectively corrects for discharge error and improves model predictions.Although but PM after the regional inverting of SCB2.5Discharge reduces 41.5%, But analog result still have it is biggish over-evaluate, error only reduces by 25.5%, and multiple city BIAS are higher than 40 μ g/m3.Pass through gas Image field assessment, discovery this area WRF 10-m wind speed simulation is seriously underestimated, in addition, to over-evaluate 37.7%, 2-m relatively wet for 2-m temperature Degree underestimates 17.5%.The wind speed underestimated is unfavorable for PM2.5Diffusion, cause the lasting accumulation of simulation, and the humidity over-evaluated promotes The moisture absorption growth of aerosol and liquid-phase conversion, the temperature over-evaluated then enhance chemical reaction rate, promote PM2.5Not medium well At.Although WRF parameter used by the application is reliable in most areas analog result, in region with a varied topography, more slightly Rough resolution ratio can may still have biggish error, this increases the uncertainty of inverting inventory to a certain extent.Cause This, needs to further increase resolution ratio to improve mode expression, especially in regions with complex terrain.
The PM of the area NCP inverting optimization2.5Discharge increases 33.4% compared with source list, which is changed to a certain extent Learn model PM2.5The incomplete influence of heterogeneous phase chemistry process of generation.
The application is to reduce a large amount of precursors to discharge error to PM of inverting2.5Influence, take while assimilating inverting PM2.5And secondary aerosol species precursor SO2、NOXThe strategy of discharge.Fig. 7 gives assimilation and does not assimilate SO2、 NO2Observation is anti- The sensitivity experiments drilled, it can be seen that there are significant differences in Middle Eastern.When the discharge of gaseous precursors object is there are when error, Also error can be generated to the concentration of the secondary inorganic matter generated in forecasting process, and then influences PM2.5Concentration, refutation process is then It can be by adjusting a PM2.5Discharge make up the error of precursor discharge, therefore very it is necessary to combine a variety of correlations of assimilation Pollutant.In addition, false increment is introduced in the analysis process to avoid observing at a distance, only in localization range scale Observation carries out inverting estimation, according to factors such as species service life and atmospheric transmission conditions.The application due to Research scale relatively Greatly, the atmospheric conditions of different zones are not quite similar, to global SO2、NO2、PM2.5The localization scale of unified setting is respectively 108, 72,108km, Fig. 8 give localization scale and are respectively set as PM when 180,108 and 180km2.5The discharge of discharge and the application Difference and scale factor.As can be seen that discharge inverting has larger sensibility to localization scale, since localization scale becomes larger, Archeus inverting less than place can also generate increment, therefore, all in all northern area discharge than the application discharge more By force, southern weaker.It is following how by model dynamical system and localization scale further combined with according to atmospheric conditions determination Different influence scale threshold values, obtaining more accurately assessment result, there is still a need for more studied and thought deeply.
The application by EnKF methodology construct set inverting assimilation system be coupled in WRF-CMAQ mode come and meanwhile it is excellent Change SO during 25 to 30 December in 20152、NOX、PM2.5Emission inventories, analysis compare the space change of inventory after inverting Change, and carries out forecast verifying with the inventory after inverting.As a result as follows:
Different pollutants discharge variation difference in different zones.In addition to northeast and the Northwest, national most area SO2Discharge is in different degrees of downward trend, wherein dropped 40.3% after the inverting of the area NCP, YRD and SCB respectively, 71.6%, 84%.For NOX, central and west regions faint increase relatively small with source list difference is discharged in discharge after inverting, And there is smaller reduction in eastern region, wherein the YRD range of decrease is obvious, reduces by 30.1% than source emission rate.For PM2.5Discharge increases Add and be mainly distributed on North China and northeast and the Northwest, and reduce be mainly distributed on Henan and on the south southern area, on the whole, The whole country averagely increases 12.5%PM2.5Discharge, wherein the area NCP increases 33.4%, and regional in YRD and SCB, PM2.5Discharge reduces 88.3%, 41.5% respectively.
In different research areas, SO is simulated using the inventory after inverting2、NO2And PM2.5It is substantially better than source list mould It is quasi-.For SO2Simulation, China's Mainland, NCP, YRD and SCB BIAS reduce 67%, 71.4%, 96.3%, 98.6%, it is bright It is aobvious to correct for system and over-evaluate.For NO2Simulation, although continent BIAS is constant, RMSE and CORR increase, in addition, three The BIAS of subregion reduces 52%, 84.1%, 59.7% respectively.In addition to SCB, other three regions are using clear after inverting Single analog PM2.5Deviation all has dropped 50% or so substantially.For SCB, although inventory has dropped than source list after inverting 41.5%, but influenced by meteorological simulation error, BIAS only reduces 25.5%, however it remains serious to over-evaluate.Generally, it removes The area YRD NO2CORR outside, analog result increases.Using the inventory of inverting after assimilation, have respectively in the whole country 85.8%, 82.5%, 85.4% city improves SO2、 NO2And PM2.5Simulation.
Utilize EnKF methodology Simultaneous Inversion SO 25 to 30 December in 20152、NOx、PM2.5Emission inventories.It analyzes excellent The Spatial Variation of the relatively primitive inventory of inventory after change.Inventory after being then based on inverting simulates PM2.5Change in time and space, point EnKF inverting emission inventories are analysed to PM2.5The influence of forecast.The result shows that the inventory after inverting can significantly improve PM2.5Forecast Effect, China's Mainland, NCP, YRD, SCB BIAS reduce by 49%, 48.4%, 53.6% and 25.5% respectively.
It is further to note that specific technical features described in the above specific embodiments, in not lance In the case where shield, it can be combined in any appropriate way.In order to avoid unnecessary repetition, the present invention is to various No further explanation will be given for possible combination.

Claims (6)

1. a kind of ground emission inventories inverting optimization method based on EnKF characterized by comprising
Step 1, building EnKF inventory Inversion System are simultaneously coupled to WRF-CMAQ modular system;
Step 2, by the inventory after daily inverting be re-applied to the same day model predictions be inverting in next day in DATA PROCESSING IN ENSEMBLE PREDICTION SYSTEM mention For initial fields, the average inventory as final optimization pass of more days inverting inventories;
Step 3, assessment inventory regional change and uncertainty;Forecast when inventory after optimization is carried out long within the assimilation period, It forecasts to compare with source list, improvement effect of the inventory after assessment optimization to forecast.
2. the ground emission inventories inverting optimization method according to claim 1 based on EnKF, which is characterized in that the structure The step of building EnKF inventory Inversion System include:
Step 11 generates source emission set sample
For source list, δ X is random perturbation field, and i is set sample number;
Step 12 considers PM2.5The influence of precursor, while to SO2、NOXDischarge carry out inverting correct, calculate the back of source emission Scape error co-variance matrix For source emission ensemble average before inverting, n is nature Number;
Step 13 is updated based on aggregate error statistics and measurement vector, the inventory after inverting by following equation:
Wherein, H Observation Operators, R are observation error covariance matrix, and K is gain matrix, determine the weight of ambient field and observation,Ensemble average for the analysis field of ith member, i.e. emission inventories after optimization, all members is optimal estimation, under Inverting provides the higher emission inventories of precision.
3. the ground emission inventories inverting optimization method based on EnKF stated according to claim 2, which is characterized in that further include step Rapid 14, inverting is carried out using average daily concentration, assimilation window is set as 1 day, to avoid model predictions diurnal variation error to the shadow of inverting It rings;Super inspection process is carried out to observation based on optimal estimation theory,
M observation y of same mesh point will be located atiPredicted value corresponding with n memberConstitute super observation and Corresponding predicted value, i=1,2 ..., m;K=1,2 ..., n, it is assumed that the observation error of different website different moments is mutually indepedent, yiCorresponding observation standard deviation is ri, then new to constitute observation ynew, observe standard deviation rnewAnd corresponding forecastUnder satisfaction Column relationship:
4. the ground emission inventories inverting optimization method according to claim 3 based on EnKF, which is characterized in that further include Step 15 successively assimilates single website using sequence assimilation mode, that is, assimilates under the analysis field conduct updated after an observation The ambient field once assimilated continues to assimilate.
5. the ground emission inventories inverting optimization method according to claim 3 based on EnKF, which is characterized in that observation misses Poor covariance R includes measurement error and representive error,
Representive error
ε0For measurement error, ε0=ermax+ermin* Π0, wherein II0For observation, ermax is elementary error;
γ is tuning scale factor, and Δ x is sizing grid, and L indicates the radius of influence of observation website, related with observation position, always Observation error be defined as
6. the ground emission inventories inverting optimization method according to claim 1 based on EnKF, which is characterized in that CMAQ mould Type system includes:
Chemical Transport module, transmission process, chemical process and infall process for simulating pollution object;
Primary condition module and boundary condition module, for providing pollutant initial fields and boundary field for CCTM;
Photochemical breakdown rate module, for calculating photochemical breakdown rate;
Meteorological-chemical interface module, is meteorologic model and the interface of CCTM, for meteorological data to be converted into what CCTM can be read Data format.
CN201910030086.5A 2019-01-14 2019-01-14 A kind of ground emission inventories inverting optimization method based on EnKF Pending CN109858686A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910030086.5A CN109858686A (en) 2019-01-14 2019-01-14 A kind of ground emission inventories inverting optimization method based on EnKF

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910030086.5A CN109858686A (en) 2019-01-14 2019-01-14 A kind of ground emission inventories inverting optimization method based on EnKF

Publications (1)

Publication Number Publication Date
CN109858686A true CN109858686A (en) 2019-06-07

Family

ID=66894600

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910030086.5A Pending CN109858686A (en) 2019-01-14 2019-01-14 A kind of ground emission inventories inverting optimization method based on EnKF

Country Status (1)

Country Link
CN (1) CN109858686A (en)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110909483A (en) * 2019-12-03 2020-03-24 河北先河环保科技股份有限公司 Point source atmospheric pollutant emission list verification method based on gridding data
CN111353634A (en) * 2020-02-12 2020-06-30 中国人民解放军国防科技大学 Sulfur dioxide emission source optimization method based on mode initial field assimilation
CN111487216A (en) * 2020-05-06 2020-08-04 北京中科锐景科技有限公司 Carbon dioxide flux inversion method and system
CN111651873A (en) * 2020-05-22 2020-09-11 上海市环境科学研究院 Urban atmosphere VOCs emission list inversion method, system, terminal and storage medium
CN111723482A (en) * 2020-06-17 2020-09-29 南京大学 Based on satellite CO2Method for observing and inverting surface carbon flux by using column concentration
CN112036038A (en) * 2020-08-31 2020-12-04 中科三清科技有限公司 Emission source emission amount adjusting and optimizing method, device, equipment and storage medium
CN114169179A (en) * 2021-12-14 2022-03-11 重庆市生态环境大数据应用中心 Air quality background concentration analysis method
CN112132336B (en) * 2020-09-22 2024-02-20 南京创蓝科技有限公司 Quarterly prediction method for PM2.5 concentration

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060229813A1 (en) * 2005-03-30 2006-10-12 Tobiska William K Ionospheric forecast system (IFS)
CN106442236A (en) * 2015-07-30 2017-02-22 中国科学院遥感与数字地球研究所 Ground PM2.5 inversion method and system based on satellite remote sensing
CN107273995A (en) * 2016-04-08 2017-10-20 株式会社日立制作所 Urban Air Pollution Methods
CN108228538A (en) * 2017-12-29 2018-06-29 航天科工智慧产业发展有限公司 SO in a kind of city2The Forecasting Methodology of concentration value
CN108519340A (en) * 2018-03-28 2018-09-11 刘诚 A kind of satellite remote sensing contamination gas scapus concentration spatial and temporal distributions inversion algorithm
CN109060617A (en) * 2018-07-06 2018-12-21 南京信息工程大学 The simulation and analysis method that different agricultural emission intensity scenes influence PM2.5 concentration

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060229813A1 (en) * 2005-03-30 2006-10-12 Tobiska William K Ionospheric forecast system (IFS)
CN106442236A (en) * 2015-07-30 2017-02-22 中国科学院遥感与数字地球研究所 Ground PM2.5 inversion method and system based on satellite remote sensing
CN107273995A (en) * 2016-04-08 2017-10-20 株式会社日立制作所 Urban Air Pollution Methods
CN108228538A (en) * 2017-12-29 2018-06-29 航天科工智慧产业发展有限公司 SO in a kind of city2The Forecasting Methodology of concentration value
CN108519340A (en) * 2018-03-28 2018-09-11 刘诚 A kind of satellite remote sensing contamination gas scapus concentration spatial and temporal distributions inversion algorithm
CN109060617A (en) * 2018-07-06 2018-12-21 南京信息工程大学 The simulation and analysis method that different agricultural emission intensity scenes influence PM2.5 concentration

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
SHUZHUANG FENG, FEI JIANG, ZIQIANG JIANG, HENGMAO WANG, ZHE CAI,: "Impact of 3DVAR assimilation of surface PM2.5 observations on PM2.5 forecasts over China during wintertime", 《ATMOSPHERIC ENVIRONMENT》 *
冯震、李兴华: "WRF-CMAQ模式在内蒙古本地化应用", 《内蒙古科技与经济》 *
吴钲、谢旻、高阳华、芦华、赵磊、高松: "利用集合均方根卡尔曼滤波反演重庆大区SO2源排放", 《环境科学研究》 *

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110909483A (en) * 2019-12-03 2020-03-24 河北先河环保科技股份有限公司 Point source atmospheric pollutant emission list verification method based on gridding data
CN110909483B (en) * 2019-12-03 2023-09-08 河北先河环保科技股份有限公司 Point source atmospheric pollutant emission list verification method based on gridding data
CN111353634B (en) * 2020-02-12 2022-09-20 中国人民解放军国防科技大学 Sulfur dioxide emission source optimization method based on mode initial field assimilation
CN111353634A (en) * 2020-02-12 2020-06-30 中国人民解放军国防科技大学 Sulfur dioxide emission source optimization method based on mode initial field assimilation
CN111487216A (en) * 2020-05-06 2020-08-04 北京中科锐景科技有限公司 Carbon dioxide flux inversion method and system
CN111487216B (en) * 2020-05-06 2023-01-24 北京中科锐景科技有限公司 Carbon dioxide flux inversion method and system
CN111651873A (en) * 2020-05-22 2020-09-11 上海市环境科学研究院 Urban atmosphere VOCs emission list inversion method, system, terminal and storage medium
CN111651873B (en) * 2020-05-22 2023-03-14 上海市环境科学研究院 Urban atmosphere VOCs emission list retrieval method, system, terminal and storage medium
CN111723482A (en) * 2020-06-17 2020-09-29 南京大学 Based on satellite CO2Method for observing and inverting surface carbon flux by using column concentration
CN111723482B (en) * 2020-06-17 2023-11-21 南京大学 Satellite CO-based 2 Method for inverting surface carbon flux by column concentration observation
CN112036038B (en) * 2020-08-31 2021-07-23 中科三清科技有限公司 Emission source emission amount adjusting and optimizing method, device, equipment and storage medium
CN112036038A (en) * 2020-08-31 2020-12-04 中科三清科技有限公司 Emission source emission amount adjusting and optimizing method, device, equipment and storage medium
CN112132336B (en) * 2020-09-22 2024-02-20 南京创蓝科技有限公司 Quarterly prediction method for PM2.5 concentration
CN114169179A (en) * 2021-12-14 2022-03-11 重庆市生态环境大数据应用中心 Air quality background concentration analysis method

Similar Documents

Publication Publication Date Title
CN109858686A (en) A kind of ground emission inventories inverting optimization method based on EnKF
Wei et al. Full-coverage mapping and spatiotemporal variations of ground-level ozone (O3) pollution from 2013 to 2020 across China
Yin et al. Surface ozone at Nam Co in the inland Tibetan Plateau: variation, synthesis comparison and regional representativeness
CN111508569B (en) Target soil property content prediction method based on soil transfer function
Monteil et al. The regional EUROpean atmospheric transport inversion COMparison, EUROCOM: first results on European wide terrestrial carbon fluxes for the period 2006–2015
Tang et al. Improvement of ozone forecast over Beijing based on ensemble Kalman filter with simultaneous adjustment of initial conditions and emissions
Itahashi et al. Seasonal source contributions of tropospheric ozone over East Asia based on CMAQ–HDDM
CN112884079A (en) Method for estimating near-surface nitrogen dioxide concentration based on Stacking integrated model
Mendoza-Dominguez et al. Estimation of emission adjustments from the application of four-dimensional data assimilation to photochemical air quality modeling
CN111487216A (en) Carbon dioxide flux inversion method and system
CN111353634B (en) Sulfur dioxide emission source optimization method based on mode initial field assimilation
Alonzo et al. Modelling the variability of the wind energy resource on monthly and seasonal timescales
Ding et al. Space-based NO x emission estimates over remote regions improved in DECSO
CN108764527B (en) Screening method for soil organic carbon library time-space dynamic prediction optimal environment variables
Trieu et al. Evaluation of summertime surface ozone in Kanto area of Japan using a semi-regional model and observation
Sun et al. Spatial-temporal prediction of air quality based on recurrent neural networks
CN114724647A (en) High spatial and temporal resolution CO2Flux inversion system and method
Feng et al. High-resolution simulation of wintertime fossil fuel CO2 in Beijing, China: Characteristics, sources, and regional transport
CN114169179A (en) Air quality background concentration analysis method
Guo et al. Construction and application of a regional kilometer-scale carbon source and sink assimilation inversion system (CCMVS-R)
CN117219183A (en) High coverage near ground NO in cloudy rain areas 2 Concentration estimation method and system
Yang et al. Estimating spatiotemporal dynamics of evapotranspiration and assessing the cause for its increase in China
Chen et al. Estimation and variation analysis of secondary inorganic aerosols across the Greater Bay Area in 2005 and 2015
Alfaro et al. CFSR-NCEP Performance for weather data forecasting in the Pernambuco Semiarid, Brazil
Mendoza-Dominguez et al. Emission strength validation using four-dimensional data assimilation: Application to primary aerosol and precursors to ozone and secondary aerosol

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