Crustal Structure of Eastern North Carolina: Piedmont and Coastal Plain

2019 ◽  
Vol 109 (6) ◽  
pp. 2288-2304 ◽  
Author(s):  
Shuai Zhao ◽  
Wenbin Guo
Keyword(s):  
North Carolina ◽  
Wave Velocity ◽  
Coastal Plain ◽  
Lower Crust ◽  
Thin Layers ◽  
P Wave ◽  
Upper Crust ◽  
S Wave ◽  
The West ◽  
Atlantic Margin

Abstract We present the results from an onshore seismic refraction and wide‐angle reflection profile, conducted in 2015, across the coastal plain and eastern Piedmont provinces of North Carolina. We use forward modeling to create 1D synthetic seismogram models and then invert first break picks to create 2D P‐ and S‐wave velocity models. The crustal thickness is 38 km beneath the Piedmont and central coastal plain, but it thins to 32 km at the coastline. The average thickness of the upper crust is 11 km with an average P‐wave velocity (VP) of 6.0  km/s and S‐wave velocity (VS) of 3.5  km/s. A prominent seismic low‐velocity zone (LVZ) (VP<6.0 and VS<3.6  km/s) exists between the depths of 6 and 11 km, beneath the western third of the seismic profile. The middle crust varies greatly in thickness, increasing from 3 km in the west (eastern Piedmont) to 13 km in the east (coastal plain), with seismic velocities of 6.5  km/s for VP and 3.8  km/s for VS. The lower crust thins significantly toward the rifted Atlantic margin, decreasing from 24 km thick in the west (Piedmont) to 8 km at the coastline, with velocities of approximately 6.9  km/s for VP and 3.9  km/s for VS. We estimate the composition of the crust by comparing the measured values of VP and Poisson’s ratio with laboratory measurements. The upper and middle crusts are in agreement with a felsic composition, while the lower crustal composition is predominately felsic to intermediate. The LVZ in the upper crust is associated with thin layers of the mylonitic rocks involved in the top and the bottom of thrusting, and the top of the lower crust could be the master detachment fault during the thin‐skinned Alleghanian orogeny. The eastward thinning of the lower crust is consistent with crustal extension during the Mesozoic rifting of the Atlantic margin.

Three-parameter prestack seismic inversion based on L1-2 minimization

Geophysics ◽  
2019 ◽  
Vol 84 (5) ◽  
pp. R753-R766 ◽  
Author(s):  
Lingqian Wang ◽  
Hui Zhou ◽  
Yufeng Wang ◽  
Bo Yu ◽  
Yuanpeng Zhang ◽  
...  
Keyword(s):  
Wave Velocity ◽  
Seismic Inversion ◽  
Thin Layers ◽  
P Wave ◽  
S Wave ◽  
Reservoir Prediction ◽  
Linear Inversion ◽  
Prestack Inversion ◽  
The Difference ◽  
Predictive Filtering

Prestack inversion has become a common approach in reservoir prediction. At present, the critical issue in the application of seismic inversion is the estimation of elastic parameters in the thin layers and weak reflectors. To improve the resolution and the accuracy of the inversion results, we introduced the difference of [Formula: see text] and [Formula: see text] norms as a nearly unbiased approximation of the sparsity of a vector, denoted as the [Formula: see text] norm, to the prestack inversion. The nonconvex penalty function of the [Formula: see text] norm can be decomposed into two convex subproblems via the difference of convex algorithm, and each subproblem can be solved efficiently by the alternating direction method of multipliers. Compared with the [Formula: see text] norm regularization, the [Formula: see text] minimization can reconstruct reflectivities more accurately. In addition, the [Formula: see text]-[Formula: see text] predictive filtering was introduced to guarantee the lateral continuity of the location and the amplitude of the reflectivity series. The generalized linear inversion and [Formula: see text]-[Formula: see text] predictive filtering are combined for stable elastic impedance inversion results, and three parameters of P-wave velocity, S-wave velocity, and density can be inverted with the Bayesian linearized amplitude variation with offset inversion. The inversion results of synthetic and real seismic data demonstrate that the proposed method can effectively improve the resolution and accuracy of the inversion results.

