KR101564094B1 - Elastic reverse-time migration system and method using absolute value function for improving the quality of subsurface structure imaging - Google Patents

Elastic reverse-time migration system and method using absolute value function for improving the quality of subsurface structure imaging Download PDF

Info

Publication number
KR101564094B1
KR101564094B1 KR1020150094460A KR20150094460A KR101564094B1 KR 101564094 B1 KR101564094 B1 KR 101564094B1 KR 1020150094460 A KR1020150094460 A KR 1020150094460A KR 20150094460 A KR20150094460 A KR 20150094460A KR 101564094 B1 KR101564094 B1 KR 101564094B1
Authority
KR
South Korea
Prior art keywords
wave
stress
displacement
back propagation
sound source
Prior art date
Application number
KR1020150094460A
Other languages
Korean (ko)
Inventor
하지호
정우근
신성렬
김영준
김원식
Original Assignee
한국지질자원연구원
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 한국지질자원연구원 filed Critical 한국지질자원연구원
Priority to KR1020150094460A priority Critical patent/KR101564094B1/en
Application granted granted Critical
Publication of KR101564094B1 publication Critical patent/KR101564094B1/en
Priority to US15/196,841 priority patent/US9798026B2/en

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. analysis, for interpretation, for correction
    • G01V1/284—Application of the shear wave component and/or several components of the seismic signal
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. analysis, for interpretation, for correction
    • G01V1/282—Application of seismic models, synthetic seismograms
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. analysis, for interpretation, for correction
    • G01V1/30—Analysis
    • G01V1/301—Analysis for determining seismic cross-sections or geostructures
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. analysis, for interpretation, for correction
    • G01V1/30—Analysis
    • G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/10—Aspects of acoustic signal generation or detection
    • G01V2210/12—Signal generation
    • G01V2210/125—Virtual source
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/20—Trace signal pre-filtering to select, remove or transform specific events or signal components, i.e. trace-in/trace-out
    • G01V2210/23—Wavelet filtering
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/50—Corrections or adjustments related to wave propagation
    • G01V2210/51—Migration
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/60—Analysis
    • G01V2210/61—Analysis by combining or comparing a seismic data set with other data
    • G01V2210/616—Data from specific type of measurement
    • G01V2210/6161—Seismic or acoustic, e.g. land or sea measurements
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/60—Analysis
    • G01V2210/66—Subsurface modeling
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/60—Analysis
    • G01V2210/67—Wave propagation modeling
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/60—Analysis
    • G01V2210/67—Wave propagation modeling
    • G01V2210/679—Reverse-time modeling or coalescence modelling, i.e. starting from receivers
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/70—Other details related to processing
    • G01V2210/74—Visualisation of seismic data

Abstract

The present invention relates to an apparatus and method for compensating an inverse temporal structure in an elastic medium to which an absolute value function is applied to enhance the quality of an underground image, and more particularly, And more particularly, to an apparatus and method for compensating inverse temporal structure in an elastic medium which improves the accuracy of an underground image by applying an absolute value function to the structure-corrected image.

Description

BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to an apparatus and method for compensating an inverse temporal structure in an elastic medium to which an absolute value function is applied to enhance the quality of an underground image,

The present invention relates to an apparatus and method for compensating an inverse temporal structure in an elastic medium to which an absolute value function is applied to enhance the quality of an underground image, and more particularly, And more particularly, to an apparatus and method for compensating inverse temporal structure in an elastic medium which improves the accuracy of an underground image by applying an absolute value function to the structure-corrected image.

Seismic exploration is a method of grasping underground structures by artificially generating seismic waves in the surface or water and analyzing the propagation behavior of the generated seismic waves. However, when a strata has a slope or a complex geological structure, the reflection point of the seismic waves has a problem of being distorted or diffracted. Therefore, for accurate underground imaging using seismic data, it is necessary to perform a structural migration, which is a process of moving a distorted reflection event to an accurate reflection position.

