Extension of forward modeling phase-screen code in isotropic and anisotropic media up to critical angle

Geophysics ◽  
2007 ◽  
Vol 72 (5) ◽  
pp. SM107-SM114 ◽  
Author(s):  
James C. White ◽  
Richard W. Hobbs

The computationally efficient phase-screen forward modeling technique is extended to allow investigation of nonnormal raypaths. The code is developed to accommodate all diffracted and converted phases up to critical angle, building on a geometric construction method. The new approach relies upon prescanning the model space to assess the complexity of each screen. The propagating wavefields are then divided as a function of horizontal wavenumber, and each subset is transformed to the spatial domain separately, carrying with it angular information. This allows both locally accurate 3D phase corrections and Zoeppritz reflection and transmission coefficients to be applied. The phase-screen code is further developed to handle simple anisotropic media. During phase-screen modeling, propagation is undertaken in the wavenumber domain where exact expressions for anisotropic phase velocities are available. Traveltimes and amplitude effects from a range of anisotropic shales are computed and compared with previous published results.

1978 ◽  
Vol 1 (16) ◽  
pp. 133
Author(s):  
Robert B. Steimer ◽  
Charles K. Sollitt

Previous attempts to analytically describe wave reflection and transmission at surface penetrating structures have neglected losses due to flow expansion, contraction, and skin drag along the structure boundaries (Black and Mei, 1970; Ijima, et al., 1972). The model described in this study includes these effects and allows for the inclusion of a dissipative medium such as rubble or closely spaced piles in the region beneath the structure. The problem of a fixed, two-dimensional structure in a train of monochromatic incident waves is modeled, as shown in Figure 1. The solution allows for 1) variable structure length and draft, 2) different depths in the regions fore, aft, and beneath the structure, 3) variable wave amplitude and period, and 4) turbulent and inertial damping in the region beneath the structure. An equivalent work technique is applied to linearize the damping beneath the structure, yielding a potential flow problem in all three regions. Amplitudes for the resulting series of eigenfunctions in each region are determined by matching pressure and horizontal mass flux at the region interfaces, orthogonalizing these expressions over the depth, and simultaneously solving the resulting equations to yield complex reflection and transmission coefficients. Complex horizontal and vertical force coefficients for the structure are also determined from the integrated Bernoulli equation. The solution technique is computationally efficient. In general, five modes in the eigen series provide satisfactory convergence for the various hydrodynamic parameters. Approximately six-tenths of a computer system second are required to solve for a single wave-structure condition. The results compare favorably with variational methods used by others.


Geophysics ◽  
1979 ◽  
Vol 44 (1) ◽  
pp. 27-38 ◽  
Author(s):  
P. F. Daley ◽  
F. Hron

The deficiency of an isotropic model of the earth in the explanation of observed traveltime phenomena has led to the mathematical investigation of elastic wave propagation in anisotropic media. A type of anisotropy dealt with in the literature (Potsma, 1955; Cerveny and Psencik, 1972; and Vlaar, 1968) is uniaxial anisotropy or transverse isotropy. A special case of transverse isotropy which assumes the wavefronts to be ellipsoids of revolution has been used by Cholet and Richard (1954) and Richards (1960) in accounting for the observed traveltimes at Berraine in the Sahara and in the foothills of Western Canada. The kinematics of this problem have been treated in a number of papers, the most notable being Gassmann (1964). However, to appreciate fully the effect of anisotropy, the dynamics of the problem must be explored. Assuming a model of the earth consisting of plane transversely isotropic layers with the axes of anisotropy perpendicular to the interfaces, a prime requisite for obtaining amplitude distance curves or synthetic seismograms is the calculation of reflection and transmission coefficients at the interfaces. In this paper the special case of ellipsoidal anisotropy will be considered. That the quasi‐shear SV wavefront is forced to be spherical by this assumption is unfortunate, but it is instructive to investigate this simple anisotropic model as it incorporates many features inherent to wave propagation in a more general anisotropic medium. A brief outline of the theory for wave propagation in an ellipsoidally anisotropic medium is given and the analytic expressions for the reflection and transmission coefficients are listed. A complete derivation of reflection and transmission coefficients in transversely isotropic media can be found in Daley and Hron (1977). Finally, all 24 possible reflection and transmission coefficients and surface conversion coefficients are displayed for varying degrees of anisotropy.


2021 ◽  
Vol 40 (4) ◽  
pp. 267-276
Author(s):  
Peter Mesdag ◽  
Leonardo Quevedo ◽  
Cătălin Tănase

Exploration and development of unconventional reservoirs, where fractures and in-situ stresses play a key role, call for improved characterization workflows. Here, we expand on a previously proposed method that makes use of standard isotropic modeling and inversion techniques in anisotropic media. Based on approximations for PP-wave reflection coefficients in orthorhombic media, we build a set of transforms that map the isotropic elastic parameters used in prestack inversion into effective anisotropic elastic parameters. When used in isotropic forward modeling and inversion, these effective parameters accurately mimic the anisotropic reflectivity behavior of the seismic data, thus closing the loop between well-log data and seismic inversion results in the anisotropic case. We show that modeling and inversion of orthorhombic anisotropic media can be achieved by superimposing effective elastic parameters describing the behavior of a horizontally stratified medium and a set of parallel vertical fractures. The process of sequential forward modeling and postinversion analysis is exemplified using synthetic data.


2008 ◽  
Vol 130 (3) ◽  
Author(s):  
J. Y. Chen ◽  
H. L. Chen ◽  
E. Pan

Reflection and transmission coefficients of plane waves with oblique incidence to a multilayered system of piezomagnetic and/or piezoelectric materials are investigated in this paper. The general Christoffel equation is derived from the coupled constitutive and balance equations, which is further employed to solve the elastic displacements and electric and magnetic potentials. Based on these solutions, the reflection and transmission coefficients in the corresponding layered structures are subsequently obtained by virtue of the propagator matrix method. Two layered examples are selected to verify and illustrate our solutions. One is the purely elastic layered system composed of aluminum and organic glass materials. The other layered system is composed of the novel magnetoelectroelastic material and the organic glass. Numerical results are presented to demonstrate the variation of the reflection and transmission coefficients with different incident angles, frequencies, and boundary conditions, which could be useful to nondestructive evaluation of this novel material structure based on wave propagations.


2016 ◽  
Vol 2016 ◽  
pp. 1-6 ◽  
Author(s):  
Yunhua Cao ◽  
Haiying Li ◽  
Zhe Wang ◽  
Zhensen Wu

Propagation characteristics of oblique incident terahertz wave from the nonuniform dusty plasma are studied using the propagation matrix method. Assuming that the electron density distribution of dusty plasma is parabolic model, variations of power reflection, transmission, and absorption coefficients with frequencies of the incident wave are calculated as the wave illuminates the nonuniform dusty plasma from different angles. The effects of incident angles, number density, and radius of the dust particles on propagation characteristics are discussed in detail. Numerical results show that the number density and radius of the dust particles have very little influences on reflection and transmission coefficients and have obvious effects on absorption coefficients. The terahertz wave has good penetrability in dusty plasma.


Sign in / Sign up

Export Citation Format

Share Document