首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 46 毫秒
1.
Shear-wave velocity is a key parameter for calibrating monitoring time-lapse 4D seismic data during CO2-EOR (Enhanced Oil Recovery) and CO2 sequestration. However, actual S-wave velocity data are lacking, especially in 4D data for CO2 sequestration because wells are closed after the CO2 injection and seismic monitoring is continued but no well log data are acquired. When CO2 is injected into a reservoir, the pressure and saturation of the reservoirs change as well as the elastic parameters of the reservoir rocks. We propose a method to predict the S-wave velocity in reservoirs at different pressures and porosities based on the Hertz–Mindlin and Gassmann equations. Because the coordination number is unknown in the Hertz–Mindlin equation, we propose a new method to predict it. Thus, we use data at different CO2 injection stages in the Gao89 well block, Shengli Oilfield. First, the sand and mud beds are separated based on the structural characteristics of the thin sand beds and then the S-wave velocity as a function of reservoir pressure and porosity is calculated. Finally, synthetic seismic seismograms are generated based on the predicted P- and S-wave velocities at different stages of CO2 injection.  相似文献   

2.
Geophysical techniques play key roles in the measuring, monitoring, and verifying the safety of CO2 sequestration and in identifying the efficiency of CO2-enhanced oil recovery. Although geophysical monitoring techniques for CO2 sequestration have grown out of conventional oil and gas geophysical exploration techniques, it takes a long time to conduct geophysical monitoring, and there are many barriers and challenges. In this paper, with the initial objective of performing CO2 sequestration, we studied the geophysical tasks associated with evaluating geological storage sites and monitoring CO2 sequestration. Based on our review of the scope of geophysical monitoring techniques and our experience in domestic and international carbon capture and sequestration projects, we analyzed the inherent difficulties and our experiences in geophysical monitoring techniques, especially, with respect to 4D seismic acquisition, processing, and interpretation.  相似文献   

3.
The spatio-temporal slip distribution of the earthquake that occurred on 8 August 2017 in Jiuzhaigou, China, was estimated from the teleseismic body wave and near-field Global Navigation Satellite System (GNSS) data (coseismic displacements and high-rate GPS data) based on a finite fault model. Compared with the inversion results from the teleseismic body waves, the near-field GNSS data can better restrain the rupture area, the maximum slip, the source time function, and the surface rupture. The results show that the maximum slip of the earthquake approaches 1.4 m, the scalar seismic moment is ~ 8.0 × 1018 N·m (Mw?≈?6.5), and the centroid depth is ~ 15 km. The slip is mainly driven by the left-lateral strike-slip and it is initially inferred that the seismogenic fault occurs in the south branch of the Tazang fault or an undetectable fault, a NW-trending left-lateral strike-slip fault, and belongs to one of the tail structures at the easternmost end of the eastern Kunlun fault zone. The earthquake rupture is mainly concentrated at depths of 5–15 km, which results in the complete rupture of the seismic gap left by the previous four earthquakes with magnitudes >?6.0 in 1973 and 1976. Therefore, the possibility of a strong aftershock on the Huya fault is low. The source duration is ~ 30 s and there are two major ruptures. The main rupture occurs in the first 10 s, 4 s after the earthquake; the second rupture peak arrives in ~ 17 s. In addition, the Coulomb stress study shows that the epicenter of the earthquake is located in the area where the static Coulomb stress change increased because of the 12 May 2017 Mw7.9 Wenchuan, China, earthquake. Therefore, the Wenchuan earthquake promoted the occurrence of the 8 August 2017 Jiuzhaigou earthquake.  相似文献   

4.
The phase velocity of seismic waves varies with the propagation frequency, and thus frequency-dependent phenomena appear when CO2 gas is injected into a reservoir. By dynamically considering these phenomena with reservoir conditions it is thus feasible to extract the frequency-dependent velocity factor with the aim of monitoring changes in the reservoir both before and after CO2 injection. In the paper, we derive a quantitative expression for the frequency-dependent factor based on the Robinson seismic convolution model. In addition, an inversion equation with a frequency-dependent velocity factor is constructed, and a procedure is implemented using the following four processing steps: decomposition of the spectrum by generalized S transform, wavelet extraction of cross-well seismic traces, spectrum equalization processing, and an extraction method for frequency-dependent velocity factor based on the damped least-square algorithm. An attenuation layered model is then established based on changes in the Q value of the viscoelastic medium, and spectra of migration profiles from forward modeling are obtained and analyzed. Frequency-dependent factors are extracted and compared, and the effectiveness of the method is then verified using a synthetic data. The frequency-dependent velocity factor is finally applied to target processing and oil displacement monitoring based on real seismic data obtained before and after CO2 injection in the G89 well block within Shengli oilfield. Profiles and slices of the frequency-dependent factor determine its ability to indicate differences in CO2 flooding, and the predicting results are highly consistent with those of practical investigations within the well block.  相似文献   

5.
At present, using Eddy Covariance (EC) method to estimate the “true value” of carbon sequestration in terrestrial ecosystem arrests more attention. However, one issue is how to solve the uncertainty of observations (especially the nighttime CO2 flux data) appearing in post-processing CO2 flux data. The ratio of effective and reliable nighttime EC CO2 flux data to all nighttime data is relatively low (commonly, less than 50%) for all the long-term and continuous observation stations in the world. Thus, the processing method of nighttime CO2 flux data and its effect analysis on estimating CO2 flux annual sums are very important. In this paper, the authors analyze and discuss the reasons for underestimating nighttime CO2 flux using EC method, and introduce the general theory and method for processing nighttime CO2 flux data. By analyzing the relationship between nighttime CO2 flux and air fraction velocity u*, we present an alternate method, Average Values Test (AVT), to determine the thresholds of fraction velocity (u*c) for screening the effective nighttime CO2 flux data. Meanwhile, taking the data observed in Yucheng and Changbai Mountains stations for an example, we analyze and discuss the effects of different methods or parameters on nighttime CO2 flux estimations. Finally, based on the data of part ChinaFLUX stations and related literatures, empirical models of nighttime respiration at different sites in ChinaFLUX are summarized.  相似文献   

6.
The sequestration of CO2 in subsurface reservoirs constitutes an immediate counter‐measure to reduce anthropogenic emissions of CO2, now recognized by international scientific panels to be the single most critical factor driving the observed global climatic warming. To ensure and verify the safe geological containment of CO2 underground, monitoring of the CO2 site is critical. In the high Arctic, environmental considerations are paramount and human impact through, for instance, active seismic surveys, has to be minimized. Efficient seismic modelling is a powerful tool to test the detectability and imaging capability prior to acquisition and thus improve the characterization of CO2 storage sites, taking both geological setting and seismic acquisition set‐up into account. The unique method presented here avoids the costly generation of large synthetic data sets by employing point spread functions to directly generate pre‐stack depth‐migrated seismic images. We test both a local‐target approach using an analytical filter assuming an average velocity and a full‐field approach accounting for the spatial variability of point spread functions. We assume a hypothetical CO2 plume emplaced in a sloping aquifer inspired by the conditions found at the University of Svalbard CO2 lab close to Longyearbyen, Svalbard, Norway, constituting an unconventional reservoir–cap rock system. Using the local‐target approach, we find that even the low‐to‐moderate values of porosity (5%–18%) measured in the reservoir should be sufficient to induce significant change in seismic response when CO2 is injected. The sensitivity of the seismic response to changes in CO2 saturation, however, is limited once a relatively low saturation threshold of 5% is exceeded. Depending on the illumination angle provided by the seismic survey, the quality of the images of five hypothetical CO2 plumes of varying volume differs depending on the steepness of their flanks. When comparing the resolution of two orthogonal 2D surveys to a 3D survey, we discover that the images of the 2D surveys contain significant artefacts, the CO2‐brine contact is misplaced and an additional reflector is introduced due to the projection of the point spread function of the unresolvable plane onto the imaging plane. All of these could easily lead to a misinterpretation of the behaviour of the injected CO2. Our workflow allows for testing the influence of geological heterogeneities in the target aquifer (igneous intrusions, faults, pervasive fracture networks) by utilizing increasingly complex and more realistic geological models as input as more information on the subsurface becomes available.  相似文献   

