Integrating S‐wave seismic‐reflection data and cone penetration test data using a multiangle multiscale approach

Geophysics ◽  
2004 ◽  
Vol 69 (2) ◽  
pp. 440-459 ◽  
Author(s):  
Ranajit Ghose ◽  
Jeroen Goudswaard

A cone penetration test (CPT) is the most common geotechnical testing method used to estimate in situ the strength properties of soil. Although CPT provides valuable information, this information is restricted to the location of the measurement. We propose a new concept to integrate shallow S‐wave reflection seismic data with CPT data in order to obtain laterally continuous subsoil information. In this vein, a valid quantitative means to relate seismic reflections to CPT data is a primary requirement. The approach proposed here is based on the characterization of the scaling behavior of the local fine‐scale S‐wave velocity information extracted from the seismic reflection data and the same behavior of the CPT cone resistance. The local velocity contrast information is extracted by linearized Zoeppritz inversion of the amplitude‐preserved prestack reflection data. We have formulated a multiscale analysis approach employing the continuous wavelet transform in order to quantitatively characterize the nature of change at an interface of the local S‐wave velocity contrast and the CPT cone resistance and to illuminate any relation between these two. The multiscale analysis estimates the singularity parameter α, which indicates the nature of the interfacial change. The application of our method to the field data has uncovered a striking relation between the nature of variation of the local S‐wave velocity contrast and that of CPT cone resistance; otherwise, such a relation was not visible. Detailed analyses of two extensive field datasets have shown that the lateral fine‐scale variation of soil strength, as seen by CPT cone resistance, has a close resemblance with the variation of the local S‐wave velocity function as seen by angle‐dependent reflection measurements. This leads to a unique possibility to integrate two very different in‐situ measurements—reflection seismic and CPT—providing laterally continuous detailed information of the soil layer boundaries.

Geophysics ◽  
1986 ◽  
Vol 51 (10) ◽  
pp. 1893-1903 ◽  
Author(s):  
Albert Tarantola

The problem of interpretation of seismic reflection data can be posed with sufficient generality using the concepts of inverse theory. In its roughest formulation, the inverse problem consists of obtaining the Earth model for which the predicted data best fit the observed data. If an adequate forward model is used, this best model will give the best images of the Earth’s interior. Three parameters are needed for describing a perfectly elastic, isotropic, Earth: the density ρ(x) and the Lamé parameters λ(x) and μ(x), or the density ρ(x) and the P-wave and S-wave velocities α(x) and β(x). The choice of parameters is not neutral, in the sense that although theoretically equivalent, if they are not adequately chosen the numerical algorithms in the inversion can be inefficient. In the long (spatial) wavelengths of the model, adequate parameters are the P-wave and S-wave velocities, while in the short (spatial) wavelengths, P-wave impedance, S-wave impedance, and density are adequate. The problem of inversion of waveforms is highly nonlinear for the long wavelengths of the velocities, while it is reasonably linear for the short wavelengths of the impedances and density. Furthermore, this parameterization defines a highly hierarchical problem: the long wavelengths of the P-wave velocity and short wavelengths of the P-wave impedance are much more important parameters than their counterparts for S-waves (in terms of interpreting observed amplitudes), and the latter are much more important than the density. This suggests solving the general inverse problem (which must involve all the parameters) by first optimizing for the P-wave velocity and impedance, then optimizing for the S-wave velocity and impedance, and finally optimizing for density. The first part of the problem of obtaining the long wavelengths of the P-wave velocity and the short wavelengths of the P-wave impedance is similar to the problem solved by present industrial practice (for accurate data interpretation through velocity analysis and “prestack migration”). In fact, the method proposed here produces (as a byproduct) a generalization to the elastic case of the equations of “prestack acoustic migration.” Once an adequate model of the long wavelengths of the P-wave velocity and of the short wavelengths of the P-wave impedance has been obtained, the data residuals should essentially contain information on S-waves (essentially P-S and S-P converted waves). Once the corresponding model of S-wave velocity (long wavelengths) and S-wave impedance (short wavelengths) has been obtained, and if the remaining residuals still contain information, an optimization for density should be performed (the short wavelengths of impedances do not give independent information on density and velocity independently). Because the problem is nonlinear, the whole process should be iterated to convergence; however, the information from each parameter should be independent enough for an interesting first solution.


