Green's function of the Lord–Shulman thermo-poroelasticity theory

2020 ◽  
Vol 221 (3) ◽  
pp. 1765-1776 ◽  
Author(s):  
Jia Wei ◽  
Li-Yun Fu ◽  
Zhi-Wei Wang ◽  
Jing Ba ◽  
José M Carcione
Keyword(s):  
Plane Wave ◽  
Green’S Function ◽  
Green's Function ◽  
Heterogeneous Media ◽  
Numerical Algorithms ◽  
Heat Diffusion ◽  
P Wave ◽  
Wave Analysis ◽  
S Wave ◽  
S Waves

SUMMARY The Lord–Shulman thermoelasticity theory combined with Biot equations of poroelasticity, describes wave dissipation due to fluid and heat flow. This theory avoids an unphysical behaviour of the thermoelastic waves present in the classical theory based on a parabolic heat equation, that is infinite velocity. A plane-wave analysis predicts four propagation modes: the classical P and S waves and two slow waves, namely, the Biot and thermal modes. We obtain the frequency-domain Green's function in homogeneous media as the displacements-temperature solution of the thermo-poroelasticity equations. The numerical examples validate the presence of the wave modes predicted by the plane-wave analysis. The S wave is not affected by heat diffusion, whereas the P wave shows an anelastic behaviour, and the slow modes present a diffusive behaviour depending on the viscosity, frequency and thermoelasticity properties. In heterogeneous media, the P wave undergoes mesoscopic attenuation through energy conversion to the slow modes. The Green's function is useful to study the physics in thermoelastic media and test numerical algorithms.

On the Green function of the Lord–Shulman thermoelasticity equations

2019 ◽  
Vol 220 (1) ◽  
pp. 393-403 ◽  
Author(s):  
Zhi-Wei Wang ◽  
Li-Yun Fu ◽  
Jia Wei ◽  
Wanting Hou ◽  
Jing Ba ◽  
...  
Keyword(s):  
Green’S Function ◽  
Green's Function ◽  
Parabolic Type ◽  
Second Order ◽  
P Wave ◽  
Thermal Mode ◽  
S Wave ◽  
Order Tensor ◽  
P Waves ◽  
Thermal Conductivities

SUMMARY Thermoelasticity extends the classical elastic theory by coupling the fields of particle displacement and temperature. The classical theory of thermoelasticity, based on a parabolic-type heat-conduction equation, is characteristic of an unphysical behaviour of thermoelastic waves with discontinuities and infinite velocities as a function of frequency. A better physical system of equations incorporates a relaxation term into the heat equation; the equations predict three propagation modes, namely, a fast P wave (E wave), a slow thermal P wave (T wave), and a shear wave (S wave). We formulate a second-order tensor Green's function based on the Fourier transform of the thermodynamic equations. It is the displacement–temperature solution to a point (elastic or heat) source. The snapshots, obtained with the derived second-order tensor Green's function, show that the elastic and thermal P modes are dispersive and lossy, which is confirmed by a plane-wave analysis. These modes have similar characteristics of the fast and slow P waves of poroelasticity. Particularly, the thermal mode is diffusive at low thermal conductivities and becomes wave-like for high thermal conductivities.

A study of earthquake mechanism using S wave data

1958 ◽  
Vol 48 (3) ◽  
pp. 201-219
Author(s):  
Wm. Mansfield Adams
Keyword(s):  
Wave Motion ◽  
P Wave ◽  
Wave Method ◽  
Wave Analysis ◽  
Poor Agreement ◽  
S Wave ◽  
Tectonic Earthquake ◽  
S Waves ◽  
The World ◽  
Earthquake Mechanism

Abstract The purpose of this paper is to determine from the seismograms of a tectonic earthquake the line of the motion which generated the observed S waves (tectonically, the A axis). By noting certain geometrical relationships between the faulting motion and the emitted S waves, it is possible to derive a method which determines the line of the generating motion from observations of the generated S waves. The results of the application of the proposed method of S wave analysis should, theoretically, make it possible to determine which of the two solutions given by the P wave method of analyzing the tectonic mechanism of earthquakes is the correct solution. The proposed procedure is applied to data collected from the original seismograms of four earthquakes as recorded at seismic observatories throughout the world. There is such poor agreement between the S wave results and the previous P wave solutions that it is necessary to conclude that one or more of the following is true: either the mechanism assumed is not the type actually occurring; the phase identified as the S wave does not correspond to the first P wave motion; the P wave method is incorrect or inadequate; or the S wave method is incorrect or inadequate. To select among the various possibilities necessitates a discussion of the relative merits, defects, and potentialities of the two methods.

