P-wave Traveltime Research Articles

Development of New Mathematical Model for Compressional and Shear Sonic Times from Wireline Log Data Using Artificial Intelligence Neural Networks (White Box)

Compressional (P-wave) and shear (S-wave) velocities are used to estimate the dynamic geomechanical properties including: Poisson’s ratio, Young’s modulus, and Lame parameters. These parameters are mainly used in estimating the static properties of the formation rocks as well as the in situ stresses. The sonic logs are not always available, epically for old wellbores. Also, in several occasions when the sonic logs are available, missing sections found in the well logs might affect the analysis results. To the authors’ knowledge, there is no single straightforward correlation that can be used to accurately estimate both P- and S-wave travel times directly from the well log data. Most of the existing correlations use the P-wave velocity to measure the S-wave velocity. The main purpose of this study is to develop accurate and simple empirical models using wireline log data (bulk density, gamma ray, and neutron porosity) to predict the sonic travel times (P-wave and S-wave). These wireline logs are slandered wireline log data that are commonly recorded in most of the wells. Three robust artificial intelligence techniques, namely: support vector machine (SVM), artificial neural network (ANN), and adaptive neurofuzzy interference systems (ANFIS), were employed and compared based on their prediction performance. Ultimately, using the weights and biases of optimized ANN model, a simple generalized empirical correlation is derived that can be used without the need of costly commercial software’s to run the AI models. The obtained results showed that ANN, ANFIS, and SVM can be used to estimate P-wave and S-wave travel times. ANN outperformed the ANFIS and SVM by yielding the lowest average absolute percentage error (AAPE) and the highest coefficient of determination $$\left( R^{2}\right) $$ for predicting P-wave and S-wave travel times. ANN model could predict the P-wave and S-wave travel times from wireline log data with high accuracy giving $$R^{2}$$ of 0.98 when compared to actual field data. In addition, the developed empirical correlations prediction completely matched the ANNs prediction. The AAPE of the predicted P and S-waves travel times was less than 5%. The developed correlations are very accurate and can help geomechanical engineers to determine the dynamic geomechanical properties (Poisson’s ratio and Young’s modulus) and propose any operation in case where sonic logs are missing.

Lateral displacement of crustal units relative to underlying mantle lithosphere: Example from the Bohemian Massif

We propose a mechanical model of deformation of the entire lithosphere of the Bohemian Massif (BM), whose core is formed by an asymmetric block of the Teplá-Barrandian (TB) unit in between the Saxothuringian (ST) and Moldanubian (MD) units. For the modelling, we have re-processed P-wave travel times recorded during the last two decades at dense networks of seismic stations installed in the BM during several passive seismic experiments. We also use previous results of anisotropic studies based on splitting of teleseismic shear waves. This allows us to refine estimates of the lithosphere thickness and delimit deep margins of the individual mantle lithosphere domains. The domains are rigid enough to preserve pre-orogenic olivine fabrics differently oriented in each of the units. Shapes and dips of the mantle boundaries, representing major zones of weakness inherited from the Variscan amalgamation of independent microplates, indicate that north-westward subductions beneath the TB unit dominated tectonic development of the core of the BM. Two mantle lithosphere domains with different fabric orientations, separated by a WSW-ENE striking shear zone, underlie the TB crust. The NW domain is the TB mantle lithosphere, while the SE domain is the MD mantle lithosphere thrust under the TB crust. Lithosphere of the north-western TB domain, compressed between early Variscan subductions of the ST continental lithosphere from the northwest and the MD continental lithosphere from the southeast, was pushed south-westward by about 50km. Though the crust of the south-westerly TB promontory is commonly attributed to the MD unit, apparently it preserves the TB mantle lithosphere. The shifted TB lithosphere provides compelling evidence in support of older views suggesting that the Zone Erbendorf-Vohenstrauss (ZEV) originally belonged to the tilted western rim of the TB unit. During the final phase of the assemblage of the BM, the rigid TB lithosphere was disrupted by the southward pushing ST lithosphere along the newly formed NW-SE striking Jáchymov Fault Zone (JFZ). This lithosphere-scale process most likely changed the tectonic regime, released subduction-related forces and started the gravity-dominated tectonics.

