Predicting subsurface CO2 movement: From laboratory to field scale

Geophysics ◽  
2012 ◽  
Vol 77 (3) ◽  
pp. M27-M37 ◽  
Author(s):  
Ranjana Ghosh ◽  
Mrinal K. Sen

Finding an appropriate model for time-lapse seismic monitoring of [Formula: see text]-sequestered carbonate reservoir poses a great challenge because carbonate-rocks have varying textures and highly reactive rock-fluid system. We introduced a frequency-dependent model based on Eshelby’s inclusion and differential effective medium (DEM) theory that can account for heterogeneity in microstructure of rocks and squirt flow. We showed that the estimated velocities from the modified DEM theory match well with the laboratory measurements (ultrasonic) of velocities of carbonate rocks saturated with [Formula: see text]-rich water. The theory predicts significant decrease in saturated P- and S-wave velocities in the seismic frequency band as a consequence of porosity and permeability enhancement by the process of chemical dissolution of carbonates with the saturating fluid. The study also showed the combined effect of chemical reaction and free [Formula: see text] saturation on the elastic properties of rock. We observed that the velocity dispersion and attenuation increased from complete gas saturation to water saturation. The proposed model can be used to invert geophysical measurements to detect changes in elastic properties of a carbonate reservoir and interpret the extent of [Formula: see text] movement with time. These are the key elements to ensure that sequestration will not damage underground geologic formation and [Formula: see text] storage is secure and environmentally acceptable.

2017 ◽  
Vol 5 (1) ◽  
pp. SB69-SB80 ◽  
Author(s):  
Jingjing Xu ◽  
Maojin Tan ◽  
Xiaochang Wang ◽  
Chunping Wu

Estimation of S-wave velocity is one of the most critical steps for prestack seismic inversion. Based on the petrophysical model of fractured carbonate rocks, theoretical methods are firstly investigated for estimating P- and S-wave velocities in the presence of fractures. Then, the methods of calculating elastic properties in fractured carbonate rocks are discussed. The mineral concentration, total porosity, and fracture porosity from core X-ray diffraction and routine core measurements or log interpretation results are used to estimate the P- and S-wave velocities. In the given carbonate rock model, the elastic properties of carbonate rocks with different porosity and fractures are calculated. Two field tests prove that the proposed new method is effective and accurate. Furthermore, the model is useful for fluid identification, which is one of the most outstanding problems for carbonate reservoir description. The simulation results suggest that the larger the fracture porosity is, the easier fluid typing. In Tahe Oilfield, the elastic properties of different fluid zones indicate that bulk modulus and Young’s modulus are more sensitive to fluid than shear modulus, the Lamé constant, and Poisson’s ratio.


2021 ◽  
Author(s):  
Kangxu Ren ◽  
Junfeng Zhao ◽  
Jian Zhao ◽  
Xilong Sun

Abstract At least three very different oil-water contacts (OWC) encountered in the deepwater, huge anticline, pre-salt carbonate reservoirs of X oilfield, Santos Basin, Brazil. The boundaries identification between different OWC units was very important to help calculating the reserves in place, which was the core factor for the development campaign. Based on analysis of wells pressure interference testing data, and interpretation of tight intervals in boreholes, predicating the pre-salt distribution of igneous rocks, intrusion baked aureoles, the silicification and the high GR carbonate rocks, the viewpoint of boundaries developed between different OWC sub-units in the lower parts of this complex carbonate reservoirs had been better understood. Core samples, logging curves, including conventional logging and other special types such as NMR, UBI and ECS, as well as the multi-parameters inversion seismic data, were adopted to confirm the tight intervals in boreholes and to predicate the possible divided boundaries between wells. In the X oilfield, hundreds of meters pre-salt carbonate reservoir had been confirmed to be laterally connected, i.e., the connected intervals including almost the whole Barra Velha Formation and/or the main parts of the Itapema Formation. However, in the middle and/or the lower sections of pre-salt target layers, the situation changed because there developed many complicated tight bodies, which were formed by intrusive diabase dykes and/or sills and the tight carbonate rocks. Many pre-salt inner-layers diabases in X oilfield had very low porosity and permeability. The tight carbonate rocks mostly developed either during early sedimentary process or by latter intrusion metamorphism and/or silicification. Tight bodies were firstly identified in drilled wells with the help of core samples and logging curves. Then, the continuous boundary were discerned on inversion seismic sections marked by wells. This paper showed the idea of coupling the different OWC units in a deepwater pre-salt carbonate play with complicated tight bodies. With the marking of wells, spatial distributions of tight layers were successfully discerned and predicated on inversion seismic sections.


Geophysics ◽  
2020 ◽  
Vol 85 (4) ◽  
pp. H51-H60
Author(s):  
Feng Zhou ◽  
Iraklis Giannakis ◽  
Antonios Giannopoulos ◽  
Klaus Holliger ◽  
Evert Slob