scholarly journals Crustal velocity structure of the Deccan Volcanic Province, Indian Peninsula, from observed surface wave dispersion

2014 ◽  
Vol 57 (4) ◽  
Author(s):  
Gaddale Suresh ◽  
Satbir S. Teotia ◽  
Sankar N. Bhattacharya
Keyword(s):  
Wave Velocity ◽  
Rayleigh Waves ◽  
Velocity Structure ◽  
Crustal Thickness ◽  
P Wave ◽  
Upper Crust ◽  
S Wave ◽  
Deccan Volcanic Province ◽  
Epicentral Zone ◽  
Volcanic Province

<p>Through inversion of fundamental mode group velocities of Love and Rayleigh waves, we study the crustal and subcrustal structure across the central Deccan Volcanic Province (DVP), which is one of the world’s largest terrestrial flood basalts. Our analysis is based on broadband seismograms recorded at seismological station Bhopal (BHPL) in the central India from earthquakes located near west coast of India, with an average epicentral distance about 768 km. The recording station and epicentral zone are situated respectively on the northern and southern edges of DVP with wave paths across central DVP. The period of group velocity data ranges from 5 to 60 s for Rayleigh waves and 5 to 45 s for Love waves. Using the genetic algorithm, the observed data have been inverted to obtain the crust and subcrustal velocity structure along the wavepaths. Using this procedure, a similar velocity structure was also obtained earlier for the northwestern DVP, which is in the west of the present study region. Comparison of results show that the crustal thickness decreases westward from central DVP (39.6 km) to northwestern DVP (37.8 km) along with the decrease of thickness of upper crust; while the thickness of lower crust remains nearly same. From east to west S-wave velocity in the upper crust decreases by 2 to 3 per cent, while P-wave velocity in the whole crust and subcrust decreases by 3 to 6 per cent. The P- and S-wave velocities are positively correlated with crustal thickness and negatively correlated with earth’s heat flow. It appears that the elevated crustal and subcrustal temperature in the western side is the main factor for low velocities on this side.</p>

Dynamics of the extension in the Fonualei Rift in the northern Lau Basin at 16 °S

2021 ◽  
Author(s):  
Anna Jegen ◽  
Anke Dannowski ◽  
Heidrun Kopp ◽  
Udo Barckhausen ◽  
Ingo Heyde ◽  
...  
Keyword(s):  
Velocity Gradient ◽  
Lower Crust ◽  
Pacific Plate ◽  
P Wave ◽  
Upper Crust ◽  
Lau Basin ◽  
Australian Plate ◽  
The Pacific ◽  
Refraction Seismic ◽  
Roll Back

