A scattered wave post-stack imaging method and device, electronic equipment and storage medium
By performing regional segmentation and scattering Fourier filtering on the reverse-time offset image, the problem of high computational cost in scattered wave imaging is solved, and fine imaging and efficient identification of small-scale slits are achieved.
Patent Information
- Authority / Receiving Office
- CN · China
- Patent Type
- Patents(China)
- Current Assignee / Owner
- CHINA PETROLEUM & CHEMICAL CORP
- Filing Date
- 2022-09-08
- Publication Date
- 2026-06-09
AI Technical Summary
In existing technologies, scattered wave imaging requires wave field separation during the reverse-time offset wave field propagation process, resulting in high computational costs and making it difficult to efficiently image small-scale slit-hole anomalies.
By dividing the reverse-time offset image into blocks in the spatial dimension and filtering the wavenumber domain using the scattering Fourier filter operator, combined with discrete Fourier transform and inverse Fourier transform, post-stack imaging of scattered waves is achieved.
This reduces computational load and hardware costs, enabling detailed imaging of small-scale slit structures and improving recognition capabilities.
Smart Images
Figure CN117706616B_ABST
Abstract
Description
Technical Field
[0001] This application belongs to the field of seismic wave imaging technology, specifically relating to a method, device, electronic device, and storage medium for post-stack imaging of scattered waves. Background Technology
[0002] Fine imaging of small-scale slots and cavities is a core challenge in high-precision ultra-deep imaging in Northwest China. Reverse-time migration can redirect observed reflected energy to the correct imaging location, but information about small-scale anomalies is primarily contained in scattered waves, making amplitude less dominant and difficult to highlight slot and cavity anomalies. Scattered wave imaging technology, by fully utilizing scattered signals, has advantages in imaging small-scale anomalies and slot and cavity structures; however, conventional scattered wave imaging techniques require wavefield separation during the reverse-time migration wavefield propagation to extract the scattered signals for imaging, incurring very high computational costs. Therefore, each solution in existing technologies has its limitations. Summary of the Invention
[0003] Based on the above technical problems, this application proposes a method, apparatus, electronic device and storage medium for post-stack imaging of scattered waves.
[0004] In a first aspect, this application proposes a method for post-stack imaging of scattered waves, including:
[0005] The reverse time-shifted image is divided into regions in the spatial dimension to obtain the reverse time-shifted image of each region.
[0006] A three-dimensional discrete Fourier transform is performed on the reverse time-shifted image of each region to obtain the reverse time-shifted image in the wavenumber domain.
[0007] The wavenumber domain reverse-time-shifted image is filtered based on the scattering Fourier filter operator to obtain the wavenumber domain scattered wave imaging result;
[0008] The wavenumber domain scattered wave imaging results are subjected to discrete inverse Fourier transform to obtain the spatial domain scattered wave imaging results;
[0009] The spatial domain scattered wave imaging results of each region are stitched together to obtain a complete scattered wave stacked imaging result.
[0010] The process of obtaining the scattering Fourier filter operator is as follows:
[0011] Extract the spatial wavenumber spectrum corresponding to the Fourier transform coordinates from the reverse time-shifted image in the wavenumber domain;
[0012] Perform an arctangent operation on the spatial wavenumber spectrum corresponding to the Fourier transform coordinates to obtain the scattering angle corresponding to the spatial wavenumber spectrum;
[0013] The scattering Fourier filter operator is determined based on the scattering angle corresponding to the spatial wavenumber spectrum.
[0014] The arctangent of the spatial wavenumber spectrum corresponding to the Fourier transform coordinates is performed using the first calculation formula to obtain the scattering angle corresponding to the spatial wavenumber spectrum. The first calculation formula is as follows:
[0015]
[0016] Where θ is the scattering angle corresponding to the spatial wavenumber distribution, α is the damping factor, kx is the spatial wavenumber spectrum in the x-direction corresponding to the Fourier transform coordinates, ky is the spatial wavenumber spectrum in the y-direction corresponding to the Fourier transform coordinates, and kz is the spatial wavenumber spectrum in the z-direction corresponding to the Fourier transform coordinates.
[0017] The scattering Fourier filter operator is determined based on the scattering angle corresponding to the spatial wavenumber spectrum using preset rules, as follows:
[0018]
[0019] Where g is the scattering Fourier filter operator, θ is the scattering angle corresponding to the space wavenumber spectrum, θ0 is the scattering angle threshold corresponding to the space wavenumber spectrum, and r is the attenuation coefficient.
[0020] The second calculation formula is used to filter the reverse time-shifted image in the wavenumber domain based on the scattering Fourier filter operator. The second calculation formula is as follows:
[0021] fp i =fv i *g
[0022] Among them, fp i The wavenumber domain scattered wave imaging result, fv i Let g be the reverse time-shifted image in the wavenumber domain, and g be the scattering Fourier filter operator.
[0023] Secondly, this application proposes a scattered wave stacking imaging device, including: a block module, a Fourier transform module, a filtering module, an inverse Fourier transform module, and a stitching module;
[0024] The segmentation module is used to divide the reverse time-shifted image into regions in the spatial dimension to obtain the reverse time-shifted image of each region.
[0025] The Fourier transform module is used to perform a three-dimensional discrete Fourier transform on the reverse time-shifted image of each region to obtain the reverse time-shifted image in the wavenumber domain.
[0026] The filtering module is used to filter the reverse time-shifted image in the wavenumber domain based on the scattering Fourier filter operator to obtain the wavenumber domain scattered wave imaging result;
[0027] The inverse Fourier transform module is used to perform discrete inverse Fourier transform on the wavenumber domain scattered wave imaging results to obtain spatial domain scattered wave imaging results.
[0028] The stitching module is used to stitch together the spatial domain scattered wave imaging results of each block area to obtain a complete scattered wave stacked imaging result.
[0029] The filtering module includes: a wavenumber spectrum acquisition unit, a scattering angle calculation unit, a filter operator determination unit, and a filtering unit;
[0030] The wavenumber spectrum acquisition unit is used to extract the spatial wavenumber spectrum corresponding to the Fourier transform coordinates from the reverse time-shifted image in the wavenumber domain.
[0031] The scattering angle calculation unit is used to perform an arctangent operation on the spatial wavenumber spectrum corresponding to the Fourier transform coordinates to obtain the scattering angle corresponding to the spatial wavenumber spectrum;
[0032] The filter operator determination unit is used to determine the scattering Fourier filter operator based on the scattering angle corresponding to the spatial wavenumber spectrum;
[0033] The filtering unit is used to filter the reverse time-shifted image in the wavenumber domain according to the scattering Fourier filter operator to obtain the wavenumber domain scattered wave imaging result.
[0034] In the scattering angle calculation unit, the arctangent of the space wavenumber spectrum corresponding to the Fourier transform coordinates is performed using a first calculation formula to obtain the scattering angle corresponding to the space wavenumber spectrum. The first calculation formula is as follows:
[0035]
[0036] Where θ is the scattering angle corresponding to the spatial wavenumber distribution, α is the damping factor, kx is the spatial wavenumber spectrum in the x-direction corresponding to the Fourier transform coordinates, ky is the spatial wavenumber spectrum in the y-direction corresponding to the Fourier transform coordinates, and kz is the spatial wavenumber spectrum in the z-direction corresponding to the Fourier transform coordinates.
[0037] Thirdly, this application proposes an electronic device, comprising: one or more processors, and a memory storing instructions that, when executed by the one or more processors, cause the one or more processors to perform the scattered wave stacking imaging method.
[0038] Fourthly, this application proposes a storage medium storing executable instructions that, when executed, cause a machine to perform the scattered wave stacking imaging method.
[0039] Beneficial effects:
[0040] This invention proposes a method, device, electronic device, and storage medium for post-stack imaging of scattered waves. Based on the idea of regional segmentation, it uses the Fourier filter operator of scattered waves to perform post-stack imaging of scattered waves in the wavenumber domain on the inverse time migration results. This solves the problem of high computational cost in the separation of scattered wave imaging, reduces the amount of computation, lowers the hardware cost, and realizes fine imaging of small-scale slit structures. Attached Figure Description
[0041] Figure 1 This is a flowchart of a method for post-stack imaging of scattered waves according to an embodiment of this application;
[0042] Figure 2 This is a flowchart illustrating the process of obtaining the scattering Fourier filter operator according to an embodiment of this application.
[0043] Figure 3 This is a two-dimensional slice of actual three-dimensional data, offset in reverse time.
[0044] Figure 4 This is a two-dimensional slice of the three-dimensional wavenumber domain scattering Fourier filter operator in an embodiment of this application;
[0045] Figure 5 This is a two-dimensional slice of the complete scattered wave post-stack imaging result of the embodiments of this application;
[0046] Figure 6 This is a schematic diagram of a scattered wave superimposed imaging device according to an embodiment of this application. Detailed Implementation
[0047] The present disclosure will be further described below with reference to the embodiments shown in the accompanying drawings.
[0048] To enable those skilled in the art to better understand the present invention and to fully understand and implement the process of how the present invention uses technical means to solve technical problems and achieve corresponding technical effects, the technical solutions of the present invention will be clearly and completely described below with reference to the accompanying drawings of the embodiments of the present invention. Obviously, the described embodiments are only some embodiments of the present invention, not all embodiments. The embodiments of the present invention and the various features therein can be combined with each other without conflict, and the resulting technical solutions are all within the protection scope of the present invention. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort should fall within the protection scope of the present invention.
[0049] It should be noted that the terms "first," "second," etc., in the specification, claims, and accompanying drawings of this invention are used to distinguish similar objects and are not necessarily used to describe a specific order or sequence. It should be understood that such data can be interchanged where appropriate so that the embodiments of the invention described herein can be implemented in orders other than those illustrated or described herein. Furthermore, the terms "comprising" and "having," and any variations thereof, are intended to cover a non-exclusive inclusion; for example, a process, method, system, product, or apparatus that comprises a series of steps or units is not necessarily limited to those steps or units explicitly listed, but may include other steps or units not explicitly listed or inherent to such processes, methods, products, or apparatus.
[0050] It should be noted that the steps shown in the flowchart in the accompanying drawings can be executed in a computer system such as a set of computer-executable instructions, and although a logical order is shown in the flowchart, in some cases the steps shown or described may be executed in a different order than that shown here.
[0051] Fine imaging of small-scale slots and holes is a core issue in high-precision ultra-deep imaging in Northwest China. Reverse-time migration (RTM) can revert observed reflected energy to the correct imaging location, but information about small-scale anomalies is primarily contained in scattered waves, making amplitude less dominant and difficult to highlight slot anomalies. Scattered wave imaging technology, by fully utilizing scattered signals, has advantages in imaging small-scale anomalies and slot structures. However, conventional scattered wave imaging techniques require wavefield separation during RTM propagation to extract scattered signals for imaging, incurring very high computational costs. To address the high computational cost of pre-stack scattered wave imaging, this application proposes a post-stack scattered wave imaging method, apparatus, electronic equipment, and storage medium. Based on RTM results, a post-stack Fourier filter operator is used to achieve rapid scattered wave imaging for small-scale slot structures. The use of the post-stack Fourier filter operator for rapid scattered wave imaging also improves the identification capability of small-scale slots and holes.
[0052] Example 1:
[0053] This embodiment proposes a method for post-stack imaging of scattered waves, such as... Figure 1 As shown, it includes:
[0054] Step S1: Divide the time-shifted image into regions in the spatial dimension to obtain the time-shifted image of each region; wherein, the time-shifted image can be obtained by the classic time-shifting technique, and the time-shifted image includes: reflected waves and scattered waves.
[0055] Step S2: Perform a three-dimensional discrete Fourier transform on the reverse time-shifted image of each region to obtain the reverse time-shifted image in the wavenumber domain;
[0056] Step S3: Filter the reverse time-shifted image in the wavenumber domain based on the scattering Fourier filter operator to obtain the wavenumber domain scattered wave imaging result; the reflected wave is filtered out by filtering to obtain the scattered wave imaging result.
[0057] Step S4: Perform discrete inverse Fourier transform on the wavenumber domain scattered wave imaging results to obtain the spatial domain scattered wave imaging results;
[0058] Step S5: Stitch together the spatial domain scattered wave imaging results in each region to obtain a complete scattered wave stacked imaging result.
[0059] This embodiment divides the reverse-time migrated image into spatial regions to obtain reverse-time migrated images for each region. A three-dimensional discrete Fourier transform is performed on the reverse-time migrated image of each region to obtain a wavenumber domain reverse-time migrated image. The wavenumber domain reverse-time migrated image is then filtered using a scattering Fourier filter operator to obtain a wavenumber domain scattered wave imaging result. A discrete inverse Fourier transform is then performed on the wavenumber domain scattered wave imaging result to obtain a spatial domain scattered wave imaging result. The spatial domain scattered wave imaging results from each region are then stitched together to obtain a complete scattered wave stacked imaging result. This embodiment achieves rapid scattered wave imaging with minimal computation, reducing hardware resource costs.
[0060] Example 2:
[0061] Based on the foregoing embodiments, the specific implementation methods of the present invention will be further described as follows:
[0062] In some implementations, the process of obtaining the scattering Fourier filter operator can be as follows: Figure 2 As shown, it includes:
[0063] Step S3.1: Extract the spatial wavenumber spectrum corresponding to the Fourier transform coordinates from the reverse time-shifted image in the wavenumber domain. The reverse time-shifted image in the wavenumber domain is the result of a three-dimensional discrete Fourier transform. The spatial wavenumber spectrum corresponding to the Fourier transform coordinates can be directly extracted from the reverse time-shifted image in the wavenumber domain.
[0064] Step S3.2: Perform arctangent operation on the spatial wavenumber spectrum corresponding to the Fourier transform coordinates to obtain the scattering angle corresponding to the spatial wavenumber spectrum;
[0065] The arctangent of the spatial wavenumber spectrum corresponding to the Fourier transform coordinates is performed using the first calculation formula to obtain the scattering angle corresponding to the spatial wavenumber spectrum. The first calculation formula is as follows:
[0066]
[0067] Where θ is the scattering angle corresponding to the spatial wavenumber distribution, α is the damping factor with a value of 0.001, kx is the spatial wavenumber spectrum in the x-direction corresponding to the Fourier transform coordinates, ky is the spatial wavenumber spectrum in the y-direction corresponding to the Fourier transform coordinates, and kz is the spatial wavenumber spectrum in the z-direction corresponding to the Fourier transform coordinates.
[0068] Step S3.3: Determine the scattering Fourier filter operator based on the scattering angle corresponding to the spatial wavenumber spectrum.
[0069] In some implementations, a preset rule is used to determine the scattering Fourier filter operator based on the scattering angle corresponding to the spatial wavenumber spectrum. The preset rule is as follows:
[0070]
[0071] Where g is the scattering Fourier filter operator, θ is the scattering angle corresponding to the space wavenumber spectrum, θ0 is the scattering angle threshold corresponding to the space wavenumber spectrum, and r is the attenuation coefficient.
[0072] when The scattering Fourier filter operator takes a value of 1, while when The scattering Fourier filter operator is based on The calculation yields , where e represents the natural constant.
[0073] In some implementations, a second calculation formula is used to filter the reverse-time-shifted image in the wavenumber domain based on the scattering Fourier filter operator. The second calculation formula is as follows:
[0074] fp i =fv i *g
[0075] Among them, fp i The wavenumber domain scattered wave imaging result, fv i Let g be the reverse time-shifted image in the wavenumber domain, and g be the scattering Fourier filter operator.
[0076] The proposed scattered wave stacking imaging method in this embodiment divides the reverse-time migrated image into spatial regions. For each region, a three-dimensional discrete Fourier transform is performed on the reverse-time migrated image to obtain a wavenumber domain reverse-time migrated image. The wavenumber domain reverse-time migrated image is then filtered using a scattering Fourier filter operator to obtain a wavenumber domain scattered wave imaging result. A discrete inverse Fourier transform is then performed on the wavenumber domain scattered wave imaging result to obtain a spatial domain scattered wave imaging result. Finally, the spatial domain scattered wave imaging results from each region are stitched together to obtain a complete scattered wave stacking imaging result. This method achieves rapid scattered wave imaging with minimal computation, reducing hardware resource costs.
[0077] This embodiment provides a specific implementation process for the scattering Fourier filter operator. The spatial wavenumber spectrum corresponding to the Fourier transform coordinates is extracted from the reverse-time offset image in the wavenumber domain. An arctangent operation is performed on the spatial wavenumber spectrum corresponding to the Fourier transform coordinates to obtain the scattering angle corresponding to the spatial wavenumber spectrum. The scattering Fourier filter operator is determined based on the angle of the scattering angle corresponding to the spatial wavenumber spectrum. The effectiveness of the filter determines the final imaging effect of this embodiment. The preset rules provided in this embodiment can accurately determine the corresponding scattering Fourier filter operator based on the angular range of the scattering angle corresponding to the spatial wavenumber spectrum. Then, the reverse-time offset image in the wavenumber domain is filtered based on the scattering Fourier filter operator to obtain the wavenumber domain scattered wave imaging result. The spatial domain scattered wave imaging results in each region are stitched together to obtain a relatively accurate complete scattered wave stacked imaging result.
[0078] Example 3:
[0079] This embodiment provides an application example of a scattered wave stacking imaging method:
[0080] (1) Input the counter-time offset image v, such as Figure 3 As shown, the horizontal axis represents horizontal distance, and the vertical axis represents depth. It can be seen that there are a large number of slits, most of which are isolated between two strong reflective layers, with some embedded at the bottom boundary of the reflective layers, resulting in insufficient recognition. The reverse-time migrated image v is divided into regions along the X, Y, and Z dimensions, for example, 5×5×5, resulting in 125 regions, yielding the reverse-time migrated image v for each region. i , where v i This is the reverse-time offset image of the i-th region.
[0081] (2) Perform a three-dimensional discrete Fourier transform on the reverse time-shifted image of each region to obtain the reverse time-shifted image fv in the wavenumber domain. i As the filtering result of the three-dimensional discrete Fourier transform, fv i Let be the wavenumber domain of the i-th region, which is a time-shifted image.
[0082] The arctangent of the spatial wavenumber spectrum corresponding to the Fourier transform coordinates is performed using the first calculation formula to obtain the scattering angle corresponding to the spatial wavenumber spectrum. The first calculation formula is as follows:
[0083]
[0084] Where θ is the scattering angle corresponding to the spatial wavenumber distribution, α is the damping factor with a value of 0.001, kx is the spatial wavenumber spectrum in the x-direction corresponding to the Fourier transform coordinates, ky is the spatial wavenumber spectrum in the y-direction corresponding to the Fourier transform coordinates, and kz is the spatial wavenumber spectrum in the z-direction corresponding to the Fourier transform coordinates.
[0085] The scattering Fourier filter operator is determined based on the scattering angle corresponding to the spatial wavenumber spectrum using preset rules, as follows:
[0086]
[0087] Where g is the scattering Fourier filter operator, θ is the scattering angle corresponding to the space wavenumber spectrum, θ0 is the scattering angle threshold corresponding to the space wavenumber spectrum, the scattering angle threshold corresponding to the space wavenumber spectrum is in the range of [0, 90] degrees, and r is the attenuation coefficient.
[0088] (3) The inverse time-shifted image fv in the wavenumber domain based on the scattering Fourier filter operator g. i Filtering is performed to obtain the wavenumber domain scattered wave imaging result fp. i ;fp i This represents the wavenumber domain scattered wave imaging result for the i-th region. The scattering Fourier filter operator g is as follows: Figure 4 As shown, the horizontal axis represents the horizontal wavenumber, and the vertical axis represents the vertical wavenumber.
[0089] (4) Based on the wavenumber domain scattered wave imaging results fp i Perform an inverse Fourier transform to obtain the spatial domain scattered wave imaging result p. i ;
[0090] The second calculation formula is used to filter the reverse time-shifted image in the wavenumber domain based on the scattering Fourier filter operator. The second calculation formula is as follows:
[0091] fp i =fv i *g
[0092] Among them, fp i The wavenumber domain scattered wave imaging result, fv i Let g be the reverse time-shifted image in the wavenumber domain, and g be the scattering Fourier filter operator.
[0093] (5) The spatial domain scattered wave imaging results within all blocks are stitched together to obtain the complete scattered wave post-stack imaging result p, such as Figure 5 As shown, the horizontal axis represents horizontal distance, and the vertical axis represents depth. It can be seen that the energy of the strong reflective layer has been largely filtered out, highlighting the fractured-vuggy reservoir and providing a clearer beaded characterization effect. This is of great significance for analyzing reservoir size and properties, and assisting in reservoir interpretation and development.
[0094] In this embodiment, a time-shifted image v is first input, and then the time-shifted image v is divided into regions along the X, Y, and Z spatial dimensions, for example, 5x5x5, resulting in 125 regions, thus obtaining the time-shifted image v for each region. i , where v iLet fv be the reverse-time migrated image of the i-th region. Perform a three-dimensional discrete Fourier transform on the reverse-time migrated image of each region to obtain the wavenumber domain reverse-time migrated image fv. i This serves as the filtering result of the three-dimensional discrete Fourier transform. The inverse time-shifted image fv in the wavenumber domain is obtained based on the scattering Fourier filter operator g. i Filtering is performed to obtain the wavenumber domain scattered wave imaging result fp. i The scattering Fourier filter operator g is based on the wavenumber domain scattered wave imaging result fp. i Perform an inverse Fourier transform to obtain the spatial domain scattered wave imaging result p. i The spatial domain scattered wave imaging results within all blocks are stitched together to obtain a complete scattered wave stacked imaging result p.
[0095] A time-shifted image is obtained using the classic time-shifting technique. Without the method described in this application, a two-dimensional time-shifted slice of actual three-dimensional data is provided, such as... Figure 3 As shown, from Figure 3 As can be seen, there are a large number of slits, most of which are isolated between two strong reflective layers, and some are embedded in the bottom boundary of the reflective layers, resulting in insufficient identification. Using the method of this application, a complete post-stack imaging result of the scattered waves is finally obtained, as shown below. Figure 5 As shown, the reflected wave of the reverse time-shifted image, i.e., the energy of the strong reflection layer, has been basically filtered out, highlighting the fractured-vuggy reservoir and providing a clearer beaded characterization effect. This is of great significance for analyzing reservoir size and properties and assisting in reservoir interpretation and development.
[0096] Example 4:
[0097] This embodiment proposes a scattered wave stacking imaging device, such as... Figure 6 As shown, it includes: a block module, a Fourier transform module, a filtering module, an inverse Fourier transform module, and a splicing module;
[0098] The block module, Fourier transform module, filtering module, inverse Fourier transform module, and splicing module are connected in sequence.
[0099] The segmentation module is used to divide the reverse time-shifted image into regions in the spatial dimension to obtain the reverse time-shifted image of each region.
[0100] The Fourier transform module is used to perform a three-dimensional discrete Fourier transform on the reverse time-shifted image of each region to obtain the reverse time-shifted image in the wavenumber domain.
[0101] The filtering module is used to filter the reverse time-shifted image in the wavenumber domain based on the scattering Fourier filter operator to obtain the wavenumber domain scattered wave imaging result;
[0102] The inverse Fourier transform module is used to perform discrete inverse Fourier transform on the wavenumber domain scattered wave imaging results to obtain spatial domain scattered wave imaging results.
[0103] The stitching module is used to stitch together the spatial domain scattered wave imaging results of each block area to obtain a complete scattered wave stacked imaging result.
[0104] The filtering module includes: a wavenumber spectrum acquisition unit, a scattering angle calculation unit, a filter operator determination unit, and a filtering unit;
[0105] The wavenumber spectrum acquisition unit is used to extract the spatial wavenumber spectrum corresponding to the Fourier transform coordinates from the reverse time-shifted image in the wavenumber domain.
[0106] The scattering angle calculation unit is used to perform an arctangent operation on the spatial wavenumber spectrum corresponding to the Fourier transform coordinates to obtain the scattering angle corresponding to the spatial wavenumber spectrum;
[0107] The filter operator determination unit is used to determine the scattering Fourier filter operator based on the scattering angle corresponding to the spatial wavenumber spectrum;
[0108] The filtering unit is used to filter the reverse time-shifted image in the wavenumber domain according to the scattering Fourier filter operator to obtain the wavenumber domain scattered wave imaging result.
[0109] In the scattering angle calculation unit, the arctangent of the space wavenumber spectrum corresponding to the Fourier transform coordinates is performed using a first calculation formula to obtain the scattering angle corresponding to the space wavenumber spectrum. The first calculation formula is as follows:
[0110]
[0111] Where θ is the scattering angle corresponding to the spatial wavenumber distribution, α is the damping factor, kx is the spatial wavenumber spectrum in the x-direction corresponding to the Fourier transform coordinates, ky is the spatial wavenumber spectrum in the y-direction corresponding to the Fourier transform coordinates, and kz is the spatial wavenumber spectrum in the z-direction corresponding to the Fourier transform coordinates.
[0112] In the filter operator determination unit, a preset rule is used to determine the scattering Fourier filter operator based on the scattering angle corresponding to the spatial wavenumber spectrum. The preset rule is as follows:
[0113]
[0114] Where g is the scattering Fourier filter operator, θ is the scattering angle corresponding to the space wavenumber spectrum, θ0 is the scattering angle threshold corresponding to the space wavenumber spectrum, and r is the attenuation coefficient.
[0115] In the filtering module, a second calculation formula is used to filter the reverse-time offset image in the wavenumber domain based on the scattering Fourier filter operator. The second calculation formula is as follows:
[0116] fp i =fv i *g
[0117] Among them, fp i The wavenumber domain scattered wave imaging result, fv i Let g be the reverse time-shifted image in the wavenumber domain, and g be the scattering Fourier filter operator.
[0118] This embodiment proposes a post-stack imaging device for scattered waves, comprising: a segmentation module, a Fourier transform module, a filtering module, an inverse Fourier transform module, and a stitching module. Each module is sequentially connected. The segmentation module divides the reverse-time migrated image into spatial regions, obtaining reverse-time migrated images for each region. The Fourier transform module performs a three-dimensional discrete Fourier transform on the reverse-time migrated image of each region, obtaining a wavenumber domain reverse-time migrated image. The filtering module filters the wavenumber domain reverse-time migrated image based on a scattering Fourier filter operator, obtaining a wavenumber domain scattered wave imaging result. The inverse Fourier transform module performs a discrete inverse Fourier transform on the wavenumber domain scattered wave imaging result, obtaining a spatial domain scattered wave imaging result. The stitching module stitches the spatial domain scattered wave imaging results from each region, obtaining a complete post-stack imaging result for scattered waves. This embodiment achieves rapid imaging of scattered waves for small-scale slot structures by using a post-stack Fourier filter operator on top of the reverse-time migrated image. Rapid imaging of scattered waves is achieved by using a post-stack Fourier filter operator, which also improves the ability to identify small-scale pores.
[0119] Example 5:
[0120] Another embodiment of this application provides an electronic device, including: one or more processors, and a memory storing instructions that, when executed by the one or more processors, cause the one or more processors to perform the scattered wave stacking imaging method of the foregoing embodiments.
[0121] The scattered wave stacking imaging method includes:
[0122] Step S1: Divide the reverse time-shifted image into regions in the spatial dimension to obtain the reverse time-shifted image of each region;
[0123] Step S2: Perform a three-dimensional discrete Fourier transform on the reverse time-shifted image of each region to obtain the reverse time-shifted image in the wavenumber domain;
[0124] Step S3: Filter the reverse time-shifted image in the wavenumber domain based on the scattering Fourier filter operator to obtain the wavenumber domain scattered wave imaging result;
[0125] The process of obtaining the scattering Fourier filter operator is as follows: Figure 2 As shown:
[0126] Step S3.1: Extract the spatial wavenumber spectrum corresponding to the Fourier transform coordinates from the reverse time-shifted image in the wavenumber domain. The reverse time-shifted image in the wavenumber domain is the result of a three-dimensional discrete Fourier transform. The spatial wavenumber spectrum corresponding to the Fourier transform coordinates can be directly extracted from the reverse time-shifted image in the wavenumber domain.
[0127] Step S3.2: Perform arctangent operation on the spatial wavenumber spectrum corresponding to the Fourier transform coordinates to obtain the scattering angle corresponding to the spatial wavenumber spectrum;
[0128] The arctangent of the spatial wavenumber spectrum corresponding to the Fourier transform coordinates is performed using the first calculation formula to obtain the scattering angle corresponding to the spatial wavenumber spectrum. The first calculation formula is as follows:
[0129]
[0130] Where θ is the scattering angle corresponding to the spatial wavenumber distribution, α is the damping factor with a value of 0.001, kx is the spatial wavenumber spectrum in the x-direction corresponding to the Fourier transform coordinates, ky is the spatial wavenumber spectrum in the y-direction corresponding to the Fourier transform coordinates, and kz is the spatial wavenumber spectrum in the z-direction corresponding to the Fourier transform coordinates.
[0131] Step S3.3: Determine the scattering Fourier filter operator based on the scattering angle corresponding to the spatial wavenumber spectrum.
[0132] The scattering Fourier filter operator is determined based on the scattering angle corresponding to the spatial wavenumber spectrum using preset rules, as follows:
[0133]
[0134] Where g is the scattering Fourier filter operator, θ is the scattering angle corresponding to the space wavenumber spectrum, θ0 is the scattering angle threshold corresponding to the space wavenumber spectrum, and r is the attenuation coefficient.
[0135] Step S4: Perform discrete inverse Fourier transform on the wavenumber domain scattered wave imaging results to obtain the spatial domain scattered wave imaging results;
[0136] The second calculation formula is used to filter the reverse time-shifted image in the wavenumber domain based on the scattering Fourier filter operator. The second calculation formula is as follows:
[0137] fp i =fv i*g
[0138] Among them, fp i The wavenumber domain scattered wave imaging result, fv i Let g be the reverse time-shifted image in the wavenumber domain, and g be the scattering Fourier filter operator.
[0139] Step S5: Stitch together the spatial domain scattered wave imaging results in each region to obtain a complete scattered wave stacked imaging result.
[0140] The electronic device may be a mobile phone, computer, or tablet computer, etc., and includes a memory and a processor. The memory stores a computer program, which, when executed by the processor, implements the scattered wave post-stack imaging method as described in the embodiments. It is understood that the electronic device may also include an input / output (I / O) interface and communication components.
[0141] The processor is used to execute all or part of the steps in the scattered wave stacking imaging method as described in the above embodiments. The memory is used to store various types of data, which may include, for example, instructions for any application or method in the electronic device, as well as application-related data.
[0142] The processor can be implemented as an application-specific integrated circuit (ASIC), a digital signal processor (DSP), a programmable logic device (PLD), a field-programmable gate array (FPGA), a controller, a microcontroller, a microprocessor, or other electronic components, and is used to execute the scattered wave stacking imaging method in the above embodiments.
[0143] Example 6:
[0144] This embodiment proposes a computer-readable storage medium storing executable instructions that, when executed, cause a machine to perform the scattered wave stacking imaging method described in the foregoing embodiment.
[0145] The scattered wave stacking imaging method includes:
[0146] Step S1: Divide the reverse time-shifted image into regions in the spatial dimension to obtain the reverse time-shifted image of each region;
[0147] Step S2: Perform a three-dimensional discrete Fourier transform on the reverse time-shifted image of each region to obtain the reverse time-shifted image in the wavenumber domain;
[0148] Step S3: Filter the reverse time-shifted image in the wavenumber domain based on the scattering Fourier filter operator to obtain the wavenumber domain scattered wave imaging result;
[0149] The process of obtaining the scattering Fourier filter operator is as follows: Figure 2 As shown:
[0150] Step S3.1: Extract the spatial wavenumber spectrum corresponding to the Fourier transform coordinates from the reverse time-shifted image in the wavenumber domain. The reverse time-shifted image in the wavenumber domain is the result of a three-dimensional discrete Fourier transform. The spatial wavenumber spectrum corresponding to the Fourier transform coordinates can be directly extracted from the reverse time-shifted image in the wavenumber domain.
[0151] Step S3.2: Perform arctangent operation on the spatial wavenumber spectrum corresponding to the Fourier transform coordinates to obtain the scattering angle corresponding to the spatial wavenumber spectrum;
[0152] The arctangent of the spatial wavenumber spectrum corresponding to the Fourier transform coordinates is performed using the first calculation formula to obtain the scattering angle corresponding to the spatial wavenumber spectrum. The first calculation formula is as follows:
[0153]
[0154] Where θ is the scattering angle corresponding to the spatial wavenumber distribution, α is the damping factor with a value of 0.001, kx is the spatial wavenumber spectrum in the x-direction corresponding to the Fourier transform coordinates, ky is the spatial wavenumber spectrum in the y-direction corresponding to the Fourier transform coordinates, and kz is the spatial wavenumber spectrum in the z-direction corresponding to the Fourier transform coordinates.
[0155] Step S3.3: Determine the scattering Fourier filter operator based on the scattering angle corresponding to the spatial wavenumber spectrum.
[0156] The scattering Fourier filter operator is determined based on the scattering angle corresponding to the spatial wavenumber spectrum using preset rules, as follows:
[0157]
[0158] Where g is the scattering Fourier filter operator, θ is the scattering angle corresponding to the space wavenumber spectrum, θ0 is the scattering angle threshold corresponding to the space wavenumber spectrum, and r is the attenuation coefficient.
[0159] Step S4: Perform discrete inverse Fourier transform on the wavenumber domain scattered wave imaging results to obtain the spatial domain scattered wave imaging results;
[0160] The second calculation formula is used to filter the reverse time-shifted image in the wavenumber domain based on the scattering Fourier filter operator. The second calculation formula is as follows:
[0161] fp i =fv i*g
[0162] Among them, fp i The wavenumber domain scattered wave imaging result, fv i Let g be the reverse time-shifted image in the wavenumber domain, and g be the scattering Fourier filter operator.
[0163] Step S5: Stitch together the spatial domain scattered wave imaging results in each region to obtain a complete scattered wave stacked imaging result.
[0164] In the various embodiments of this application, the functional units can be integrated into one processing unit, or each unit can exist physically separately, or two or more units can be integrated into one unit. If the functions are implemented as software functional units and sold or used as independent products, they can be stored in a computer-readable storage medium.
[0165] Based on this understanding, the technical solution of this application, or the part that contributes to the prior art, or part of the technical solution, can be embodied in the form of a software product. The computer software product is stored in a storage medium and includes several instructions to cause a computer device (which may be a personal computer, a server, or a network device, etc.) to execute all or part of the steps of the scattered wave-based post-stack imaging method described in the various embodiments of this invention.
[0166] The aforementioned storage media include: flash memory, hard disk, multimedia card, card-type memory (e.g., SD (Secure Digital Memory Card) or DX (Memory Data Register, MDR) memory, random access memory (RAM), static random access memory (SRAM), read-only memory (ROM), electrically erasable programmable read-only memory (EEPROM), programmable read-only memory (PROM), magnetic memory, disk, optical disk, server, APP (Application) application store, and other media that can store program verification codes, on which computer programs are stored. When the computer programs are executed by the processor, they can implement the various steps of the above-mentioned scattered wave superimposed imaging method.
[0167] Those skilled in the art will recognize that the units and algorithm steps of the various examples described in conjunction with the embodiments disclosed herein can be implemented in electronic hardware, or a combination of computer software and electronic hardware. Whether these functions are implemented in hardware or software depends on the specific application and design constraints of the technical solution. Those skilled in the art can use different methods to implement the described functions for each specific application, but such implementation should not be considered beyond the scope of this invention. It will be clearly understood by those skilled in the art that, for the sake of convenience and brevity, the specific working processes of the systems, devices, and units described above can be referred to the corresponding processes in the foregoing method embodiments, and will not be repeated here.
[0168] In the several embodiments provided in this application, it should be understood that the disclosed systems, apparatuses, and methods can be implemented in other ways. For example, the apparatus embodiments described above are merely illustrative. For instance, the division of units is only a logical functional division, and in actual implementation, there may be other division methods. For example, multiple units or components may be combined or integrated into another system, or some features may be ignored or not executed. Furthermore, the shown or discussed mutual couplings, direct couplings, or communication connections may be through some interfaces; indirect couplings or communication connections between devices or units may be electrical, mechanical, or other forms. The units described as separate components may or may not be physically separate. The components shown as units may or may not be physical units, i.e., they may be located in one place or distributed across multiple network units. Some or all of the units can be selected to achieve the purpose of this embodiment according to actual needs.
[0169] In the description of this application, "multiple" means at least two, such as two, three, etc., unless otherwise explicitly specified. All directional indications (such as up, down, left, right, front, back, top, bottom, etc.) in the embodiments of this application are only used to explain the relative positional relationships and movement of the components in a specific posture (as shown in the figures). If the specific posture changes, the directional indication will also change accordingly. Furthermore, the terms "comprising" and "having," and any variations thereof, are intended to cover non-exclusive inclusion. For example, a process, method, system, product, or device that includes a series of steps or units is not limited to the listed steps or units, but may optionally include steps or units not listed, or may optionally include other steps or units inherent to these processes, methods, products, or devices.
[0170] Furthermore, the reference to "embodiment" herein means that a particular feature, structure, or characteristic described in connection with an embodiment may be included in at least one embodiment of this application. The appearance of this phrase in various places throughout the specification does not necessarily refer to the same embodiment, nor is it a separate or alternative embodiment mutually exclusive with other embodiments. It will be explicitly and implicitly understood by those skilled in the art that the embodiments described herein can be combined with other embodiments.
[0171] The above description is merely a specific embodiment of the present invention, but the scope of protection of the present invention is not limited thereto. Any variations or substitutions that can be easily conceived by those skilled in the art within the technical scope disclosed in the present invention should be included within the scope of protection of the present invention. Therefore, the scope of protection of the present invention should be determined by the scope of the claims. The above description is merely a preferred embodiment of the present invention and is not intended to limit the present invention. Any modifications, equivalent substitutions, etc., made within the spirit and principles of the present invention should be included within the scope of protection of the present invention.
Claims
1. A method for imaging after scattering wave stacking, characterized in that, include: The reverse time-shifted image is divided into regions in the spatial dimension to obtain the reverse time-shifted image of each region. A three-dimensional discrete Fourier transform is performed on the reverse time-shifted image of each region to obtain the reverse time-shifted image in the wavenumber domain. The wavenumber domain reverse-time-shifted image is filtered based on the scattering Fourier filter operator to obtain the wavenumber domain scattered wave imaging result; The wavenumber domain scattered wave imaging results are subjected to discrete inverse Fourier transform to obtain the spatial domain scattered wave imaging results; The spatial domain scattered wave imaging results of each region are stitched together to obtain a complete scattered wave stacked imaging result. The process of obtaining the scattering Fourier filter operator is as follows: Extract the spatial wavenumber spectrum corresponding to the Fourier transform coordinates from the reverse time-shifted image in the wavenumber domain; Perform an arctangent operation on the spatial wavenumber spectrum corresponding to the Fourier transform coordinates to obtain the scattering angle corresponding to the spatial wavenumber spectrum; The scattering Fourier filter operator is determined based on the scattering angle corresponding to the spatial wavenumber spectrum.
2. The scattered wave stacking imaging method as described in claim 1, characterized in that, The arctangent of the spatial wavenumber spectrum corresponding to the Fourier transform coordinates is performed using the first calculation formula to obtain the scattering angle corresponding to the spatial wavenumber spectrum. The first calculation formula is as follows: in, It is the damping factor. The wavenumber spectrum in the spatial x-direction corresponding to the Fourier transform coordinates. The wavenumber spectrum in the spatial y-direction corresponding to the Fourier transform coordinates. This represents the wavenumber spectrum in the spatial z-direction corresponding to the Fourier transform coordinates.
3. The scattered wave stacking imaging method as described in claim 1, characterized in that, The scattering Fourier filter operator is determined based on the scattering angle corresponding to the spatial wavenumber spectrum using preset rules, as follows: Where g is the scattering Fourier filter operator. for Threshold, r is the attenuation coefficient.
4. The scattered wave stacking imaging method as described in claim 1, characterized in that, The second calculation formula is used to filter the reverse time-shifted image in the wavenumber domain based on the scattering Fourier filter operator. The second calculation formula is as follows: in, The results are from wavenumber domain scattered wave imaging. This is a reverse-time migration image in the wavenumber domain. This is the scattering Fourier filter operator.
5. A device for imaging after scattering wave stacking, characterized in that, include: Blocking module, Fourier transform module, filtering module, inverse Fourier transform module, splicing module; The segmentation module is used to divide the reverse time-shifted image into regions in the spatial dimension to obtain the reverse time-shifted image of each region. The Fourier transform module is used to perform a three-dimensional discrete Fourier transform on the reverse time-shifted image of each region to obtain the reverse time-shifted image in the wavenumber domain. The filtering module is used to filter the reverse time-shifted image in the wavenumber domain based on the scattering Fourier filter operator to obtain the wavenumber domain scattered wave imaging result; The inverse Fourier transform module is used to perform discrete inverse Fourier transform on the wavenumber domain scattered wave imaging results to obtain spatial domain scattered wave imaging results. The stitching module is used to stitch together the spatial domain scattered wave imaging results of each block area to obtain a complete scattered wave stacked imaging result. The filtering module includes: a wavenumber spectrum acquisition unit, a scattering angle calculation unit, a filter operator determination unit, and a filtering unit; The wavenumber spectrum acquisition unit is used to extract the spatial wavenumber spectrum corresponding to the Fourier transform coordinates from the reverse time-shifted image in the wavenumber domain. The scattering angle calculation unit is used to perform an arctangent operation on the spatial wavenumber spectrum corresponding to the Fourier transform coordinates to obtain the scattering angle corresponding to the spatial wavenumber spectrum; The filter operator determination unit is used to determine the scattering Fourier filter operator based on the scattering angle corresponding to the spatial wavenumber spectrum; The filtering unit is used to filter the reverse time-shifted image in the wavenumber domain according to the scattering Fourier filter operator to obtain the wavenumber domain scattered wave imaging result.
6. The scattered wave superposition imaging device as described in claim 5, characterized in that, In the scattering angle calculation unit, the arctangent of the space wavenumber spectrum corresponding to the Fourier transform coordinates is performed using a first calculation formula to obtain the scattering angle corresponding to the space wavenumber spectrum. The first calculation formula is as follows: in, It is the damping factor. The wavenumber spectrum in the spatial x-direction corresponding to the Fourier transform coordinates. The wavenumber spectrum in the spatial y-direction corresponding to the Fourier transform coordinates. This represents the wavenumber spectrum in the spatial z-direction corresponding to the Fourier transform coordinates.
7. An electronic device, characterized in that, include: One or more processors, and a memory storing instructions that, when executed by the one or more processors, cause the one or more processors to perform the scattered wave stacking imaging method according to any one of claims 1 to 4.
8. A computer-readable storage medium, characterized in that, It stores executable instructions that, when executed, cause the machine to perform the scattered wave stacking imaging method according to any one of claims 1 to 4.