CN110118993B - Diffracted wave imaging method and device - Google Patents
Diffracted wave imaging method and device Download PDFInfo
- Publication number
- CN110118993B CN110118993B CN201910443218.7A CN201910443218A CN110118993B CN 110118993 B CN110118993 B CN 110118993B CN 201910443218 A CN201910443218 A CN 201910443218A CN 110118993 B CN110118993 B CN 110118993B
- Authority
- CN
- China
- Prior art keywords
- angle
- gather
- module
- imaging
- zero
- Prior art date
- Legal status (The legal status 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 status listed.)
- Expired - Fee Related
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 86
- 238000001228 spectrum Methods 0.000 claims abstract description 34
- 238000013508 migration Methods 0.000 claims abstract description 32
- 230000005012 migration Effects 0.000 claims abstract description 31
- 230000006870 function Effects 0.000 claims description 25
- 238000000034 method Methods 0.000 claims description 17
- 230000009466 transformation Effects 0.000 claims description 11
- 238000007781 pre-processing Methods 0.000 claims description 7
- 238000004590 computer program Methods 0.000 claims description 5
- 238000010183 spectrum analysis Methods 0.000 claims description 4
- 230000001629 suppression Effects 0.000 abstract description 7
- 238000004364 calculation method Methods 0.000 abstract description 4
- 230000002159 abnormal effect Effects 0.000 description 5
- 238000012545 processing Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 238000004891 communication Methods 0.000 description 3
- 230000014509 gene expression Effects 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 239000003245 coal Substances 0.000 description 2
- 230000002238 attenuated effect Effects 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 239000002360 explosive Substances 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000002310 reflectometry Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 150000003839 salts Chemical class 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
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. for interpretation or for event detection
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention provides a diffracted wave imaging method and a diffracted wave imaging device, which relate to the technical field of seismic exploration, wherein an angle gather is extracted from seismic data, and Fourier transform is carried out on the angle gather in the angle direction to obtain an amplitude spectrum and a phase spectrum; enabling the phase spectrum to be zero, and then carrying out Fourier inversion, so that the angle gather is migrated according to the dip angle by taking the true dip angle migration of the reflected wave to the zero dip angle as a reference, and the migrated angle gather is obtained, so that the energy of the reflected wave is distributed near the zero dip angle, and the energy of the diffracted wave is distributed in a wide angle range; designing a truncation function, and setting the energy near the zero dip angle to zero, so that reflected waves are effectively suppressed, and diffracted waves are obtained; finally, performing superposition imaging to obtain a diffracted wave imaging result; the invention avoids the calculation of the dip angle field, realizes the effective suppression of the reflected wave, and has high imaging quality, simplicity and high efficiency.
Description
Technical Field
The invention relates to the technical field of seismic exploration, in particular to a diffracted wave imaging method and a diffracted wave imaging device.
Background
Correctly identifying geological discontinuities, such as faults, pinch-outs, and small-sized scatterers, is a significant challenge in seismic exploration. The seismic response of these objects is present in diffracted waves. Seismic diffraction is an ideal carrier of information for these small scale discrete objects and can therefore be used to detect these geological objects. However, it is difficult to identify because diffracted energy is weak relative to specular reflected energy. The positioning of scatterers such as faults, pinch-out, sharp changes in reflectivity, salt side, etc. is a significant challenge in seismic exploration. The wavefields produced by these objects are characterized by the presence of diffracted energy. Conventional seismic processing tends to enhance reflections, treating other waves, including diffracted waves, as noise. Therefore, the diffracted waves are attenuated during conventional processing, and it is necessary to provide an effective imaging method for the diffracted wave characteristics.
The angle gather is an imaging gather in which the energy of reflected waves is concentrated near the true dip and diffracted waves are dispersed over a wide angle range. The conventional reflected wave suppression method is to remove the energy of the reflected wave on the premise of knowing the inclination angle and then realize diffracted wave imaging. However, this method requires the calculation of a complex tilt angle field, and the accuracy error of the tilt angle field will affect the imaging quality of the diffraction body.
Disclosure of Invention
The invention aims to provide a diffracted wave imaging method and a diffracted wave imaging device, which do not need to calculate an inclination angle field, avoid the error of imaging results of a diffraction body caused by the error of the inclination angle field, have high imaging quality, can realize the effective suppression of reflected waves, and have simple and efficient method.
In a first aspect, the present invention provides a diffracted wave imaging method, including:
determining an angle gather from the seismic data, the angle gather comprising reflected wave energy and diffracted wave energy;
taking the migration of the true inclination angle of the reflected wave to a zero inclination angle as a reference, and migrating the angle gather according to the inclination angle to obtain a migrated angle gather;
removing energy near the zero dip angle in the corner gather after the migration to obtain the corner gather after the removal;
and carrying out superposition imaging on the removed angle gather to obtain a diffracted wave imaging result.
In a second aspect, the invention provides a diffracted wave imaging apparatus, which includes a preprocessing module, a migration module, a removal module, and an imaging module;
the preprocessing module is used for determining an angle gather according to the seismic data, wherein the angle gather comprises reflected wave energy and diffracted wave energy;
the migration module is used for migrating the angle gather according to the inclination angle by taking the true inclination angle of the reflected wave as a reference to be migrated to a zero inclination angle, so as to obtain the migrated angle gather;
the removing module is used for removing energy near a zero dip angle in the corner gather after the migration to obtain the corner gather after the removal;
and the imaging module is used for carrying out superposition imaging on the removed angle gather to obtain a diffracted wave imaging result.
In a third aspect, the present invention provides an electronic device, which includes a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor implements the steps of the method of the first aspect when executing the computer program.
In a fourth aspect, the present invention provides a computer readable medium having non-volatile program code executable by a processor, the program code causing the processor to perform the method of the first aspect.
According to the diffracted wave imaging method and device, angle gather data are extracted from seismic data, reflected wave energy in the angle gather data is transferred to a zero inclination angle from a true inclination angle position, reflected wave energy is removed, diffracted wave angle gather data are obtained, and diffracted wave imaging is finally carried out, so that the position of a geological abnormal body is determined according to a diffracted wave imaging result; the invention can realize effective suppression of reflected waves without calculating an inclination angle field, improves the imaging quality and has simple and efficient method.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, and it is obvious that the drawings in the following description are some embodiments of the present invention, and other drawings can be obtained by those skilled in the art without creative efforts.
FIG. 1 is a flowchart of a diffracted wave imaging method according to an embodiment of the present invention;
FIG. 2 is a schematic diagram of an angle gather of a diffracted wave imaging method according to an embodiment of the present invention;
FIG. 3 is a schematic diagram of a diffracted wave imaging apparatus according to an embodiment of the present invention;
fig. 4 is a schematic diagram of an electronic device according to an embodiment of the present invention.
Icon: 301-a pre-processing module; 302-a migration module; 303-removal module; 304-an imaging module; 400-an electronic device; 401 — a communication interface; 402-a processor; 403-a memory; 404-bus.
Detailed Description
The technical solutions of the present invention will be described clearly and completely with reference to the following embodiments, and it should be understood that the described embodiments are some, but not all, embodiments of the present invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
At present, the existing diffracted wave imaging method mainly has two ideas, namely, the diffracted wave imaging is directly carried out by utilizing the kinematics and dynamics characteristics (such as hyperbolic characteristic and amplitude attenuation rule of a travel curve) of diffracted waves and reflected waves; secondly, aiming at the difference between the diffracted wave and the reflected wave, the diffracted wave field is firstly separated, and then imaging is carried out. The plane wave decomposition method is a commonly used diffraction wave imaging method, and the suppression is performed on reflected waves under the condition that the inclination angle of the reflected waves is known, so that a complex inclination angle field needs to be calculated, and the precision error of the inclination angle field can influence the imaging quality of a diffraction body. Based on the method and the device, the diffracted wave imaging method and the device do not need to calculate an inclination angle field, can effectively suppress reflected waves and obtain a high-quality diffracted wave imaging result.
Referring to fig. 1, a diffracted wave imaging method includes:
s101, determining an angle gather according to seismic data, wherein the angle gather comprises reflection wave energy and diffracted wave energy;
specifically, the seismic data are acquired by a seismic prospecting instrument, the seismic prospecting instrument comprises a seismic source generating device, a seismic wave receiving device and a recording device, the seismic wave receiving device is mainly a seismic detector, and the recording device mainly adopts a seismic information acquisition system. The seismic source is mainly generated by explosive charge discharge, manual swinging of a sledge hammer or a seismic source instrument and the like.
Seismic data are obtained through a seismic prospecting instrument, migration processing is carried out on the seismic data, and an angle gather is obtained. Commonly used migration methods are Kirchhoff integration, finite difference method, and frequency-wavenumber domain migration.
An angular gather is a gather of imaging points, the full name of which is an angular domain common imaging points gather (ADCIGs).
Seismic waves are vibrations that propagate around a seismic source, and refer to elastic waves that radiate from the source to the surroundings.
The reflected wave is formed under conditions in which the subterranean formation presents a wave impedance interface. When seismic waves pass through discontinuous elastic discontinuities such as discontinuities of strata, pinch points or unconformity contact points of strata, corner points of faults and the like, as long as the sizes of geologic bodies are approximately equal to the wavelength of the seismic waves, the discontinuous discontinuities can be regarded as a new seismic source, the new seismic source generates a new disturbance and transmits the new disturbance to the periphery of an elastic space, and the waves are called diffracted waves in seismic exploration. The reflected wave energy is related to the reflected wave amplitude, and similarly, the diffracted wave energy is related to the diffracted wave amplitude.
It should be noted that the seismic waves include other waves besides the reflected waves and the diffracted waves, and only the reflected waves and the diffracted waves are referred to in this embodiment, and other waveforms are not described in detail.
And S102, taking the true inclination angle of the reflected wave as a reference, and migrating the angle gather according to the inclination angle to obtain the migrated angle gather.
In particular, the true dip is the maximum angle of the horizontal plane to the formation level, i.e. the true dip is the angle between the inclined plane and a horizontal reference plane measured on a cross section running perpendicular to the inclined plane.
Referring to FIG. 2, a schematic diagram of an angle gather is shown, where curve 201 is a diffracted wave and curves 202 and 203 are reflected waves. As can be seen from fig. 2, the images of the diffracted waves with different inclinations are aligned, the curve of the reflected wave is represented by a pseudo-hyperbolic curve, and the inclination corresponding to the vertex of the reflected wave is the true inclination of the reflecting interface. That is, the reflective interface corresponding to curve 202 is horizontal and the reflective interface corresponding to curve 203 is inclined.
And migrating the angle gather by taking the migration of the true dip angle of the reflected wave to the zero dip angle as a reference, wherein after migration, the diffracted wave energy is distributed in a wide angle range, and the reflected wave energy is distributed near the zero dip angle.
S103, removing energy near the zero dip angle in the corner gather after the migration to obtain the removed corner gather.
Specifically, a cut-off function is designed by utilizing the Fresnel zone range under the angle domain, and energy near a zero dip angle is removed to remove reflected wave energy, so that the reflected wave is effectively suppressed.
And S104, carrying out superposition imaging on the removed angle gather to obtain a diffracted wave imaging result.
Specifically, the angle gathers only protecting the diffracted wave energy are obtained in step S103, and are subjected to superposition imaging, so as to obtain the diffracted wave imaging result.
In the embodiment, address data is subjected to offset processing to extract an angle gather, the angle gather is migrated to migrate reflected wave energy to be near a zero dip angle, then the energy near the zero dip angle is set to be zero by setting a truncation function, so that the reflected wave energy is suppressed, the angle gather only containing diffracted wave energy is obtained, and the angle gather is subjected to superposition imaging, so that a diffracted wave imaging result is obtained; according to the method, the calculation of the inclination angle field is avoided, the effective suppression of the reflected wave is realized, the imaging quality is effectively improved, and the method is simple and efficient.
Optionally, the step S102 includes:
carrying out Fourier transform on the angle gather to obtain first time-frequency domain data;
obtaining a corresponding amplitude spectrum and a corresponding phase spectrum according to the first time-frequency domain data;
setting the phase spectrum to be zero, and obtaining second time-frequency domain data according to the amplitude spectrum and the phase spectrum;
and performing Fourier inverse transformation on the second time-frequency domain data to obtain the corner gather after the migration.
Specifically, the relationship between the frequency of the seismic signal and the amplitude and phase may be referred to as an amplitude spectrum and a phase spectrum, and in fourier transform, the amplitude of each component becomes the amplitude spectrum of the signal with the change in frequency, and the phase of each component with the change in angular frequency is referred to as the phase spectrum of the signal.
Setting the angle gather subjected to the offset processing as theta (x, y, t, α), wherein (x, y) represents geological imaging point coordinates, t represents imaging point corresponding time, and α represents an inclination angle;
fourier transforming the angle gather θ (x, y, t, α) according to equation (1):
θ(x,y,t,f)＝∫θ(x,y,t,α)e-jfαdα (1)，
obtaining first time-frequency domain data theta (x, y, t, f), wherein f represents the frequency of seismic data, and j represents an imaginary unit;
obtaining an amplitude map and a phase spectrum according to the equations (2) and (3):
wherein A (x, y, t, f) is an amplitude spectrum, P (x, y, t, f) is a phase spectrum, real represents a real part of an imaginary number, and imag represents the imaginary part;
setting the phase spectrum of the expression (3) to zero, and obtaining second time-frequency domain data theta according to the expressions (4) and (5)p(x,y,t,f)，
real(θp(x,y,t,f))＝A(x,y,t,f) (4)，
imag(θp(x,y,t,f))＝A(x,y,t,f) (5)，
According to equation (6), time-frequency domain data thetap(x, y, t, f) carrying out inverse Fourier transform to obtain the angle gather data theta after the migrationp(x,y,t,α)，
In the embodiment, the angle gather is subjected to fourier transform, that is, fourier transform is performed according to the angle direction to obtain an amplitude spectrum and a phase spectrum, the phase spectrum is made to be zero, second time-frequency domain data is obtained through expressions (2) and (3), and fourier inverse transform is performed on the second time-frequency domain data to obtain the transferred angle gather, so that the reflected wave is transferred from a true dip angle position to a zero dip angle position, the energy of the transferred reflected wave is mainly concentrated near the zero dip angle, and the energy of the diffracted wave is in a wide angle range; therefore, the reflected wave energy is transferred from the true inclination angle position to the zero inclination angle position through a time shifting method, and preparation is made for enabling the energy near the zero inclination angle to be returned to zero through a stage function in the next step.
Optionally, the step S103 includes:
setting a truncation function with a zero dip angle as a center;
and multiplying the truncation function by the angle gather after migration to obtain the angle gather after the energy near the zero dip angle is removed.
Optionally, the truncation function in step S103 is:
wherein w (x, y, t, α) represents a truncation function, coordinates of (x, y) tabular geological imaging point, t represents corresponding time of the imaging point, α represents an inclination angle, fmRepresenting the seismic data dominant frequency.
Specifically, the angle gathers containing only diffracted waves are obtained according to the following equation:
θpw(x,y,t,α)＝w(x,y,t,α)×θp(x,y,t,α) (8)，
wherein, thetap(x, y, t, α) represents the corner gather after migration, θpw(x, y, t, α) represents an angle gather containing only diffracted waves.
And (3) performing superposition imaging according to the formula (8) to obtain a diffracted wave imaging result I (x, y, t):
I(x,y,t)＝∫θpw(x,y,t,α)dα (9)，
after the reflected wave energy is concentrated near the zero inclination angle through the previous step, the diffracted wave energy is distributed in a wide angle range, and thus the zero inclination angle is taken as the center; the step sets a truncation function to zero the energy near the zero inclination angle, and the reflected wave energy is concentrated and the diffracted wave energy is distributed widely, so that the reflected wave can be removed and the diffracted wave can be reserved.
According to the diffracted wave imaging method provided by the embodiment, the time shifting method is utilized to shift the energy of the reflected wave from the true dip angle position to the zero dip angle position, after the shift, the energy of the diffracted wave is distributed in a wide angle range, and the energy of the reflected wave is distributed near the zero dip angle; and then designing a truncation function by utilizing the Fresnel zone range under the angular domain to remove the energy of the reflected wave. The method avoids the calculation of an inclination field, can realize the effective suppression of reflected waves, obtains a diffracted wave imaging result, improves the imaging quality, and is simple and efficient. Accurate identification of geological abnormal bodies (such as breakpoints, stratum sharp vanishing points and the like) can be realized by utilizing the diffracted wave imaging result, and the geological abnormal bodies are closely related to oil and gas migration and coal mining safety. The embodiment can be used for detecting complex structures, can accurately position the positions of the diffraction body-geological abnormal bodies, such as fault break points, coal seam point vanishing points and the like, and realizes high-resolution imaging of the underground geological abnormal bodies.
The embodiment of the invention also provides a diffracted wave imaging device, which refers to fig. 3 and comprises a preprocessing module 301, a migration module 302, a removal module 303 and an imaging module 304;
the preprocessing module 301 is configured to determine an angle gather according to the seismic data, where the angle gather includes reflected wave energy and diffracted wave energy;
the migration module 302 is configured to migrate the angle gather according to the inclination angle based on migration of the true inclination angle of the reflected wave to a zero inclination angle, so as to obtain a migrated angle gather;
the removing module 303 is configured to remove energy near a zero dip in the migrated angle gather to obtain a removed angle gather;
the imaging module 304 is configured to perform superposition imaging on the removed angle gathers to obtain a diffracted wave imaging result.
Optionally, the migration module 302 includes a first transformation module, a spectrum analysis module, a zeroing module, and a second transformation module;
the first transformation module is used for carrying out Fourier transformation on the angle gather to obtain first time-frequency domain data;
the spectrum analysis module is used for obtaining a corresponding amplitude spectrum and a corresponding phase spectrum according to the first time-frequency domain data;
the zero setting module is used for setting the phase spectrum to be zero and obtaining second time-frequency domain data according to the amplitude spectrum and the phase spectrum;
and the second transformation module is used for carrying out Fourier inversion on the second time-frequency domain data to obtain the corner gather after the migration.
Optionally, the removal module 303 comprises a truncation function module and an energy removal module;
the truncation function module is used for setting a truncation function by taking a zero inclination angle as a center;
and the energy removing module is used for multiplying the truncation function by the angle gather after migration to obtain the angle gather after the energy near the zero dip angle is removed.
Optionally, the truncation function is:
wherein w (x, y, t, α) represents a truncation function, coordinates of (x, y) tabular geological imaging points, t represents corresponding time of the imaging points, α represents an inclination angle, and fm represents a seismic data main frequency.
Referring to fig. 4, an embodiment of the present invention further provides an apparatus, and an embodiment of the present invention further provides an electronic apparatus 400, which includes a communication interface 401, a processor 402, a memory 403, and a bus 404, where the processor 402, the communication interface 401, and the memory 403 are connected by the bus 404; the memory 403 is used for storing computer programs that support the processor 402 to execute the diffracted wave imaging method, and the processor 402 is configured to execute the programs stored in the memory 403.
Optionally, an embodiment of the present invention further provides a computer readable medium having a non-volatile program code executable by a processor, where the program code causes the processor to execute the diffracted wave imaging method in the above embodiment.
Finally, it should be noted that: the above embodiments are only used to illustrate the technical solution of the present invention, and not to limit the same; while the invention has been described in detail and with reference to the foregoing embodiments, it will be understood by those skilled in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some or all of the technical features may be equivalently replaced; and the modifications or the substitutions do not make the essence of the corresponding technical solutions depart from the scope of the technical solutions of the embodiments of the present invention.
Claims (8)
1. A diffracted wave imaging method, comprising:
determining an angle gather from the seismic data, the angle gather comprising reflected wave energy and diffracted wave energy;
carrying out Fourier transform on the angle gather to obtain first time-frequency domain data;
obtaining a corresponding amplitude spectrum and a corresponding phase spectrum according to the first time-frequency domain data;
setting the phase spectrum to be zero, and obtaining second time-frequency domain data according to the amplitude spectrum and the phase spectrum;
performing Fourier inversion on the second time-frequency domain data to obtain an angle gather after migration;
removing energy near the zero dip angle in the corner gather after the migration to obtain the corner gather after the removal;
and carrying out superposition imaging on the removed angle gather to obtain a diffracted wave imaging result.
2. The diffracted wave imaging method of claim 1, wherein removing energy in the vicinity of zero dip in the migrated angular gathers, resulting in removed angular gathers comprising:
setting a truncation function with a zero dip angle as a center;
and multiplying the truncation function by the angle gather after migration to obtain the angle gather after the energy near the zero dip angle is removed.
3. The diffracted wave imaging method of claim 2, wherein the truncation function is:
wherein w (x, y, t, α) represents a truncation function, coordinates of (x, y) tabular geological imaging point, t represents corresponding time of the imaging point, α represents an inclination angle, fmRepresenting the seismic data dominant frequency.
4. A diffraction wave imaging device is characterized by comprising a preprocessing module, a migration module, a removal module and an imaging module;
the preprocessing module is used for determining an angle gather according to the seismic data, wherein the angle gather comprises reflected wave energy and diffracted wave energy;
the migration module comprises a first transformation module, a spectrum analysis module, a zero setting module and a second transformation module;
the first transformation module is used for carrying out Fourier transformation on the angle gather to obtain first time-frequency domain data;
the spectrum analysis module is used for obtaining a corresponding amplitude spectrum and a corresponding phase spectrum according to the first time-frequency domain data;
the zero setting module is used for setting the phase spectrum to be zero and obtaining second time-frequency domain data according to the amplitude spectrum and the phase spectrum;
the second transformation module is used for carrying out Fourier inversion on the second time-frequency domain data to obtain an angle gather after the migration;
the removing module is used for removing energy near a zero dip angle in the corner gather after the migration to obtain the corner gather after the removal;
and the imaging module is used for carrying out superposition imaging on the removed angle gather to obtain a diffracted wave imaging result.
5. The diffracted wave imaging apparatus of claim 4, wherein the removal module comprises a truncation function module and an energy removal module;
the truncation function module is used for setting a truncation function by taking a zero inclination angle as a center;
and the energy removing module is used for multiplying the truncation function and the angle gather after migration to obtain the angle gather after the energy near the zero dip angle is removed.
6. The diffracted wave imaging apparatus of claim 5, wherein the truncation function is:
wherein w (x, y, t, α) represents a truncation function, coordinates of (x, y) tabular geological imaging point, t represents corresponding time of the imaging point, α represents an inclination angle, fmRepresenting the seismic data dominant frequency.
7. An electronic device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the steps of the method of any of the preceding claims 1 to 3 are implemented when the computer program is executed by the processor.
8. A computer-readable medium having non-volatile program code executable by a processor, wherein the program code causes the processor to perform the method of any of claims 1 to 3.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910443218.7A CN110118993B (en) | 2019-05-22 | 2019-05-22 | Diffracted wave imaging method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910443218.7A CN110118993B (en) | 2019-05-22 | 2019-05-22 | Diffracted wave imaging method and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110118993A CN110118993A (en) | 2019-08-13 |
| CN110118993B true CN110118993B (en) | 2020-05-15 |
Family
ID=67523258
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910443218.7A Expired - Fee Related CN110118993B (en) | 2019-05-22 | 2019-05-22 | Diffracted wave imaging method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110118993B (en) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114114420B (en) * | 2020-09-01 | 2024-02-23 | 中国石油化工股份有限公司 | Diffraction identification imaging method, diffraction identification imaging device, electronic equipment and medium |
| CN114910951A (en) * | 2021-02-10 | 2022-08-16 | 中国石油天然气股份有限公司 | Diffracted wave imaging method and device |
| CN114910952A (en) * | 2021-02-10 | 2022-08-16 | 中国石油天然气股份有限公司 | Scattered wave imaging method and device |
| CN113608261B (en) * | 2021-07-30 | 2021-12-28 | 中国矿业大学(北京) | Diffracted wave imaging method and device and electronic equipment |
| CN113687416B (en) * | 2021-08-19 | 2022-08-30 | 吉林大学 | Imaging method for inclined metal ore body |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102455439A (en) * | 2010-11-02 | 2012-05-16 | 中国石油大学(北京) | Diffracted wave field separation method based on kirchhoff integration method |
| CN103984012A (en) * | 2014-04-16 | 2014-08-13 | 孙赞东 | Diffracted wave field separation method based on pre-stack gaussian beam depth migration |
| CN109143343A (en) * | 2018-11-14 | 2019-01-04 | 中国矿业大学(北京) | Diffraction wave imaging method, device and electronic equipment |
| CN110031900A (en) * | 2019-05-27 | 2019-07-19 | 中国矿业大学(北京) | Diffraction wave imaging method, device, electronic equipment and storage medium |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2673662B1 (en) * | 2011-02-08 | 2014-12-17 | Total SA | Method of analyzing seismic data |
-
2019
- 2019-05-22 CN CN201910443218.7A patent/CN110118993B/en not_active Expired - Fee Related
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102455439A (en) * | 2010-11-02 | 2012-05-16 | 中国石油大学(北京) | Diffracted wave field separation method based on kirchhoff integration method |
| CN103984012A (en) * | 2014-04-16 | 2014-08-13 | 孙赞东 | Diffracted wave field separation method based on pre-stack gaussian beam depth migration |
| CN109143343A (en) * | 2018-11-14 | 2019-01-04 | 中国矿业大学(北京) | Diffraction wave imaging method, device and electronic equipment |
| CN110031900A (en) * | 2019-05-27 | 2019-07-19 | 中国矿业大学(北京) | Diffraction wave imaging method, device, electronic equipment and storage medium |
Non-Patent Citations (1)
| Title |
|---|
| 基于Kirchhoff积分的绕射波分离与成像方法研究;郭威;《中国优秀硕士学位论文全文数据库 基础科学辑》;20150615(第6期);第3章 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110118993A (en) | 2019-08-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110118993B (en) | Diffracted wave imaging method and device | |
| US8451682B2 (en) | Method and apparatus for deghosting seismic data | |
| US9250341B2 (en) | Device and method for extrapolating specular energy of reverse time migration three dimensional angle gathers | |
| EP3749987B1 (en) | Systems and methods to enhance 3-d prestack seismic data based on non-linear beamforming in the cross-spread domain | |
| US8437218B2 (en) | Correcting geometry-related time and amplitude errors | |
| EA032186B1 (en) | Seismic adaptive focusing | |
| EA023300B1 (en) | Method for dynamic aperture determination for three-dimensional surface-related multiple elimination | |
| US20150362608A1 (en) | Combined interpolation and primary estimation | |
| CN113640872A (en) | Diffracted wave separation method and device and electronic equipment | |
| US11125900B2 (en) | Estimating an earth response | |
| US11662492B2 (en) | Seismic random noise attenuation | |
| WO2014132023A2 (en) | Enhanced visualisation of geologic features in 3d seismic survey data | |
| CN106610507B (en) | A kind of method and apparatus of removal ghosting interference | |
| Li et al. | A causal imaging condition for reverse time migration using the Discrete Hilbert transform and its efficient implementation on GPU | |
| CN111474580A (en) | Azimuth angle gather extraction method and system based on offset vector piece | |
| CN111691876B (en) | Method, device and storage medium for imaging adjacent well by using acoustic logging | |
| CN111983682B (en) | Seismic profile imaging method and device and electronic equipment | |
| CN108802822B (en) | The direct prestack time migration method of guarantor's width and device in direction anisotropy medium | |
| US9052408B2 (en) | Method and system for detecting geometry variations of seismic reflections | |
| US10739483B2 (en) | Absolute strength and absolute sensitivity in seismic data | |
| US10067254B2 (en) | Removal of an estimated acquisition effect from a marine survey measurement | |
| Tan et al. | Combined adaptive multiple subtraction based on optimized event tracing and extended wiener filtering | |
| CN116774293B (en) | Method, system, electronic equipment and medium for automatically picking up same phase shaft | |
| CN115184987B (en) | Stratum inclination angle information extraction method and device, server and storage medium | |
| CN112925025B (en) | Seismic data static correction method and device based on surface waves and electronic equipment |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200515 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |