Source-rock seismic-velocity models: Gassmann versus Backus

Geophysics ◽  
2011 ◽  
Vol 76 (5) ◽  
pp. N37-N45 ◽  
Author(s):  
José M. Carcione ◽  
Hans B. Helle ◽  
Per Avseth

Source rocks are described by a porous transversely isotropic medium composed of illite and organic matter (kerogen, oil, and gas). The bulk modulus of the oil/gas mixture is calculated by using a model of patchy saturation. Then, the moduli of the kerogen/fluid mixture are obtained with the Kuster and Toksöz model, assuming that oil is the inclusion in a kerogen matrix. To obtain the seismic velocities of the shale, we used Backus averaging and Gassmann equations generalized to the anisotropic case with a solid-pore infill. In the latter case, the dry-rock elastic constants are calculated with a generalization of Krief equations to the anisotropic case. We considered 11 samples of the Bakken-shale data set, with a kerogen pore infill. The Backus model provides lower and upper bounds of the velocities, whereas the Krief/Gassmann model provides a good match to the data. Alternatively, we obtain the dry-rock elastic moduli by using the inverse Gassmann equation, instead of using Krief equations. Four cases out of 11 yielded physically unstable results. We also considered samples of the North Sea Kimmeridge shale. In this case, Backus performed as well as the Krief/Gassmann model. If there is gas and oil in the shale, we found that the wave velocities are relatively constant when the amount of kerogen is kept constant. Varying kerogen content implies significant velocity changes versus fluid (oil) saturation.

2006 ◽  
Vol 46 (1) ◽  
pp. 215 ◽  
Author(s):  
C.I. Uruski ◽  
B.D. Field ◽  
R. Funnell ◽  
C. Hollis ◽  
A. Nicol ◽  
...  

Oil production in the East Coast Basin began in the late 19th century from wildcat wells near oil seeps. By the mid-20th century, geology was being applied to oil exploration, but with little success. In the late 20th century, seismic techniques were added to the exploration arsenal and several gas discoveries were made. At each stage it was recognised that exploration in this difficult but tantalising basin required more information than was available. Continuing work by exploration companies, as well as by the Institute of Geological & Nuclear Sciences (GNS), has begun to reduce the risk of exploration. Source rocks have been identified and sophisticated thermal models show that petroleum is being generated and expelled from them as shown by numerous oil and gas seeps onshore. Many potential reservoir sequences have been recognised from outcrop studies and depositional models are being refined. All components of petroleum systems have been demonstrated to be present. The most important deficiency to date is the general lack of high-quality seismic data which would allow recognition of reservoir facies in the subsurface.During early 2005, Crown Minerals, the New Zealand government group charged with promoting and regulating oil and gas exploration, commissioned a high specification regional 2D survey intended to address some of the main data gaps in the offshore East Coast Basin. A broad grid was planned with several regional lines to be acquired with a 12,000 m streamer and infill lines to be acquired with a streamer 8,000 m long. It was expected that the long streamer would increase resolution of Paleogene and Cretaceous units. Several of the lines were actually acquired with a 4,000 m streamer due to unexpectedly high rates of unserviceability. The resulting 2,800 km data set consists of a series of northwest–southeast lines approximately orthogonal to the coast at a spacing of about 10 km as well as several long strike lines.GNS was contracted to produce a series of reports covering source rock distribution, a catalogue of reservoir rocks, a regional seismic interpretation, thermal models and structural reconstruction. The data package and reports are available free of charge to any interested exploration company to accompany the licensing round that was announced on 1 September 2005. The new data set has confirmed the existence of a large, little-deformed basin to the north of North Island and the Bay of Plenty; it has elucidated the complex structure of a large part of the East Coast Basin and has enabled generation of a general sequence stratigraphic model which assists in delineating reservoir targets. On 1 September 2005, the New Zealand government launched a licensing round covering about 43,000 km2 of the East Coast Basin, from the far offshore East Cape Ridge in the north to the northern Wairarapa coast in the south. Four blocks (I, J, K and L) were on offer for a competitive staged work programme bid, closing on 17 February 2006.


