Influence of lithology and fluid flow on the temperature distribution in a sedimentary basin: A case study from the Cold Lake area, Alberta, Canada

A detailed study of the groundwater and terrestrial heat flow was carried out over an area of 23,700 km² west of Cold Lake, Alberta, which is part of the western Canada sedimentary basin. The information for the study was provided from data from 3100 wells drilled in the area. The screening and processing of thousands of stratigraphic picks, drillstem test data, bottom hole temperatures and formation water chemistry data was performed mainly using a specially designed software package. As a result, every stratigraphic unit is characterized by appropriate hydraulic and thermal parameters.

A sequence of aquifers, aquitards and aquicludes was differentiated. The groundwater flow in the Paleozoic aquifers is regional in nature and mainly horizontal. The flow in the Cretaceous aquifers is of intermediate type, mainly downward oriented. In general, the permeability of the Cretaceous and Paleozoic strata has such low values that the fluid velocity is less than 1 cm/yr.

The convective heat transport in the hydrostratigraphic sequence is negligible with respect to the conductive heat transfer, as shown by the Peclet number of the fluid and heat flow in porous media. The flow of the terrestrial heat flux from the Precambrian basement of the sedimentary basin to the atmosphere is controlled by the variability in the thermal properties of the formations in the basin.

The geothermal gradients were computed by hydrostratigraphic unit using a linear regression fit to the temperature data. As expected, they show higher values for the less conductive layers, and lower values for the more conductive ones. The weighted average, or the integral geothermal gradient of the whole sedimentary column, was computed by considering the difference between the temperature measured at the Precambrian basement and the annual average temperature at the surface. The areal distribution of the integral geothermal gradient (with an average of 22.0 mK/m) shows a strong correlation with the lithology.

The areal temperature distribution for each hydrostratigraphic unit was analyzed by mapping the deviation of the measured value from the computed geothermal gradient. The lateral heat flow from warmer to colder areas is one order of magnitude smaller than the vertical heat flow. In the more homogeneous units, the lateral heat flow presents a trend that seems to reflect the geometry and lithology.     

Effects of in-situ conditions on relative permeability characteristics of CO2-brine systems

Carbon dioxide capture and geological storage (CCGS) is an emerging technology that is increasingly being considered for reducing greenhouse gas emissions to the atmosphere. Deep saline aquifers provide a very large capacity for CO₂ storage and, unlike hydrocarbon reservoirs and coal beds, are immediately accessible and are found in all sedimentary basins. Proper understanding of the displacement character of CO₂-brine systems at in-situ conditions is essential in ascertaining CO₂ injectivity, migration and trapping in the pore space as a residual gas or supercritical fluid, and in assessing the suitability and safety of prospective CO₂ storage sites. Because of lack of published data, the authors conducted a program of measuring the relative permeability and other displacement characteristics of CO₂-brine systems for sandstone, carbonate and shale formations in central Alberta in western Canada. The tested formations are representative of the in-situ characteristics of deep saline aquifers in compacted on-shore North American sedimentary basins. The results show that the capillary pressure, interfacial tension, relative permeability and other displacements characteristics of CO₂-brine systems depend on the in-situ conditions of pressure, temperature and water salinity, and on the pore size distribution of the sedimentary rock. This paper presents a synthesis and interpretation of the results.     

On the effective thermal and hydraulic conductivity of binary heterogeneous sediments

Heterogeneity is present in geological sedimentary structures at all scales from pore to basin, and its distribution has an impact on transport processes such as heat and fluid flow. The rock masses at any scale need to be characterized by their effective properties at that scale, based on the individual characteristics of the heterogeneous porous medium. The effective thermal and hydraulic conductivity of sediments characterized by a stochastic distribution of heterogeneity is studied using an inverse approach based on numerical experiments. The simulations, covering a large range of conductivity contrasts, are carried out for actual core-scale cases of shale clasts in a sand matrix, and on a diagrammatic cross-section through a clastic sedimentary group at the basin scale.

The effective conductivity depends primarily on the heterogeneity fraction and on the conductivity contrast between heterogeneities and the embedding matrix, a dependency which can be described by a generalized weighted mean model. This model is better suited to estimate the effective conductivity at any scale than other models like the self-consistent, or any of the arithmetic, geometric or harmonic averages. The effective conductivity has an anisotropic character although the individual components are homogeneous and isotropic. The variation in effective conductivity is significant even for small conductivity contrasts, like in heat flow processes, and exhibits an asymptotic behaviour for large conductivity contrasts characteristic of fluid flow processes. The effective conductivity has a second-order dependence on such heterogeneity characteristics as shape, aspect ratio, orientation, and distribution. Depending on these characteristics, the bounds of effective conductivity values can be narrowed further from the extreme bounds expressed by the arithmetic and harmonic averages.     

Spatial characterization of the location of potentially leaky wells penetrating a deep saline aquifer in a mature sedimentary basin

This work was motivated by considerations of potential leakage pathways for CO₂ injected into deep geological formations for the purpose of carbon sequestration. Because existing wells represent a potentially important leakage pathway, a spatial analysis of wells that penetrate a deep aquifer in the Alberta Basin was performed and various statistical measures to quantify the spatial distribution of these wells were presented. The data indicate spatial clustering of wells, due to oil and gas production activities. The data also indicate that the number of wells that could be impacted by CO₂ injection, as defined by the spread of an injected CO₂ plume, varies from several hundred in high well-density areas to about 20 in low-density areas. These results may be applied to other mature continental sedimentary basins in North America and elsewhere, where detailed information on well location and status may not be available.     

New insights into the origin and migration of brines in deep Devonian aquifers, Alberta, Canada

Analysis of hydraulic heads and chemical compositions of Devonian formation waters in the west central part of the Alberta Basin, Canada, characterizes the origin of formation waters and migration of brines. The Devonian succession in the study area lies 2000–5000 m below the ground surface, and has an approximate total thickness of 1000 m and an average slope of 15 m/km. Four Devonian aquifers are present in the study area, which form two aquifer systems [i.e., a Middle–Upper Devonian aquifer system (MUDAS) consisting of the Elk Point and Woodbend–Beaverhill Lake aquifers, and an Upper Devonian aquifer system (UDAS) consisting of the Winterburn and Wabamun aquifers]. The Ireton is an effective aquitard between these two systems in the eastern parts of the study area. The entire Devonian succession is confined below by efficient aquitards of the underlying Cambrian shales and/or the Precambrian basement, and above by overlying Carboniferous shales of the Exshaw and Lower Banff Formations.The formation water chemistry shows that the Devonian succession contains two distinct brine types: a ‘heavy brine,’ located updip, defined approximately by TDS >200 g/l, and a ‘light brine’ with TDS <200 g/l. Hydraulic head distributions suggest that, presently, the ‘light brine’ attempts to flow updip, thereby pushing the ‘heavy brine’ ahead. The interface between the two brines is lobate and forms large-scale tongues that are due to channeled flow along high-permeability pathways. Geological and hydrogeochemical data suggest that the following processes determined the present composition of the ‘light’ and ‘heavy’ brines: original seawater, evaporation beyond gypsum but below halite saturation, dolomitization, clay dehydration, gypsum dewatering, thermochemical sulfate reduction (TSR), and halite dissolution. The influx of meteoric (from the south) and metamorphic (from the west) waters can be recognized only in the ‘light brine.’ Albitization can be unequivocally identified only in the ‘heavy brine.’ The ‘heavy brine’ may be residual Middle Devonian evaporitic brine from the Williston Basin or the Elk Point Basin, or it may have originated from partial dissolution of thick, laterally extensive Middle Devonian evaporite deposits to the east of the study area. The ‘light brine’ most probably originated from dilution of ‘heavy brine’ in post-Laramide times.     

An assessment of the potential impact of climate change on flood risk in Mumbai

Managing risks from extreme events will be a crucial component of climate change adaptation. In this study, we demonstrate an approach to assess future risks and quantify the benefits of adaptation options at a city-scale, with application to flood risk in Mumbai. In 2005, Mumbai experienced unprecedented flooding, causing direct economic damages estimated at almost two billion USD and 500 fatalities. Our findings suggest that by the 2080s, in a SRES A2 scenario, an ??upper bound?? climate scenario could see the likelihood of a 2005-like event more than double. We estimate that total losses (direct plus indirect) associated with a 1-in-100 year event could triple compared with current situation (to $690?C$1,890 million USD), due to climate change alone. Continued rapid urbanisation could further increase the risk level. The analysis also demonstrates that adaptation could significantly reduce future losses; for example, estimates suggest that by improving the drainage system in Mumbai, losses associated with a 1-in-100 year flood event today could be reduced by as much as 70%.,We show that assessing the indirect costs of extreme events is an important component of an adaptation assessment, both in ensuring the analysis captures the full economic benefits of adaptation and also identifying options that can help to manage indirect risks of disasters. For example, we show that by extending insurance to 100% penetration, the indirect effects of flooding could be almost halved. We conclude that, while this study explores only the upper-bound climate scenario, the risk-assessment core demonstrated in this study could form an important quantitative tool in developing city-scale adaptation strategies. We provide a discussion of sources of uncertainty and risk-based tools could be linked with decision-making approaches to inform adaptation plans that are robust to climate change.