Shallow and deep lithosphere slabs beneath the Dinarides from teleseismic tomography as the result of the Adriatic lithosphere downwelling

The study area covers the Dinarides and southwestern part of the Pannonian basin as the marginal zone between the Adriatic microplate (African plate) and the Pannonian tectonic segment (Eurasian plate). We created a three-dimensional seismic velocity model to 450km depth using teleseismic tomography. Our travel-time dataset was collected by means of 40 seismic stations from the ORFEUS database and Croatian Seismological Survey database. A set of 90 teleseismic earthquakes were selected in the time range 2014–2015, and relative P-wave travel-time residuals were calculated. For the first time the seismic P-wave velocity model of a relatively high resolution on the entire Dinaridic mountain belt was obtained. Based on this model, a more reliable insight in the relations of the lithosphere plates has been achieved.We imaged a fast velocity anomaly extending underneath the entire Dinaridic mountain belt which indicates cold, rigid materials. The anomaly is steeply sloping towards the northeast and directly indicates the sinking of the Adriatic microplate underneath the Pannonian tectonic segment. In the Northern Dinarides the anomaly extends to the depth of 250km, whereas in the Southern Dinarides it covers greater depths, up to 450km. The shallow Adriatic slab extends along the External Dinarides, while the deep Adriatic slab extends beneath the Internal Dinarides and ophiolite zones in the area of central and southern Dinarides. Different slab depths are interpreted as the faster convergence of the plate in the southern Dinarides than in the northern, or the convergence of the plates had started in the southern part and systematically developed to the north.

P-wave anisotropy, mantle wedge flow and olivine fabrics beneath Japan

We present a new 3-D anisotropic P-wave velocity (Vp) model for the crust and upper mantle of the Japan subduction zone obtained by inverting a large number of high-quality P-wave traveltime data of local earthquakes and teleseismic events. By assuming orthorhombic anisotropy with a vertical symmetry axis existing in the modeling space, isotropic Vp tomography and 3-D Vp azimuthal and radial anisotropies are determined simultaneously. According to a simple flow field and the obtained Vp anisotropic tomography, we estimate the distribution of olivine fabrics in the mantle wedge. Our results show that the forearc mantle wedge above the subducting Pacific slab beneath NE Japan exhibits an azimuthal anisotropy with trench-parallel fast velocity directions (FVDs) and Vhf > Vv > Vhs (here Vv is Vp in the vertical direction, Vhf and Vhs are P-wave velocities in the fast and slow directions in the horizontal plane), where B-type olivine fabric with vertical trench-parallel flow may dominate. Such an anisotropic feature is not obvious in the forearc mantle wedge above the Philippine Sea (PHS) slab under SW Japan, probably due to higher temperatures and more fluids there associated with the young and warm PHS slab subduction. Trench-normal FVDs and Vhf > Vv > Vhs are generally revealed in the mantle wedge beneath the arc and backarc in Japan, where E-type olivine fabric with FVD-parallel horizontal flow may dominate. Beneath western Honshu, however, the mantle wedge exhibits an anisotropy of Vv > Vhf > Vhs and so C-type olivine fabric may dominate, suggesting that the water content is the highest there, because both the PHS and Pacific slabs exist there and their dehydration reactions release abundant fluids to the overlying mantle wedge.

An acoustic eikonal equation for attenuating orthorhombic media

Attenuating orthorhombic models are often used to describe the azimuthal variation of the seismic wave velocity and attenuation in finely layered hydrocarbon reservoirs with vertical fractures. In addition to the P-wave related medium parameters, S-wave parameters are also present in the complex eikonal equation needed to describe the P-wave complex-valued traveltime in an attenuating orthorhombic medium, which increases the complexity of using the P-wave traveltime to invert for the medium parameters in practice. We have used the acoustic assumption to derive an acoustic eikonal equation that approximately governs the complex-valued traveltime of P-waves in an attenuating orthorhombic medium. For a homogeneous attenuating orthorhombic media, we solve the eikonal equation using a combination of the perturbation method and Shanks transform. For a horizontal attenuating orthorhombic layer, the real and imaginary parts of the complex-valued reflection traveltime have nonhyperbolic behaviors in terms of the source-receiver offset. Similar to the roles of normal moveout (NMO) velocity and anellipticity, the attenuation NMO velocity and the attenuation anellipticity characterize the variation of the imaginary part of the complex-valued reflection traveltime around zero source-receiver offset.

Crustal structure and evolution of the Arctic Caledonides: Results from controlled-source seismology

The continuation of the Caledonides into the Barents Sea has long been a subject of discussion, and two major orientations of the Caledonian deformation fronts have been suggested: NNW-SSE striking and NE-SW striking. A regional NW-SE oriented ocean bottom seismic profile across the western Barents Sea was acquired in 2014. In this paper we map the crust and upper mantle structure along this profile in order to discriminate between different interpretations of Caledonian structural trends and orientation of rift basins in the western Barents Sea. Modeling of P-wave travel times has been done using a ray-tracing method, and combined with gravity modeling. The results show high P-wave velocities (4km/s) close to the seafloor, as well as localized sub-horizontal high velocity zones (6.0km/s and 6.9km/s) at shallow depths which are interpreted as magmatic sills. Refractions from the top of the crystalline basement together with reflections from the Moho give basement velocities from 6.0km/s at the top to 6.7km/s at the base of the crust. P-wave travel time modeling of the OBS profile indicate an eastwards increase in velocities from 6.4km/s to 6.7km/s at the base of the crystalline crust, and the western part of the profile is characterized by a higher seismic reflectivity than the eastern part. This change in seismic character is consistent with observations from vintage reflection seismic data and is interpreted as a Caledonian suture extending through the Barents Sea, separating Barentsia and Baltica. Local deepening of Moho (from 27km to 33km depth) creates “root structures” that can be linked to the Caledonian compressional deformation or a suture zone imprinted in the lower crust. Our model supports a separate NE-SW Caledonian trend extending into the central Barents Sea, branching off from the northerly trending Svalbard Caledonides, implying the existence of Barentsia as an independent microcontinent between Laurentia and Baltica.

Velocity variations and uncertainty from transdimensional P-wave tomography of North America

High-resolution models of seismic velocity variations constructed using body-wave tomography inform the study of the origin, fate and thermochemical state of mantle domains. In order to reliably relate these variations to material properties including temperature, composition and volatile content, we must accurately retrieve both the patterns and amplitudes of variations and quantify the uncertainty associated with the estimates of each. For these reasons, we image the mantle beneath North America with P-wave traveltimes from USArray using a novel method for 3-D probabilistic body-wave tomography. The method uses a Transdimensional Hierarchical Bayesian framework with a reversible-jump Markov Chain Monte Carlo algorithm in order to generate an ensemble of possible velocity models. We analyse this ensemble solution to obtain the posterior probability distribution of velocities, thereby yielding error bars and enabling rigorous hypothesis testing. Overall, we determine that the average uncertainty (1σ) of compressional wave velocity estimates beneath North America is ∼0.25 per cent dVP/VP, increasing with proximity to complex structure and decreasing with depth. The addition of USArray data reduces the uncertainty beneath the Eastern US by over 50 per cent in the upper mantle and 25–40 per cent below the transition zone and ∼30 per cent throughout the mantle beneath the Western US. In the absence of damping and smoothing, we recover amplitudes of variations 10–80 per cent higher than a standard inversion approach. Accounting for differences in data coverage, we infer that the length scale of heterogeneity is ∼50 per cent longer at shallow depths beneath the continental platform than beneath tectonically active regions. We illustrate the model trade-off analysis for the Cascadia slab and the New Madrid Seismic Zone, where we find that smearing due to the limitations of the illumination is relatively minor.

