Storage of CO2 in saline aquifers: Effects of gravity, viscous, and capillary forces on amount and timing of trapping

CO₂ can be effectively immobilized during CO₂ injection into saline aquifers by residual trapping – also known as capillary trapping – a process resulting from capillary snap-off of isolated CO₂ bubbles. Simulations of CO₂ injection were performed to investigate the interplay of viscous and gravity forces and capillary trapping of CO₂. Results of those simulations show that gas injection processes in which gravitational forces are weak compared to viscous forces (low gravity number N_gv) trap significantly more CO₂ than do flows with strong gravitational forces relative to the viscous forces (high N_gv). The results also indicate that over a wide range of gravity numbers (N_gv), significant fractions of the trapping of CO₂ can occur relatively quickly. The amount of CO₂ that is trapped after injection ceases is demonstrated to correlate with N_gv. For some simulated displacements, effects of capillary pressure and aquifer dip angle on the amount and the rate of trapping are reported. Trapping increases when effects of capillary pressure and aquifer inclination are included in the model. Finally we show that injection schemes such as alternating injection of brine and CO₂ or brine injection after CO₂ injection can also enhance the trapping behavior.     

An intergrated waste survey and environmental effects of COGIDO, a bleached pulp and paper mill in Vietnam, on the receiving waterbody

An integrated investigation on wastewater characterization and the environmental effects from the COGIDO pulp and paper mill in Bien Hoa Industrial Estate, Vietnam, a chlorine bleached soda integrated pulp and paper mill operating without a chemical recovery system, on the receiving water body was conducted during the rainy and dry seasons in 1993 and 1995. The pollution load from the mill was very high in terms of BOD, COD and SS (COD_m: 58.7 t/d; BOD: 33.3 t/d and SS: 25.1 t/d). The effluent toxicity was determined using four toxicity tests: the green micro-alga, Selenastrum capricornutum, Microtox (marine bacteria: Photobacterium phosphoreum), the duckweed, Lemna aequinoctialis, and fish (silver barb: Puntius gonionotus, and Tilapia: Tilapia nilotica). Selenastrum capricornutum was the most sensitive among the tested organisms. The mill toxicity emission rate (TER) was as high as 338 610 (Selenastrum test). The bleaching-pulp and semi-chemical pulp plants which contributed the largest pollution load to the total COGIDO effluent, therefore, were targeted for abatement measures. Physico-chemical parameters as well as qualitative and quantitative aquatic organism composition for the river water were established. The BOD₅ and COD values exceeded the potable surface water standard by a factor of 2 to 4. The species diversity and abundance of the phytoplankton, zooplankton and zoobenthos were found to be lower (20–40%) than that of unpolluted rivers in Vietnam, whereas pollution-indicator species increased up to four times during the dry season 1995.     

Design of a representative and cost-effective sampling programs for industrial wastewater with examples from bleached kraft pulp and paper mills

The diagram constructed for selection of sampling methods indicates that, for a systematic error, E ≤ 13%, grab sampling (GS) may be used to characterize effluents with variation coefficient of flow ≤ 120% and of contaminant concentrations ≤ 10%. For the whole studied range of variation coefficient of contaminant concentrations (2–82%), time-proportional compositing (TC) method may be applied with E ≤ 10% for effluents characterized by variation coefficients in flow < 90%. The more complicated flow-proportional compositing (FC) method is required only for effluents with larger variation coefficients or to produce more precise results. The diagram constructed for selection of sampling frequencies indicates that sequential sampling at intervals of approximately 60 min may be applied with E ≤ 10% for effluents characterized by variation coefficients ≤ 30%. Practical application of the diagrams, constructed using normal series, was checked against monitoring data from two pulp and paper mills in Vietnam. The two diagrams provided results on sampling methods and frequencies in good agreement with those obtained from actual monitoring data with percentages of agreement cases of 80 and 75%, respectively. The approach was applied in design of a monitoring program at the Bai Bang integrated pulp and paper mill in Vietnam.     

Long-term variations of CO2 trapped in different mechanisms in deep saline formations: A case study of the Songliao Basin,China