2015 ◽  
Vol 55 (2) ◽  
pp. 400 ◽  
Author(s):  
Catherine Belgarde ◽  
Gianreto Manatschal ◽  
Nick Kusznir ◽  
Sonia Scarselli ◽  
Michal Ruder

Acquisition of long-offset (8–10 km), long-record length (12–18 sec), 2D reflection seismic and ship-borne potential fields data (WestraliaSpan by Ion/GXT and New Dawn by PGS) on the North West Shelf of Australia provide the opportunity to study rift processes in the context of modern models for rifted margins (Manatschal, 2004). Basement and Moho surfaces were interpreted on seismic reflection data. Refraction models from Geoscience Australia constrain Moho depth and initial densities for gravity modelling through standard velocity-density transformation. 2D joint inversion of seismic reflection and gravity data for Moho depth and basement density constrain depth to basement on seismic. 2D gravity and magnetic intensity forward modelling of key seismic lines constrain basement thickness, type and density. Late Permian and Jurassic-Early Cretaceous rift zones were mapped on seismic reflection data and constrained further by inversion and forward modelling of potential fields data. The Westralian Superbasin formed as a marginal basin in Eastern Gondwana during the Late Permian rifting of the Sibumasu terrane. Crustal necking was localised along mechanically-weak Proterozoic suture belts or Early Paleozoic sedimentary basins (such as Paterson and Canning). Mechanically-strong cratons (such as Pilbara and Kimberley) remained intact, resulting in necking and hyper-extension at their edges. Late Permian hyper-extended areas (such as Exmouth Plateau) behaved as mechanically-strong blocks during the Jurassic to Early Cretaceous continental break-up. Late Permian necking zones were reactivated as failed-rift basins and localised the deposition of the Jurassic oil-prone source rocks that have generated much of the oil discovered on the North West Shelf.


Solid Earth ◽  
2019 ◽  
Vol 10 (5) ◽  
pp. 1651-1662 ◽  
Author(s):  
Juan Alcalde ◽  
Clare E. Bond ◽  
Gareth Johnson ◽  
Armelle Kloppenburg ◽  
Oriol Ferrer ◽  
...  

Abstract. The use of conceptual models is essential in the interpretation of reflection seismic data. It allows interpreters to make geological sense of seismic data, which carries inherent uncertainty. However, conceptual models can create powerful anchors that prevent interpreters from reassessing and adapting their interpretations as part of the interpretation process, which can subsequently lead to flawed or erroneous outcomes. It is therefore critical to understand how conceptual models are generated and applied to reduce unwanted effects in interpretation results. Here we have tested how interpretation of vertically exaggerated seismic data influenced the creation and adoption of the conceptual models of 161 participants in a paper-based interpretation experiment. Participants were asked to interpret a series of faults and a horizon, offset by those faults, in a seismic section. The seismic section was randomly presented to the participants with different horizontal–vertical exaggeration (1:4 or 1:2). Statistical analysis of the results indicates that early anchoring to specific conceptual models had the most impact on interpretation outcome, with the degree of vertical exaggeration having a subdued influence. Three different conceptual models were adopted by participants, constrained by initial observations of the seismic data. Interpreted fault dip angles show no evidence of other constraints (e.g. from the application of accepted fault dip models). Our results provide evidence of biases in interpretation of uncertain geological and geophysical data, including the use of heuristics to form initial conceptual models and anchoring to these models, confirming the need for increased understanding and mitigation of these biases to improve interpretation outcomes.


Geophysics ◽  
2000 ◽  
Vol 65 (1) ◽  
pp. 286-292 ◽  
Author(s):  
Lingxiu Jiao ◽  
Wooil M. Moon