Hunting for the Tristan mantle plume – An upper mantle tomography around the volcanic island of Tristan da Cunha

The active volcanic island Tristan da Cunha, located at the southwestern and youngest end of the Walvis Ridge – Tristan/Gough hotspot track, is believed to be the surface expression of a huge thermal mantle anomaly. While several criteria for the diagnosis of a classical hotspot track are met, the Tristan region also shows some peculiarities. Consequently, it is vigorously debated if the active volcanism in this region is the expression of a deep mantle plume, or if it is caused by shallow plate tectonics and the interaction with the nearby Mid-Atlantic Ridge. Because of a lack of geophysical data in the study area, no model or assumption has been completely confirmed.We present the first amphibian P-wave finite-frequency travel time tomography of the Tristan da Cunha region, based on cross-correlated travel time residuals of teleseismic earthquakes recorded by 24 ocean-bottom seismometers. The data can be used to image a low velocity structure southwest of the island. The feature is cylindrical with a radius of ∼100 km down to a depth of 250 km. We relate this structure to the origin of Tristan da Cunha and name it the Tristan conduit. Below 250 km the low velocity structure ramifies into narrow veins, each with a radius of ∼50 km. Furthermore, we imaged a linkage between young seamounts southeast of Tristan da Cunha and the Tristan conduit.

An acoustic eikonal equation for attenuating transversely isotropic media with a vertical symmetry axis

Seismic-wave attenuation is an important component of describing wave propagation. Certain regions, such as gas clouds inside the earth, exert highly localized attenuation. In fact, the anisotropic nature of the earth induces anisotropic attenuation because the quasi P-wave dispersion effect should be profound along the symmetry direction. We have developed a 2D acoustic eikonal equation governing the complex-valued traveltime of quasi P-waves in attenuating, transversely isotropic media with a vertical-symmetry axis (VTI). This equation is derived under the assumption that the complex-valued traveltime of quasi P-waves in attenuating VTI media are independent of the S-wave velocity parameter [Formula: see text] in Thomsen’s notation and the S-wave attenuation coefficient [Formula: see text] in Zhu and Tsvankin’s notation. We combine perturbation theory and Shanks transform to develop practical approximations to the acoustic attenuating eikonal equation, capable of admitting an analytical description of the attenuation in homogeneous media. For a horizontal-attenuating VTI layer, we also derive the nonhyperbolic approximations for the real and imaginary parts of the complex-valued reflection traveltime. These equations reveal that (1) the quasi SV-wave velocity and the corresponding quasi SV-wave attenuation coefficient given as part of Thomsen-type notation barely affect the ray velocity and ray attenuation of quasi P-waves in attenuating VTI media; (2) combining the perturbation method and Shanks transform provides an accurate analytic eikonal solution for homogeneous attenuating VTI media; (3) for a horizontal attenuating VTI layer with weak attenuation, the real part of the complex-valued reflection traveltime may still be described by the existing nonhyperbolic approximations developed for nonattenuating VTI media, and the imaginary part of the complex-valued reflection traveltime still has the shape of nonhyperbolic curves. In addition, we have evaluated the possible extension of the proposed eikonal equation to realistic attenuating media, an alternative perturbation solution to the proposed eikonal equation, and the feasibility of applying the proposed nonhyperbolic equation for the imaginary part of the complex-valued traveltime to invert for interval attenuation parameters.

Regional travel-time residual studies and station correction from 1-D velocity models for some stations around Peninsular Malaysia and Singapore

