Method for computing an exact impulse response of a plane acoustic reflector at zero offset due to a point acoustic source
First Claim
1. A method for determining impulse response [us(t)] and reflection response [r(t)] of a plane acoustic reflector at zero offset, due to a point acoustic source by varying parameters consisting of reference wave speed about which speed perturbation takes place (c), change in wave speed across a reflector (Δ
- c), plane-wave-reflection coefficient of a reflector for normal incidence (R), depth to the reflector from the source/receiver point (h) and source time function s(t), employing Fourier method, wherein perturbation parameter is substituted in terms of reference wave speed and change in wave speed at an infinite plane interface across which the wave speed changes from the reference wave speed to a sum of the reference wave speed and the change in wave speed, followed by removing the restrictions of Born approximation, the method being useful for seismological applications, the method comprising the steps of;
i) inputting distance of the reflector (h), reference wave speed (c), change in wave speed (Δ
c), source time function s(t) and truncation time of the impulse response of a reflector (T) in equation A;
us(t)=R/8π
h[δ
(t−
{hacek over (c)})−
2{hacek over (c)}2H(t−
{hacek over (c)})/t3]
(Equation A);
ii) obtaining impulse response uS (t) of a plane acoustic reflector at zero offset due to a point acoustic source according to
us(t)=R/8π
h[δ
(t−
{hacek over (c)})−
2{hacek over (c)}2H(t−
{hacek over (c)})/t3],wherein;
R=Δ
c/(2c+Δ
c);
t=observation time that begins at the onset of the source;
δ
(t)=Dirac delta function;
c=Reference wave speed;
Δ
c=change in wave speed;
h=shortest distance between the source and the reflector;
{hacek over (c)}=2h/c;
H=Heaviside function, denoted by H (t−
{hacek over (c)}) and is defined as
H(t−
{hacek over (c)})=1,t>
{hacek over (c)}=0,t<
{hacek over (c)}iii) obtaining reflection response [r(t)] of a plane acoustic reflector at zero offset due to a point acoustic source described by source time function s (t), from the expression as obtained in Equation B;
r(t)=R/8π
h{s(t)−
2{hacek over (c)}2s(t)*[H(t−
{hacek over (c)})/t3]}
(Equation B);
wherein;
R=Δ
c/(2c+Δ
c);
{hacek over (c)}=2h/c; and
* denotes convolution;
H=Heaviside function;
t=observation time that begins at the onset of the source;
s(t)=source time function;
c=Reference wave speed;
Δ
c=change in wave speed; and
h=shortest distance between the source and the reflector, wherein a visualization of the transformed data is output.
1 Assignment
0 Petitions
Accused Products
Abstract
Originating from a novel and an exact algebraic formula for the impulse response of a plane acoustic reflector at zero offset due to a point acoustic source the present invention provides a method for computing an exact impulse response of a plane acoustic reflector at zero offset due to a point acoustic source; and originating from the method, methods for testing and validating algorithms for numerical modeling of seismic reflection, seismic migration and seismic inversion; a method for testing the efficacy of ray-theoretical solution for a given source-reflector configuration; another method for computing zero-offset reflection response of a circular reflector at its central axis; yet another method for validating an interpretation of a reflector as a planar structure; still yet another method for estimating the seismic source-time function when the zero-offset reflection response of a plane reflector is given. Although the algebraic formula and the methods originating from it are, in a strict sense, valid for an acoustic earth and an acoustic source, these would also be of immense utility in the seismic industry where the earth is successfully approximated as an acoustic medium and a seismic source as an acoustic source.
-
Citations
5 Claims
-
1. A method for determining impulse response [us(t)] and reflection response [r(t)] of a plane acoustic reflector at zero offset, due to a point acoustic source by varying parameters consisting of reference wave speed about which speed perturbation takes place (c), change in wave speed across a reflector (Δ
- c), plane-wave-reflection coefficient of a reflector for normal incidence (R), depth to the reflector from the source/receiver point (h) and source time function s(t), employing Fourier method, wherein perturbation parameter is substituted in terms of reference wave speed and change in wave speed at an infinite plane interface across which the wave speed changes from the reference wave speed to a sum of the reference wave speed and the change in wave speed, followed by removing the restrictions of Born approximation, the method being useful for seismological applications, the method comprising the steps of;
i) inputting distance of the reflector (h), reference wave speed (c), change in wave speed (Δ
c), source time function s(t) and truncation time of the impulse response of a reflector (T) in equation A;
us(t)=R/8π
h[δ
(t−
{hacek over (c)})−
2{hacek over (c)}2H(t−
{hacek over (c)})/t3]
(Equation A);ii) obtaining impulse response uS (t) of a plane acoustic reflector at zero offset due to a point acoustic source according to
us(t)=R/8π
h[δ
(t−
{hacek over (c)})−
2{hacek over (c)}2H(t−
{hacek over (c)})/t3],wherein; R=Δ
c/(2c+Δ
c);t=observation time that begins at the onset of the source; δ
(t)=Dirac delta function;c=Reference wave speed; Δ
c=change in wave speed;h=shortest distance between the source and the reflector; {hacek over (c)}=2h/c; H=Heaviside function, denoted by H (t−
{hacek over (c)}) and is defined as
H(t−
{hacek over (c)})=1,t>
{hacek over (c)}=0,t<
{hacek over (c)}iii) obtaining reflection response [r(t)] of a plane acoustic reflector at zero offset due to a point acoustic source described by source time function s (t), from the expression as obtained in Equation B;
r(t)=R/8π
h{s(t)−
2{hacek over (c)}2s(t)*[H(t−
{hacek over (c)})/t3]}
(Equation B);wherein; R=Δ
c/(2c+Δ
c);{hacek over (c)}=2h/c; and
* denotes convolution;H=Heaviside function; t=observation time that begins at the onset of the source; s(t)=source time function; c=Reference wave speed; Δ
c=change in wave speed; andh=shortest distance between the source and the reflector, wherein a visualization of the transformed data is output. - View Dependent Claims (2, 3, 4, 5)
- c), plane-wave-reflection coefficient of a reflector for normal incidence (R), depth to the reflector from the source/receiver point (h) and source time function s(t), employing Fourier method, wherein perturbation parameter is substituted in terms of reference wave speed and change in wave speed at an infinite plane interface across which the wave speed changes from the reference wave speed to a sum of the reference wave speed and the change in wave speed, followed by removing the restrictions of Born approximation, the method being useful for seismological applications, the method comprising the steps of;
Specification