VTI medium imaging method, device, equipment, medium and system

Abstract

The method comprises the following steps: acquiring a point scattering model and a reverse time migration image of a VTI medium in a target area, wherein the VTI medium is a transverse isotropic medium of a vertical symmetry axis in the target area; determining a point spread function data volume of the VTI medium according to the point scattering model; transforming the point spread function data body and the reverse time migration image to a frequency wavenumber domain by using Fourier transform to obtain a wavenumber domain reverse time migration image block and a wavenumber domain point spread function data block; and determining a least square image of the VTI medium according to the wave number domain reverse time migration image block and the wave number domain point spread function data block. According to the method, VTI medium least square migration imaging is realized, the influence of anisotropy on construction imaging positioning is corrected, the influence on the imaging position and amplitude is eliminated, and the imaging precision is improved.

Description

Technical Field

[0001] The invention belongs to the technical field of seismic data imaging, and specifically relates to an imaging method, device, equipment, medium and system of a VTI medium. Background Art

[0002] Least squares migration imaging iteratively estimates the underground reflectivity model by matching seismic observation data, and is currently the most accurate inversion imaging theory. The imaging domain least squares migration method calculates the point spread function and approximates the inverse operator of the Hessian. It does not require an iterative step to obtain the least squares migration result, which alleviates the algorithm's requirements for high-precision velocity and wavelets. It has demonstrated good application results in actual seismic data processing. At present, this technology still remains in the problem of high-precision imaging of isotropic media, and insufficient consideration is given to VTI (vertical anisotropy) media. Anisotropy parameters affect the travel time and phase of seismic wave propagation. Therefore, in areas where anisotropy is developed, it is very necessary to develop imaging algorithms for anisotropic media. Summary of the invention

[0003] In response to the above technical problems, the present invention proposes an imaging method for VTI media. This application obtains a point scattering model and a reverse time migration image of the VTI medium in the target area, wherein the VTI medium is a transversely isotropic medium with a vertical symmetry axis in the target area; determines the point spread function data body of the VTI medium according to the point scattering model; uses Fourier transform to transfer the point spread function data body and the reverse time migration image to the frequency wavenumber domain to obtain a reverse time migration image block in the wavenumber domain and a point spread function data block in the wavenumber domain; determines the least squares image of the VTI medium according to the reverse time migration image block in the wavenumber domain and the point spread function data block in the wavenumber domain. Least squares migration imaging of VTI media is realized, the influence of anisotropy on structural imaging positioning is corrected, the influence on imaging position and amplitude is eliminated, and imaging accuracy is improved.

[0004] To solve the above technical problems, the technical solution adopted by the present invention includes four aspects.

[0005] In a first aspect, a method for imaging a VTI medium is provided, comprising acquiring a point scattering model and a reverse time migration image of the VTI medium in a target area, wherein the VTI medium is a transversely isotropic medium with a vertical symmetry axis of a vertically anisotropic medium in the target area; determining a point spread function data volume of the VTI medium according to the point scattering model; transferring the point spread function data volume and the reverse time migration image to a frequency wavenumber domain by Fourier transform to obtain a wavenumber domain reverse time migration image block and a wavenumber domain point spread function data block; and determining a least squares image of the VTI medium according to the wavenumber domain reverse time migration image block and the wavenumber domain point spread function data block.

[0006] In some embodiments, determining the point spread function data volume of the VTI medium according to the point scattering model includes: acquiring the migration velocity corresponding to the reverse time migration image and various anisotropic parameters of the VTI medium; reversely migrating the point scattering model according to the migration velocity and the anisotropic parameters to obtain a disturbed wave field; and reverse time migrating the disturbed wave field into a point spread function data volume.

[0007] In some embodiments, the point spread function data volume and the reverse time migration image are transferred to the frequency wave number domain by Fourier transform to obtain the wave number domain reverse time migration image block and the wave number domain point spread function data block. The point spread function data volume and the reverse time migration image are transferred to the frequency wave number domain by Fourier transform to obtain the wave number domain reverse time migration image block and the wave number domain point spread function data block, including: obtaining the model parameters of the point scattering model; dividing the point spread function data volume and the reverse time migration image into regional blocks according to the model parameters to obtain multiple point spread function data blocks and reverse time migration image blocks; converting each of the point spread function data volume blocks into a wave number domain point spread function data block by Fourier transform; converting each of the reverse time migration image blocks into a wave number domain reverse time migration image block by Fourier transform.