In the structural correction, Kirchhoff migration and reverse-time migration using wave equation are mainly used. In particular, inverse time structure correction using a two-way wave equation Since the inverse temporal structure correction using the bidirectional wave equation requires a large amount of computation and relatively long computation time, it is possible to obtain an image with high accuracy even in complicated terrain, media).

However, recently, as computing devices such as supercomputers have become widespread, inverse time structure corrections have been applied to elastic media that can be reflected more closely to actual conditions. In the elastic medium, various waves such as P wave, S wave, and PS wave are propagated unlike in the acoustic medium, and radio waves between various waves act as noise between individual signals. .

Therefore, it is necessary to separate the wave fields of various components in the inverse time structure correction. Particularly, many studies have been made on the separation of the P wave wave field capable of propagating to deep depths. For this purpose, wave field separation using a divergence operator has been performed.

 Since then, research has proceeded in the direction of anisotropy and three-dimensional data processing. However, in the case of the wave field separation technique using the diversion operator, the spatial differential of the stress component is required, so that the waveform of the wave field is deformed. It is also necessary to study for a clear interpretation of underground strata. Especially, the lower stratum structure of the high-speed layer such as rock salt, which is a representative form of the reservoir with oil and gas resources, is an important factor that influences the shape and size of the reservoir. There is a limit. Therefore, there is a need for a way to overcome these problems.

Korean Patent No. 10-1413751 (registered on June 24, 2014)

It is an object of the present invention to provide an apparatus and method for compensating inverse time structure in an elastic medium which can minimize waveform deformation of a wave field in performing separation of a wave field in structural correction, Method.

It is still another object of the present invention to provide an apparatus and method for compensating inverse time structure in an elastic medium to which an absolute value function capable of further improving the quality of a structure corrected image of an underground structure is applied by applying an absolute value function.

The apparatus for correcting inverse temporal structure in an elastic medium according to the present invention includes: a virtual source calculation unit for calculating a virtual source by receiving seismic exploration data and an underground structure model; A virtual sound source wave length calculation unit for calculating a virtual sound source wave length from the calculated virtual transmission source; A back propagation wave length calculation unit for receiving the seismic wave exploration data and processing back propagation and calculating a back propagation wave length; A virtual sound source wave - length waveform demultiplexer for separating the P - wave sound field from the virtual sound source wave field by using a relationship between stress and displacement; A backwave wave long wave form separator for separating a P wave wave field from the back propagation wave field using a relationship between stress and displacement; And a convolution unit for convoluting the two wave fields separated by the virtual sound source wave length separation unit and the back propagation wave length separation unit.

The apparatus may further include an image quality enhancing unit that applies an absolute value function to the structure-corrected seismic image acquired through convolution in the convolution unit.

Also, the inverse temporal structure correction method in the elastic medium of the present invention includes the steps of: a) receiving seismic exploration data and underground structure model data; b) calculating a virtual transmission source based on the input data and calculating a virtual sound source wave field; c) calculating a P wave back propagation wave length by back propagation processing the seismic wave survey data; d) separating the P-wave wave field from the virtual sound source wave field by using a relationship between stress and displacement; e) separating the wave field from the back propagation wave field using a relationship between stress and displacement; And f) convoluting the two wave fields separated in the step c) and the step c), respectively.

G) applying an absolute value function to the structure-corrected seismic image acquired through convolution in step f).

The present invention utilizes the relationship between stress and displacement when separating a P wave wave field, and it has an advantage of minimizing a numerical error that can be generated by using a variation operator when a conventional P wave wave field is separated .

In addition, the present invention has an advantage of obtaining an image of an improved underground structure in which an interlayer boundary is more clearly displayed by applying an absolute value function.