&lt;p&gt;The Lau Basin is a young back-arc basin steadily forming at the Indo-Australian-Pacific plate boundary, where the Pacific plate is subducting underneath the Australian plate along the Tonga-Kermadec island arc. Roughly 25 Ma ago, roll-back of the Kermadec-Tonga subduction zone commenced, which lead to break up of the overriding plate and thus the formation of the western Lau Ridge and the eastern Tonga Ridge separated by the emerging Lau Basin.&lt;/p&gt;&lt;p&gt;As an analogue to the asymmetric roll back of the Pacific plate, the divergence rates decline southwards hence dictating an asymmetric, V-shaped basin opening. Further, the decentralisation of the extensional motion over 11 distinct spreading centres and zones of active rifting has led to the formation of a composite crust formed of a microplate mosaic. A simplified three plate model of the Lau Basin comprises the Tonga plate, the Australian plate and the Niuafo'ou microplate. The northeastern boundary of the Niuafo'ou microplate is given by two overlapping spreading centres (OLSC), the southern tip of the eastern axis of the Mangatolu Triple Junction (MTJ-S) and the northern tip of the Fonualei Rift spreading centre (FRSC) on the eastern side. Slow to ultraslow divergence rates were identified along the FRSC (8-32 mm/a) and slow divergence at the MTJ (27-32 mm/a), both decreasing southwards. However, the manner of divergence has not yet been identified. Additional regional geophysical data are necessary to overcome this gap of knowledge.&lt;/p&gt;&lt;p&gt;Research vessel RV Sonne (cruise SO267) set out to conduct seismic refraction and wide-angle reflection data along a 185 km long transect crossing the Lau Basin at ~16 &amp;#176;S from the Tonga arc in the east, the overlapping spreading centres, FRSC1 and MTJ-S2, and extending as far as a volcanic ridge in the west. The refraction seismic profile consisted of 30 ocean bottom seismometers. Additionally, 2D MCS reflection seismic data as well as magnetic and gravimetric data were acquired.&lt;/p&gt;&lt;p&gt;The results of our P-wave traveltime tomography show a crust that varies between 4.5-6 km in thickness. Underneath the OLSC the upper crust is 2-2.5 km thick and the lower crust 2-2.5 km thick. The velocity gradients of the upper and lower crust differ significantly from tomographic models of magmatically dominated oceanic ridges. Compared to such magmatically dominated ridges, our final P-wave velocity model displays a decreased velocity gradient in the upper crust and an increased velocity gradient in the lower crust more comparable to tectonically dominated rifts with a sparse magmatic budget.&lt;/p&gt;&lt;p&gt;The dominance of crustal stretching in the regional rifting process leads to a tectonical stretching, thus thinning of the crust under the OLSC and therefore increasing the lower crust&amp;#8217;s velocity gradient. Due to the limited magmatic budget of the area, neither the magnetic anomaly nor the gravity data indicate a magmatically dominated spreading centre. We conclude that extension in the Lau Basin at the OLSC at 16 &amp;#176;S is dominated by extensional processes with little magmatism, which is supported by the distribution of seismic events concentrated at the northern tip of the FRSC.&lt;/p&gt;

Reconnaissance seismic refraction-reflection surveys in northwestern New Mexico

1984 ◽  
Vol 74 (4) ◽  
pp. 1263-1274
Author(s):  
Lawrence H. Jaksha ◽  
David H. Evans
Keyword(s):  
New Mexico ◽  
Wave Velocity ◽  
Lower Crust ◽  
Seismic Waves ◽  
Seismic Refraction ◽  
Velocity Model ◽  
P Wave ◽  
Uppermost Mantle ◽  
P Wave Velocity ◽  
Crustal Thinning

Abstract A velocity model of the crust in northwestern New Mexico has been constructed from an interpretation of direct, refracted, and reflected seismic waves. The model suggests a sedimentary section about 3 km thick with an average P-wave velocity of 3.6 km/sec. The crystalline upper crust is 28 km thick and has a P-wave velocity of 6.1 km/sec. The lower crust below the Conrad discontinuity has an average P-wave velocity of about 7.0 km/sec and a thickness near 17 km. Some evidence suggests that velocity in both the upper and lower crust increases with depth. The P-wave velocity in the uppermost mantle is 7.95 ± 0.15 km/sec. The total crustal thickness near Farmington, New Mexico, is about 48 km (datum = 1.6 km above sea level), and there is evidence for crustal thinning to the southeast.

Numerical and Theoretical Investigation on the Origin of the Multiple Peaks of the H/V Ratio Curve

2021 ◽  
Author(s):  
Wanbo Xiao ◽  
Siqi Lu ◽  
Yanbin Wang
Keyword(s):  
Rayleigh Wave ◽  
Wave Velocity ◽  
Subsurface Layer ◽  
P Wave ◽  
Theoretical Formula ◽  
S Wave ◽  
Multiple Peaks ◽  
Ratio Curve ◽  
Multiple Layer ◽  
Wave Resonance