We have investigated the average P-wave travel-time residuals for some stations around Southern Thailand, Peninsular Malaysia and Singapore at regional distances. Six years (January, 2010–December, 2015) record of events from central and northern Sumatra was obtained from the digital seismic archives of Integrated Research Institute for Seismology (IRIS). The criteria used for the data selection are designed to be above the magnitude of mb 4.5, depth less than 200km and an epicentral distance shorter than 1000km. Within this window a total number of 152 earthquakes were obtained. Furthermore, data were filtered based on the clarity of the seismic phases that are manually picked. A total of 1088 P-wave arrivals and 962 S-wave arrivals were hand-picked from 10 seismic stations around the Peninsula. Three stations IPM, KUM, and KOM from Peninsular Malaysia, four stations BTDF, NTU, BESC and KAPK from Singapore and three stations SURA, SRIT and SKLT located in the southern part of Thailand are used. Station NTU was chosen as the Ref. station because it recorded the large number of events. Travel-times were calculated using three 1-D models (Preliminary Ref. Earth Model PREM (Dziewonski and Anderson, 1981, IASP91, and Lienert et al., 1986) and an adopted two-point ray tracing algorithm. For the three models, we corroborate our calculated travel-times with the results from the use of TAUP travel-time calculation software. Relative to station NTU, our results show that the average P wave travel-time residual for PREM model ranges from −0.16 to 0.45s for BESC and IPM respectively. For IASP91 model, the average residual ranges from −0.25 to 0.24s for SRIT and SKLT respectively, and ranges from −0.22 to 0.30s for KAPK and IPM respectively for Lienert et al. (1986) model. Generally, most stations have slightly positive residuals relative to station NTU. These corrections reflect the difference between actual and estimated model velocities along ray paths to stations and can compensate for heterogeneous velocity structure near individual stations. The computed average travel-time residuals can reduce errors attributable to station correction in the inversion of hypocentral parameters around the Peninsula. Due to the heterogeneity occasioned by the numerous fault systems, a better 1-D velocity model for the Peninsula is desired for more reliable hypocentral inversion and other seismic investigations.

Application of perturbation theory to a P-wave eikonal equation in orthorhombic media

The P-wave eikonal equation for orthorhombic (ORT) anisotropic media is a highly nonlinear partial differential equation requiring the solution of a sixth-order polynomial to obtain traveltimes, resulting in complex and time-consuming numerical solutions. To alleviate this complexity, we approximate the solution of this equation by applying a multiparametric perturbation approach. We also investigated the sensitivity of traveltime surfaces in ORT media with respect to three anelliptic parameters. As a result, a simple and accurate P-wave traveltime approximation valid for ORT media was derived. Two different possible anelliptic parameterizations were compared. One of the parameterizations includes anelliptic parameters defined at zero offset: [Formula: see text], [Formula: see text], and [Formula: see text]. Another parameterization includes anelliptic parameters defined for all symmetry planes: [Formula: see text], [Formula: see text], and [Formula: see text]. The azimuthal behavior of sensitivity coefficients with different parameterizations was used to analyze the crosstalk between anelliptic parameters.

Descending lithosphere slab beneath the Northwest Dinarides from teleseismic tomography

The area of study covers the marginal zone between the Adriatic microplate (African plate) and the Pannonian segment (Eurasian plate). We present a tomography model for this area, with special emphasis on the northwest Dinarides. A dense distribution of temporary seismic stations in the area of the Northern Dinarides along with permanent seismic stations located in the area, allowed us to construct this P-wave tomographic model. We assembled our travel-time dataset based on 26 seismic stations were used to collect the dataset. Teleseismic events were recorded for a period of 18 months and a set of 76 distant earthquakes were used to calculate the P-wave travel-time residuals. We calculated relative rather than absolute arrival-time residuals in the inversion to obtain depths of 0–400km.We imaged a pronounced fast velocity anomaly below the NW Dinarides which directly indicates a lithosphere slab downgoing beneath the Dinarides. This fast anomaly extends towards the NW direction to at least 250km depth, and we interpreted it as a descending lithosphere slab. The thrusting of the Adriatic microplate may be brought about by sub-lithosphere rising movement beneath the Pannonian region, along with a push from African plate. In our interpretation, the Adriatic lower lithosphere has been detached from the crust, and steeply sinks beneath the Dinarides. A lithosphere model of the contact between the Adriatic microplate and Pannonian tectonic segment was constructed based on the tomographic velocity model and results of previous crustal studies.

Velocity model building and imaging in the presence of shallow gas

Multiple shallow layers of low saturation gas charged sands cause poor seismic imaging and hamper exploitation of many oil and gas discoveries in SE Asia. These gas ‘clouds’ are characterized by bright high amplitude reflectors of low P-wave velocity beneath which the seismic reflections have delayed P-wave travel times, reduced bandwidth, reduced reflected and transmitted amplitudes and associated phase distortion (Figure 1). All gas clouds are not equal, the appearance of gas clouds varies worldwide and within regions. A single gas-charged sand may sometimes have almost no effect on P-wave propagation; a shallow region of limited spatial extent may sometimes be undershot resulting in an interpretable image or multiple shallow layers can result in complete ‘wipeout’ of the final image. We believe that both anelastic losses and multiple scattering phenomena (indicated by rays on Figure 1) contribute to weak reflected energy beneath the gas-charged sands. A similar problem in the North Sea led to the development of multi-component seabed seismic solutions utilizing converted waves to penetrate beneath the gas (Berg et al., 1994). These solutions have been successfully applied to the largest of the Malaysian producing fields and are the PETRONAS preferred technical solution (e.g. Radzi et al., 2015). As with many problems, a mixture of acquisition and processing related solutions produces the optimal result but sometimes processing alone can produce a cost-effective result. Despite cost reductions and technical advances, converted wave imaging is still expensive, difficult and time consuming. In this paper we attempt to build an efficient and accurate workflow for velocity model building of P-wave data in order to best exploit the latest imaging algorithms. We use a mixture of synthetic and real data examples from offshore Malaysia to illustrate this difficult and unresolved problem.

Mapping b-values beneath Abu Dabbab from June to August 2004 earthquake swarm

Abu Dabbab area is considered as one of the most active earthquake sources in Egypt. It is defined by its swarm type activity, and complicated stress pattern. This study was conducted to evaluate the two and three dimensional spatial distribution of b-value at Abu Dabbab area (Margin of the northern Red Sea Rift, Egypt). The gridding technique of Wiemer and Wyss (1997) was used to compute b-value using ZMAP software. The b-value is calculated from a catalog consisting of 850 well-located earthquakes, which were recorded from 1st June to August 2004, using the maximum likelihood method. These earthquakes were recorded by temporary digital seismic network, with magnitudes ranging from −1 to 3.4ML. It is important to mention that the variations of b-value with time cannot be easily detected for a short period. Hence, this study has been carried out to examine the variations of b-value in space. The computed b-value in the Abu Dabbab area does not follow a uniform distribution. A small volume of anomalously high b-value (b>1.8) exists in the central part of the area at a depth between 6 and 9km. This seems to agree with the reported low velocity value derived from previous P-wave travel time tomography studies (Hosny et al., 2009) and the low Q value (Abdel-Fattah et al., 2008). The existence of an anomalously high b-value region may be attributed to the presence of a magma reservoir or dyke zone beneath the northern Red Sea Rift that causes an intensively heterogeneous fractured crust or unusually high pore pressure.

Downhole microseismic monitoring for low signal-to-noise ratio events

Microseismic monitoring plays an important role in the process of hydraulic fracturing for shale gas/oil production. The accuracy of event location is an essential issue in microseismic monitoring. The data obtained from downhole monitoring system usually show a higher signal-to-noise ratio (SNR) than the recorded data from the surface. For small microseismic events, however, P waves recorded in a downhole array may be very weak, while S waves are generally dominant and strong. Numerical experiments suggest that inverting S-wave arrival times alone is not sufficient to constrain event locations. In this study, we perform extensive location tests with various noise effects using a grid search method that matches the travel time data of the S wave across a recording array. We conclude that fitting S-wave travel time data along with at least one P-wave travel time of the same event can significantly improve location accuracy. In practice, picking S-wave arrival time data and at least one P-wave pick is possible for many small events. We demonstrate that fitting the combination of the travel time data is a robust approach, which can help increase the number of microseismic events to be located accurately during hydraulic fracturing.

Density and P‐wave velocity structure beneath the Paraná Magmatic Province: Refertilization of an ancient lithospheric mantle

AbstractWe estimate density and P‐wave velocity perturbations in the mantle beneath the southeastern South America plate from geoid anomalies and P‐wave traveltime residuals to constrain the structure of the lithosphere underneath the Paraná Magmatic Province (PMP) and conterminous geological provinces. Our analysis shows a consistent correlation between density and velocity anomalies. The P‐wave speed and density are 1% and 15 kg/m3 lower, respectively, in the upper mantle under the Late Cretaceous to Cenozoic alkaline provinces, except beneath the Goiás Alkaline Province (GAP), where density (+20 kg/m3) and velocity (+0.5%) are relatively high. Underneath the PMP, the density is higher by about 50 kg/m3 in the north and 25 kg/m3 in the south, to a depth of 250 − 300 km. These values correlate with high‐velocity perturbations of +0.5% and +0.3%, respectively. Profiles of density perturbation versus depth in the upper mantle are different for the PMP and the adjacent Archean São Francisco (SFC) and Amazonian (AC) cratons. The Paleoproterozoic PMP basement has a high‐density root. The density is relatively low in the SFC and AC lithospheres. A reduction of density is a typical characteristic of chemically depleted Archean cratons. A more fertile Proterozoic and Phanerozoic subcontinental lithospheric mantle has a higher density, as deduced from density estimates of mantle xenoliths of different ages and composition. In conjunction with Re‐Os isotopic studies of the PMP basalts, chemical and isotopic analyses of peridodite xenoliths from the GAP in the northern PMP, and electromagnetic induction experiments of the PMP lithosphere, our density and P‐wave speed models suggest that the densification of the PMP lithosphere and flood basalt generation are related to mantle refertilization. Metasomatic refertilization resulted from the introduction of asthenospheric components from the mantle wedge above Proterozoic subduction zones, which surrounded the Paraná lithosphere. The high‐density PMP lithosphere is presently gravitationally unstable and prone to delamination.

Time-lapse monitoring of fluid-induced geophysical property variations within an unstable earthwork using P-wave refraction

A significant portion of the UK’s transportation system relies on a network of geotechnical earthworks (cuttings and embankments) that were constructed more than 100 years ago, whose stability is affected by the change in precipitation patterns experienced over the past few decades. The vulnerability of these structures requires a reliable, cost- and time-effective monitoring of their geomechanical condition. We have assessed the potential application of P-wave refraction for tracking the seasonal variations of seismic properties within an aged clay-filled railway embankment, located in southwest England. Seismic data were acquired repeatedly along the crest of the earthwork at regular time intervals, for a total period of 16 months. P-wave first-break times were picked from all available recorded traces, to obtain a set of hodocrones referenced to the same spatial locations, for various dates along the surveyed period of time. Traveltimes extracted from each acquisition were then compared to track the pattern of their temporal variability. The relevance of such variations over time was compared with the data experimental uncertainty. The multiple set of hodocrones was subsequently inverted using a tomographic approach, to retrieve a time-lapse model of [Formula: see text] for the embankment structure. To directly compare the reconstructed [Formula: see text] sections, identical initial models and spatial regularization were used for the inversion of all available data sets. A consistent temporal trend for P-wave traveltimes, and consequently for the reconstructed [Formula: see text] models, was identified. This pattern could be related to the seasonal distribution of precipitation and soil-water content measured on site.

Quantitative evaluation of rock brittleness and fracability based on elastic-wave velocity variation around borehole

Abstract Brittleness and fracability are two important rock properties in hydraulic fracturing of unconventional reservoirs. Based on the variation of compressional and shear velocity around borehole using acoustic measurement, an effective technique is developed to estimate these parameters to guide reservoir-fracturing. During drilling, when the rock is broken, a significant amount of drilling induced cracks will occur in the formation around the borehole, reducing elastic wave velocity around borehole and causing the wave velocity variation away from borehole. The radial variation of compressional and shear velocities of formation rocks surrounding a borehole were respectively obtained from P-wave travel time tomography and dipole shear-wave dispersion inversion. By integrating the two variation profiles along the radial direction, the brittleness-fracability index is obtained to estimate the brittleness and fracability of formation rocks. The index shows fairly good consistency and correlation with rock brittleness and fracability, which demonstrates the practicability and effectiveness of the proposed approach. Well log data analysis examples are presented to demonstrate the effectiveness of our technique.

Subcrustal structure of the black sea basin from seismological data

The P-wave travel time data from the earthquakes offshore and onshore around the Black Sea are used for the tomographic reconstruction of the three-dimensional (3D) velocity distribution in the lithosphere of the region. The preliminary refinement of the foci parameters (the coordinates and origin time) has reduced the random errors in the travel-time data. The earthquake data were supplemented by the previous deep seismic sounding (DSS) data on the profiles in Crimea and offshore off the Black Sea. The dataset included more than 4000 travel times overall. In order to eliminate the crustal effect, the travel times were reduced to a surface at a depth of 35 km corresponding to the mean Moho depth in the region. The improved crustal model was used for removing the contribution of the crust from the initial data. The new tomography method, which was recently developed by one of the authors and which relies on the assumption of smoothness of the lateral velocity variations, was applied for reconstructing the velocity structure of the upper mantle beneath the Black Sea up to a depth of 95 km. The lateral velocity variation maps at different depths and the vertical velocity distributions along the meridional and sublatitudinal cross sections across the Black Sea were constructed. High velocities were revealed in the subcrustal lithosphere, and the structural difference below two subbasins—the West Black Sea (WBS) and the East Black Sea (EBS) ones—was established. It shows that the high-velocity body below the WBS is located deeper than below the EBS and is distinguished by higher velocities. Based on these results, it is concluded that the lithosphere beneath the Black Sea has a continental origin.

C13 and Thomsen anisotropic parameter distributions for hydraulic fracture monitoring

Most shale reservoirs are described assuming transverse isotropy. A full characterization of such reservoirs requires the knowledge of six elastic parameters comprising, for instance, the density; P- and S-wave vertical velocities; and the three anisotropic Thomsen parameters [Formula: see text], [Formula: see text], and [Formula: see text]. For hydraulic fracturing monitoring purposes, the parameters [Formula: see text] and [Formula: see text] are typically estimated from perforation microseismic data by analyzing the variation with phase angles of the traveltimes of direct P-waves and direct SH-waves. The estimation of the Thomsen parameter [Formula: see text] is more challenging because most perforation microseismic data do not exhibit direct SV-waves to which the parameter [Formula: see text] is the most sensitive. A Bayesian inference method is proposed to estimate [Formula: see text] from traveltimes of secondary SV-waves observed on microseismic data. Because the interpretation of the physics underlying the observation of these secondary SV-waves is often not readily available for wellsite operations, a marginalization scheme is used to account for uncertainties associated with not knowing, or only partially knowing, the locations where they originated from. I determined the benefits of the proposed method by applying it to data sets collected in the Jonah field using a single monitoring well. The results revealed a significant variance reduction of more than 15-fold on the posterior probability distribution associated with [Formula: see text] when using the secondary-source SV-waves. Moreover, a more accurate estimate of the microseismically derived fractured volume is achieved with variance reduced by more than a factor of 75.