Geostatistical seismic inversion for frontier exploration

2017 ◽  
Vol 5 (4) ◽  
pp. T477-T485 ◽  
Author(s):  
Ângela Pereira ◽  
Rúben Nunes ◽  
Leonardo Azevedo ◽  
Luís Guerreiro ◽  
Amílcar Soares
Keyword(s):  
Acoustic Impedance ◽  
Seismic Reflection ◽  
Seismic Inversion ◽  
Well Log ◽  
Data Set ◽  
Seismic Reflection Data ◽  
Log Data ◽  
The Real ◽  
Reflection Data ◽  
Impedance Models

Numerical 3D high-resolution models of subsurface petroelastic properties are key tools for exploration and production stages. Stochastic seismic inversion techniques are often used to infer the spatial distribution of the properties of interest by integrating simultaneously seismic reflection and well-log data also allowing accessing the spatial uncertainty of the retrieved models. In frontier exploration areas, the available data set is often composed exclusively of seismic reflection data due to the lack of drilled wells and are therefore of high uncertainty. In these cases, subsurface models are usually retrieved by deterministic seismic inversion methodologies based exclusively on the existing seismic reflection data and an a priori elastic model. The resulting models are smooth representations of the real complex geology and do not allow assessing the uncertainty. To overcome these limitations, we have developed a geostatistical framework that allows inverting seismic reflection data without the need of experimental data (i.e., well-log data) within the inversion area. This iterative geostatistical seismic inversion methodology simultaneously integrates the available seismic reflection data and information from geologic analogs (nearby wells and/or analog fields) allowing retrieving acoustic impedance models. The model parameter space is perturbed by a stochastic sequential simulation methodology that handles the nonstationary probability distribution function. Convergence from iteration to iteration is ensured by a genetic algorithm driven by the trace-by-trace mismatch between real and synthetic seismic reflection data. The method was successfully applied to a frontier basin offshore southwest Europe, where no well has been drilled yet. Geologic information about the expected impedance distribution was retrieved from nearby wells and integrated within the inversion procedure. The resulting acoustic impedance models are geologically consistent with the available information and data, and the match between the inverted and the real seismic data ranges from 85% to 90% in some regions.

scholarly journals Geophysical imaging of porosity variations in the Danish North Sea chalk

2008 ◽  
Vol 15 ◽  
pp. 17-20 ◽  
Author(s):  
Tanni Abramovitz
Keyword(s):  
North Sea ◽  
Seismic Reflection ◽  
Seismic Inversion ◽  
Rock Properties ◽  
Well Log ◽  
Reservoir Properties ◽  
Seismic Reflection Data ◽  
Log Data ◽  
The North ◽  
The North Sea

More than 80% of the present-day oil and gas production in the Danish part of the North Sea is extracted from fields with chalk reservoirs of late Cretaceous (Maastrichtian) and early Paleocene (Danian) ages (Fig. 1). Seismic reflection and in- version data play a fundamental role in mapping and characterisation of intra-chalk structures and reservoir properties of the Chalk Group in the North Sea. The aim of seismic inversion is to transform seismic reflection data into quantitative rock properties such as acoustic impedance (AI) that provides information on reservoir properties enabling identification of porosity anomalies that may constitute potential reservoir compartments. Petrophysical analyses of well log data have shown a relationship between AI and porosity. Hence, AI variations can be transformed into porosity variations and used to support detailed interpretations of porous chalk units of possible reservoir quality. This paper presents an example of how the chalk team at the Geological Survey of Denmark and Greenland (GEUS) integrates geological, geophysical and petrophysical information, such as core data, well log data, seismic 3-D reflection and AI data, when assessing the hydrocarbon prospectivity of chalk fields.

Prediction of Petrophysical Properties from Seismic Inversion and Neural Network: A case study

2021 ◽  
Author(s):  
Siddharth Garia ◽  
Arnab Kumar Pal ◽  
Karangat Ravi ◽  
Archana M Nair
Keyword(s):  
Neural Network ◽  
Bulk Density ◽  
Low Frequency ◽  
Seismic Inversion ◽  
Well Log ◽  
Petrophysical Properties ◽  
Seismic Section ◽  
Frequency Model ◽  
Seismic Reflection Data ◽  
Reflection Data