In oil drilling, mud filtrate penetrates into porous formations and alters the compositions and properties of the pore fluids. This disturbs the logging signals and brings errors to reservoir evaluation. Drilling and logging engineers therefore deem mud invasion as undesired and attempt to eliminate its adverse effects. However, the mud-contaminated formation carries valuable information, notably with regard to its hydraulic properties. Typically, the invasion depth critically depends on the formation porosity and permeability. Therefore, if adequately characterized, mud invasion effects could be used for reservoir evaluation. To pursue this objective, we have applied borehole radar to measure mud invasion depth considering its high radial spatial resolution compared with conventional logging tools, which then allows us to estimate the reservoir permeability based on the acquired invasion depth. We investigate the feasibility of this strategy numerically through coupled electromagnetic and fluid modeling in an oil-bearing layer drilled using freshwater-based mud. Time-lapse logging is simulated to extract the signals reflected from the invasion front, and a dual-offset downhole antenna mode enables time-to-depth conversion to determine the invasion depth. Based on drilling, coring, and logging data, a quantitative interpretation chart is established, mapping the porosity, permeability, and initial water saturation into the invasion depth. The estimated permeability is in a good agreement with the actual formation permeability. Our results therefore suggest that borehole radar has significant potential to estimate permeability through mud invasion effects.


Geophysics ◽  
1993 ◽  
Vol 58 (1) ◽  
pp. 30-38 ◽  
Author(s):  
Gang Yu ◽  
Keeva Vozoff ◽  
David W. Durney

Laboratory measurements are described on Permian coals from Wollongong, New South Wales, Australia related to the dependence of ultrasonic P‐ and S‐wave velocities, attenuation, anisotropy and the dynamic elastic moduli on confining pressure, water saturation, and pore pressure. Five independent stiffness constants are used to represent the elastic anisotropy of the specimens as a function of confining pressure and water saturation. The anisotropy is believed to be controlled mainly by the internal structure of the coals, while the pressure dependence of the constants is controlled mainly by randomly oriented cracks. P‐ and S‐wave dispersions were measured on water‐saturated specimens as confining pressures increased from 2 MPa to 40 MPa. The samples represented cores taken both parallel and perpendicular to bedding planes. Velocities along bedding planes are marginally higher than those across bedding planes. This anisotropy is insensitive to confining pressure. Attenuation was also measured, both normal and parallel to bedding planes, on dry and water‐saturated specimens from 2 MPa to 40 MPa confining pressures. The experimental results show that dynamic elastic properties are potential indicators of the states of stress and saturation in coal seams, and provide necessary information for computer modeling and interpreting seismic surveys carried out to assist mine development.


2019 ◽  
Vol 17 (2) ◽  
Author(s):  
M. Wahdanadi Haidar ◽  
Reza Wardhana ◽  
M. Iskan ◽  
M. Syamsu Rosid

The pore systems in carbonate reservoirs are more complex than the pore systems in clastic rocks. There are three types of pores in carbonate rocks: interparticle pores, stiff pores and cracks. The complexity of the pore types can lead to changes in the P-wave velocity by up to 40%, and carbonate reservoir characterization becomes difficult when the S-wave velocity is estimated using the dominant interparticle pore type only. In addition, the geometry of the pores affects the permeability of the reservoir. Therefore, when modelling the elastic modulus of the rock it is important to take into account the complexity of the pore types in carbonate rocks. The Differential Effective Medium (DEM) is a method for modelling the elastic modulus of the rock that takes into account the heterogeneity in the types of pores in carbonate rocks by adding pore-type inclusions little by little into the host material until the required proportion of the material is reached. In addition, the model is optimized by calculating the bulk modulus of the fluid filler porous rock under reservoir conditions using the Adaptive Batzle-Wang method. Once a fluid model has been constructed under reservoir conditions, the model is entered as input for the P-wave velocity model, which is then used to estimate the velocity of the S-wave and the proportion of primary and secondary pore types in the rock. Changes in the characteristics of the P-wave which are sensitive to the presence of fluid lead to improvements in the accuracy of the P-wave model, so the estimated S-wave velocity and the calculated ratio of primary and secondary pores in the reservoir are more reliable.


Author(s):  
Handoyo Handoyo ◽  
M Rizki Sudarsana ◽  
Restu Almiati