Seismic signals in deep crustal surveys are often contaminated with various types of noise, mainly caused by the low signal‐to‐noise (S/N) earth environment. A variance fractal dimension (VFD) technique is investigated and tested with real data sets for detection of seismic refraction signals from background noise. The data tested in this study were collected during the 1992 Lithoprobe Abitibi‐Grenville Transect high‐resolution refraction and wide‐angle reflection seismic experiments. The sharpness of transition features on the VFD trajectory is used as a criterion for distinguishing specific seismic phases. The window size and window interval applied in the application of VFD technique were determined using synthetic seismic data for generation of the optimum VFD trajectory. The window size of 48 samples and the window interval of 8 sample intervals were chosen to calculate the fractal dimension values and create the trajectories for detecting phases Pg, Pn, PmP, and ground roll. The VFD technique was also tested and applied for automatic detection of the first breaks in the high‐resolution seismic reflection data collected during the 1990 Lithoprobe regional and high‐resolution seismic surveys. The sharp transition features corresponding to the first arrivals in the seismic reflection data are distinct and provide us with a robust and powerful tool for separating the seismic signals from noise.


2019 ◽  
Vol 149 (3) ◽  
pp. 279 ◽  
Author(s):  
Gábor Csaba Tari ◽  
Viktória Németh ◽  
Ferenc Horváth † ◽  
Viktor Wesztergom

The so-called Transdanubian Conductivity Anomaly (TCA) of the Hungarian part of the NW Pannonian Basin has been well known for more than five decades. The exceptionally low resistivity (i.e. 1–2 Ωm) zone has a very large areal extent (on the order a few thousand km2) and it is an entirely subsurface anomaly occurring at depth between circa 3–15 km, with no corresponding outcrops. Various geological explanations of this enigmatic crustal-scale geophysical anomaly range from invoking sub-horizontal Alpine nappe contacts to sub-vertical dikes with graphite and/or saline fluid content. Only one possible analogue outcrop area was considered for the high conductivity anomaly so far, namely the Drauzug/Gailtal area of the Eastern Alps in Austria, some 300 km to the West from the TCA area. Previous attempts to find correspondence between the TCA and prominent seismic reflectors seen on 2D seismic reflection profiles were based on data acquired by research institutions. This study systematically correlates, for the first time, the TCA with 2D industry seismic reflection data in the same area. Our new results show a very strong correlation between the subsurface extent and location of the TCA with various sub-horizontally oriented Cretaceous Alpine nappe surfaces. In addition, we draw on the latest structural correlation of the Alpine nappe stack of the Transdanubian Range with its proper tectonic counterpart in the Eastern Alps.At the southern edge of the Upper Austroalpine units in northern Styria, in the Veitsch Nappe of the Greywacke Zone, numerous graphite localities are known historically. These laterally extensive graphite units in NW Styria formed as the result of greenschist-grade metamorphism of a Carboniferous coal sequence during the Cretaceous. For the first time, we describe here one well penetration of possibly age-equivalent graphitic units in NW Hungary. Correlation of the magnetotelluric anomaly with the distinct reflection seismic signature suggests that the same Palaeozoic graphitebearing Upper Austroalpine units should be present at 3–15 km depth in our study area.Therefore we propose that the best explanation for the observed extent and geometry of the TCA is the presence of graphite in subhorizontal, tectonically thinned detachment surfaces at the base of the Upper Austroalpine nappe edifice of NW Hungary


Geophysics ◽  
2009 ◽  
Vol 74 (6) ◽  
pp. G35-G45 ◽  
Author(s):  
Laura Valentina Socco ◽  
Daniele Boiero ◽  
Sebastiano Foti ◽  
Roger Wisén

Seismic reflection data contain surface waves that can be processed and interpreted to supply shear-wave velocity models along seismic reflection lines. The coverage of seismic reflection data allows the use of automated multifold processing to extract high-quality dispersion curves and experimental uncertainties in a moving spatial window. The dispersion curves are then inverted using a deterministic, laterally constrained inversion to obtain a pseudo-2D model of the shear-wave velocity. A Monte Carlo global search inversion algorithm optimizes the parameterization. When the strategy is used with synthetic and field data, consistent final models with smooth lateral variations are successfully retrieved. This method constitutes an improvement over the individual inversion of single dispersion curves.


Sign in / Sign up

Export Citation Format

Share Document