1 is a configuration diagram of a structure correcting apparatus according to the present invention;
2 is a flow chart of a structure correction method according to the present invention;
Fig. 3 is a graph showing the relationship between the horizontal component wave field (a), the vertical component wave field (b), and the P wave field (c) separated by using the divergence in the elastic medium. (g), (h), and (i) on the 150th channel.
4 is a P wave field (c) separated by using a relational expression of the displacement and stress of the present invention, and the wave component field (G), (h), and (i) in the 150th channel, respectively, of the corresponding elastic wave record (d), (e), and (f).
Fig. 5 is a view showing a P wave wave field (a) separated using a conventional divergence and a P wave wave field (b) separated using a relational expression of displacement and stress according to the present invention;
6 shows a signal (b) which takes an arbitrary sine signal (a) and an envelope and absolute value function;
7 is a diagram showing SEG / EAGE rock salt P-wave model (a), S wave model (b) and density model (c).
Fig. 8 is a diagram showing SEG / EAGE saline-pom P wave model (a), S wave model (b) and density model (c) smoothed for application to the invention.
9 is an inverse time structure correction image using a conventional structure correction method.
10 is an inverse time structure correction image using the structure correction method according to the present invention.

Hereinafter, the technical idea of the present invention will be described more specifically with reference to the accompanying drawings.

BRIEF DESCRIPTION OF THE DRAWINGS The accompanying drawings, which are included to provide a further understanding of the technical concept of the present invention, are incorporated in and constitute a part of the specification, and are not intended to limit the scope of the present invention.

Fig. 1 is a configuration diagram of a structure correcting apparatus according to the present invention, and Fig. 2 is a flowchart of a structure correcting method according to the present invention.

The apparatus for correcting inverse temporal structure in an elastic medium according to the present invention includes a structure correcting unit, and the structure correcting unit performs inverse time structure correction using a bi-directional wave equation to the elastic medium. 1, the structure modification unit 100 includes a virtual transmission source calculation unit 110, a virtual sound source wave field calculation unit 120, a back propagation wave field calculation unit 130, a virtual transmission source wave- A back propagation wave long-wave separator 150, and a convolution unit 160. The back- At this time, the virtual transmission source calculating unit 110, the virtual sound source wave length calculating unit 120, and the back propagation wave length calculating unit 130 are well-known technologies, and a detailed description is omitted here.

First, the seepage wave survey data and the underground structure model are input to the virtual transmission source calculation unit 110. (S100) The virtual transmission source calculation unit 110 calculates a virtual transmission source based on the input data, (Virtual transmission source) wave field from the calculated virtual transmission source (S200). Also, the back propagation wave length calculation unit 130 receives the seismic wave exploration data, processes the back propagation wave, and calculates the back propagation wave length from the back propagation processed data (S300).

Thereafter, the virtual sound source wave-length waveform separator 140 separates the P-wave wave field from the calculated virtual sound source wave field (S400), and the back-propagation wave-length waveform separator 150 separates the P- Wave field separation is performed (S500). At this time, in the present invention, the wave field separation is performed by using the relationship between the stress and the displacement of each wave field.

As described above, in the inverse temporal structure correction, a process (S400, S500) of separating the wave fields of various components is required in order to minimize the influence of noise and obtain an accurate underground structure image. A method of using a divergence operator as a conventional wave field separation technique has been mainly used, but there has been a problem that a waveform is deformed due to application of a differential to a space. In the present invention, in order to solve the above-mentioned problem, waveform separation is performed using the relationship between stress and displacement instead of the die version operator. That is, since the present invention uses a stress tensor to which a spatial differential is applied, there is an advantage that no deformation of the waveform occurs after the wave field separation. A specific method for this will be described later.

Thereafter, the virtual source P wave field and the backward propagation P wave field separated by the virtual transmission source wave length separator 140 and the backward wave long wave form separator 150 are convoluted through the convolution unit 160, (S600), thereby obtaining a structure-corrected image.

Meanwhile, the structure correcting apparatus according to the present invention may further include an image quality improving unit 200 for improving the accuracy of the image obtained through the structure correcting unit 100. The image quality enhancement unit 200 applies an absolute value function to the structure-compensated seismic image acquired through convolution in the convolution unit 160 (S700), thereby obtaining a clear inter- There is an advantage in obtaining images of underground structures.