[0008] In some embodiments,

[0009] The method of determining the least squares image of the VTI medium based on the wavenumber domain reverse time migration image block and the wavenumber domain point spread function data block includes: performing high-dimensional spatial deconvolution on each of the wavenumber domain point spread function data blocks and the corresponding wavenumber domain reverse time migration image block to obtain a wavenumber domain least squares image block; performing Fourier transformation on each of the wavenumber domain least squares image blocks to obtain the least squares image block; and merging multiple least squares image blocks to obtain the least squares image of the VTI medium.

[0010] In some embodiments, demigrating the point scattering model according to the migration velocity and the anisotropy parameter to obtain a disturbed wave field includes: acquiring an acoustic seismic wave field of the target area; demigrating the point scattering model according to the acoustic seismic wave field, the migration velocity and the anisotropy parameter to obtain a disturbed wave field.

[0011] In some embodiments, the de-migration of the point scattering model according to the acoustic seismic wave field, the migration velocity and the anisotropy parameter to obtain the disturbance wave field is achieved by the following formula: Where P and Q are acoustic seismic wave fields, ε and δ are anisotropy parameters, d represents the disturbance term, and v p represents the migration velocity, dP represents the disturbance wave field, dv prepresents the disturbance parameter, x represents the x-axis direction, y represents the y-axis direction, z represents the z-axis direction, and t represents the time axis direction.

[0012] In some embodiments, the step of performing high-dimensional spatial deconvolution on the plurality of wavenumber domain point spread function data blocks and the plurality of wavenumber domain reverse time migration image blocks to obtain the wavenumber domain least squares image block is implemented according to the following formula: represents the wavenumber domain reverse time migration image block, represents the wavenumber domain point spread function data block, represents the least squares image block in the wavenumber domain, α is the regularization factor, and is 5% of maximum value.

[0013] In a second aspect, the present application proposes an imaging device for a VTI medium, comprising: a first acquisition module, used to acquire a point scattering model and a reverse time migration image of the VTI medium in a target area, wherein the VTI medium is a transversely isotropic medium or a vertically anisotropic medium with a vertical symmetry axis in the target area; a first determination module, used to determine a point spread function data volume of the VTI medium according to the point scattering model; a first execution module, used to transfer the point spread function data volume and the reverse time migration image to a frequency wave number domain by using Fourier transform, and obtain a wave number domain reverse time migration image block and a wave number domain point spread function data block; a second determination module, used to determine a least squares image of the VTI medium according to the wave number domain reverse time migration image block and the wave number domain point spread function data block.

[0014] In a third aspect, the present application proposes an electronic device, comprising: a memory and a processor, wherein the memory stores a computer program, and when the computer program is executed by the processor, the method described in the first aspect is performed.

[0015] In a fourth aspect, the present application proposes a storage medium, which stores a computer program that can be executed by one or more processors, and the computer program can be used to implement the method described in the first aspect.

[0016] Beneficial effects created by the present invention: The present application obtains the point scattering model and reverse time migration image of the VTI medium in the target area, wherein the VTI medium is a transversely isotropic medium with a vertical symmetry axis in the target area; determines the point spread function data body of the VTI medium according to the point scattering model; uses Fourier transform to transfer the point spread function data body and the reverse time migration image to the frequency wave number domain to obtain the wave number domain reverse time migration image block and the wave number domain point spread function data block; determines the least square image of the VTI medium according to the wave number domain reverse time migration image block and the wave number domain point spread function data block. Least square migration imaging of VTI medium is realized, the influence of anisotropy on structural imaging positioning is corrected, the influence on imaging position and amplitude is eliminated, and imaging accuracy is improved. BRIEF DESCRIPTION OF THE DRAWINGS

[0017] The scope of the present disclosure may be better understood by reading the following detailed description of exemplary embodiments in conjunction with the accompanying drawings. The drawings included are:

[0018] Figure 1 An overall flow chart of an imaging method of a VTI medium provided in an embodiment of the present application;

[0019] Figure 2 This is a structural block diagram of an imaging device of a VTI medium provided in an embodiment of the present application. DETAILED DESCRIPTION

