首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
Geological heterogeneities affect the dynamics of carbon dioxide (CO2) plumes in subsurface environments in important ways. Previously, we showed how the dynamics of CO2 plumes are influenced by the multiscaled sedimentary architecture in deep brine fluvial-type reservoirs. The results confirm that representing small-scale features and the corresponding heterogeneity in saturation functions, along with hysteresis in saturation functions, are all critical to understanding capillary trapping processes. Here, we show that when heterogeneity and hysteresis are represented, the two conventional approaches for defining saturation functions, Brooks-Corey and van Genuchten, represent fundamentally different physical systems. The Brooks-Corey approach represents heterogeneity in entry pressures, and leads to trapping by capillary pinning. The van Genuchten approach represents a network of pores transporting the nonwetting fluid, across rock types, with negligible capillary entry pressure, and leads to capillary retardation. These differences significantly affect the large-scale characteristics of CO2 plumes (i.e., their mass, shape, and position).  相似文献   

2.
The typical shape of a capillary-pressure curve is either convex (e.g., Brooks–Corey model) or S-shaped (e.g., van Genuchten model). It is not universally agreed which model reflects natural rocks better. The difference between the two models lies in the representation of the capillary entry pressure. This difference does not lead to significantly different simulation results for modeling CO2 sequestration in aquifers without considering CO2 dissolution. However, we observe that the van-Genuchten-type capillary-pressure model accelerates CO2 solubility trapping significantly compared with the Brooks–Corey-type model. We also show that the simulation results are very sensitive to the slope of the van-Genuchten-type curve around the entry-pressure region. For the representative examples we study, the differences can be so large as to have complete dissolution of the CO2 plume versus persistence of over 50% of the plume over a 5000-year period.The cause of such sensitivity to the capillary-pressure model is studied. Particularly, we focus on how the entry pressure is represented in each model. We examine the mass-transfer processes under gravity-capillary equilibrium, molecular diffusion, convective mixing, and in the presence of small-scale heterogeneities. Laboratory measurement of capillary-pressure curves and some important implementation issues of capillary-pressure models in numerical simulators are also discussed. Most CO2 sequestration simulations in the literature employ one of the two capillary-pressure models. It is important to recognize that these two representations lead to very different predictions of long-term CO2 sequestration.  相似文献   

3.
Integrated modeling of basin- and plume-scale processes induced by full-scale deployment of CO2 storage was applied to the Mt. Simon Aquifer in the Illinois Basin. A three-dimensional mesh was generated with local refinement around 20 injection sites, with approximately 30 km spacing. A total annual injection rate of 100 Mt CO2 over 50 years was used. The CO2-brine flow at the plume scale and the single-phase flow at the basin scale were simulated. Simulation results show the overall shape of a CO2 plume consisting of a typical gravity-override subplume in the bottom injection zone of high injectivity and a pyramid-shaped subplume in the overlying multilayered Mt. Simon, indicating the important role of a secondary seal with relatively low-permeability and high-entry capillary pressure. The secondary-seal effect is manifested by retarded upward CO2 migration as a result of multiple secondary seals, coupled with lateral preferential CO2 viscous fingering through high-permeability layers. The plume width varies from 9.0 to 13.5 km at 200 years, indicating the slow CO2 migration and no plume interference between storage sites. On the basin scale, pressure perturbations propagate quickly away from injection centers, interfere after less than 1 year, and eventually reach basin margins. The simulated pressure buildup of 35 bar in the injection area is not expected to affect caprock geomechanical integrity. Moderate pressure buildup is observed in Mt. Simon in northern Illinois. However, its impact on groundwater resources is less than the hydraulic drawdown induced by long-term extensive pumping from overlying freshwater aquifers.  相似文献   

4.
  相似文献   