<p>Seismic inversion method is widely used to characterize reservoirs and detect zones of interest, i.e., hydrocarbon-bearing zone in the subsurface by transforming seismic reflection data into quantitative subsurface rock properties. The primary aim of seismic inversion is to transform the 3D seismic section/cube into an acoustic impedance (AI) cube. The integration of this elastic attribute, i.e., AI cube with well log data, can thereafter help to establish correlations between AI and different petrophysical properties. The seismic inversion algorithm interpolates and spatially populates data/parameters of wells to the entire seismic section/cube based on the well log information. The case study presented here uses machine learning-neural network based algorithm to extract the different petrophysical properties such as porosity and bulk density from the seismic data of the Upper Assam basin, India. We analyzed three different stratigraphic  units that are established to be producing zones in this basin.</p><p> AI model is generated from the seismic reflection data with the help of colored inversion operator. Subsequently, low-frequency model is generated from the impedance data extracted from the well log information. To compensate for the band limited nature of the seismic data, this low-frequency model is added to the existing acoustic model. Thereafter, a feed-forward neural network (NN) is trained with AI as input and porosity/bulk density as target, validated with NN generated porosity/bulk density with actual porosity/bulk density from well log data. The trained network is thus tested over the entire region of interest to populate these petrophysical properties.</p><p>Three seismic zones were identified from the seismic section ranging from 681 to 1333 ms, 1528 to 1575 ms and 1771 to 1814 ms. The range of AI, porosity and bulk density were observed to be 1738 to 6000 (g/cc) * (m/s), 26 to 38% and 1.95 to 2.46 g/cc respectively. Studies conducted by researchers in the same basin yielded porosity results in the range of 10-36%. The changes in acoustic impedance, porosity and bulk density may be attributed to the changes in lithology. NN method was prioritized over other traditional statistical methods due to its ability to model any arbitrary dependency (non-linear relationships between input and target values) and also overfitting can be avoided. Hence, the workflow presented here provides an estimation of reservoir properties and is considered useful in predicting petrophysical properties for reservoir characterization, thus helping to estimate reservoir productivity.</p>

Shallow seismic reflection study of a glaciated valley

Geophysics ◽  
1998 ◽  
Vol 63 (4) ◽  
pp. 1395-1407 ◽  
Author(s):  
Frank Büker ◽  
Alan G. Green ◽  
Heinrich Horstmeyer
Keyword(s):  
Seismic Reflection ◽  
High Amplitude ◽  
Data Sets ◽  
Seismic Section ◽  
Data Set ◽  
Seismic Reflection Data ◽  
Northern Switzerland ◽  
Reflection Data ◽  
Shallow Seismic ◽  
Glaciated Valley

Shallow seismic reflection data were recorded along two long (>1.6 km) intersecting profiles in the glaciated Suhre Valley of northern Switzerland. Appropriate choice of source and receiver parameters resulted in a high‐fold (36–48) data set with common midpoints every 1.25 m. As for many shallow seismic reflection data sets, upper portions of the shot gathers were contaminated with high‐amplitude, source‐generated noise (e.g., direct, refracted, guided, surface, and airwaves). Spectral balancing was effective in significantly increasing the strength of the reflected signals relative to the source‐generated noise, and application of carefully selected top mutes ensured guided phases were not misprocessed and misinterpreted as reflections. Resultant processed sections were characterized by distributions of distinct seismic reflection patterns or facies that were bounded by quasi‐continuous reflection zones. The uppermost reflection zone at 20 to 50 ms (∼15 to ∼40 m depth) originated from a boundary between glaciolacustrine clays/silts and underlying glacial sands/gravels (till) deposits. Of particular importance was the discovery that the deepest part of the valley floor appeared on the seismic section at traveltimes >180 ms (∼200 m), approximately twice as deep as expected. Constrained by information from boreholes adjacent to the profiles, the various seismic units were interpreted in terms of unconsolidated glacial, glaciofluvial, and glaciolacustrine sediments deposited during two principal phases of glaciation (Riss at >100 000 and Würm at ∼18 000 years before present).

Migration with reduced artifacts from internal multiple reflections