S waves: Alaska and other earthquakes

1960 ◽  
Vol 50 (4) ◽  
pp. 581-597 ◽  
Author(s):  
William Stauder
Keyword(s):  
Focal Mechanism ◽  
Fault Plane ◽  
P Wave ◽  
Wave Analysis ◽  
Point Model ◽  
S Wave ◽  
S Waves ◽  
Wave Solutions ◽  
Single Force

ABSTRACT Techniques of S wave analysis are used to investigate the focal mechanism of four earthquakes. In all cases the results of the S wave analysis agree with previously determined P wave solutions and conform to a dipole with moment or single couple as the point model of the focus. Further, the data from S waves select one of the two nodal planes of P as the fault plane. Small errors in the determination of the angle of polarization of S are shown to result in scatter in the data of a peculiar character which might lead to misinterpretation. The same methods of analysis which in the present instances show excellent agreement with a dipole with moment source are the methods which in a previous paper required a single force type mechanism for a different group of earthquakes.

Separation and enhancement of split S‐waves on multicomponent shot records from the BIRPS WISPA experiment

Geophysics ◽  
1994 ◽  
Vol 59 (1) ◽  
pp. 131-139 ◽  
Author(s):  
M. Boulfoul ◽  
D. R. Watts
Keyword(s):  
Synthetic Data ◽  
P Wave ◽  
Moving Window ◽  
Wave Analysis ◽  
S Wave ◽  
Near Surface ◽  
S Waves ◽  
Wave Velocities ◽  
Wave Splitting ◽  
Explosive Source

Instantaneous rotations are combined with f-k filtering to extract coherent S‐wave events from multicomponent shot records recorded by British Institutions Reflection Profiling Syndicate (BIRPS) Weardale Integrated S‐wave and P‐wave analysis (WISPA) experiment. This experiment was an attempt to measure the Poisson’s ratio of the lower crest by measuring P‐wave and S‐wave velocities. The multihole explosive source technique did generate S‐waves although not of opposite polarization. Attempts to produce stacks of the S‐wave data are unsuccessful because S‐wave splitting in the near surface produced random polarizations from receiver group to receiver group. The delay between the split wavelets varies but is commonly between 20 to 40 ms for 10 Hz wavelets. Dix hyperbola are produced on shot records after instantaneous rotations are followed by f-k filtering. To extract the instantaneous polarization, the traces are shifted back by the length of a moving window over which the calculation is performed. The instantaneous polarization direction is computed from the shifted data using the maximum eigenvector of the covariance matrix over the computation window. Split S‐waves are separated by the instantaneous rotation of the unshifted traces to the directions of the maximum eigenvectors determined for each position of the moving window. F-K filtering is required because of the presence of mode converted S‐waves and S‐waves produced by the explosive source near the time of detonation. Examples from synthetic data show that the method of instantaneous rotations will completely separate split S‐waves if the length of the moving window over which the calculation is performed is the length of the combined split wavelets. Separation may be achieved on synthetic data for wavelet delays as small as two sample intervals.

Theory and modeling of constant-Q P- and S-waves using fractional time derivatives

Geophysics ◽  
2009 ◽  
Vol 74 (1) ◽  
pp. T1-T11 ◽  
Author(s):  
José M. Carcione
Keyword(s):  
Heterogeneous Media ◽  
Fourier Method ◽  
P Wave ◽  
Difference Operators ◽  
S Wave ◽  
Central Difference ◽  
S Waves ◽  
The Time Domain ◽  
Derivatives Of ◽  
Modeling Algorithm

I have developed and solved the constant-[Formula: see text] model for the attenuation of P- and S-waves in the time domain using a new modeling algorithm based on fractional derivatives. The model requires time derivatives of order [Formula: see text] applied to the strain components, where [Formula: see text] and [Formula: see text], with [Formula: see text] the P-wave or S-wave quality factor. The derivatives are computed with the Grünwald-Letnikov and central-difference fractional approximations, which are extensions of the standard finite-difference operators for derivatives of integer order. The modeling uses the Fourier method to compute the spatial derivatives, and therefore can handle complex geometries and general material-property variability. I verified the results by comparison with the 2D analytical solution obtained for wave propagation in homogeneous Pierre Shale. Moreover, the modeling algorithm was used to compute synthetic seismograms in heterogeneous media corresponding to a crosswell seismic experiment.

Fast plane-wave reverse time migration

Geophysics ◽  
2018 ◽  
Vol 83 (6) ◽  
pp. S549-S556 ◽  
Author(s):  
Xiongwen Wang ◽  
Xu Ji ◽  
Hongwei Liu ◽  
Yi Luo
Keyword(s):  
Time Delay ◽  
Plane Wave ◽  
Green’S Function ◽  
Green's Function ◽  
Time Domain ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Wave Source ◽  
Computation Cost ◽  
Time Migration

Plane-wave reverse time migration (RTM) could potentially provide quick subsurface images by migrating fewer plane-wave gathers than shot gathers. However, the time delay between the first and the last excitation sources in the plane-wave source largely increases the computation cost and decreases the practical value of this method. Although the time delay problem is easily overcome by periodical phase shifting in the frequency domain for one-way wave-equation migration, it remains a challenge for time-domain RTM. We have developed a novel method, referred as to fast plane-wave RTM (FP-RTM), to eliminate unnecessary computation burden and significantly reduce the computational cost. In the proposed FP-RTM, we assume that the Green’s function has finite-length support; thus, the plane-wave source function and its responding data can be wrapped periodically in the time domain. The wrapping length is the assumed total duration length of Green’s function. We also determine that only two period plane-wave source and data after the wrapping process are required for generating the outcome with adequate accuracy. Although the computation time for one plane-wave gather is twice as long as a normal shot gather migration, a large amount of computation cost is saved because the total number of plane-wave gathers to be migrated is usually much less than the total number of shot gathers. Our FP-RTM can be used to rapidly generate RTM images and plane-wave domain common-image gathers for velocity model building. The synthetic and field data examples are evaluated to validate the efficiency and accuracy of our method.

Numerical simulation of azimuthal acoustic logging in a borehole penetrating a rock formation boundary

Geophysics ◽  
2016 ◽  
Vol 81 (3) ◽  
pp. D283-D291 ◽  
Author(s):  
Peng Liu ◽  
Wenxiao Qiao ◽  
Xiaohua Che ◽  
Xiaodong Ju ◽  
Junqiang Lu ◽  
...  
Keyword(s):  
Domain Model ◽  
P Wave ◽  
S Wave ◽  
Angle Variation ◽  
P Waves ◽  
Acoustic Logging ◽  
S Waves ◽  
Rock Formation ◽  
Logging Tool ◽  
Difference Time

We have developed a new 3D acoustic logging tool (3DAC). To examine the azimuthal resolution of 3DAC, we have evaluated a 3D finite-difference time-domain model to simulate a case in which the borehole penetrated a rock formation boundary when the tool worked at the azimuthal-transmitting-azimuthal-receiving mode. The results indicated that there were two types of P-waves with different slowness in waveforms: the P-wave of the harder rock (P1) and the P-wave of the softer rock (P2). The P1-wave can be observed in each azimuthal receiver, but the P2-wave appears only in the azimuthal receivers toward the softer rock. When these two types of rock are both fast formations, two types of S-waves also exist, and they have better azimuthal sensitivity compared with P-waves. The S-wave of the harder rock (S1) appears only in receivers toward the harder rock, and the S-wave of the softer rock (S2) appears only in receivers toward the softer rock. A model was simulated in which the boundary between shale and sand penetrated the borehole but not the borehole axis. The P-wave of shale and the S-wave of sand are azimuthally sensitive to the azimuth angle variation of two formations. In addition, waveforms obtained from 3DAC working at the monopole-transmitting-azimuthal-receiving mode indicate that the corresponding P-waves and S-waves are azimuthally sensitive, too. Finally, we have developed a field example of 3DAC to support our simulation results: The azimuthal variation of the P-wave slowness was observed and can thus be used to reflect the azimuthal heterogeneity of formations.

Scalar reverse‐time depth migration of prestack elastic seismic data

Geophysics ◽  
2001 ◽  
Vol 66 (5) ◽  
pp. 1519-1527 ◽  
Author(s):  
Robert Sun ◽  
George A. McMechan
Keyword(s):  
Seismic Waves ◽  
Velocity Model ◽  
P Wave ◽  
Common Source ◽  
Reverse Time ◽  
S Wave ◽  
Wave Source ◽  
Reflected Waves ◽  
S Waves ◽  
Elastic Data