Carbonate rock are important hydrocarbon reservoir rocks with complex texture and petrophysical properties (porosity and permeability). These complexities make the prediction reservoir characteristics (e.g. porosity and permeability) from their seismic properties more difficult. The goal of this paper are to understanding the relationship of physical properties and to see the signature carbonate initial rock and shally-carbonate rock from the reservoir. To understand the relationship between the seismic, petrophysical and geological properties, we used rock physics modeling from ultrasonic P- and S- wave velocity that measured from log data. The measurements obtained from carbonate reservoir field (gas production). X-ray diffraction and scanning electron microscope studies shown the reservoir rock are contain wackestone-packstone content. Effective medium theory to rock physics modeling are using Voigt, Reuss, and Hill.  It is shown the elastic moduly proposionally decrease with increasing porosity. Elastic properties and wave velocity are decreasing proporsionally with increasing porosity and shally cemented on the carbonate rock give higher elastic properties than initial carbonate non-cemented. Rock physics modeling can separated zones which rich of shale and less of shale.


2021 ◽  
pp. 60-74
Author(s):  
G. A. Smolyakov ◽  
N. V. Gilmanova ◽  
A. V. Sivkova

The article deals with the determination of the reservoir properties of Permian-age carbonate rocks. There is a section dissection technique, taking into account the fossil organisms prevailing for a particular geological age. It was noted a high content of silica in the lower Artinskian deposits of Toravey and Varandey fields of the Komi Republic. The presence of silicon is associated with an increase in the population of siliceous sponges during this period of sedimentation; this fact caused the maximum values of porosity and permeability in the corresponding intervals. However, this was also the reason for the high values of the residual water saturation factor and, as a result, low oil flow rates from the lower Artinskian stage. The need for detailed correlation and accounting for the content of fossils in the rocks when dissecting the well section made it necessary to systematize the available actual material on core and well testing. It became obvious that the separation of reservoirs and stages at the qualitative level isn't possible, so quantitative estimates of the parameters that are significantly different for the upper and lower Artinskian deposits and allow you to dissect the well section were proposed.


Geophysics ◽  
2019 ◽  
Vol 84 (1) ◽  
pp. MR45-MR59 ◽  
Author(s):  
Jean-Baptiste Regnet ◽  
Jérôme Fortin ◽  
Aurélien Nicolas ◽  
Matthieu Pellerin ◽  
Yves Guéguen

We have provided new insights into the controlling factors of elastic properties in continental carbonate rocks and introduced an applicable model for acoustic-velocity predictions in such a medium. Petrophysical properties (porosity, permeability, P- and S-wave velocities) from laboratory measurements have been coupled with thin-section observations and characterizations, and X-ray diffraction (XRD) analyses. A major achievement is the establishment of the link between the mineralogical composition and the P- and S-wave velocity dispersion at a given porosity. This reflects the subtle interplay between physicochemical and biological precipitation of continental carbonates, which can also be associated with a strong influence of detrital mineralogical inputs. The result is a mineralogical commixture, coupled to a wide array of pore types inherited from the strong ability of carbonate rocks to undergo diagenetic alteration. The proposed model takes into account the elastic moduli of the minerals, porosity, and pore shape, and it is based on the effective medium theory. We have considered the case in which the medium contained randomly oriented pores with different aspect ratios. Overall, the fit between the predicted trends and the experimental data is fairly good, especially for calcite and quartz matrix mineralogy. The results are even better when considering mineralogy inferred from XRD data, although in some case, and despite the aspect ratio variation in both simulations, the model fails to accurately predict the P-wave velocities. This probably means that another factor is at stake beside mineralogy. This can be explained by the limitation of the effective medium approach, which oversimplifies the reality and fails to account for the variability of some aspect ratio from one inclusion to another.


2021 ◽  
Author(s):  
ahmed wattan ◽  
Mohammed AL‑Jawad

Abstract Shear and compressional wave velocities are useful for drilling operations, the exploration of reservoirs, stimulation processes, and hydraulic fracturing. An ultrasonic device will be used in this investigation to anticipate and analyze the elastic characteristics of carbonate rocks. At the summit of the field, the well WQ1-20 obtained samples of the Mishrif formation from a variety of various depths. The number of samples taken from the well is nine from different units whereas the number of samples taken from the main unit (MB2) was five. The relations between the elastic properties for the carbonate rocks with P-and S-waves were defined. The relations between Vp and Vs with elastic properties were defined by applied Regression analysis. The results showed that a linear relationship between P-and S-wave velocities with the elastic properties of the carbonate rocks. It is found that the relationship between Vp and Young's modulus (E) is R2 equal to 0.979 while the relationship between Vs and Young's modulus (E) is R2 equal to 0.925. The relationship between shear modulus and Vs is good in comparison with Vp where the values of R2 were 0.985 and 0.94 respectively. R2 values for the Bulk modulus and Lame's constant of Vp are 0.925 and 0.6, respectively, while the values for Vs are 0.925 and 0.6 for the latter. The relation between Vp/Vs ratio with Poisson’s ratio showed a good R2 with a value of 0.97. When it comes to predicting the dynamic elastic characteristics of a material, the ultrasonic approach may be regarded as a cost-effective, easy, and non-destructive method.


Sign in / Sign up

Export Citation Format

Share Document