Scattered wave high-precision imaging method and device, electronic equipment and storage medium

Abstract

The invention discloses a scattered wave high-precision imaging method and device, electronic equipment and a storage medium, and the method comprises the steps: obtaining a least square imaging result, carrying out the processing of the least square imaging result based on a pre-constructed scattering Fourier filtering operator, determining a scattered wave imaging result, and outputting the scattered wave imaging result. And scanning a preset scattered wave target area based on the scattered wave imaging result to obtain an energy compensation factor, and performing weighted stacking on the scattered wave imaging result and the least square imaging result based on the energy compensation factor to generate a scattered wave high-precision imaging result. Scattering energy is separated by utilizing a post-stack scattering wave imaging technology, scattering wave enhanced imaging is realized by utilizing an energy compensation factor, and the high-resolution and high-display imaging effect of the fracture-cavity reservoir is improved.

Description

Technical Field

[0001] The present application relates to the technical field of seismic wave imaging, and in particular to a scattered wave high-precision imaging method, device, electronic equipment and storage medium. Background Art

[0002] At present, the least squares migration imaging technology can effectively eliminate the problems caused by seismic observation, wavelet band limitation, and spatial stretching by utilizing the inverse Hessian, and is an important technical means for high-precision seismic imaging.

[0003] However, among the above technologies, for the fine imaging problem of fracture-vuggy reservoirs, scattered wave imaging can effectively solve the high-quality imaging of small-scale anomalies by stripping off the reflected energy in the seismic wave field. However, separation imaging strips off the reflected imaging information, which is not conducive to the subsequent reservoir inversion interpretation.

[0004] Therefore, how to improve the resolution and display of small-scale abnormal body imaging is a problem that needs to be solved. Summary of the invention

[0005] In view of the above problems, the present application provides a scattered wave high-precision imaging method, device, electronic device and storage medium to at least solve the problems existing in the related art.

[0006] In a first aspect, an embodiment of the present application provides a scattered wave high-precision imaging method, comprising:

[0007] Obtain least squares imaging results;

[0008] Processing the least squares imaging result based on a pre-constructed scattering Fourier filter operator to determine a scattered wave imaging result;

[0009] Scanning a preset scattered wave target area based on the scattered wave imaging result to obtain an energy compensation factor;

[0010] The scattered wave imaging result and the least square imaging result are weightedly superimposed based on the energy compensation factor to generate a scattered wave high-precision imaging result.

[0011] In some embodiments, obtaining the least squares imaging result includes:

[0012] Obtain reverse time migration images and global space-variant scattering functions from seismic data;

[0013] The least squares imaging result is determined by performing high-dimensional spatial deconvolution processing based on the reverse time migration image and the global space-variant scattering function.

[0014] In some embodiments, the energy compensation factor is calculated as:

[0015]

[0016] Where U is the mean amplitude of the scattered wave imaging result, U0 is the mean amplitude of the least squares imaging result, rms is the mean square error, sigmoid is the activation function, θ s is the opening angle, θ r is the azimuth.

[0017] In some embodiments, the scattered wave imaging result and the least square imaging result are weightedly superimposed based on the energy compensation factor to generate a scattered wave high-precision imaging result, and the calculation formula used is:

[0018] m enh =m lsm +F·m svm ;

[0019] In the formula, m lsm is the least squares imaging result, m svm is the scattered wave imaging result, m enh This is the high-precision imaging result of the scattered wave.

[0020] In some embodiments, the calculation formula of the pre-constructed scattering Fourier filter operator is:

[0021]

[0022] Where θ is the scattering angle corresponding to the spatial wave number distribution, and its expression is α is the damping factor, k is the spatial wave number, x, y, z are the spatial coordinate directions, θ0 is the scattering angle of the imaging point, and r is the attenuation smoothing function.

[0023] In some embodiments, the processing of the least square imaging result based on a pre-constructed scattering Fourier filter operator to determine the scattered wave imaging result includes:

[0024] Performing filtering processing on the least squares imaging result based on the scattering Fourier filter operator to obtain a wavenumber domain scattering wave imaging result;

[0025] The scattered wave imaging results in the wavenumber domain are spliced ​​and processed to determine the scattered wave imaging results.

[0026] In a second aspect, an embodiment of the present application provides a scattered wave high-precision imaging method and device, comprising:

[0027] A first acquisition module, used for acquiring least squares imaging results;

[0028] A determination module, configured to determine a scattered wave imaging result by processing the least square imaging result based on a pre-constructed scattering Fourier filter operator;

[0029] A second acquisition module, configured to scan a preset scattered wave target area based on the scattered wave imaging result to acquire an energy compensation factor;