5.
Previous experimental studies of capillary barriers have identified highly hysteretic soil moisture retention characteristics in the materials used. In this study, numerical modelling is used to analyse the role of soil moisture hysteresis in capillary barrier functioning. Comparisons between published experimental results and model simulations indicate that soil moisture hysteresis was a necessary inclusion in the modelling approach to adequately reproduce pore water pressure distributions and the timing of breakthrough occurrences. Under hypothetical intermittent infiltration and evaporation conditions, the predicted volumetric water content in the moisture retention layer was significantly different for hysteretic and non‐hysteretic models. The hysteresis effect was found to be dependent on the nature of infiltration–evaporation cycling, although the predicted volume of flow through the hysteretic barrier was lower than that of the non‐hysteretic case, regardless of the nature of the cyclic upper boundary conditions. For practical engineering designs, where the water leakage through the barrier is the primary concern, the inclusion of soil moisture hysteresis in numerical modelling is needed. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

6.
Geological storage of carbon dioxide (CO2) is a promising technology for reducing atmospheric emissions. The large discrepancy in the time- and length-scales between up-dip migration of buoyant supercritical CO2 and the sinking fingers of dissolved CO2 poses a challenge for numerical simulations aimed at describing the fate of the plume. Hence, several investigators have suggested methods to simplify the problem, but to date there has been no reference solution with which these simplified models can be compared. We investigate the full problem of Darcy-based two-phase flow with gravity-current propagation and miscible convective mixing, using high-resolution numerical simulations. We build on recent developments of the Automatic Differentiation - General Purpose Research Simulator (AD-GPRS) at Stanford. The results show a CO2 plume that travels for 5000 years reaching a final distance of 14 km up-dip from the injection site. It takes another 2000 years before the CO2 is completely trapped as residual (40%) and dissolved (60%) CO2. Dissolution causes a significant reduction of the plume speed. While fingers of dissolved CO2 appear under the propagating gravity current, the resident brine does not become fully saturated with CO2 anywhere under the plume. The overall mass transfer of CO2 into the brine under the plume remains practically constant for several thousands of years. These results can be used as a benchmark for verification, or improvements, of simplified (reduced-dimensionality, upscaled) models. Our results indicate that simplified models need to account for: (i) reduced dissolution due to interaction with the plume, and (ii) gradual reduction of the local dissolution rate after the fingers begin to interact with the bottom of the aquifer.  相似文献   

7.
During geologic storage of carbon dioxide (CO2), trapping of the buoyant CO2 after injection is essential in order to minimize the risk of leakage into shallower formations through a fracture or abandoned well. Models for the subsurface behavior of the CO2 are useful for the design, implementation, and long-term monitoring of injection sites, but traditional reservoir-simulation tools are currently unable to resolve the impact of small-scale trapping processes on fluid flow at the scale of a geologic basin. Here, we study the impact of solubility trapping from convective dissolution on the up-dip migration of a buoyant gravity current in a sloping aquifer. To do so, we conduct high-resolution numerical simulations of the gravity current that forms from a pair of miscible analogue fluids. Our simulations fully resolve the dense, sinking fingers that drive the convective dissolution process. We analyze the dynamics of the dissolution flux along the moving CO2–brine interface, including its decay as dissolved buoyant fluid accumulates beneath the buoyant current. We show that the dynamics of the dissolution flux and the macroscopic features of the migrating current can be captured with an upscaled sharp-interface model.  相似文献   

8.
This study presents the impact of fractures on CO2 transport, capillary pressure and storage capacity by conducting both experimental and numerical studies. A series of laboratory experiment tests was designed with both a homogeneous and a fractured core under CO2 storage conditions. The experimental results reveal a piston-like brine displacement with gravity override effects in the homogeneous core regardless of CO2 injection rates. In the fractured core, however, two distinctive types of brine displacements were observed; one showing brine displacement only in the fracture whereas the other shows brine displacement both in the fracture and matrix with different rates, which were dependent on the magnitude of the pressure build-up in the matrix. The injectivity in the fractured core was twice of the homogeneous core, while the amount of calculated CO2 in the homogeneous core was over 1.5 times greater than the fractured core. Salt precipitation, which is likely to occur near injection wells, was observed in the experiments; X-ray images enabled the observation of salt-precipitation during CO2-flooding tests. Finally, numerical simulations predict free-phase CO2 transfer between fracture and matrix in a fracture-matrix system. Pressure gradients between the fracture and matrix enforced CO2 to transfer from the fracture into matrix at the front of the CO2 plume, whereas, the reversal of pressure gradients at the rear zone of the CO2 plume reversed the transfer process. The variation of CO2 saturation within the fracture was caused by fracture aperture variations, and local variations of fracture permeability control the free-phase CO2 transfer between the fracture and matrix.  相似文献   

