scholarly journals Seismic characterization of multiple fracture sets at Rulison Field, Colorado

Geophysics ◽  
2007 ◽  
Vol 72 (2) ◽  
pp. B19-B30 ◽  
Author(s):  
Ivan Vasconcelos ◽  
Vladimir Grechka
Keyword(s):  
Velocity Model ◽  
Multiple Fracture ◽  
Fracture Characterization ◽  
Seismic Reflection Data ◽  
Characterization Methods ◽  
S Waves ◽  
Reflection Data ◽  
Seismic Characterization ◽  
Inversion Procedure

Conventional fracture-characterization methods assume the presence of a single set of aligned, vertical cracks in the subsur-face. We relax this assumption and demonstrate the feasibility of seismic characterization of multiple fracture sets. Our technique relies on recent numerical findings indicating that multiple, differently oriented, possibly intersecting planar cracks embedded in an otherwise isotropic host rock result in a nearly orthorhombic (or orthotropic) effective medium. Here, the governing parameters of crack-induced orthotropy are estimated from 3D, wide-azimuth, multicomponent seismic reflection data acquired over the tight-gas Rulison Field in Colorado. We translate strong azimuthal variations of the normal-moveout velocities intointerval crack densities, fracture orientations, type of fluid infill, and velocities of P- and S-waves in an unfractured rock. Our inversion procedure identifies a set of cracks aligned in approximately west northwest-east southeast direction in the western part of the study area and multiple, likely intersecting fractures in its eastern part. We validate both our underlying theoretical model and the obtained estimates by two independent measurements: (1) the estimated fluid-infill parameter indicates dry cracks as expected for the gas-producing sandstones at Rulison; and (2) the obtained crack orientations are supported by well observations. As a by-product of fracture characterization, we build an anisotropic velocity model of the Rulison reservoir which, we believe, is the first orthorhombic velocity field constructed from surface seismic data.

scholarly journals Study on the limitations of travel-time inversion applied to sub-basalt imaging

Solid Earth ◽  
2013 ◽  
Vol 4 (2) ◽  
pp. 543-554 ◽  
Author(s):  
I. Flecha ◽  
R. Carbonell ◽  
R. W. Hobbs
Keyword(s):  
Oil Industry ◽  
Velocity Model ◽  
Time Inversion ◽  
Wide Angle ◽  
Seismic Reflection Data ◽  
Additional Information ◽  
Velocity Models ◽  
Reflection Data ◽  
Travel Time Inversion ◽  
Tomographic Inversions

Abstract. The difficulties of seismic imaging beneath high velocity structures are widely recognised. In this setting, theoretical analysis of synthetic wide-angle seismic reflection data indicates that velocity models are not well constrained. A two-dimensional velocity model was built to simulate a simplified structural geometry given by a basaltic wedge placed within a sedimentary sequence. This model reproduces the geological setting in areas of special interest for the oil industry as the Faroe-Shetland Basin. A wide-angle synthetic dataset was calculated on this model using an elastic finite difference scheme. This dataset provided travel times for tomographic inversions. Results show that the original model can not be completely resolved without considering additional information. The resolution of nonlinear inversions lacks a functional mathematical relationship, therefore, statistical approaches are required. Stochastic tests based on Metropolis techniques support the need of additional information to properly resolve sub-basalt structures.

A multioffset, three‐component VSP study in the Sudbury Basin

Geophysics ◽  
1995 ◽  
Vol 60 (2) ◽  
pp. 341-353 ◽  
Author(s):  
Xiao‐Gui Miao ◽  
Wooil M. Moon ◽  
B. Milkereit
Keyword(s):  
Data Processing ◽  
Seismic Reflection ◽  
Velocity Structure ◽  
Least Square ◽  
Seismic Reflection Data ◽  
S Waves ◽  
Velocity Models ◽  
Reflection Data ◽  
Wavefield Separation ◽  
Vsp Data

