Claprood, Maxime

Duchesne, Mathieu J.

Giroux, Bernard

Gloaguen, Erwan

Konstantinovskaya, Elena

Malo, Michel

Espace INRS

Combining borehole logs and 2D post-stack seismic profiles by kriging with external drift to build a 3D geological model for CO  2storage in the Bécancour area, Québec, Canada1 1 1Maxime Claprood , Elena Konstantinovskaya , and Michel Malo1Institut National de la Recherche Scientifique, Centre Eau Terre et Environnement (INRS-ETE), 490, rue de la Couronne, Québec, Qc, Canada;2Geological Survey of Canada (Québec), 490, rue de la Couronne, Québec, Qc, Canada;1 1 1,2Erwan Gloaguen , Bernard Giroux , Mathieu Duchesne , SUMMARYDeep saline aquifers are identified in the Bécancour area (Québec, Canada), where their potential for the geological storage of CO  is studied . The first step for times, to build a 3D model of the geological units in Bécancour. The results obtained on the geological horizons show a good compromise between the small 2wavelength variations reproduced by the small scale variograms and the long wavelength reproduced by the external drift. The deterministic 3D geological model evaluating and forecasting the injection of CO  in geological units is to build a 3D model. In contrast to Enhanced Oil Recovery (EOR) and other projects where CO  2 2will be used as a starting point to geostatistically model the distribution of lithofacies and petro-physical properties within each formation unit.is injected in operating oil or gas fields or after their closure, the availability of new and coherent data is very limited due to financial considerations. In this study, we krige the tops of the geological formations recorded in the depth domain at 11 wells, using an external trend interpreted from 2D seismic horizons picked in two-way1 - INTRODUCTIONDeep saline aquifers are identified in the sedimentary successions of the St. Lawrence Lowlands in the province of Québec, Canada (Figure 1a for location). Their potential for CO  storage is 2currently studied in the Bécancour area, half-way between Montréal and Québec on the south shore of the St. Lawrence River.As opposed to EOR and others projects where CO  is injected in operating or after-closure oil or gas 2fields, the availability of new and coherent data above deep saline aquifers is limited due to financial considerations. The challenge of this project is to optimally use all kind of existing data to better estimate the reservoir storage capacity and geological variability. This poster presents the first step of the project: building a 3D model of the geological units of the St. Lawrence Lowlands in the Bécancour area.7 - CONCLUSIONThe results obtained by the kriging of formation tops with modeled horizons as an external trend are very encouraging, considering the limited amount of available  seismic and well logs. The horizons modeled in depth exactly match the formation tops identified at the 11 wells located within the limits of the model. The spatial variations of the horizons are evaluated from the seismic data in TWT and transferred to the depth domain by using an external trend during the kriging process. We think that this approach is a very promising workflow in non-EOR CO  storage projects where the amount of funds and thus 2available existing data sets will be a limiting factor.The next step will consist in integrating all kriged horizons representing the sedimentary successions to geostatistically model the distribution of lithofacies and petro-physical properties (porosity, permeability) within each formation unit that will later be used for flow simulations in the deep saline reservoirs in the Bécancour area.Delhomme, J.P., Boucher, M., Meunier, G. and Jensen, F., 1981. Apport de la géostatistique à la description de gaz en aquifère. Revue de l ’Institut Français du Pétrole, 36, 309-327 (in French).Doveton, J.H., 1986. Log analysis of subsurface geology - Concepts and computer methods: John Wiley & Sons, New York, USA, 273pp.Dubrule, O., 2003. Geostatistics for seismic data integration in earth models: Distinguished Instructor Series no.6, sponsored by Society of Exploration Geophysicists, European Association of Geoscientists & Engineers, Tulsa, USA.Konstantinovskaya, E.A., Rodriguez, D., Kirkwood, D., Harris, L.B., and Thériault, R., 2009. Effects of basement structure, sedimentation and erosion on thrust wedge geometry: An example from the Quebec Appalachians and analogue models: Bulletin of Canadian Petroleum Geology, 57, 34-62.Mallet, J.-L., 2002. Geomodeling: Oxford University Press, New York, USA, 599pp.REFERENCESFigure 1. (a) Location of Bécancour, Québec, Canada. Red square is the Bécancour area. (b) Map of the Grenville basement in the Bécancour area in two-way time (TWT) in ms. Green dots are wells location. Blue lines are 2D seismic lines. White star is location of well A198.AcknowledgementsThis research was realized under financial support of Ministère du Développement Durable, de l'Environnement et des Parcs du Québec, JUNEX Inc. is greatly acknowledged for giving access to seismic and well log data sets. Seismic Micro-Technology kindly provided the seismic interpretation software used in this study.3 - CREATION OF TD CHARTSFigure 3. TD Charts at well A198 (Figure 1b for location) created with the help of the gamma-ray log (GR, last column in yellow). Dashed-dotted lines are formations tops identified on logs, with same color code used in Figure 2. Vertical axis is TWT (ms), except for GR log which is depth (m). TD charts are used to tie formation tops identified in depth with seismic horizons interpreted in TWT.Apply running average on sonic logs and compute the seismic velocityApply running average on density logMultiply velocity and density logs to obtain the acoustic impedanceUse of the coefficient of correlation to find best tie between synthetic and observed seismic tracesCompute the reflection coefficients from the acoustic impedanceConvolve the reflection coefficients with a seismic wavelet to obtain the synthetic seismic tracesTie the synthetic traces with the appropriate reflections on observed seismic tracesUse formation tops identified with other borehole logs (GR) to tie synthetic traces to seismic signatures of some geological units on observed seismic2 possibilities- many solutions possible when DT log recorded on limited depth interval.- imprecision in final estimation of TD charts- better identification of reflections on seismic lines.- improved precision in estimation of TD charts.Bécancour(a)2 - DATA SETS??????????30 post-stack 2D seismic lines acquired between 1970 and 2008 (25 of them are used to build the 3D model, for a total of 100km of 2D seismic):use to delineate the lateral and vertical extents of 9 horizons corresponding to major changes in geology.borehole logs (sonic, gamma-ray, density, porosity, electric) at 18 wells (11 with sonic log to build TD charts):use to locate the formation tops with precision.Some horizons have characteristic signatures on seismic lines, corresponding to geological units identified on well logs. The sequence of reflexions picked on Figure 2 represents the following geological units:Utica Shale (dark green);limestone of the Trenton Group (dark blue, low coherency and low amplitude reflections);Black River - Chazy (BRChazy, light blue) and Beekmantown (Bmt, massive dolomite in purple, dolomitic sandstone in pink) Groups, and sandstone of the Cairnside Formation (yellow) as 3 strong amplitude reflections.Sonic logs are used to generate TD charts, to tie seismic data acquired in TWT to well logs recorded in depth. An example of a joint sonic log - 2D seismic analysis completed at 5 wells located close to seismic line L1 (Figure 1) is presented in Figure 2.Synthetic traces (blue wiggles) are generated by convolving the reflection coefficients obtained from the product of the sonic (red curves) and density (not shown) logs with a seismic wavelet. The wavelet was modeled by extracting the 50 nearest shot points from each well from time-depths between 0.1 s to 1.2 s.4 - MODELING HORIZONS IN TWT??2D surfaces corresponding to the geological horizons are interpolated between seismic picks using the discrete smooth interpolator (DSI) (Mallet, 2002) implemented in the software SKUA. To avoid repetitions, only the modeled horizon of the Trenton Group is presented in this poster.Figure 4a presents the picked geological horizons and structural elements identified in TWT on the seismic profiles.Figure 4b presents the modeled horizon of the Trenton Group in TWT.Figure 6. Experimental (blue squares) and modeled (blue lines) areal variograms at azimuths 95 , 125 , 140 , 155 , 5 , 35 , 50 , and 65° computed from the seismic picks of the Trenton Group. X-axis is the range in meters. Modeled variograms is a gaussian model with geometrical anisotropy, with maximum range a  = 1482m, minimum range a  = 1007m, azimuth of 1 2maximum axis = 48°, and sill = 882.° ° ° ° ° ° °Figure 2. Seismic line identified L1 in Figure 1b, with sonic (red) and gamma-ray (yellow and green) logs at wells located at proximity to seismic line. 9 seismic horizons corresponding to geological units are identified. Vertical axis is two-way time (TWT).(b)A198L1lateral extents of model5 - KRIGING FORMATION TOPS WITH EXTERNAL DRIFT?Geological horizons are modeled in depth by kriging the formation tops identified on well logs (Figure 5a), using the horizons modeled in TWT (such as the example of the Trenton Group in Figure 4b) as an external drift to guide the conversion to the depth domain between the wells (Delhomme et al., 1981; Dubrule, 2003).Figure 5. (a) Formation tops identified at 11 wells in depth. (b) Kriged geological horizon of the Trenton Group in depth. White cube defining the model is 13.6km x 10.5km x 3200m.Figure 4. (a) Geological horizons and structural elements picked on 2D seismic lines in TWT. (b) Modeled horizon of the Trenton Group. Color bar is TWT in ms from short time (red) to late time (white). White cube defining the model is 13.6km x 10.5km x 1750ms. Vertical extension not to scale.  6 - BUILDING 3D GEOLOGICAL MODEL IN DEPTH?The methodology presented in Sections 4 & 5 is applied to all 9 horizons to construct a deterministic 3D geological model of the sedimentary succession of the St. Lawrence Lowlands from the top of the Grenville basement (red in Figure 7a) to the top of the Utica Shale (cap rock, dark green in Figure 7a).?Formation tops are identified by Konstantinoskaya et al. (2009) from well logs analyses, mostly using gamma-ray (GR), sonic (DT), density porosity (DPHI), neutron porosity ( NPHI), and photo-electric factor (PEF) from the methodology described by Doveton (1986).Figure 7. Two different representations of the 3D geological model of the St. Lawrence Lowlands with a color code representing (a) each modeled horizon and (b) the depth below sea level. Thin red lines are limits of 2 faults included in the model. White cube defining the model is 13.6km x 10.5km x 3200m.Potsdam (Cairnside) sandstonePotsdam (Covey Hill) sandstoneGrenville basementBmt (Theresa) dolomitic sandstoneBmt (Beauharnois) dolomiteBlack River - Chazy GroupsTrenton limestoneLow Utica - Up Trenton shale & limestoneUtica ShaleA198(a) (b)gray spheres are faults pickedcolor code for horizons is same as Figure 2color code is TWT to surface:450ms (red) to 1100ms (white)vertical axis is TWT vertical axis is TWT??Formation tops identified from the interpretation of selected logs at 11 wells are the only information available in depth in the Bécancour model.2D variograms (Figure 6) defining the geometric structure of the primary variable (formation tops) are computed on the interpolated 2D seismic horizons due to limited amount of hard data available for the formation tops from well logs.Model variogram is fitted to the experimental variogram over a limited distance (< 2000m). Small scale variograms will reproduce the small wavelength variations, and the external drift will support the long wavelength of the extrapolation.(a) (b)color code is depth to surface: 500m (red) to 2500m (white)blue dots are formation tops of the Trenton Groupcolor code for formation tops is same as Figure 2vertical axis is depth (m)vertical exaggeration x2vertical axis is depth (m)vertical exaggeration x2(a)(b)vertical axis is depth (m)vertical exaggeration x2vertical axis is depth (m)vertical exaggeration x2color code is depth to surface: 500m (red) to 2500m (white)color code for horizons is same as Figure 2Natural ResourcesCanadaRessources naturellesCanadaVelocity ties are mostly based on qualitative fit between 2D post-stack seismic lines and sonic logs using synthetic seismograms. While there is a lack of quantitative validation of the obtained TD charts, the continuity of seismic signatures suggests that the velocity analyses are reliable. 

Combining borehole logs and 2D post-stack seismic profiles by kriging with external drift to build a 3D geological model for CO2 storage in the Bécancour area, Québec, Canada.

Combining borehole logs and 2D post-stack seismic profiles by kriging with external drift to build a 3D geological model for CO2 storage in the Bécancour area, Québec, Canada.

Abstract

Similar works

Full text

Available Versions

Espace INRS