A robust workflow for performing joint impedance inversion and its applications

2020 ◽  
Vol 8 (1) ◽  
pp. T141-T149
Author(s):  
Ritesh Kumar Sharma ◽  
Satinder Chopra ◽  
Larry R. Lines
Keyword(s):  
Time Domain ◽  
Seismic Data ◽  
Sedimentary Basin ◽  
Vertical Component ◽  
Data Sets ◽  
Data Set ◽  
Matching Process ◽  
Multicomponent Seismic ◽  
Impedance Inversion ◽  
Straightforward Approach

Multicomponent seismic data offer several advantages for characterizing reservoirs with the use of the vertical component (PP) and mode-converted (PS) data. Joint impedance inversion inverts both of these data sets simultaneously; hence, it is considered superior to simultaneous impedance inversion. However, the success of joint impedance inversion depends on how accurately the PS data are mapped on the PP time domain. Normally, this is attempted by performing well-to-seismic ties for PP and PS data sets and matching different horizons picked on PP and PS data. Although it seems to be a straightforward approach, there are a few issues associated with it. One of them is the lower resolution of the PS data compared with the PP data that presents difficulties in the correlation of the equivalent reflection events on both the data sets. Even after a few consistent horizons get tracked, the horizon matching process introduces some artifacts on the PS data when mapped into PP time. We have evaluated such challenges using a data set from the Western Canadian Sedimentary Basin and then develop a novel workflow for addressing them. The importance of our workflow was determined by comparing data examples generated with and without its adoption.

Addressing artifacts in PP-PS registration prior to performing joint impedance inversion

2020 ◽  
Vol 39 (1) ◽  
pp. 47-52
Author(s):  
Satinder Chopra ◽  
Ritesh Kumar Sharma
Keyword(s):  
Seismic Data ◽  
The United States ◽  
Synthetic Seismograms ◽  
Data Sets ◽  
Well Log ◽  
Data Set ◽  
Anadarko Basin ◽  
Matching Process ◽  
Impedance Inversion ◽  
Integrated Interpretation

Multicomponent seismic data analysis enhances confidence in interpretation by providing mode-converted PS data for imaging the subsurface. Integrated interpretation of PP and PS data begins with the identification of reflections corresponding to similar geologic events on both data sets. This identification is accomplished by carrying out well-log correlation through the generation of PP and PS synthetic seismograms. There are a few issues associated with the approach. One of the issues is that PS data have lower resolution than PP data. This presents difficulties in the correlation of equivalent reflection events on both data sets. Even if few consistent horizons are tracked, the horizon-matching process introduces artifacts on the PS data mapped in PP time. In this paper, we elaborate on such challenges with a data set from the Anadarko Basin in the United States. We then propose a novel workflow to address the challenges.

Sparse reflectivity inversion for nonstationary seismic data with surface-related multiples: Numerical and field-data experiments

Geophysics ◽  
2017 ◽  
Vol 82 (3) ◽  
pp. R199-R217 ◽  
Author(s):  
Xintao Chai ◽  
Shangxu Wang ◽  
Genyang Tang
Keyword(s):  
Seismic Data ◽  
Resolution Enhancement ◽  
Synthetic Data ◽  
Data Sets ◽  
Data Set ◽  
Anelastic Attenuation ◽  
Seismic Resolution ◽  
Text Filtering ◽  
The Stability ◽  
Reflectivity Inversion