&lt;p&gt;Despite the popularity of the horizontal to vertical spectral ratio (HVSR) method in site effect studies, the origin of the H/V peaks has been controversial since this method was proposed. Many previous studies mainly focused on the explanation of the first or single peak of the H/V ratio, trying to distinguish between the two hypotheses &amp;#8212; the S-wave resonance and ellipticity of Rayleigh wave. However, it is common both in numerical simulations and practical experiments that the H/V ratio exhibits multiple peaks, which is essential to explore the origin of the H/V peaks.&lt;/p&gt;&lt;p&gt;The cause for the multiple H/V peaks has not been clearly figured out, and once was simply explained as the result of multi subsurface layers. Therefore, we adopted numerical method to simulate the ambient noise in various layered half-space models and calculated the H/V ratio curves for further comparisons. The peak frequencies of the H/V curves accord well with the theoretical frequencies of S-wave resonance in two-layer models, whose frequencies only depend on the S wave velocity and the thickness of the subsurface layer. The same is true for models with varying model parameters. Besides, the theoretical formula of the S-wave resonance in multiple-layer models is proposed and then supported by numerical investigations as in the cases of two-layer models. We also extended the S-wave resonance to P-wave resonance and found that its theoretical frequencies fit well with the V/H peaks, which could be an evidence to support the S-wave resonance theory from a new perspective. By contrast, there are obvious differences between the higher orders of the H/V ratio peaks and the higher orders of Rayleigh wave ellipticity curves both in two-layer and multiple-layer models. The Rayleigh wave ellipticity curves are found to be sensitive to the Poisson&amp;#8217;s ratio and the thickness of the subsurface layer, so the variation of the P wave velocity can affect the peak frequencies of the Rayleigh wave ellipticity curves while the H/V peaks show slight change. The Rayleigh wave ellipticity theory is thus proved to be inappropriate for the explanation of the multiple H/V peaks, while the possible effects of the Rayleigh wave on the fundamental H/V peak still cannot be excluded.&lt;/p&gt;&lt;p&gt;Based on the analyses above, we proposed a new evidence to support the claim that the peak frequencies of the H/V ratio curve, except the fundamental peaks, are caused by S-wave resonance. The relationship between the P-wave resonance and the V/H peaks may also find further application.&lt;/p&gt;

Using Seismic Data to Predict Shale Pore Pressure and Overpressure Sweet Spots: A Case Study from the Lower Silurian Longmaxi Formation in Sichuan Basin, China

2021 ◽  
Author(s):  
Sheng Chen ◽  
Qingcai Zeng ◽  
Xiujiao Wang ◽  
Qing Yang ◽  
Chunmeng Dai ◽  
...  
Keyword(s):  
Prediction Model ◽  
Wave Velocity ◽  
Shale Gas ◽  
Seismic Data ◽  
P Wave ◽  
Formation Pressure ◽  
S Wave ◽  
P Wave Velocity ◽  
New Formation ◽  
Sweet Spots

Abstract Practices of marine shale gas exploration and development in south China have proved that formation overpressure is the main controlling factor of shale gas enrichment and an indicator of good preservation condition. Accurate prediction of formation pressure before drilling is necessary for drilling safety and important for sweet spots predicting and horizontal wells deploying. However, the existing prediction methods of formation pore pressures all have defects, the prediction accuracy unsatisfactory for shale gas development. By means of rock mechanics analysis and related formulas, we derived a formula for calculating formation pore pressures. Through regional rock physical analysis, we determined and optimized the relevant parameters in the formula, and established a new formation pressure prediction model considering P-wave velocity, S-wave velocity and density. Based on regional exploration wells and 3D seismic data, we carried out pre-stack seismic inversion to obtain high-precision P-wave velocity, S-wave velocity and density data volumes. We utilized the new formation pressure prediction model to predict the pressure and the spatial distribution of overpressure sweet spots. Then, we applied the measured pressure data of three new wells to verify the predicted formation pressure by seismic data. The result shows that the new method has a higher accuracy. This method is qualified for safe drilling and prediction of overpressure sweet spots for shale gas development, so it is worthy of promotion.

Constraints on the composition of the crust and uppermost mantle in northwestern Canada: Vp/Vs variations along Lithoprobe's SNorCLE transect