Geophysics ◽  
2020 ◽  
Vol 85 (4) ◽  
pp. A25-A29
Author(s):  
Lele Zhang
Keyword(s):  
Seismic Reflection ◽  
Computational Cost ◽  
Data Sets ◽  
Square Root ◽  
Multiple Reflections ◽  
Data Set ◽  
Seismic Reflection Data ◽  
Reflection Data ◽  
Recent Developments ◽  
Multiple Elimination

Migration of seismic reflection data leads to artifacts due to the presence of internal multiple reflections. Recent developments have shown that these artifacts can be avoided using Marchenko redatuming or Marchenko multiple elimination. These are powerful concepts, but their implementation comes at a considerable computational cost. We have derived a scheme to image the subsurface of the medium with significantly reduced computational cost and artifacts. This scheme is based on the projected Marchenko equations. The measured reflection response is required as input, and a data set with primary reflections and nonphysical primary reflections is created. Original and retrieved data sets are migrated, and the migration images are multiplied with each other, after which the square root is taken to give the artifact-reduced image. We showed the underlying theory and introduced the effectiveness of this scheme with a 2D numerical example.

Shallow 3D seismic-reflection imaging of fracture zones in crystalline rock

Geophysics ◽  
2007 ◽  
Vol 72 (6) ◽  
pp. B149-B160 ◽  
Author(s):  
Cedric Schmelzbach ◽  
Heinrich Horstmeyer ◽  
Christopher Juhlin
Keyword(s):  
Seismic Reflection ◽  
Crystalline Rock ◽  
Three Dimensions ◽  
Disposal Site ◽  
3D Seismic ◽  
Fracture Zones ◽  
Data Set ◽  
Seismic Reflection Data ◽  
Reflection Data ◽  
Receiver Array

A limited 3D seismic-reflection data set was used to map fracture zones in crystalline rock for a nuclear waste disposal site study. Seismic-reflection data simultaneously recorded along two roughly perpendicular profiles (1850 and [Formula: see text] long) and with a [Formula: see text] receiver array centered at the intersection of the lines sampled a [Formula: see text] area in three dimensions. High levels of source-generated noise required a processing sequence involving surface-consistent deconvolution, which effectively increased the strength of reflected signals, and a linear [Formula: see text] filtering scheme to suppress any remaining direct [Formula: see text]-wave energy. A flexible-binning scheme significantly balanced and increased the CMP fold, but the offset and azimuth distributions remain irregular; a wide azimuth range and offsets [Formula: see text] are concentrated in the center of the survey area although long offsets [Formula: see text] are only found at the edges of the site. Three-dimensional dip moveout and 3D poststack migration were necessary to image events with conflicting dips up to about 40°. Despite the irregular acquisition geometry and the high level of source-generated noise, we obtained images rich in structural detail. Seven continuous to semicontinuous reflection events were traced through the final data volume to a maximum depth of around [Formula: see text]. Previous 2D seismic-reflection studies and borehole data indicate that fracture zones are the most likely cause of the reflections.

scholarly journals Multidimensional scaling for the evaluation of a geostatistical seismic elastic inversion methodology

Geophysics ◽  
2014 ◽  
Vol 79 (1) ◽  
pp. M1-M10 ◽  
Author(s):  
Leonardo Azevedo ◽  
Ruben Nunes ◽  
Pedro Correia ◽  
Amílcar Soares ◽  
Luis Guerreiro ◽  
...  
Keyword(s):  
Multidimensional Scaling ◽  
Metric Space ◽  
Synthetic Data ◽  
Real Data ◽  
Seismic Inversion ◽  
Inversion Algorithm ◽  
Data Set ◽  
The Real ◽  
Geostatistical Inversion ◽  
Impedance Models

