Air Quality Assessment And Control Of Emission Rates

Mathematical methods based on the adjoint model approach are given for the air-pollution estimation and control in an urban region. A simple advection–diffusion-reaction model and its adjoint are used to illustrate the application of the methods. Dual pollution concentration estimates in ecologically important zones are derived and used to develop two non-optimal strategies and one optimal strategy for controlling the emission rates of enterprises. A linear convex combination of these strategies represents a new sufficient strategy. A method for detecting the enterprises, which violate the emission rates prescribed by a control, is given. A method for determining an optimal position for a new enterprise in the region is also described.     

Industrial pollution transport. Part 1. Formulation of the problem and air pollution estimates

A pollution transport problem is formulated in a limited area. As the pollution sources we take emissions from industrial plants. Physically and mathematically suitable conditions are prescribed on the open boundaries. We show that the problem (as well as its adjoint) is well posed in the sense that a weak solution exists, is unique and depends continuously on its data. Direct and adjoint estimates of the average pollution concentration in an ecologically important zone are given, and the sensitivity of these estimates to perturbations in model parameters is analyzed.     

Industrial pollution transport. Part 2. Control of industrial emissions

Solutions of the pollution transport problem and its adjoint are used to monitor mean pollution concentration in an ecologically important zone. Four strategies of control over pollutants released into the atmosphere by industrial plants are suggested. They differ by the restrictions imposed on the emission rate of each plant. All the strategies use solutions of the adjoint transport problem and assure the fulfillment of the sanitary norm in the zone. A linear interpolation of these strategies also brings pollution level in the zone down to the sanitary norm. A method of detecting the plants violating the prescribed emission rates is also given. A simple example is given to illustrate the strategies suggested.     

Direct and Adjoint Oil Spill Estimates

Propagation of the oil spilling from a damaged oil tanker is considered in a limited sea area. The accident consequences are evaluated by means of direct and adjoint oil concentration estimates in ecologically sensitive zones. While the direct estimates are preferable to get a comprehensive idea of the oil spill impact on the whole area, the adjoint ones are useful and economical in studying the sensitivity of the oil concentration in some zones to variations in the accident site and oil spill rate from the tanker. Thanks to special boundary conditions set at the inflow and outflow parts of the open boundary, the main and adjoint oil transport problems are both well-posed according to Hadamard (1923). The estimates obtained in Skiba (1996a) are generalized to the three dimensions. Balanced, absolutely stable 2nd-order finite-difference schemes based on the splitting method are constructed for the two- and three-dimensional cases, both. The main and adjoint schemes are compatible in the sense that at every fractional step of the splitting algorithm, the one-dimensional split operators of both the schemes satisfy a discrete form of the Lagrange identity (Marchuk, 1995). In the special unforced and non-dissipative case, the schemes have two conservation laws each. Each split problem is solved by the factorization method.     

Application of Adjoint Approach to Oil Spill Problems

Adjoint method is applied to various oil spill problems. A three-dimensional model for describing the dispersion of a quasi-passive substance (a pollutant or a nutrient) and its adjoint model are considered in a limited sea region. Direct and adjoint estimates are used to get dual (equivalent) estimates of the mean concentration of the substance in important zones of the region. The role of dual estimates is illustrated with a few examples. They include such oil spill problems as the search of the most dangerous point of the oil tanker route, the oil dispersion with a climatic velocity, and the dependence of the oil concentration estimates on the oil spill rate. One more example is the application of optimal bioremediation strategy for cleaning a few zones polluted by oil. In this case, instead of oil, the model describes the dispersion of a nutrient released to marine environment. Balanced, unconditionally stable second-order finite-difference schemes based on the splitting method for the solution of the dispersion model and its adjoint are suggested. The main and adjoint difference schemes are compatible in the sense that at every fractional step of the splitting algorithm, the one-dimensional split operators of both schemes satisfy a discrete form of Lagrange identity. In the special unforced and non-dissipative case, each scheme has two conservation laws. Every split one-dimensional problem is solved by Thomas’ factorization method.     

Dual oil concentration estimates in ecologically sensitive zones

