Improved imaging of the Strzelecki Formation by the reprocessing of 3D seismic reflection data: onshore Gippsland Basin, Australia

2015 ◽  
Vol 55 (2) ◽  
pp. 465
Author(s):  
Hadi Nourollah ◽  
Javad Aliemrani

The Wombat 3D seismic survey was recorded by Lakes Oil in 2008 following the drilling of the Wombat–3 well. The survey was aimed at identifying the structures better inside the Strzelecki Formation and the underlying Rintoul Creek Sandstone, and to give a better understanding of the provenance of the gas encountered in the first three Wombat wells and the oil encountered at depth in the Wombat–3 well. Seismic acquisition design and processing in this area are challenged by the presence of coal measures in the Latrobe Group, which overlies the zone of interest, absorbs most of the incident seismic energy and complicates seismic imaging through the generation of multiple reflection events. A novel approach, described as multi-line decomposition, modelling and synthesis, was developed to model critical imaging parameters and in particular the velocity model for midpoint stacking and migration. This approach is demonstrated through the reprocessing of the Wombat 3D seismic survey and is compared with a conventional 3D processing approach and original survey processing examples. The reprocessed seismic section provides a significantly improved image of the deeper structures, delineating a number of continuous reflectors in the Early Cretaceous formations. It also provides the opportunity to build a more accurate model of the basement.

Geophysics ◽  
1986 ◽  
Vol 51 (1) ◽  
pp. 12-19 ◽  
Author(s):  
James F. Mitchell ◽  
Richard J. Bolander

Subsurface structure can be mapped using refraction information from marine multichannel seismic data. The method uses velocities and thicknesses of shallow sedimentary rock layers computed from refraction first arrivals recorded along the streamer. A two‐step exploration scheme is described which can be set up on a personal computer and used routinely in any office. It is straightforward and requires only a basic understanding of refraction principles. Two case histories from offshore Peru exploration demonstrate the scheme. The basic scheme is: step (1) shallow sedimentary rock velocities are computed and mapped over an area. Step (2) structure is interpreted from the contoured velocity patterns. Structural highs, for instance, exhibit relatively high velocities, “retained” by buried, compacted, sedimentary rocks that are uplifted to the near‐surface. This method requires that subsurface structure be relatively shallow because the refracted waves probe to depths of one hundred to over one thousand meters, depending upon the seismic energy source, streamer length, and the subsurface velocity distribution. With this one requirement met, we used the refraction method over a wide range of sedimentary rock velocities, water depths, and seismic survey types. The method is particularly valuable because it works well in areas with poor seismic reflection data.


2020 ◽  
pp. 3024-3035
Author(s):  
Kamal K. Ali ◽  
Gassak F. Kadhim

This study includes structural and stratigraphic interpretation of 3D seismic reflection data for Zubair Formation (L. Cretaceous) within Al-Akhadeir area, southwestern Iraq (Karbala Governorate). Depending on the 3D seismic reflection interpretation process, and based on the synthetic seismogram  and well logs data, two horizons were identified and selected (top and base Zubair  reflectors). These horizons were followed up over the entire area in order to obtain structural and stratigraphic settings. TWT, depth, and velocity maps for the base and top Zubair Formation were constructed. From the interpretation of these maps and based on the seismic section, the study concluded that there are some enclosures that represent anticline in the NW of the horizon and syncline in the NE, while the nose structure appears in the middle of the horizon and trends N-S. The horizon represents a progradational with sigmoid configuration.. Other seismic structural phenomena were recognized in this part of the area, such as flat spot, down lap, and top lap, which give indicators of potential hydrocarbon accumulations


2020 ◽  
Vol 28 ◽  
pp. 1-27 ◽  
Author(s):  
David R. Cox ◽  
Paul C. Knutz ◽  
D. Calvin Campbell ◽  
John R. Hopper ◽  
Andrew M. W. Newton ◽  
...  

Abstract. A geohazard assessment workflow is presented that maximizes the use of 3D seismic reflection data to improve the safety and success of offshore scientific drilling. This workflow has been implemented for International Ocean Discovery Program (IODP) Proposal 909 that aims to core seven sites with targets between 300 and 1000 m below seabed across the north-western Greenland continental shelf. This glaciated margin is a frontier petroleum province containing potential drilling hazards that must be avoided during drilling. Modern seismic interpretation techniques are used to identify, map and spatially analyse seismic features that may represent subsurface drilling hazards, such as seabed structures, faults, fluids and challenging lithologies. These hazards are compared against the spatial distribution of stratigraphic targets to guide site selection and minimize risk. The 3D seismic geohazard assessment specifically advanced the proposal by providing a more detailed and spatially extensive understanding of hazard distribution that was used to confidently select eight new site locations, abandon four others and fine-tune sites originally selected using 2D seismic data. Had several of the more challenging areas targeted by this proposal only been covered by 2D seismic data, it is likely that they would have been abandoned, restricting access to stratigraphic targets. The results informed the targeted location of an ultra-high-resolution 2D seismic survey by minimizing acquisition in unnecessary areas, saving valuable resources. With future IODP missions targeting similarly challenging frontier environments where 3D seismic data are available, this workflow provides a template for geohazard assessments that will enhance the success of future scientific drilling.


2021 ◽  
Author(s):  
Siddharth Garia ◽  
Arnab Kumar Pal ◽  
Karangat Ravi ◽  
Archana M Nair

<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>


Solid Earth ◽  
2013 ◽  
Vol 4 (2) ◽  
pp. 543-554 ◽  
Author(s):  
I. Flecha ◽  
R. Carbonell ◽  
R. W. Hobbs

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.


Geophysics ◽  
1998 ◽  
Vol 63 (4) ◽  
pp. 1395-1407 ◽  
Author(s):  
Frank Büker ◽  
Alan G. Green ◽  
Heinrich Horstmeyer

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).


2020 ◽  
Vol 177 (5) ◽  
pp. 1039-1056
Author(s):  
Thomas B. Phillips ◽  
Craig Magee

Intraplate volcanism is widely distributed across the continents, but the controls on the 3D geometry and longevity of individual volcanic systems remain poorly understood. Geophysical data provide insights into magma plumbing systems, but, as a result of the relatively low resolution of these techniques, it is difficult to evaluate how magma transits highly heterogeneous continental interiors. We use borehole-constrained 2D seismic reflection data to characterize the 3D geometry of the Tuatara Volcanic Field located offshore New Zealand's South Island and investigate its relationship with the pre-existing structure. This c. 270 km2 field is dominated by a dome-shaped lava edifice, surrounded and overlain by c. 69 volcanoes and >70 sills emplaced over 40 myr from the Late Cretaceous to Early Eocene (c. 85–45 Ma). The Tuatara Volcanic Field is located above a basement terrane boundary represented by the Livingstone Fault; the recently active Auckland Volcanic Field is similarly located along-strike on North Island. We suggest that the Livingstone Fault controlled the location of the Tuatara Volcanic Field by producing relief at the base of the lithosphere, thereby focussing lithospheric detachment over c. 40 myr, and provided a pathway that facilitated the ascent of magma. We highlight how observations from ancient intraplate volcanic systems may inform our understanding of active intraplate volcanic systems, including the Auckland Volcanic Field.Supplementary material: Interpreted seismic section showing well control on stratigraphic interpretation is available at https://doi.org/10.6084/m9.figshare.c.5004464


Geophysics ◽  
1998 ◽  
Vol 63 (3) ◽  
pp. 1062-1065 ◽  
Author(s):  
Thomas Gruber ◽  
Stewart A. Greenhalgh

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.


Sign in / Sign up

Export Citation Format

Share Document