GeoArabia ◽  
2009 ◽  
Vol 14 (3) ◽  
pp. 199-228 ◽  
Author(s):  
Mohammad Faqira ◽  
Martin Rademakers ◽  
AbdulKader M. Afifi

ABSTRACT During the past decade, considerable improvements in the seismic imaging of the deeper Paleozoic section, along with data from new well penetrations, have significantly improved our understanding of the mid-Carboniferous deformational event. Because it occurred at the same time as the Hercynian Orogeny in Europe, North Africa and North America it has been commonly referred to by the same name in the Middle East. This was the main tectonic event during the late Paleozoic, which initiated or reactivated many of the N-trending block uplifts that underlie the major hydrocarbon accumulations in eastern Arabia. The nature of the Hercynian deformation away from these structural features was poorly understood due to inadequate seismic imaging and insufficient well control, along with the tectonic overprint of subsequent deformation events. Three Hercynian NE-trending arches are recognized in the Arabian Plate (1) the Levant Arch, which extended from Egypt to Turkey along the coast of the Mediterranean Sea, (2) the Al-Batin Arch, which extended from the Arabian Shield through Kuwait to Iran, and (3) the Oman-Hadhramaut Arch, which extended along the southeast coast of Oman and Yemen. These arches were initiated during the mid-Carboniferous Hercynian Orogeny, and persisted until they were covered unconformably by the Khuff Formation during the Late Permian. Two Hercynian basins separate these arches: the Nafud-Ma’aniya Basin in the north and Faydah-Jafurah Basin in the south. The pre-Hercynian Paleozoic section was extensively eroded over the arches, resulting in a major angular unconformity, but generally preserved within the basins. Our interpretation suggests that most of the Arabian Shield, except the western highlands along the Red Sea, is the exhumed part of the Al-Batin Arch. The Hercynian structural fabric of regional arches and basins continue in northern Africa, and in general appear to be oriented orthogonal to the old margin of the Gondwana continent. The Hercynian structure of arches and basins was partly obliterated by subsequent Mesozoic and Cenozoic tectonic events. In eastern Saudi Arabia, Qatar, and Kuwait, regional extension during the Triassic formed N-trending horsts and graben that cut across the NE-trending Hercynian mega-structures, which locally inverted them. Subsequent reactivation during the Cretaceous and Neogene resulted in additional growth of the N-trending structures. The Hercynian Arches had major impact on the Paleozoic hydrocarbon accumulations. The Silurian source rocks are generally preserved in the basins and eroded over the arches, which generally confined Silurian-sourced hydrocarbons either within the basins or along their flanks. Furthermore, the relict Hercynian paleo-topography generally confined the post-Hercynian continental clastics of the Unayzah Formation and equivalents to the Hercynian basins. These clastics contain the main Paleozoic oil and gas reservoirs, particularly along the basin margins where they overlie the sub-crop of the Silurian section with angular unconformity, thus juxtaposing reservoir and source rock.


2018 ◽  
Vol 36 (5) ◽  
pp. 1229-1244
Author(s):  
Xiao-Rong Qu ◽  
Yan-Ming Zhu ◽  
Wu Li ◽  
Xin Tang ◽  
Han Zhang

The Huanghua Depression is located in the north-centre of Bohai Bay Basin, which is a rift basin developed in the Mesozoic over the basement of the Huabei Platform, China. Permo-Carboniferous source rocks were formed in the Huanghua Depression, which has experienced multiple complicated tectonic alterations with inhomogeneous uplift, deformation, buried depth and magma effect. As a result, the hydrocarbon generation evolution of Permo-Carboniferous source rocks was characterized by discontinuity and grading. On the basis of a detailed study on tectonic-burial history, the paper worked on the burial history, heating history and hydrocarbon generation history of Permo-Carboniferous source rocks in the Huanghua Depression combined with apatite fission track testing and fluid inclusion analyses using the EASY% Ro numerical simulation. The results revealed that their maturity evolved in stages with multiple hydrocarbon generations. In this paper, we clarified the tectonic episode, the strength of hydrocarbon generation and the time–spatial distribution of hydrocarbon regeneration. Finally, an important conclusion was made that the hydrocarbon regeneration of Permo-Carboniferous source rocks occurred in the Late Cenozoic and the subordinate depressions were brought forward as advantage zones for the depth exploration of Permo-Carboniferous oil and gas in the middle-northern part of the Huanghua Depression, Bohai Bay Basin, China.


