Seismic parameters for media with elliptical velocity dependencies

Geophysics ◽  
1982 ◽  
Vol 47 (12) ◽  
pp. 1621-1626 ◽  
Author(s):  
Bok S. Byun

The concept of wavefront curvature has been discussed extensively in the literature to relate surface seismic reflection data to subsurface geologic parameters. Developed initially for the case of homogeneous, isotropic, but arbitrarily dipping layered media, this technique has been extended to the inhomogeneous case. Now with the advent of new seismic techniques, such as vertical seismic profiling, three‐dimensional seismic methods, and shear‐wave techniques, the problem of velocity anisotropy is of growing concern to exploration seismologists. The essence of this paper is to extend the method of wavefront curvature to the “elliptically anisotropic” case in which the ray velocity varies elliptically with the direction of propagation. A fundamental feature of wave propagation in the anisotropic medium is that the direction of propagation of the disturbance (or the ray velocity direction) generally differs from that of the wavefront (or the phase velocity direction). Based on the assumption of two‐dimensional dipping layers in which velocity is “elliptically dependent” on the angle of propagation, relationships are developed between important seismic properties and model parameters. First, a relationship between the incident and refracted rays across the interface is expressed through the ray parameter. A geometrical divergence law is then developed relating the radii of the elliptical wave surfaces of the incident and refracted rays. A Dix‐type formula is finally derived which relates the normal moveout (NMO) velocity to the subsurface parameters. An example is shown to compare the radius of the wave surface and the NMO velocity for the elliptically anisotropic case with those for the equivalent isotropic case.

2021 ◽  
Author(s):  
Eric Roots ◽  
Graham Hill ◽  
Ben M. Frieman ◽  
James A. Craven ◽  
Richard S. Smith ◽  
...  

<p>The role of melts and magmatic/metamorphic fluids in mineralization processes is well established. However, the role of crustal architecture in defining source and sink zones in the middle to lower crust remains enigmatic. Integration of three dimensional magnetotelluric (MT) modelling and seismic reflection data across the Archean Abitibi greenstone belt of the Superior Province, Canada, reveals a ‘whole-of-crust’ mineralizing system and highlights the controls by crustal architecture on metallogenetic processes. Electrically conductive conduits in an otherwise resistive upper crust are coincident with truncations and offsets of seismic reflections that are mostly interpreted as major brittle-ductile fault zones. The spatial association between these features and low resistivity zones imaged in the 3D models suggest that these zones acted as pathways through which fluids and melts ascended toward the surface. At mid-crustal levels, these ‘conduit’ zones connect to ~50 km long, north-south striking conductors, and are inferred to represent graphite and/or sulphide deposited from cooling fluids. At upper mantle to lower crustal depths, east-west trending conductive zones dominate and display shallow dips. The upper mantle features are broadly coincident with the surface traces of the major deformation zones with which a large proportion of the gold endowment is associated. We suggest that these deep conductors represent interconnected graphitic zones perhaps augmented by sulphides that are relicts from metamorphic fluid and melt emplacement associated primarily with the later stages of regional deformation.  Thus, from the combined MT and seismic data, we develop a crustal-scale architectural model that is consistent with existing geological and deformational models, providing constraints on the sources for and signatures of fluid and magma emplacement that resulted in widespread metallogenesis in the Abitibi Subprovince.</p>


2021 ◽  
Author(s):  
Hongdan Deng ◽  
Ken McClay

<div>Basement fault reactivation, and the growth, interaction, and linkage with new fault segments are fundamentally three-dimensional and critical for understanding the evolution of fault network development in sedimentary basins. This paper analyses the evolution of a complex, basement-involved extensional fault network on the Enderby Terrace on the eastern margin of the Dampier sub-basin, NW Shelf of Australia. A high-resolution, depth-converted, 3D seismic reflection data volume is used to show that multiphase, oblique extensional reactivation of basement-involved faults controlled the development of the fault network in the overlying strata. Oblique reactivation of the pre-existing faults initially led to the formation of overlying, en échelon Late Triassic – Middle Jurassic fault segments that, as WNW–directed rifting progressed on the margin, linked by breaching of relay ramp to form two intersecting fault systems (F1 and F2-F4). Further reactivation in the Latest Jurassic – Early Cretaceous (NNW–SSE extension) produced an additional set of en échelon fault arrays in the cover strata. The final fault network consists of main or principal faults and subordinate or splay faults, together with branch lines that link the various components. Our study shows that breaching of relay ramps and/or vertical linkages produces vertical and horizontal branch lines giving complex final fault geometries. We find that repeated activity of the basement-involved faults tends to form continuous and planar fault architectures that favor displacement transfer between the main constituent segments along strike and with depth.</div>


1989 ◽  
Vol 26 (2) ◽  
pp. 305-321 ◽  
Author(s):  
François Marillier ◽  
Charlotte E. Keen ◽  
Glen S. Stockmal ◽  
Garry Quinlan ◽  
Harold Williams ◽  
...  