[0030] A generating module is used to perform weighted superposition on the scattered wave imaging result and the least square imaging result based on the energy compensation factor to generate a scattered wave high-precision imaging result.

[0031] In some embodiments, the first acquisition module further includes:

[0032] Sub-acquisition module, used to obtain reverse time migration image and global space-variant scattering function in seismic data;

[0033] The sub-confirmation module is used to determine the least squares imaging result by performing high-dimensional spatial deconvolution processing based on the reverse time migration image and the global space-variant scattering function.

[0034] In a third aspect, an embodiment of the present application provides an electronic device, comprising: at least one processor and a memory; the processor is used to execute a computer program stored in the memory to implement the scattered wave high-precision imaging method as described in any implementation of the first aspect.

[0035] In a fourth aspect, an embodiment of the present application provides a computer storage medium storing one or more programs, which can be executed by an electronic device as described in the third aspect to implement a high-precision imaging method for scattered waves as described in any embodiment of the first aspect.

[0036] The embodiments of the present application provide a scattered wave high-precision imaging method, device, electronic device and storage medium, comprising: obtaining a least squares imaging result, processing the least squares imaging result based on a pre-constructed scattering Fourier filter operator to determine the scattered wave imaging result, scanning a preset scattering wave target area based on the scattered wave imaging result to obtain an energy compensation factor, weightedly superimposing the scattered wave imaging result and the least squares imaging result based on the energy compensation factor to generate a scattered wave high-precision imaging result, separating the scattered energy by using post-stack scattered wave imaging technology under the least squares inversion imaging framework, and realizing scattered wave enhanced imaging by using the energy compensation factor, thereby improving the high-resolution and high-visibility imaging effect of fracture-cavity reservoirs.

[0037] It should be understood that the content described in this section is not intended to identify the key or important features of the embodiments of the present application, nor is it intended to limit the scope of the present application. Other features of the present application will become easily understood through the following description. BRIEF DESCRIPTION OF THE DRAWINGS

[0038] In order to more clearly illustrate the technical solutions in the embodiments of the present application, the drawings required for use in the description of the embodiments will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present application. For those skilled in the art, other drawings can be obtained based on these drawings without creative work.

[0039] Figure 1 A schematic diagram of a scattered wave high-precision imaging method proposed in an embodiment of the present application is shown;

[0040] Figure 2 A schematic diagram of a least squares scattered wave enhanced imaging process proposed in an embodiment of the present application is shown;

[0041] Figure 3 A schematic diagram of a two-dimensional slice of a global space-varying scattering function of three-dimensional actual data proposed in one embodiment of the present application is shown;

[0042] Figure 4 A schematic diagram of a two-dimensional slice of a three-dimensional actual data scattered wave least square imaging proposed in an embodiment of the present application is shown;

[0043] Figure 5 It shows another schematic diagram of a two-dimensional slice of least square imaging of three-dimensional actual data scattered waves proposed in an embodiment of the present application;

[0044] Figure 6 A schematic diagram of a two-dimensional slice of least squares scattered wave enhanced imaging of three-dimensional actual data proposed in one embodiment of the present application is shown;

[0045] Figure 7 A structural block diagram of a scattered wave high-precision imaging method device proposed in an embodiment of the present application is shown;

[0046] Figure 8 A structural block diagram of an electronic device for executing a scattered wave high-precision imaging method according to an embodiment of the present application is shown;

[0047] Fig. 9 A computer-readable storage medium for storing or carrying a scattered wave high-precision imaging method according to an embodiment of the present application is shown. DETAILED DESCRIPTION

[0048] In order to make the objectives, technical solutions and advantages of the present invention more clearly understood, the present invention is further described in detail below in conjunction with embodiments and drawings. The exemplary embodiments of the present invention and their description are only used to explain the present invention and are not intended to limit the present invention.

[0049] Embodiments of the present disclosure will be described in more detail below with reference to the accompanying drawings. Although certain embodiments of the present disclosure are shown in the accompanying drawings, it should be understood that the present disclosure can be implemented in various forms and should not be construed as being limited to the embodiments described herein, which are instead provided for a more thorough and complete understanding of the present disclosure. It should be understood that the drawings and embodiments of the present disclosure are only for exemplary purposes and are not intended to limit the scope of protection of the present disclosure.

[0050] It should be understood that the various steps described in the method embodiments of the present disclosure may be performed in different orders and / or in parallel. In addition, the method embodiments may include additional steps and / or omit the steps shown. The scope of the present disclosure is not limited in this respect.

[0051] The term "including" and its variations used herein are open inclusions, i.e., "including but not limited to". The term "based on" means "based at least in part on". The term "one embodiment" means "at least one embodiment"; the term "another embodiment" means "at least one additional embodiment"; the term "some embodiments" means "at least some embodiments". The relevant definitions of other terms will be given in the following description.