7.
Time‐lapse seismics is the methodology of choice for remotely monitoring changes in oil/gas reservoir depletion, reservoir stimulation or CO2 sequestration, due to good sensitivity and resolving power at depths up to several kilometres. This method is now routinely applied offshore, however, the use of time‐lapse methodology onshore is relatively rare. The main reason for this is the relatively high cost of commercial seismic acquisition on land. A widespread belief of a relatively poor repeatability of land seismic data prevents rapid growth in the number of land time‐lapse surveys. Considering that CO2 sequestration on land is becoming a necessity, there is a great need to evaluate the feasibility of time‐lapse seismics for monitoring. Therefore, an understanding of the factors influencing repeatability of land seismics and evaluating limitations of the method is crucially important for its application in many CO2 sequestration projects. We analyse several repeated 2D and 3D surveys acquired within the Otway CO2 sequestration pilot project (operated by the Cooperative Research Centre for Greenhouse Technologies, CO2CRC) in Australia, in order to determine the principal limitations of land time‐lapse seismic repeatability and investigate the influence of the main factors affecting it. Our findings are that the intrinsic signal‐to‐noise ratio (S/N, signal to coherent and background noise levels) and the normalized‐root‐mean‐square (NRMS) difference are controlled by the source strength and source type. However, the post‐stack S/N ratio and corresponding NRMS residuals are controlled mainly by the data fold. For very high‐fold data, the source strength and source type are less critical.  相似文献   

8.
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.  相似文献   

9.
Lake Albano (Alban Hills volcanic complex, Central Italy) is located in a densely populated area near Rome. The deep lake waters have significant dissolved CO2 concentrations, probably related to sub-lacustrine fluid discharges fed by a pressurized CO2-rich reservoir. The analytical results of geochemical surveys carried out in 1989–2010 highlight the episodes of CO2 removal from the lake. The total mass of dissolved CO2 decreased from ∼5.8 × 107 kg in 1989 to ∼0.5 × 107 kg in 2010, following an exponential decreasing trend. Calculated values of both dissolved inorganic carbon and CO2 concentrations along the vertical profile of the lake indicate that this decrease is caused by CO2 release from the epilimnion, at depth <9 m, combined with (1) water circulation at depth <95 m and (2) CO2 diffusion from the deeper lake layers. According to this model, Lake Albano was affected by a large CO2 input that coincided with the last important seismic swarm at Alban Hills in 1989, suggesting an intimate relationship between the addition of deep-originated CO2 to the lake and seismic activity. In the case of a CO2 degassing event of an order of magnitude larger than the one that occurred in 1989, the deepest part of Lake Albano would become CO2-saturated, resulting in conditions compatible with the occurrence of a gas outburst. These results reinforce the idea that a sudden CO2 input into the lake may cause the release of a dense gas cloud, presently representing the major volcanic threat for this densely populated area.  相似文献   

10.
Geological sequestration of CO2 (carbon dioxide) shows great potential to reduce Greenhouse gas emissions. However, CO2 injection into geological formations may give rise to a variety of coupled chemical and physical processes. The thermo-hydro-mechanical (THM) impact of CO2 injection can induce fault instability, even possibly lead to seismic activities in and around the disposal reservoir. A sequential coupling approach under some assumptions was proposed in the numerical study to investigate the THM behavior of the CO2 sequestration system concerning the temperature, initial geological stress, injection pressure and CO2 buoyancy. The fault was treated as a flexible contact model. The effects of CO2 injection on the mechanical behavior of the faults were investigated. The Drucker-Prager model and the cap model were used to model the constitutive relationship of formations. The numerical results show that injection pressure sensitively affects the relative slip change of the fault. At the initial stage of the sequestration process, the injection pressure plays a key role in affecting the pore pressure of the formations. However, as time continues, the influence of CO2-induced buoyancy becomes obvious on the pore pressure of the formations. In general, The THM effects of CO2 geosequestration do not affect the mechanical stability of formations and faults.  相似文献   