[0020] In order to make the purpose, technical solutions and advantages of the present application clearer, the present application will be further described in detail below in conjunction with the accompanying drawings. The described embodiments should not be regarded as limiting the present application. All other embodiments obtained by ordinary technicians in the field without making creative work are within the scope of protection of this application.

[0021] In the following description, reference is made to “some embodiments”, which describe a subset of all possible embodiments, but it will be understood that “some embodiments” may be the same subset or different subsets of all possible embodiments and may be combined with each other without conflict.

[0022] If similar descriptions of "first\second\third" appear in the application documents, the following instructions will be added. In the following description, the terms "first\second\third" involved are merely used to distinguish similar objects and do not represent a specific ordering of the objects. It can be understood that "first\second\third" can be interchanged in a specific order or sequence where permitted, so that the embodiments of the present application described herein can be implemented in an order other than that illustrated or described herein.

[0023] Unless otherwise defined, all technical and scientific terms used herein have the same meaning as those commonly understood by those skilled in the art to which this application belongs. The terms used herein are only for the purpose of describing the embodiments of this application and are not intended to limit this application.

[0024] Embodiment 1:

[0025] Least squares migration imaging iteratively estimates the underground reflectivity model by matching seismic observation data, and is currently the most accurate inversion imaging theory. The imaging domain least squares migration method calculates the point spread function and approximates the inverse operator of the Hessian. It does not require an iterative step to obtain the least squares migration result, which alleviates the algorithm's requirements for high-precision velocity and wavelets. It has demonstrated good application results in actual seismic data processing. At present, this technology still remains in the problem of high-precision imaging of isotropic media, and insufficient consideration is given to VTI (vertical anisotropy) media. Anisotropy parameters affect the travel time and phase of seismic wave propagation. Therefore, in areas where anisotropy is developed, it is very necessary to develop imaging algorithms for anisotropic media.

[0026] In view of the problems existing in the prior art, such as Figure 1 As shown, the present application provides an imaging method of a VTI medium, and the method is applied to an electronic device, and the electronic device can be a server, a mobile terminal, a computer, a cloud platform, etc. The function implemented by the device data processing provided in the embodiment of the present application can be implemented by calling a program code by a processor of the electronic device, wherein the program code can be stored in a computer storage medium, and the imaging method of the VTI medium includes:

[0027] Step S1: Acquire a point scattering model and a reverse time migration image of a VTI medium in a target area, wherein the VTI medium is a transversely isotropic medium with a vertical symmetry axis in the target area.

[0028] Reverse time migration imaging is obtained by reverse time migration of various anisotropic parameters of VTI media at a preset migration velocity. The point scattering model is obtained by uniformly sparsely sampling the seismic data in the target area.

[0029] Step S2: Determine the point spread function data volume of the VTI medium according to the point scattering model.

[0030] In some embodiments, step S2 of “determining a point spread function data volume of the VTI medium according to the point scattering model” includes:

[0031] Step S21: obtaining the migration velocity corresponding to the reverse time migration image and various anisotropic parameters of the VTI medium.

[0032] Step S22: De-migrate the point scattering model according to the migration velocity and the anisotropy parameter to obtain a disturbance wave field.

[0033] In some embodiments, step S22 of "de-migrating the point scattering model according to the migration velocity and the anisotropy parameter to obtain a disturbance wave field" includes:

[0034] Step S221: Acquire the acoustic seismic wave field of the target area.

[0035] Step S222: De-migrate the point scattering model according to the acoustic seismic wave field, the migration velocity and the anisotropy parameter to obtain a disturbance wave field.

[0036] In some embodiments, step S222 is implemented by the following formula:

[0037]

[0038] Where P and Q are acoustic seismic wave fields, ε and δ are anisotropy parameters, d represents the disturbance term, and v p represents the migration velocity, dP represents the disturbance wave field, dv p represents the disturbance parameter, x represents the x-axis direction, y represents the y-axis direction, z represents the z-axis direction, and t represents the time axis direction.

[0039] Step S23: reverse time migration of the disturbance wave field into a point spread function data volume.

[0040] Since it is impossible to directly obtain the point spread function data volume through the point scattering model, it is necessary to reverse migrate the point spread model through the acoustic wave migration equation to convert the point spread model into waveform data, that is, the disturbed wave field in this application, and then reverse time migrate the disturbed wave field to obtain a point spread function data volume that can be processed.

[0041] Step S3: using Fourier transform to convert the point spread function data volume and the reverse time migration image into the frequency wave number domain, so as to obtain a wave number domain reverse time migration image block and a wave number domain point spread function data block.

[0042] Since the obtained reverse time migration image and point spread function data body are in the time domain, it is not convenient to fit them. Therefore, in this application, they need to be transferred to the wave number domain to obtain the wave number domain reverse time migration image block and wave number domain point spread function data block.

[0043] Therefore, in some embodiments, step S3 "using Fourier transform to convert the point spread function data volume and the reverse time migration image into the frequency wave number domain to obtain a wave number domain reverse time migration image block and a wave number domain point spread function data block" includes:

[0044] Step S31: obtaining model parameters of the point scattering model.

[0045] Step S32: dividing the point spread function data volume and the reverse time migration image into blocks according to the model parameters to obtain a plurality of point spread function data blocks and reverse time migration image blocks.

[0046] Step S33: converting each of the point spread function data blocks into a wave number domain point spread function data block through Fourier transform.

[0047] Step S34: convert each of the reverse time migrated image blocks into a reverse time migrated image block in the wave number domain through Fourier transformation.

[0048] Since the point scattering model is obtained by uniform sparse sampling, the meaningful data in the point spread function data volume obtained by the point scattering model is discontinuous. In order to improve the accuracy of the least squares image, it is necessary to perform regional block division on the point spread function data volume. When performing regional block division, in order to make the data volume in each regional block meaningful, it is necessary to obtain the model parameters of the point scattering module, and then perform regional block division according to the model parameters. Correspondingly, since the least squares image needs to be determined by the reverse time migration image and the point spread function data volume, the reverse time migration image also needs to be correspondingly regionally blocked, that is, there is a regional correspondence between the obtained point spread function data block and the reverse time migration image block.

[0049] Then, each point spread function data block and reverse time migration image block is subjected to Fourier transformation to obtain the corresponding wave number domain point spread function data block and wave number domain reverse time migration block, wherein since the wave number domain point spread function data block is obtained by converting a point spread function data block, it is actually a data block. Similarly, the wave number domain reverse time migration image block is also a data block, so multiple wave number domain point spread data blocks and wave number domain reverse time migration image blocks will be obtained in the end.

[0050] Step S4: determining the least square image of the VTI medium according to the wavenumber domain reverse time migration image block and the wavenumber domain point spread function data block.

[0051] In some embodiments, step S4 "determining the least squares image of the VTI medium according to the wavenumber domain reverse time migration image block and the wavenumber domain point spread function data block" includes:

[0052] Step S41: performing high-dimensional spatial deconvolution on each of the wavenumber domain point spread function data blocks and the corresponding wavenumber domain reverse time migration image block to obtain a wavenumber domain least squares image block.

[0053] In some embodiments, step S41 is implemented by the following formula:

[0054]

[0055] represents the reverse time migration in the wavenumber domain, represents the wavenumber domain point spread function data block, represents the least squares image block in the wavenumber domain, α is the regularization factor, and is 5% of maximum value.

[0056] Step S42: Perform Fourier transformation on each of the wavenumber domain least squares image blocks to obtain the least squares image block.

[0057] Step S43: merging the plurality of least square image blocks to obtain the least square image of the VTI medium.

[0058] Since the point spread data volume and the reverse time migration image have been regionally divided in the previous step, a plurality of wave number domain point spread data blocks and wave number domain reverse time migration image blocks are obtained. Therefore, when determining the least squares image, it is necessary to first perform high-dimensional spatial deconvolution on each wave number domain point spread data block and the reverse time migration image block belonging to the same region as the wave number domain point spread data block, wherein the same region refers to being in the same region in the regional division. Since the wave number domain data is used when performing high-dimensional spatial deconvolution, and they are all data blocks, the wave number domain least squares image blocks are also obtained, so there are also multiple wave number domain least squares image blocks existing at this time.

[0059] Then, each wave number domain least squares image block is Fourier transformed to obtain multiple time domain least squares images. Then, multiple time domain least squares images are fused to obtain the least squares image of the VTI medium.