9.
In this study we performed three categories of steady- and unsteady-state core-flooding experiments to investigate capillary trapping, relative permeability, and capillary pressure, in a scCO2 + SO2/brine/limestone system at elevated temperature and pressure conditions, i.e., 60 °C and 19.16 MPa. We used a Madison limestone core sample acquired from the Rock Springs Uplift in southwest Wyoming. We carried out two sets of steady-state drainage-imbibition relative permeability experiments with different initial brine saturations to study hysteresis. We found that the final scCO2 + SO2 drainage relative permeability was very low, i.e., 0.04. We also observed a rapid reduction in the scCO2-rich phase imbibition relative permeability curve, which resulted in a high residual trapping. The results showed that between 62.8% and more than 76% of the initial scCO2 + SO2 at the end of drainage was trapped by capillary trapping mechanism (trapping efficiency). We found that at higher initial brine saturations, the trapping efficiency was higher. The maximum initial and residual scCO2-rich phase saturations at the end of primary drainage and imbibition were 0.525 and 0.329, respectively. Each drainage-imbibition cycle was followed by a dissolution process to re-establish Sw = 1. The dissolution brine relative permeabilities for both cycles were also obtained. We characterized the scCO2 + SO2/brine capillary pressure hysteresis behavior through unsteady-state primary drainage, imbibition, and secondary drainage experiments. We observed negative imbibition capillary pressure curve indicative of possible wettability alteration throughout the experiments due to contact with scCO2 + SO2/brine fluid system. The trapping results were compared to those reported in literature for other carbonate core samples. We noticed slightly more residual trapping in our sample, which might be attributed to heterogeneity, different viscosity ratio, and pore-space topologies. The impact of dynamic effects, i.e., high brine flow rate imbibition tests, on trapping of the scCO2-rich phase was also explored. We performed two imbibition experiments with relatively high brine flow rates. The residual scCO2 saturation dropped to 0.291 and 0.262 at the end of the first and second imbibition tests, i.e., 11.5% and 20.4%, respectively, compared to 0.329 under capillary-dominated regime.  相似文献   

10.
This study aims to evaluate the application of 222Rn in groundwater as a tracer for monitoring CO2 plume migration in a shallow groundwater system, which is important to detect potential CO2 leakage in the carbon capture and storage (CCS) project. For this research, an artificial CO2-infused water injection experiment was performed in a shallow aquifer by monitoring hydrogeochemical parameters, including 222Rn. Radon in groundwater can be a useful tracer because of its sensitivity to sudden changes in subsurface environment. To monitor the CO2 plume migration, the data were analysed based on (a) the influence of mixing processes on the distribution of 222Rn induced by the artificial injection experiment and (b) the influence of a carrier gas role by CO2 on the variation of 222Rn. The spatio-temporal distributions of radon concentrations were successfully explained in association with horizontal and vertical mixing processes by the CO2-infused water injection. Additionally, the mixing ratios of each monitoring well were calculated, quantitatively confirming the influence of these mixing processes on the distribution of radon concentrations. Moreover, one monitoring well showed a high positive relationship between 222Rn and Total dissolved inorganic carbon (TIC) by the carrier gas effect of CO2 through volatilization from the CO2 plume. It indicated the applicability of 222Rn as a sensitive tracer to directly monitor CO2 leakage. When with a little effect of carrier gas, natural 222Rn in groundwater can be used to compute mixing ratio of CO2-infused water indicative of CO2 migration pathways. CO2 carrier gas effect can possibly increase 222Rn concentration in groundwater and, if fully verified with more field tests, will pose a great potential to be used as a natural tracer for CO2.  相似文献   