11.
Long-term measurement of carbon metabolism of old-growth forests is critical to predict their behaviors and to reduce the uncertainties of carbon accounting under changing climate. Eddy covariance technology was applied to investigate the long-term carbon exchange over a 200 year-old Chinese broad-leaved Korean pine mixed forest in the Changbai Mountains (128°28′E and 42°24′N, Jilin Province, P. R. China) since August 2002. On the data obtained with open-path eddy covariance system and CO2 profile measurement system from Jan. 2003 to Dec. 2004, this paper reports (i) annual and seasonal variation of F NEE, F GPP and R E; (ii) regulation of environmental factors on phase and amplitude of ecosystem CO2 uptake and release Corrections due to storage and friction velocity were applied to the eddy carbon flux.LAI and soil temperature determined the seasonal and annual dynamics of FGPP and RE separately. VPD and air temperature regulated ecosystem photosynthesis at finer scales in growing seasons. Water condition at the root zone exerted a significant influence on ecosystem maintenance carbon metabolism of this forest in winter.The forest was a net sink of atmospheric CO2 and sequestered ?449 g C·m?2 during the study period; ?278 and ?171 gC·m?2 for 2003 and 2004 respectively. F GPP and F RE over 2003 and 2004 were ?1332, ?1294 g C·m?2. and 1054, 1124 g C·m?2 respectively. This study shows that old-growth forest can be a strong net carbon sink of atmospheric CO2.There was significant seasonal and annual variation in carbon metabolism. In winter, there was weak photosynthesis while the ecosystem emitted CO2. Carbon exchanges were active in spring and fall but contributed little to carbon sequestration on an annual scale. The summer is the most significant season as far as ecosystem carbon balance is concerned. The 90 days of summer contributed 66.9, 68.9% of F GPP, and 60.4, 62.1% of R E of the entire year.  相似文献   

12.
We present a new petro-elastical and numerical-simulation methodology to compute synthetic seismograms for reservoirs subject to CO2 sequestration. The petro-elastical equations model the seismic properties of reservoir rocks saturated with CO2, methane, oil and brine. The gas properties are obtained from the van der Waals equation and we take into account the absorption of gas by oil and brine, as a function of the in situ pore pressure and temperature. The dry-rock bulk and shear moduli can be obtained either by calibration from real data or by using rock-physics models based on the Hertz-Mindlin and Hashin-Shtrikman theories. Mesoscopic attenuation due to fluids effects is quantified by using White's model of patchy saturation, and the wet-rock velocities are calculated with Gassmann equations by using an effective fluid modulus to describe the velocities predicted by White's model. The simulations are performed with a poro-viscoelastic modeling code based on Biot's theory, where viscoelasticity is described by generalizing the solid/fluid coupling modulus to a relaxation function. Using the pseudo-spectral method, which allows general material variability, a complete and accurate characterization of the reservoir can be obtained. A simulation, that considers the Utsira sand of the North Sea, illustrates the methodology.  相似文献   

13.
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.  相似文献   

14.
We present an uncertainty analysis of ecological process parameters and CO2 flux components (R eco, NEE and gross ecosystem exchange (GEE)) derived from 3 years’ continuous eddy covariance measurements of CO2 fluxes at subtropical evergreen coniferous plantation, Qianyanzhou of ChinaFlux. Daily-differencing approach was used to analyze the random error of CO2 fluxes measurements and bootstrapping method was used to quantify the uncertainties of three CO2 flux components. In addition, we evaluated different models and optimization methods in influencing estimation of key parameters and CO2 flux components. The results show that: (1) Random flux error more closely follows a double-exponential (Laplace), rather than a normal (Gaussian) distribution. (2) Different optimization methods result in different estimates of model parameters. Uncertainties of parameters estimated by the maximum likelihood estimation (MLE) are lower than those derived from ordinary least square method (OLS). (3) The differences between simulated Reco, NEE and GEE derived from MLE and those derived from OLS are 12.18% (176 g C·m−2·a−1), 34.33% (79 g C·m−2·a−1) and 5.4% (92 g C·m−2·a−1). However, for a given parameter optimization method, a temperature-dependent model (T_model) and the models derived from a temperature and water-dependent model (TW_model) are 1.31% (17.8 g C·m−2·a−1), 2.1% (5.7 g C·m−2·a−1), and 0.26% (4.3 g C·m−2·a−1), respectively, which suggested that the optimization methods are more important than the ecological models in influencing uncertainty in estimated carbon fluxes. (4) The relative uncertainty of CO2 flux derived from OLS is higher than that from MLE, and the uncertainty is related to timescale, that is, the larger the timescale, the smaller the uncertainty. The relative uncertainties of Reco, NEE and GEE are 4%−8%, 7%−22% and 2%−4% respectively at annual timescale. Supported by the National Natural Science Foundation of China (Grant No. 30570347), Innovative Research International Partnership Project of the Chinese Academy of Sciences (Grant No. CXTD-Z2005-1) and National Basic Research Program of China (Grant No. 2002CB412502)  相似文献   