2005 ◽  
Vol 42 (6) ◽  
pp. 1205-1222 ◽  
Author(s):  
Gabriela Fernández-Viejo ◽  
Ron M Clowes ◽  
J Kim Welford
Keyword(s):  
Wave Velocity ◽  
Upper Mantle ◽  
Laboratory Data ◽  
High Ratio ◽  
P Wave ◽  
S Wave ◽  
Uppermost Mantle ◽  
Crust And Upper Mantle ◽  
Slave Craton ◽  
Crustal Composition

Shear-wave seismic data recorded along four profiles during the SNoRE 97 (1997 Slave – Northern Cordillera Refraction Experiment) refraction – wide-angle reflection experiment in northwestern Canada are analyzed to provide S-wave velocity (Vs) models. These are combined with previous P-wave velocity (Vp) models to produce cross sections of the ratio Vp/Vs for the crust and upper mantle. The Vp/Vs values are related to rock types through comparisons with published laboratory data. The Slave craton has low Vp/Vs values of 1.68–1.72, indicating a predominantly silicic crustal composition. Higher values (1.78) for the Great Bear and eastern Hottah domains of the Wopmay orogen imply a more mafic than average crustal composition. In the western Hottah and Fort Simpson arc, values of Vp/Vs drop to ∼1.69. These low values continue westward for 700 km into the Foreland and Omineca belts of the Cordillera, providing support for the interpretation from coincident seismic reflection studies that much of the crust from east of the Cordilleran deformation front to the Stikinia terrane of the Intermontane Belt consists of quartzose metasedimentary rocks. Stikinia shows values of 1.78–1.73, consistent with its derivation as a volcanic arc terrane. Upper mantle velocity and ratio values beneath the Slave craton indicate an ultramafic peridotitic composition. In the Wopmay orogen, the presence of low Vp/Vs ratios beneath the Hottah – Fort Simpson transition indicates the presence of pyroxenite in the upper mantle. Across the northern Cordillera, low Vp values and a moderate-to-high ratio in the uppermost mantle are consistent with the region's high heat flow and the possible presence of partial melt.

Nonlinear two‐dimensional elastic inversion of multioffset seismic data

Geophysics ◽  
1987 ◽  
Vol 52 (9) ◽  
pp. 1211-1228 ◽  
Author(s):  
Peter Mora
Keyword(s):  
Wave Equation ◽  
Wave Velocity ◽  
A Priori ◽  
Gaussian Model ◽  
Wave Impedance ◽  
P Wave ◽  
S Wave ◽  
Elastic Wave Equation ◽  
Gradient Direction ◽  
P Wave Velocity

The treatment of multioffset seismic data as an acoustic wave field is becoming increasingly disturbing to many geophysicists who see a multitude of wave phenomena, such as amplitude‐offset variations and shearwave events, which can only be explained by using the more correct elastic wave equation. Not only are such phenomena ignored by acoustic theory, but they are also treated as undesirable noise when they should be used to provide extra information, such as S‐wave velocity, about the subsurface. The problems of using the conventional acoustic wave equation approach can be eliminated via an elastic approach. In this paper, equations have been derived to perform an inversion for P‐wave velocity, S‐wave velocity, and density as well as the P‐wave impedance, S‐wave impedance, and density. These are better resolved than the Lamé parameters. The inversion is based on nonlinear least squares and proceeds by iteratively updating the earth parameters until a good fit is achieved between the observed data and the modeled data corresponding to these earth parameters. The iterations are based on the preconditioned conjugate gradient algorithm. The fundamental requirement of such a least‐squares algorithm is the gradient direction which tells how to update the model parameters. The gradient direction can be derived directly from the wave equation and it may be computed by several wave propagations. Although in principle any scheme could be chosen to perform the wave propagations, the elastic finite‐ difference method is used because it directly simulates the elastic wave equation and can handle complex, and thus realistic, distributions of elastic parameters. This method of inversion is costly since it is similar to an iterative prestack shot‐profile migration. However, it has greater power than any migration since it solves for the P‐wave velocity, S‐wave velocity, and density and can handle very general situations including transmission problems. Three main weaknesses of this technique are that it requires fairly accurate a priori knowledge of the low‐ wavenumber velocity model, it assumes Gaussian model statistics, and it is very computer‐intensive. All these problems seem surmountable. The low‐wavenumber information can be obtained either by a prior tomographic step, by the conventional normal‐moveout method, by a priori knowledge and empirical relationships, or by adding an additional inversion step for low wavenumbers to each iteration. The Gaussian statistics can be altered by preconditioning the gradient direction, perhaps to make the solution blocky in appearance like well logs, or by using large model variances in the inversion to reduce the effect of the Gaussian model constraints. Moreover, with some improvements to the algorithm and more parallel computers, it is hoped the technique will soon become routinely feasible.