Finally, the scaling unit 300 scales the underground structure image (S800) to finally obtain the structure-corrected seismic image.

Hereinafter, this will be described in more detail with reference to mathematical equations. In the inverse time structure correction in the elastic medium, the process of separating the waves should be given priority. Generally, a method using Helmholtz decomposition is mainly used for wave separation, and the displacement of an acoustic wave is expressed as Equation 1 below.

Figure 112015064257149-pat00001

here

Figure 112015064257149-pat00002
Is a displacement vector,
Figure 112015064257149-pat00003
and
Figure 112015064257149-pat00004
Are curl and divergence components, respectively. Using this, the P-wave and the S-wave wave length can be expressed by the following equations (2) and (3).

Figure 112015064257149-pat00005

Figure 112015064257149-pat00006

In other words, since the wave field can be separated by taking the divergence operator or the curl operator in the displacement of the acoustic wave, the waveform separation is performed using a die version operator in the past. However, in the present invention, waveform separation is performed by using the relationship between stress and displacement as shown in Equations (4) and (5) below (S400, S500) instead of using a die version or a curl operator.

Figure 112015064257149-pat00007

Figure 112015064257149-pat00008

here

Figure 112015064257149-pat00009
,
Figure 112015064257149-pat00010
Represents a Lame constant,
Figure 112015064257149-pat00011
,
Figure 112015064257149-pat00012
Is the stress term,
Figure 112015064257149-pat00013
,
Figure 112015064257149-pat00014
Represents the displacement vector for each x and z axis. In general, the separated P-wave wave field is a displacement vector
Figure 112015064257149-pat00015
,
Figure 112015064257149-pat00016
Can be expressed by the sum of partial derivatives for each axis. This can be expressed as a sum of stress terms as shown in Equation (6) below using the relationship between stress and displacement as shown in Equations (4) and (5). That is, the wave field separation is performed using Equation (6) for the relationship between stress and displacement in the elastic medium of the present invention.

Figure 112015064257149-pat00017

3 and 4, the vertical component and the horizontal component of each wave field in the elastic medium and the separated P wave wave in the elastic medium are obtained by using the waveform separation method using the conventional divergence and the waveform separation method using the stress- You can check the chapter. In particular, when the 150th trace of the elastic wave recording is examined, when waveform separation is performed using the die version, the waveforms of the horizontal and vertical components (g, h) and the waveforms of the separated wave field (i) On the other hand, in the waveform separation using the stress-displacement relationship of the present invention, the waveform (i) does not change.

5 shows waveforms of P wave waves (a) separated by waveforms using conventional diverting and waveforms P wave waves (b) obtained by using the relationship of stress-displacement of the present invention. It can be seen that the waveform separation can be performed by using the relationship between the displacement and the stress rather than the die version which has been conventionally utilized through the coincidence of the results.

That is, the stress-displacement relationship can be used in the separation of the P-wave wave field by the virtual sound source wave-length separator 140 and the back-propagation wave-field separator 150 of the present invention through FIGS. 3 to 5, It can be seen that the waveform change due to the spatial differential generated in the method using the conventional diversity operator is minimized.

On the other hand, in order to apply the waveform separation, the inverse time structure correction in the elastic medium can be expressed as shown in Equation 7 in the time domain.

Figure 112015064257149-pat00018

Where k is the lattice position,

Figure 112015064257149-pat00019
Modeling wave field,
Figure 112015064257149-pat00020
Lt; RTI ID = 0.0 > k <
Figure 112015064257149-pat00021
Is a partial derivative wave field vector,
Figure 112015064257149-pat00022
Is the seismic exploration data acquired in the field. The partial derivative wave field is a convolution of the green function and the virtual sound source vector, and is represented by the following equation (8).

Figure 112015064257149-pat00023