15.
Distribution of parameters characterizing soil response during the 1999 Chi-Chi, Taiwan, earthquake (M w = 7.6) around the fault plane is studied. The results of stochastic finite-fault simulations performed in Pavlenko and Wen (2008) and constructed models of soil behavior at 31 soil sites were used for the estimation of amplification of seismic waves in soil layers, average stresses, strains, and shear moduli reduction in the upper 30 m of soil, as well as nonlinear components of soil response during the Chi-Chi earthquake. Amplification factors were found to increase with increasing distance from the fault (or, with decreasing the level of “input” motion to soil layers), whereas average stresses and strains, shear moduli reduction, and nonlinear components of soil response decrease with distance as ~ r ?1 . The area of strong nonlinearity, where soil behavior is substantially nonlinear (the content of nonlinear components in soil response is more than ~40–50% of the intensity of the response), and spectra of oscillations on the surface take the smoothed form close to E(f) ~ f ?n , is located within ~20–25 km from the fault plane (~ 1/4 of its length). Nonlinearity decreases with increasing distance from the fault, and at ~40–50 km from the fault (~ 1/2 of the fault length), soil response becomes virtually linear. Comparing soil behavior in near-fault zones during the 1999 Chi-Chi, the 1995 Kobe (M w = 6.8), and the 2000 Tottori (Japan) (M w = 6.7) earthquakes, we found similarity in the behavior of similar soils and predominance of the hard type of soil behavior. Resonant phenomena in upper soil layers were observed at many studied sites; however, during the Chi-Chi earthquake they involved deeper layers (down to ~ 40–60 m) than during lesser-magnitude Kobe and Tottori earthquakes.  相似文献   

16.
Between the years 2008 and 2013, approximately 67 kilotons of CO2 have been injected at the Ketzin site, Germany. As part of the geophysical monitoring programme, time‐lapse electrical resistivity tomography has been applied using crosshole and surface‐downhole measurements of electrical resistivity tomography. The data collection of electrical resistivity tomography is partly based on electrodes that are permanently installed in three wells at the site (one injection well and two observation wells). Both types of ERT measurements consistently show the build‐up of a CO2‐related resistivity signature near the injection point. Based on the imaged resistivity changes and a petrophysical model, CO2 saturation levels are estimated. These CO2 saturations are interpreted in conjunction with CO2 saturations inferred from neutron‐gamma loggings. Apart from the CO2–brine substitution response in the observed resistivity changes, significant imprints from the dynamic behaviour of the CO2 in the reservoir are observed.  相似文献   

17.
The effects of elevated CO2 on vegetation dynamics and the hydrological cycle have been widely studied at the site level. However, quantitative assessments of these effects on a regional scale remain a challenge. We conducted numerical simulations to predict the possible responses of vegetation and the hydrological cycle in the Sino-Mongolia arid and semi-arid region (SMASR) to doubled CO2 and its associated climate change using the Community Earth System Model in tandem with a dynamic global vegetation model. The results showed that the doubled CO2 had a positive effect on the leaf area index of the SMASR, but its associated climate change exerted a negative effect in most parts of the SMASR. Although climate change had a weak negative effect on ground runoff at the regional scale, a 4.74 mm increase was predicted under the combined effect of doubled CO2 and climate change, largely due to the positive effect of doubled CO2. Spatially, the evident increase in ground runoff, which primarily occurred in the southeastern part of the SMASR, resulted from decreased ground evaporation and canopy transpiration under the doubled CO2 condition. A negative effect was predicted in the central west as a result of increased temperature and a changed precipitation under doubled CO2. These findings implied that the condition of water resources would be improved slightly under a doubled CO2 condition, whereas there would be a larger spatial heterogeneity in relation to different sensitivities of vegetation and hydrological variables to doubled CO2 and associated climate change.  相似文献   