11.
The injection of CO2 at the Ketzin pilot CO2 storage site started in June 2008 and ended in August 2013. During the 62 months of injection, a total amount of about 67 kt of CO2 was injected into a saline aquifer. A third repeat three‐dimensional seismic survey, serving as the first post‐injection survey, was acquired in 2015, aiming to investigate the recent movement of the injected CO2. Consistent with the previous two time‐lapse surveys, a predominantly west–northwest migration of the gaseous CO2 plume in the up‐dip direction within the reservoir is inferred in this first post‐injection survey. No systematic anomalies are detected through the reservoir overburden. The extent of the CO2 plume west of the injection site is almost identical to that found in the 2012 second repeat survey (after injection of 61 kt); however, there is a significant decrease in its size east of the injection site. Assessment of the CO2 plume distribution suggests that the decrease in the size of the anomaly may be due to multiple factors, such as limited vertical resolution, CO2 dissolution, and CO2 migration into thin layers, in addition to the effects of ambient noise. Four‐dimensional seismic modelling based on dynamic flow simulations indicates that a dynamic balance between the newly injected CO2 after the second repeat survey and the CO2 migrating into thin layers and being dissolved was reached by the time of the first post‐injection survey. In view of the significant uncertainties in CO2 mass estimation, both patchy and non‐patchy saturation models for the Ketzin site were taken into consideration.  相似文献   

12.
Geologic carbon sequestration (GCS) is being considered as a climate change mitigation option in many future energy scenarios. Mathematical modeling is routinely used to predict subsurface CO2 and resident brine migration for the design of injection operations, to demonstrate the permanence of CO2 storage, and to show that other subsurface resources will not be degraded. Many processes impact the migration of CO2 and brine, including multiphase flow dynamics, geochemistry, and geomechanics, along with the spatial distribution of parameters such as porosity and permeability. In this article, we review a set of multiphase modeling approaches with different levels of conceptual complexity that have been used to model GCS. Model complexity ranges from coupled multiprocess models to simplified vertical equilibrium (VE) models and macroscopic invasion percolation models. The goal of this article is to give a framework of conceptual model complexity, and to show the types of modeling approaches that have been used to address specific GCS questions. Application of the modeling approaches is shown using five ongoing or proposed CO2 injection sites. For the selected sites, the majority of GCS models follow a simplified multiphase approach, especially for questions related to injection and local‐scale heterogeneity. Coupled multiprocess models are only applied in one case where geomechanics have a strong impact on the flow. Owing to their computational efficiency, VE models tend to be applied at large scales. A macroscopic invasion percolation approach was used to predict the CO2 migration at one site to examine details of CO2 migration under the caprock.  相似文献   

13.
Growing evidence suggests microbial respiration of dissolved organic carbon (DOC) may be a principal driver of subsurface dissolution and cave formation in eogenetic carbonate rock. Analyses of samples of vadose zone gasses, and geochemical and hydrological data collected from shallow, uncased wells on San Salvador Island, Bahamas, suggest tidally varying water tables may help fuel microbial respiration and dissolution through oxygenation. Respiration of soil organic carbon transported to water tables generates dysaerobic to anaerobic groundwater, limiting aerobic microbial processes. Positive correlations of carbon dioxide (CO2), radon-222 (222Rn) and water table elevation indicate, however, that tidal pumping of water tables pulls atmospheric air that is rich in oxygen, and low in CO2 and 222Rn, into contact with the tidal capillary fringe during falling tides. Ratios of CO2 and O2 in vadose gas relative to the atmosphere indicate this atmospheric oxygen fuels respiration within newly-exposed, wetted bedrock. Deficits of expected CO2 relative to O2 concentrations indicate some respired CO2 is likely removed by carbonate mineral dissolution. Tidal pumping also appears capable of transferring oxygen to the freshwater lens, where it could also contribute to respiration and dissolution; dissolved oxygen concentrations at the water table are at least 5% saturated and decline to anaerobic conditions 1–2 m below. Our results demonstrate how tidal pumping of air to vadose zones can drive mineral dissolution reactions that are focused near water tables and may contribute to the formation of laterally continuous vuggy horizons and potentially caves. © 2020 John Wiley & Sons, Ltd.  相似文献   