Reflected P‐to‐P and P‐to‐S converted seismic waves in a two‐component elastic common‐source gather generated with a P‐wave source in a two‐dimensional model can be imaged by two independent scalar reverse‐time depth migrations. The inputs to migration are pure P‐ and S‐waves that are extracted by divergence and curl calculations during (shallow) extrapolation of the elastic data recorded at the earth’s surface. For both P‐to‐P and P‐to‐S converted reflected waves, the imaging time at each point is the P‐wave traveltime from the source to that point. The extracted P‐wave is reverse‐time extrapolated and imaged with a P‐velocity model, using a finite difference solution of the scalar wave equation. The extracted S‐wave is reverse‐time extrapolated and imaged similarly, but with an S‐velocity model. Converted S‐wave data requires a polarity correction prior to migration to ensure constructive interference between data from adjacent sources. Synthetic examples show that the algorithm gives satisfactory results for laterally inhomogeneous models.

On plane‐wave decomposition: Alias removal

Geophysics ◽  
1989 ◽  
Vol 54 (10) ◽  
pp. 1339-1343 ◽  
Author(s):  
S. C. Singh ◽  
G. F. West ◽  
C. H. Chapman
Keyword(s):  
Plane Wave ◽  
Seismic Data ◽  
P Wave ◽  
S Wave ◽  
Plane Wave Decomposition ◽  
Time Migration ◽  
Time Distance ◽  
P Domain ◽  
Wave Decomposition ◽  
Rapid Modeling

The delay‐time (τ‐p) parameterization, which is also known as the plane‐wave decomposition (PWD) of seismic data, has several advantages over the more traditional time‐distance (t‐x) representation (Schultz and Claerbout, 1978). Plane‐wave seismograms in the (τ, p) domain can be used for obtaining subsurface elastic properties (P‐wave and S‐wave velocities and density as functions of depth) from inversion of the observed oblique‐incidence seismic data (e.g., Yagle and Levy, 1985; Carazzone, 1986; Carrion, 1986; Singh et al., 1989). Treitel et al. (1982) performed time migration of plane‐wave seismograms. Diebold and Stoffa (1981) used plane‐wave seismograms to derive a velocity‐depth function. Decomposing seismic data also allows more rapid modeling, since it is faster to compute synthetic seismograms in the (τ, p) than in the (t, x) domain. Unfortunately, the transformation of seismic data from the (t, x) to the (τ, p) domain may produce artifacts, such as those caused by discrete sampling, of the data in space.

scholarly journals A new acoustic assumption for orthorhombic media

2020 ◽  
Vol 223 (2) ◽  
pp. 1118-1129
Author(s):  
Mohammad Mahdi Abedi ◽  
Alexey Stovas
Keyword(s):  
Wave Propagation ◽  
Conventional Approach ◽  
P Wave ◽  
Model Parameters ◽  
Governing Equations ◽  
S Wave ◽  
S Waves ◽  
Exploration Seismology ◽  
And Inversion ◽  
Approximate Equations

SUMMARY In exploration seismology, the acquisition, processing and inversion of P-wave data is a routine. However, in orthorhombic anisotropic media, the governing equations that describe the P-wave propagation are coupled with two S waves that are considered as redundant noise. The main approach to free the P-wave signal from the S-wave noise is the acoustic assumption on the wave propagation. The conventional acoustic assumption for orthorhombic media zeros out the S-wave velocities along three orthogonal axes, but leaves significant S-wave artefacts in all other directions. The new acoustic assumption that we propose mitigates the S-wave artefacts by zeroing out their velocities along the three orthogonal symmetry planes of orthorhombic media. Similar to the conventional approach, our method reduces the number of required model parameters from nine to six. As numerical experiments on multiple orthorhombic models show, the accuracy of the new acoustic assumption also compares well to the conventional approach. On the other hand, while the conventional acoustic assumption simplifies the governing equations, the new acoustic assumption further complicates them—an issue that emphasizes the necessity of simple approximate equations. Accordingly, we also propose simpler rational approximate phase-velocity and eikonal equations for the new acoustic orthorhombic media. We show a simple ray tracing example and find out that the proposed approximate equations are still highly accurate.

Sign in / Sign up

Export Citation Format

Share Document