18.
We estimate the corner frequencies of 20 crustal seismic events from mainshock–aftershock sequences in different tectonic environments (mainshocks 5.7 < M W < 7.6) using the well-established seismic coda ratio technique (Mayeda et al. in Geophys Res Lett 34:L11303, 2007; Mayeda and Malagnini in Geophys Res Lett, 2010), which provides optimal stability and does not require path or site corrections. For each sequence, we assumed the Brune source model and estimated all the events’ corner frequencies and associated apparent stresses following the MDAC spectral formulation of Walter and Taylor (A revised magnitude and distance amplitude correction (MDAC2) procedure for regional seismic discriminants, 2001), which allows for the possibility of non-self-similar source scaling. Within each sequence, we observe a systematic deviation from the self-similar \( M_{0} \propto \mathop f\nolimits_{\text{c}}^{ - 3} \) line, all data being rather compatible with \( M_{0} \propto \mathop f\nolimits_{\text{c}}^{ - (3 + \varepsilon )} \) , where ε > 0 (Kanamori and Rivera in Bull Seismol Soc Am 94:314–319, 2004). The deviation from a strict self-similar behavior within each earthquake sequence of our collection is indicated by a systematic increase in the estimated average static stress drop and apparent stress with increasing seismic moment (moment magnitude). Our favored physical interpretation for the increased apparent stress with earthquake size is a progressive frictional weakening for increasing seismic slip, in agreement with recent results obtained in laboratory experiments performed on state-of-the-art apparatuses at slip rates of the order of 1 m/s or larger. At smaller magnitudes (M W < 5.5), the overall data set is characterized by a variability in apparent stress of almost three orders of magnitude, mostly from the scatter observed in strike-slip sequences. Larger events (M W > 5.5) show much less variability: about one order of magnitude. It appears that the apparent stress (and static stress drop) does not grow indefinitely at larger magnitudes: for example, in the case of the Chi–Chi sequence (the best sampled sequence between M W 5 and 6.5), some roughly constant stress parameters characterize earthquakes larger than M W ~ 5.5. A representative fault slip for M W 5.5 is a few tens of centimeters (e.g., Ide and Takeo in J Geophys Res 102:27379–27391, 1997), which corresponds to the slip amount at which effective lubrication is observed, according to recent laboratory friction experiments performed at seismic slip velocities (V ~ 1 m/s) and normal stresses representative of crustal depths (Di Toro et al. in Nature in press, 2011, and references therein). If the observed deviation from self-similar scaling is explained in terms of an asymptotic increase in apparent stress (Malagnini et al. in Pure Appl Geophys, 2014, this volume), which is directly related to dynamic stress drop on the fault, one interpretation is that for a seismic slip of a few tens of centimeters (M W ~ 5.5) or larger, a fully lubricated frictional state may be asymptotically approached.  相似文献   

19.
In the context of geological carbon sequestration (GCS), carbon dioxide (CO2) is often injected into deep formations saturated with a brine that may contain dissolved light hydrocarbons, such as methane (CH4). In this multicomponent multiphase displacement process, CO2 competes with CH4 in terms of dissolution, and CH4 tends to exsolve from the aqueous into a gaseous phase. Because CH4 has a lower viscosity than injected CO2, CH4 is swept up into a ‘bank’ of CH4‐rich gas ahead of the CO2 displacement front. On the one hand, this may provide a useful tracer signal of an approaching CO2 front. On the other hand, the emergence of gaseous CH4 is undesirable because it poses a leakage risk of a far more potent greenhouse gas than CO2 if the cap rock is compromised. Open fractures or faults and wells could result in CH4 contamination of overlying groundwater aquifers as well as surface emissions. We investigate this process through detailed numerical simulations for a large‐scale GCS pilot project (near Cranfield, Mississippi) for which a rich set of field data is available. An accurate cubic‐plus‐association equation‐of‐state is used to describe the non‐linear phase behavior of multiphase brine‐CH4‐CO2 mixtures, and breakthrough curves in two observation wells are used to constrain transport processes. Both field data and simulations indeed show the development of an extensive plume of CH4‐rich (up to 90 mol%) gas as a consequence of CO2 injection, with important implications for the risk assessment of future GCS projects.  相似文献   

20.
Carbon capture and storage is a viable greenhouse gas mitigation technology and the Sleipner CO2 sequestration site in the North Sea is an excellent example. Storage of CO2 at the Sleipner site requires monitoring over large areas, which can successfully be accomplished with time lapse seismic imaging. One of the main goals of CO2 storage monitoring is to be able to estimate the volume of the stored CO2 in the reservoir. This requires a parametrization of the subsurface as exact as possible. Here we use elastic 2D time‐domain full waveform inversion in a time lapse manner to obtain a P‐wave velocity constrain directly in the depth domain for a base line survey in 1994 and two post‐injection surveys in 1999 and 2006. By relating velocity change to free CO2 saturation, using a rock physics model, we find that at the considered location the aquifer may have been fully saturated in some places in 1999 and 2006.  相似文献   

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

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