[0052] It should be noted that the concepts such as "first" and "second" mentioned in the present disclosure are only used to distinguish different devices, modules or units, and are not used to limit the order or interdependence of the functions performed by these devices, modules or units.

[0053] It should be noted that the modifications of "one" and "plurality" mentioned in the present disclosure are illustrative rather than restrictive, and those skilled in the art should understand that unless otherwise clearly indicated in the context, it should be understood as "one or more".

[0054] The names of the messages or information exchanged between multiple devices in the embodiments of the present disclosure are only used for illustrative purposes and are not used to limit the scope of these messages or information.

[0055] It is understandable that before using the technical solutions disclosed in the embodiments of the present disclosure, the types, scope of use, usage scenarios, etc. of the personal information involved in the present disclosure should be informed to the user and the user's authorization should be obtained in an appropriate manner in accordance with relevant laws and regulations.

[0056] For example, in response to receiving an active request from a user, a prompt message is sent to the user to clearly prompt the user that the operation requested to be performed will require obtaining and using the user's personal information. Thus, the user can autonomously choose whether to provide personal information to software or hardware such as an electronic device, application, server, or storage medium that performs the operation of the technical solution of the present disclosure according to the prompt message.

[0057] As an optional but non-limiting implementation, in response to receiving an active request from the user, the prompt information may be sent to the user in the form of a pop-up window, in which the prompt information may be presented in text form. In addition, the pop-up window may also carry a selection control for the user to choose "agree" or "disagree" to provide personal information to the electronic device.

[0058] It is understandable that the above notification and the process of obtaining user authorization are merely illustrative and do not constitute a limitation on the implementation of the present disclosure. Other methods that meet the relevant laws and regulations may also be applied to the implementation of the present disclosure.

[0059] At the same time, it is understandable that the data involved in this technical solution (including but not limited to the data itself, the acquisition or use of the data) shall comply with the requirements of relevant laws, regulations and relevant provisions.

[0060] The following will be combined with the drawings in the embodiments of the present application to clearly and completely describe the technical solutions in the embodiments of the present application. Obviously, the described embodiments are only part of the embodiments of the present application, not all of the embodiments. Based on the embodiments in the present application, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of this application.

[0061] At present, the least squares migration imaging technology can effectively eliminate the problems caused by seismic observation, wavelet band limitation, and spatial stretching by utilizing the inverse Hessian, and is an important technical means for high-precision seismic imaging.

[0062] In the face of the problem of fine imaging of fracture-vuggy reservoirs, scattered wave imaging can effectively solve the problem of high-quality imaging of small-scale anomalies by stripping off the reflected energy in the seismic wave field.

[0063] The scattered wave least squares imaging technology proposed earlier has realized scattered wave separation imaging under the least squares inversion imaging framework, and has shown certain advantages in imaging small-scale anomalies. However, the separation imaging strips off the reflection imaging information, which is not conducive to the subsequent reservoir inversion interpretation.

[0064] In view of the fact that the scattered wave least squares imaging in the previous study can only separate the scattered wave information, but cannot provide the full wave field imaging result, which is not conducive to the inversion interpretation of the reservoir, the applicant proposed the scattered wave high-precision imaging method, device, electronic device and storage medium provided in the embodiment of the present application, by obtaining the least squares imaging result, processing the least squares imaging result based on the pre-constructed scattering Fourier filter operator to determine the scattered wave imaging result, scanning the preset scattered wave target area based on the scattered wave imaging result to obtain the energy compensation factor, and weightedly superimposing the scattered wave imaging result and the least squares imaging result based on the energy compensation factor to generate the scattered wave high-precision imaging result, by using the post-stack scattered wave imaging technology to separate the scattered energy under the least squares inversion imaging framework, using the energy compensation factor to realize the scattered wave enhanced imaging, and improving the high-resolution and high-display imaging effect of the fracture-cavity reservoir. Among them, the scattered wave high-precision imaging method is described in detail in the subsequent embodiments.

[0065] Embodiment 1

[0066] The following is an introduction to the application scenarios of the scattered wave high-precision imaging method provided in the embodiments of the present application:

[0067] See also Figure 1 , Figure 1 The figure is a flow chart of a scattered wave high-precision imaging method provided in the embodiment of the present application, so as to improve the high-resolution and high-visibility imaging effect of the fracture-cavity reservoir and provide a technical means for high-precision imaging of the fracture-cavity structure. In the embodiment, the scattered wave high-precision imaging method can be applied to the following: Figure 7 The scattered wave high-precision imaging method and device 300 shown in FIG. Figure 8 In the electronic device 200 shown, the electronic device 200 can be a smart terminal such as a desktop computer, a tablet computer, a smart phone, etc., wherein the electronic device 200 can include one or more, and information can be transmitted between multiple electronic devices through wireless and / or wired methods. Multiple electronic devices can collaborate to complete the scattered wave high-precision imaging method. By way of example, various types of data can be obtained through the smart terminal to complete the scattered wave high-precision imaging method process.