The geological storage of CO₂ in deep saline formations is increasing seen as a viable strategy to reduce the release of greenhouse gases to the atmosphere. There are numerous sedimentary basins in China, in which a number of suitable CO₂ geologic reservoirs are potentially available. To identify the multi-phase processes, geochemical changes and mineral alteration, and CO₂ trapping mechanisms after CO₂ injection, reactive geochemical transport simulations using a simple 2D model were performed. Mineralogical composition and water chemistry from a deep saline formation of Songliao Basin were used. Results indicate that different storage forms of CO₂ vary with time. In the CO₂ injection period, a large amount of CO₂ remains as a free supercritical phase (gas trapping), and the amount dissolved in the formation water (solubility trapping) gradually increases. Later, gas trapping decrease, solubility trapping increases significantly due to the migration and diffusion of CO₂ plume and the convective mixing between CO₂-saturated water and unsaturated water, and the amount trapped by carbonate minerals increases gradually with time. The residual CO₂ gas keeps dissolving into groundwater and precipitating carbonate minerals. For the Songliao Basin sandstone, variations in the reaction rate and abundance of chlorite, and plagioclase composition affect significantly the estimates of mineral alteration and CO₂ storage in different trapping mechanisms. The effect of vertical permeability and residual gas saturation on the overall storage is smaller compared to the geochemical factors. However, they can affect the spatial distribution of the injected CO₂ in the formations. The CO₂ mineral trapping capacity could be in the order of 10 kg/m³ medium for the Songliao Basin sandstone, and may be higher depending on the composition of primary aluminosilicate minerals especially the content of Ca, Mg, and Fe.     

Measurements of non-wetting phase trapping applied to carbon dioxide storage

We measure the trapped non-wetting phase saturation as a function of the initial saturation in sand packs. The application of the work is for CO₂ storage in aquifers where capillary trapping is a rapid and effective mechanism to render injected CO₂ immobile. The CO₂ is injected into the formation followed by chase brine injection, or natural groundwater flow, which displaces and traps CO₂ on the pore scale as a residual immobile phase. Current models to predict the amount of trapping are based on experiments in consolidated media, while CO₂ may be stored in relatively shallow, poorly consolidated systems. We use analogue fluids at ambient conditions. The trapped saturation initially rises linearly with initial saturation to a value of approximately 0.13 for oil/water systems and 0.14 for gas/water systems. There then follows a region where the residual saturation is constant with further increases in initial saturation. This behaviour is not predicted by the traditional literature trapping models, but is physically consistent with unconsolidated media where most of the larger pores can easily be invaded at relatively low saturation and there is, overall, relatively little trapping. A good match to our experimental data was achieved with the trapping model proposed by Aissaoui.     

Numerical simulations on the influence of matrix diffusion to carbon sequestration in double porosity fissured aquifers

The double porosity model for fissured rocks, such as limestones and dolomites, has some features that may be relevant for carbon sequestration. Numerical simulations were conducted to study the influence of matrix diffusion on the trapping mechanisms relevant for the long-term fate of CO₂ injected in fissured rocks. The simulations show that, due to molecular diffusion of CO₂ into the rock matrix, dissolution trapping and hydrodynamic trapping are more effective in double porosity aquifers than in an equivalent porous media. Mineral trapping, although assessed indirectly, is also probably more relevant in double porosity aquifers due to the larger contact surface and longer contact time between dissolved CO₂ and rock minerals. However, stratigraphic/structural trapping is less efficient in double porosity media, because at short times CO₂ is stored only in the fissures, requiring large aquifer volumes and increasing the risk associated to the occurrence of imperfections in the cap-rock through which leakage can occur. This increased risk is also a reality when considering storage in aquifers with a regional flow gradient, since the CO₂ free-phase will move faster due to the higher flow velocities in fissured media and discharge zones may be reached sooner.     

Natural resources and 'gradual' reform in Uzbekistan and Turkmenistan

