APPARATUS AND METHOD FOR CALCULATING EFFICIENT 3D TRAVELTIME BY USING COARSE-GRID MESH FOR SHALLOW DEPTH SOURCE

• US 20180246240A1
• Filed: 09/26/2017
• Published: 08/30/2018
• Est. Priority Date: 02/28/2017
• Status: Active Grant

• Current Assignee: Korea Institute of Geoscience and Mineral Resources
• Original Assignee: Korea Institute of Geoscience and Mineral Resources
• Inventors: Son, Woohyun, Pyun, Sukjoon, Lee, Ho-Young, Koo, Nam-Hyung, Shin, Changsoo
• Granted Patent: US 10,451,756 B2
• Time in Patent Office: Days
• CPC Class Codes: G01V 1/305 Travel times
• CPC Class Codes: G01V 2210/6222 Velocity; travel time
• CPC Class Codes: G01V 2210/675 Wave equation; Green's functions
• CPC Class Codes: G06F 30/00 Computer-aided design [CAD]
• CPC Class Codes: G06F 30/20 Design optimisation, verification or simulation optimisation, verification or simulation of circuit designs G06F30/30

Abstract

The present invention relates generally to an apparatus and method for calculating efficient 3-dimensional (3D) traveltime by using coarse-grid mesh for a shallow depth source. More particularly, the present invention relates to an efficient 3D traveltime calculation method for a shallow depth source by combining a suppressed wave equation estimation of traveltime (SWEET) algorithm and an equivalent source distribution (ESD) algorithm, wherein the SWEET algorithm is a traveltime calculation algorithm using an damped wave equation and the ESD algorithm is for equivalently distributed sources; and to an apparatus and method for calculating efficient 3D traveltime by using coarse-grid mesh for a shallow depth source which may need less calculation time compared with that of a conventional SWEET algorithm.

Claims

1. An apparatus for calculating efficient 3D traveltime by using coarse-grid mesh for a shallow depth source, the apparatus comprising:

a Green'"'"'s function calculation unit configured to calculate a Green'"'"'s function for a homogeneous half-space medium;

an equivalent source vector calculation unit configured to calculate an equivalent source vector equivalent to an original point source vector by using a wavefield vector sampled at coarse-grid points calculated by the Green'"'"'s function calculation unit;

a wavefield vector calculation unit configured to calculate the wavefield vector by using the equivalent source vector calculated by the equivalent source vector calculation unit;

a wavefield vector'"'"'s partial derivative generation unit configured to generate a partial derivative of the wavefield vector calculated by the wavefield vector calculation unit; and

a first-arrival traveltime calculation unit configured to calculate a first-arrival traveltime by a suppressed wave equation estimation of traveltime (SWEET) algorithm with both the wavefield vector calculated by the wavefield vector calculation unit and the partial derivative of the wavefield vector generated by the wavefield vector'"'"'s partial derivative generation unit.

2. The apparatus of claim 1, wherein the Green'"'"'s function calculation unit calculates the Green'"'"'s function by the following equation 7:

3. The apparatus of claim 2, wherein the equivalent source vector calculation unit calculates the equivalent source vector by the following equation 8:

f^equi=Sũ

, 

 

(Equation

     8)

where f^equi is a new equivalent source vector for ahomogeneous half-space, S is an impedance matrix, and ũ

is the wavefield vector sampled at the coarse-grid points from the analytical solution of the equation 7.

4. The apparatus of claim 3, wherein the wavefield vector calculation unit calculates the wavefield vector by the following equation 9:

u^equi=S<sup>−

1</sup>f^equi, 

 

(Equation

     9)

where u^equi is the wavefield vector generated from the equivalent source vector.

5. The apparatus of claim 4, wherein the wavefield vector'"'"'s partial derivative generation unit generates the partial derivative of the wavefield vector by the following equation 10:

6. The apparatus of claim 5, wherein the first-arrival traveltime calculation unit calculates the first-arrival traveltime by the following equation 11:

7. A method of calculating efficient 3D traveltime by using coarse-grid mesh for a shallow depth source by using the apparatus of claim 1, the method comprising:

calculating, by a Green'"'"'s function calculation unit, a Green'"'"'s function for a homogeneous half-space medium;

calculating, by an equivalent source vector calculation unit, an equivalent source vector equivalent to an original point source vector by using a wavefield vector sampled at coarse-grid points calculated at the calculating the Green'"'"'s function;

calculating, by a wavefield vector calculation unit, a wavefield vector by using the equivalent source vector calculated at the calculating the equivalent source vector;

generating, by the wavefield vector'"'"'s partial derivative generation unit, a partial derivative of the wavefield vector calculated at the calculating the wavefield vector; and

calculating, by a first-arrival traveltime calculation unit, a first-arrival traveltime by using the SWEET algorithm with both the wavefield vector calculated at the calculating the wavefield vector and the partial derivative of the wavefield vector generated at the generating the partial derivative of the wavefield vector.

8. The method of claim 7, wherein at the calculating the Green'"'"'s function, the Green'"'"'s function is calculated by the following equation 7:

9. The method of claim 8, wherein at the calculating the equivalent source vector, the equivalent source vector is calculated by the following equation 8:

f^equi=Sũ

, 

 

(Equation

     8)

where f^equi is a new equivalent source vector for the homogeneous half-space, S is an impedance matrix, and ũ

is the wavefield vector sampled at the coarse-grid points from the analytical solution of the equation 7.

10. The method of claim 9, wherein at the calculating the wavefield vector, the wavefield vector is calculated by the following equation 9:

u^equi=S<sup>−

1</sup>f^equi, 

 

(Equation

     9)

where u^equi is the wavefield vector generated from the equivalent source vector.

11. The method of claim 10, wherein at the generating the partial derivative of the wavefield vector, the partial derivative of the wavefield vector is generated by the following equation 10:

12. The method of claim 11, wherein at the calculating the first-arrival traveltime, the first-arrival traveltime is calculated by the following equation 11: