Using orthorhombic depth imaging to image fractures in tight reservoirs of the RDG Field, Barmer Basin, India

2019 ◽  
Vol 38 (4) ◽  
pp. 268-273
Author(s):  
Maheswara Phani ◽  
Sushobhan Dutta ◽  
Kondal Reddy ◽  
Sreedurga Somasundaram

Raageshwari Deep Gas (RDG) Field is situated in the southern part of the Barmer Basin in Rajasthan, India, at a depth of 3000 m. With both clastic and volcanic lithologies, the main reservoirs are tight, and hydraulic fracturing is required to enhance productivity, especially to improve permeability through interaction of induced fractures with natural fractures. Therefore, optimal development of the RDG Field reservoirs requires characterization of faults and natural fractures. To address this challenge, a wide-azimuth 3D seismic data set over the RDG Field was processed to sharply define faults and capture anisotropy related to open natural fractures. Anisotropy was indicated by the characteristic sinusoidal nature of gather reflection events processed using conventional tilted transverse imaging (TTI); accordingly, we used orthorhombic imaging to correct for these, to quantify fracture-related anisotropy, and to yield a more correct subsurface image. During prestack depth migration (PSDM) processing of the RDG data, TTI and orthorhombic velocity modeling was done with azimuthal sectoring of reflection arrivals. The resultant PSDM data using this velocity model show substantial improvement in image quality and vertical resolution at the reservoir level compared to vintage seismic data. The improved data quality enabled analysis of specialized seismic attributes like curvature and thinned fault likelihood for more reliable characterization of faults and fractures. These attributes delineate the location and distribution of probable fracture networks within the volcanic reservoirs. Interpreted subtle faults, associated with fracture zones, were validated with microseismic, production, and image log data.

Geophysics ◽  
2003 ◽  
Vol 68 (6) ◽  
pp. 1782-1791 ◽  
Author(s):  
M. Graziella Kirtland Grech ◽  
Don C. Lawton ◽  
Scott Cheadle

We have developed an anisotropic prestack depth migration code that can migrate either vertical seismic profile (VSP) or surface seismic data. We use this migration code in a new method for integrated VSP and surface seismic depth imaging. Instead of splicing the VSP image into the section derived from surface seismic data, we use the same migration algorithm and a single velocity model to migrate both data sets to a common output grid. We then scale and sum the two images to yield one integrated depth‐migrated section. After testing this method on synthetic surface seismic and VSP data, we applied it to field data from a 2D surface seismic line and a multioffset VSP from the Rocky Mountain Foothills of southern Alberta, Canada. Our results show that the resulting integrated image exhibits significant improvement over that obtained from (a) the migration of either data set alone or (b) the conventional splicing approach. The integrated image uses the broader frequency bandwidth of the VSP data to provide higher vertical resolution than the migration of the surface seismic data. The integrated image also shows enhanced structural detail, since no part of the surface seismic section is eliminated, and good event continuity through the use of a single migration–velocity model, obtained by an integrated interpretation of borehole and surface seismic data. This enhanced migrated image enabled us to perform a more robust interpretation with good well ties.


Geophysics ◽  
2021 ◽  
pp. 1-45
Author(s):  
Guofeng Liu ◽  
Xiaohong Meng ◽  
Johanes Gedo Sea

Seismic reflection is a proven and effective method commonly used during the exploration of deep mineral deposits in Fujian, China. In seismic data processing, rugged depth migration based on wave-equation migration can play a key role in handling surface fluctuations and complex underground structures. Because wave-equation migration in the shot domain cannot output offset-domain common-image gathers in a straightforward way, the use of traditional tools for updating the velocity model and improving image quality can be quite challenging. To overcome this problem, we employed the attribute migration method. This worked by sorting the migrated stack results for every single-shot gather into the offset gathers. The value of the offset that corresponded to each image point was obtained from the ratio of the original migration results to the offset-modulated shot-data migration results. A Gaussian function was proposed to map every image point to a certain range of offsets. This helped improve the signal-to-noise ratio, which was especially important in handing low quality seismic data obtained during mineral exploration. Residual velocity analysis was applied to these gathers to update the velocity model and improve image quality. The offset-domain common-image gathers were also used directly for real mineral exploration seismic data with rugged depth migration. After several iterations of migration and updating the velocity, the proposed procedure achieved an image quality better than the one obtained with the initial velocity model. The results can help with the interpretation of thrust faults and deep deposit exploration.


Geophysics ◽  
1993 ◽  
Vol 58 (6) ◽  
pp. 873-882 ◽  
Author(s):  
Roelof Jan Versteeg

To get a correct earth image from seismic data acquired over complex structures it is essential to use prestack depth migration. A necessary condition for obtaining a correct image is that the prestack depth migration is done with an accurate velocity model. In cases where we need to use prestack depth migration determination of such a model using conventional methods does not give satisfactory results. Thus, new iterative methods for velocity model determination have been developed. The convergence of these methods can be accelerated by defining constraints on the model in such a way that the method only looks for those components of the true earth velocity field that influence the migrated image. In order to determine these components, the sensitivity of the prestack depth migration result to the velocity model is examined using a complex synthetic data set (the Marmousi data set) for which the exact model is known. The images obtained with increasingly smoothed versions of the true model are compared, and it is shown that the minimal spatial wavelength that needs to be in the model to obtain an accurate depth image from the data set is of the order of 200 m. The model space that has to be examined to find an accurate velocity model from complex seismic data can thus be constrained. This will increase the speed and probability of convergence of iterative velocity model determination methods.


2019 ◽  
Vol 38 (11) ◽  
pp. 872a1-872a9 ◽  
Author(s):  
Mauricio Araya-Polo ◽  
Stuart Farris ◽  
Manuel Florez

Exploration seismic data are heavily manipulated before human interpreters are able to extract meaningful information regarding subsurface structures. This manipulation adds modeling and human biases and is limited by methodological shortcomings. Alternatively, using seismic data directly is becoming possible thanks to deep learning (DL) techniques. A DL-based workflow is introduced that uses analog velocity models and realistic raw seismic waveforms as input and produces subsurface velocity models as output. When insufficient data are used for training, DL algorithms tend to overfit or fail. Gathering large amounts of labeled and standardized seismic data sets is not straightforward. This shortage of quality data is addressed by building a generative adversarial network (GAN) to augment the original training data set, which is then used by DL-driven seismic tomography as input. The DL tomographic operator predicts velocity models with high statistical and structural accuracy after being trained with GAN-generated velocity models. Beyond the field of exploration geophysics, the use of machine learning in earth science is challenged by the lack of labeled data or properly interpreted ground truth, since we seldom know what truly exists beneath the earth's surface. The unsupervised approach (using GANs to generate labeled data)illustrates a way to mitigate this problem and opens geology, geophysics, and planetary sciences to more DL applications.


2017 ◽  
Vol 5 (3) ◽  
pp. SJ81-SJ90 ◽  
Author(s):  
Kainan Wang ◽  
Jesse Lomask ◽  
Felix Segovia

Well-log-to-seismic tying is a key step in many interpretation workflows for oil and gas exploration. Synthetic seismic traces from the wells are often manually tied to seismic data; this process can be very time consuming and, in some cases, inaccurate. Automatic methods, such as dynamic time warping (DTW), can match synthetic traces to seismic data. Although these methods are extremely fast, they tend to create interval velocities that are not geologically realistic. We have described the modification of DTW to create a blocked dynamic warping (BDW) method. BDW generates an automatic, optimal well tie that honors geologically consistent velocity constraints. Consequently, it results in updated velocities that are more realistic than other methods. BDW constrains the updated velocity to be constant or linearly variable inside each geologic layer. With an optimal correlation between synthetic seismograms and surface seismic data, this algorithm returns an automatically updated time-depth curve and an updated interval velocity model that still retains the original geologic velocity boundaries. In other words, the algorithm finds the optimal solution for tying the synthetic to the seismic data while restricting the interval velocity changes to coincide with the initial input blocking. We have determined the application of the BDW technique on a synthetic data example and field data set.


Geophysics ◽  
2011 ◽  
Vol 76 (4) ◽  
pp. S157-S164 ◽  
Author(s):  
Robert Sun ◽  
George A. McMechan

We have extended prestack parsimonious Kirchhoff depth migration for 2D, two-component, reflected elastic seismic data for a P-wave source recorded at the earth’s surface. First, we separated the P-to-P reflected (PP-) waves and P-to-S converted (PS-) waves in an elastic common-source gather into P-wave and S-wave seismograms. Next, we estimated source-ray parameters (source p values) and receiver-ray parameters (receiver p values) for the peaks and troughs above a threshold amplitude in separated P- and S-wavefields. For each PP and PS reflection, we traced (1) a source ray in the P-velocity model in the direction of the emitted ray angle (determined by the source p value) and (2) a receiver ray in the P- or S-velocity model back in the direction of the emergent PP- or PS-wave ray angle (determined by the PP- or PS-wave receiver p value), respectively. The image-point position was adjusted from the intersection of the source and receiver rays to the point where the sum of the source time and receiver-ray time equaled the two-way traveltime. The orientation of the reflector surface was determined to satisfy Snell’s law at the intersection point. The amplitude of a P-wave (or an S-wave) was distributed over the first Fresnel zone along the reflector surface in the P- (or S-) image. Stacking over all P-images of the PP-wave common-source gathers gave the stacked P-image, and stacking over all S-images of the PS-wave common-source gathers gave the stacked S-image. Synthetic examples showed acceptable migration quality; however, the images were less complete than those produced by scalar reverse-time migration (RTM). The computing time for the 2D examples used was about 1/30 of that for scalar RTM of the same data.