14.
Tracking the migration of the CO2 plume is essential in order to better manage the operation of geologic sequestration of CO2. However, the large cost of most monitoring technologies, such as time-lapse seismic, limits their application. We investigated the application of a probabilistic history matching methodology using routine measurements of injection data, which are affected by the presence of large-scale heterogeneities, as an inexpensive alternative to track the migration of CO2 plume in an aquifer. The approach is demonstrated first through a synthetic example in which the ability to detect the presence of flow barriers was tested. In a second example, we applied our method to the In Salah field, one of the largest geological sequestration projects in the world, where the main direction of high permeability features was inferred. The accuracy and reproducibility of the results show that our approach for assisted history matching is an economic and viable option for plume monitoring during geologic CO2 sequestration.  相似文献   

15.
Geophysical methods have been used experimentally during the last decade, a period of strong development, being adopted as complementary techniques for characterizing and monitoring hydrocarbon and gas reservoirs. In this study, we evaluated the ability of the controlled source electromagnetic (CSEM) method to monitor the storage of CO2 at the Research Laboratory on Geological Storage of CO2 at Hontomín (Burgos, Spain). Two aspects of the CSEM monitoring were examined considering the geoelectrical structure at the site, the technological constraints and the noise conditions of the Hontomín area. Borehole-to-surface simulations were performed to evaluate the detectability of the resistivity changes in the reservoir and the capacity to determine the location of the CO2 plume. The synthetic time-lapse study explores the possibilities of CSEM monitoring with a deep electric source. Three depths of the source are analyzed: above the plume, inside the plume, and beneath the stored CO2. In terms of the Hontomín storage site, the study confirmed that a deep electric source located beneath the injection depth can provide valuable information on the behavior of the stored CO2.  相似文献   

16.
This paper presents an overview of earlier and current methods of modelling atmospheric dispersion, and proposes and evaluates a screening model for operation over the World Wide Web. The On Line Atmospheric Dispersion Model (OLADMO) is a quasi boundary layer parameterised Gaussian plume model with additional algorithms to account for plume rise, building wake effects and deposition processes. The Monin–Obukhov length boundary layer parameter is utilised to define six stability classes in order to determine atmospheric turbulence and stability, whilst new equations, derived from an intercomparison study of old and next generation dispersion models, are used to calculate the horizontal and vertical dispersion coefficients σ y and σ z . Using data from two field experiments in Copenhagen, Denmark and Lillestrøm, Norway, the model results from OLADMO are found to compare favourably with the results from several old and next generation dispersion models. As a consequence of the unique nature of the meteorological and location factors of the Lillestrøm experiment, all models struggled to represent the concentrations observed during the field study adequately. However, OLADMO was the best performing model in this case, with a mean normalised crosswind integrated concentration 13% closer to the mean observed concentration than its nearest competitor. Because the evaluation of the model was conducted with a limited dataset, several limitations and improvements to both the model and experimental procedure are suggested.  相似文献   

17.
On 27 February 2007, a new eruption occurred on Stromboli which lasted until 2 April. It was characterized by effusive activity on the Sciara del Fuoco and by a paroxysmal event (15 March). This crisis represented an opportunity for us to refine the model that had been developed previously (2002–2003 eruption) and to improve our understanding of the relationship between the magmatic dynamics of the volcano and the geochemical variations in the fluids. In particular, the evaluation of the dynamic equilibrium between the volatiles (CO2 and SO2) released from the magma and the corresponding fluids discharged from the summit area allowed us to evaluate the level of criticality of the volcanic activity. One of the major accomplishments of this study is a 4-year database of summit soil CO2 flux on the basis of which we define the thresholds (low–medium–high) for this parameter that are empirically based on the natural volcanological evolution of Stromboli. The SO2 fluxes of the degassing plume and the CO2 fluxes emitted from the soil at Pizzo Sopra la Fossa are also presented. It is noteworthy that geochemical signals of volcanic unrest have been clearly identified before, during and after the effusive activity. These signals were found almost simultaneously in the degassing plume (SO2 flux) and in soil degassing (CO2 flux) at the summit, although the two degassing processes are shown to be clearly different. The interpretation of the results will be useful for future volcanic surveillance at Stromboli.  相似文献   