A multioffset, three‐component vertical seismic profiling (VSP) experiment was carried out in the Sudbury Basin, Ontario, as a part of the LITHOPROBE Sudbury Transect. The main objectives were determination of the shallow velocity structure in the middle of the Sudbury Basin, development of an effective VSP data processing flow, correlation of the VSP survey results with the surface seismic reflection data, and demonstration of the usefulness of the VSP method in a crystalline rock environment. The VSP data processing steps included rotation of the horizontal component data, traveltime inversion for velocity analysis, Radon transform for wavefield separation, and preliminary analysis of shear‐wave data. After wavefield separation, the flattened upgoing wavefields for both P‐waves and S‐waves display consistent reflection events from three depth levels. The VSP-CDP transformed section and corridor stacked section correlate well with the high‐resolution surface reflection data. In addition to obtaining realistic velocity models for both P‐ and S‐waves through least‐square inversion and synthetic seismic modeling for the Chelmsford area, the VSP experiment provided an independent estimation for the reflector dip using three component hodogram analysis, which indicates that the dip of the contact between the Chelmsford and Onwatin formations, at an approximate depth of 380 m in the Chelmsford borehole, is approximately 10.5° southeast. This study demonstrates that multioffset, three‐component VSP experiments can provide important constraints and auxiliary information for shallow crustal seismic studies in crystalline terrain. Thus, the VSP technique bridges the gap between the surface seismic‐reflection technique and well‐log surveys.

Precision analysis of first‐break times in grid models

Geophysics ◽  
1998 ◽  
Vol 63 (3) ◽  
pp. 1062-1065 ◽  
Author(s):  
Thomas Gruber ◽  
Stewart A. Greenhalgh
Keyword(s):  
Numerical Models ◽  
Grid Point ◽  
Tomographic Reconstruction ◽  
Velocity Model ◽  
Model Parameters ◽  
Seismic Reflection Data ◽  
Arrival Times ◽  
Velocity Models ◽  
Reflection Data ◽  
Inversion Techniques

Rectangular grid velocity models and their derivatives are widely used in geophysical inversion techniques. Specifically, seismic tomographic reconstruction techniques, whether they be based on raypath methods (Bregman et al., 1989; Moser, 1991; Schneider et al., 1992; Cao and Greenhalgh, 1993; Zhou, 1993) or full wave equation methods (Vidale, 1990; Qin and Schuster, 1993; Cao and Greenhalgh, 1994) for calculating synthetic arrival times, involve propagation through a grid model. Likewise, migration of seismic reflection data, using asymptotic ray theory or finite difference/pseudospectral methods (Stolt and Benson, 1986; Zhe and Greenhalgh, 1997) involve assigning traveltimes to upward and downward propagating waves at every grid point in the model. The traveltimes in both cases depend on the grid specification. However, the precision level of such numerical models and their dependence on the model parameters is often unknown. In this paper, we describe a two‐dimensional velocity model and derive an error bound for first‐break times calculated with such a model. The analysis provides clear guidelines for grid specifications.

scholarly journals Reconstructing the primary reflections in seismic data by Marchenko redatuming and convolutional interferometry

Geophysics ◽  
2016 ◽  
Vol 81 (2) ◽  
pp. Q15-Q26 ◽  
Author(s):  
Giovanni Angelo Meles ◽  
Kees Wapenaar ◽  
Andrew Curtis
Keyword(s):  
Velocity Model ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Data Set ◽  
Seismic Reflection Data ◽  
Surface Reflection ◽  
Reflection Data ◽  
Reference Velocity ◽  
Time Migration ◽  
Recorded Data

State-of-the-art methods to image the earth’s subsurface using active-source seismic reflection data involve reverse time migration. This and other standard seismic processing methods such as velocity analysis provide best results only when all waves in the data set are primaries (waves reflected only once). A variety of methods are therefore deployed as processing to predict and remove multiples (waves reflected several times); however, accurate removal of those predicted multiples from the recorded data using adaptive subtraction techniques proves challenging, even in cases in which they can be predicted with reasonable accuracy. We present a new, alternative strategy to construct a parallel data set consisting only of primaries, which is calculated directly from recorded data. This obviates the need for multiple prediction and removal methods. Primaries are constructed by using convolutional interferometry to combine the first-arriving events of upgoing and direct-wave downgoing Green’s functions to virtual receivers in the subsurface. The required upgoing wavefields to virtual receivers are constructed by Marchenko redatuming. Crucially, this is possible without detailed models of the earth’s subsurface reflectivity structure: Similar to the most migration techniques, the method only requires surface reflection data and estimates of direct (nonreflected) arrivals between the virtual subsurface sources and the acquisition surface. We evaluate the method on a stratified synclinal model. It is shown to be particularly robust against errors in the reference velocity model used and to improve the migrated images substantially.

scholarly journals Study on the limitations of traveltime inversion in the presence of extreme velocity anomalies

2013 ◽  
Vol 5 (1) ◽  
pp. 189-226
Author(s):  
I. Flecha ◽  
R. Carbonell ◽  
R. W. Hobbs
Keyword(s):  
Oil Industry ◽  
Velocity Model ◽  
Wide Angle ◽  
Seismic Reflection Data ◽  
Traveltime Inversion ◽  
Additional Information ◽  
Velocity Models ◽  
Reflection Data ◽  
Velocity Anomalies ◽  
Tomographic Inversions

Abstract. The difficulties of seismic imaging beneath high velocity structures are widely recognised. In this setting, theoretical analysis of synthetic wide-angle seismic reflection data indicates that velocity models are not well constrained. A two-dimensional velocity model was built to simulate a simplified structural geometry given by a basaltic wedge placed within a sedimentary sequence. This model reproduces the geological setting in areas of special interest for the oil industry as the Faroe-Shetland Basin. A wide-angle synthetic dataset was calculated on this model using an elastic finite difference scheme. This dataset provided travel times for tomographic inversions. Results show that the original model can not be completely resolved without considering additional information. The resolution of nonlinear inversions lacks a functional mathematical relationship, therefore, statistical approaches are required. Stochastical tests based on Metropolis techniques support the need of additional information to properly resolve subbasalt structures.

Application of waveform tomography to marine seismic reflection data from the Queen Charlotte Basin of western Canada

Geophysics ◽  
2011 ◽  
Vol 76 (2) ◽  
pp. B55-B70 ◽  
Author(s):  
E. M. Takam Takougang ◽  
A. J. Calvert
Keyword(s):  
Seismic Reflection ◽  
Field Data ◽  
Velocity Model ◽  
P Wave ◽  
Low Velocity ◽  
Seismic Line ◽  
Seismic Reflection Data ◽  
Reflection Data ◽  
Sonic Log ◽  
Waveform Tomography

To obtain a higher resolution quantitative P-wave velocity model, 2D waveform tomography was applied to seismic reflection data from the Queen Charlotte sedimentary basin off the west coast of Canada. The forward modeling and inversion were implemented in the frequency domain using the visco-acoustic wave equation. Field data preconditioning consisted of f-k filtering, 2D amplitude scaling, shot-to-shot amplitude balancing, and time windowing. The field data were inverted between 7 and 13.66 Hz, with attenuation introduced for frequencies ≥ 10.5 Hz to improve the final velocity model; two different approaches to sampling the frequencies were evaluated. The limited maximum offset of the marine data (3770 m) and the relatively high starting frequency (7 Hz) were the main challenges encountered during the inversion. An inversion strategy that successively recovered shallow-to-deep structures was designed to mitigate these issues. The inclusion of later arrivals in the waveform tomography resulted in a velocity model that extends to a depth of approximately 1200 m, twice the maximum depth of ray coverage in the ray-based tomography. Overall, there is a good agreement between the velocity model and a sonic log from a well on the seismic line, as well as between modeled shot gathers and field data. Anomalous zones of low velocity in the model correspond to previously identified faults or their upward continuation into the shallow Pliocene section where they are not readily identifiable in the conventional migration.

Nonlinear inversion of seismic reflection data in a laterally invariant medium

Geophysics ◽  
1990 ◽  
Vol 55 (3) ◽  
pp. 284-292 ◽  
Author(s):  
A. Pica ◽  
J. P. Diet ◽  
A. Tarantola
Keyword(s):  
Conjugate Gradient ◽  
Velocity Model ◽  
Stratified Medium ◽  
Single Shot ◽  
Nonlinear Inversion ◽  
Reverse Time ◽  
Data Set ◽  
Seismic Reflection Data ◽  
Seismic Waveforms ◽  
Reflection Data

Interpretation of seismic waveforms can be expressed as an optimization problem based on a non‐linear least‐squares criterion to find the model which best explains the data. An initial model is corrected iteratively using a gradient method (conjugate gradient). At each iteration, computation of the direction of the model perturbation requires the forward propagation of the actual sources and the reverse‐time propagation of the residuals (misfit between the data and the synthetics); the two wave fields thus obtained are then correlated. An extra forward propagation is required to compute the amplitude of the perturbation along the conjugate‐gradient direction. The number of propagations to be simulated numerically in each iteration equals three times the number of shots. Since a 2-D finite‐difference code is employed to solve forward‐ and backward‐propagation problems, the method is general and can handle arbitrary 2-D source‐receiver configurations and lateral heterogeneities. Using conventional velocity analysis to derive an initial velocity model, the method is implemented on a real marine data set. The data set which has been selected corresponds approximately to a horizontally stratified medium. Consequently, a single‐shot gather has been used for inversion. In spite of some simplifying assumptions used for wave‐field propagation (acoustic approximation, point source), and using synthetics generated by a nearby sonic log to calibrate amplitudes, our final synthetics match the input data very well and the inversion result has clear similarities to the log.

Sign in / Sign up

Export Citation Format

Share Document