Figure 112015064257149-pat00024
Is a greens function, * is a convolution operator,
Figure 112015064257149-pat00025
Is a virtual sound source vector. Applying Equation (8) to Equation (7), it can be expanded to Equation (9) as follows.

Figure 112015064257149-pat00026

here

Figure 112015064257149-pat00027
Is a zero-lag cross correlation,
Figure 112015064257149-pat00028
Total recording time,
Figure 112015064257149-pat00029
Represents the recorded time. Back propagated wave field (
Figure 112015064257149-pat00030
) Is expressed as Equation (10) below,

Figure 112015064257149-pat00031

Equation (9) can be expressed as Equation (11) using Equation (10).

Figure 112015064257149-pat00032

In addition, Chung, W., Shin, J., Shin, C. and Shin, S. 2012, Elastic reverse-time migration using Helmholtz decomposition in the time domain, SEG Expanded abstract, SEG Las Vegas 2012 Annual Meeting. , A method of expressing Equation (11) using a virtual sound source vector has been proposed. Using the virtual sound source vector, the inverse time structure correction is expressed as Equation (12).

Figure 112015064257149-pat00033

here

Figure 112015064257149-pat00034
Is a virtual sound source vector. Considering the blasting positions in all the sound sources, the structure-corrected image can be expressed as shown in Equation (13), and the waveform separation using the diversion operator can be expressed as Equation (14).
Figure 112015064257149-pat00035
Is the number of times the sound source is blasted.

Figure 112015064257149-pat00036

Figure 112015064257149-pat00037

Equation (14) is utilized in the inverse time structure correction in a conventional elastic medium. On the other hand, in the present invention, since waveform separation is performed using Equation 6, which is a relational expression of stress and displacement, the inverse time structure correction using the relational expression of stress and displacement corresponding to Equation 14 is expressed as Equation 15 below .

Figure 112015064257149-pat00038

here

Figure 112015064257149-pat00039
Is a wave field vector that separates the P wave from the virtual sound source vector. At this time
Figure 112015064257149-pat00040
,
Figure 112015064257149-pat00041
Lt;
Figure 112015064257149-pat00042
Wow
Figure 112015064257149-pat00043
Can be expressed by the following expression (16).

Figure 112015064257149-pat00044

Figure 112015064257149-pat00045
The
Figure 112015064257149-pat00046
, And
Figure 112015064257149-pat00047
Wow
Figure 112015064257149-pat00048
Represents the stress term of the wave field subjected to the back propagation process, respectively.

Also, in the present invention, an absolute value function is applied to the structure-corrected image defined by Equation (15) to obtain an improved structure-corrected image.

In general, envelope functions are used to analyze structural information of underground media such as boundaries and faults between strata during elastic wave complex analysis. The seismic data using this envelope technique is very useful for analyzing the structure of strata because the low frequency components are strengthened and expressed only by the reflection intensity.

 In the present invention, an absolute value function, which has similar characteristics to the envelope signal used in the above-described acoustic wave complex analysis and has a polarity change of the signal more clearly, is applied to the structure correction technique.

Fig. 6 is a diagram showing a signal (b) which takes an arbitrary sine signal (a) and an envelope and absolute value function. As shown in the figure, when the envelope signal and the absolute value function are applied to an arbitrary sine function signal, the amplitude and the wavelength of the envelope signal are similar to each other. However, the amplitude of the signal changes more remarkably at the point where the polarity changes . For this reason, when using the absolute value function, it is advantageous to obtain an image of an improved underground structure in which the interlayer boundary is more clearly displayed.

Here, the structure correction image to which the absolute value function is applied can be expressed as Equation (17).

Figure 112015064257149-pat00049

By using the above equation, a structure-corrected seismic image with an absolute value function can be obtained.

In order to confirm the effect of the present invention, the structure correction method according to the present invention was applied to a SEG / EASE rock salt dome model. The SEG / EASE armillary dome model is widely used for validating related studies and techniques and is considered to be sufficient to confirm the effects of the present invention. The SEG / EAGE armillary dome model used in the present invention is a P wave velocity model (a), an S wave velocity model (b), and a density model (c) of FIG. 7, which are obtained by smoothing the model of FIG. The velocity model was used for the inverse time structure correction according to the present invention.