Shear‐wave velocity and density estimation from PS-wave AVO analysis: Application to an OBS dataset from the North Sea

Geophysics ◽  
2000 ◽  
Vol 65 (5) ◽  
pp. 1446-1454 ◽  
Author(s):  
Side Jin ◽  
G. Cambois ◽  
C. Vuillermoz
Keyword(s):  
Wave Velocity ◽  
North Sea ◽  
Random Noise ◽  
P Wave ◽  
S Wave ◽  
Data Set ◽  
Avo Analysis ◽  
The North ◽  
Avo Inversion ◽  
The North Sea

S-wave velocity and density information is crucial for hydrocarbon detection, because they help in the discrimination of pore filling fluids. Unfortunately, these two parameters cannot be accurately resolved from conventional P-wave marine data. Recent developments in ocean‐bottom seismic (OBS) technology make it possible to acquire high quality S-wave data in marine environments. The use of (S)-waves for amplitude variation with offset (AVO) analysis can give better estimates of S-wave velocity and density contrasts. Like P-wave AVO, S-wave AVO is sensitive to various types of noise. We investigate numerically and analytically the sensitivity of AVO inversion to random noise and errors in angles of incidence. Synthetic examples show that random noise and angle errors can strongly bias the parameter estimation. The use of singular value decomposition offers a simple stabilization scheme to solve for the elastic parameters. The AVO inversion is applied to an OBS data set from the North Sea. Special prestack processing techniques are required for the success of S-wave AVO inversion. The derived S-wave velocity and density contrasts help in detecting the fluid contacts and delineating the extent of the reservoir sand.

scholarly journals Discussion on the Theoretical Basis for Cross-Over Method Applied to Downhole Wave Velocity Test

2021 ◽  
Vol 2021 ◽  
pp. 1-9
Author(s):  
Qing Dong ◽  
Zheng-hua Zhou ◽  
Su Jie ◽  
Bing Hao ◽  
Yuan-dong Li
Keyword(s):  
Numerical Simulation ◽  
Numerical Simulations ◽  
Wave Velocity ◽  
Theoretical Basis ◽  
Influence Factors ◽  
P Wave ◽  
Engineering Practice ◽  
S Wave ◽  
Simulation Results ◽  
The Cross

At engineering practice, the theoretical basis for the cross-over method, used to obtain shear wave arrival time in the downhole method of the wave velocity test by surface forward and backward strike, is that the polarity of P-wave keeps the same, while the polarity of S-wave transforms when the direction of strike inverted. However, the characteristics of signals recorded in tests are often found to conflict with this theoretical basis for the cross-over method, namely, the polarity of the P-wave also transforms under the action of surface forward and backward strike. Therefore, 3D finite element numerical simulations were conducted to study the validity of the theoretical basis for the cross-over method. The results show that both shear and compression waves are observed to be in 180° phase difference between horizontal signal traces, consistent with the direction of excitation generated by reversed impulse. Furthermore, numerical simulation results prove to be reliable by the analytic solution; it shows that the theoretical basis for the cross-over method applied to the downhole wave velocity test is improper. In meanwhile, numerical simulations reveal the factors (inclining excitation, geophone deflection, inclination, and background noise) that may cause the polarity of the P-wave not to reverse under surface forward and backward strike. Then, as to reduce the influence factors, we propose a method for the downhole wave velocity test under surface strike, the time difference of arrival is based between source peak and response peak, and numerical simulation results show that the S-wave velocity by this method is close to the theoretical S-wave velocity of soil.

Sign in / Sign up

Export Citation Format

Share Document