[0068] The following is for Figure 1 The process shown is described in detail, and the scattered wave high-precision imaging method may include S110 to S140.

[0069] S110: Obtaining least squares imaging results.

[0070] In the embodiment of the present application, the least squares imaging result is an image generated by minimizing the difference between the observed data and the predicted data. It is used to reconstruct the spatial distribution of the target from the collected data. Least squares imaging uses an optimization process to find the image that best matches the observed data, thereby improving the accuracy and resolution of imaging.

[0071] In the field of seismic exploration, least squares imaging can be used to reconstruct underground structures to understand the distribution and properties of underground strata. By collecting seismic data and applying least squares imaging technology, seismic profiles or seismic volume maps can be generated to reveal the characteristics of underground strata, such as rock types, oil and gas reservoir locations, and structural changes.

[0072] S120: Processing the least square imaging result based on the pre-constructed scattering Fourier filter operator to determine the scattered wave imaging result.

[0073] In the embodiment of the present application, the scattering Fourier filter operator is used to improve the quality of seismic data and interpret underground structures. The operator performs filtering operations on seismic imaging data in the frequency domain to enhance or suppress the energy of different frequency components.

[0074] It should be noted that the role of the scattering Fourier filter operator is to process the data in the wavenumber domain according to the spectral characteristics of the seismic data, which can change the wavenumber components of the data, thereby highlighting or suppressing specific underground flat reflection structure information. By selecting appropriate filtering parameters and applying the scattering Fourier filter operator, the signal in the seismic data can be enhanced, noise interference can be reduced, and the resolution and clarity of the underground structure can be improved.

[0075] In some aspects, the scattering Fourier filter operator is usually designed based on prior information and the characteristics of seismic data, and it can select appropriate filtering parameters according to the scattering characteristics of the underground medium and the characteristics of the exploration target. By filtering the seismic data in the frequency domain, the scattering Fourier filter operator can improve the spectral distribution of the data, making the energy of the underground target more prominent, thereby improving the accuracy of interpretation and imaging.

[0076] In some aspects, scattered wave imaging is usually based on the propagation direction information of scattered waves. When seismic waves propagate in the underground medium, they scatter with the inhomogeneities in the medium to produce scattered waves. These scattered waves can carry information about the microstructure of the underground medium, such as the pores, cracks, and stratum interfaces of the rock.

[0077] Among them, during the scattered wave imaging process, the propagation characteristics of the scattered waves of the anomaly are analyzed, and the information such as the location, shape and intensity of the underground anomaly is reflected in the imaging results.

[0078] S130: Scanning a preset scattered wave target area based on the scattered wave imaging result to obtain an energy compensation factor.

[0079] In the embodiment of the present application, by analyzing the scattered wave imaging results, the preset scattered wave target area is found, and the corresponding energy compensation factor is calculated. In the scattered wave imaging results, the scattered target is usually manifested as an area where local energy is concentrated. By analyzing the imaging results, the position and shape of the scattered target can be determined and matched with the preset scattered wave target area.

[0080] In some embodiments, an energy compensation factor can be calculated in a preset scattered wave target area by scanning the imaging results. The energy compensation factor represents the relative compensation degree of the energy of the scattered waves in the area relative to other areas. By calculating the energy compensation factor, the energy in the scattered wave target area can be adjusted to improve the energy resolution and clarity of the target area.

[0081] S140: Performing weighted superposition on the scattered wave imaging result and the least square imaging result based on the energy compensation factor to generate a scattered wave high-precision imaging result.

[0082] The scattered wave high-precision imaging result may include a least squares scattering enhancement imaging result.

[0083] In the embodiment of the present application, the energy compensation factor is used to adjust the weight between the scattered wave imaging result and the least squares imaging result. The energy compensation factor is obtained by analyzing the scattered wave imaging result and is used to indicate the relative energy contribution of the scattered wave in different regions. The energy compensation factor is used as a weight factor to perform weighted averaging on the scattered wave imaging result and the least squares imaging result. The weighted averaging can retain more energy in the scattered wave target region, thereby improving the resolution and clarity of the imaging result in the target region.

[0084] In this embodiment, under the framework of least squares migration theory, scattered wave separation and energy compensation factors are used to realize separated imaging and adaptive amplitude compensation of small-scale anomalies, and energy compensation factors are used to realize scattered wave enhanced imaging, thereby improving the high-resolution and high-visibility imaging effect of fracture-vuggy reservoirs, and effectively improving the imaging resolution and visibility of small-scale anomalies.