The image using the conventional inverse time structure correction method is shown in FIG. 9, and the inverse time structure corrected image obtained by the structure correction method of the present invention is shown in FIG. Comparing the two images, the boundaries of the strata are more apparent in the images using the method according to the present invention than in the conventional method. Particularly, in the case of the conventional method, the structure of the lower part of the cancer dome is hardly shown, and in the case of the method according to the present invention, accurate image can be confirmed even under the lower part of the cancer dome. Thus, it can be seen that the structure correction method according to the present invention provides a more improved image than the conventional method.

It will be understood by those skilled in the art that various changes in form and details may be made therein without departing from the spirit and scope of the invention as defined by the appended claims.

100:
110: virtual transmission source calculation unit 120: virtual sound source wave length calculation unit
130: back propagation wave length calculation unit 140: virtual sound source wave length waveform separation unit
150: backward wave long wave form separator 160: convolution unit
200: Image quality enhancement unit using an absolute value function
300: Scaling unit

Claims (4)

A virtual transmission source calculation unit for receiving the seismic wave exploration data and the underground structure model and calculating a virtual transmission source;
A virtual sound source wave length calculation unit for calculating a virtual sound source wave length from the calculated virtual transmission source;
A back propagation wave length calculation unit for receiving the seismic wave exploration data and processing back propagation and calculating a back propagation wave length;
A virtual sound source wave - length waveform demultiplexer for separating the P - wave sound field from the virtual sound source wave field by using a relationship between stress and displacement;
A backwave wave long wave form separator for separating a P wave wave field from the back propagation wave field using a relationship between stress and displacement; And
A convolution unit for convoluting the two wave fields separated by the virtual sound source wave length separation unit and the back propagation wave length separation unit;
And,
The relationship between the stress and the displacement,
Figure 112015084809975-pat00060

Figure 112015084809975-pat00061

(At this time
Figure 112015084809975-pat00062
,
Figure 112015084809975-pat00063
Represents a Lame constant,
Figure 112015084809975-pat00064
,
Figure 112015084809975-pat00065
Is the stress term,
Figure 112015084809975-pat00066
,
Figure 112015084809975-pat00067
Represents the displacement vector for each x and z axis.)
Of the inverse time structure in the elastic medium.
The method according to claim 1,
An image quality enhancement unit applying an absolute value function to the structure-corrected seismic image acquired through convolution in the convolution unit;
And an inverse temporal structure compensating device in the elastic medium.
a) receiving the seismic wave exploration data and the underground structure model data;
b) calculating a virtual transmission source based on the input data and calculating a virtual sound source wave field;
c) calculating a P wave back propagation wave length by back propagation processing the seismic wave survey data;
d) separating the P-wave wave field from the virtual sound source wave field by using a relationship between stress and displacement;
e) separating the wave field from the back propagation wave field using a relationship between stress and displacement; And
f) convoluting the two wave fields separated in the steps c) and e), respectively;
And,
The relationship between the stress and the displacement,
Figure 112015084809975-pat00068

Figure 112015084809975-pat00069

(At this time
Figure 112015084809975-pat00070
,
Figure 112015084809975-pat00071
Represents a Lame constant,
Figure 112015084809975-pat00072
,
Figure 112015084809975-pat00073
Is the stress term,
Figure 112015084809975-pat00074
,
Figure 112015084809975-pat00075
Represents the displacement vector for each x and z axis.)
Of the inverse time structure in the elastic medium.
The method of claim 3,
g) applying an absolute value function to the structure-corrected seismic image acquired through convolution in step f);
Wherein the inverse temporal structure correction method further comprises:
KR1020150094460A 2015-07-02 2015-07-02 Elastic reverse-time migration system and method using absolute value function for improving the quality of subsurface structure imaging KR101564094B1 (en)