[0060] This application obtains the point scattering model and reverse time migration image of the VTI medium in the target area, wherein the VTI medium is a transversely isotropic medium with a vertical symmetry axis in the target area; determines the point spread function data body of the VTI medium according to the point scattering model; uses Fourier transform to transfer the point spread function data body and the reverse time migration image to the frequency wave number domain to obtain the wave number domain reverse time migration image block and the wave number domain point spread function data block; determines the least square image of the VTI medium according to the wave number domain reverse time migration image block and the wave number domain point spread function data block. Least square migration imaging of VTI medium is realized, the influence of anisotropy on structural imaging positioning is corrected, the influence on imaging position and amplitude is eliminated, and imaging accuracy is improved.

[0061] Embodiment 2:

[0062] Based on the foregoing embodiments, the embodiments of the present application provide an imaging device of a VTI medium. The modules included in the device and the units included in the modules can be implemented by a processor in a computer device; of course, they can also be implemented by a specific logic circuit; in the implementation process, the processor can be a central processing unit (CPU, Central Processing Unit), a microprocessor (MPU, Microprocessor Unit), a digital signal processor (DSP, Digital Signal Processing) or a field programmable gate array (FPGA, Field Programmable Gate Array), etc.

[0063] like Figure 2 As shown, a VTI medium imaging device includes: a first acquisition module 1, a first determination module 2, a first execution module 3 and a second execution module 4.

[0064] The first acquisition module 1 is used to acquire the point scattering model and reverse time migration image of the VTI medium in the target area, wherein the VTI medium is a transversely isotropic medium with a vertical symmetry axis in the target area. The first determination module 2 is used to determine the point spread function data volume of the VTI medium according to the point scattering model. The first execution module 3 is used to convert the point spread function data volume and the reverse time migration image into the frequency wave number domain by Fourier transform to obtain the wave number domain reverse time migration image block and the wave number domain point spread function data block. The second determination module is used to determine the least squares image of the VTI medium according to the wave number domain reverse time migration image block and the wave number domain point spread function data block.

[0065] In some embodiments, the first determination module 2 includes: a second acquisition module, a third execution module and a fourth execution module.

[0066] The second acquisition module is used to obtain the migration velocity corresponding to the reverse time migration image and the anisotropic parameters of the VTI medium. The third execution module is used to reverse migrate the point scattering model according to the migration velocity and the anisotropic parameters to obtain a disturbance wave field. The fourth execution module is used to reverse time migrate the disturbance wave field into a point spread function data volume.

[0067] In some embodiments, the first execution module 3 includes: a third acquisition module, a fifth execution module, a sixth execution module and a seventh execution module.

[0068] The third acquisition module is used to obtain the model parameters of the point scattering model. The fifth execution module is used to regionally block the point spread function data volume and the reverse time migration image according to the model parameters to obtain multiple point spread function data blocks and reverse time migration image blocks. The sixth execution module is used to convert each of the point spread function data blocks into a wave number domain point spread function data block through Fourier transformation. The seventh execution module is used to convert each of the reverse time migration image blocks into a wave number domain reverse time migration image block through Fourier transformation.

[0069] In some embodiments, the second execution module 4 includes: an eighth execution module, a ninth execution module and a tenth execution module.

[0070] The eighth execution module is used to perform high-dimensional spatial deconvolution on each of the wavenumber domain point spread function data blocks and the corresponding wavenumber domain reverse time migration image blocks to obtain the wavenumber domain least squares image blocks. The ninth execution module is used to perform Fourier transformation on each of the wavenumber domain least squares image blocks to obtain the least squares image blocks. The tenth execution module is used to merge multiple least squares image blocks to obtain the least squares image of the VTI medium.

[0071] In some embodiments, the third execution module includes: a fourth acquisition module and an eleventh execution module.

[0072] The fifth acquisition module is used to acquire the acoustic seismic wave field of the target area. The eleventh execution module is used to perform de-migration on the point scattering model according to the acoustic seismic wave field, the migration velocity and the anisotropy parameter to obtain a disturbance wave field.

[0073] Each module in the imaging device of the VTI medium described above can be implemented in whole or in part by software, hardware, or a combination thereof. Each module described above can be embedded in or independent of the processor in the device in the form of hardware, or can be stored in the memory in the processing device in the form of software, so that the processor can call and execute the operations corresponding to each of the above modules. It should be noted that the division of modules in the embodiment of the present application is schematic, and is only a logical function division. There may be other division methods in actual implementation.

[0074] Embodiment 3:

[0075] A third aspect provides an electronic device, including a storage device and a processor, wherein the storage device stores a computer program, and the processor implements the steps of a VTI medium imaging method when executing the computer program.

[0076] Embodiment 4:

[0077] A fourth aspect provides a storage medium, wherein a computer program stored in the storage medium can be executed by one or more processors, and the computer program can be used to implement the steps of any one of the VTI medium imaging methods in the first aspect.

[0078] Those of ordinary skill in the art can understand that all or part of the processes in the above-mentioned embodiment methods can be completed by instructing the relevant hardware through a computer program, and the computer program can be stored in a non-volatile computer-readable storage medium. When the computer program is executed, it can include the processes of the embodiments of the above-mentioned methods. Among them, any reference to memory, storage, database or other media used in the embodiments provided in this application can include at least one of non-volatile and volatile memory. Non-volatile memory can include read-only memory (ROM), magnetic tape, floppy disk, flash memory or optical memory, etc. Volatile memory can include random access memory (RAM) or external cache memory. As an illustration and not limitation, RAM can be in various forms, such as static random access memory (SRAM) or dynamic random access memory (DRAM).

[0079] It should be understood that "one embodiment" or "an embodiment" mentioned throughout the specification means that specific features, structures or characteristics related to the embodiment are included in at least one embodiment of the present application. Therefore, "in one embodiment" or "in an embodiment" appearing throughout the specification does not necessarily refer to the same embodiment. In addition, these specific features, structures or characteristics can be combined in one or more embodiments in any suitable manner. It should be understood that in various embodiments of the present application, the size of the sequence number of the above-mentioned processes does not mean the order of execution, and the execution order of each process should be determined by its function and internal logic, and should not constitute any limitation on the implementation process of the embodiments of the present application. The above-mentioned sequence numbers of the embodiments of the present application are only for description and do not represent the advantages and disadvantages of the embodiments.

[0080] 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 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 existence of other identical elements in the process, method, article or device including the element.

[0081] In the several embodiments provided in the present application, it should be understood that the disclosed devices and methods can be implemented in other ways. The device embodiments described above are only schematic. For example, the division of the units is only a logical function division. There may be other division methods in actual implementation, such as: multiple units or components can be combined, or can be integrated into another system, or some features can be ignored or not executed. In addition, the coupling, direct coupling, or communication connection between the components shown or discussed can be through some interfaces, and the indirect coupling or communication connection of the devices or units can be electrical, mechanical or other forms.

[0082] The units described above as separate components may or may not be physically separated, and the components displayed as units may or may not be physical units; they may be located in one place or distributed on multiple network units; some or all of the units may be selected according to actual needs to achieve the purpose of the present embodiment.

[0083] In addition, all functional units in the embodiments of the present application may be integrated into one processing unit, or each unit may be a separate unit, or two or more units may be integrated into one unit; the above-mentioned integrated units may be implemented in the form of hardware or in the form of hardware plus software functional units.

[0084] A person of ordinary skill in the art can understand that: all or part of the steps of implementing the above-mentioned method embodiment can be completed by hardware related to program instructions, and the aforementioned program can be stored in a computer-readable storage medium. When the program is executed, it executes the steps of the above-mentioned method embodiment; and the aforementioned storage medium includes: mobile storage devices, read-only memories (ROM, Read Only Memory), disks or optical disks, etc. Various media that can store program codes.

[0085] Alternatively, if the above-mentioned integrated unit of the present application is implemented in the form of a software function module and sold or used as an independent product, it can also be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the embodiment of the present application can essentially or in other words, the part that contributes to the prior art can be embodied in the form of a software product, which is stored in a storage medium and includes a number of instructions for a controller to execute all or part of the methods described in each embodiment of the present application. The aforementioned storage medium includes: various media that can store program codes, such as mobile storage devices, ROMs, magnetic disks, or optical disks.

[0086] The above is only an implementation method of the present application, but the protection scope of the present application is not limited thereto. Any person skilled in the art who is familiar with the present technical field can easily think of changes or substitutions within the technical scope disclosed in the present application, which should be included in the protection scope of the present application. Therefore, the protection scope of the present application should be based on the protection scope of the claims.

Claims

