Multiple migration imaging method and device based on iterative deconvolution imaging condition

Abstract

The invention provides a multiple migration imaging method and device based on an iterative deconvolution imaging condition, and relates to the technical field of seismic data imaging and migration of geophysical exploration for petroleum, and the method comprises the steps: obtaining a seismic wave forward propagation wave field and a seismic wave backward propagation wave field; performing multiple wave imaging on the seismic wave forward transmission wave field and the seismic wave backward transmission wave field by adopting a deconvolution imaging condition to obtain an initial imaging result; taking the seismic wave forward propagation wave field, the seismic wave backward propagation wave field and the initial imaging result as input, and performing multiple imaging by using an iterative deconvolution imaging method to obtain an iterative deconvolution multiple imaging result; in the multiple wave imaging process, a deconvolution imaging objective function value is determined; and if the objective function value does not meet the preset condition, iterating the imaging result until the deconvolution imaging objective function value meets the preset condition, and generating a final iterative deconvolution multiple imaging result according to the iterated imaging result.

Description

Technical Field

[0001] The invention relates to the field of petroleum geophysical exploration seismic data imaging and migration technology, and in particular to a multiple wave migration imaging method and device based on iterative deconvolution imaging conditions. Background Art

[0002] This section is intended to provide a background or context to the embodiments of the invention recited in the claims. No admission is made that the description herein is prior art by inclusion in this section.

[0003] In traditional marine exploration seismic processing, free surface multiple wave signals are generally removed as interference information, but this part of the information containing underground medium parameters has not been fully utilized and mined, which is a problem and bottleneck in current seismic processing and imaging. Therefore, based on free surface multiple wave imaging, free surface multiple wave information can be fully utilized to perform seismic migration of shallow and middle layers of underground media, realize high-precision imaging of the seabed and its vicinity, and solve the problem of insufficient shallow imaging range of traditional imaging methods;

[0004] Berkhout and Verschuur (1994) proposed free surface multiple wave migration based on the up and down wave framework. Subsequently, Schuster and Rickett (2000) converted multiple waves into primary waves through the interference principle and migrated according to the traditional imaging framework. Shan (2003) perfected this quasi-primary wave imaging theory and applied it to actual production. Lu et al. (2021) systematically studied and improved this type of method and applied it to industrial processing procedures. Multiple wave imaging can be realized in the framework of one-way wave equation and wave equation. Due to the existence of multi-order multiple wave information, the traditional multiple wave inversion method causes multiple wave signals of different orders to image each other, form interference, and reduce the imaging quality. This is also an important problem leading to its industrial application.

[0005] At present, the methods for suppressing crosstalk in the process of multiple wave imaging can be divided into the following categories: the first is to suppress crosstalk by using deconvolution imaging conditions. This method is relatively simple to implement, but this method can only be used as an approximate elimination method to eliminate the influence of crosstalk, and can only suppress or eliminate most of the crosstalk. The second is to image multiple waves of different orders during imaging to avoid imaging crosstalk caused by different orders, but this method cannot avoid the crosstalk problem caused by multiple wave signals of the same order and requires additional separation of multiple wave signals of different orders in seismic data. The third is multiple wave imaging based on the migration algorithm. The core of this type of method is to use the migration algorithm to estimate the phase axis of the crosstalk and then eliminate it. This type of method requires adaptive subtraction processing and requires the separation of primary waves and multiple waves. The fourth method is the least squares type of multiple wave imaging method, which can perfectly solve the imaging crosstalk problem and can be implemented using the one-way wave and wave equation method. However, the calculation amount of this method will increase by several orders of magnitude and have a greater impact on the noise in the data.

[0006] At present, the industrial application of multiple wave imaging mainly adopts the multiple wave migration method under the deconvolution imaging condition, and the imaging data result has low clarity and poor imaging effect. Therefore, a technical solution that can overcome the above defects and improve the imaging effect is urgently needed. Summary of the invention

[0007] In order to solve the problems existing in the prior art, the present invention proposes a multiple wave migration imaging method and device based on iterative deconvolution imaging conditions. The patent of the present invention is different from the multiple wave imaging under the traditional deconvolution imaging conditions. It introduces regularization terms in the frequency and space methods, and realizes a wider frequency of the multiple wave imaging results and a clearer imaging result.

[0008] In a first aspect of an embodiment of the present invention, a multiple wave migration imaging method based on iterative deconvolution imaging conditions is proposed, comprising:

[0009] Obtaining seismic wave forward wave field and seismic wave reverse wave field;

[0010] Performing multiple wave imaging on the forward wave field of the seismic wave and the reverse wave field of the seismic wave using deconvolution imaging conditions to obtain an initial imaging result;

[0011] Taking the seismic wave forward wave field, the seismic wave reverse wave field and the initial imaging result as input, performing multiple wave imaging by using an iterative deconvolution imaging method, and obtaining an iterative deconvolution multiple wave imaging result;

[0012] In the process of multiple wave imaging, the objective function value of deconvolution imaging is determined;

[0013] If the objective function value does not meet the preset condition, the imaging result is iterated until the deconvolution imaging objective function value meets the preset condition, and the final iterative deconvolution multiple wave imaging result is generated according to the iterative imaging result.