Propagation of the oil spilling from a damaged oil tanker is considered in a limited sea area. Direct and adjoint estimates of the average oil concentration in special zones are derived by using solutions of the 2-D main and adjoint oil transport problems, respectively. The dual estimates complement each other nicely in studying the oil spill consequences. While the direct estimates are preferable to get a comprehensive idea of the oil spill impact on the whole area, the adjoint ones are specially useful and economical when the accident site-dependence or/and the oil spill rate-dependence of the oil concentration is/are studied only in a few ecologically sensitive zones. Indeed, each adjoint estimate explicitly relates the average oil concentration in a zone to the oil spill rate using the adjoint solution values at the accident site. Being independent of the two parameters (the accident site and oil spill rate), the adjoint solution can be found for each zone regardless of a concrete accident and used repeatedly for various possible values of these parameters. Several examples explain how to decide between two estimates.Thanks to special boundary conditions, the main and adjoint problems are both well-posed according to Hadamard (1923). The dual estimates can be generalized to the three dimensions. The balanced, absolutely stable and compatibie main and adjoint 3-D numerical algorithms by Skiba (1993) can easily be adapted to the problem discussed here.     

Mathematical Modeling and Numerical Algorithms for Simulation of Oil Pollution

This paper deals with the mathematical modeling and algorithms for the problem of oil pollution. For solving this task, we derive the adjoint problem for the advection–diffusion equation describing the propagation of oil slick after an accident, which we call the main problem. We prove a fundamental equality between the solutions of the main and the adjoint problems. Based on this equality, we propose a novel method for the identification of the pollution source location and the accident time of oil emission. This approach is illustrated on an example for an accident in the offshore of the central part of the Vietnamese coast. Numerical simulations demonstrate the effectiveness of the proposed method. Besides, the method is verified for 1D model of substance propagation.     

A Linear-Programming-Based Strategy for Bioremediation of Oil-Polluted Marine Environments

A linear programming problem is considered with the aim to determine the optimal discharge point and the optimal discharge rate of a nutrient to be released to a marine environment polluted with oil. The objective is to minimize the total discharge of nutrient into the system provided that the concentrations of nutrient will reach critical values sufficient to eliminate oil residuals in certain affected zones through bioremediation. An initial boundary-value 3D problem for the advection–diffusion equation and its adjoint problems are considered to model, estimate, and control the dispersion of nutrient in a limited region. It is shown that the advection–diffusion problem is well posed, and its solution satisfies the mass balance equation. In each oil-polluted zone, the mean concentration of nutrient is determined by means of an integral formula in which the adjoint model solution serves as a weight function. Critical values of these mean concentrations are used as the constraints of linear programming problem. Some additional constraints are posed in order to limit not only the local discharge of the nutrient, but also the mean concentration of this substance in the whole region. Both constraints serve for environmental protection. The ability of the new method is demonstrated by numerical experiments on the remediation in oil-polluted channel using three control zones. The experiments show that the optimal discharge rate can always be got with a simple combination of step functions.     

基于数值模拟的青岛市O3快速来源解析

为了建立臭氧污染快速来源解析方法,提高解析结果的时空分辨率,以2017、2018年青岛市环境空气质量精细化管理实践为契机,利用CAMx-OSAT模型的污染源识别与追踪技术,预测解析未来时段特别是污染期间不同区域、行业的排放源对目标站点O3浓度的贡献量和贡献率。结果显示:模拟的春、夏季2个时段青岛市的O3来源主要为工业、电厂、交通、生活源,2个时段的模拟结果本地和外来源中工业源的贡献分别占62. 0%和65. 0%,交通源分别占24. 5%和16. 0%,生活源分别占8. 4%和8. 0%,电厂源分别占5. 0%和11. 0%,O3高浓度污染时,工业源的贡献上升明显,青岛市的O3浓度外来源的贡献约占40%~80%。该数值模型O3来源解析及预测方法已成功用于青岛市的空气质量精细化管理及应急管控。     

A quantitative methodology for indexing environmental sensitivity and pollution potential

A methodology for rating the suitability of sites forthe location of industrial facilities is formulated and appliedto the case of a coal-fired power plant location. The methodologycomprises two major interlinked components: the environmentalplant location indexing component, which involves theidentification, scaling and weighting of environmentalsensitivity factors; and the impact analysis component, whichinvolves the superimposition of the pollution generation impactsof an industrial facility on spatially gridded zones of variousenvironmental sensitivities. For each rectangular areal unitdefined by a square grid, the Unit Pollution Potential Index isdetermined by the severity and distribution of key environmentalsensitivity factors and the coverage of superimposed pollutanteffects as determined by contaminant fate and transport models.For any alternative site of a planned facility, the summation ofthe unit indices over the area of influence of the facilityprovides the quantitative Pollution Index, which can be used as abasis for comparison of alternative sites for planned facilities.For this paper, this methodology is applied to the hypotheticalcase of the siting of a coal-fired power plant in thenortheastern region of the United States, in which threealternative sites are considered. The three sites: A, B and Cyielded indices of 47.83, 47.91 and 47.6, respectively, indicatingthat site C is the most suitable for siting the power plant.     

沱江流域主要污染负荷预测及时空分布特征

耦合社会-经济因子探究工业点源和生活污染源污染负荷未来变化趋势,可为优化水环境规划和管理方案提供理论依据。选取沱江流域为研究区域,采用经济增长预测法、工业点源传统统计法、人口趋势灰色模型预测法和排污系数法分别计算了2020-2025年该区域28个县(市、区)的工业GDP值,工业点源废水排放量及主要污染负荷(COD、NH₃-N、TN、TP),农村与城镇人口及生活污染源的主要污染负荷,并利用ArcGIS技术探究了工业点源和生活污染源主要污染负荷空间分布特征。结果表明:2020-2025年,工业GDP值总体呈逐年增加趋势,而工业废水排放量总体呈逐年减少趋势,预计到2025年,流域工业GDP值将增加至2.52×10¹²元,而工业废水排放量将减少至0.64×10⁸ t。工业点源主要污染负荷表现为COD>NH₃-N>TN>TP。沱江流域总人口数与生活污染源污染负荷呈逐年增加趋势,其中城镇人口与生活污染源污染负荷呈逐年增加趋势,农村人口与生活污染源污染负荷呈逐年减少趋势,且城镇人口及生活污染源污染负荷增加量大于农村人口及生活污染源污染负荷减少量。城镇、农村生活污染源的主要污染负荷表现为COD>NH₃-N>TN>TP。工业点源和生活污染源主要污染负荷在空间上存在高度异质性。2025年,来自工业点源的主要污染负荷均呈上游较少,中、下游较多的特征;来自城市生活污染源的主要污染负荷均呈中、上游较多,下游较少的特征;来自农村生活污染源的主要污染负荷均呈中游较多,上、下游较少的特征。笔者提出耦合社会-经济因子预测流域污染负荷的方法可以推广到其他与社会经济指标相关联的流域工业点源、生活污染源污染负荷的预测研究中,以期为未来流域水环境管理与治理提供科学参考。     

Source Release-Rate Estimation of Atmospheric Pollution from a Non-Steady Point Source at a Known Location

The goal is to build up an inverse model capable of finding the release history of atmospheric pollution by using measured gas concentration data at just one location on the ground and identify the factors which affects the accuracy of the model predictions. The problem involves a non-steady point source of pollution at a known location in the atmosphere. This problem of finding the release rate is an ill-posed inverse problem and its solution is extremely sensitive to errors in the measurement data. Special regularisation methods, which stabilise the process of the solution, must be used to solve the problem. The method described in this paper is based on linear least-squares regression and Tikhonov regularisation, coupled with the solution of an advection-diffusion equation for a non-steady point source. The accuracy of the method is examined by imposing normally-distributed relative noise into the concentration data generated by the forward model as well as some real experimental data.     

Preparation and Validation of Gridded Emission Inventory of Criteria Air Pollutants and Identification of Emission Hotspots for Megacity Delhi

Delhi is one of the many megacities struggling with punishing levels of pollution from industrial, residential, and transportation sources. Over the years, pollution abatement in Delhi has become an important constituent of state policies. In the past one decade a lot of policies and regulations have been implemented which have had a noticeable effect on pollution levels. In this context, air quality models provide a powerful tool to study the impact of development plans on the expected air pollution levels and thus aid the regulating and planning authorities in decision-making process. In air quality modeling, emissions in the modeling domain at regular interval are one of the most important inputs. From the annual emission data of over a decade (1990–2000), emission inventory is prepared for the megacity Delhi. Four criteria pollutants namely, CO, SO₂, PM, and NO _x are considered and a gridded emission inventory over Delhi has been prepared taking into account land use pattern, population density, traffic density, industrial areas, etc. A top down approach is used for this purpose. Emission isopleths are drawn and annual emission patterns are discussed mainly for the years 1990, 1996 and 2000. Primary and secondary areas of emission hotspots are identified and emission variations discussed during the study period. Validation of estimated values is desired from the available data. There is a direct relationship of pollution levels and emission strength in a given area. Hence, an attempt has been made to validate the emission inventory for all criteria pollutants by analyzing emissions in various sampling zones with the ambient pollution levels. For validation purpose, the geographical region encompassing the study area (Delhi) has been divided into seven emission zones as per the air quality monitoring stations using Voronoi polygon concept. Dispersion modeling is also used for continuous elevated sources to have the contributing emissions at the ground level to facilitate validation. A good correlation between emission estimates and concentration has been found. Correlation coefficient of 0.82, 0.77, 0.58 and 0.68 for CO, SO₂, PM and NO _x respectively shows a reasonably satisfactory performance of the present estimates.     

环境统计与排污许可衔接的工业源排放核算研究

环境统计面临依据排污许可数据估算工业污染排放总量的挑战.以某市4个行业6类污染物的排放数据为案例,讨论了在非重点源污染排放量估算问题上衔接排污许可和环境统计的可行性和存在的问题.研究结果表明:排污许可重点管理单位与环境统计重点调查单位存在显著差异,不能采用比率估算方法估算工业非重点源排污量;排污许可简化和登记管理单位排...     

A Simple Model for Urban Background Pollution

A simple urban background pollution model is presented. Contributions from the individual area sources, subdivided into a grid net of a resolution of 2km × 2km, are integrated along the wind direction path assuming linear dispersion with the distance to the receptor point. Horizontal dispersion is accounted for by averaging the calculated concentrations over a certain, wind speed dependent wind direction sector, centred on the average wind direction. Formation of the nitrogen dioxide due to oxidation of nitrogen monoxide by ozone is calculated using a simple chemical model based on assumption of a photochemical equilibrium on the time scale of the pollution transport across the city area. The rate of entrainment of fresh rural ozone is governed by this time scale. The model is suitable for calculations of urban background when the dominating source is the road traffic. For this source the emissions take place at ground level, and a good approximation is to treat the emissions as area sources, but with an initial vertical dispersion determined by the height of the buildings.     

陕西省重金属污染特征分析

以陕西省重金属工业污染源为例,利用等标污染负荷法和单位产值等标污染负荷法,从污染物种类、时间变化、空间分布、行业排放4个方面,分析了陕西省工业污染源重金属的排放特征,重点比较了1991~2009年近20年以来重金属排放特征的变化。结果表明,近20年陕西省工业污染源重金属等标污染负荷大小排序为铅>砷>镉>六价铬>汞; 1991~2009年,重金属排放总量呈较大下降趋势,单位产值重金属排放总量呈显著下降趋势; 1991年重金属等标污染负荷最大的行政区和流域分别为渭南市和渭河流域,2009年则转移至西安市和嘉陵江流域; 1991年和2009年,污染负荷最大的行业均为有色金属矿采选业、有色金属冶炼及压延加工业,其铅排放应作为整治的重点。     

Magnetic mapping of fly-ash pollution and heavy metals from soil samples around a point source in a dry tropical environment

The Singrauli region in the southeastern part of Uttar Pradesh, India is one of the most polluted industrial sites of Asia. It encompasses 11 open cast coalmines and six thermal power stations that generate about 7,500 MW (about 10% of India’s installed generation capacity) electricity. Thermal power plants represent the main source of pollution in this region, emitting six million tonnes of fly-ash per annum. Fly-ash is deposited on soils over a large area surrounding thermal power plants. Fly-ashes have high surface concentrations of several toxic elements (heavy metals) and high atmospheric mobility. Fly ash is produced through high-temperature combustion of fossil fuel rich in ferromagnetic minerals. These contaminants can be identified using rock-magnetic methods. Magnetic susceptibility is directly linked to the concentration of ferromagnetic minerals, primarily high values of magnetite. In this study, magnetic susceptibility of top soil samples collected from surrounding areas of a bituminous-coal-fired power plant were measured to identify areas of high emission levels and to chart the spatial distribution of airborne solid particles. Sites close to the power plant have shown higher values of susceptibility that decreases with increasing distance from the source. A significant correlation between magnetic susceptibility and heavy metal content in soils is found. A comparison of the spatial distribution of magnetic susceptibility with heavy-metal concentrations in soil samples suggests that magnetic measurements can be used as a rapid and inexpensive method for proxy mapping of air borne pollution due to industrial activity.     

Waste management modeling by MARKAL model: A case study for Basilicata Region

Air pollution is one of the most pressing environmental problems which affects likewise urban, industrial and rural areas. Environmental planners, regulators and decision makers need reliable, scientifically based tools to find out strategies for controlling air pollution in a cost-effective way, taking into account the whole productive system. In this framework the basic elements of energy-environmental planning have to be extended to include also waste processing technologies amongst the usually considered pollution sources. Bottom-up optimizing models, based on linear programming techniques and customized for specific cases, represent a powerful tool in energy-environmental management. This paper focuses on the integrated modeling of material flows and energy system performed on a local scale case study (Basilicata Region, Southern Italy) using the linear programming model IEA-MARKAL. We have evaluated the feasibility of the model in representing the waste management system to estimate the environmental impact of the waste processing technologies in the context of the whole productive system. A sensitivity analysis has been carried out to emphasize the connections between tariffs, waste disposal technologies assessment and atmospheric emissions. This revised version was published online in July 2006 with corrections to the Cover Date.     

基于空间计量模型的四川省大气污染特征及影响因素分析

四川省细颗粒物污染问题越来越受到重视,为有效识别四川省大气污染空间分布情况及影响因素,利用2015—2020年PM_2.5监测数据,综合分析了四川省大气污染时空分布特征,选取同期气象要素观测数据和社会经济数据,区分出全省及省内不同经济区大气污染的主要影响因素。结果表明:2015—2020年四川省的PM_2.5浓度逐年下降,日变化存在明显的双峰双谷趋势,且具有明显季节性特征,空间分布上具有明显的空间聚集现象;PM_2.5的排放与人口密度、经济水平和气温呈显著正相关,与城市绿化、风速呈显著负相关。该研究为经济增长方式优化、产业结构调整、绿化水平改善等提供了政策建议,可为污染防治、优化人居环境提供参考。     

Siting analyses for water quality sampling in a catchment

Pollution loads discharged from upstream development or human activities significantly degrade the water quality of a reservoir. The design of an appropriate water quality sampling network is therefore important for detecting potential pollution events and monitoring pollution trends. However, under a limited budgetary constraint, how to site an appropriate number of sampling stations is a challenging task. A previous study proposed a method applying the simulated annealing algorithm to design the sampling network based on three cost factors including the number of reaches, bank length, and subcatchment area. However, these factors are not directly related to the distribution of possible pollution. Thus, this study modified the method by considering three additional factors, i.e. total phosphorus, nitrogen, and sediment loads. The larger the possible load, the higher the probability of a pollution event may occur. The study area was the Derchi reservoir catchment in Taiwan. Pollution loads were derived from the AGNPS model with rainfall intensity estimated using the Thiessen method. Analyses for a network with various numbers of sampling sites were implemented. The results obtained based on varied cost factors were compared and discussed. With the three additional factors, the chosen sampling network is expected to properly detect pollution events and monitor pollution distribution and temporal trends.