CN113219525B - Offset imaging defuzzification method and device - Google Patents
Offset imaging defuzzification method and device Download PDFInfo
- Publication number
- CN113219525B CN113219525B CN202010081486.1A CN202010081486A CN113219525B CN 113219525 B CN113219525 B CN 113219525B CN 202010081486 A CN202010081486 A CN 202010081486A CN 113219525 B CN113219525 B CN 113219525B
- Authority
- CN
- China
- Prior art keywords
- vector
- wave number
- point
- target scattering
- 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.)
- Active
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 209
- 238000000034 method Methods 0.000 title claims abstract description 49
- 239000013598 vector Substances 0.000 claims abstract description 255
- 230000006870 function Effects 0.000 claims description 136
- 239000006185 dispersion Substances 0.000 claims description 28
- 238000004590 computer program Methods 0.000 claims description 16
- 238000003860 storage Methods 0.000 claims description 8
- 238000004378 air conditioning Methods 0.000 claims description 4
- 230000009466 transformation Effects 0.000 claims description 3
- 230000005012 migration Effects 0.000 abstract description 15
- 238000013508 migration Methods 0.000 abstract description 15
- 238000010586 diagram Methods 0.000 description 13
- 238000012545 processing Methods 0.000 description 10
- 238000004364 calculation method Methods 0.000 description 5
- 230000008569 process Effects 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 3
- 238000005286 illumination Methods 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005755 formation reaction Methods 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 230000001131 transforming effect Effects 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000004519 manufacturing process Methods 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
- 239000003208 petroleum Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000005070 sampling Methods 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/50—Corrections or adjustments related to wave propagation
- G01V2210/51—Migration
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention discloses an offset imaging defuzzification method and device, the method comprises the steps of carrying out ray tracing on a seismic source and a receiving point corresponding to each seismic channel, and determining an imaging wave number vector at a target scattering point; determining a wave number domain point spread function of the seismic data according to the imaging wave number vector and a seismic source wavelet extracted from the seismic data; determining a spatial domain point spread function of the seismic data according to the wave number domain point spread function of the seismic data; and deblurring the primary imaging result of the seismic data by using a spatial domain point spread function of the seismic data to form a final imaging result of the seismic data. According to the method, forward modeling and migration are not needed, the spatial domain point spread function of the seismic data can be rapidly determined only on the basis of the imaging wave number vector and the seismic source wavelet at the target scattering point, and then the imaging result is deblurred by using the obtained spatial domain point spread function, so that the efficiency of seismic migration imaging deblurring can be improved.
Description
Technical Field
The invention relates to the technical field of petroleum geophysical exploration, in particular to an offset imaging defuzzification method and device.
Background
This section is intended to provide a background or context to the embodiments of the invention that are recited in the claims. The description herein is not admitted to be prior art by inclusion in this section.
The seismic migration imaging technology is one of key technical links of exploration seismology, has important roles in the research of seismic exploration theory and practical application, and plays an increasingly important role in the imaging of underground media which are increasingly complex and refined. The traditional offset imaging method is only the conjugate transpose of the Born forward operator, and can only generate a fuzzy structural imaging result for the band-limited seismic data acquired by the limited observation system. Furthermore, complex subsurface medium formations and irregularities in the spatial sampling of seismic data can also lead to migration artifacts and irregular imaging illumination, severely affecting the reliability of imaging amplitude. The least squares migration technique based on the linearized seismic inversion theory is an effective approach to solve the above-described problems. The technology solves the linear inverse problem corresponding to the least square offset by using an optimization method, and can theoretically eliminate the influence of factors such as acquisition illumination on an imaging result.
However, the classical data domain least square offset imaging method requires a conventional forward modeling and offset operation once for each iteration, and the complete inversion process often requires several tens of times of the calculation amount of the conventional offset, so the huge calculation amount severely restricts the data processing efficiency of the offset imaging method. The imaging domain least square offset imaging method needs to calculate and store huge Hessian matrix, even if the corresponding calculation process can be optimized through data coding and other technologies, the calculation and storage cost is still unacceptable, and the data processing efficiency of offset imaging is also low.
The point spread function is a localized Hessian matrix that describes the filtering and blurring effects of the seismic imaging system on the scatter rate of a single point, reflecting the illumination capability of the imaging system for a certain imaging point. In the case of known point spread functions, the conventional seismic imaging results may be deblurred using the point spread function, thereby improving the resolution and amplitude fidelity of the imaging results. However, when offset imaging deblurring is performed by using the point spread function, forward modeling and offset are required in the process of calculating the point spread function, and the calculation amount is also large, so that the data processing efficiency of seismic offset imaging deblurring is still low.
Therefore, the existing seismic offset imaging defuzzification has the problem of low data processing efficiency.
Disclosure of Invention
The embodiment of the invention provides an offset imaging defuzzification method, which is used for improving the data processing efficiency of seismic offset imaging defuzzification, and comprises the following steps of:
ray tracing is carried out on the seismic source and the receiving point corresponding to each seismic channel, and an imaging wave number vector at a target scattering point is determined;
determining a wave number domain point spread function of the seismic data formed by accumulating wave number domain point spread functions of each seismic channel according to an imaging wave number vector corresponding to the target scattering point and a seismic source wavelet extracted from the seismic data;
Determining a spatial domain point spread function of the seismic data according to the wave number domain point spread function of the seismic data;
deblurring the preliminary imaging result of the seismic data by using a spatial domain point spread function of the seismic data to form a final imaging result of the seismic data;
determining a wave number domain point spread function of each seismic trace based on the following formula:
wherein ,representing the wavenumber domain point spread function of each trace, F (ω) representing the source wavelet extracted from the seismic data, ω representing the seismic frequency, +.>Wavenumber vector representing the point spread function of the wavenumber domain, is->Representing the imaging wavenumber vector at the scattering point of the target.
The embodiment of the invention also provides an offset imaging defuzzification device for improving the data processing efficiency of seismic offset imaging defuzzification, which comprises the following steps:
the ray tracing module is used for carrying out ray tracing on the seismic source and the receiving point corresponding to each seismic channel and determining an imaging wave number vector at the target scattering point;
the wave number domain determining module is used for determining a wave number domain point spread function of the seismic data formed by accumulating wave number domain point spread functions of each seismic channel according to an imaging wave number vector corresponding to a target scattering point and a seismic source wavelet extracted from the seismic data;
The space domain determining module is used for determining a space domain point spread function of the seismic data according to the wave number domain point spread function of the seismic data;
the defuzzification module is used for defuzzifying the primary imaging result of the seismic data by utilizing the spatial domain point spread function of the seismic data to form a final imaging result of the seismic data;
determining a wave number domain point spread function of each seismic trace based on the following formula:
wherein ,representing the wavenumber domain point spread function of each trace, F (ω) representing the source wavelet extracted from the seismic data, ω representing the seismic frequency, +.>Wavenumber vector representing the point spread function of the wavenumber domain, is->Representing the imaging wavenumber vector at the scattering point of the target.
The embodiment of the invention also provides computer equipment, which comprises a memory, a processor and a computer program stored on the memory and capable of running on the processor, wherein the processor realizes the offset imaging defuzzification method when executing the computer program.
Embodiments of the present invention also provide a computer-readable storage medium storing a computer program for executing the offset imaging defuzzification method described above.
In the embodiment of the invention, the imaging wave number vector at the target scattering point is determined based on ray tracing, and then the spatial domain point spread function of the seismic data can be obtained according to the imaging wave number vector at the target scattering point and the seismic source wavelet, and finally the primary imaging result of the seismic data is defuzzified by utilizing the obtained spatial domain point spread function of the seismic data. According to the embodiment of the invention, forward modeling and migration are not needed, and the spatial domain point spread function of the seismic data can be rapidly determined only based on the imaging wavenumber vector and the seismic source wavelet at the target scattering point, so that the imaging result is deblurred by using the obtained spatial domain point spread function, and the efficiency of seismic migration imaging deblurring can be improved.
Drawings
In order to more clearly illustrate the embodiments of the invention or the technical solutions in the prior art, the drawings that are required in the embodiments or the description of the prior art will be briefly described, it being obvious that the drawings in the following description are only some embodiments of the invention, and that other drawings may be obtained according to these drawings without inventive effort for a person skilled in the art. In the drawings:
FIG. 1 is a flowchart of an implementation of an offset imaging defuzzification method provided by an embodiment of the present invention;
FIG. 2 is a flowchart of step 101 in the offset imaging defuzzification method according to an embodiment of the present invention;
FIG. 3-1 is a schematic diagram of a seismic imaging ray pair in an offset imaging defuzzification method according to an embodiment of the present invention;
FIG. 3 is a flowchart of another implementation of step 101 in the offset imaging defuzzification method according to the embodiment of the present invention;
FIG. 4 is a flowchart illustrating a step 102 in the offset imaging defuzzification method according to an embodiment of the present invention;
FIG. 5 is a flowchart illustrating a step 103 in the offset imaging defuzzification method according to an embodiment of the present invention;
FIG. 6 is a functional block diagram of an offset imaging defuzzification apparatus according to an embodiment of the present invention;
FIG. 7 is a block diagram illustrating a ray tracing module 601 in an offset imaging defuzzification apparatus according to an embodiment of the present invention;
FIG. 8 is another block diagram illustrating a ray tracing module 601 in an offset imaging defuzzification apparatus according to an embodiment of the present invention;
fig. 9 is a block diagram illustrating a wave number domain determining module 602 in an offset imaging defuzzification apparatus according to an embodiment of the present invention;
FIG. 10 is a block diagram illustrating a spatial domain determination module 603 in an offset imaging defuzzification apparatus according to an embodiment of the present invention;
FIG. 11 is a wavenumber domain point spread function of seismic data obtained by the offset imaging defuzzification method of the present invention, provided by an embodiment of the present invention;
FIG. 12 is a spatial domain point spread function of seismic data obtained using the offset imaging defuzzification method of the present invention, in accordance with an embodiment of the present invention.
Detailed Description
For the purpose of making the objects, technical solutions and advantages of the embodiments of the present invention more apparent, the embodiments of the present invention will be described in further detail with reference to the accompanying drawings. The exemplary embodiments of the present invention and their descriptions herein are for the purpose of explaining the present invention, but are not to be construed as limiting the invention.
Fig. 1 shows a flow of implementation of the offset imaging defuzzification method according to the embodiment of the present invention, and for convenience of description, only the portions relevant to the embodiment of the present invention are shown, which are described in detail below:
As shown in fig. 1, the offset imaging defuzzification method includes:
102, determining a wave number domain point spread function of the seismic data formed by accumulating wave number domain point spread functions of each seismic channel according to an imaging wave number vector corresponding to a target scattering point and a seismic source wavelet extracted from the seismic data;
In an embodiment of the invention, the seismic data includes multiple seismic traces, i.e., multiple seismic data. And when the imaging wave number vector at the target scattering point of each seismic trace is determined, ray tracing is carried out at the seismic source and receiving point corresponding to each seismic trace, so that the imaging wave number vector at the target scattering point is determined.
In addition, the seismic wavelets can be extracted from the seismic data, the extracted seismic wavelets are used as initial seismic wavelets, the initial seismic wavelets are further updated through regular inversion optimization, updated seismic wavelets are obtained, and the updated seismic wavelets are used as the seismic wavelets in the embodiment of the invention. It will be appreciated by those skilled in the art that the source wavelet may be extracted from the seismic data by other methods, and embodiments of the invention are not particularly limited.
After the imaging wavenumber vector corresponding to the target scattering point is obtained and the source sub-wave is extracted from the seismic data, the wavenumber domain point spread function of each seismic trace can be obtained based on the imaging wavenumber vector corresponding to the target scattering point and the source sub-wave extracted from the seismic data, and then the wavenumber domain point spread functions of each seismic trace are accumulated, so that the wavenumber domain point spread function of the seismic data is formed.
After the wave number domain point spread function of the seismic data is determined, the wave number domain point spread function of the seismic data is converted to obtain a spatial domain point spread function of the seismic data. The conversion of the wavenumber domain point spread function of the seismic data into the spatial domain point spread function of the seismic data is a well-known prior art for those skilled in the art, and any method in the prior art may be adopted by those skilled in the art, so that detailed description of the embodiment of the present invention is omitted.
And (3) after the spatial domain point spread function of the seismic data is obtained, the spatial domain point spread function of the seismic data can be utilized to deblur the primary imaging result of the seismic data. Any seismic data migration imaging method in the prior art can be adopted by a person skilled in the art to obtain a primary imaging result of the seismic data. The primary imaging result of the seismic data is deblurred by utilizing the spatial domain point spread function of the seismic data, so that the effect and resolution of seismic migration imaging can be improved.
In the embodiment of the invention, the imaging wave number vector at the target scattering point is determined based on ray tracing, and then the spatial domain point spread function of the seismic data can be obtained according to the imaging wave number vector at the target scattering point and the seismic source wavelet, and finally the primary imaging result of the seismic data is defuzzified by utilizing the obtained spatial domain point spread function of the seismic data. According to the embodiment of the invention, forward modeling and migration are not needed, and the spatial domain point spread function of the seismic data can be rapidly determined only based on the imaging wavenumber vector and the seismic source wavelet at the target scattering point, so that the imaging result is deblurred by using the obtained spatial domain point spread function, and the efficiency of seismic migration imaging deblurring can be improved.
Fig. 2 shows a flow chart of step 101 in the offset imaging defuzzification method according to the embodiment of the present invention, and for convenience of description, only the portions relevant to the embodiment of the present invention are shown, which are described in detail below:
in an embodiment of the present invention, in order to improve the efficiency of determining the imaging wavenumber vector, as shown in fig. 2, step 101 of performing ray tracing on the source and receiving points corresponding to each seismic trace, determining the imaging wavenumber vector at the target scattering point includes:
When the imaging wave number vector at the target scattering point is determined, firstly, ray tracing is carried out at the seismic source and the receiving point corresponding to each seismic channel so as to determine the propagation vector at the target scattering point. Further, after the propagation vector at the target scattering point is determined, the wave number vector at the target scattering point can be determined based on the propagation vector at the target scattering point and the seismic wave dispersion relation. The seismic wave frequency dispersion can cause propagation energy attenuation and phase velocity frequency dispersion of the seismic wave so as to influence the resolution of seismic imaging, so that the problem of energy attenuation and phase velocity frequency dispersion is required to be paid attention to in forward modeling and inversion of the seismic wave. Seismic dispersion refers to the relationship of wave number to frequency. After obtaining the wave number vector at the target scattering point, the imaging wave number vector at the target scattering point may be determined based on the wave number vector at the target scattering point.
Specifically, the following formula may be used to represent the seismic dispersion relationship:
wherein ,representing the wavenumber vector at the scattering point of the target, ω representing the seismic wave frequency, +.>Representing the propagation vector at the scattering point of the target.
In the embodiment of the invention, firstly, ray tracing is carried out on the seismic source and the receiving point corresponding to each seismic trace, the propagation vector at the target scattering point is determined, then, the wave number vector at the target scattering point is determined according to the propagation vector at the target scattering point and the seismic wave frequency dispersion relation, and finally, the imaging wave number vector at the target scattering point is determined according to the wave number vector at the target scattering point, so that the efficiency of determining the imaging wave number vector can be improved, and further, the offset imaging efficiency is improved.
Fig. 3-1 shows a schematic diagram of a seismic imaging ray pair in the offset imaging defuzzification method provided by the embodiment of the present invention, fig. 3 shows another implementation flow of step 101 in the offset imaging defuzzification method provided by the embodiment of the present invention, and for convenience of description, only the portion relevant to the embodiment of the present invention is shown in detail as follows:
in one embodiment of the present invention, to further improve the efficiency of determining the imaging wavenumber vector, the propagation vector at the target scattering point includes a source propagation vector at the target scattering point and a receiver propagation vector at the target scattering point, and the wavenumber vector at the target scattering point includes a source wavenumber vector at the target scattering point and a receiver wavenumber vector at the target scattering point. As shown in fig. 3-1 and 3, step 202, determining a wave number vector at the target scattering point according to the propagation vector at the target scattering point and the seismic wave dispersion relation, includes:
In an embodiment of the present invention, the propagation vector at the target scattering point comprises a source propagation vector at the target scattering pointThe amount and the received point propagation vector at the target scattering point. Accordingly, the wave number vector at the target scattering point includes a source wave number vector at the target scattering point and a receiver wave number vector at the target scattering point. For convenience of description, we respectively adopt Is->Representing a source propagation vector at the target scattering point, a receiver propagation vector at the target scattering point, a source wavenumber vector at the target scattering point, and a receiver wavenumber vector at the target scattering point.
Specifically, when determining an imaging wave number vector at a target scattering point, firstly, performing ray tracing on a seismic source and a receiving point corresponding to a seismic channel, and determining a seismic source propagation vector at the target scattering point based on the ray tracingAnd a receiving point propagation vector at the target scattering point +.>/>
Source propagation vector at the resulting target scattering pointThen, the propagation vector of the seismic source at the target scattering point can be determined>And the above-mentioned seismic wave dispersion relation, determining the source wavenumber vector +.>Specifically, the following formula can be used:
wherein ,representing the source wavenumber vector at the target scattering point, ω represents the seismic wave frequency, +.>Representing the source propagation vector at the target scattering point.
Receiving point propagation vectors at the resulting target scattering pointsAfter that, the propagation vector can be determined according to the receiving point at the target scattering point>And the above-mentioned seismic wave dispersion relation, and determining the wave number vector of the receiving point at the target scattering point +.>Specifically, the following formula can be used:
wherein ,representing the wave number vector of the receiving point at the scattering point of the target, ω representing the seismic frequency, +.>Representing the received point propagation vector at the scattering point of the object.
Source wavenumber vector at each of the determined target scattering points Wave number vector of receiving point at scattering point of targetThen, the wave number vector of the seismic source at the scattering point of the target can be based on +.>And the wave number vector of the receiving point at the scattering point of the object +.>An imaging wavenumber vector at the target scattering point is determined. Here we use +.>Representing the imaging wavenumber vector at the scattering point of the target.
In one embodiment of the present invention, to further increase the efficiency of determining the imaging wavenumber vector, step 303, is performed based on the source wavenumber vector at the target scattering pointAnd the wave number vector of the receiving point at the scattering point of the object +.>Determining the imaging wavenumber vector at the scattering point of the object>Comprising the following steps:
wavenumber vector of receiving point at scattering point of targetSeismic source wavenumber vector at scattering point with target +.>The result of the subtraction is taken as the formation at the target scattering pointImage wavenumber vector +.>
Specifically, the imaging wavenumber vector at the target scattering point can be determined using the following formula
wherein ,representing the imaging wavenumber vector at the scattering point of the target, is->A wave number vector representing the receiving point at the scattering point of the object, is->Representing the source wavenumber vector at the scattering point of the target.
In the embodiment of the invention, a seismic source wave number vector at a target scattering point is determined according to a seismic source propagation vector and a seismic wave frequency dispersion relation at the target scattering point; according to the receiving point propagation vector and the seismic wave frequency dispersion relation at the target scattering point, the receiving point wave number vector at the target scattering point is determined, and finally, according to the source wave number vector at the target scattering point and the receiving point wave number vector at the target scattering point, the imaging wave number vector at the target scattering point is determined, so that the efficiency of determining the imaging wave number vector can be further improved, and further the offset imaging efficiency is improved.
Fig. 4 shows a flow chart of step 102 in the offset imaging defuzzification method according to the embodiment of the present invention, and for convenience of description, only the portions relevant to the embodiment of the present invention are shown, which is described in detail below:
in an embodiment of the present invention, in order to improve the efficiency of determining the wave number domain point spread function, as shown in fig. 4, step 102, determining the wave number domain point spread function of the seismic data formed by accumulating the wave number domain point spread functions of each seismic trace according to the imaging wave number vector corresponding to the target scattering point and the source wavelet extracted from the seismic data, includes:
In determining the wavenumber domain point spread function of the seismic data, in view of the seismic data comprising a plurality of seismic traces, a wavenumber domain point spread function of each trace is first determined based on an imaging wavenumber vector at a target scattering point and a source wavelet extracted from the seismic data.
Specifically, the wavenumber domain point spread function for each trace may be determined based on the following formula:
wherein ,representing the wavenumber domain point spread function of each trace, F (ω) representing the source wavelet extracted from the seismic data, ω representing the seismic frequency, +.>Wavenumber vector representing the point spread function of the wavenumber domain, is->Representing the imaging wavenumber vector at the scattering point of the target. Wherein the wavenumber vector of the wavenumber domain point spread function +.>To the extent known, embodiments of the present invention are not specifically described.
After the wave number domain point spread function of each seismic channel is obtained, the wave number domain point spread function of each seismic channel is further accumulated, and the wave number domain point spread function of the seismic data can be obtained.
In the embodiment of the invention, firstly, the wave number domain point spread function of each seismic channel is determined according to the imaging wave number vector at the target scattering point and the seismic source wavelet extracted from the seismic data, and then the wave number domain point spread function of each seismic channel is accumulated, so that the wave number domain point spread function of the seismic data can be obtained, and the efficiency of determining the wave number domain point spread function can be improved.
Fig. 5 shows a flow chart of step 103 in the offset imaging defuzzification method according to the embodiment of the present invention, and for convenience of description, only the portions relevant to the embodiment of the present invention are shown, which are described in detail below:
In one embodiment of the present invention, to improve the efficiency of determining the spatial domain point spread function, as shown in fig. 5, step 103, determining the spatial domain point spread function of the seismic data according to the wavenumber domain point spread function of the seismic data includes:
The wave number domain point spread function of the seismic data is converted into the spatial domain point spread function of the seismic data, and two-dimensional inverse Fourier transform can be carried out on the wave number domain point spread function of the seismic data. The two-dimensional inverse fourier transform is performed on the wavenumber domain point spread function of the seismic data to obtain the spatial domain point spread function of the seismic data, which is a part of the prior art known to those skilled in the art, and the detailed description of the embodiment of the present invention is omitted.
In the embodiment of the invention, the wave number domain point spread function of the seismic data is subjected to two-dimensional inverse Fourier transform, and the spatial domain point spread function of the seismic data is determined, so that the efficiency of determining the spatial domain point spread function can be improved, and further the offset imaging efficiency is improved.
The embodiment of the invention also provides an offset imaging defuzzification device, which is described in the following embodiment. Since the principle of solving the problems of these devices is similar to that of the offset imaging defuzzification method, the implementation of these devices can be referred to as the implementation of the method, and the repetition is omitted.
Fig. 6 shows functional modules of the offset imaging defuzzification apparatus according to an embodiment of the present invention, and for convenience of explanation, only the portions relevant to the embodiment of the present invention are shown in detail as follows:
referring to fig. 6, each module included in the offset imaging defuzzification apparatus is configured to perform each step in the corresponding embodiment of fig. 1, and detailed descriptions of the corresponding embodiments of fig. 1 are omitted herein. In the embodiment of the present invention, the offset imaging defuzzification apparatus includes a ray tracing module 601, a wave number domain determining module 602, a spatial domain determining module 603, and a defuzzification module 604.
The ray tracing module 601 is configured to perform ray tracing on the source and receiving point corresponding to each seismic trace, and determine an imaging wavenumber vector at the target scattering point.
The wave number domain determining module 602 is configured to determine a wave number domain point spread function of the seismic data formed by accumulating wave number domain point spread functions of each seismic trace according to an imaging wave number vector corresponding to the target scattering point and a source wavelet extracted from the seismic data.
The spatial domain determining module 603 is configured to determine a spatial domain point spread function of the seismic data according to the wave number domain point spread function of the seismic data.
The defuzzification module 604 is configured to defuzzify a preliminary imaging result of the seismic data by using a spatial domain point spread function of the seismic data, so as to form a final imaging result of the seismic data.
In the embodiment of the present invention, the ray tracing module 601 determines the imaging wavenumber vector at the target scattering point based on ray tracing, and then the wave number domain determining module 602 and the spatial domain determining module 603 obtain the spatial domain point spread function of the seismic data according to the imaging wavenumber vector at the target scattering point and the source wavelet, and finally the defuzzification module 604 defuzzifies the primary imaging result of the seismic data by using the obtained spatial domain point spread function of the seismic data. According to the embodiment of the invention, forward modeling and migration are not needed, and the spatial domain point spread function of the seismic data can be rapidly determined only based on the imaging wavenumber vector and the seismic source wavelet at the target scattering point, so that the imaging result is deblurred by using the obtained spatial domain point spread function, and the efficiency of seismic migration imaging deblurring can be improved.
Fig. 7 shows a schematic structure of a ray tracing module 601 in an offset imaging defuzzification apparatus according to an embodiment of the present invention, and for convenience of explanation, only the portions related to the embodiment of the present invention are shown in detail as follows:
In an embodiment of the present invention, in order to improve the efficiency of determining the imaging wavenumber vector, referring to fig. 7, each unit included in the ray tracing module 601 is configured to perform each step in the corresponding embodiment of fig. 2, and detailed descriptions in fig. 2 and the corresponding embodiment of fig. 2 are omitted herein. In the embodiment of the present invention, the ray tracing module 601 includes a propagation vector determining unit 701, a wave number vector determining unit 702, and an imaging wave number vector determining unit 703.
The propagation vector determining unit 701 is configured to perform ray tracing on the source and receiving point corresponding to each seismic trace, and determine a propagation vector at the target scattering point.
The wave number vector determining unit 702 is configured to determine a wave number vector at the scattering point of the target according to the propagation vector at the scattering point of the target and the seismic wave dispersion relation.
An imaging wavenumber vector determination unit 703 for determining an imaging wavenumber vector at the target scattering point from the wavenumber vector at the target scattering point.
In the embodiment of the present invention, firstly, the propagation vector determining unit 701 performs ray tracing on the seismic source and the receiving point corresponding to each seismic trace to determine the propagation vector at the target scattering point, then the wave number vector determining unit 702 determines the wave number vector at the target scattering point according to the propagation vector at the target scattering point and the seismic wave frequency dispersion relation, and finally the imaging wave number vector determining unit 703 determines the imaging wave number vector at the target scattering point according to the wave number vector at the target scattering point, so that the efficiency of determining the imaging wave number vector can be improved, and further the offset imaging efficiency is improved.
Fig. 8 shows another schematic structure of a ray tracing module 601 in an offset imaging defuzzification apparatus according to an embodiment of the present invention, and for convenience of explanation, only the portions related to the embodiment of the present invention are shown in detail as follows:
in one embodiment of the present invention, to further improve the efficiency of determining the imaging wavenumber vector, the propagation vector at the target scattering point includes a source propagation vector at the target scattering point and a receiver propagation vector at the target scattering point, and the wavenumber vector at the target scattering point includes a source wavenumber vector at the target scattering point and a receiver wavenumber vector at the target scattering point. Referring to fig. 6, each unit included in the ray tracing module 601 is configured to perform each step in the corresponding embodiment of fig. 3, and detailed descriptions of fig. 3 and the corresponding embodiment of fig. 3 are omitted herein. In the embodiment of the present invention, the wave number vector determining unit 702 includes a source wave number vector determining subunit 801 and a receiving point wave number vector determining subunit 802.
The source wave number vector determination subunit 801 is configured to determine a source wave number vector at the target scattering point according to the source propagation vector and the seismic wave dispersion relationship at the target scattering point. A kind of electronic device with high-pressure air-conditioning system
A receiving point wave number vector determining subunit 802, configured to determine a receiving point wave number vector at the target scattering point according to the receiving point propagation vector and the seismic wave frequency dispersion relationship at the target scattering point.
The imaging wave number vector determination unit 703 includes an imaging wave number vector determination unit 803.
An imaging wavenumber vector determination unit 803 for determining an imaging wavenumber vector at the target scattering point from the source wavenumber vector at the target scattering point and the receiving point wavenumber vector at the target scattering point.
In an embodiment of the present invention, to further improve the efficiency of determining the imaging wavenumber vector, the imaging wavenumber vector determination subunit is specifically configured to use a result of subtracting the source wavenumber vector at the target scattering point from the receiving point wavenumber vector at the target scattering point as the imaging wavenumber vector at the target scattering point.
In the embodiment of the present invention, the source wave number vector determination subunit 801 determines a source wave number vector at a target scattering point according to a source propagation vector and a seismic wave dispersion relationship at the target scattering point; the receiving point wave number vector determining subunit 802 determines a receiving point wave number vector at the target scattering point according to the receiving point propagation vector at the target scattering point and the seismic wave frequency dispersion relation, and finally the imaging wave number vector determining subunit 803 determines an imaging wave number vector at the target scattering point according to the source wave number vector at the target scattering point and the receiving point wave number vector at the target scattering point, which can further improve the efficiency of determining the imaging wave number vector and further improve the efficiency of offset imaging.
Fig. 9 shows a schematic structure of a wave number domain determining module 602 in an offset imaging defuzzification apparatus according to an embodiment of the present invention, and for convenience of explanation, only the portions related to the embodiment of the present invention are shown in detail as follows:
in an embodiment of the present invention, in order to improve the efficiency of determining the point spread function in the wave number domain, referring to fig. 9, each unit included in the wave number domain determining module 602 is configured to perform each step in the corresponding embodiment of fig. 4, and detailed descriptions in fig. 4 and the corresponding embodiment of fig. 4 are omitted herein. In the embodiment of the present invention, the wave number domain determining module 602 includes a single-channel wave number domain determining unit 901 and an accumulating unit 902.
The single-channel wave number domain determining unit 901 is configured to determine a wave number domain point spread function of each channel of the seismic channel according to an imaging wave number vector at the target scattering point and a source wavelet extracted from the seismic data.
And the accumulating unit 902 is configured to accumulate the wave number domain point spread functions of each seismic trace, and determine the wave number domain point spread function of the seismic data.
In the embodiment of the present invention, firstly, the single-channel wave number domain determining unit 901 determines the wave number domain point spread function of each channel according to the imaging wave number vector at the target scattering point and the source wavelet extracted from the seismic data, and then the accumulating unit 902 accumulates the wave number domain point spread function of each channel, so as to obtain the wave number domain point spread function of the seismic data, and the efficiency of determining the wave number domain point spread function can be improved.
Fig. 10 shows a schematic structure of the spatial domain determining module 603 in the offset imaging defuzzification apparatus according to the embodiment of the present invention, and for convenience of explanation, only the portions related to the embodiment of the present invention are shown, which is described in detail below:
in an embodiment of the present invention, in order to improve the efficiency of determining the spatial domain point spread function, referring to fig. 10, each unit included in the spatial domain determining module 603 is configured to perform each step in the corresponding embodiment of fig. 5, and detailed descriptions of fig. 5 and the corresponding embodiment of fig. 5 are omitted herein. In the embodiment of the present invention, the spatial domain determining module 603 includes a transforming unit 1001.
The transformation unit 1001 is configured to perform two-dimensional inverse fourier transform on the wave number domain point spread function of the seismic data, and determine a spatial domain point spread function of the seismic data.
In the embodiment of the present invention, the transforming unit 1001 performs two-dimensional inverse fourier transform on the wave number domain point spread function of the seismic data, and determines the spatial domain point spread function of the seismic data, so that the efficiency of determining the spatial domain point spread function can be improved, and further the efficiency of offset imaging is improved.
Fig. 11 shows a schematic representation of a wave number domain point spread function of a seismic data obtained by using the offset imaging defuzzification method according to the embodiment of the present invention, and fig. 12 shows a schematic representation of a spatial domain point spread function of a seismic data obtained by using the offset imaging defuzzification method according to the embodiment of the present invention, where, for convenience of explanation, only the portions related to the embodiment of the present invention are shown, and the details are as follows:
FIG. 11 is a graph showing a wave number domain point spread function of seismic data obtained by the offset imaging defuzzification method of the present invention; FIG. 12 is a schematic representation of a spatial domain point spread function of seismic data obtained using the offset imaging defuzzification method of the present invention.
The embodiment of the invention also provides computer equipment, which comprises a memory, a processor and a computer program stored on the memory and capable of running on the processor, wherein the processor realizes the moving image defuzzification method when executing the computer program.
Embodiments of the present invention also provide a computer-readable storage medium storing a computer program for executing the above-described shift image defuzzification method.
In summary, the embodiment of the invention determines the imaging wavenumber vector at the target scattering point based on ray tracing, and then obtains the spatial domain point spread function of the seismic data according to the imaging wavenumber vector at the target scattering point and the source wavelet, and finally deblurs the primary imaging result of the seismic data by using the obtained spatial domain point spread function of the seismic data. According to the embodiment of the invention, forward modeling and migration are not needed, and the spatial domain point spread function of the seismic data can be rapidly determined only based on the imaging wavenumber vector and the seismic source wavelet at the target scattering point, so that the imaging result is deblurred by using the obtained spatial domain point spread function, and the efficiency of seismic migration imaging deblurring can be improved.
It will be appreciated by those skilled in the art that embodiments of the present invention may be provided as a method, system, or computer program product. Accordingly, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present invention may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems) and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flowchart illustrations and/or block diagrams, and combinations of flows and/or blocks in the flowchart illustrations and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
The foregoing description of the embodiments has been provided for the purpose of illustrating the general principles of the invention, and is not meant to limit the scope of the invention, but to limit the invention to the particular embodiments, and any modifications, equivalents, improvements, etc. that fall within the spirit and principles of the invention are intended to be included within the scope of the invention.
Claims (12)
1. An offset imaging defuzzification method, comprising:
ray tracing is carried out on the seismic source and the receiving point corresponding to each seismic channel, and an imaging wave number vector at a target scattering point is determined;
determining a wave number domain point spread function of the seismic data formed by accumulating wave number domain point spread functions of each seismic channel according to an imaging wave number vector corresponding to the target scattering point and a seismic source wavelet extracted from the seismic data;
determining a spatial domain point spread function of the seismic data according to the wave number domain point spread function of the seismic data;
deblurring the preliminary imaging result of the seismic data by using a spatial domain point spread function of the seismic data to form a final imaging result of the seismic data;
determining a wave number domain point spread function of each seismic trace based on the following formula:
wherein ,representing the wavenumber domain point spread function of each trace, F (ω) representing the source wavelet extracted from the seismic data, ω representing the seismic frequency, +.>Wavenumber vector representing the point spread function of the wavenumber domain, is->Representing the imaging wavenumber vector at the scattering point of the target.
2. The offset imaging defuzzification method of claim 1, wherein performing ray tracing at source and receiver points corresponding to each seismic trace to determine an imaging wavenumber vector at a target scattering point, comprises:
Ray tracing is carried out on the seismic source and the receiving point corresponding to each seismic channel, and the propagation vector of the target scattering point is determined;
determining a wave number vector at the target scattering point according to the propagation vector at the target scattering point and the seismic wave frequency dispersion relation;
and determining the imaging wave number vector at the target scattering point according to the wave number vector at the target scattering point.
3. The offset imaging defuzzification method of claim 2, wherein the propagation vectors at the target scattering points comprise a source propagation vector at the target scattering points and a receiving point propagation vector at the target scattering points, the wave number vectors at the target scattering points comprise a source wave number vector at the target scattering points and a receiving point wave number vector at the target scattering points, and determining the wave number vector at the target scattering points from the propagation vectors at the target scattering points and the seismic wave dispersion relationship comprises:
determining a seismic source wave number vector at the target scattering point according to the seismic source propagation vector and the seismic wave frequency dispersion relation at the target scattering point; a kind of electronic device with high-pressure air-conditioning system
According to the receiving point propagation vector and the seismic wave frequency dispersion relation at the target scattering point, determining the receiving point wave number vector at the target scattering point;
determining an imaging wavenumber vector at the target scattering point from the wavenumber vector at the target scattering point, comprising:
And determining an imaging wave number vector at the target scattering point according to the source wave number vector at the target scattering point and the receiving wave number vector at the target scattering point.
4. The offset imaging defuzzification method of claim 3, wherein determining the imaging wavenumber vector at the target scattering point from the source wavenumber vector at the target scattering point and the receiving wavenumber vector at the target scattering point comprises:
the result of subtracting the source wavenumber vector at the target scattering point from the receiving point wavenumber vector at the target scattering point is taken as the imaging wavenumber vector at the target scattering point.
5. The offset imaging defuzzification method of claim 1, wherein determining a spatial domain point spread function of the seismic data from a wavenumber domain point spread function of the seismic data comprises:
and performing two-dimensional inverse Fourier transform on the wave number domain point spread function of the seismic data, and determining the spatial domain point spread function of the seismic data.
6. An offset imaging defuzzification apparatus, comprising:
the ray tracing module is used for carrying out ray tracing on the seismic source and the receiving point corresponding to each seismic channel and determining an imaging wave number vector at the target scattering point;
The wave number domain determining module is used for determining a wave number domain point spread function of the seismic data formed by accumulating wave number domain point spread functions of each seismic channel according to an imaging wave number vector corresponding to a target scattering point and a seismic source wavelet extracted from the seismic data;
the space domain determining module is used for determining a space domain point spread function of the seismic data according to the wave number domain point spread function of the seismic data;
the defuzzification module is used for defuzzifying the primary imaging result of the seismic data by utilizing the spatial domain point spread function of the seismic data to form a final imaging result of the seismic data;
determining a wave number domain point spread function of each seismic trace based on the following formula:
wherein ,representing the wavenumber domain point spread function of each trace, F (ω) representing the source wavelet extracted from the seismic data, ω representing the seismic frequency, +.>Wavenumber vector representing the point spread function of the wavenumber domain, is->Representing the imaging wavenumber vector at the scattering point of the target.
7. The offset imaging defuzzification apparatus of claim 6, wherein the ray tracing module comprises:
the propagation vector determining unit is used for carrying out ray tracing on the seismic source and the receiving point corresponding to each seismic channel and determining a propagation vector at the target scattering point;
The wave number vector determining unit is used for determining the wave number vector at the target scattering point according to the propagation vector at the target scattering point and the seismic wave frequency dispersion relation;
and the imaging wave number vector determining unit is used for determining the imaging wave number vector at the target scattering point according to the wave number vector at the target scattering point.
8. The offset imaging defuzzification apparatus of claim 7, wherein the propagation vector at the target scattering point comprises a source propagation vector at the target scattering point and a receiving point propagation vector at the target scattering point, the wave number vector at the target scattering point comprises a source wave number vector at the target scattering point and a receiving point wave number vector at the target scattering point, the wave number vector determination unit comprising:
the focus wave number vector determination subunit is used for determining the focus wave number vector at the target scattering point according to the focus propagation vector and the seismic wave frequency dispersion relation at the target scattering point; a kind of electronic device with high-pressure air-conditioning system
The receiving point wave number vector determining subunit is used for determining the receiving point wave number vector at the target scattering point according to the receiving point propagation vector and the seismic wave frequency dispersion relation at the target scattering point;
an imaging wave number vector determination unit comprising:
and the imaging wave number vector determining subunit is used for determining the imaging wave number vector at the target scattering point according to the source wave number vector at the target scattering point and the receiving point wave number vector at the target scattering point.
9. The offset imaging defuzzification apparatus of claim 8, wherein the imaging wavenumber vector determination subunit is operable to determine a result of subtracting the source wavenumber vector at the target scattering point from the receive point wavenumber vector at the target scattering point as the imaging wavenumber vector at the target scattering point.
10. The offset imaging defuzzification apparatus of claim 6, wherein the spatial domain determination module comprises:
the transformation unit is used for carrying out two-dimensional inverse Fourier transformation on the wave number domain point spread function of the seismic data and determining the spatial domain point spread function of the seismic data.
11. A computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor implements the offset imaging defuzzification method of any of claims 1 to 5 when the computer program is executed by the processor.
12. A computer-readable storage medium, wherein the computer-readable storage medium stores a computer program for performing the offset imaging defuzzification method of any one of claims 1 to 5.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010081486.1A CN113219525B (en) | 2020-02-06 | 2020-02-06 | Offset imaging defuzzification method and device |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010081486.1A CN113219525B (en) | 2020-02-06 | 2020-02-06 | Offset imaging defuzzification method and device |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113219525A CN113219525A (en) | 2021-08-06 |
CN113219525B true CN113219525B (en) | 2023-04-25 |
Family
ID=77085580
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010081486.1A Active CN113219525B (en) | 2020-02-06 | 2020-02-06 | Offset imaging defuzzification method and device |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113219525B (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115166817B (en) * | 2022-06-30 | 2023-02-17 | 哈尔滨工程大学 | Ice sound positioning method based on ice layer modal group slowness difference characteristics |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108562937A (en) * | 2018-03-15 | 2018-09-21 | 东北石油大学 | A kind of seismic imaging method |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
NO322089B1 (en) * | 2003-04-09 | 2006-08-14 | Norsar V Daglig Leder | Procedure for simulating local preamp deep-migrated seismic images |
US7860333B2 (en) * | 2007-01-09 | 2010-12-28 | University Of Utah Research Foundation | Systems and methods for deblurring data corrupted by shift variant blurring |
US10677949B2 (en) * | 2017-01-05 | 2020-06-09 | China Petroleum & Chemical Corporation | Anisotropy matching filtering for attenuation of seismic migration artifacts |
CN107422379B (en) * | 2017-07-27 | 2019-01-11 | 中国海洋石油集团有限公司 | Multiple dimensioned seismic full-field shape inversion method based on local auto-adaptive convexification method |
CN108445532B (en) * | 2018-02-12 | 2019-11-08 | 中国石油天然气集团有限公司 | A kind of Depth Domain inverse migration method and device |
-
2020
- 2020-02-06 CN CN202010081486.1A patent/CN113219525B/en active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108562937A (en) * | 2018-03-15 | 2018-09-21 | 东北石油大学 | A kind of seismic imaging method |
Also Published As
Publication number | Publication date |
---|---|
CN113219525A (en) | 2021-08-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Zhu et al. | Integrating deep neural networks with full-waveform inversion: Reparameterization, regularization, and uncertainty quantification | |
CN110806602B (en) | Intelligent seismic data random noise suppression method based on deep learning | |
US9612351B2 (en) | System and method for estimating and attenuating noise in seismic data | |
CN105425301A (en) | Frequency domain three-dimensional irregular earthquake data reconstruction method | |
JP2020508454A (en) | Geophysical image generation using directional wavefield imaging | |
Chen et al. | Seismic inversion by hybrid machine learning | |
CN103105623A (en) | Data waveform processing method in seismic exploration | |
CN113219525B (en) | Offset imaging defuzzification method and device | |
WO2021055152A1 (en) | Noise attenuation methods applied during simultaneous source deblending and separation | |
Ferreira et al. | Internal multiple removal in offshore Brazil seismic data using the inverse scattering series | |
CN105259575A (en) | Method for fast predicting 3D surface-related multiples | |
Wang et al. | Fast 3D time-domain airborne EM forward modeling using random under-sampling | |
Sun et al. | Full-waveform Inversion Using A Learned Regularization | |
CN113031072B (en) | Multiple wave pressing method, device and equipment between virtual phase axis layers | |
Geetha et al. | An improved variational mode decomposition for seismic random noise attenuation using grasshopper optimization via shape dynamic time warping | |
Lee et al. | Zero-offset data estimation using CNN for applying 1D full waveform inversion | |
CN115170428A (en) | Noise reduction method for acoustic wave remote detection imaging graph | |
CN105319594A (en) | Fourier domain seismic data reconstruction method on the basis of least-square parametric inversion | |
CN106405504A (en) | Combined shear wave transformation and singular value decomposition ground penetrating radar data denoising method | |
Li et al. | Flexibility-residual BiSeNetV2 for GPR image decluttering | |
Qiu et al. | Mitigating the cycle-skipping of full-waveform inversion: An optimal transport approach with exponential encoding | |
CN114063159A (en) | Seismic surface wave velocity determination method and device | |
CN113311485A (en) | Seismic sedimentary feature enhanced filtering method and device | |
EP2304625A2 (en) | System and method for seismic trace analysis | |
Aghamiry et al. | Large-scale highly-accurate extended full waveform inversion using convergent Born series |
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 |