Priority Applications (2)

Application Number Priority Date Filing Date Title
KR1020150094460A KR101564094B1 (en) 2015-07-02 2015-07-02 Elastic reverse-time migration system and method using absolute value function for improving the quality of subsurface structure imaging
US15/196,841 US9798026B2 (en) 2015-07-02 2016-06-29 Elastic reverse-time migration system and method using absolute value function for improving the quality of subsurface structure imaging

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR1020150094460A KR101564094B1 (en) 2015-07-02 2015-07-02 Elastic reverse-time migration system and method using absolute value function for improving the quality of subsurface structure imaging

Publications (1)

Publication Number Publication Date
KR101564094B1 true KR101564094B1 (en) 2015-10-29

Family

ID=54430677

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020150094460A KR101564094B1 (en) 2015-07-02 2015-07-02 Elastic reverse-time migration system and method using absolute value function for improving the quality of subsurface structure imaging

Country Status (2)

Country Link
US (1) US9798026B2 (en)
KR (1) KR101564094B1 (en)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106970416A (en) * 2017-03-17 2017-07-21 中国地质科学院地球物理地球化学勘查研究所 Elastic wave least square reverse-time migration system and method based on wave field separation
KR101864314B1 (en) * 2017-11-07 2018-06-29 한국지질자원연구원 Method for signal processing of raw data acquired using chirp sub-bottom profiler and readable storage medium stroing the method
KR101903475B1 (en) 2016-02-04 2018-10-02 한국해양대학교 산학협력단 Analysis method and apparatus for rock properties with vertically transverse isotropy media
KR101907635B1 (en) 2016-05-30 2018-10-12 한양대학교 산학협력단 Seismic imaging apparatus and method using a reverse time migration
KR20210103339A (en) * 2020-02-13 2021-08-23 한양대학교 산학협력단 Image interpolation method by deep learning and device thereof
CN113406698A (en) * 2021-05-24 2021-09-17 中国石油大学(华东) Dual-phase medium elastic wave reverse time migration imaging method based on longitudinal and transverse wave decoupling

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109307890A (en) * 2017-07-28 2019-02-05 中国石油化工股份有限公司 Reverse-time migration method and system based on uplink and downlink wavefield decomposition
CN108802813B (en) * 2018-06-13 2019-07-23 中国石油大学(华东) A kind of multi-component seismic data offset imaging method and system
CN109143339B (en) * 2018-08-14 2020-06-09 中国石油天然气集团有限公司 Elastic reverse time migration imaging method and device based on transverse wave stress invariant
US11086036B2 (en) 2019-01-22 2021-08-10 Saudi Arabian Oil Company AVO imaging condition in elastic reverse time migration
CN111158047B (en) * 2020-03-04 2021-05-11 中国石油大学(北京) Three-dimensional elastic wave field vector decomposition method, device and computer storage medium
CN112327359B (en) * 2020-10-14 2022-06-14 山东省科学院海洋仪器仪表研究所 Elastic reverse time migration method based on imaging energy flow vector
CN112904426B (en) * 2021-03-27 2022-09-30 中国石油大学(华东) Decoupling elastic wave reverse time migration method, system and application

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101219746B1 (en) * 2010-08-24 2013-01-10 서울대학교산학협력단 Apparatus and method for imaging a subsurface using frequency domain reverse time migration in an elastic medium
US20130242693A1 (en) 2012-03-13 2013-09-19 Seoul National University R&Db Foundation Seismic imaging system using a reverse time migration algorithm

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101903475B1 (en) 2016-02-04 2018-10-02 한국해양대학교 산학협력단 Analysis method and apparatus for rock properties with vertically transverse isotropy media
KR101907635B1 (en) 2016-05-30 2018-10-12 한양대학교 산학협력단 Seismic imaging apparatus and method using a reverse time migration
CN106970416A (en) * 2017-03-17 2017-07-21 中国地质科学院地球物理地球化学勘查研究所 Elastic wave least square reverse-time migration system and method based on wave field separation
KR101864314B1 (en) * 2017-11-07 2018-06-29 한국지질자원연구원 Method for signal processing of raw data acquired using chirp sub-bottom profiler and readable storage medium stroing the method
US10488514B2 (en) 2017-11-07 2019-11-26 Korea Institute Of Geoscience And Mineral Resources Signal processing method of raw data obtained by using chirp sub-bottom profiler and readable storage medium storing the method
KR20210103339A (en) * 2020-02-13 2021-08-23 한양대학교 산학협력단 Image interpolation method by deep learning and device thereof
KR102324784B1 (en) 2020-02-13 2021-11-10 한양대학교 산학협력단 Image interpolation method by deep learning and device thereof
CN113406698A (en) * 2021-05-24 2021-09-17 中国石油大学(华东) Dual-phase medium elastic wave reverse time migration imaging method based on longitudinal and transverse wave decoupling
CN113406698B (en) * 2021-05-24 2023-04-21 中国石油大学(华东) Double-phase medium elastic wave reverse time migration imaging method based on longitudinal and transverse wave decoupling