Among low‐income transition reformers, natural resource rents are an important initial condition that helps explain choice of reform strategy. Resource‐rich Uzbekistan and Turkmenistan and resource‐poor China and Vietnam all claim to pursue gradual reform, but their strategies differ. In China and Vietnam, low resource rents have nurtured developmental political conditions and encouraged efficient resource use, which initially promoted agriculture as a dynamic market sector, capable of absorbing labour from the lagging state sector. In contrast, the scale and ease of natural resource rent extraction in the Central Asian countries has consolidated authoritarian governments that postpone reform. Despite high energy rents, Uzbekistan and Turkmenistan still extract agricultural rents in ways that repress farm incentives, perpetuate environmental degradation and liquidate irrigation assets. Uzbekistan uses its rents to subsidize a manufacturing sector, that is neither dynamic nor competitive. As its dynamic sector, Turkmenistan promotes natural gas exports that depend on volatile markets. Resource‐driven development models suggest that reform is required in both countries to avert a growth collapse. Turkmenistan's large energy rent‐stream may postpone a collapse for some years, but Uzbekistan's position is already precarious: it has run down its rural infrastructure and accumulated sizeable foreign debts and will require external assistance to recover from a growth collapse. Such assistance should be made conditional on accelerated economic and political reform.     

Weed manipulation for insect pest management in corn

Populations of insect pests and associated predaceous arthropods were sampled by direct observation and other relative methods in simple and diversified corn habitats at two sites in north Florida during 1978 and 1979. Through various cultural manipulations, characteristic weed communities were established selectively in alternate rows within corn plots.Fall armyworm (Spodoptera frugiperda J. E. Smith) incidence was consistently higher in the weed-free habitats than in the corn habitats containing natural weed complexes or selected weed associations. Corn earworm (Heliothis zea Boddie) damage was similar in all weed-free and weedy treatments, suggesting that this insect is not affected greatly by weed diversity. Only the diversification of corn with a strip of soybean significantly reduced corn earworm damage.In one site, distance between plots was reduced. Because predators moved freely between habitats, it was difficult to identify between-treatment differences in the composition of predator communities. In the other site, increased distances between plots minimized such migrations, resulting in greater population densities and diversity of common foliage insect predators in the weed-manipulated corn systems than in the weed-free plots. Trophic relationships in the weedy habitats were more complex than food webs in monocultures.Predator diversity (measured as mean number of species per area) and predator density was higher in com plots surrounded by mature, complex vegetation than at those surrounded by annual crops. This suggests that diverse adjacent areas to crops provide refuge for predators, thus acting as colonization sources.     

Predicting δCDIC dynamics in CCS: A scheme based on a review of inorganic carbon chemistry under elevated pressures and temperatures

Stable carbon isotopes are important tools to assess potential storage sites for CO₂, as they allow the quantification of ionic trapping via isotope mass balances. In deep geological formations high p/T conditions need to be considered, because CO₂ dissolution, equilibrium constants and isotope fractionation of dissolved inorganic carbon (DIC) depend on temperature, pressure and solute composition. After reviewing different approaches to account for these dependencies, an expanded scheme is presented for speciation and carbon isotope fractionation of DIC and dissolution of CaCO₃ for pCO₂ up to 100 bar, pH down to 3 and temperatures of up to 200 °C. The scheme evaluates the influence of respective parameters on isotope ratios during CO₂ sequestration. The pCO₂ and pH are the dominant controlling factors in the DIC/δ¹³C/pH system. The fugacity of CO₂ has major impact on DIC concentrations at temperatures below 100 °C at high pCO₂. Temperature dependency of activities and equilibrium dominates at temperatures above 100 °C. Isotope ratios of DIC are expected to be about 1–2‰ more depleted in ¹³C compared to the free CO₂ at pCO₂ values above 10 bar. This depletion is controlled by carbon isotope fractionation between CO₂ and H₂CO₃* which is the dominant species of DIC at the resulting pH below 5.     

Comparison of Indexes for Prioritizing Placement of Water Quality Buffers in Agricultural Watersheds1

Dosskey, Michael G. and Zeyuan Qiu, 2011. Comparison of Indexes for Prioritizing Placement of Water Quality Buffers in Agricultural Watersheds. Journal of the American Water Resources Association (JAWRA) 47(4):662‐671. DOI: 10.1111/j.1752‐1688.2011.00532.x Abstract: Five physically based, spatially distributed, empirical indexes were compared for the degree to which they identified the same or different locations in watersheds where vegetative buffers would function better for reducing agricultural nonpoint source pollution. All five indexes were calculated on a 10 m × 10 m digital elevation grid on agricultural land in the 144‐km² Neshanic River watershed in New Jersey. The indexes included the topography‐based Wetness Index (WI) and Topographic Index (TI) and three soil survey‐based indexes (sediment trapping efficiency [STE], water trapping efficiency [WTE], and groundwater interaction [GI]). Results showed that each index associated higher pollution risk and mitigation potential to a different part of the landscape. The WI and TI identified swales and riparian areas where runoff converges, whereas STE and WTE identified upland sites. The STE and WTE lack the fine scale of slope resolution and the accounting for convergent runoff patterns that can be important for properly locating buffers in some watersheds. The GI index indicates the existence of a shallow water table but the correspondence with WI‐ and TI‐identified sites was only modest. For watersheds where pollutant loading is generated by both saturation‐excess (emphasized by TI and WI) and infiltration‐excess processes (emphasized by STE and WTE), the indexes could be complementary. However, techniques would be needed for properly apportioning priority among sites identified by each index.     

Dealing with uncertainty: a case study of controlling insect populations in natural ecosystems

This paper analyses the political process for handling an outbreak of an insect causing human allergic reactions. In the southern part of the Swedish island Gotland in the Baltic Sea, an outbreak of the northern pine processionary moth, Thaumetopoea pinivora, has occurred. With regard to the human nuisance and long-term effects on the tourist industry, demands have been raised for intervention to reduce and control the insect population. At the same time, there have been warnings against treating the insect population because there are knowledge gaps concerning the wider ecological consequences, including effects on biodiversity. This paper analyses the political process and its problem-solving efforts. Of particular interest is how it dealt with the issue of uncertainty. It was found that the uncertainty associated with the issue hindered the development of a shared understanding of the problem and a possible solution. There seems to be a growing need in society to develop the institutional capacity to handle complex issues that cross different sectors, regulatory frameworks and policy targets.     

Small farm ponds: overlooked features with important impacts on watershed sediment transport

Despite their size, small farm ponds are important features in many landscapes. Yet hydrographical databases often fail to capture these ponds, and their impacts on watershed processes remain unclear. For a 230‐km² portion of central Texas, United States (U.S.), we created a historical inventory of ponds and quantified the accuracy of automated detection methods under varying drought conditions. In addition, we documented pond dredging/enlargement events and identified sites that had been abandoned. We also analyzed sediment cores from downstream reservoirs to track changes in watershed sediment transport. Over 75 years, pond densities increased more than 350% — to among the highest documented in the U.S. — and the ability of automated methods to detect these ponds was highly dependent on drought severity (R² = 0.96). Approximately 5% of ponds present in the 1950s were no longer present in 2012, while 33% were dredged between 1937 and 2012. Downstream reservoir sedimentation has decreased by an average of 55% as ponds have increased in number. These findings suggest that small ponds and the maintenance of trapping efficiency have large‐scale impacts on sediment dynamics. Accurately accounting for these storage effects is vital to water resource planning efforts.     

Barriers and enablers to climate change adaptation in hierarchical governance systems: the case of Vietnam

Governments fulfil important roles in increasing the adaptive capacity of local communities to respond to climate change impacts, particularly in developing countries. Existing studies on how governments enable and constrain the ways in which local level communities learn and build their adaptive capacity, however, generally adopt network or market-oriented types of governance. However, the most vulnerable regions to climate change impact in the world are generally governed through hierarchical policy systems. This research aims to understand how the hierarchical policy system in Vietnam creates enables and/or constrains the policy capacity of policy actors to contribute to effective climate change adaptation. We conducted interviews (n?=?26) with key actors at multiple levels of government. Our findings show the importance of clear legal institutions, available financing for implementing policies, and the training of governmental staff, particularly at district and commune levels where the policy capacities are generally too low to deal with climate change impacts. We conclude that any efforts to support local actors (i.e. smallholder farmers) should include investments in policy capacity to ensure uptake and upscaling of adaptation actions more broadly.     

CO2 storage capacity estimation: Methodology and gaps

Implementation of CO₂ capture and geological storage (CCGS) technology at the scale needed to achieve a significant and meaningful reduction in CO₂ emissions requires knowledge of the available CO₂ storage capacity. CO₂ storage capacity assessments may be conducted at various scales—in decreasing order of size and increasing order of resolution: country, basin, regional, local and site-specific. Estimation of the CO₂ storage capacity in depleted oil and gas reservoirs is straightforward and is based on recoverable reserves, reservoir properties and in situ CO₂ characteristics. In the case of CO₂-EOR, the CO₂ storage capacity can be roughly evaluated on the basis of worldwide field experience or more accurately through numerical simulations. Determination of the theoretical CO₂ storage capacity in coal beds is based on coal thickness and CO₂ adsorption isotherms, and recovery and completion factors. Evaluation of the CO₂ storage capacity in deep saline aquifers is very complex because four trapping mechanisms that act at different rates are involved and, at times, all mechanisms may be operating simultaneously. The level of detail and resolution required in the data make reliable and accurate estimation of CO₂ storage capacity in deep saline aquifers practical only at the local and site-specific scales. This paper follows a previous one on issues and development of standards for CO₂ storage capacity estimation, and provides a clear set of definitions and methodologies for the assessment of CO₂ storage capacity in geological media. Notwithstanding the defined methodologies suggested for estimating CO₂ storage capacity, major challenges lie ahead because of lack of data, particularly for coal beds and deep saline aquifers, lack of knowledge about the coefficients that reduce storage capacity from theoretical to effective and to practical, and lack of knowledge about the interplay between various trapping mechanisms at work in deep saline aquifers.     

Evaluation of the efficiency of some sediment trapping methods after a Mediterranean forest fire

Forest fires are common in Mediterranean environments and may become increasingly more frequent as the climate changes. Destruction of the forest cover and litter layer leads to greater overland flow and increased erosion rates. The greatest risk occurs during the first rainstorms following a major fire, so local authorities must act quickly to put erosion control methods in place in order to avoid excessive post-fire sediment loads in river channels. Deciding on which methods to use requires accurate knowledge of their impact on sediment load and an estimate of their cost efficiency. The objective of this study was to evaluate the efficiency of Log Debris Dams (LDDs) and a sedimentation basin for their effectiveness in trapping sediments. Paired sub-catchments were studied to quantify the amount of sediments trapped in stream channels by a series of LDDs and a sedimentation basin. Cost efficiency was evaluated for each of the measures as a function of the cost per unit volume of sediments trapped. In addition, grain size analyses were performed to characterise the nature of the sediments trapped. A third sediment trapping method, Log Erosion Barriers (LEBs) was evaluated more superficially than the first two and conclusions regarding this method are tentative. LDDs trapped a mean volume of 1.57 m³ per unit (median = 1.28 m³); mean LDD height was 105.4 cm (std. dev. = 21.9 cm), and mean height of trapped sediments was only 50.0 cm (std. dev. = 22.9 cm), showing that the traps were only half filled. Sediment height was limited by the presence of gaps between logs or branches that allowed runoff to flow through. Comparison of the textural characteristics of slope and trapped sediments showed distinct sorting: particles greater than 20 mm were not mobilised from the slopes during the study period, sediments in the medium to coarse sand size fractions were trapped preferentially by the LDDs, and sediments in the sedimentation basin were enriched by clay and silt sized (<0.050 mm) particles as coarser sediments were trapped upstream by the LDDs. Cost efficiency of LDDs was estimated at about 143 € m⁻³ for the LDDs and 217 € m⁻³ for the sedimentation basin at the time of sampling. The LDDs are therefore a cost effective method of trapping sediments, but they can only be used when pine trees or straight-trunked trees are locally available. In this case, they should be combined with LEBs, which had a cost efficiency estimated at about 250 € m⁻³. Installation of the LEBs had not been optimised and they have the advantage of trapping sediments on the slopes where they can continue to play an ecological role, so this method can give better results with more care. Sedimentation basins can be emptied if necessary and are useful in areas where pine trees are not available and where the site can be secured.     

NONLINEAR PROGRAMMING IN RIVER BASIN MODELING1

ABSTRACT: This paper explores the use of nonlinear programming in river basin water quality modelling. Applications recently reported in the literature, along with the author's experience with nonlinear programming, are reviewed. Results obtained using nonlinear programming are compared with the results obtained by other researchers using linear and dynamic programming to solve river basin water quality optimization problems. These water quality models have objective functions with continuous first partial derivatives, several inequality and variable bound constraints, and are of the form: minizie Σ^j=n_j=1Y_j(X_j) subject to Σ^j=n_j=1a_ijX_jb_i, i=1,2, …, m c_jX_jd_j, j= 1,2, …, n The variable X_i is the maximum allowable ratio of the BOD (biochemical oxygen demand) of the effluent outflow to the BOD of the wastewater inflow for treatment plant j, in the range c_j to d_j. The a_ijd and b_i are constants in the DO (dissolved oxygen) and BOD constraints. The resuks show, given certain assumptions about the data, that nonlinear programming is a better solution method for these problems than is either linear programming or dynamic programming.     

Volatilization characteristics of organic solutes in stirred solution

The effects of turbulence intensity (velocity gradient, G (s⁻¹)), Henry's law constant (H), and molecular weight (M) on the volatilization rates of organic compounds are examined using changes in the mass transfer coefficients (K_OL (cm/min)) under specific liquid-mixing intensities. The selected compounds were divided into three groups according to their H values (mole in gas/mole in liquid, dimensionless), which ranged from 10² to 10⁻⁵. The relationship of the K_OL relative to G, H and M was obtained via multiple regression. The obtained values of these parameters indicate that the primary factor affecting the K_OL values of the high H compounds is their M values. The effects of the H values on the K_OL values of the high H compounds can be neglected. On the other hand, the H value is the major factor determining the K_OL values of the low H compounds. The changes in the K_OL values of the different H compounds exhibit different profiles as the liquid-mixing intensity increases. The M and H values of middle H compounds possibly affect their K_OL values. The effects of the liquid-mixing intensity on the K_OL values of the organic compounds increase with increasing H values. The variation in the K_OL values might be a result of the concentration of the organic compounds at the interface between the liquid and gas films. The empirical relationship between K_OL and some selected parameters, G, H and M, is examined in this study. The obtained results can help to estimate volatilization loss of organic solutes in wastewater treatment facilities.     

The determination and matching analysis of pinch point temperature difference in evaporator and condenser of organic rankine cycle for mixed working fluid

Exergo-economic analysis of the pinch point temperature difference (PPTD) in both evaporator and condenser of sub-critical organic Rankine cycle system (ORCs) are performed based on the first and second laws of thermodynamics. Taking mixture R13I1/R601a as a working fluid and the annual total cost per net output power Z as exergo-economic performance evaluation criterion, the effects of PPTD in evaporator ΔT_e, and the PPTD ratio of condenser to evaporator y, on the exergo-economic performance of ORCs are analyzed. Moreover, how some other parameters influence the optimal PPTD in evaporator ΔT_e,opt and the optimal PPTD ratio of condenser to evaporator y_opt are also discussed. It has been found that the exergo-economic performance of ORCs is remarkably influenced by ΔT_e and y, and there exists ΔT_e,opt and y_opt. In addition, ΔT_e,opt and y_opt are affected by heat transfer coefficient ratio of condenser to evaporator ß, the temperature of working fluid at dew point in condenser T_1a, and composition of R13I1/R601a: larger ß and T_1a lead to lower ΔT_e,opt and y_opt; by contraries, larger mass fraction of R13I1 makes ΔT_e,opt and y_opt increase, and y_opt increases linearly. The effects of the temperature of working fluid at bubble point in evaporator T_3a, mass flow rate of exhaust flue gas m_g, and inlet temperature of exhaust flue gas T_gi on ΔT_e,opt and y_opt are very slight. For comparison, three additional working fluids, namely R601a, R245fa, and 0.32R245fa/0.68R601a, are also taken into account.