Seismic data are nonstationary due to subsurface anelastic attenuation and dispersion effects. These effects, also referred to as the earth’s [Formula: see text]-filtering effects, can diminish seismic resolution. We previously developed a method of nonstationary sparse reflectivity inversion (NSRI) for resolution enhancement, which avoids the intrinsic instability associated with inverse [Formula: see text] filtering and generates superior [Formula: see text] compensation results. Applying NSRI to data sets that contain multiples (addressing surface-related multiples only) requires a demultiple preprocessing step because NSRI cannot distinguish primaries from multiples and will treat them as interference convolved with incorrect [Formula: see text] values. However, multiples contain information about subsurface properties. To use information carried by multiples, with the feedback model and NSRI theory, we adapt NSRI to the context of nonstationary seismic data with surface-related multiples. Consequently, not only are the benefits of NSRI (e.g., circumventing the intrinsic instability associated with inverse [Formula: see text] filtering) extended, but also multiples are considered. Our method is limited to be a 1D implementation. Theoretical and numerical analyses verify that given a wavelet, the input [Formula: see text] values primarily affect the inverted reflectivities and exert little effect on the estimated multiples; i.e., multiple estimation need not consider [Formula: see text] filtering effects explicitly. However, there are benefits for NSRI considering multiples. The periodicity and amplitude of the multiples imply the position of the reflectivities and amplitude of the wavelet. Multiples assist in overcoming scaling and shifting ambiguities of conventional problems in which multiples are not considered. Experiments using a 1D algorithm on a synthetic data set, the publicly available Pluto 1.5 data set, and a marine data set support the aforementioned findings and reveal the stability, capabilities, and limitations of the proposed method.

Common-focus point-based target-oriented imaging approach for continuous seismic reservoir monitoring

Geophysics ◽  
2018 ◽  
Vol 83 (4) ◽  
pp. M41-M48 ◽  
Author(s):  
Hongwei Liu ◽  
Mustafa Naser Al-Ali
Keyword(s):  
Seismic Data ◽  
Rock Physics ◽  
Time Lapse ◽  
Velocity Estimation ◽  
Estimation Methods ◽  
Data Sets ◽  
Data Set ◽  
Velocity Models ◽  
Reservoir Monitoring ◽  
Focus Point

The ideal approach for continuous reservoir monitoring allows generation of fast and accurate images to cope with the massive data sets acquired for such a task. Conventionally, rigorous depth-oriented velocity-estimation methods are performed to produce sufficiently accurate velocity models. Unlike the traditional way, the target-oriented imaging technology based on the common-focus point (CFP) theory can be an alternative for continuous reservoir monitoring. The solution is based on a robust data-driven iterative operator updating strategy without deriving a detailed velocity model. The same focusing operator is applied on successive 3D seismic data sets for the first time to generate efficient and accurate 4D target-oriented seismic stacked images from time-lapse field seismic data sets acquired in a [Formula: see text] injection project in Saudi Arabia. Using the focusing operator, target-oriented prestack angle domain common-image gathers (ADCIGs) could be derived to perform amplitude-versus-angle analysis. To preserve the amplitude information in the ADCIGs, an amplitude-balancing factor is applied by embedding a synthetic data set using the real acquisition geometry to remove the geometry imprint artifact. Applying the CFP-based target-oriented imaging to time-lapse data sets revealed changes at the reservoir level in the poststack and prestack time-lapse signals, which is consistent with the [Formula: see text] injection history and rock physics.

Recording the Aurora at Seismometers across Alaska

2020 ◽  
Vol 91 (6) ◽  
pp. 3039-3053 ◽  
Author(s):  
Carl Tape ◽  
Adam T. Ringler ◽  
Don L. Hampton
Keyword(s):  
Magnetic Field ◽  
Seismic Data ◽  
Vertical Component ◽  
Seismic Array ◽  
Magnetic Shielding ◽  
Data Sets ◽  
S Waves ◽  
Auroral Activity ◽  
Seismic Sensors ◽  
Magnetometer Data

