CN115061206A - Method and system for separating and imaging diffracted waves of remote detection acoustic logging - Google Patents
Method and system for separating and imaging diffracted waves of remote detection acoustic logging Download PDFInfo
- Publication number
- CN115061206A CN115061206A CN202210397636.9A CN202210397636A CN115061206A CN 115061206 A CN115061206 A CN 115061206A CN 202210397636 A CN202210397636 A CN 202210397636A CN 115061206 A CN115061206 A CN 115061206A
- Authority
- CN
- China
- Prior art keywords
- data
- diffracted wave
- diffracted
- rank
- imaging
- 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.)
- Pending
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 61
- 238000000034 method Methods 0.000 title claims abstract description 56
- 238000001514 detection method Methods 0.000 title claims abstract description 21
- 238000000926 separation method Methods 0.000 claims abstract description 41
- 238000002447 crystallographic data Methods 0.000 claims abstract description 32
- 238000004458 analytical method Methods 0.000 claims abstract description 18
- 238000012937 correction Methods 0.000 claims abstract description 10
- 238000012545 processing Methods 0.000 claims abstract description 10
- 239000011159 matrix material Substances 0.000 claims description 32
- 238000004088 simulation Methods 0.000 claims description 22
- 238000004364 calculation method Methods 0.000 claims description 20
- 238000000354 decomposition reaction Methods 0.000 claims description 15
- 238000010521 absorption reaction Methods 0.000 claims description 9
- 238000004422 calculation algorithm Methods 0.000 claims description 7
- 238000005516 engineering process Methods 0.000 claims description 7
- 238000005457 optimization Methods 0.000 claims description 7
- 230000008569 process Effects 0.000 claims description 7
- 238000009499 grossing Methods 0.000 claims description 6
- 230000009467 reduction Effects 0.000 claims description 5
- 238000010276 construction Methods 0.000 claims description 4
- 238000012795 verification Methods 0.000 claims description 4
- 230000009466 transformation Effects 0.000 claims description 3
- 230000007246 mechanism Effects 0.000 abstract description 4
- 230000009021 linear effect Effects 0.000 description 6
- 239000013598 vector Substances 0.000 description 4
- 230000001747 exhibiting effect Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000002349 favourable effect Effects 0.000 description 2
- 238000005755 formation reaction Methods 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 239000011435 rock Substances 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 238000010183 spectrum analysis Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 239000000945 filler Substances 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 230000000750 progressive 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/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
- G01V1/44—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
- G01V1/48—Processing data
-
- 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/52—Move-out correction
-
- 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/52—Move-out correction
- G01V2210/522—Dip move-out [DMO]
Abstract
The invention relates to a method and a system for separating and imaging diffracted waves of a remote detection acoustic logging, wherein the method comprises the following steps: processing the reflection data and the diffraction data based on a multi-channel singular analysis method to obtain separated diffraction wave signal data; and respectively carrying out high-precision travel time estimation and dip angle domain travel time correction on the diffracted wave signal data to obtain diffracted wave imaging data. The diffracted wave separation method under low-rank constraint can be researched based on multi-channel singular analysis, the problem of diffracted wave separation under the condition of strong reflected wave energy of a sound wave field is solved, and the separated diffracted waves have waveform consistency and integrity; secondly, by utilizing the dynamic characteristics of the diffracted wave such as polarity, amplitude and the like, a two-dimensional angle domain diffracted wave dynamic imaging mechanism is researched, and the diffracted wave imaging resolution is improved.
Description
Technical Field
The invention relates to the technical field of geological exploration, in particular to a method and a system for separating and imaging diffracted waves of a remote detection acoustic logging.
Background
Unconventional formations tend to be enriched with abundant geological resources. At present, research on exploration and development of unconventional strata is carried out successively at home and abroad. However, the unconventional stratum has complicated lithology, compactness and strong heterogeneity, which makes the implementation of fracturing difficult, and it is important to clear microstructures such as fractures and holes in the unconventional stratum. The remote detection acoustic imaging logging can detect stratum information of fracture and hole type reservoirs within tens of meters around the well, provides a new high-precision identification means for the fine description of unconventional strata, and how to remotely detect the acoustic imaging logging and evaluate the structure of the periphery of the unconventional strata is a work which has wide prospect and is very challenging.
The remote detection acoustic logging has important significance for evaluating unconventional formation information beside the stratum well, but the remote detection acoustic logging imaging method mainly takes reflected waves as main parts and is insufficient for the separation and imaging research of diffracted waves. Since diffracted waves are weak in energy and rapid in amplitude attenuation along with the continuation of propagation time, borehole mode waves, reflected waves and converted waves in the far detection acoustic wave data are strong in energy, and diffracted wave fields with weak energy are easily covered, so that diffracted wave research is extremely challenging. Therefore, how to separate the diffracted wave with weak energy from the mirror reflected wave with strong energy is an urgent problem to be solved.
Disclosure of Invention
In order to overcome the defects of the prior art, the invention aims to provide a method and a system for separating and imaging a far-detection acoustic logging diffracted wave.
In order to achieve the purpose, the invention provides the following scheme:
a method for separating and imaging a far detection acoustic logging diffracted wave comprises the following steps:
processing the reflection data and the diffraction data based on a multi-channel singular analysis method to obtain separated diffraction wave signal data;
and respectively carrying out high-precision travel time estimation and dip angle domain travel time correction on the diffracted wave signal data to obtain diffracted wave imaging data.
Preferably, the method further comprises the following steps:
and optimizing the rank of the reflection data and the rank of the diffraction data according to a singular value curve in the process of obtaining the separated diffraction wave signal data.
Preferably, the method further comprises the following steps:
and performing numerical simulation on the diffracted wave imaging data by adopting a method of two-dimensional staggered grid finite difference and complex frequency shift complete matching layer absorption boundary to obtain a simulation result.
Preferably, the processing the reflection data and the diffraction data based on the multi-channel singular analysis method to obtain separated diffracted wave signal data includes:
fourier transformation is carried out on the time domain pre-stack data or post-stack data in the reflection data and the diffraction data, and a Hankel matrix is constructed in a frequency domain;
performing effective decomposition on the Hankel matrix based on a singular value decomposition technology;
estimating thresholds of a reflected wave and a diffracted wave rank by combining a singular value curve curvature calculation formula;
reconstructing a diffracted wave Hankel matrix through a rank reduction algorithm;
and obtaining the diffracted wave signal data according to the reconstructed diffracted wave Hankel matrix.
Preferably, the optimizing the rank of the reflection data and the rank of the diffraction data according to the singular value curve includes:
smoothing and curvature calculation are carried out on the singular value curve to obtain a singular value curve curvature model;
and determining the rank of the reflection data and the position of the rank of the diffraction data according to the singular value curve curvature model.
A remote sounding acoustic logging diffracted wave separation and imaging system, comprising:
the separation module is used for processing the reflection data and the diffraction data based on a multi-channel singular analysis method to obtain separated diffracted wave signal data;
and the reflection module is used for respectively carrying out high-precision travel time estimation and inclination angle domain travel time correction on the diffracted wave signal data to obtain diffracted wave imaging data.
Preferably, the method further comprises the following steps:
and the optimization module is used for optimizing the rank of the reflection data and the rank of the diffraction data according to a singular value curve in the process of obtaining the separated diffracted wave signal data.
Preferably, the method further comprises the following steps:
and the verification module is used for performing numerical simulation on the diffracted wave imaging data by adopting a method of two-dimensional staggered grid finite difference and complex frequency shift complete matching layer absorption boundary to obtain a simulation result.
Preferably, the separation module specifically includes:
the matrix construction unit is used for performing Fourier transform on the time domain pre-stack or post-stack data in the reflection data and the diffraction data and constructing a Hankel matrix in a frequency domain;
the decomposition unit is used for effectively decomposing the Hankel matrix based on a singular value decomposition technology;
the estimation unit is used for estimating the threshold values of the reflected wave and the diffracted wave rank by combining a singular value curve curvature calculation formula;
the reconstruction unit is used for reconstructing a diffracted wave Hankel matrix through a rank subtraction algorithm;
and the separation unit is used for obtaining the diffracted wave signal data according to the reconstructed diffracted wave Hankel matrix.
Preferably, the optimization module specifically includes:
the calculation unit is used for carrying out smoothing and curvature calculation on the singular value curve to obtain a singular value curve curvature model;
and a determining unit configured to determine a rank of the reflection data and a position of the rank of the diffraction data based on the singular value curve curvature model.
According to the specific embodiment provided by the invention, the invention discloses the following technical effects:
the invention provides a method and a system for separating and imaging diffracted waves of a remote detection acoustic logging, wherein the method comprises the following steps: processing the reflection data and the diffraction data based on a multi-channel singular analysis method to obtain separated diffraction wave signal data; and respectively carrying out high-precision travel time estimation and dip angle domain travel time correction on the diffracted wave signal data to obtain diffracted wave imaging data. The diffracted wave separation method under low-rank constraint can be researched based on multi-channel singular analysis, the problem of diffracted wave separation under the condition of strong reflected wave energy of a sound wave field is solved, and the separated diffracted waves have waveform consistency and integrity; secondly, by utilizing the dynamic characteristics of the diffracted wave such as polarity, amplitude and the like, a two-dimensional angle domain diffracted wave dynamic imaging mechanism is researched, and the diffracted wave imaging resolution is improved.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings needed to be used in the embodiments will be briefly described below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art to obtain other drawings without inventive exercise.
FIG. 1 is a flow chart of a method in an embodiment provided by the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be obtained by a person skilled in the art without making any creative effort based on the embodiments in the present invention, belong to the protection scope of the present invention.
The invention aims to provide a method and a system for separating and imaging diffracted waves of a remote detection acoustic logging, which can solve the problem of diffracted wave separation under the condition of stronger reflected wave energy of an acoustic wave field, so that the separated diffracted waves have waveform consistency and integrity, and the imaging resolution of the diffracted waves is improved.
In order to make the aforementioned objects, features and advantages of the present invention comprehensible, embodiments accompanied with figures are described in further detail below.
Fig. 1 is a flowchart and a technical route schematic diagram of a method in an embodiment provided by the present invention, and as shown in fig. 1, the present invention provides a method for separating and imaging a far-detected acoustic logging diffracted wave, including:
step 100: processing the reflection data and the diffraction data based on a multi-channel singular analysis method to obtain separated diffraction wave signal data;
step 200: and respectively carrying out high-precision travel time estimation and dip angle domain travel time correction on the diffracted wave signal data to obtain diffracted wave imaging data.
According to the method, firstly, a diffracted wave separation method under low-rank constraint is researched based on multi-channel singular analysis, the problem of diffracted wave separation under the condition that the energy of a reflected wave of a sound wave field is strong is solved, and the separated diffracted waves have waveform consistency and integrity; secondly, researching a two-dimensional angle domain diffraction wave dynamic imaging mechanism by using dynamic characteristics such as the polarity and the amplitude of diffraction waves, and improving the imaging resolution of the diffraction waves; and finally, carrying out numerical simulation experiments and application of practical logging information of unconventional strata (such as hot dry rocks), and verifying the reliability of the method.
Preferably, the step 100 comprises:
fourier transformation is carried out on the time domain pre-stack data or post-stack data in the reflection data and the diffraction data, and a Hankel matrix is constructed in a frequency domain;
performing effective decomposition on the Hankel matrix based on a singular value decomposition technology;
estimating the threshold values of the reflected wave and the diffracted wave rank by combining a singular value curve curvature calculation formula;
reconstructing a diffracted wave Hankel matrix through a rank reduction algorithm;
and obtaining the diffracted wave signal data according to the reconstructed diffracted wave Hankel matrix.
Specifically, the first step in this embodiment is low-rank constrained rank subtraction and separation, and currently existing methods, such as a PWD filtering method and a CRS wavefront property method, mostly consider only the kinematic difference between the reflected wave and the diffracted wave, ignore the kinetic difference, and are not favorable for diffracted wave imaging. The kinematic characteristics are as follows: the reflected wave obeys snell's law and the diffracted wave obeys the huygens principle, so the reflected wave generally exhibits linear characteristics, while the diffracted wave generally exhibits hyperbolic characteristics, both of which exhibit completely different kinematic characteristics. Dynamic characteristics: the reflected wave is a wave field with higher energy formed by multiple interference superposition of the element diffraction wave field, and the diffraction wave field is generated by a secondary seismic source, and the energy of the diffraction wave field is generally one to two orders of magnitude weaker than that of the reflected wave. The MSSA-based diffraction wave high-quality separation technical process comprises the following steps:
a. carrying out Fourier transform on the time domain pre-stack or post-stack data, and constructing a Hankel matrix in a frequency domain;
b. carrying out effective decomposition on the Hankel matrix based on a Singular Value Decomposition (SVD) technology;
c. estimating the threshold values of the reflected wave and the diffracted wave rank by combining a singular value curve curvature calculation formula;
d. reconstructing a diffracted wave Hankel matrix through a rank reduction algorithm by using a rank reduction operator (truncation SVD);
e. and reconstructing time domain high-quality diffracted wave signal data according to the reconstructed diffracted wave Hankel matrix.
The method comprises the steps that firstly, a multi-channel singular spectrum analysis (MSSA) wave field signal generally has low-rank property, the signal is expressed by the acquaintance of singular values and singular vectors on the basis of the linear property and the energy intensity of the signal, and the rank is closely related to the number of the linear signals. Assuming df represents frequency domain data, it is at a certain frequency w 0 When, it can be expressed as:
hankel matrix construction:
data decomposition:
wherein Ω r is signal singular value, and Ω n is noise singular value.
Signal separation:
and secondly, reflected wave signals separated by the MSSA diffracted waves are in a linear form in a common offset distance or a post-stack domain and have strong energy, while diffracted waves are in a hyperbolic form and have weak energy, so that a theoretical basis is laid for the separation of the MSSA diffracted waves.
The mathematical model is as follows: d ═ D r +D d +D n 。
And (3) decomposing reflection and diffraction data:
wherein Ω r is a singular value of the reflected wave, Ω d is a singular value of the reflected wave, and Ω n is a singular value of the noise.
Diffraction wave separation:
preferably, the method further comprises the following steps:
and optimizing the rank of the reflection data and the rank of the diffraction data according to a singular value curve in the process of obtaining the separated diffraction wave signal data.
Specifically, the rank in this embodiment is preferably that in the diffracted wave separation process, the rank estimation of the reflected wave signal has a direct influence on the separation result. Since the singular values are directly related to the energy of the linear signal, the singular values change sharply around the threshold, i.e. the rank r1 for the reflected wave will appear at the maximum curvature of the singular value curve and the rank r2 threshold for the reflected wave + diffracted wave will appear at the second largest curvature of the curve.
Smoothed singular value curve: l ═ s 1 ，s 2 ，…s n ]。
Singular value curve curvature calculation:
further, in this embodiment, the high-precision diffracted wave imaging includes two steps, specifically as follows:
high-precision travel time estimation: slowness vectors are key parameters of ray tracing travel time calculation, single anisotropic parameters are estimated through low-order Taylor expansion in traditional travel time calculation, and then slowness vectors are estimated, but the low-order expansion causes insufficient travel time estimation accuracy, and high-precision imaging is not facilitated. Therefore, the travel time estimation precision can be effectively improved by adopting high-order Taylor expansion, and the imaging effect is improved.
Low-order Taylor expansion:
high-order Taylor expansion:
and secondly, a key parameter xi in the conventional correction formula for correcting the travel time of the dip angle domain represents the ratio of the horizontal distance between an observed imaging point and an underground diffraction point to the depth of the diffraction point, and under the actual condition, the position of the underground diffraction point is difficult to judge and know, and xi is a geological information parameter which cannot be obtained. And each underground point is assumed to be the position of the small-scale diffraction point, so that the method is simpler in form, and related parameters are reduced, so that the method has more practical application value and popularization.
Preferably, the method further comprises the following steps:
and performing numerical simulation on the diffracted wave imaging data by adopting a method of two-dimensional staggered grid finite difference and complex frequency shift complete matching layer absorption boundary to obtain a simulation result.
Specifically, a two-dimensional staggered grid finite difference method and a complex frequency shift complete matching layer absorption boundary are adopted, numerical simulation of remote detection acoustic logging in micro-structured unconventional strata with different depths, different scales and fillers is achieved, and diffracted wave separation and angle domain dynamic imaging methods are tested.
The embodiment also discloses feasibility analysis of the technical scheme, which comprises three parts of contents, specifically as follows:
(1) the feasibility analysis of the diffracted wave separation method most of the diffracted wave separation methods at present take the kinematic difference of reflected waves and diffracted waves as the separation basis, and the kinetic characteristic difference is ignored. In the post-stack domain, the reflected wave can be generally considered as a plane wave incident, with the in-phase axis exhibiting a linear morphology, while the diffracted wave incident as a non-plane wave, exhibiting a hyperbolic morphology, both exhibiting completely different kinematic properties. And (4) suppressing the reflected wave field by using the difference of the kinematic characteristics of the wave fields, and separating a diffracted wave field. However, the reflected wave and the diffracted wave are used as wave field responses of the geologic body, carry rich geological information representing the geologic body, and besides kinematic information, the dynamic information such as amplitude, phase and the like of the wave has an important role in researching characteristics such as geologic body structure, lithology and the like. In addition, since the diffracted waves are weak in energy carried by themselves and are easily influenced by other wave field information or noise, only the kinematic characteristics of the wave field are considered, which is not favorable for accurate imaging of the diffracted waves. The diffracted wave separation method based on low-rank constraint considers the kinematics and dynamics full wave field characteristic difference of the reflected wave and the diffracted wave field, and realizes high-quality separation of the diffracted wave field by utilizing multi-channel singular spectrum analysis (MSSA). Due to the difference of the kinematic characteristics of the same-phase axis forms of the reflected wave and the diffracted wave and the difference of the energy dynamics characteristics of the reflected wave and the diffracted wave, the reflected wave and the diffracted wave have different spatial distributions on a singular spectrum. By analyzing the mutual relation between the signal energy and the singular value, threshold value optimization is carried out by utilizing the curvature difference of the singular spectrum, the diffracted wave signal is reconstructed only by considering the singular value and the singular vector of the diffracted wave, and a high-quality diffracted wave field can be obtained.
(2) Diffraction wave imaging method feasibility analysis in the imaging domain, reflected waves and diffraction have different response degrees to disturbance of the offset speed, the reflected waves are basically not influenced by the speed disturbance and show weaker sensitivity, the diffraction wave imaging result is greatly influenced by the disturbance of the offset speed and show stronger sensitivity, and when the offset speed is higher or lower, the diffraction waves show a smiling face shape with a downward or upward opening, and when and only when the offset speed is accurate, the diffraction waves can be completely converged to the position of the diffraction point. In the dip domain, the reflected wave and the diffracted wave have the performance similar to that of the imaging domain, and the largest difference is that when the diffracted wave has an offset speed error in the dip domain and the imaging domain, the opening directions of the same-phase axes are opposite, so that mutual verification of the diffracted wave in the dip domain and the imaging domain is facilitated. Therefore, by adopting an approximate inclination angle domain diffracted wave travel time correction formula, high-precision diffracted wave imaging can be realized.
(3) The feasibility analysis finite difference method of the far detection sound wave numerical simulation method is to approximately replace the differential in a partial differential equation by the difference. The adopted differential formats are various, and for numerical simulation of the wave equation, the differential formats comprise a staggered grid, a rotary staggered grid and a Lebedev grid, wherein the staggered grid is the most mature and widely used differential format at present. In numerical simulation of acoustic logging, acoustic waves generate a tube wave at a solid-liquid interface, and the existence of the tube wave puts higher requirements on a perfectly matched absorption layer (PML) method. Several PML methods currently in use are Split PML (SPML), multi-axis PML (MPML), non-split PML (NPML) and complex frequency-shifted PML (CFS-PML). Compared with SPML, CFS-PML can absorb reflected guided waves from the calculation boundary more effectively; when the simulation time is longer, the SPML has the phenomenon of unstable numerical value; CFS-PML is faster than SPML operation time when numerical simulation parameters are the same. Among many PML methods, CFS-PML is the best choice for array acoustic wave FDTD simulation. Therefore, the two-dimensional staggered grid finite difference method and the complex frequency shift complete matching layer absorption boundary are adopted, and the remote detection acoustic logging numerical simulation and the diffracted wave separation and imaging method inspection can be realized.
Preferably, the optimizing the rank of the reflection data and the rank of the diffraction data according to the singular value curve includes:
smoothing and curvature calculation are carried out on the singular value curve to obtain a singular value curve curvature model;
and determining the rank of the reflection data and the position of the rank of the diffraction data according to the singular value curve curvature model.
The embodiment also discloses a far detection acoustic logging diffracted wave separation and imaging system, which comprises:
the separation module is used for processing the reflection data and the diffraction data based on a multi-channel singular analysis method to obtain separated diffracted wave signal data;
and the reflection module is used for respectively carrying out high-precision travel time estimation and inclination angle domain travel time correction on the diffracted wave signal data to obtain diffracted wave imaging data.
Preferably, the method further comprises the following steps:
and the optimization module is used for optimizing the rank of the reflection data and the rank of the diffraction data according to a singular value curve in the process of obtaining the separated diffracted wave signal data.
Preferably, the method further comprises the following steps:
and the verification module is used for performing numerical simulation on the diffracted wave imaging data by adopting a method of two-dimensional staggered grid finite difference and complex frequency shift complete matching layer absorption boundary to obtain a simulation result.
Preferably, the separation module specifically includes:
the matrix construction unit is used for carrying out Fourier transform on the time domain pre-stack or post-stack data in the reflection data and the diffraction data and constructing a Hankel matrix in a frequency domain;
the decomposition unit is used for effectively decomposing the Hankel matrix based on a singular value decomposition technology;
the estimation unit is used for estimating the threshold values of the reflected wave and the diffracted wave rank by combining a singular value curve curvature calculation formula;
the reconstruction unit is used for reconstructing a diffracted wave Hankel matrix through a rank subtraction algorithm;
and the separation unit is used for obtaining the diffracted wave signal data according to the reconstructed diffracted wave Hankel matrix.
Preferably, the optimization module specifically includes:
the calculation unit is used for carrying out smoothing and curvature calculation on the singular value curve to obtain a singular value curve curvature model;
and the determining unit is used for determining the rank of the reflection data and the position of the rank of the diffraction data according to the singular value curve curvature model.
The invention has the following beneficial effects:
based on multi-channel singular analysis, a diffracted wave separation method under low-rank constraint is researched, the problem of diffracted wave separation under the condition that the energy of a reflected wave of a sound wave field is strong is solved, and the separated diffracted waves have waveform consistency and integrity; secondly, researching a two-dimensional angle domain diffraction wave dynamic imaging mechanism by using dynamic characteristics such as the polarity and the amplitude of diffraction waves, and improving the imaging resolution of the diffraction waves; and finally, developing a numerical simulation experiment and the application of the actual data of the dry hot rock well logging, and verifying the reliability of the method.
The embodiments in the present description are described in a progressive manner, each embodiment focuses on differences from other embodiments, and the same and similar parts among the embodiments are referred to each other.
The principles and embodiments of the present invention have been described herein using specific examples, which are provided only to help understand the method and the core concept of the present invention; meanwhile, for a person skilled in the art, according to the idea of the present invention, the specific embodiments and the application range may be changed. In view of the above, the present disclosure should not be construed as limiting the invention.
Claims (10)
1. A method for separating and imaging a far-detection acoustic logging diffracted wave is characterized by comprising the following steps:
processing the reflection data and the diffraction data based on a multi-channel singular analysis method to obtain separated diffraction wave signal data;
and respectively carrying out high-precision travel time estimation and dip angle domain travel time correction on the diffracted wave signal data to obtain diffracted wave imaging data.
2. The method for borehole diffraction separation and imaging for remote sensing acoustic logging according to claim 1, further comprising:
and optimizing the rank of the reflection data and the rank of the diffraction data according to a singular value curve in the process of obtaining the separated diffraction wave signal data.
3. The method of claim 1, further comprising:
and performing numerical simulation on the diffracted wave imaging data by adopting a method of two-dimensional staggered grid finite difference and complex frequency shift complete matching layer absorption boundary to obtain a simulation result.
4. The method for separating and imaging diffracted waves from a far-sounding acoustic log according to claim 1, wherein the processing the reflection data and the diffraction data based on a multi-channel singular analysis method to obtain separated diffracted wave signal data comprises:
fourier transformation is carried out on the time domain pre-stack data or post-stack data in the reflection data and the diffraction data, and a Hankel matrix is constructed in a frequency domain;
performing effective decomposition on the Hankel matrix based on a singular value decomposition technology;
estimating the threshold values of the reflected wave and the diffracted wave rank by combining a singular value curve curvature calculation formula;
reconstructing a diffracted wave Hankel matrix through a rank reduction algorithm;
and obtaining the diffracted wave signal data according to the reconstructed diffracted wave Hankel matrix.
5. The method of claim 2, wherein optimizing the rank of the reflection data and the rank of the diffraction data according to a singular value curve comprises:
smoothing and curvature calculation are carried out on the singular value curve to obtain a singular value curve curvature model;
and determining the rank of the reflection data and the position of the rank of the diffraction data according to the singular value curve curvature model.
6. A system for remotely detecting acoustic logging diffracted wave separation and imaging, comprising:
the separation module is used for processing the reflection data and the diffraction data based on a multi-channel singular analysis method to obtain separated diffracted wave signal data;
and the reflection module is used for respectively carrying out high-precision travel time estimation and inclination angle domain travel time correction on the diffracted wave signal data to obtain diffracted wave imaging data.
7. The far detection acoustic logging diffracted wave separation and imaging system of claim 6, further comprising:
and the optimization module is used for optimizing the rank of the reflection data and the rank of the diffraction data according to a singular value curve in the process of obtaining the separated diffraction wave signal data.
8. The far detection acoustic logging diffracted wave separation and imaging system of claim 6, further comprising:
and the verification module is used for performing numerical simulation on the diffracted wave imaging data by adopting a method of two-dimensional staggered grid finite difference and complex frequency shift complete matching layer absorption boundary to obtain a simulation result.
9. The far detection acoustic logging diffracted wave separation and imaging system of claim 6, wherein the separation module specifically comprises:
the matrix construction unit is used for carrying out Fourier transform on the time domain pre-stack or post-stack data in the reflection data and the diffraction data and constructing a Hankel matrix in a frequency domain;
the decomposition unit is used for effectively decomposing the Hankel matrix based on a singular value decomposition technology;
the estimation unit is used for estimating the threshold values of the reflected wave and the diffracted wave rank by combining a singular value curve curvature calculation formula;
the reconstruction unit is used for reconstructing a diffracted wave Hankel matrix through a rank subtraction algorithm;
and the separation unit is used for obtaining the diffracted wave signal data according to the reconstructed diffracted wave Hankel matrix.
10. The far detection acoustic logging diffracted wave separation and imaging system of claim 7, wherein the optimization module specifically comprises:
the calculation unit is used for carrying out smoothing and curvature calculation on the singular value curve to obtain a singular value curve curvature model;
and the determining unit is used for determining the rank of the reflection data and the position of the rank of the diffraction data according to the singular value curve curvature model.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210397636.9A CN115061206A (en) | 2022-04-15 | 2022-04-15 | Method and system for separating and imaging diffracted waves of remote detection acoustic logging |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210397636.9A CN115061206A (en) | 2022-04-15 | 2022-04-15 | Method and system for separating and imaging diffracted waves of remote detection acoustic logging |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN115061206A true CN115061206A (en) | 2022-09-16 |
Family
ID=83196482
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210397636.9A Pending CN115061206A (en) | 2022-04-15 | 2022-04-15 | Method and system for separating and imaging diffracted waves of remote detection acoustic logging |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115061206A (en) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20180292553A1 (en) * | 2017-01-10 | 2018-10-11 | Institute Of Geology And Geophysics, Chinese Academy Of Sciences | Method and Apparatus for Separating Seismic Diffracted Wave |
| CN113945968A (en) * | 2021-10-19 | 2022-01-18 | 中国矿业大学(北京) | Diffracted wave imaging method and device for discontinuous geologic body and electronic equipment |
| CN113945967A (en) * | 2021-10-14 | 2022-01-18 | 中国矿业大学(北京) | Diffraction wave separation method and device |
| CN114200525A (en) * | 2021-12-10 | 2022-03-18 | 河北地质大学 | Self-adaptive multi-channel singular spectrum analysis seismic data denoising method |
-
2022
- 2022-04-15 CN CN202210397636.9A patent/CN115061206A/en active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20180292553A1 (en) * | 2017-01-10 | 2018-10-11 | Institute Of Geology And Geophysics, Chinese Academy Of Sciences | Method and Apparatus for Separating Seismic Diffracted Wave |
| CN113945967A (en) * | 2021-10-14 | 2022-01-18 | 中国矿业大学(北京) | Diffraction wave separation method and device |
| CN113945968A (en) * | 2021-10-19 | 2022-01-18 | 中国矿业大学(北京) | Diffracted wave imaging method and device for discontinuous geologic body and electronic equipment |
| CN114200525A (en) * | 2021-12-10 | 2022-03-18 | 河北地质大学 | Self-adaptive multi-channel singular spectrum analysis seismic data denoising method |
Non-Patent Citations (5)
| Title |
|---|
| 张亮: "基于紧致交错网格的井间地震波场数值模拟", 中国优秀硕士学位论文全文数据库, no. 03, pages 011 - 403 * |
| 杨茜娜: "地震波数值模拟中的完全匹配层参数取值研究", 中国优秀硕士学位论文全文数据库, no. 01, pages 011 - 300 * |
| 林朋: "地震绕射波分离与高精度成像方法研究", 中国博士学位论文全文数据库, no. 01, pages 011 - 146 * |
| 田坤 等: "实现复频移完全匹配层吸收边界条件的递推卷积方法", 吉林大学学报(地球科学版), vol. 43, no. 3, pages 1022 - 1032 * |
| 谢俊法 等: "基于无分裂复频移卷积完全匹配层边界的黏弹介质勒夫波模拟", 地震学报, vol. 38, no. 22, pages 244 - 259 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN113759425B (en) | Method and system for evaluating filling characteristics of deep paleo-karst reservoir stratum by well-seismic combination | |
| CN111239802B (en) | Deep learning speed modeling method based on seismic reflection waveform and velocity spectrum | |
| CN113759424B (en) | Karst reservoir filling analysis method and system based on spectral decomposition and machine learning | |
| CN106405651B (en) | Full waveform inversion initial velocity model construction method based on logging matching | |
| CN104237940A (en) | Diffracted wave imaging method and diffracted wave imaging device based on dynamical features | |
| CN103758511B (en) | A kind of method and device of the down-hole reverse-time migration hidden reservoir of imaging identification | |
| CN104570072A (en) | Method for modeling reflection coefficient of spherical PP wave in viscoelastic medium | |
| US9952341B2 (en) | Systems and methods for aligning a monitor seismic survey with a baseline seismic survey | |
| CN102053267A (en) | Method for separating VSP (vertical seismic profiling) wave field based on parametric inversion during seismic profile data processing | |
| CN102937720A (en) | Method for improving seismic data resolution ratio for well control | |
| Li et al. | Identification of fractured carbonate vuggy reservoirs in the S48 well area using 3D 3C seismic technique: A case history from the Tarim Basin | |
| Li et al. | Estimation of porosity, fluid bulk modulus, and stiff-pore volume fraction using a multitrace Bayesian amplitude-variation-with-offset petrophysics inversion in multiporosity reservoirs | |
| CN102565852B (en) | Angle domain pre-stack offset data processing method aiming to detect oil-gas-bearing property of reservoir | |
| Ning et al. | High-frequency surface-wave imaging from traffic-induced noise by selecting in-line sources | |
| Eisen et al. | Age estimates of isochronous reflection horizons by combining ice core, survey, and synthetic radar data | |
| CN103558637B (en) | Based on the detection method far away of three component sensor | |
| Zhou et al. | Passive surface-wave waveform inversion for source-velocity joint imaging | |
| CN112462427A (en) | Multi-component seismic data amplitude-preserving angle domain common imaging point gather extraction method and system | |
| CN104570090A (en) | Method for extracting full-waveform inversion noise filtering operator and performing noise filtering through full-waveform inversion noise filtering operator | |
| CN115061206A (en) | Method and system for separating and imaging diffracted waves of remote detection acoustic logging | |
| CN101937101A (en) | Method for identifying whether time-lapse seism is implemented or not | |
| CN105929453B (en) | State estimation method for infinite distribution time lag neural network system with channel fading | |
| Ren et al. | Pre-stack elastic reverse time migration in tunnels based on cylindrical coordinates | |
| Zhao et al. | The 3D global pre-stack seismic inversion in the time-frequency mixed domain | |
| Tang et al. | Numerical Simulation and Characteristic Analysis of Remote Detection Acoustic Reflection Imaging Logging of Near-Borehole Fracture and Vuggy Reservoirs |
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 |