In 1986, 1181 km of marine seismic reflection data was collected to 18–20 s of two-way traveltime in the Gulf of St. Lawrence area. The seismic profiles sample all major surface tectono-stratigraphic zones of the Canadian Appalachians. They complement the 1984 deep reflection survey northeast of Newfoundland. Together, the seismic profiles reveal the regional three-dimensional geometry of the orogen.Three lower crustal blocks are distinguished on the seismic data. They are referred to as the Grenville, Central, and Avalon blocks, from west to east. The Grenville block is wedge shaped in section, and its subsurface edge follows the form of the Appalachian structural front. The Grenville block abuts the Central block at mid-crustal to mantle depths. The Avalon block meets the Central block at a steep junction that penetrates the entire crust.Consistent differences in the seismic character of the Moho help identify boundaries of the deep crustal blocks. The Moho signature varies from uniform over extended distances to irregular with abrupt depth changes. In places the Moho is offset by steep reflections that cut the lower crust and upper mantle. In other places, the change in Moho elevation is gradual, with lower crustal reflections following its form. In all three blocks the crust is generally highly reflective, with no distinction between a transparent upper crust and reflective lower crust.In general, Carboniferous and Mesozoic basins crossed by the seismic profiles overlie thinner crust. However, a deep Moho is found at some places beneath the Carboniferous Magdalen Basin.The Grenville block belongs to the Grenville Craton; the Humber Zone is thrust over its dipping southwestern edge. The Dunnage Zone is allochthonous above the opposing Grenville and Central blocks. The Gander Zone may be the surface expression of the Central block or may be allochthonous itself. There is a spatial analogy between the Avalon block and the Avalon Zone. Our profile across the Meguma Zone is too short to seismically distinguish this zone from the Avalon Zone.


2018 ◽  
Vol 477 (1) ◽  
pp. 223-240 ◽  
Author(s):  
D. M. Hodgson ◽  
H. L. Brooks ◽  
A. Ortiz-Karpf ◽  
Y. Spychala ◽  
D. R. Lee ◽  
...  

AbstractMany mass transport complexes (MTCs) contain up to kilometre-scale (mega)clasts encased in a debritic matrix. Although many megaclasts are sourced from the headwall areas, the irregular basal shear surfaces of many MTCs indicate that megaclast entrainment during the passage of flows into the deeper basin is also common. However, the mechanisms responsible for the entrainment of large blocks of substrate, and their influence on the longitudinal behaviour of the associated flows, have not been widely considered. We present examples of megaclasts from exhumed MTCs (the Neuquén Basin, Argentina and the Karoo Basin, South Africa) and MTCs imaged in three-dimensional seismic reflection data (Magdalena Fan, offshore Colombia and Santos Basin, offshore Brazil) to investigate these process–product interactions. We show that highly sheared basal surfaces are well developed in distal locations, sometimes extending beyond their associated deposit. This points to deformation and weakening of the substrate ahead of the flow, suggesting that preconditioning of the substrate by distributed shear ahead of, and to the side of, a mass flow could result in the entrainment of large fragments. An improved understanding of the interactions between flow evolution, seabed topography, and the entrainment and abrasion of megaclasts will help to refine estimates of run-out distances, and therefore the geohazard potential of submarine landslides.


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.


1999 ◽  
Vol 36 (2) ◽  
pp. 161-173 ◽  
Author(s):  
H John Broome ◽  
David Viljoen

The NATMAP Shield Margin Project, which began in 1991, straddles the Manitoba-Saskatchewan border and studied a rectangular area with east-west and north-south extents of approximately 250 and 150 km, respectively. Among the principal objectives of the NATMAP program were a compilation of a digital geoscience database for the study area and development of digital methodology to build this database and its utilization to aid in reaching the geological goals of the project. It was anticipated that these initiatives would encourage integrated interpretation of data, improve the effectiveness of the geological mapping process, and accelerate publication of results. One of the keys to effectively applying digital methods to a geological mapping project is maintaining data in digital form throughout the project, from the field through to publication. Transcription errors are eliminated and the production of both preliminary and final maps and digital products is accelerated. Access to digital data also facilitates application of digital tools for analysis and visualization. Integrated image products generated from geophysical and geological data sets were used throughout the project and assisted in development of geological models and their visualization. Three-dimensional visualization methods were used to combine NATMAP surface and subsurface mapping with Lithoprobe interpretations of seismic reflection data. The digital database was used to prepare a digital archive of the project which will be published in CD-ROM after completion of the project.


2005 ◽  
Vol 42 (4) ◽  
pp. 403-419 ◽  
Author(s):  
Z Hajnal ◽  
J Lewry ◽  
D White ◽  
K Ashton ◽  
R Clowes ◽  
...  