[0085] In some embodiments, obtaining the least squares imaging result includes: may include S210 to S220.

[0086] S210: Obtain the reverse time migration image and global space-variant scattering function in the seismic data.

[0087] In the embodiment of the present application, the global space-variant scattering function may be a function that describes the scattering characteristics of underground media. The reverse time migration is first used to obtain the reverse time migration image in the seismic data, which can be obtained by reverse propagating the seismic data and stacking them.

[0088] S220: Perform high-dimensional spatial deconvolution processing based on the reverse time migration image and the global space-variant scattering function to determine the least squares imaging result.

[0089] In the embodiment of the present application, the reverse time migration image and the global space-varying scattering function are subjected to high-dimensional spatial deconvolution processing in order to correct or improve the scattering effect in the seismic data to improve the imaging quality and accuracy.

[0090] S110: Obtaining least squares imaging results.

[0091] In some possible implementations, obtaining the least squares imaging result may be implemented by the following steps:

[0092] S112: Input the reverse time migration image p and the global space-variant scattering function k.

[0093] S114: Divide the reverse time migration image p and the global space-variant scattering function k into regional blocks.

[0094] S116: For any block, use Fourier transform to transform the reverse time offset image p i and the global space-variant scattering function k i Switch to the frequency wavenumber domain and perform spatial deconvolution to obtain the least squares imaging result l i .

[0095] Among them, the expression of S116 can be:

[0096]

[0097] In the formula, α is the regularization factor, α is k i k i * 1% of the maximum value, k i * k i The conjugate function of .

[0098] S120: Processing the least square imaging result based on the pre-constructed scattering Fourier filter operator to determine the scattered wave imaging result.

[0099] In some possible implementations, S120 may include the following steps:

[0100] S122: Construct the scattering Fourier filter operator g in the wave number domain to perform the least squares imaging result l iFilter and obtain the wavenumber domain scattered wave imaging result s i .

[0101] S124: The least squares imaging results of scattered waves for all blocks s i Perform stitching to obtain the final scattered wave least squares imaging result s.

[0102] S130: Scanning a preset scattered wave target area based on the scattered wave imaging result to obtain an energy compensation factor.

[0103] S140: Performing weighted superposition on the scattered wave imaging result and the least square imaging result based on the energy compensation factor to generate a scattered wave high-precision imaging result.

[0104] In some embodiments, the energy compensation factor is calculated as:

[0105]

[0106] Where U is the mean amplitude of the scattered wave imaging result, U0 is the mean amplitude of the least squares imaging result, rms is the mean square error, sigmoid is the activation function, θ s is the opening angle, θ r is the azimuth.

[0107] In some embodiments, the scattered wave imaging result and the least squares imaging result are weightedly superimposed based on the energy compensation factor to generate a scattered wave high-precision imaging result, using the calculation formula:

[0108] m enh =m lsm +F·m svm ;

[0109] In the formula, m lsm is the least squares imaging result, m svm is the scattered wave imaging result, m enh It is the result of high-precision imaging of scattered waves.

[0110] In some embodiments, the calculation formula of the pre-constructed scattering Fourier filter operator is:

[0111]

[0112] Where θ is the scattering angle corresponding to the spatial wave number distribution, and its expression is α is the damping factor, k is the spatial wave number, x, y, z are the spatial coordinate directions, θ0 is the scattering angle of the imaging point, and r is the attenuation smoothing function.

[0113] Among them, θ0 is the scattering angle of the imaging point, which is adjusted according to the filtering requirements and ranges from [0,90] degrees. In order to avoid the truncation effect of Fourier transform, the scattering Fourier filter operator g needs to use an attenuation smoothing function r.

[0114] In some embodiments, according to the scattering Fourier filter operator g, the reverse time migration image l in the wave number domain is i Filter and obtain the wavenumber domain scattered wave imaging result s i , which is expressed as follows:

[0115] s i = l i g;

[0116] In some embodiments, the least squares imaging result is processed based on a pre-constructed scattering Fourier filter operator to determine the scattered wave imaging result in S120, including:

[0117] The least square imaging result is filtered based on the scattering Fourier filter operator to obtain the scattering wave imaging result in the wave number domain;

[0118] The scattered wave imaging results in the wavenumber domain are spliced ​​and processed to determine the scattered wave imaging results.

[0119] In summary, the present embodiment provides a scattered wave high-precision imaging method, which obtains a least squares imaging result, processes the least squares imaging result based on a pre-constructed scattering Fourier filter operator to determine the scattered wave imaging result, scans a preset scattering wave target area based on the scattering wave imaging result to obtain an energy compensation factor, performs weighted superposition of the scattering wave imaging result and the least squares imaging result based on the energy compensation factor to generate a scattered wave high-precision imaging result, separates the scattered energy by using post-stack scattered wave imaging technology under the least squares inversion imaging framework, and uses the energy compensation factor to achieve scattered wave enhanced imaging, thereby improving the high-resolution and high-visibility imaging effect of fracture-vuggy reservoirs.

[0120] Embodiment 2

[0121] Based on Example 1, the following uses the actual data processing effect as an example to illustrate the technical features of this application in combination with the compliance and specific implementation methods:

[0122] See also Figure 2 , Figure 2 A schematic diagram of a least squares scattered wave enhanced imaging process provided in an embodiment of the present application, the specific implementation process may include the following steps:

[0123] The reverse time migration results and point spread function can be first input into the model.

[0124] S1: Perform high-dimensional spatial deconvolution on the reverse time migration image and the global space-variant scattering function to obtain the least squares imaging result.

[0125] In this step, the reverse time migration image can be referred to as follows Figure 3 A schematic diagram of a two-dimensional slice of a global space-varying scattering function of three-dimensional actual data is shown in FIG. The processing of the global space-varying scattering function can be referred to as follows: Figure 4 A schematic diagram of a two-dimensional slice of least squares imaging of three-dimensional actual data scattered waves is shown.

[0126] S2: Based on the least squares imaging results, the scattered energy is extracted using the scattering Fourier filter operator to obtain the scattered wave imaging results.

[0127] In this step, the scattered wave imaging results can be found in Figure 5 Another schematic diagram of two-dimensional slice of three-dimensional actual data scattered wave least squares imaging is shown.

[0128] S3: Based on the energy compensation factor, the scattered wave imaging result and the least squares imaging result are weighted superimposed to obtain the least squares scattering enhanced imaging result.

[0129] In this step, the least squares scattering enhanced imaging results can be found in Figure 6 The figure shows a schematic diagram of a two-dimensional slice of least squares scattered wave enhanced imaging of three-dimensional actual data.

[0130] In this step, it can be seen from the above figure that the least squares scattered wave enhanced imaging results not only improve the imaging resolution, enhance the energy of the anomaly, but also comprehensively improve the small-scale crack and hole recognition capability.

[0131] Embodiment 3

[0132] See also Figure 7 , Figure 7 A scattered wave high-precision imaging method and device 300 provided in the present application includes: a first acquisition module 310, a determination module 320, a second acquisition module 330 and a generation module 340, wherein:

[0133] The first acquisition module 310 is used to acquire the least squares imaging result.

[0134] The determination module 320 is used to determine the scattered wave imaging result by processing the least square imaging result based on the pre-constructed scattered wave Fourier filter operator.

[0135] The second acquisition module 330 is used to scan a preset scattered wave target area based on the scattered wave imaging result to acquire an energy compensation factor.

[0136] The generating module 340 is used to perform weighted superposition on the scattered wave imaging result and the least square imaging result based on the energy compensation factor to generate a scattered wave high-precision imaging result.

[0137] Optionally, the first acquisition module 310 may further include: a sub-acquisition module and a sub-confirmation module, wherein:

[0138] The sub-acquisition module is used to obtain the reverse time migration image and the global space-variant scattering function in the seismic data.

[0139] The sub-confirmation module is used to determine the least squares imaging result by performing high-dimensional spatial deconvolution processing based on the reverse time migration image and the global space-variant scattering function.

[0140] It should be noted that the device embodiments in the present application correspond to each other with the aforementioned method embodiments. The specific principles in the device embodiments can be found in the contents of the aforementioned method embodiments and will not be repeated here.

[0141] In several embodiments provided in this embodiment, the coupling between modules may be electrical, mechanical or other forms of coupling.

[0142] In addition, each functional module in each embodiment of the present invention may be integrated into one processing module, or each module may exist physically separately, or two or more modules may be integrated into one module. The above-mentioned integrated module may be implemented in the form of hardware or in the form of software functional modules.

[0143] The flow chart and block diagram in the accompanying drawings illustrate the possible architecture, function and operation of the system, method and computer program product according to various embodiments of the present disclosure. In this regard, each square box in the flow chart or block diagram can represent a module, a program segment or a part of a code, and the module, the program segment or a part of the code contains one or more executable instructions for realizing the specified logical function. It should also be noted that in some implementations as replacements, the functions marked in the square box can also occur in a sequence different from that marked in the accompanying drawings. For example, two square boxes represented in succession can actually be executed substantially in parallel, and they can sometimes be executed in the opposite order, depending on the functions involved. It should also be noted that each square box in the block diagram and / or flow chart, and the combination of the square boxes in the block diagram and / or flow chart can be implemented with a dedicated hardware-based system that performs a specified function or operation, or can be implemented with a combination of dedicated hardware and computer instructions.

[0144] The modules involved in the embodiments described in the present disclosure may be implemented by software or hardware, wherein the name of a module does not, in some cases, limit the module itself.

[0145] Embodiment 4

[0146] See also Figure 8 , Figure 8 The present application provides a structural block diagram of an electronic device 200 that can execute the above-mentioned scattered wave high-precision imaging method. The electronic device 200 can be a computer, a tablet computer, a smart phone, a portable computer or other devices.

[0147] The electronic device 200 further includes a processor 202 and a memory 204. The memory 204 stores a program that can execute the contents of the aforementioned embodiments, and the processor 202 can execute the program stored in the memory 204.

[0148] Among them, the processor 202 may include one or more cores and message matrix units for processing data. The processor 202 uses various interfaces and lines to connect various parts of the entire electronic device 200, and executes various functions and processes data of the electronic device 200 by running or executing instructions, programs, code sets or instruction sets stored in the memory 204, and calling data stored in the memory 204. Optionally, the processor 202 can be implemented in at least one hardware form of digital signal processing (Digital Signal Processing, DSP), field programmable gate array (Field-Programmable Gate Array, FPGA), and programmable logic array (Programmable Logic Array, PLA). The processor 202 can integrate one or a combination of a central processing unit (Central Processing Unit, CPU), a graphics processing unit (Graphics Processing Unit, GPU) and a modulation decoder. Among them, the CPU mainly processes the operating system, user interface and application programs; the GPU is responsible for rendering and drawing display content; and the modem is used to process wireless communications. It can be understood that the above-mentioned modulation decoder may not be integrated into the processor, but may be implemented separately through a communication chip.

[0149] The memory 204 may include a random access memory (RAM) or a read-only memory (ROM). The memory 204 may be used to store instructions, programs, codes, code sets or instruction sets. The memory 204 may include a program storage area and a data storage area, wherein the program storage area may store instructions for implementing an operating system, instructions for implementing at least one function (e.g., instructions for a user to obtain a random number), instructions for implementing the following various method embodiments, etc. The data storage area may also store data (e.g., random numbers) created by the terminal during use, etc.

[0150] The electronic device 200 may also include a network module and a screen. The network module is used to receive and send electromagnetic waves, realize the mutual conversion between electromagnetic waves and electrical signals, and thus communicate with a communication network or other devices, such as communicating with an audio playback device. The network module may include various existing circuit components for performing these functions, such as antennas, radio frequency transceivers, digital signal processors, encryption / decryption chips, user identity modules (SIM) cards, memories, etc. The network module may communicate with various networks such as the Internet, corporate intranets, wireless networks, or communicate with other devices via wireless networks. The above-mentioned wireless networks may include cellular telephone networks, wireless local area networks, or metropolitan area networks. The screen may display interface content and perform data interaction.

[0151] Embodiment 5

[0152] Please refer to Fig. 9 , which shows a block diagram of a computer-readable storage medium provided in an embodiment of the present application. The computer-readable storage medium 400 stores a program code 410, which can be called by a processor to execute the method described in the above method embodiment.

[0153] The computer-readable storage medium 400 may be an electronic memory such as a flash memory, an EEPROM (electrically erasable programmable read-only memory), an EPROM, a hard disk, or a ROM. Optionally, the computer-readable storage medium includes a non-transitory computer-readable storage medium. The computer-readable storage medium 400 has storage space for program codes 410 that perform any method steps in the above method. These program codes 410 can be read from or written to one or more computer program products. The program codes 410 can be compressed, for example, in an appropriate form.

[0154] The embodiment of the present application also provides a computer program product or a computer program, which includes a computer instruction stored in a computer-readable storage medium. A processor of a computer device reads the computer instruction from the computer-readable storage medium, and the processor executes the computer instruction, so that the computer device executes the scattered wave high-precision imaging method described in the above various optional implementations.

[0155] In the several embodiments provided in the present disclosure, it should be understood that the disclosed method can also be implemented in other ways. The method embodiments described above are merely illustrative.

[0156] It should be noted that, in this article, the terms "include", "comprises" or any other variations thereof are intended to cover non-exclusive inclusion, so that a process, method, article or device including a series of elements includes not only those elements, but also includes other elements not explicitly listed, or also includes elements inherent to such process, method, article or device. In the absence of further restrictions, an element defined by the sentence "comprises a ..." does not exclude the presence of other identical elements in the process, method, article or device including the element.