2008 ◽  
Vol 48 (1) ◽  
pp. 53 ◽  
Author(s):  
Chris Uruski ◽  
Callum Kennedy ◽  
Rupert Sutherland ◽  
Vaughan Stagpoole ◽  
Stuart Henrys

The East Coast of North Island, New Zealand, is the site of subduction of the Pacific below the Australian plate, and, consequently, much of the basin is highly deformed. An exception is the Raukumara Sub-basin, which forms the northern end of the East Coast Basin and is relatively undeformed. It occupies a marine plain that extends to the north-northeast from the northern coast of the Raukumara Peninsula, reaching water depths of about 3,000 m, although much of the sub-basin lies within the 2,000 m isobath. The sub-basin is about 100 km across and has a roughly triangular plan, bounded by an east-west fault system in the south. It extends about 300 km to the northeast and is bounded to the east by the East Cape subduction ridge and to the west by the volcanic Kermadec Ridge. The northern seismic lines reveal a thickness of around 8 km increasing to 12–13 km in the south. Its stratigraphy consists of a fairly uniformly bedded basal section and an upper, more variable unit separated by a wedge of chaotically bedded material. In the absence of direct evidence from wells and samples, analogies are drawn with onshore geology, where older marine Cretaceous and Paleogene units are separated from a Neogene succession by an allochthonous series of thrust slices emplaced around the time of initiation of the modern plate boundary. The Raukumara Sub-basin is not easily classified. Its location is apparently that of a fore-arc basin along an ocean-to-ocean collision zone, although its sedimentary fill must have been derived chiefly from erosion of the New Zealand land mass. Its relative lack of deformation introduces questions about basin formation and petroleum potential. Although no commercial discoveries have been made in the East Coast Basin, known source rocks are of marine origin and are commonly oil prone, so there is good potential for oil as well as gas in the basin. New seismic data confirm the extent of the sub-basin and its considerable sedimentary thickness. The presence of potential trapping structures and direct hydrocarbon indicators suggest that the Raukumara Sub-basin may contain large volumes of oil and gas.


2018 ◽  
Vol 58 (2) ◽  
pp. 871 ◽  
Author(s):  
Melissa Thompson ◽  
Fred Wehr ◽  
Jack Woodward ◽  
Jon Minken ◽  
Gino D'Orazio ◽  
...  

Commencing in 2014, Quadrant Energy and partners have undertaken an active exploration program in the Bedout Sub-basin with a 100% success rate, discovering four hydrocarbon accumulations with four wells. The primary exploration target in the basin, the Middle Triassic Lower Keraudren Formation, encompasses the reservoirs, source rocks and seals that have trapped hydrocarbons in a self-contained petroleum system. This petroleum system is older than the traditional plays on the North-West Shelf and before recent activity was very poorly understood and easily overlooked. Key reservoirs occur at burial depths of 3500–5500 m, deeper than many of the traditional plays on the North-West Shelf and exhibit variable reservoir quality. Oil and gas-condensate discovered in the first two wells, Phoenix South-1 and Roc-1, raised key questions on the preservation of effective porosity and productivity sufficient to support a commercial development. With the acquisition and detailed interpretation of 119 m of core over the Caley Member reservoir in Roc-2 and a successful drill stem test that was surface equipment constrained to 55 MMscf/d, the productive potential of this reservoir interval has been confirmed. The results of the exploration program to date, combined with acquisition of new 3D/2D seismic data, have enabled a deeper understanding of the potential of the Bedout Sub-basin. A detailed basin model has been developed and a large suite of prospects and leads are recognised across a family of hydrocarbon plays. Two key wells currently scheduled for 2018 (Phoenix South-3 and Dorado-1) will provide critical information about the scale of this opportunity.


2020 ◽  
Author(s):  
Rebecca Bell

<p>The discovery of slow slip events (SSEs) at subduction margins in the last two decades has changed our understanding of how stress is released at subduction zones. Fault slip is now viewed as a continuum of different slip modes between regular earthquakes and aseismic creep, and an appreciation of seismic hazard can only be realised by understanding the full spectrum of slip. SSEs may have the potential to trigger destructive earthquakes and tsunami on faults nearby, but whether this is possible and why SSEs occur at all are two of the most important questions in earthquake seismology today. Laboratory and numerical models suggest that slow slip can be spontaneously generated under conditions of very low effective stresses, facilitated by high pore fluid pressure, but it has also been suggested that variations in frictional behaviour, potentially caused by very heterogeneous fault zone lithology, may be required to promote slow slip.</p><p>Testing these hypotheses is difficult as it requires resolving rock properties at a high resolution many km below the seabed sometimes in km’s of water, where drilling is technically challenging and expensive. Traditional geophysical methods like travel-time tomography cannot provide fine-scale enough velocity models to probe the rock properties in fault zones specifically. In the last decade, however, computational power has improved to the point where 3D full-waveform inversion (FWI) methods make it possible to use the full wavefield rather than just travel times to produce seismic velocity models with a resolution an order of magnitude better than conventional models. Although the hydrocarbon industry have demonstrated many successful examples of 3D FWI the method requires extremely high density arrays of instruments, very different to the 2D transect data collection style which is still commonly employed at subduction zones.</p><p> The north Hikurangi subduction zone, New Zealand is special, as it hosts the world’s most well characterised shallow SSEs (<2 km to 15 km below the seabed).  This makes it an ideal location to collect 3D data optimally for FWI to resolve rock properties in the slow slip zone. In 2017-2018 an unprecedentedly large 3D experiment including 3D multi-channel seismic reflection, 99 ocean bottom seismometers and 194 onshore seismometers was conducted along the north Hikurangi margin in an 100 km x 15 km area, with an average 2 km instrument spacing. In addition, IODP Expeditions 372 and 375 collected logging-while drilling and core data, and deployed two bore-hole observatories to target slow slip in the same area. In this presentation I will introduce you to this world class 3D dataset and preliminary results, which will enable high resolution 3D models of physical properties to be made to bring slow slip processes into focus.  </p>


Geophysics ◽  
2005 ◽  
Vol 70 (5) ◽  
pp. U51-U65 ◽  
Author(s):  
Stig-Kyrre Foss ◽  
Bjørn Ursin ◽  
Maarten V. de Hoop

We present a method of reflection tomography for anisotropic elastic parameters from PP and PS reflection seismic data. The method is based upon the differential semblance misfit functional in scattering angle and azimuth (DSA) acting on common-image-point gathers (CIGs) to find fitting velocity models. The CIGs are amplitude corrected using a generalized Radon transform applied to the data. Depth consistency between the PP and PS images is enforced by penalizing any mis-tie between imaged key reflectors. The mis-tie is evaluated by means of map migration-demigration applied to the geometric information (times and slopes) contained in the data. In our implementation, we simplify the codepthing approach to zero-scattering-angle data only. The resulting measure is incorporated as a regularization in the DSA misfit functional. We then resort to an optimization procedure, restricting ourselves to transversely isotropic (TI) velocity models. In principle, depending on the available surface-offset range and orientation of reflectors in the subsurface, by combining the DSA with codepthing, the anisotropic parameters for TI models can be determined, provided the orientation of the symmetry axis is known. A proposed strategy is applied to an ocean-bottom-seismic field data set from the North Sea.


2021 ◽  
Author(s):  
Bryant Chow