1. A method for imaging a VTI medium, characterized in that: include: Acquire a point scattering model and a reverse time migration image of a VTI medium in a target area, wherein the VTI medium is a transversely isotropic medium with a vertical symmetry axis in the target area; determining a point spread function data volume of the VTI medium according to the point scattering model; Using Fourier transform to transfer the point spread function data volume and the reverse time migration image to the frequency wave number domain, so as to obtain a wave number domain reverse time migration image block and a wave number domain point spread function data block; A least square image of the VTI medium is determined based on the wavenumber domain reverse time migration image block and the wavenumber domain point spread function data block.

2. The method according to claim 1, characterized in that: Determining the point spread function data volume of the VTI medium according to the point scattering model includes: Acquiring a migration velocity corresponding to the reverse time migration image and various anisotropic parameters of the VTI medium; De-migrating the point scattering model according to the migration velocity and the anisotropy parameter to obtain a disturbance wave field; The disturbance wave field is reverse time migrated into a point spread function data volume.

3. The method according to claim 1, characterized in that The method of transferring the point spread function data volume and the reverse time migration image to the frequency wave number domain by Fourier transform to obtain the wave number domain reverse time migration image block and the wave number domain point spread function data block comprises: Acquiring model parameters of the point scattering model; Dividing the point spread function data volume and the reverse time migration image into blocks according to the model parameters to obtain a plurality of point spread function data blocks and reverse time migration image blocks; Convert each of the point spread function data blocks into a wave number domain point spread function data block by Fourier transform; Each of the reverse time migrated image blocks is converted into a wave number domain reverse time migrated image block by Fourier transformation.

4. The method according to claim 3, characterized in that The determining of the least square image of the VTI medium according to the wavenumber domain reverse time migration image block and the wavenumber domain point spread function data block comprises: Performing high-dimensional spatial deconvolution on each wavenumber domain point spread function data block and the corresponding wavenumber domain reverse time migration image block to obtain a wavenumber domain least squares image block; Performing Fourier transformation on each of the wavenumber domain least squares image blocks to obtain the least squares image block; The plurality of least square image blocks are merged to obtain the least square image of the VTI medium.

5. The method according to claim 2, characterized in that: De-migrating the point scattering model according to the migration velocity and the anisotropy parameter to obtain a disturbance wave field includes: Acquiring the acoustic seismic wave field of the target area; The point scattering model is de-migrated according to the acoustic seismic wave field, the migration velocity and the anisotropy parameter to obtain a disturbance wave field.

6. The method according to claim 5, characterized in that De-migrating the point scattering model according to the acoustic seismic wave field, the migration velocity and the anisotropy parameter to obtain the disturbance wave field is achieved by the following formula: Where P and Q are acoustic seismic wave fields, ε and δ are anisotropy parameters, d represents the disturbance term, and v p represents the migration velocity, dP represents the disturbance wave field, dv p represents the disturbance parameter, x represents the x-axis direction, y represents the y-axis direction, z represents the z-axis direction, and t represents the time axis direction.

7. The method according to claim 4, characterized in that The high-dimensional spatial deconvolution of the plurality of wavenumber domain point spread function data blocks and the plurality of wavenumber domain reverse time migration image blocks is performed according to the following formula to obtain the wavenumber domain least squares image block: represents the wavenumber domain reverse time migration image block, represents the wavenumber domain point spread function data block, represents the least squares image block in the wavenumber domain, α is the regularization factor, and is 5% of maximum value.

8. An imaging device for VTI medium, characterized in that: include: A first acquisition module is used to acquire a point scattering model and a reverse time migration image of a VTI medium in a target area, wherein the VTI medium is a transversely isotropic medium with a vertical symmetry axis in the target area; A first determination module, configured to determine a point spread function data volume of the VTI medium according to the point scattering model; A first execution module is used to convert the point spread function data volume and the reverse time migration image into a frequency wave number domain by Fourier transform to obtain a wave number domain reverse time migration image block and a wave number domain point spread function data block; The second determination module is used to determine the least square image of the VTI medium according to the wavenumber domain reverse time migration image block and the wavenumber domain point spread function data block.

9. An electronic device, characterized in that: include: A memory and a processor, wherein a computer program is stored in the memory, and when the computer program is executed by the processor, the method according to any one of claims 1 to 7 is executed.

10. A storage medium, characterized in that: The computer program stored in the storage medium can be executed by one or more processors, and the computer program can be used to implement the method according to any one of claims 1 to 7.

Images