Abstract We examine three continuously recording data sets related to the aurora: all-sky camera images, three-component magnetometer data, and vertical-component, broadband seismic data as part of the EarthScope project (2014 to present). Across Alaska there are six all-sky cameras, 13 magnetometers, and >200 seismometers. The all-sky images and magnetometers have the same objective, which is to monitor space weather and improve our understanding of auroral activity, including the influence on magnetic fields in the ground. These variations in the magnetic field are also visible on seismometers, to the extent that during an auroral event, the long-period (40–800 s) waves recorded by a seismometer are magnetic field variations, not true ground motion. Although this is a problem—one that can be rectified with magnetic shielding at each seismometer site—it is also an opportunity because the present seismic array in Alaska is much broader than the coverage by magnetometers and all-sky cameras. Here we focus on three aurora events and document a direct link between aurora images in the night sky and seismometer recordings on ground. Simultaneous recordings by magnetometers provide a critical link between the sky images and the seismometer recordings. We document qualitative correlations among sky, magnetic, and seismic data. The findings suggest that the signature of auroral activity is widespread across seismometers in Alaska, implying that the seismic array could be used to enhance the spatial resolution of the existing network of all-sky cameras and magnetometers. Future efforts to improve the multisensor seismic stations in Alaska, for the purpose of monitoring seismic and auroral activity, should consider installation of all-sky cameras, installation of magnetometers, and magnetic shielding of seismic sensors.

Accurate poststack acoustic-impedance inversion by well-log calibration

Geophysics ◽  
2009 ◽  
Vol 74 (5) ◽  
pp. R59-R67 ◽  
Author(s):  
Igor B. Morozov ◽  
Jinfeng Ma
Keyword(s):  
Seismic Data ◽  
Acoustic Impedance ◽  
Joint Inversion ◽  
Low Frequency ◽  
Real Data ◽  
Data Sets ◽  
Frequency Model ◽  
Physical Constraints ◽  
Impedance Inversion ◽  
Reliable Solution

The seismic-impedance inversion problem is underconstrained inherently and does not allow the use of rigorous joint inversion. In the absence of a true inverse, a reliable solution free from subjective parameters can be obtained by defining a set of physical constraints that should be satisfied by the resulting images. A method for constructing synthetic logs is proposed that explicitly and accurately satisfies (1) the convolutional equation, (2) time-depth constraints of the seismic data, (3) a background low-frequency model from logs or seismic/geologic interpretation, and (4) spectral amplitudes and geostatistical information from spatially interpolated well logs. The resulting synthetic log sections or volumes are interpretable in standard ways. Unlike broadly used joint-inversion algorithms, the method contains no subjectively selected user parameters, utilizes the log data more completely, and assesses intermediate results. The procedure is simple and tolerant to noise, and it leads to higher-resolution images. Separating the seismic and subseismic frequency bands also simplifies data processing for acoustic-impedance (AI) inversion. For example, zero-phase deconvolution and true-amplitude processing of seismic data are not required and are included automatically in this method. The approach is applicable to 2D and 3D data sets and to multiple pre- and poststack seismic attributes. It has been tested on inversions for AI and true-amplitude reflectivity using 2D synthetic and real-data examples.

3D ray+Born migration/inversion—Part 2: Application to the SEG/EAGE overthrust experiment

Geophysics ◽  
2003 ◽  
Vol 68 (4) ◽  
pp. 1357-1370 ◽  
Author(s):  
Stéphane Operto ◽  
Gilles Lambaré ◽  
Pascal Podvin ◽  
Philippe Thierry
Keyword(s):  
Spatial Resolution ◽  
Seismic Data ◽  
Velocity Model ◽  
Nonuniform Distribution ◽  
Data Sets ◽  
Data Set ◽  
Quantitative Estimates ◽  
3D Processing ◽  
Type Data ◽  
New Formula

The SEG/EAGE overthrust model is a synthetic onshore velocity model that was used to generate several large synthetic seismic data sets using acoustic finite‐difference modeling. From this database, several realistic subdata sets were extracted and made available for testing 3D processing methods. For example, classic onshore‐type data‐acquisition geometries are available such as a swath acquisition, which is characterized by a nonuniform distribution of long offsets with azimuth and midpoints. In this paper, we present an application of 2.5D and 3D ray‐Born migration/inversion to several classical data sets from the SEG/EAGE overthrust experiment. The method is formulated as a linearized inversion of the scattered wavefield. The method allows quantitative estimates of short wavelength components of the velocity model. First, we apply a 3D migration/inversion formula formerly developed for marine acquisitions to the swath data set. The migrated sections exhibit significant amplitude artifacts and acquisition footprints, also revealed by the shape of the local spatial resolution filters. From the analysis of these spatial resolution filters, we propose a new formula significantly improving the migrated dip section. We also present 3D migrated results for the strike section and a small 3D target containing a channel. Finally, the applications demonstrate, that the ray+Born migration formula must be adapted to the acquisition geometry to obtain reliable estimates of the true amplitude of the model perturbations. This adaptation is relatively straightforward in the frame of the ray+Born formalism and can be guided by the analysis of the resolution operator.

Structurally constrained impedance inversion

2016 ◽  
Vol 4 (4) ◽  
pp. T577-T589 ◽  
Author(s):  
Haitham Hamid ◽  
Adam Pidlisecky
Keyword(s):  
Objective Function ◽  
Seismic Data ◽  
Field Data ◽  
A Priori ◽  
Structural Constraints ◽  
Inversion Algorithm ◽  
Data Set ◽  
Impedance Inversion ◽  
Single Trace ◽  
Complex Geology

In complex geology, the presence of highly dipping structures can complicate impedance inversion. We have developed a structurally constrained inversion in which a computationally well-behaved objective function is minimized subject to structural constraints. This approach allows the objective function to incorporate structural orientation in the form of dips into our inversion algorithm. Our method involves a multitrace impedance inversion and a rotation of an orthogonal system of derivative operators. Local dips used to constrain the derivative operators were estimated from migrated seismic data. In addition to imposing structural constraints on the inversion model, this algorithm allows for the inclusion of a priori knowledge from boreholes. We investigated this algorithm on a complex synthetic 2D model as well as a seismic field data set. We compared the result obtained with this approach with the results from single trace-based inversion and laterally constrained inversion. The inversion carried out using dip information produces a model that has higher resolution that is more geologically realistic compared with other methods.

scholarly journals Characterization of the umbra–penumbra boundary by the vertical component of the magnetic field

2020 ◽  
Vol 638 ◽  
pp. A25
Author(s):  
P. Lindner ◽  
R. Schlichenmaier ◽  
N. Bello González
Keyword(s):  
Magnetic Field ◽  
Near Infrared ◽  
Statistical Study ◽  
Vertical Component ◽  
Fundamental Property ◽  
Field Vector ◽  
Data Sets ◽  
Data Set ◽  
Geometric Difference ◽  
The Magnetic Field

Context. The vertical component of the magnetic field was found to reach a constant value at the boundary between penumbra and umbra of stable sunspots in a recent statistical study of Hinode/SP data. This finding has profound implications as it can serve as a criterion to distinguish between fundamentally different magneto-convective modes operating in the sun. Aims. The objective of this work is to verify the existence of a constant value for the vertical component of the magnetic field (B⊥) at the boundary between umbra and penumbra from ground-based data in the near-infrared wavelengths and to determine its value for the GREGOR Infrared Spectrograph (GRIS@GREGOR) data. This is the first statistical study on the Jurčák criterion with ground-based data, and we compare it with the results from space-based data (Hinode/SP and SDO/HMI). Methods. Eleven spectropolarimetric data sets from the GRIS@GREGOR slit-spectograph containing fully-fledged stable sunspots were selected from the GRIS archive. SIR inversions including a polarimetric straylight correction are used to produce maps of the magnetic field vector using the Fe I 15648 Å and 15662 Å lines. Averages of B⊥ along the contours between penumbra and umbra are analyzed for the 11 data sets. In addition, contours at the resulting B⊥const are drawn onto maps and compared to intensity contours. The geometric difference between these contours, ΔP, is calculated for each data set. Results. Averaged over the 11 sunspots, we find a value of B⊥const = (1787 ± 100) gauss. The difference from the values previously derived from Hinode/SP and SDO/HMI data is explained by instrumental differences and by the formation characteristics of the respective lines that were used. Contours at B⊥ = B⊥const and contours calculated in intensity maps match from a visual inspection and the geometric distance ΔP was found to be on the order of 2 pixels. Furthermore, the standard deviation between different data sets of averages along umbra–penumbra contours is smaller for B⊥ than for B∥ by a factor of 2.4. Conclusions. Our results provide further support to the Jurčák criterion with the existence of an invariable value B⊥const at the umbra–penumbra boundary. This fundamental property of sunspots can act as a constraining parameter in the calibration of analysis techniques that calculate magnetic fields. It also serves as a requirement for numerical simulations to be realistic. Furthermore, it is found that the geometric difference, ΔP, between intensity contours and contours at B⊥ = B⊥const acts as an index of stability for sunspots.

