CN112904426A - Decoupling elastic wave reverse time migration method, system and application - Google Patents
Decoupling elastic wave reverse time migration method, system and application Download PDFInfo
- Publication number
- CN112904426A CN112904426A CN202110329092.8A CN202110329092A CN112904426A CN 112904426 A CN112904426 A CN 112904426A CN 202110329092 A CN202110329092 A CN 202110329092A CN 112904426 A CN112904426 A CN 112904426A
- Authority
- CN
- China
- Prior art keywords
- wave
- imaging
- elastic
- vector
- longitudinal
- 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.)
- Granted
Links
- 230000005012 migration Effects 0.000 title claims abstract description 68
- 238000013508 migration Methods 0.000 title claims abstract description 68
- 238000000034 method Methods 0.000 title claims abstract description 46
- 238000003384 imaging method Methods 0.000 claims abstract description 142
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 34
- 238000001514 detection method Methods 0.000 claims abstract description 34
- 238000000926 separation method Methods 0.000 claims abstract description 32
- 230000006870 function Effects 0.000 claims description 23
- 238000012633 nuclear imaging Methods 0.000 claims description 14
- 239000011159 matrix material Substances 0.000 claims description 13
- 238000004364 calculation method Methods 0.000 claims description 8
- 238000005286 illumination Methods 0.000 claims description 6
- 238000004590 computer program Methods 0.000 claims description 5
- 230000005284 excitation Effects 0.000 claims description 3
- 239000004576 sand Substances 0.000 claims description 3
- 238000010168 coupling process Methods 0.000 claims 2
- 238000005859 coupling reaction Methods 0.000 claims 2
- 230000008878 coupling Effects 0.000 claims 1
- 230000017105 transposition Effects 0.000 claims 1
- 239000000523 sample Substances 0.000 description 10
- 238000010586 diagram Methods 0.000 description 9
- 238000004422 calculation algorithm Methods 0.000 description 3
- 230000007547 defect Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 230000010287 polarization Effects 0.000 description 2
- 150000003839 salts Chemical class 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
- 238000004321 preservation Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 230000001568 sexual effect Effects 0.000 description 1
- 239000007787 solid Substances 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
-
- 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
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 belongs to the technical field of exploration seismology data processing, and discloses a decoupling elastic wave reverse time migration method, a system and application, wherein an input seismic record, seismic wavelets, a density model and a velocity model are obtained; obtaining a seismic source and wave detection point continuation wave field by solving an elastic medium wave equation; elastic wave field separation is carried out by using vector Helmholtz decomposition, and separated amplitude-preserving pure longitudinal wave and transverse wave vector wave fields are obtained; calculating imaging conditions by using an elastic wave impedance sensitive kernel function to obtain high signal-to-noise ratio and high precision longitudinal wave-longitudinal wave (PP) and longitudinal wave-transverse wave (PS) imaging results; and superposing all shot integrated image results to obtain a final offset imaging section. The invention can obtain the separation result of pure longitudinal wave and pure transverse wave vector wave fields with true amplitude, improve the fidelity of imaging amplitude, automatically avoid low-frequency noise generated by cross-correlation of double-pass wave fields, improve the signal-to-noise ratio and resolution of imaging, and assist in high-precision longitudinal and transverse wave combined seismic interpretation.
Description
Technical Field
The invention belongs to the technical field of exploration seismology data processing, and particularly relates to a decoupling elastic wave reverse time migration method, a system and application.
Background
At present, in the processing of multi-wave multi-component exploration seismic data, elastic wave reverse time migration is a migration method with the highest imaging precision at present. However, because the double-pass wave information is used for wave field extension, strong low-frequency noise and migration false images are often generated in complex areas such as salt domes, fault blocks, retrograde faults and the like, the spatial resolution and the imaging signal-to-noise ratio of migration results are reduced, and the accuracy of subsequent seismic interpretation is seriously influenced. In addition, the traditional elastic wave imaging method does not usually adopt longitudinal and transverse wave field separation, and the extended transverse and vertical wave field components are directly used for cross-correlation imaging, so that the physical meaning of an imaging result is not clear, and serious longitudinal and transverse wave crosstalk noise is easily caused.
Through the above analysis, the problems and defects of the prior art are as follows:
(1) in the conventional multiwave multicomponent elastic wave reverse-time migration method, because the wave field continuation is carried out by solving a two-way wave equation, strong low-frequency noise and migration false images are often generated in complex areas such as salt domes, fault blocks, reverse-impact faults and the like, the imaging resolution and the signal-to-noise ratio are reduced, and the precision of subsequent seismic interpretation is seriously influenced.
(2) The traditional elastic wave imaging method does not usually carry out longitudinal and transverse wave field separation, and directly uses the extended transverse and vertical elastic wave field components to carry out cross-correlation imaging, so that the physical meaning of an imaging result is not clear, and serious longitudinal and transverse wave crosstalk noise is easily caused.
The difficulty in solving the above problems and defects is:
(1) for low-frequency noise occurring in elastic wave reverse time migration, a conventional processing mode is to firstly use elastic wave reverse time migration to obtain an imaging result with low-frequency noise, and then use a post-stack laplacian high-pass filter to filter the low-frequency noise, and the processing flow easily causes phase distortion of the filtered imaging result and inaccurate imaging amplitude. And imaging processing methods to remove these low frequency noise are less studied directly in wavefield continuation and application imaging conditions.
(2) In order to reduce the cross-wave noise in elastic wave imaging, the extended wave fields of the seismic source and the wave detection point need to be subjected to longitudinal and transverse wave separation before the imaging conditions are applied, and although pure longitudinal wave and pure transverse wave fields can be extracted by using simple Helmholtz decomposition, the phase and amplitude of the separated wave fields are inconsistent with those of the original wave fields, so that the final imaging amplitude is inaccurate.
The significance of solving the problems and the defects is as follows:
(1) the invention aims at the problem of elastic wave reverse time migration low-frequency noise, a brand-new impedance sensitive kernel function is constructed through linear combination of wave fields of a seismic source and a wave detection point, and the sensitive kernel function is used for calculating a zero-delay cross-correlation imaging condition so as to achieve the effect of directly removing the low-frequency noise in the imaging process.
(2) In addition, the invention provides a novel true amplitude vector wave Helmholtz wave field decomposition method, compared with the conventional simple scalar Helmholtz decomposition, the method firstly obtains an auxiliary wave field by solving a Poisson equation, and then obtains separated true amplitude vector pure longitudinal wave and pure transverse wave fields by applying divergence and rotation operations to the auxiliary wave field twice.
Disclosure of Invention
Aiming at the problems in the prior art, the invention provides a decoupling elastic wave reverse time migration theoretical method, a processing algorithm and a technical process, in particular relates to a decoupling elastic wave reverse time migration theoretical method, an efficient algorithm and a processing process based on vector Helmholtz decomposition and impedance sensitive nuclear imaging conditions, and aims to solve the problems of low-frequency noise and cross wave crosstalk artifact in elastic wave reverse time migration.
The invention is realized in such a way that a decoupling elastic wave reverse time migration method comprises the following steps: obtaining input seismic records, seismic wavelets, density and velocity models; obtaining a seismic source and wave detection point continuation wave field by solving an elastic wave fluctuation equation; performing wave field separation by using vector Helmholtz decomposition to obtain separated pure longitudinal (P) wave and pure transverse (S) wave vector wave fields with true amplitude; calculating imaging conditions by using an elastic wave impedance sensitive kernel function to obtain high-precision PP and PS imaging results; and superposing all shot integrated image results to obtain a final offset imaging section.
Further, the decoupling elastic wave reverse time migration method comprises the following steps:
step one, acquiring input data; wherein the input data comprises: longitudinal wave offset velocity model vp(x) Transverse wave migration velocity model vs(x) Density model rho (x), source function wavelet f (t), and multi-wave multi-component observation data d (x)rT), providing a model and seismic data basis for subsequent wavefield continuation;
secondly, according to the input seismic source wavelet f (t), the density model rho (x) and the longitudinal wave velocity model vp(x) And transverse wave velocity model vs(x) Calculating the extended vector wave field u at one side of the seismic source by solving the first-order velocity-stress elastic wave equations(x, t) providing a required source wavefield for subsequent wavefield separation and application of imaging conditions;
thirdly, according to the input density model rho (x) and the longitudinal wave velocity model vp(x) Transverse wave velocity model vs(x) And observed data d (x)rT), calculating the extended wave field u of the wave detection point vector by solving the elastic wave equationr(x, t) providing a required demodulator point wavefield for subsequent wavefield separation and application of imaging conditions;
step four, adopting vector Helmholtz decomposition to carry out longitudinal and transverse wave separation on extended seismic source and wave detection point wave fields, and extracting vector pure longitudinal wave and pure transverse wave fields, wherein the innovative processing step can eliminate crosstalk noise of longitudinal waves and transverse waves in an imaging result, improve the imaging quality and enable the imaging result to have clear physical significance;
step five, after the separated seismic source and wave detection point wave fields are obtained, cross-correlation imaging is carried out by using elastic impedance sensitive nuclear imaging conditions to obtain PP and PS imaging results, and the innovative processing step can automatically remove low-frequency noise in the imaging results when the imaging conditions are applied, and improve the signal-to-noise ratio and the resolution of imaging;
and step six, adding the imaging results of all shot gathers to obtain a final offset imaging section, and improving the balance of deep imaging amplitude by using seismic source illumination.
Further, in the second step, the first-order velocity-stress elasticity fluctuation equation is in the form of:
where x is the imaging spatial position, xsFor the seismic source location, t is the wavefield travel time, us＝[us,x us,y us,z]TFor the seismic source polarization velocity wave field, σs＝[σs,xx σs,yy σs,zz σs,xy σs,xz σs,yz]TFor the source stress wavefield, T is the transposed symbol and δ is the kronecker function. L is a partial derivative matrix, C (x) is a parameter matrix, having the following expression:
further, in step three, the expression of the elastic wave equation is as follows:
wherein x isrFor the position of the detector point, L is a partial derivative matrix, C (x) is a parameter matrix, and the expression is:
further, in step four, the longitudinal and transverse wave separation is performed on the wave fields of the source and the demodulator probe by using vector helmholtz decomposition, and the following expression is provided:
wherein,
in order to perform the gradient operation, the method comprises the following steps,
in order to calculate the divergence, the method comprises the following steps,
in order to calculate the rotation degree,
for the separated longitudinal wave source wavefield,
is the shear wave source wavefield and,
is a longitudinal wave demodulator probe wave field,
is a shear wave demodulator probe wavefield. w is asAnd wrTo assist the vector wavefield, it can be obtained by solving the following poisson equation:
where Δ is the laplace operator.
Further, in the fifth step, after the separated seismic source and geophone point wave fields are obtained, cross-correlation imaging is performed by using the elastic impedance sensitive nuclear imaging condition to obtain PP and PS imaging results, which have the following expression:
wherein,
in order to be able to obtain a shear modulus,
in order to be the bulk modulus,
and I is an identity matrix.
Another object of the present invention is to provide a decoupling elastic wave reverse time migration system using the decoupling elastic wave reverse time migration method, wherein the decoupling elastic wave reverse time migration system includes:
the input data acquisition module is used for acquiring input data; wherein the input data comprises: longitudinal wave offset velocity model vp(x) Transverse wave migration velocity model vs(x) Density model rho (x), source function wavelet f (t), and multi-wave multi-component observation data d (x)r,t)；
A seismic source vector continuation wave field calculation module for calculating the vector continuation wave field of the seismic source according to the input seismic source wavelets f (t), the density model rho (x) and the longitudinal wave velocity model vp(x) And transverse wave velocity model vs(x) Calculating the vector continuation wave field u of the seismic source by solving a first-order velocity-stress elastic wave equations(x,t)；
A wave detection point vector continuation wave field calculation module for calculating the wave field according to the input density model rho (x) and the longitudinal wave velocity model vp(x) Transverse wave velocity model vs(x) And observed data d (x)rT) by solving for the bulletCalculating a wave detection point vector continuation wave field u by using a sexual wave equationr(x,t)；
The longitudinal and transverse wave separation module is used for performing longitudinal and transverse wave separation on wave fields of the seismic source and the wave detection point by adopting vector Helmholtz decomposition;
the imaging result acquisition module is used for performing cross-correlation imaging by using elastic impedance sensitive nuclear imaging conditions after acquiring the separated seismic source and wave detection point wave fields to acquire PP and PS imaging results;
and the offset imaging section acquisition module is used for adding the imaging results of all shot gathers to obtain a final offset imaging section and improving the balance of deep imaging amplitude by using seismic source illumination.
It is a further object of the invention to provide a computer device comprising a memory and a processor, the memory storing a computer program which, when executed by the processor, causes the processor to perform the steps of:
obtaining input seismic records, seismic wavelets, density and velocity models; obtaining a seismic source and wave detection point continuation wave field by solving an elastic wave fluctuation equation; performing wave field separation by using vector Helmholtz decomposition to obtain separated pure longitudinal (P) wave and pure transverse (S) wave vector wave fields with amplitude preserved; calculating imaging conditions by using an elastic wave impedance sensitive kernel function to obtain high-precision PP and PS imaging results; and superposing all shot integrated image results to obtain a final offset imaging section.
It is another object of the present invention to provide a computer-readable storage medium storing a computer program which, when executed by a processor, causes the processor to perform the steps of:
obtaining input seismic records, seismic wavelets, density and velocity models; obtaining a seismic source and wave detection point continuation wave field by solving an elastic wave fluctuation equation; performing wave field separation by using vector Helmholtz decomposition to obtain separated pure longitudinal (P) wave and pure transverse (S) wave vector wave fields with amplitude preserved; calculating imaging conditions by using an elastic wave impedance sensitive kernel function to obtain high-precision PP and PS imaging results; and superposing all shot integrated image results to obtain a final offset imaging section.
Another objective of the present invention is to provide an information data processing terminal, which is used for implementing the decoupled elastic wave reverse time migration system.
By combining all the technical schemes, the invention has the advantages and positive effects that: according to the decoupling elastic wave reverse time migration method provided by the invention, firstly, vector Helmholtz decomposition is used for wave field separation to obtain pure longitudinal (P) wave and pure transverse (S) wave vector wave fields with amplitude preservation, and then, an elastic wave impedance sensitive kernel function is used for calculating imaging conditions to obtain high-precision PP and PS imaging results.
According to the scheme of the invention, after input seismic data, a seismic source wavelet and model parameters are obtained, extended seismic source and demodulator probe vector wave fields are obtained by solving an elastic wave field equation, then vertical and horizontal wave field separation is carried out by using vector Helmholtz decomposition, P wave and S wave vector wave fields of the seismic source and the demodulator probe are obtained, and finally cross-correlation imaging is carried out by using imaging conditions based on an elastic impedance sensitive core to obtain a final imaging result. Compared with the existing elastic wave reverse time migration imaging technology, the technical method of the invention adopts vector Helmholtz decomposition, can obtain the separation result of pure P wave field and S wave field of true amplitude, further avoid crosstalk noise in a migration profile and improve the imaging quality; in addition, the technical method uses the elastic impedance sensitive kernel function to calculate the imaging condition, can well avoid background noise in reverse time migration, and improve the imaging resolution and the signal-to-noise ratio so as to assist high-precision longitudinal and transverse wave combined seismic interpretation.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings needed to be used in the embodiments of the present invention will be briefly described below, and it is obvious that the drawings described below are only some embodiments of the present invention, and it is obvious for those skilled in the art that other drawings can be obtained according to the drawings without creative efforts.
Fig. 1 is a flowchart of a decoupling elastic wave reverse time migration method according to an embodiment of the present invention.
Fig. 2 is a schematic diagram of a decoupling elastic wave reverse time migration method provided in the embodiment of the present invention.
FIG. 3 is a block diagram of a decoupled elastic wave reverse time migration system according to an embodiment of the present invention;
in the figure: 1. an input data acquisition module; 2. a seismic source vector continuation wave field calculation module; 3. a wave detection point vector continuation wave field calculation module; 4. a longitudinal and transverse wave separation module; 5. an imaging result acquisition module; 6. an offset imaging profile acquisition module.
FIG. 4 is a longitudinal wave velocity model v provided by an embodiment of the present inventionp(x) Schematic representation of (a).
FIG. 5 shows a shear velocity model v provided by an embodiment of the present inventions(x) Schematic representation of (a).
FIG. 6 is a diagram of the results of conventional elastic wave reverse time migration PP imaging based on divergence and rotation wavefield separation and conventional cross-correlation imaging conditions, provided by an embodiment of the present invention.
FIG. 7 is a diagram of the results of conventional elastic wave reverse time migration PS imaging based on divergence and rotation wavefield separation and conventional cross-correlation imaging conditions, provided by an embodiment of the present invention.
Fig. 8 is a diagram illustrating elastic wave reverse time migration PP imaging results based on vector helmholtz decomposition and conventional cross-correlation imaging conditions according to an embodiment of the present invention.
Fig. 9 is a diagram illustrating elastic wave reverse time migration PS imaging results based on vector helmholtz decomposition and conventional cross-correlation imaging conditions according to an embodiment of the present invention.
FIG. 10 is a graphical illustration of decoupled elastic wave reverse time migration PP imaging results based on vector Helmholtz decomposition and impedance sensitive nuclear imaging conditions as provided by an embodiment of the present invention.
Fig. 11 is a diagram illustrating the result of decoupled elastic wave reverse time migration PS imaging based on vector helmholtz decomposition and impedance sensitive nuclear imaging conditions, according to an embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is further described in detail with reference to the following embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the invention and are not intended to limit the invention.
Aiming at the problems in the prior art, the invention provides a decoupling elastic wave reverse time migration method, an efficient algorithm and a processing flow, and the invention is described in detail below with reference to the accompanying drawings.
As shown in fig. 1, the decoupling elastic wave reverse time migration method provided in the embodiment of the present invention includes the following steps:
s101, obtaining input seismic records, seismic wavelets, density and velocity models;
s102, obtaining a seismic source and demodulator probe continuation wave field by solving an elastic wave fluctuation equation;
s103, obtaining separated true-amplitude pure longitudinal wave and pure transverse wave fields by using vector Helmholtz decomposition;
s104, performing offset imaging by using elastic impedance sensitive nuclear imaging conditions;
and S105, overlapping all shot integrated image results to obtain a final offset imaging section.
A schematic diagram of a decoupling elastic wave reverse time migration method provided by the embodiment of the invention is shown in fig. 2.
As shown in fig. 3, the decoupled elastic wave reverse time migration system provided in the embodiment of the present invention includes:
an input data acquisition module 1 for acquiring input data; wherein the input data comprises: longitudinal wave offset velocity model vp(x) Transverse wave migration velocity model vs(x) Density model rho (x), source function wavelet f (t), and multi-wave multi-component observation data d (x)rT), these inputs enhance the model and data base for subsequent wavefield computations and offset imaging;
a seismic source vector continuation wave field calculation module 2 for calculating the velocity model v according to the input seismic source wavelets f (t), the density model rho (x) and the longitudinal wave velocity model vp(x) And transverse wave velocity model vs(x) Calculating a seismic source vector continuation wave field u by solving a first-order velocity-stress elastic wave equation in a point source excitation modes(x,t)；
Wave detection point vector continuation wave field calculation module3, a density model rho (x) and a longitudinal wave velocity model v which are input according to the inputp(x) Transverse wave velocity model vs(x) And observed data d (x)rT), calculating the extended wave field u of the wave detection point vector by solving the elastic wave equationr(x,t)；
The longitudinal and transverse wave separation module 4 is used for performing longitudinal and transverse wave separation on wave fields of a seismic source and a wave detection point by adopting vector Helmholtz decomposition, and the innovative technical steps can obtain pure longitudinal wave and pure transverse wave vector wave fields with true amplitude, ensure the amplitude validity of an imaging result and avoid cross-interference noise of longitudinal and transverse waves;
the imaging result acquisition module 5 is used for performing cross-correlation imaging by using elastic impedance sensitive nuclear imaging conditions after acquiring separated seismic source and wave detection point wave fields, and the innovative processing step can automatically remove low-frequency noise when the imaging conditions are applied so as to obtain PP and PS imaging results with high signal-to-noise ratio and high resolution;
and the offset imaging section acquisition module 6 is used for adding the imaging results of all shot gathers to obtain a final offset imaging section, and improving the balance of deep imaging amplitude by using the seismic source illumination.
The technical solution of the present invention will be further described with reference to the following examples.
Example 1
The invention provides a decoupling elastic reverse time migration imaging method based on vector Helmholtz decomposition and impedance sensitive nuclear imaging conditions, aiming at solving the problems of low-frequency noise and cross wave crosstalk artifact in elastic wave reverse time migration. According to the method, firstly, vector Helmholtz decomposition is used for wave field separation to obtain amplitude-preserved pure longitudinal (P) wave and pure transverse (S) wave vector wave fields, then, an elastic wave impedance sensitive kernel function is used for calculating imaging conditions, and high-precision PP and PS imaging results are obtained.
The invention relates to a decoupling elastic reverse time migration imaging method based on an impedance sensitive nucleus, which comprises the following steps:
(1) obtaining input data, the input data comprising: longitudinal wave offset velocity model vp(x) Transverse wave migration velocity model vs(x) Density model rho (x), source function wavelet f (t), and multi-wave multi-component observation data d (x)r,t)。
(2) According to the input seismic source wavelets f (t), the density model rho (x) and the longitudinal wave velocity model vp(x) And transverse wave velocity model vs(x) Calculating the vector continuation wave field u of the seismic source by solving a first-order velocity-stress elastic wave equations(x, t), said first order velocity-stress elastic wave equation being of the form:
where x is the imaging spatial position, xsFor the seismic source location, t is the wavefield travel time, us＝[us,x us,y us,z]TFor the polarization velocity wave field, σs＝[σs,xx σs,yy σs,zz σs,xy σs,xz σs,yz]TFor the stress wavefield, T is the transposed symbol and δ is the kronecker function. L is a partial derivative matrix, C (x) is a parameter matrix, having the following expression:
(3) according to the input density model rho (x) and longitudinal wave velocity model vp(x) Transverse wave velocity model vs(x) And observed data d (x)rT), calculating the extended wave field u of the wave detection point vector by solving the elastic wave equationr(x, t), the expression of the elastic wave equation is:
wherein x isrFor the detector position, L and C (x) have the same expression as equations (2) and (3).
(4) And performing longitudinal and transverse wave separation on the wave fields of the seismic source and the wave detection point by adopting vector Helmholtz decomposition, wherein the specific expression is as follows:
wherein,
in order to perform the gradient operation, the method comprises the following steps,
in order to calculate the divergence, the method comprises the following steps,
in order to calculate the rotation degree,
for the separated longitudinal wave source wavefield,
is the shear wave source wavefield and,
is a longitudinal wave demodulator probe wave field,
is a shear wave demodulator probe wavefield. w is asAnd wrTo assist the vector wavefield, it can be obtained by solving the poisson equation:
where Δ is the laplace operator.
(5) After the separated seismic source and wave detection point wave fields are obtained, performing cross-correlation imaging by using elastic impedance sensitive nuclear imaging conditions to obtain PP and PS imaging results, wherein the specific expression is as follows:
wherein,
in order to be able to obtain a shear modulus,
in order to be the bulk modulus,
and I is an identity matrix.
(6) And adding the imaging results of all shot gathers to obtain a final offset imaging section, and improving the balance of deep imaging amplitude by using the source illumination.
The embodiment of the invention adopts at least one technical scheme which can achieve the following beneficial effects:
according to the scheme of the invention, after input seismic data, a seismic source wavelet and model parameters are obtained, extended seismic source and demodulator probe vector wave fields are obtained by solving an elastic wave field equation, then vertical and horizontal wave field separation is carried out by using vector Helmholtz decomposition, P wave and S wave vector wave fields of the seismic source and the demodulator probe are obtained, and finally cross-correlation imaging is carried out by using imaging conditions based on an elastic impedance sensitive core to obtain a final imaging result. Compared with the existing elastic wave reverse time migration imaging technology, the technical method adopts vector Helmholtz decomposition, can obtain an amplitude-preserved wave field separation result, further avoids crosstalk noise in a migration profile, and improves imaging quality; in addition, the technical method uses the elastic impedance sensitive kernel function to calculate the imaging condition, can well avoid background noise in reverse time migration, and improve the imaging resolution and the signal-to-noise ratio so as to assist high-precision longitudinal and transverse wave combined seismic interpretation.
Example 2
The following description is given as an explanation of the practical effects of the embodiment in the model.
The method provided by the invention is applied to the imaging of the international standard Marmousi model, and the ideal PP and PS imaging results are obtained. A longitudinal wave velocity model (as shown in fig. 4), a shear wave velocity model (as shown in fig. 5); adopting an observation aperture of 3km, and carrying out full waveform forward modeling by using a real velocity model to obtain observation seismic records, wherein the total number of the observation seismic records is 103 shot gather data; and inputting the source wavelet and the migration model to carry out source wave field continuation, reversely transmitting and observing the seismic record to obtain a wave detection point continuation wave field, and applying a sensitive kernel function imaging condition to obtain an imaging result (as shown in figures 10 and 11). For contrast imaging effects, imaging results obtained using the conventional elastic wave reverse time migration method based on both the rotation and divergence wave field decomposition and conventional cross-correlation imaging conditions are shown in fig. 6-7, and imaging results obtained using the vector helmholtz decomposition and conventional cross-correlation imaging conditions are shown in fig. 8-9. Compared with the conventional elastic wave reverse time migration imaging result (as shown in fig. 6-9), the scheme of the embodiment of the specification adopts elastic wave reverse time migration based on vector helmholtz decomposition and impedance sensitive nuclear imaging conditions, can automatically eliminate low-frequency noise in the reverse time migration, improves imaging section resolution and signal-to-noise ratio (as shown in fig. 10 and 11), and can greatly improve the accuracy of subsequent seismic interpretation. In the foregoing schematic diagram, Distance corresponds to the abscissa x and Depth corresponds to the ordinate z.
In the above embodiments, the implementation may be wholly or partially realized by software, hardware, firmware, or any combination thereof. When used in whole or in part, can be implemented in a computer program product that includes one or more computer instructions. When loaded or executed on a computer, cause the flow or functions according to embodiments of the invention to occur, in whole or in part. The computer may be a general purpose computer, a special purpose computer, a network of computers, or other programmable device. The computer instructions may be stored in a computer readable storage medium or transmitted from one computer readable storage medium to another, for example, the computer instructions may be transmitted from one website site, computer, server, or data center to another website site, computer, server, or data center via wire (e.g., coaxial cable, fiber optic, Digital Subscriber Line (DSL), or wireless (e.g., infrared, wireless, microwave, etc.)). The computer-readable storage medium can be any available medium that can be accessed by a computer or a data storage device, such as a server, a data center, etc., that includes one or more of the available media. The usable medium may be a magnetic medium (e.g., floppy Disk, hard Disk, magnetic tape), an optical medium (e.g., DVD), or a semiconductor medium (e.g., Solid State Disk (SSD)), among others.
The above description is only for the purpose of illustrating the present invention and the appended claims are not to be construed as limiting the scope of the invention, which is intended to cover all modifications, equivalents and improvements that are within the spirit and scope of the invention as defined by the appended claims.
Claims (10)
1. A decoupling elastic wave reverse time migration method is characterized in that the elastic wave reverse time migration method of vector wave field decoupling comprises the following steps: obtaining input seismic records, seismic wavelets, density and longitudinal and transverse wave velocity models; obtaining a seismic source and wave detection point continuation wave field by solving an elastic medium wave equation; elastic wave field separation is carried out by using vector Helmholtz decomposition, and separated true-amplitude pure longitudinal (P) wave and pure transverse (S) wave vector wave fields are obtained; calculating imaging conditions by using an elastic wave impedance sensitive kernel function to obtain high-precision longitudinal wave-longitudinal wave (PP) and longitudinal wave-transverse wave (PS) imaging results; and superposing all the single-shot imaging results to obtain a final offset imaging section.
2. A decoupled elastic wave reverse time migration method according to claim 1, wherein said vector wavefield decoupled elastic wave reverse time migration method comprises the steps of:
step one, acquiring input data; wherein the input data comprises: longitudinal wave offset velocity model vp(x) Transverse wave migration velocity model vs(x) Density model rho (x), source function wavelet f (t), multi-wave multi-componentObserved data d (x)r,t)；
Secondly, according to the input seismic source wavelet f (t), the density model rho (x) and the longitudinal wave velocity model vp(x) And transverse wave velocity model vs(x) Calculating an elastic vector wave field u extended at one side of the seismic source by solving a first-order velocity-stress elastic medium wave equation in a point source excitation modes(x,t)；
Thirdly, according to the input density model rho (x) and the longitudinal wave velocity model vp(x) Transverse wave velocity model vs(x) Using the observed data d (x)rT) as boundary condition, calculating extended elastic vector wave field u at one side of detection point by solving elastic wave equationr(x,t)；
Step four, performing vertical and horizontal wave field separation on the coupling elastic wave fields at the seismic source side and the wave detection point side by adopting vector Helmholtz decomposition;
step five, after the separated seismic source and wave detection point wave fields are obtained, performing zero-delay cross-correlation imaging by using elastic impedance sensitive nuclear imaging conditions to obtain PP and PS imaging results;
and step six, adding the imaging results of all shot gathers to obtain a final offset imaging section, and improving the balance of deep imaging amplitude by using seismic source illumination as a preconditioner.
3. The method for decoupling reverse time migration of an elastic wave according to claim 2, wherein in the second step, the first order velocity-stress elastic wave equation is in the form of:
where x is the imaging spatial position, xsFor the seismic source location, t is the wavefield travel time, us＝[us,x us,y us,z]TFor the seismic source side deflection velocity wave field, σs＝[σs,xx σs,yy σs,zz σs,xy σs,xz σs,yz]TAs seismic sourceA side stress wave field, T is a transposition symbol, and delta is a kronecker function; l is a first order partial derivative matrix, C (x) is a model parameter matrix, having the following expression:
4. the method for decoupling reverse time migration of an elastic wave according to claim 2, wherein in step three, the expression of the elastic wave equation is:
wherein x isrFor the position of the detection point, L is a partial derivative matrix, C (x) is a rigidity matrix formed by model parameters, and the specific expression is as follows:
5. the method for de-coupling elastic wave reverse time migration according to claim 2, wherein in step four, the source and demodulator point wavefields are separated by using vector Helmholtz decomposition, and the method has the following expression:
wherein,
in order to perform the gradient operation, the method comprises the following steps,
in order to calculate the divergence, the method comprises the following steps,
in order to calculate the rotation degree,
for the separated longitudinal wave source wavefield,
for the separated shear wave source wavefield,
for the separated longitudinal geophone point wavefield,
the separated transverse wave detection point wave field is obtained; w is asAnd wrTo assist the vector wavefield, it can be obtained by solving the following poisson equation:
where Δ is the laplace operator.
6. The method for decoupling elastic wave reverse time migration according to claim 2, wherein in step five, after the separated source and demodulator wave fields are obtained, cross-correlation imaging is performed by using elastic impedance sensitive nuclear imaging conditions, and PP and PS imaging results are obtained, wherein the expression is as follows:
wherein,
in order to be able to obtain a shear modulus,
in order to be the bulk modulus,
and I is an identity matrix.
7. An elastic wave reverse time migration system for implementing the vector wave field decoupling elastic wave reverse time migration method according to any one of claims 1 to 6, wherein the vector wave field decoupling elastic wave reverse time migration method comprises the following steps:
the input data acquisition module is used for acquiring input data; wherein the input data comprises: longitudinal wave offset velocity model vp(x) Transverse wave migration velocity model vs(x) Density model rho (x), source function wavelet f (t), and multi-wave multi-component observation data d (x)r,t)；
A seismic source vector continuation wave field calculation module for calculating the vector continuation wave field of the seismic source according to the input seismic source wavelets f (t), the density model rho (x) and the longitudinal wave velocity model vp(x) And transverse wave velocity model vs(x) Calculating a seismic source vector continuation wave field u by solving a first-order velocity-stress elastic wave equation in a point source excitation modes(x,t)；
A wave detection point vector continuation wave field calculation module for calculating the wave field according to the input density model rho (x) and the longitudinal wave velocity model vp(x) Transverse wave velocity model vs(x) Using the observed data d (x)rT) as boundary condition, calculating the extended wave field u of the wave detection point vector by solving the elastic wave equationr(x,t)；
The longitudinal and transverse wave field separation module is used for performing longitudinal and transverse wave separation on wave fields of the seismic source and the wave detection point by adopting vector Helmholtz decomposition;
the imaging result acquisition module is used for performing cross-correlation imaging by using elastic impedance sensitive nuclear imaging conditions after acquiring the separated seismic source and wave detection point wave fields to acquire PP and PS imaging results;
and the offset imaging section acquisition module is used for adding all the single-shot imaging results to obtain a final offset imaging section, and using the seismic source illumination as a preconditioner to improve the balance of the deep imaging amplitude.
8. A computer device, characterized in that the computer device comprises a memory and a processor, the memory storing a computer program which, when executed by the processor, causes the processor to carry out the steps of:
obtaining input seismic records, seismic wavelets, density and velocity models; obtaining a seismic source and wave detection point continuation wave field by solving an elastic wave fluctuation equation; performing vertical and horizontal wave field separation by using vector Helmholtz decomposition to obtain separated true-amplitude pure longitudinal (P) wave and pure transverse (S) wave vector wave fields; calculating imaging conditions by using an elastic wave impedance sensitive kernel function to obtain high-precision PP and PS imaging results; and superposing all shot integrated image results to obtain a final offset imaging section.
9. A computer-readable storage medium storing a computer program which, when executed by a processor, causes the processor to perform the steps of:
obtaining input seismic records, seismic wavelets, density and velocity models; obtaining a seismic source and wave detection point continuation wave field by solving an elastic wave fluctuation equation; performing wave field separation by using vector Helmholtz decomposition to obtain separated pure longitudinal (P) wave and pure transverse (S) wave vector wave fields with amplitude preserved; calculating imaging conditions by using an elastic wave impedance sensitive kernel function to obtain high-precision PP and PS imaging results; and superposing all shot integrated image results to obtain a final offset imaging section.
10. An information data processing terminal, characterized in that the information data processing terminal is used for implementing the decoupled elastic wave reverse time migration system according to claim 7.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110329092.8A CN112904426B (en) | 2021-03-27 | 2021-03-27 | Decoupling elastic wave reverse time migration method, system and application |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110329092.8A CN112904426B (en) | 2021-03-27 | 2021-03-27 | Decoupling elastic wave reverse time migration method, system and application |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112904426A true CN112904426A (en) | 2021-06-04 |
| CN112904426B CN112904426B (en) | 2022-09-30 |
Family
ID=76109272
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110329092.8A Active CN112904426B (en) | 2021-03-27 | 2021-03-27 | Decoupling elastic wave reverse time migration method, system and application |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112904426B (en) |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113885079A (en) * | 2021-08-23 | 2022-01-04 | 中国石油大学(华东) | Elastic wave field decoupling-based high-precision multi-azimuth reverse-time seismic source imaging method |
| CN114660659A (en) * | 2022-03-30 | 2022-06-24 | 中国矿业大学 | Longitudinal and transverse wave decoupling high-precision microseism positioning method |
| CN114814944A (en) * | 2022-03-15 | 2022-07-29 | 长江大学 | Elastic longitudinal and transverse wave field separation and reverse time migration imaging method based on divergence and rotation |
| CN114910962A (en) * | 2022-05-11 | 2022-08-16 | 中国石油大学(华东) | Decoupling method for transverse wave, fast longitudinal wave and slow longitudinal wave in two-phase medium |
| CN115373022A (en) * | 2022-01-11 | 2022-11-22 | 长江大学 | Elastic wave field Helmholtz decomposition method based on amplitude phase correction |
| CN117075197A (en) * | 2023-10-12 | 2023-11-17 | 中国石油大学(华东) | Transverse wave decoupling equation construction method for transverse isotropic dielectric wave field separation |
| CN117233838A (en) * | 2023-09-20 | 2023-12-15 | 长江大学 | Elastic quasi-longitudinal and transverse wave field separation and reverse time migration imaging method in two-dimensional VTI medium |
| CN117647838A (en) * | 2024-01-29 | 2024-03-05 | 东北石油大学三亚海洋油气研究院 | Angle-constrained reverse time offset imaging method |
Citations (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20120051182A1 (en) * | 2010-08-24 | 2012-03-01 | Snu R&B Foundation | Apparatus and method for imaging a subsurface using frequency-domain elastic reverse-time migration |
| CN105807315A (en) * | 2016-03-14 | 2016-07-27 | 中国石油大学(华东) | Elastic vector reverse time migration imaging method |
| CN105974470A (en) * | 2016-07-04 | 2016-09-28 | 中国石油大学(华东) | Multi-component seismic data least squares reverse time migration imaging method and system |
| US20170003410A1 (en) * | 2015-07-02 | 2017-01-05 | Korea Institute Of Geoscience And Mineral Resource | Elastic reverse-time migration system and method using absolute value function for improving the quality of subsurface structure imaging |
| CN106842300A (en) * | 2016-12-21 | 2017-06-13 | 中国石油大学(华东) | A kind of high efficiency multi-component seismic data true amplitude migration imaging method |
| CN106970416A (en) * | 2017-03-17 | 2017-07-21 | 中国地质科学院地球物理地球化学勘查研究所 | Elastic wave least square reverse-time migration system and method based on wave field separation |
| CN108802813A (en) * | 2018-06-13 | 2018-11-13 | 中国石油大学(华东) | A kind of multi-component seismic data offset imaging method and system |
| CN109143339A (en) * | 2018-08-14 | 2019-01-04 | 中国石油天然气集团有限公司 | Elastic reverse-time migration imaging method and device |
| CN109557582A (en) * | 2018-12-17 | 2019-04-02 | 中国石油大学(华东) | A kind of two dimension multi-component seismic data offset imaging method and system |
| CN110857998A (en) * | 2018-08-23 | 2020-03-03 | 中国石油化工股份有限公司 | Elastic reverse time migration method and system based on lowrank finite difference |
| CN111221037A (en) * | 2020-01-21 | 2020-06-02 | 中国石油大学(华东) | Decoupling elastic reverse time migration imaging method and device |
| CN111239804A (en) * | 2020-02-12 | 2020-06-05 | 中国石油大学(华东) | Elastic energy reverse time migration imaging method, device, equipment and system |
| US20200233112A1 (en) * | 2019-01-22 | 2020-07-23 | Saudi Arabian Oil Company | Avo imaging condition in elastic reverse time migration |
| CN112327359A (en) * | 2020-10-14 | 2021-02-05 | 山东省科学院海洋仪器仪表研究所 | Elastic reverse time migration method based on imaging energy flow vector |
-
2021
- 2021-03-27 CN CN202110329092.8A patent/CN112904426B/en active Active
Patent Citations (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20120051182A1 (en) * | 2010-08-24 | 2012-03-01 | Snu R&B Foundation | Apparatus and method for imaging a subsurface using frequency-domain elastic reverse-time migration |
| US20170003410A1 (en) * | 2015-07-02 | 2017-01-05 | Korea Institute Of Geoscience And Mineral Resource | Elastic reverse-time migration system and method using absolute value function for improving the quality of subsurface structure imaging |
| CN105807315A (en) * | 2016-03-14 | 2016-07-27 | 中国石油大学(华东) | Elastic vector reverse time migration imaging method |
| CN105974470A (en) * | 2016-07-04 | 2016-09-28 | 中国石油大学(华东) | Multi-component seismic data least squares reverse time migration imaging method and system |
| CN106842300A (en) * | 2016-12-21 | 2017-06-13 | 中国石油大学(华东) | A kind of high efficiency multi-component seismic data true amplitude migration imaging method |
| CN106970416A (en) * | 2017-03-17 | 2017-07-21 | 中国地质科学院地球物理地球化学勘查研究所 | Elastic wave least square reverse-time migration system and method based on wave field separation |
| CN108802813A (en) * | 2018-06-13 | 2018-11-13 | 中国石油大学(华东) | A kind of multi-component seismic data offset imaging method and system |
| CN109143339A (en) * | 2018-08-14 | 2019-01-04 | 中国石油天然气集团有限公司 | Elastic reverse-time migration imaging method and device |
| CN110857998A (en) * | 2018-08-23 | 2020-03-03 | 中国石油化工股份有限公司 | Elastic reverse time migration method and system based on lowrank finite difference |
| CN109557582A (en) * | 2018-12-17 | 2019-04-02 | 中国石油大学(华东) | A kind of two dimension multi-component seismic data offset imaging method and system |
| US20200233112A1 (en) * | 2019-01-22 | 2020-07-23 | Saudi Arabian Oil Company | Avo imaging condition in elastic reverse time migration |
| CN111221037A (en) * | 2020-01-21 | 2020-06-02 | 中国石油大学(华东) | Decoupling elastic reverse time migration imaging method and device |
| CN111239804A (en) * | 2020-02-12 | 2020-06-05 | 中国石油大学(华东) | Elastic energy reverse time migration imaging method, device, equipment and system |
| CN112327359A (en) * | 2020-10-14 | 2021-02-05 | 山东省科学院海洋仪器仪表研究所 | Elastic reverse time migration method based on imaging energy flow vector |
Non-Patent Citations (2)
| Title |
|---|
| 任丽娟等: "基于矢量波场分离的TI介质逆时偏移", 《物探与化探》 * |
| 李闯等: "平面波最小二乘逆时偏移编码策略分析", 《石油物探》 * |
Cited By (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113885079A (en) * | 2021-08-23 | 2022-01-04 | 中国石油大学(华东) | Elastic wave field decoupling-based high-precision multi-azimuth reverse-time seismic source imaging method |
| CN115373022A (en) * | 2022-01-11 | 2022-11-22 | 长江大学 | Elastic wave field Helmholtz decomposition method based on amplitude phase correction |
| CN115373022B (en) * | 2022-01-11 | 2024-04-19 | 长江大学 | Amplitude phase correction-based elastic wave field Helmholtz decomposition method |
| CN114814944A (en) * | 2022-03-15 | 2022-07-29 | 长江大学 | Elastic longitudinal and transverse wave field separation and reverse time migration imaging method based on divergence and rotation |
| CN114814944B (en) * | 2022-03-15 | 2022-10-21 | 长江大学 | Elastic longitudinal and transverse wave field separation and reverse time migration imaging method based on divergence and rotation |
| CN114660659A (en) * | 2022-03-30 | 2022-06-24 | 中国矿业大学 | Longitudinal and transverse wave decoupling high-precision microseism positioning method |
| CN114910962A (en) * | 2022-05-11 | 2022-08-16 | 中国石油大学(华东) | Decoupling method for transverse wave, fast longitudinal wave and slow longitudinal wave in two-phase medium |
| CN117233838B (en) * | 2023-09-20 | 2024-04-05 | 长江大学 | Elastic quasi-longitudinal and transverse wave field separation and reverse time migration imaging method in two-dimensional VTI medium |
| CN117233838A (en) * | 2023-09-20 | 2023-12-15 | 长江大学 | Elastic quasi-longitudinal and transverse wave field separation and reverse time migration imaging method in two-dimensional VTI medium |
| CN117075197A (en) * | 2023-10-12 | 2023-11-17 | 中国石油大学(华东) | Transverse wave decoupling equation construction method for transverse isotropic dielectric wave field separation |
| CN117075197B (en) * | 2023-10-12 | 2024-02-06 | 中国石油大学(华东) | Transverse wave decoupling equation construction method for transverse isotropic dielectric wave field separation |
| CN117647838A (en) * | 2024-01-29 | 2024-03-05 | 东北石油大学三亚海洋油气研究院 | Angle-constrained reverse time offset imaging method |
| CN117647838B (en) * | 2024-01-29 | 2024-04-12 | 东北石油大学三亚海洋油气研究院 | Angle-constrained reverse time offset imaging method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112904426B (en) | 2022-09-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112904426B (en) | Decoupling elastic wave reverse time migration method, system and application | |
| Gao et al. | Irregular seismic data reconstruction based on exponential threshold model of POCS method | |
| CN111221037B (en) | Decoupling elastic reverse time migration imaging method and device | |
| CN102156296A (en) | Elastic reverse time migration imaging method by combining seismic multi-component | |
| CN112487604B (en) | Long-time nonlinear drift compensation method for output data of marine gravimeter | |
| Li et al. | Wavelet-based higher order correlative stacking for seismic data denoising in the curvelet domain | |
| MX2011003850A (en) | Image domain signal to noise estimate. | |
| Zhang et al. | 3D simultaneous seismic data reconstruction and noise suppression based on the curvelet transform | |
| CN113885077A (en) | Multi-seismic-source seismic data separation method based on deep learning | |
| Ma et al. | Multichannel block sparse Bayesian learning reflectivity inversion with lp-norm criterion-based Q estimation | |
| CN111239805B (en) | Block constraint time-lapse seismic difference inversion method and system based on reflectivity method | |
| Saengduean et al. | Multi-source wavefield reconstruction combining interferometry and compressive sensing: application to a linear receiver array | |
| CN111538088B (en) | Offshore inclined cable wave field correction method | |
| CN112149614B (en) | Pre-stack well-seismic combined intelligent denoising method | |
| CN114740528A (en) | Pre-stack multi-wave joint inversion method based on ultramicro Laplace block constraint | |
| CN105866838B (en) | Seismic data low-frequency information compensation method and device | |
| Wang et al. | Inversion-based non-stationary normal moveout correction along with prestack high-resolution processing | |
| CN113655534A (en) | Nuclear magnetic resonance FID signal noise suppression method based on multi-linear singular value tensor decomposition | |
| Li et al. | Enhancement of the seismic data resolution through Q-compensated denoising based on dictionary learning | |
| Yang et al. | Deep learning with soft attention mechanism for small-scale ground roll attenuation | |
| CN114609671A (en) | Ghost wave attenuation method and device, computer equipment and readable storage medium | |
| Bing et al. | Crosshole acoustic velocity imaging with full-waveform spectral data: 2.5-D numerical simulations | |
| CN112799133B (en) | Wave field separation and time shift correction method for staggered grid simulation seismic data | |
| Liu et al. | A method for denoising active source seismic data via Fourier transform and spectrum reconstruction | |
| WO2023123951A1 (en) | Method and apparatus for improving das signal-to-noise ratio by means of local fk transform |
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 |