[0014] In the multiple wave migration imaging method based on iterative deconvolution imaging conditions in the embodiment of the present invention, obtaining the seismic wave forward propagation wave field includes:

[0015] Get input data; in forward propagation, the input data is seismic full-wave field shot record data containing direct wave, primary reflection wave and free surface multiple wave information.

[0016] According to the input data, a one-way wave equation algorithm is used to perform forward propagation of the wave field to obtain the seismic wave forward propagation wave field.

[0017] In the multiple wave migration imaging method based on iterative deconvolution imaging conditions in the embodiment of the present invention, obtaining the seismic wave back propagation wave field includes:

[0018] Get input data; in reverse propagation, the input data is multiple wave information or seismic full wave field shot records;

[0019] According to the input data, a one-way wave equation algorithm is used to perform reverse propagation of the wave field to obtain the seismic wave reverse propagation wave field.

[0020] In the multiple wave migration imaging method based on iterative deconvolution imaging conditions in the embodiment of the present invention, during the multiple wave imaging process, determining the deconvolution imaging objective function value includes:

[0021] The relationship of the deconvolution imaging objective function is:

[0022]

[0023] Where E(R) is the objective function value;

[0024] U is the back propagation wave field of seismic waves;

[0025] D is the positive wave field of seismic waves;

[0026] R is the imaging result;

[0027] Indicates obtaining the spatial derivative of the imaging result;

[0028] a and β are dimensionless parameters;

[0029] ω is the circular frequency;

[0030] It means to obtain the frequency derivative of the imaging result.

[0031] In the multiple wave migration imaging method based on iterative deconvolution imaging conditions in an embodiment of the present invention, if the objective function value does not meet the preset condition, the imaging result is iterated until the deconvolution imaging objective function value meets the preset condition, and the final iterative deconvolution multiple wave imaging result is generated according to the imaging result after the iteration, including:

[0032] The formula for iterative deconvolution imaging conditions is:

[0033]

[0034] Among them, R n+1 is the imaging result after the n+1th iteration;

[0035] R n is the imaging result after the nth iteration;

[0036] i, j, k, and l are the spatial and frequency coordinates of the imaging results (x, y, z, ω);

[0037] It is the complex conjugate of the forward propagation wave field of the seismic wave.

[0038] In the multiple wave migration imaging method based on iterative deconvolution imaging conditions in an embodiment of the present invention, if the objective function value does not meet the preset condition, the imaging result is iterated until the deconvolution imaging objective function value meets the preset condition, and the final iterative deconvolution multiple wave imaging result is generated according to the imaging result after the iteration, including:

[0039] According to the imaging result R after iteration n+1 , recalculate the objective function value; wherein the objective function value gradually approaches 0 during calculation;

[0040] When the deconvolution imaging objective function value meets the preset condition, a final iterative deconvolution multiple wave imaging result is generated according to the imaging result after iteration.

[0041] In a second aspect of an embodiment of the present invention, a multiple wave migration imaging device based on iterative deconvolution imaging conditions is proposed, comprising:

[0042] A data acquisition module, used to acquire the forward wave field of seismic waves and the reverse wave field of seismic waves;

[0043] A multiple wave imaging processing module, used for performing multiple wave imaging on the seismic wave forward propagation wave field and the seismic wave reverse propagation wave field by using deconvolution imaging conditions to obtain an initial imaging result;

[0044] Taking the seismic wave forward wave field, the seismic wave reverse wave field and the initial imaging result as input, performing multiple wave imaging by using an iterative deconvolution imaging method, and obtaining an iterative deconvolution multiple wave imaging result;

[0045] In the process of multiple wave imaging, the objective function value of deconvolution imaging is determined;

[0046] If the objective function value does not meet the preset condition, the imaging result is iterated until the deconvolution imaging objective function value meets the preset condition, and the final iterative deconvolution multiple wave imaging result is generated according to the iterative imaging result.

[0047] In the multiple wave migration imaging device based on iterative deconvolution imaging conditions in the embodiment of the present invention, the data acquisition module is specifically used for:

[0048] Get input data; in forward propagation, the input data is seismic full-wave field shot record data containing direct wave, primary reflection wave and free surface multiple wave information.

[0049] According to the input data, a one-way wave equation algorithm is used to perform forward propagation of the wave field to obtain the seismic wave forward propagation wave field.

[0050] In the multiple wave migration imaging device based on iterative deconvolution imaging conditions in the embodiment of the present invention, the data acquisition module is specifically used for:

[0051] Get input data; in reverse propagation, the input data is multiple wave information or seismic full wave field shot records;

[0052] According to the input data, a one-way wave equation algorithm is used to perform reverse propagation of the wave field to obtain the seismic wave reverse propagation wave field.

[0053] In the multiple wave migration imaging device based on iterative deconvolution imaging conditions in the embodiment of the present invention, the multiple wave imaging processing module determines the deconvolution imaging objective function value during the multiple wave imaging process, including:

[0054] The relationship of the deconvolution imaging objective function is:

[0055]

[0056] Where E(R) is the objective function value;

[0057] U is the back propagation wave field of seismic waves;

[0058] D is the positive wave field of seismic waves;

[0059] R is the imaging result;

[0060] Indicates obtaining the spatial derivative of the imaging result;

[0061] a and β are dimensionless parameters;

[0062] ω is the circular frequency;

[0063] It means to obtain the frequency derivative of the imaging result.

[0064] In the multiple wave migration imaging device based on iterative deconvolution imaging conditions of the embodiment of the present invention, the multiple wave imaging processing module, if the objective function value does not meet the preset condition, iterates the imaging result until the deconvolution imaging objective function value meets the preset condition, and generates a final iterative deconvolution multiple wave imaging result according to the iterative imaging result, including:

[0065] The formula for iterative deconvolution imaging conditions is:

[0066]

[0067] Among them, R n+1 is the imaging result after the n+1th iteration;

[0068] R n is the imaging result after the nth iteration;

[0069] i, j, k, and l are the spatial and frequency coordinates of the imaging results (x, y, z, ω);

[0070] It is the complex conjugate of the forward propagation wave field of the seismic wave.

[0071] In the multiple wave migration imaging device based on iterative deconvolution imaging conditions in the embodiment of the present invention, the multiple wave imaging processing module is specifically used for:

[0072] According to the imaging result R after iteration n+1 , recalculate the objective function value; wherein the objective function value gradually approaches 0 during calculation;

[0073] When the deconvolution imaging objective function value meets the preset condition, a final iterative deconvolution multiple wave imaging result is generated according to the imaging result after iteration.

[0074] In a third aspect of an embodiment of the present invention, a computer device is proposed, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein when the processor executes the computer program, a multiple wave migration imaging method based on iterative deconvolution imaging conditions is implemented.

[0075] In a fourth aspect of an embodiment of the present invention, a computer-readable storage medium is provided, wherein the computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, a multiple-wave migration imaging method based on iterative deconvolution imaging conditions is implemented.

[0076] In a fifth aspect of the embodiments of the present invention, a computer program product is proposed. The computer program product includes a computer program. When the computer program is executed by a processor, a multiple wave migration imaging method based on iterative deconvolution imaging conditions is implemented.

[0077] In the implementation process of the multiple wave migration imaging method and device based on iterative deconvolution imaging conditions proposed by the present invention, the iterative deconvolution imaging conditions have constraints in the imaging space and frequency directions, and the imaging controllable flexibility is higher. The iterative deconvolution imaging conditions are an extension of the conventional deconvolution imaging method, and the physical basis is more complete. The iterative deconvolution imaging conditions can make the frequency expansion of the imaging results wider and achieve better imaging resolution, providing strong technical support for oil exploration. BRIEF DESCRIPTION OF THE DRAWINGS

[0078] In order to more clearly illustrate the technical solutions of 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 some embodiments of the present application. For ordinary technicians in this field, other drawings can be obtained based on these drawings without creative work.

[0079] Figure 1 It is a schematic flow chart of a multiple wave migration imaging method based on iterative deconvolution imaging conditions according to an embodiment of the present invention.

[0080] Figure 2 It is a schematic diagram of the multiple wave imaging process under iterative deconvolution conditions according to a specific embodiment of the present invention.

[0081] Figure 3A It is a schematic diagram of a full-wavefield seismic record of a single-layer model according to a specific embodiment of the present invention.

[0082] Figure 3B It is a schematic diagram of multiple wave recording according to a specific embodiment of the present invention.

[0083] Figure 4A It is a schematic diagram of the multiple wave migration result of the cross-correlation imaging condition according to a specific embodiment of the present invention.

[0084] Figure 4B yes Figure 4A A partial enlarged schematic diagram of .

[0085] Figure 4CIt is a schematic diagram of the multiple wave migration result of conventional deconvolution imaging conditions according to a specific embodiment of the present invention.

[0086] Figure 4D yes Figure 4C A partial enlarged schematic diagram of .

[0087] Figure 4E It is a schematic diagram of the iterative deconvolution imaging conditional multiple wave migration result according to a specific embodiment of the present invention.

[0088] Figure 4F yes Figure 4E A partial enlarged schematic diagram of .

[0089] Figure 5 It is a schematic diagram of spectrum comparison results under cross-correlation, conventional deconvolution and iterative deconvolution conditions according to a specific embodiment of the present invention.

[0090] Figure 6 It is a schematic diagram of VSP acquisition and observation using the Marmousi model according to a specific embodiment of the present invention.

[0091] Fig. 7A It is a schematic diagram of seismic recording at a receiving point depth of 200m according to a specific embodiment of the present invention.

[0092] Figure 7B It is a schematic diagram of seismic recording at a receiving point depth of 1700m according to a specific embodiment of the present invention.

[0093] Figure 7C It is a schematic diagram of seismic records at a receiving point depth of 3200m according to a specific embodiment of the present invention.

[0094] Fig. 8A It is a schematic diagram of conventional deconvolution migration results according to a specific embodiment of the present invention.

[0095] Figure 8B It is a schematic diagram of iterative deconvolution migration results according to a specific embodiment of the present invention.

[0096] Fig. 9 It is a schematic diagram of the architecture of a multiple wave migration imaging device based on iterative deconvolution imaging conditions according to an embodiment of the present invention.

[0097] Fig.10 It is a schematic diagram of the structure of a computer device according to an embodiment of the present invention. DETAILED DESCRIPTION

[0098] The principles and spirit of the present invention will be described below with reference to several exemplary embodiments. It should be understood that these embodiments are provided only to enable those skilled in the art to better understand and implement the present invention, and are not intended to limit the scope of the present invention in any way. On the contrary, these embodiments are provided to make the present disclosure more thorough and complete, and to fully convey the scope of the present disclosure to those skilled in the art.

[0099] Those skilled in the art will appreciate that the embodiments of the present invention may be implemented as a system, device, apparatus, method or computer program product. Therefore, the present disclosure may be specifically implemented in the following forms, namely: complete hardware, complete software (including firmware, resident software, microcode, etc.), or a combination of hardware and software.

[0100] According to an embodiment of the present invention, a multiple wave migration imaging method and device based on iterative deconvolution imaging conditions are proposed, which relates to the technical field of petroleum geophysical exploration seismic data imaging and migration. The present invention proposes an objective function of regularization constraints of space and frequency domain constraints, and provides an iterative solution method. When performing multiple wave migration, the frequency band range of the migration imaging can be increased as much as possible, and the imaging quality can be improved.

[0101] The principle and spirit of the present invention are explained in detail below with reference to several representative embodiments of the present invention.

[0102] Figure 1 FIG. 1 is a flow chart of a multiple wave migration imaging method based on iterative deconvolution imaging conditions according to an embodiment of the present invention. Figure 1 As shown, the method includes:

[0103] Step 101, obtaining a seismic wave forward propagation wave field and a seismic wave reverse propagation wave field;

[0104] Step 102, performing multiple wave imaging on the forward propagation wavefield of the seismic wave and the reverse propagation wavefield of the seismic wave using deconvolution imaging conditions to obtain an initial imaging result;

[0105] Step 103, taking the seismic wave forward wave field, the seismic wave reverse wave field and the initial imaging result as input, performing multiple wave imaging by using an iterative deconvolution imaging method, and obtaining an iterative deconvolution multiple wave imaging result;

[0106] Step 104, determining a deconvolution imaging objective function value during the multiple wave imaging process;

[0107] Step 105: if the objective function value does not meet the preset condition, the imaging result is iterated until the deconvolution imaging objective function value meets the preset condition, and a final iterative deconvolution multiple wave imaging result is generated according to the iterated imaging result.

[0108] In actual application scenarios, the multiple wave migration imaging method based on iterative deconvolution imaging conditions has the following working principles: obtaining a seismic wave forward wave field and a seismic wave reverse wave field; performing multiple wave imaging on the seismic wave forward wave field and the seismic wave reverse wave field using deconvolution imaging conditions to obtain an initial imaging result; using the seismic wave forward wave field, the seismic wave reverse wave field and the initial imaging result as input, performing multiple wave imaging using an iterative deconvolution imaging method to obtain an iterative deconvolution multiple wave imaging result; during the multiple wave imaging process, determining a deconvolution imaging objective function value; if the objective function value does not meet a preset condition, iterating the imaging result until the deconvolution imaging objective function value meets the preset condition, and generating a final iterative deconvolution multiple wave imaging result based on the iterative imaging result. The present invention proposes an objective function of regularization constraints constrained by space and frequency domains, and provides an iterative solution method. When performing multiple wave migration, the frequency band range of the migration imaging can be increased as much as possible, and the imaging quality can be improved.

[0109] In order to explain the multiple wave migration imaging method based on iterative deconvolution imaging conditions more clearly, each step is described in detail below.

[0110] In one embodiment, step 101, obtaining a seismic wave forward propagation wave field, includes:

[0111] Get input data; in forward propagation, the input data is seismic full-wave field shot record data containing direct wave, primary reflection wave and free surface multiple wave information.

[0112] According to the input data, a one-way wave equation algorithm is used to perform forward propagation of the wave field to obtain the seismic wave forward propagation wave field.

[0113] The forward wavelength of seismic waves can be represented by D.

[0114] In one embodiment, step 101, obtaining a seismic wave back propagation wave field, includes:

[0115] Get input data; in reverse propagation, the input data is multiple wave information or seismic full wave field shot records;

[0116] According to the input data, a one-way wave equation algorithm is used to perform reverse propagation of the wave field to obtain the seismic wave reverse propagation wave field.

[0117] The wavelength of the seismic wave can be represented by U.

[0118] In one embodiment, step 102, multiple wave imaging is performed on the forward wave field and the reverse wave field of the seismic wave using deconvolution imaging conditions to obtain an initial imaging result R 0 ;

[0119] In one embodiment, step 103, the seismic wave forward wave field, the seismic wave reverse wave field and the initial imaging result R 0 As input, multiple wave imaging is performed using an iterative deconvolution imaging method to obtain an iterative deconvolution multiple wave imaging result;

[0120] In one embodiment, step 104, determining a deconvolution imaging objective function value during multiple wave imaging, includes:

[0121] The relational expression of the deconvolution imaging objective function is:

[0122]

[0123] Where, E(R) is the objective function value;

[0124] U is the back propagation wave field of seismic waves;

[0125] D is the positive wave field of seismic waves;

[0126] R is the imaging result;

[0127] Indicates obtaining the spatial derivative of the imaging result;

[0128] a and β are dimensionless parameters;

[0129] ω is the circular frequency;

[0130] It means to obtain the frequency derivative of the imaging result.

[0131] In one embodiment, in step 105, if the objective function value does not meet the preset condition, the imaging result is iterated until the deconvolution imaging objective function value meets the preset condition, and a final iterative deconvolution multiple wave imaging result is generated according to the iterated imaging result, including:

[0132] The formula for iterative deconvolution imaging conditions is:

[0133]

[0134] Among them, R n+1 is the imaging result after the n+1th iteration;

[0135] R n is the imaging result after the nth iteration;

[0136] i, j, k, l are the spatial and frequency coordinates (x, y, z, ω) of the imaging results; i, j, k correspond to the spatial x, y, z, and l corresponds to the frequency ω;

[0137] It is the complex conjugate of the forward propagation wave field of the seismic wave.

[0138] According to the imaging result R after iteration n+1 , recalculate the objective function value; wherein the objective function value gradually approaches 0 during calculation;

[0139] When the deconvolution imaging objective function value meets the preset condition, a final iterative deconvolution multiple wave imaging result is generated according to the imaging result after iteration.

[0140] In order to explain the multiple wave migration imaging method based on iterative deconvolution imaging conditions more clearly, a detailed description is given below in conjunction with a specific embodiment.

[0141] In the context of the prior art, the objective function E(R) that needs to be solved for multiple wave migration under the traditional deconvolution imaging condition is:

[0142]

[0143] The multiple wave imaging under the traditional deconvolution imaging condition is to solve formula (1) and derive the following specific imaging formula:

[0144]

[0145] in, For imaging results;

[0146] The existing industrial applications of multiple wave imaging mainly use the multiple wave migration method under the deconvolution imaging condition. The patent of this invention is different from the multiple wave imaging under the traditional deconvolution imaging condition. It introduces regularization terms in the frequency and space methods, and realizes a wider frequency of the multiple wave imaging results and a clearer imaging result. Specifically, compared with the objective function of the traditional deconvolution imaging condition, the present invention proposes a new deconvolution imaging objective function:

[0147]

[0148] Formula (3) cannot use the traditional method of obtaining gradients to separate the imaging results. It is necessary to use the optical flow method to solve the objective function. The calculation method is:

[0149]

[0150] Formula (4) is the formula for iterative deconvolution imaging conditions. By using it to iteratively solve the imaging results, iterative deconvolution multiple wave imaging can be achieved.

[0151] The multiple wave migration imaging method based on iterative deconvolution imaging conditions proposed in the present invention has at least the following advantages compared with the prior art:

[0152] 1. Iterative deconvolution imaging conditions have constraints in both the spatial and frequency directions of imaging, and the imaging is more controllable and flexible.

[0153] 2. Iterative deconvolution imaging conditions are an extension of conventional deconvolution imaging methods, and their physical basis is more complete.

[0154] 3. Iterative deconvolution imaging conditions can make the frequency of imaging results wider and achieve better imaging resolution.

[0155] Combined with a practical scenario, reference Figure 2 , is a schematic diagram of the multiple wave imaging process of iterative deconvolution conditions according to a specific embodiment of the present invention. Figure 2 As shown, first, the seismic shot record is input; the input data in the forward propagation is the seismic full-wavefield shot record data containing the direct wave, the primary reflection wave and the free surface multiple wave information; the input data in the reverse propagation can be only the multiple wave information or the seismic full-wavefield shot record.

[0156] Then, using the input shot records, the one-way wave equation algorithm is used to perform forward and reverse propagation of the wave field to obtain the forward propagation wave field D and the reverse propagation wave field U. The conventional deconvolution imaging conditions are used to perform multiple wave imaging on the forward and reverse propagation wave fields, and the imaging result R is calculated. 0 .

[0157] Then the wave fields U, D and the initial imaging results R 0 As input, the iterative deconvolution imaging formula is used to perform multiple wave imaging. After n iterations, the iterative deconvolution multiple wave imaging result R can be obtained. n After n+1 iterations, the iterative deconvolution multiple wave imaging result R can be obtained. n+1 .

[0158] To verify the effectiveness of this method, the full wavefield seismic records of a single-layer model ( Figure 3A ) and multiple wave records ( Figure 3B ) as input, the multiple wave imaging under the cross-correlation imaging condition, the conventional deconvolution imaging condition and the iterative deconvolution imaging condition is calculated, and the calculation results are shown in 4A to 4F As shown, Figure 4A , 4C , 4E, the horizontal axis range is 0-135, and the vertical axis range is 0-123; the corresponding partial enlargement diagram ( Figure 4B , Figure 4D , Figure 4F )The horizontal axis range is 56-77, and the vertical axis range is 0-100.

[0159] The spectra of the three results are further compared, such as Figure 5As shown in the figure, the iterative deconvolution imaging result has the widest spectrum and the highest resolution, which proves the effectiveness of the method. To further verify the stability of the method, the Marmousi model is used for VSP seismic acquisition, as shown in Figure 6 As shown in Figure 2, the model is relatively complex, and relatively complex common detection point gathers are received at different depths, such as 7A to 7C shown.

[0160] Finally, conventional deconvolution imaging and iterative deconvolution imaging were used to realize multiple wave migration. The results are shown in Figure 2. Fig. 8A and Figure 8B As shown in the figure, it can be seen from the comparison that the resolution and lateral continuity of the iterative deconvolution imaging results have been significantly improved, which illustrates the advantages and effectiveness of the method of the present application.

[0161] It should be noted that, although the operations of the method of the present invention are described in a specific order in the above embodiments and the accompanying drawings, this does not require or imply that the operations must be performed in the specific order, or that all the operations shown must be performed to achieve the desired results. Additionally or alternatively, some steps may be omitted, multiple steps may be combined into one step, and / or one step may be decomposed into multiple steps.

[0162] After introducing the method of the exemplary embodiment of the present invention, next, refer to Fig. 9 A multiple wave migration imaging device based on iterative deconvolution imaging conditions according to an exemplary embodiment of the present invention is introduced.

[0163] The implementation of the multiple wave migration imaging device based on iterative deconvolution imaging conditions can refer to the implementation of the above method, and the repeated parts will not be repeated. The term "module" or "unit" used below can be a combination of software and / or hardware that implements a predetermined function. Although the device described in the following embodiments is preferably implemented in software, the implementation of hardware, or a combination of software and hardware, is also possible and conceivable.

[0164] Based on the same inventive concept, the present invention also proposes a multiple wave migration imaging device based on iterative deconvolution imaging conditions, such as Fig. 9 As shown, the device comprises:

[0165] A data acquisition module 910 is used to acquire a seismic wave forward propagation wave field and a seismic wave reverse propagation wave field;

[0166] The multiple wave imaging processing module 920 is used to perform multiple wave imaging on the seismic wave forward propagation wave field and the seismic wave reverse propagation wave field by using deconvolution imaging conditions to obtain an initial imaging result;

[0167] Taking the seismic wave forward wave field, the seismic wave reverse wave field and the initial imaging result as input, performing multiple wave imaging by using an iterative deconvolution imaging method, and obtaining an iterative deconvolution multiple wave imaging result;

[0168] In the process of multiple wave imaging, the objective function value of deconvolution imaging is determined;

[0169] If the objective function value does not meet the preset condition, the imaging result is iterated until the deconvolution imaging objective function value meets the preset condition, and the final iterative deconvolution multiple wave imaging result is generated according to the iterative imaging result.

[0170] In one embodiment, the data acquisition module 910 is specifically used for:

[0171] Get input data; in forward propagation, the input data is seismic full-wave field shot record data containing direct wave, primary reflection wave and free surface multiple wave information.

[0172] According to the input data, a one-way wave equation algorithm is used to perform forward propagation of the wave field to obtain the seismic wave forward propagation wave field.

[0173] In one embodiment, the data acquisition module 910 is specifically used for:

[0174] Get input data; in reverse propagation, the input data is multiple wave information or seismic full wave field shot records;

[0175] According to the input data, a one-way wave equation algorithm is used to perform reverse propagation of the wave field to obtain the seismic wave reverse propagation wave field.

[0176] In one embodiment, the multiple wave imaging processing module 920 determines the deconvolution imaging objective function value during the multiple wave imaging process, including:

[0177] The relationship of the deconvolution imaging objective function is:

[0178]

[0179] Where E(R) is the objective function value;

[0180] U is the back propagation wave field of seismic waves;

[0181] D is the positive wave field of seismic waves;

[0182] R is the imaging result;

[0183] Indicates obtaining the spatial derivative of the imaging result;

[0184] a and β are dimensionless parameters;

[0185] ω is the circular frequency;

[0186] It means to obtain the frequency derivative of the imaging result.

[0187] In one embodiment, if the objective function value does not satisfy a preset condition, the multiple wave imaging processing module 920 iterates the imaging result until the deconvolution imaging objective function value satisfies the preset condition, and generates a final iterative deconvolution multiple wave imaging result according to the iterated imaging result, including:

[0188] The formula for iterative deconvolution imaging conditions is:

[0189]

[0190] Among them, R n+1 is the imaging result after the n+1th iteration;

[0191] R n is the imaging result after the nth iteration;

[0192] i, j, k, and l are the spatial and frequency coordinates of the imaging results (x, y, z, ω);

[0193] It is the complex conjugate of the forward propagation wave field of the seismic wave.

[0194] In one embodiment, the multiple wave imaging processing module 920 is specifically used for:

[0195] According to the imaging result R after iteration n+1 , recalculate the objective function value; wherein the objective function value gradually approaches 0 during calculation;

[0196] When the deconvolution imaging objective function value meets the preset condition, a final iterative deconvolution multiple wave imaging result is generated according to the imaging result after iteration.

[0197] It should be noted that, although several modules of the multiple wave migration imaging device based on iterative deconvolution imaging conditions are mentioned in the above detailed description, such division is only exemplary and not mandatory. In fact, according to an embodiment of the present invention, the features and functions of two or more modules described above can be embodied in one module. Conversely, the features and functions of one module described above can be further divided into multiple modules for embodiment.

[0198] Based on the above invention concept, Fig.10As shown, the present invention also proposes a computer device 1000, including a memory 1010, a processor 1020, and a computer program 1030 stored in the memory 1010 and executable on the processor 1020, wherein the processor 1020 implements the aforementioned multiple wave migration imaging method based on iterative deconvolution imaging conditions when executing the computer program 1030.

[0199] Based on the aforementioned inventive concept, the present invention proposes a computer-readable storage medium, wherein the computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the multiple-wave migration imaging method based on iterative deconvolution imaging conditions is implemented.

[0200] Based on the aforementioned inventive concept, the present invention proposes a computer program product, wherein the computer program product comprises a computer program, and when the computer program is executed by a processor, a multiple wave migration imaging method based on iterative deconvolution imaging conditions is implemented.

[0201] In the implementation process of the multiple wave migration imaging method and device based on iterative deconvolution imaging conditions proposed by the present invention, the iterative deconvolution imaging conditions have constraints in the imaging space and frequency directions, and the imaging controllable flexibility is higher. The iterative deconvolution imaging conditions are an extension of the conventional deconvolution imaging method, and the physical basis is more complete. The iterative deconvolution imaging conditions can make the frequency expansion of the imaging results wider and achieve better imaging resolution, providing strong technical support for oil exploration.

[0202] The acquisition, storage, use, and processing of data in the technical solution of this application comply with the relevant provisions of laws and regulations.

[0203] It will be appreciated by those skilled in the art that embodiments of the present invention may be provided as methods, devices, or computer program products. Therefore, the present invention may take the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware. Moreover, the present invention may take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

[0204] The present invention is described with reference to flowcharts and / or block diagrams of methods and computer program products according to embodiments of the present invention. It should be understood that each process and / or block in the flowchart and / or block diagram, as well as the combination of processes and / or blocks in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to generate a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowchart and / or block diagram. Figure 1 A process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0205] These computer program instructions may also be stored in a computer-readable memory capable of directing a computer or other programmable data processing device to operate in a specific manner, so that the instructions stored in the computer-readable memory produce an article of manufacture comprising an instruction device, which implements the process Figure 1 A process or multiple processes and / or boxes Figure 1 A function specified in one or more boxes.

[0206] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operating steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing instructions for implementing the process. Figure 1 A process or multiple processes and / or boxes Figure 1 The steps for the functions specified in one or more boxes.

[0207] Finally, it should be noted that the above-described embodiments are only specific implementations of the present invention, which are used to illustrate the technical solutions of the present invention, rather than to limit them. The protection scope of the present invention is not limited thereto. Although the present invention is described in detail with reference to the above-described embodiments, ordinary technicians in the field should understand that any technician familiar with the technical field can still modify the technical solutions recorded in the above-described embodiments within the technical scope disclosed by the present invention, or can easily think of changes, or make equivalent replacements for some of the technical features therein; and these modifications, changes or replacements do not make the essence of the corresponding technical solutions deviate from the spirit and scope of the technical solutions of the embodiments of the present invention, and should be included in the protection scope of the present invention. Therefore, the protection scope of the present invention shall be subject to the protection scope of the claims.

Claims

1. A multiple wave migration imaging method based on iterative deconvolution imaging conditions, It is characterized in that include: Obtaining seismic wave forward wave field and seismic wave reverse wave field; Performing multiple wave imaging on the forward wave field of the seismic wave and the reverse wave field of the seismic wave using deconvolution imaging conditions to obtain an initial imaging result; Taking the seismic wave forward wave field, the seismic wave reverse wave field and the initial imaging result as input, performing multiple wave imaging by using an iterative deconvolution imaging method, and obtaining an iterative deconvolution multiple wave imaging result; In the process of multiple wave imaging, the objective function value of deconvolution imaging is determined; If the objective function value does not meet the preset condition, the imaging result is iterated until the deconvolution imaging objective function value meets the preset condition, and the final iterative deconvolution multiple wave imaging result is generated according to the iterative imaging result.

2. The multiple wave migration imaging method based on iterative deconvolution imaging conditions according to claim 1, It is characterized in that Obtain seismic wave forward wave field, including: Get input data; in forward propagation, the input data is seismic full-wave field shot record data containing direct wave, primary reflection wave and free surface multiple wave information. According to the input data, a one-way wave equation algorithm is used to perform forward propagation of the wave field to obtain the seismic wave forward propagation wave field.

3. The multiple wave migration imaging method based on iterative deconvolution imaging conditions according to claim 1, It is characterized in that Obtain seismic wave back propagation field, including: Get input data; in reverse propagation, the input data is multiple wave information or seismic full wave field shot records; According to the input data, a one-way wave equation algorithm is used to perform reverse propagation of the wave field to obtain the seismic wave reverse propagation wave field.

4. The multiple wave migration imaging method based on iterative deconvolution imaging conditions according to claim 1, It is characterized in that In the process of multiple wave imaging, the deconvolution imaging objective function value is determined, including: The relational expression of the deconvolution imaging objective function is: Where, E(R) is the objective function value; U is the back propagation wave field of seismic waves; D is the positive wave field of seismic waves; R is the imaging result; Indicates obtaining the spatial derivative of the imaging result; a and β are dimensionless parameters; ω is the circular frequency; It means to obtain the frequency derivative of the imaging result.

5. The multiple wave migration imaging method based on iterative deconvolution imaging conditions according to claim 4, It is characterized in that If the objective function value does not meet the preset condition, the imaging result is iterated until the deconvolution imaging objective function value meets the preset condition, and a final iterative deconvolution multiple wave imaging result is generated according to the iterative imaging result, including: The formula for iterative deconvolution imaging conditions is: Among them, R n+1 is the imaging result after the n+1th iteration; R n is the imaging result after the nth iteration; i, j, k, and l are the spatial and frequency coordinates of the imaging results (x, y, z, ω); It is the complex conjugate of the forward propagation wave field of the seismic wave.

6. The multiple wave migration imaging method based on iterative deconvolution imaging conditions according to claim 5, It is characterized in that If the objective function value does not meet the preset condition, the imaging result is iterated until the deconvolution imaging objective function value meets the preset condition, and a final iterative deconvolution multiple wave imaging result is generated according to the iterative imaging result, including: According to the imaging result R after iteration n+1 , recalculate the objective function value; wherein the objective function value gradually approaches 0 during calculation; When the deconvolution imaging objective function value meets the preset condition, a final iterative deconvolution multiple wave imaging result is generated according to the imaging result after iteration.

7. A multiple wave migration imaging device based on iterative deconvolution imaging conditions, It is characterized in that include: A data acquisition module, used to acquire the forward wave field of seismic waves and the reverse wave field of seismic waves; A multiple wave imaging processing module, used for performing multiple wave imaging on the seismic wave forward propagation wave field and the seismic wave reverse propagation wave field by using deconvolution imaging conditions to obtain an initial imaging result; Taking the seismic wave forward wave field, the seismic wave reverse wave field and the initial imaging result as input, performing multiple wave imaging by using an iterative deconvolution imaging method, and obtaining an iterative deconvolution multiple wave imaging result; In the process of multiple wave imaging, the objective function value of deconvolution imaging is determined; If the objective function value does not meet the preset condition, the imaging result is iterated until the deconvolution imaging objective function value meets the preset condition, and the final iterative deconvolution multiple wave imaging result is generated according to the iterative imaging result.

8. The multiple wave migration imaging device based on iterative deconvolution imaging conditions according to claim 7, It is characterized in that The data acquisition module is specifically used for: Get input data; in forward propagation, the input data is seismic full-wave field shot record data containing direct wave, primary reflection wave and free surface multiple wave information. According to the input data, a one-way wave equation algorithm is used to perform forward propagation of the wave field to obtain the seismic wave forward propagation wave field.

9. The multiple wave migration imaging device based on iterative deconvolution imaging conditions according to claim 7, It is characterized in that The data acquisition module is specifically used for: Get input data; in reverse propagation, the input data is multiple wave information or seismic full wave field shot records; According to the input data, a one-way wave equation algorithm is used to perform reverse propagation of the wave field to obtain the seismic wave reverse propagation wave field.

10. The multiple wave migration imaging device based on iterative deconvolution imaging conditions according to claim 7, It is characterized in that The multiple wave imaging processing module determines the deconvolution imaging objective function value during the multiple wave imaging process, including: The relational expression of the deconvolution imaging objective function is: Where, E(R) is the objective function value; U is the back propagation wave field of seismic waves; D is the positive wave field of seismic waves; R is the imaging result; Indicates obtaining the spatial derivative of the imaging result; a and β are dimensionless parameters; ω is the circular frequency; It means to obtain the frequency derivative of the imaging result.

11. The multiple wave migration imaging device based on iterative deconvolution imaging conditions according to claim 10, It is characterized in that The multiple wave imaging processing module, if the objective function value does not meet the preset condition, iterates the imaging result until the deconvolution imaging objective function value meets the preset condition, and generates a final iterative deconvolution multiple wave imaging result according to the iterated imaging result, including: The formula for iterative deconvolution imaging conditions is: Among them, R n+1 is the imaging result after the n+1th iteration; R n is the imaging result after the nth iteration; i, j, k, and l are the spatial and frequency coordinates of the imaging results (x, y, z, ω); It is the complex conjugate of the forward propagation wave field of the seismic wave.

12. The multiple wave migration imaging device based on iterative deconvolution imaging conditions according to claim 11, It is characterized in that The multiple wave imaging processing module is specifically used for: According to the imaging result R after iteration n+1 , recalculate the objective function value; wherein the objective function value gradually approaches 0 during calculation; When the deconvolution imaging objective function value meets the preset condition, a final iterative deconvolution multiple wave imaging result is generated according to the imaging result after iteration.

13. A computer device comprising a memory, a processor and a computer program stored in the memory and executable on the processor, It is characterized in that When the processor executes the computer program, the method according to any one of claims 1 to 6 is implemented.

14. A computer-readable storage medium, It is characterized in that The computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the method according to any one of claims 1 to 6 is implemented.

15. A computer program product, It is characterized in that The computer program product comprises a computer program, and when the computer program is executed by a processor, the method according to any one of claims 1 to 6 is implemented.

Images