 Library Home /
 Search Collections /
 Open Collections /
 Browse Collections /
 UBC Theses and Dissertations /
 Truncated asymptotic solution of the onedimensional...
Open Collections
UBC Theses and Dissertations
UBC Theses and Dissertations
Truncated asymptotic solution of the onedimensional inhomogeneous wave equation Zelt, Barry Curtis
Abstract
I present a new timedomain method for solving for the stress and particle velocity of normally incident plane waves propagating in a smoothly varying onedimensional medium. Both the Young's modulus E and the density ⍴ are allowed to vary smoothly with depth. The restriction of geometrical optics, that the wavelength be much less than the material stratification length, is not required in this new method. The infinite geometrical optics expansion is truncated after n terms, imposing a condition on the acoustic impedance I for exact solutions to exist. For the case ռ = 2 there are three general classes of impedance functions for which the resultant expansion is uniform and exact. To check the numerical validity of the "truncated asymptotic" (TA) solution results are calculated for the case of a medium with a linear velocity gradient for which there is an exact solution in the frequency domain. Since a linear velocity gradient is not one of the foregoing classes of impedance functions, a curvefitting approach is necessary. The TA method compares favourably to the exact solution and is accurate to within the error of the required curve fit. Two classes of synthetic seismograms are calculated for smooth velocity and density variations. The same impedance as a function of traveltime is used for both classes. In the first class the principal variation in impedance is due to velocity, while in the second it is mainly due to density. The amplitudes in both classes of synthetic seismograms are very similar, but, as expected, the traveltime curves for each class are widely separated. For the case ⍴ = constant the TA solution is used as a benchmark to analyse a twoterm WKBJ approximation for three classes of velocity functions. The velocity functions are such that the TA solutions are exact. For two of the classes the WKBJ solution performs well when the length of the transition zone is of the same order, or larger, than the length of the applied wavelet. For steeper velocity gradients the WKBJ solution begins to differ significantly from the exact TA solution. The WKBJ solution for the third class performs extremely well even for steep gradients. Equations governing the validity of the WKBJ solution are examined to explain the above results. Equations are derived to describe the distortion of a stress pulse propagating through a transition zone. For small velocity gradients (relative to the length of the applied pulse) the wavelet changes amplitude but its phase is not effected. As the gradient increases and the velocity function becomes a discontinuity at z = 0 the wavelet travels through undistorted. Only when the transition zone width is of the order of the length of the wavelet is there any visible phase distortion. Reflection and transmission coefficients as functions of time are calculated for low, intermediate and high gradient transition zones. The transmission coefficient is a delta function in each case. The reflection coefficient has the shape of a Hilbert transform for low gradients. For higher gradients the reflection coefficient approaches the shape of a delta function.
Item Metadata
Title 
Truncated asymptotic solution of the onedimensional inhomogeneous wave equation

Creator  
Publisher 
University of British Columbia

Date Issued 
1987

Description 
I present a new timedomain method for solving for the stress and particle velocity of normally incident plane waves propagating in a smoothly varying onedimensional medium. Both the Young's modulus E and the density ⍴ are allowed to vary smoothly with depth. The restriction of geometrical optics, that the wavelength be much less than the material stratification length, is not required in this new method. The infinite geometrical optics expansion is truncated after n terms, imposing a condition on the acoustic impedance I for exact solutions to exist. For the case ռ = 2 there are three general classes of impedance functions for which the resultant expansion is uniform and exact.
To check the numerical validity of the "truncated asymptotic" (TA) solution results are calculated for the case of a medium with a linear velocity gradient for which there is an exact solution in the frequency domain. Since a linear velocity gradient is not one of the foregoing classes of impedance functions, a curvefitting approach is necessary. The TA method compares favourably to the exact solution and is accurate to within the error of the required curve fit.
Two classes of synthetic seismograms are calculated for smooth velocity and density variations. The same impedance as a function of traveltime is used for both classes. In the first class the principal variation in impedance is due to velocity, while in the second it is mainly due to density. The amplitudes in both classes of synthetic seismograms are very similar, but, as expected, the traveltime curves for each class are widely separated.
For the case ⍴ = constant the TA solution is used as a benchmark to analyse a twoterm WKBJ approximation for three classes of velocity functions. The velocity functions are such that the TA solutions are exact. For two of the classes the WKBJ solution performs well when the length of the transition zone is of the same order, or larger, than the length of the applied wavelet. For steeper velocity gradients the WKBJ solution begins to differ significantly from the exact TA solution. The WKBJ solution for the third class performs extremely well even for steep gradients. Equations governing the validity of the WKBJ solution are examined to explain the above results.
Equations are derived to describe the distortion of a stress pulse propagating through a transition zone. For small velocity gradients (relative to the length of the applied pulse) the wavelet changes amplitude but its phase is not effected. As the gradient increases and the velocity function becomes a discontinuity at z = 0 the wavelet travels through undistorted. Only when the transition zone width is of the order of the length of the wavelet is there any visible phase distortion.
Reflection and transmission coefficients as functions of time are calculated for low, intermediate and high gradient transition zones. The transmission coefficient is a delta function in each case. The reflection coefficient has the shape of a Hilbert transform for low gradients. For higher gradients the reflection coefficient approaches the shape of a delta function.

Genre  
Type  
Language 
eng

Date Available 
20100720

Provider 
Vancouver : University of British Columbia Library

Rights 
For noncommercial purposes only, such as research, private study and education. Additional conditions apply, see Terms of Use https://open.library.ubc.ca/terms_of_use.

DOI 
10.14288/1.0052971

URI  
Degree  
Program  
Affiliation  
Degree Grantor 
University of British Columbia

Campus  
Scholarly Level 
Graduate

Aggregated Source Repository 
DSpace

Item Media
Item Citations and Data
Rights
For noncommercial purposes only, such as research, private study and education. Additional conditions apply, see Terms of Use https://open.library.ubc.ca/terms_of_use.