Also Published As

Publication number Publication date
US20170003410A1 (en) 2017-01-05
US9798026B2 (en) 2017-10-24

Similar Documents

Publication Publication Date Title
KR101564094B1 (en) Elastic reverse-time migration system and method using absolute value function for improving the quality of subsurface structure imaging
CN108139499B (en) Q-compensated full wavefield inversion
US8437998B2 (en) Hybrid method for full waveform inversion using simultaneous and sequential source method
KR101549388B1 (en) Prestack elastic generalized-screen migration method for seismic multicomponent data
US8352190B2 (en) Method for analyzing multiple geophysical data sets
US10228476B2 (en) Method for survey data processing compensating for visco-acoustic effects in tilted transverse isotropy reverse time migration
EP3803469B1 (en) Method for generating an image of a subsurface of an area of interest from seismic data
EP2491428B1 (en) Full-waveform inversion in the traveltime domain
EP3259619B1 (en) Method and system of processing seismic data by providing surface aperture common image gathers
CN102695970B (en) An improved process for characterising the evolution of an oil or gas reservoir over time
EP3167315B1 (en) Method of constraining an inversion in the characterisation of the evolution of a subsurface volume
AU2011312806B2 (en) Hybrid method for full waveform inversion using simultaneous and sequential source method
GB2588488A (en) A full waveform inversion of seismic data using partial match filtering
EP3243089B1 (en) Method for obtaining estimates of a model parameter so as to characterise the evolution of a subsurface volume over a period of time
CN110187382A (en) A kind of diving Wave and back wave wave equation Travel Time Inversion method
Kamei et al. Application of waveform tomography to a crooked-line 2D land seismic data set
Alam Near-surface characterization using traveltime and full-waveform inversion with vertical and horizontal component seismic data
CN113376689B (en) Elastic reflection wave travel time inversion method considering interlayer multiples
US20200348434A1 (en) 4d time shift and amplitude joint inversion for velocity perturbation
US8699298B1 (en) 3D multiple prediction and removal using diplets
CN111480097A (en) Sub-salt imaging tool for interpreters
Giustiniani et al. Imaging subsurface structures using wave equation datuming advanced seismic techniques
Luo et al. Salt and subsalt structure recovery—An envelope-based waveform inversion with local angle domain illumination compensation and L-BFGS
Chen et al. Joint data and model-driven simultaneous inversion of velocity and density
CN117665965A (en) Method and system for positioning sika type iron-rich ore body based on full waveform inversion

Legal Events

Date Code Title Description
E701 Decision to grant or registration of patent right
GRNT Written decision to grant
FPAY Annual fee payment

Payment date: 20180927

Year of fee payment: 4

FPAY Annual fee payment

Payment date: 20191001

Year of fee payment: 5