18.
Stream temperature (Ts) is a key water quality parameter that controls several biological, ecological, and chemical processes in aquatic systems. In forested headwaters, exchanges of energy across air-water-streambed interfaces may influence Ts regimes, especially during storm events as the sources of runoff change over space and time. Analysis of the hysteretic behaviour of Ts during storm events may provide insights into rainfall-runoff responses, but such relationships have not been thoroughly investigated. As such, our objectives were to (a) quantify the variability of stream temperature hysteresis across seasons in different sub-regions and (b) investigate the relationship between the hysteretic response and catchment characteristics. Ts hysteresis during storm events was assessed based on the hysteresis index (HI), which describes the directionality of hysteresis loops, and the temperature response index (TRI), which indicates whether Ts increased or decreased during a storm event. We analysed Ts data from 10 forested headwater reaches in two sub-regions (McGarvey and West Fork Tectah) in Northern California. We also performed a clustering analysis to examine the relationship amongst HI, TRI, topographic metrics, and meteorological characteristics of the study areas. Overall, the hysteretic behaviour of Ts varied across seasons—the greatest HI occurred during spring and summer. Interestingly, in the McGarvey streams the variability in Ts hysteresis co-varied strongly with topographic metrics (i.e., upslope accumulative area, average channel slope, topographic wetness index). Comparatively, in West Fork Tectah the variability of Ts hysteresis co-varied most strongly with meteorological metrics (i.e., antecedent rainfall events, solar radiation, and air temperature). Variables such as the gradient between stream and air temperatures, slope, and wetted width were significant for both sub-regional hysteretic patterns. We posit that the drivers of Ts response during storms are likely dependent on catchment physiographic characteristics. Our study also illustrated the potential utility of stream temperature as a tracer for improving the understanding of hydrologic connectivity and shifts in the dominant runoff contributions to streamflow during storm events.  相似文献   

19.
Here, we report the first continuous data of geochemical parameters acquired directly from the active summit crater of Vulcano. This approach provides a means to better investigate deep geochemical processes associated with the degassing system of Vulcano Island. In particular, we report on soil CO2 fluxes from the upper part of Vulcano, a closed-conduit volcano, from September 2007 to October 2010. Large variations in the soil CO2 and plume SO2 fluxes (order of magnitude), coinciding with other discontinuous geochemical parameters (CO2 concentrations in fumarole gas) and physical parameters (increase of shallow seismic activity and fumarole temperatures) have been recorded. The results from this work suggest new prospects for strengthening geochemical monitoring of volcanic activity and for improving the constraints in the construction of a “geochemical model”, this being a necessary condition to better understand the functioning of volcanic systems.  相似文献   

20.
Soil erosion, transport and deposition by water drastically affect the distribution of soil organic carbon (SOC) within a landscape. Moreover, soil redistribution may have a large impact on the exchange of carbon (C) between the pedosphere and the atmosphere. One of the large information gaps within this research domain, concerns the fate of SOC after erosion by water. According to different (mainly laboratory) studies, soil redistribution leads to aggregate breakdown, thereby exposing the contained SOC to mineralization. Our study aims to quantify the extent to which such increased mineralization occurs in a real field situation. Carbon dioxide (CO2)‐efflux was measured in the field after an important erosion event for a continuous period of 112 days. The specific situation on the field ensured that almost none of eroded SOC was exported from the field. Measurements of CO2‐efflux were done in areas with sediment deposition, as well as in comparable areas without sedimentation. Comparison of these measurements allowed the net effect of soil deposition on CO2‐efflux to be assessed. Field data were complemented by measurements on incubated, undisturbed soil core samples, in order to disentangle the contribution of environmental factors (moisture, temperature) from any erosional effect on CO2‐efflux. Results of these measurements on the field showed that CO2‐efflux was regulated by a complex interplay of different factors (mostly soil porosity, soil moisture and soil temperature). In combination with the incubation measurements, it could be concluded that the processes of erosion and transport indeed led to an increased mineralization of SOC, as a result of aggregate breakdown and exposure of previously encapsulated SOC. This effect was, however, much smaller than observed in previous laboratory studies. Moreover, it was only important in the first weeks, immediately after the erosion event. The calculated net erosional effect on CO2‐efflux represented a mere 1·6% of total SOC, originally present in the soil. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号