[0157] The above description is only a preferred embodiment of the present disclosure and an explanation of the technical principles used. Those skilled in the art should understand that the scope of disclosure involved in the present disclosure is not limited to the technical solutions formed by a specific combination of the above technical features, but should also cover other technical solutions formed by any combination of the above technical features or their equivalent features without departing from the above disclosed concept. For example, the above features are replaced with the technical features with similar functions disclosed in the present disclosure (but not limited to) by each other to form a technical solution.

[0158] In addition, although each operation is described in a specific order, this should not be understood as requiring these operations to be performed in the specific order shown or in a sequential order. Under certain circumstances, multitasking and parallel processing may be advantageous. Similarly, although some specific implementation details are included in the above discussion, these should not be interpreted as limiting the scope of the present disclosure. Some features described in the context of a separate embodiment can also be implemented in a single embodiment in combination. On the contrary, the various features described in the context of a single embodiment can also be implemented in multiple embodiments individually or in any suitable sub-combination mode.

[0159] The above embodiments are only used to illustrate the technical solutions of the present application, rather than to limit it. Although the present application has been described in detail with reference to the aforementioned embodiments, a person skilled in the art should understand that the technical solutions described in the aforementioned embodiments can still be modified, or some of the technical features can be replaced by equivalents. However, these modifications or replacements do not cause the essence of the corresponding technical solutions to deviate from the spirit and scope of the technical solutions of the embodiments of the present application.

Claims

1. A scattered wave high-precision imaging method, characterized in that: The method comprises: Obtain least squares imaging results; Processing the least squares imaging result based on a pre-constructed scattering Fourier filter operator to determine a scattered wave imaging result; Scanning a preset scattered wave target area based on the scattered wave imaging result to obtain an energy compensation factor; The scattered wave imaging result and the least square imaging result are weightedly superimposed based on the energy compensation factor to generate a scattered wave high-precision imaging result.

2. The method according to claim 1, characterized in that The obtaining of the least squares imaging result comprises: Obtain reverse time migration images and global space-variant scattering functions from seismic data; The least squares imaging result is determined by performing high-dimensional spatial deconvolution processing based on the reverse time migration image and the global space-variant scattering function.

3. The method according to claim 1, characterized in that The calculation formula of the energy compensation factor is: Where U is the mean amplitude of the scattered wave imaging result, U0 is the mean amplitude of the least squares imaging result, rms is the mean square error, sigmoid is the activation function, θ s is the opening angle, θ r is the azimuth.

4. The method according to claim 3, characterized in that The scattered wave imaging result and the least square imaging result are weightedly superimposed based on the energy compensation factor to generate a scattered wave high-precision imaging result, and the calculation formula used is: m enh =m lsm +F·m svm ; In the formula, m lsm is the least squares imaging result, m svm is the scattered wave imaging result, m enh This is the high-precision imaging result of the scattered wave.

5. The method according to claim 2, characterized in that: The calculation formula of the pre-constructed scattering Fourier filter operator is: Where θ is the scattering angle corresponding to the spatial wave number distribution, and its expression is α is the damping factor, k is the spatial wave number, x, y, z are the spatial coordinate directions, θ0 is the scattering angle of the imaging point, and r is the attenuation smoothing function.

6. The method according to claim 5, characterized in that The method of processing the least square imaging result based on the pre-constructed scattering Fourier filter operator to determine the scattered wave imaging result includes: Performing filtering processing on the least squares imaging result based on the scattering Fourier filter operator to obtain a wavenumber domain scattering wave imaging result; The scattered wave imaging results in the wavenumber domain are spliced ​​and processed to determine the scattered wave imaging results.

7. A scattered wave high-precision imaging device, characterized in that: The device comprises: A first acquisition module, used for acquiring least squares imaging results; A determination module, configured to determine a scattered wave imaging result by processing the least square imaging result based on a pre-constructed scattering Fourier filter operator; A second acquisition module, configured to scan a preset scattered wave target area based on the scattered wave imaging result to acquire an energy compensation factor; A generating module is used to perform weighted superposition on the scattered wave imaging result and the least square imaging result based on the energy compensation factor to generate a scattered wave high-precision imaging result.

8. The scattered wave high-precision imaging device according to claim 7, characterized in that: The first acquisition module also includes: Sub-acquisition module, used to obtain reverse time migration image and global space-variant scattering function in seismic data; The sub-confirmation module is used to determine the least squares imaging result by performing high-dimensional spatial deconvolution processing based on the reverse time migration image and the global space-variant scattering function.

9. An electronic device, characterized in that: include: one or more processors; Memory; One or more programs, wherein the one or more programs are stored in the memory and configured to be executed by the one or more processors, and the one or more programs are configured to execute the scattered wave high-precision imaging method according to any one of claims 1 to 6.

10. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores program codes, and the program codes can be called by one or more processors to execute the scattered wave high-precision imaging method according to any one of claims 1 to 6.

Images