<p><b>Seismic tomography is a powerful tool for understanding Earth structure. In New Zealand, velocity models derived using ray-based tomography have been used extensively to characterize the complex plate boundary between the Australian and Pacific plates. Advances in computational capabilities now allow us to improve these velocity models using adjoint tomography, an imaging method which minimizes differences between observed and simulated seismic waveforms. We undertake the first application of adjoint tomography in New Zealand to improve a ray-based New Zealand velocity model containing the Hikurangi subduction zone and the North Island of New Zealand.</b></p> <p>In support of this work we deployed the Broadband East Coast Network (BEACON), a temporary seismic network aimed at improving coverage of the New Zealand permanent network, along the east coast of the North Island. We concurrently develop an automated, open-source workflow for full-waveform inversion using spectral element and adjoint methods. We employ this tool to assess a candidate velocity model’s suitability for adjoint tomography. Using a 3D ray-based traveltime tomography model of New Zealand, we generate synthetic seismic waveforms for more than 10 000 source–receiver pairs and evaluate waveform misfits. We subsequently perform synthetic checkerboard inversions with a realistic New Zealand source–receiver distribution. Reasonable systematic time shifts and satisfactory checkerboard resolution in synthetic inversions indicate that the candidate model is appropriate as an initial model for adjoint tomography. This assessment also demonstrates the relative ease of use and reliability of the automated tools.</p> <p>We then undertake a large-scale adjoint tomography inversion for the North Island of New Zealand using up to 1 800 unique source–receiver pairs to fit waveforms with periods 4–30 s, relating to minimum waveform sensitivities on the order of 5 km. Overall, 60 geographically well-distributed earthquakes and as many as 88 broadband station locations are included. Using a nonlinear optimization algorithm, we undertake 28 model updates of Vp and Vs over six distinct inversion legs which progressively increase resolution. The total inversion incurred a computational cost of approximately 500 000 CPU-hours. The overall time shift between observed and synthetic seismograms is reduced, and updated velocities show as much as ±30% change with respect to initial values. A formal resolution analysis using point spread tests highlights that velocity changes are strongly resolved onland and directly offshore, at depths above 30 km, with low-amplitude changes (> 1%) observed down to 100 km depth. The most striking velocity changes coincide with areas related to the active Hikurangi subduction zone.</p> <p>We interpret the updated velocity model in terms of New Zealand tectonics and geology, and observe good agreement with known basement terranes, and major structural elements such as faults, sedimentary basins, broad-scale subduction related features. We recover increased spatial heterogeneity in seismic velocities along the strike of the Hikurangi subduction zone with respect to the initial model. Below the East Coast, we interpret two localized high-velocity anomalies as previously unidentified subducted seamounts. We corroborate this interpretation with other work, and discuss the implications of deeply subducted seamounts on slip behavior along the Hikurangi margin. In the Cook Strait we observe a low-velocity zone that we interpret as a deep sedimentary basin. Strong velocity gradients bounding this low-velocity zone support hypotheses of a structural boundary here separating the North and South Islands of New Zealand. In the central North Island, low-velocity anomalies are linked to surface geology, and we relate seismic velocities at depth to crustal magmatic activity below the Taupo Volcanic Zone.</p> <p>This new velocity model provides more accurate synthetic seismograms and additional constraints on enigmatic tectonic processes related to the North Island of New Zealand. Both the velocity model itself, and the underpinning methodological contributions, improve our ever-expanding understanding of the North Island of New Zealand, the Hikurangi subduction zone, and the broader Australian-Pacific plate boundary.</p>


Geophysics ◽  
2020 ◽  
Vol 85 (5) ◽  
pp. R437-R446
Author(s):  
Taylor Dahlke ◽  
Biondo Biondi ◽  
Robert Clapp

As oil and gas extraction becomes more advanced, deep-water exploration becomes increasingly focused on imaging near or under complex salt geology, which necessitates detailed velocity models with strong contrast interfaces. These interfaces can be elegantly tracked using the level sets of an implicit surface. One can invert for the velocity model that best fits the recorded data in a full-waveform inversion (FWI) style objective function by reparameterizing the model in terms of an implicit surface representation of the salt interface. With this parameterization of the FWI objective function, we find the Hessian and solve a conjugate gradient system for the Newton step at every nonlinear iteration. We sparsify the representation of the implicit surface using radial basis functions, which can hasten convergence of the inner inversion by reducing the number of model parameters. We have developed a guided inversion approach that embeds information about the certainty of different salt boundary regions by the initialization of the implicit surface slope at the salt interface. This can help guide the inversion away from perceived local minima. The results of testing this inversion workflow on a 3D Gulf of Mexico data set show that it can be a useful tool for refining salt models because the initial and final seismic images show clearer and more consistent features below the updated salt area.


Sign in / Sign up

Export Citation Format

Share Document