A three-dimensional model of the regional crustal architecture of the western Trans-Hudson Orogen, based on the interpretation of 590 km of deep-sounding seismic reflection data and a comparable length of existing seismic reflection information, is presented. The seismic images identify the regional geometry of the basal detachment zone (Pelican thrust) that separates juvenile allochthonous terranes from the underlying Archean microcontinent (Sask craton). The Sask Craton is inferred to have a minimum spatial extent of over 100 000 km2 with an associated crustal root that extends for 200 km along strike. During terminal collision, complete convergence of the Rae–Hearne and Superior continental blocks was precluded by the presence of the Sask Craton, resulting in the preservation of anomalous amounts of oceanic and associated sedimentary juvenile material. Along regional tectonic strike, consistency of crustal structure across the Rae–Hearne margin – Reindeer zone boundary is established. Several phases of tectonic development, including multistage subduction and continent–continent collision, are inferred for the western margin of the orogen. A bright, shallow (2–3.5 s two-way traveltime) band of reflectivity (Wollaston Lake reflector) imaged over ~150 000 km2 area is inferred to be a large post-orogenic mafic intrusion. A highly reflective, well-defined and structurally disturbed Moho discontinuity is mapped throughout the western Trans-Hudson Orogen. The present-day crustal architecture of the western Trans-Hudson Orogen is described in terms of the tectonic evolution within the region.


Author(s):  
Hongdan Deng ◽  
Ken McClay

Basement fault reactivation, and the growth, interaction, and linkage with new fault segments are fundamentally three-dimensional and critical for understanding the evolution of fault network development in sedimentary basins. This paper analyzes the evolution of a complex, basement-involved extensional fault network on the Enderby Terrace on the eastern margin of the Dampier sub-basin, North West Shelf of Australia. A high-resolution, depth-converted, 3-D seismic reflection data volume is used to show that multiphase, oblique extensional reactivation of basement-involved faults controlled the development of the fault network in the overlying strata. Reactivation of the pre-existing faults initially led to the formation of overlying, en échelon Late Triassic−Middle Jurassic fault segments that, as WNW-directed rifting progressed on the margin, linked by breaching of relay zones to form two intersecting fault systems (F1 and F2−F4). Further reactivation in the latest Jurassic−Early Cretaceous (NNW-SSE extension) produced an additional set of en échelon fault arrays in the cover strata. The final fault network consists of main or principal faults and subordinate or splay faults, together with branch lines that link the various components. Our study shows that breaching of relay ramps and/or vertical linkages produces vertical and horizontal branch lines giving complex final fault geometries. We find that repeated activity of the basement-involved faults tends to form continuous and planar fault architectures that favor displacement transfer between the main constituent segments along strike and with depth.


Geophysics ◽  
1984 ◽  
Vol 49 (12) ◽  
pp. 2117-2131 ◽  
Author(s):  
Guan Cheng ◽  
Shimon Coen

The relationship between Born inversion and wave‐equation migration of common‐midpoint (CMP) stacked seismic reflection data is analytically determined. The three‐dimensional (3-D) velocity distribution obtained by Born inversion is shown to be directly related to the 3-D reflectivity function obtained by wave‐equation migration for full bandwidth or band‐limited data. The relationship is obtained by the reformulations of migration and Born inversion methods as inverse source problems for the 3-D wave equation. The reformulation leads to a definition of the reflectivity function as the source function for the wave equation. It also leads to determination of the Born inversion results by applying the algorithm for wave‐equation migration to modified surface data. The modified data are simply related to the CMP stacked data. Alternatively, Born inversion results may be obtained directly from the migrated section. Results from synthetic and recorded data are presented and found to be consistent with the theoretical developments.


Geophysics ◽  
1989 ◽  
Vol 54 (8) ◽  
pp. 991-1000 ◽  
Author(s):  
R. Snieder ◽  
M. Y. Xie ◽  
A. Pica ◽  
A. Tarantola

Recorded seismic reflection waveforms contain information as to the small‐scale variations of impedance and the large‐scale variations of velocity. This information can be retrieved by minimizing the misfit between the recorded waveforms and synthetic seismograms as a function of the model parameters. Because of the different physical characters of the velocity and the impedance, we update these parameters in an alternating fashion, which amounts to a relaxation approach to the minimization of the waveform misfit. As far as the impedance is concerned, this minimization can be performed efficiently using gradient algorithms. For the inversion for seismic velocities, gradient methods do not work nearly as well; therefore, we use different minimization methods for determining impedances and velocities. However, the determination of the impedance and the determination of the velocity are strongly coupled; relaxation is most effective when this coupling is as weak as possible. Weak coupling can be achieved partially by parameterizing the impedances not as a function of depth but as a function of traveltime. A nonlinear, nonlocal method is presented for determining the smooth reference velocity from seismic reflection data. This technique is applied both to synthetic seismograms and to real marine data. In both cases, the velocity information implicitly contained in the curvature of the reflection hyperbolas was fully retrieved using nonlinear waveform optimization. In this way, it is possible to reconstruct both the impedance contrast and the smooth reference velocity from band‐limited seismic reflection data using a single waveform‐fit criterion.


Sign in / Sign up

Export Citation Format

Share Document