Land FWI in the Delaware Basin, west Texas: A case study

2020 ◽  
Vol 39 (5) ◽  
pp. 324-331
Author(s):  
Gary Murphy ◽  
Vanessa Brown ◽  
Denes Vigh
Keyword(s):  
Seismic Data ◽  
Waveform Inversion ◽  
Surface Velocity ◽  
High Signal ◽  
Data Sets ◽  
Data Set ◽  
Low Frequencies ◽  
Near Surface ◽  
Delaware Basin ◽  
West Texas

As part of a wide-reaching full-waveform inversion (FWI) research program, FWI is applied to an onshore seismic data set collected in the Delaware Basin, west Texas. FWI is routinely applied on typical marine data sets with high signal-to-noise ratio (S/N), relatively good low-frequency content, and reasonably long offsets. Land seismic data sets, in comparison, present significant challenges for FWI due to low S/N, a dearth of low frequencies, and limited offsets. Recent advancements in FWI overcome limitations due to poor S/N and low frequencies making land FWI feasible to use to update the shallow velocities. The chosen area has contrasting and variable near-surface conditions providing an excellent test data set on which to demonstrate the workflow and its challenges. An acoustic FWI workflow is used to update the near-surface velocity model in order to improve the deeper image and simultaneously help highlight potential shallow drilling hazards.

Time-lapse acoustic impedance variations during CO2 injection in Weyburn oilfield, Canada

Geophysics ◽  
2019 ◽  
Vol 85 (1) ◽  
pp. M1-M13 ◽  
Author(s):  
Yichuan Wang ◽  
Igor B. Morozov
Keyword(s):  
Oil Recovery ◽  
Seismic Data ◽  
Acoustic Impedance ◽  
Time Lapse ◽  
Time Shift ◽  
Co2 Injection ◽  
Inversion Method ◽  
Data Sets ◽  
Storage Site ◽  
Data Set

For seismic monitoring injected fluids during enhanced oil recovery or geologic [Formula: see text] sequestration, it is useful to measure time-lapse (TL) variations of acoustic impedance (AI). AI gives direct connections to the mechanical and fluid-related properties of the reservoir or [Formula: see text] storage site; however, evaluation of its subtle TL variations is complicated by the low-frequency and scaling uncertainties of this attribute. We have developed three enhancements of TL AI analysis to resolve these issues. First, following waveform calibration (cross-equalization) of the monitor seismic data sets to the baseline one, the reflectivity difference was evaluated from the attributes measured during the calibration. Second, a robust approach to AI inversion was applied to the baseline data set, based on calibration of the records by using the well-log data and spatially variant stacking and interval velocities derived during seismic data processing. This inversion method is straightforward and does not require subjective selections of parameterization and regularization schemes. Unlike joint or statistical inverse approaches, this method does not require prior models and produces accurate fitting of the observed reflectivity. Third, the TL AI difference is obtained directly from the baseline AI and reflectivity difference but without the uncertainty-prone subtraction of AI volumes from different seismic vintages. The above approaches are applied to TL data sets from the Weyburn [Formula: see text] sequestration project in southern Saskatchewan, Canada. High-quality baseline and TL AI-difference volumes are obtained. TL variations within the reservoir zone are observed in the calibration time-shift, reflectivity-difference, and AI-difference images, which are interpreted as being related to the [Formula: see text] injection.

Sign in / Sign up

Export Citation Format

Share Document