Due to the nature of seismic inversion problems, there are multiple possible solutions that can equally fit the observed seismic data while diverging from the real subsurface model. Consequently, it is important to assess how inverse-impedance models are converging toward the real subsurface model. For this purpose, we evaluated a new methodology to combine the multidimensional scaling (MDS) technique with an iterative geostatistical elastic seismic inversion algorithm. The geostatistical inversion algorithm inverted partial angle stacks directly for acoustic and elastic impedance (AI and EI) models. It was based on a genetic algorithm in which the model perturbation at each iteration was performed recurring to stochastic sequential simulation. To assess the reliability and convergence of the inverted models at each step, the simulated models can be projected in a metric space computed by MDS. This projection allowed distinguishing similar from variable models and assessing the convergence of inverted models toward the real impedance ones. The geostatistical inversion results of a synthetic data set, in which the real AI and EI models are known, were plotted in this metric space along with the known impedance models. We applied the same principle to a real data set using a cross-validation technique. These examples revealed that the MDS is a valuable tool to evaluate the convergence of the inverse methodology and the impedance model variability among each iteration of the inversion process. Particularly for the geostatistical inversion algorithm we evaluated, it retrieves reliable impedance models while still producing a set of simulated models with considerable variability.

Minimizing field operations in shallow 3‐D seismic reflection surveying

Geophysics ◽  
2001 ◽  
Vol 66 (6) ◽  
pp. 1761-1773 ◽  
Author(s):  
Roman Spitzer ◽  
Alan G. Green ◽  
Frank O. Nitsche
Keyword(s):  
Seismic Reflection ◽  
Cost Effective ◽  
Unconsolidated Sediments ◽  
Data Set ◽  
Seismic Reflection Data ◽  
Northern Switzerland ◽  
Shallow Subsurface ◽  
Reflection Data ◽  
Effective Manner ◽  
Original Survey

By appropriately decimating a comprehensive shallow 3‐D seismic reflection data set recorded across unconsolidated sediments in northern Switzerland, we have investigated the potential and limitations of four different source‐receiver acquisition patterns. For the original survey, more than 12 000 shots and 18 000 receivers deployed on a [Formula: see text] grid resulted in common midpoint (CMP) data with an average fold of ∼40 across a [Formula: see text] area. A principal goal of our investigation was to determine an acquisition strategy capable of producing reliable subsurface images in a more efficient and cost‐effective manner. Field efforts for the four tested acquisition strategies were approximately 50%, 50%, 25%, and 20% of the original effort. All four data subsets were subjected to a common processing sequence. Static corrections, top‐mute functions, and stacking velocities were estimated individually for each subset. Because shallow reflections were difficult to discern on shot and CMP gathers generated with the lowest density acquisition pattern (20% field effort) such that dependable top‐mute functions could not be estimated, data resulting from this acquisition pattern were not processed to completion. Of the three fully processed data subsets, two (50% field effort and 25% field effort) yielded 3‐D migrated images comparable to that derived from the entire data set, whereas the third (50% field effort) resulted in good‐quality images only in the shallow subsurface because of a lack of far‐offset data. On the basis of these results, we concluded that all geological objectives associated with our particular study site, which included mapping complex lithological units and their intervening shallow dipping boundaries, would have been achieved by conducting a 3‐D seismic reflection survey that was 75% less expensive than the original one.

scholarly journals Structural Interpretation of Yamama and Naokelekan Formations in Tuba Oil Field Using 2D Seismic Data

2021 ◽  
Vol 54 (2B) ◽  
pp. 55-64
Author(s):  
Belal M. Odeh
Keyword(s):  
Seismic Data ◽  
Seismic Reflection ◽  
Synthetic Seismogram ◽  
Oil Field ◽  
Well Log ◽  
Seismic Reflection Data ◽  
Depth Maps ◽  
Reflection Data ◽  
Interpretation Process ◽  
Southern Iraq

This research includes structure interpretation of the Yamama Formation (Lower Cretaceous) and the Naokelekan Formation (Jurassic) using 2D seismic reflection data of the Tuba oil field region, Basrah, southern Iraq. The two reflectors (Yamama and Naokelekan) were defined and picked as peak and tough depending on the 2D seismic reflection interpretation process, based on the synthetic seismogram and well log data. In order to obtain structural settings, these horizons were followed over all the regions. Two-way travel-time maps, depth maps, and velocity maps have been produced for top Yamama and top Naokelekan formations. The study concluded that certain longitudinal enclosures reflect anticlines in the east and west of the study area representing Zubair and Rumaila fold confined between them a fold consist of two domes represents Tuba fold with the same trending of Zubair and Rumaila structures. The study confirmed the importance of this field as a reservoir of the accumulation of hydrocarbons.

Sign in / Sign up

Export Citation Format

Share Document