Geophysics ◽  
1998 ◽  
Vol 63 (2) ◽  
pp. 546-556 ◽  
Author(s):  
Herman Chang ◽  
John P. VanDyke ◽  
Marcelo Solano ◽  
George A. McMechan ◽  
Duryodhan Epili

Portable, production‐scale 3-D prestack Kirchhoff depth migration software capable of full‐volume imaging has been successfully implemented and applied to a six‐million trace (46.9 Gbyte) marine data set from a salt/subsalt play in the Gulf of Mexico. Velocity model building and updates use an image‐driven strategy and were performed in a Sun Sparc environment. Images obtained by 3-D prestack migration after three velocity iterations are substantially better focused and reveal drilling targets that were not visible in images obtained from conventional 3-D poststack time migration. Amplitudes are well preserved, so anomalies associated with known reservoirs conform to the petrophysical predictions. Prototype development was on an 8-node Intel iPSC860 computer; the production version was run on an 1824-node Intel Paragon computer. The code has been successfully ported to CRAY (T3D) and Unix workstation (PVM) environments.


Geophysics ◽  
2014 ◽  
Vol 79 (3) ◽  
pp. S105-S111 ◽  
Author(s):  
Sheng Xu ◽  
Feng Chen ◽  
Bing Tang ◽  
Gilles Lambare

When using seismic data to image complex structures, the reverse time migration (RTM) algorithm generally provides the best results when the velocity model is accurate. With an inexact model, moveouts appear in common image gathers (CIGs), which are either in the surface offset domain or in subsurface angle domain; thus, the stacked image is not well focused. In extended image gathers, the strongest energy of a seismic event may occur at non-zero-lag in time-shift or offset-shift gathers. Based on the operation of RTM images produced by the time-shift imaging condition, the non-zero-lag time-shift images exhibit a spatial shift; we propose an approach to correct them by a second pass of migration similar to zero-offset depth migration; the proposed approach is based on the local poststack depth migration assumption. After the proposed second-pass migration, the time-shift CIGs appear to be flat and can be stacked. The stack enhances the energy of seismic events that are defocused at zero time lag due to the inaccuracy of the model, even though the new focused events stay at the previous positions, which might deviate from the true positions of seismic reflection. With the stack, our proposed approach is also able to attenuate the long-wavelength RTM artifacts. In the case of tilted transverse isotropic migration, we propose a scheme to defocus the coherent noise, such as migration artifacts from residual multiples, by applying the original migration velocity model along the symmetry axis but with different anisotropic parameters in the second pass of migration. We demonstrate that our approach is effective to attenuate the coherent noise at subsalt area with two synthetic data sets and one real data set from the Gulf of Mexico.


Geophysics ◽  
2003 ◽  
Vol 68 (4) ◽  
pp. 1357-1370 ◽  
Author(s):  
Stéphane Operto ◽  
Gilles Lambaré ◽  
Pascal Podvin ◽  
Philippe Thierry

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.


2005 ◽  
Vol 45 (1) ◽  
pp. 421
Author(s):  
P. Bocca ◽  
L. Fava ◽  
E. Stolf

3D pre-stack depth migration (PSDM) reprocessing was conducted in 2003 on a portion of the Onnia 3D seismic cube, located in exploration permit AC/P-21, Timor Sea.The main objective of the reprocessing was to obtain the best seismic depth image and the most realistic structural reconstruction of the sub-surface to mitigate the risk factors associated with trap definition (trap retention and trap efficiency). This represents one of the main challenges for oil exploration in the area.The 3D PSDM methodology was chosen as the most appropriate imaging tool to define the correct sub-surface geometry and fault imaging through the use of an appropriate velocity field. An integrated approach to building the final velocity model was adopted, with a substantial contribution from the regional geological model.Several examples are given to demonstrate that the 3D PSDM reprocessing significantly improved the seismic image and thus the confidence in the interpretation, contributing strongly to the definition of the exploration targets.The interpretation of the new seismic data has resulted in a new structural picture in which higher confidence in seismic imaging has improved fault correlation. This has enabled better structural definition at the Middle Jurassic Plover Formation level that has reduced the complexity of the large Vesta Prospect, in the centre of the Swan Graben to the northwest of East Swan–1. Improved understanding of the fault reactivation mechanism and the structural elements of the trap (trap integrity) were eventually incorporated in the prospect risking.In the Swan Graben 3D PSDM has proved to be a very powerful instrument capable of producing significant impact on the exploration even in an area with a complex geological setting and a fairly poor seismic data quality.


Sign in / Sign up

Export Citation Format

Share Document