An adaptive partial angle superimposed data volume interpolation hole filling method and device

Abstract

The present application relates to the field of seismic signal processing and discloses an adaptive method and device for interpolating and filling holes in partial angle stacked data volumes. The method comprises: determining seismic traces with holes, deleting the data of the seismic traces with holes from the original seismic data, and obtaining an input sample space; determining interpolated data based on the sample space and the original seismic data; determining the absolute amplitude value of the original seismic data, determining the area with an absolute amplitude value equal to a second preset threshold as a hole area, and generating a hole factor to identify whether it is a hole area; determining interpolated hole filling data for the hole area based on the hole factor and the interpolated data; determining complete interpolated hole filling data based on the original seismic data and the interpolated hole filling data; and denoising the complete interpolated hole filling data to obtain the final adaptive interpolated hole filling data. This method overcomes the boundary effect caused by holes, incomplete data, and the resulting unsatisfactory inversion effect.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of seismic signal processing, in particular to an adaptive partial angle stack data body interpolation hole filling method and device, a storage medium and an electronic equipment. BACKGROUND

[0002] The purpose of the background description provided herein is to generally present the context of the application. The statements described in this part are only provided to give a background of the application and do not necessarily constitute the prior art.

[0003] The pre-stack inversion combines the advantages of AVO (Amplitude variation with offset, AVO) inversion and post-stack wave impedance inversion, and the elastic wave impedance inversion method using multiple partial angle stack data bodies has become a research hotspot of geophysicists.

[0004] Compared with conventional stack data, partial angle stack data has its own characteristics and advantages:

[0005] Firstly, multiple partial angle stack data bodies still contain AVO information, have more information than conventional stack data, and have the potential to obtain underground elastic parameters through inversion methods;

[0006] Secondly, since multiple partial angle stack data bodies are the result of stacking partial data within a certain angle range, they have higher signal-to-noise ratio than conventional pre-stack gather data, thus helping to obtain reliable inversion results.

[0007] Unlike conventional partial offset stack data bodies, there will be holes in the data body below the first non-zero sample in the angle partial stack data body, such as Figure 2 Part a of the above figure, when the target layer is deep, the holes in the data will seriously affect the inversion results of the seismic data. The existing pre-stack inversion interpretive processing does not consider the influence of data holes, especially when there are no holes in the target layer data, the hole problem is basically ignored, but in actual inversion, even if there are no holes in the target layer, due to the boundary effect, the application effect of the inversion method and the integrity of the inversion result will also be affected.

[0008] Therefore, an adaptive partial angle stack data body interpolation hole filling method is needed to solve the above-mentioned defects. SUMMARY

[0009] To solve the above problems, the application provides an adaptive partial angle stack data body interpolation hole filling method and device, a storage medium and an electronic device. The adaptive interpolation hole filling processing method according to the data itself can adaptively monitor the holes in the data, and combines the matched pursuit Fourier interpolation method to perform hole filling processing, so that the data below the first non-zero sample point value of the seismic data is complete, thereby reducing the boundary effect of the data in seismic inversion and increasing the stability and integrity of the inversion result.

[0010] In a first aspect, the application provides an adaptive partial angle stack data body interpolation hole filling method, which comprises the following steps:

[0011] determining a seismic trace with a hole and deleting the data of the seismic trace with a hole from the original seismic data to obtain an input sample space for interpolation processing;

[0012] determining interpolation data according to the sample space and the original seismic data;

[0013] determining the absolute amplitude value of the original seismic data, determining the region with the absolute amplitude value equal to a second preset threshold as a hole region, and generating a hole factor for identifying whether it is a hole region;

[0014] determining interpolation hole filling data of the hole region according to the hole factor and the interpolation data;

[0015] determining complete interpolation hole filling data according to the original seismic data and the interpolation hole filling data;

[0016] performing denoising processing on the complete interpolation hole filling data to obtain final adaptive interpolation hole filling data.

[0017] Further, the original seismic data comprises:

[0018] partial angle stack seismic data.

[0019] Further, the determination of the seismic trace with a hole comprises:

[0020] determining the number of non-zero sample points of the seismic trace within a preset time length;

[0021] if the number of non-zero sample points is not less than a first preset threshold, determining that the seismic trace is a seismic trace with a hole.

[0022] Further, the determination of the interpolation data according to the sample space and the original seismic data comprises:

[0023] taking the position information of the original seismic data as the expected output space of the interpolation processing;

[0024] performing matching pursuit Fourier interpolation processing according to the sample space and the expected output space to obtain interpolation data.

[0025] Further, the interpolation data determination module is configured to determine interpolation data according to the sample space and the original seismic data.

[0026] Further, the interpolation data determination module is configured to determine interpolation data according to the sample space and the original seismic data.

[0027] Further, the interpolation data determination module is configured to determine interpolation data according to the sample space and the original seismic data.

[0028] Further, the interpolation data determination module is configured to determine interpolation data according to the sample space and the original seismic data.

[0029] Further, the interpolation data determination module is configured to determine interpolation data according to the sample space and the original seismic data.

[0030] Further, the interpolation data determination module is configured to determine interpolation data according to the sample space and the original seismic data.

[0031] A second aspect of the present application provides an adaptive partial angle stack data body interpolation hole filling device, the device comprising:

[0032] A sample space determination module is configured to determine seismic traces with holes and delete data of seismic traces with holes from original seismic data to obtain an input sample space for interpolation processing.

[0033] An interpolation data determination module is configured to determine interpolation data according to the sample space and the original seismic data.

[0034] A hole region determination module is configured to determine absolute amplitude values of the original seismic data, determine a region with an absolute amplitude value equal to a second preset threshold as a hole region, and generate a hole factor for identifying whether the region is a hole region.

[0035] An interpolation hole filling data determination module is configured to determine interpolation hole filling data of the hole region according to the hole factor and the interpolation data.

[0036] A complete interpolation hole filling data determination module is configured to determine complete interpolation hole filling data according to the original seismic data and the interpolation hole filling data.

[0037] A denoising processing module is configured to perform denoising processing on the complete interpolation hole filling data to obtain final adaptive interpolation hole filling data.

[0038] In a third aspect, the present application provides a computer readable storage medium storing a computer program, which can be executed by one or more processors to implement the steps of the method as described above.

[0039] In a fourth aspect, the present application provides an electronic device comprising a memory and one or more processors, wherein the memory stores a computer program, and the memory and the one or more processors are communicatively connected, and the computer program is executed by the one or more processors to implement the steps of the method as described above.

[0040] Compared with the prior art, the technical scheme of the present application has the advantages or beneficial effects as follows:

[0041] According to the adaptive interpolation hole-filling method based on data itself disclosed in the present application, the holes in the data can be adaptively monitored, and the matching pursuit Fourier interpolation method is combined to perform hole-filling processing, so that the data below the first non-zero sample point value of the seismic data is complete, thereby reducing the boundary effect of the data in seismic inversion and increasing the stability and integrity of the inversion result. The situation that the inversion effect is not ideal due to the boundary effect and incomplete data caused by the existence of holes in the data during the calculation of seismic inversion and the interpretative processing of the basic inversion data is overcome. BRIEF DESCRIPTION OF DRAWINGS

[0042] In order to more clearly illustrate the technical scheme in the embodiments of the present application or the prior art, the drawings needed to be used in the embodiments or the prior art description will be briefly introduced as follows. Obviously, the drawings in the following description are only embodiments of the present application, and those skilled in the art can obtain other drawings according to the provided drawings without any creative effort.

[0043] In addition, it should be further pointed out that, for the convenience of description, only the parts related to the present disclosure are shown in the drawings. The drawings accompanying the specification form part of the present application and are used to provide a further understanding of the present application. The schematic embodiments and their descriptions in the present application are used to explain the present application and do not constitute an improper limitation on the present application. In the drawings:

[0044] Figure 1 A flowchart of an adaptive partial-angle stack data body interpolation hole-filling method provided by the embodiments of the present application is shown in the figure;

[0045] Figure 2 A schematic diagram of the effect of an adaptive partial-angle stack data body interpolation hole-filling method provided by the embodiments of the present application is shown in the figure;

[0046] Figure 3A schematic diagram of an adaptive partial angle stack data body interpolation hole filling before and after explanatory processing effect provided by an embodiment of the present application.

[0047] Figure 4 A flowchart of another adaptive partial angle stack data body interpolation hole filling method provided by an embodiment of the present application.

[0048] Figure 5 A structural schematic diagram of an adaptive partial angle stack data body interpolation hole filling device provided by an embodiment of the present application.

[0049] Figure 6 A connection block diagram of an electronic device provided by an embodiment of the present application. DETAILED DESCRIPTION

[0050] The embodiments of the present application will be described in detail below with reference to the accompanying drawings and embodiments, so that how the present application applies technical means to solve technical problems and achieves corresponding technical effects can be fully understood and implemented. The embodiments of the present application and each feature in the embodiments can be combined with each other on the premise of no conflict, and the formed technical solutions are all within the protection scope of the present application.

[0051] It should be clear that the embodiments described below are only part of the embodiments of the present application, not all the embodiments. Based on the embodiments in the present application, all other embodiments obtained by those skilled in the art without creative labor are within the scope of protection of the present application.

[0052] As known from the background art, pre-stack inversion combines the advantages of AVO (Amplitude variation with offset, AVO for short) inversion and post-stack wave impedance inversion, and the elastic wave impedance inversion method using multiple partial angle stack data bodies becomes a research hotspot of geophysicists.

[0053] Compared with conventional stack data, partial angle stack data has its own characteristics and advantages:

[0054] Firstly, multiple partial angle stack data bodies still contain AVO information, have more information than conventional stack data, and have the potential to obtain underground elastic parameters through inversion methods;

[0055] Secondly, since multiple partial angle stack data bodies are the result of stacking partial data within a certain angle range, they have higher signal-to-noise ratio than conventional pre-stack gather data, thus helping to obtain reliable inversion results.

[0056] Different from the conventional offset-based partial stacking data volume, the angle-based partial stacking data volume will have holes in the data volume below the first non-zero sample point, such as Figure 2 In part a of Figure 1, when the target layer is deep, holes in the data can seriously affect the inversion results of the seismic data. Existing prestack inversion interpretive processing does not consider the impact of holes in the data. In particular, when the target layer data does not have holes, the hole problem is basically ignored. However, in actual inversion, even if there are no holes in the target layer, boundary effects can still affect the application of the inversion method and the integrity of the inversion results.

[0057] In view of this, the present application proposes an adaptive method for interpolation and filling holes in partial angle stacked data. According to the adaptive interpolation and filling hole processing method based on the data itself disclosed in the present application, holes in the data can be adaptively monitored, and the matching pursuit Fourier interpolation method can be combined to perform hole filling processing, so that the data below the first non-zero sample value of the seismic data is complete, so as to reduce the boundary effect of the data during seismic inversion, thereby increasing the stability and integrity of the inversion results. It overcomes the situation in which the inversion effect is not ideal due to the boundary effect caused by holes in the data and the incomplete data during seismic inversion calculation and explanatory processing of the inversion basic data.

[0058] Example 1

[0059] This embodiment provides an adaptive method for filling holes by interpolating data volumes with partial angle superposition. Figure 1 This is a flow chart of an adaptive partial angle superposition data volume interpolation method for filling holes provided in an embodiment of the present application, such as Figure 1 As shown, the method disclosed in this embodiment includes the following steps:

[0060] Step 110: Determine the seismic traces with voids, and delete the data of the seismic traces with voids from the original seismic data to obtain an input sample space for interpolation processing;

[0061] Step 120: determining interpolation data according to the sample space and the original seismic data;

[0062] Step 130: determining an absolute amplitude value of the original seismic data, determining a region where the absolute amplitude value is equal to a second preset threshold as a cavity region, and generating a cavity factor for identifying whether the region is a cavity region;

[0063] Step 140: Determine interpolation hole filling data for the hole area according to the hole factor and the interpolation data;

[0064] Step 150: determining complete interpolation void-filling data based on the original seismic data and the interpolation void-filling data;

[0065] Step 160, denoising the complete interpolation data for filling the holes to obtain final adaptive interpolation data for filling the holes.

[0066] In some embodiments, the original seismic data comprises:

[0067] Partial angle stack seismic data.

[0068] In the present embodiment, the purpose is to fill the hole part in the partial angle stack data body, to provide a regular and complete basis data for eliminating the hard boundary of the hole, and to obtain a good inversion effect and the application of subsequent inversion and other explanatory processing methods.

[0069] In some embodiments, the determining of the seismic trace with the hole comprises:

[0070] Determining the number of non-zero sample values of the seismic trace within a preset time length;

[0071] In the case that the number of non-zero sample values is not less than a first preset threshold, the seismic trace is determined as the seismic trace with the hole.

[0072] As an example, the number of non-zero sample values of the seismic trace within a given time window is counted, and a first preset threshold is set. When the number of non-zero sample values within the given time window is greater than the set first preset threshold, the seismic trace is defined as the seismic trace with the hole, i.e. the trace that needs to be interpolated.

[0073] It should be noted that the first preset threshold can be set according to actual needs and / or experience.

[0074] In some embodiments, the determining of the interpolation data according to the sample space and the original seismic data comprises:

[0075] Taking the position information of the original seismic data as the expected output space of the interpolation processing;

[0076] Performing matching pursuit Fourier interpolation processing according to the sample space and the expected output space to obtain the interpolation data.

[0077] As an example, taking the position information of the original data as the expected output space of the interpolation processing, performing matching pursuit Fourier interpolation processing on the input sample data and the expected output space to obtain the interpolation data after calculation and processing.

[0078] In some embodiments, the determining of the interpolation data for filling the holes of the hole region according to the hole factor and the interpolation data comprises:

[0079] Taking the product of the hole factor and the interpolation data as the interpolation data for filling the holes of the hole region.

[0080] As an example, first, absolute amplitude values of the original seismic data are determined, and then regions with absolute amplitude values equal to a second preset threshold are taken as cavity regions, and regions with absolute amplitude values not equal to the second preset threshold are taken as non-cavity regions. Wherein, a cavity factor of the cavity region is represented by 1, and a cavity factor of the non-cavity region is represented by 0.

[0081] Further, a product obtained by multiplying the cavity factor with the interpolation data is taken as interpolation cavity filling data of the cavity region.

[0082] Optionally, the second preset threshold can be set as 0.

[0083] In some embodiments, the determining of the complete interpolation cavity filling data according to the original seismic data and the interpolation cavity filling data comprises:

[0084] Taking a sum of the original seismic data and the interpolation cavity filling data as the complete interpolation cavity filling data.

[0085] As an example, the original data and the interpolation cavity filling data of the cavity region are added to obtain the complete interpolation cavity filling data.

[0086] In some embodiments, the de-noising processing of the complete interpolation cavity filling data comprises:

[0087] The complete interpolation cavity filling data is subjected to matched pursuit Fourier interpolation de-noising processing.

[0088] As an example, the complete interpolation cavity filling data is subjected to matched pursuit Fourier interpolation de-noising processing to obtain final adaptive interpolation cavity filling data.

[0089] The adaptive partial angle stacking data volume interpolation method for filling holes provided in this embodiment can adaptively monitor holes in the data and perform hole filling processing in combination with the matching pursuit Fourier interpolation method, so that the data below the first non-zero sample value of the seismic data is complete, thereby reducing the boundary effect of the data during seismic inversion and increasing the stability and integrity of the inversion results. Specifically, the method includes the following steps: Step 110, determining the seismic traces with holes, and deleting the data of the seismic traces with holes from the original seismic data to obtain an input sample space for interpolation processing; Step 120, determining interpolation data based on the sample space and the original seismic data; Step 130, determining the absolute amplitude value of the original seismic data, determining the area where the absolute amplitude value is equal to a second preset threshold as a hole area, and generating a hole factor for identifying whether it is a hole area; Step 140, determining interpolation hole-filling data for the hole area based on the hole factor and the interpolation data; Step 150, determining complete interpolation hole-filling data based on the original seismic data and the interpolation hole-filling data; Step 160, denoising the complete interpolation hole-filling data to obtain final adaptive interpolation hole-filling data. The data-based adaptive interpolation hole-filling method disclosed in this application can adaptively monitor holes in the data and, combined with the matching pursuit Fourier interpolation method, fill the holes. This ensures that the data below the first non-zero sample value in the seismic data is complete, thereby reducing the boundary effects of the data during seismic inversion and increasing the stability and integrity of the inversion results. This overcomes the problem of unsatisfactory inversion results caused by boundary effects and incomplete data caused by holes in the data during seismic inversion calculations and interpretive processing of inversion basic data.

[0090] Example 2

[0091] Based on the first embodiment, this embodiment further illustrates the adaptive hole filling method of partial angle stacking data volume interpolation disclosed in the first embodiment.

[0092] As an example, see Figure 1 , Figure 1 This is a flow chart of an adaptive partial angle superposition data volume interpolation method for filling holes provided in an embodiment of the present application, such as Figure 1 As shown, the method disclosed in this embodiment includes the following steps:

[0093] Step 110: Determine the seismic traces with voids, and delete the data of the seismic traces with voids from the original seismic data to obtain an input sample space for interpolation processing.

[0094] In some embodiments, the raw seismic data includes:

[0095] Seismic data with partial angle stacking.

[0096] In the embodiment, the purpose is to fill the hollow part in the partial angle superposition data body, to provide a complete and regular data without hollow hard boundary for obtaining good inversion effect and subsequent inversion and application of other explanatory processing methods.

[0097] In some embodiments, the determining of the seismic trace with the hollow includes:

[0098] Determining the number of non-zero sample values of the seismic trace within a preset time length;

[0099] In a case where the number of non-zero sample values is not less than a first preset threshold, the seismic trace is determined as the seismic trace with the hollow.

[0100] As an example, the number of non-zero sample values of the seismic trace within a given time window is counted, and a first preset threshold is set. When the number of non-zero sample values within the given time window is greater than the set first preset threshold, the seismic trace is defined as the seismic trace with the hollow, i.e., the trace that needs to be interpolated.

[0101] It should be noted that the first preset threshold can be set according to actual requirements and / or experience.

[0102] Step 120, determining interpolation data according to the sample space and the original seismic data.

[0103] In some embodiments, the determining of the interpolation data according to the sample space and the original seismic data includes:

[0104] Taking the position information of the original seismic data as an expected output space of the interpolation processing;

[0105] Performing matched pursuit Fourier interpolation processing according to the sample space and the expected output space to obtain the interpolation data.

[0106] As an example, taking the position information of the original data as the expected output space of the interpolation processing, performing matched pursuit Fourier interpolation processing on the input sample data and the expected output space to obtain the interpolation data after calculation and processing.

[0107] Step 130, determining an absolute amplitude value of the original seismic data, determining a region with an absolute amplitude value equal to a second preset threshold as a hollow region, and generating a hollow factor for identifying whether the region is the hollow region.

[0108] Step 140, determining interpolation hole-filling data of the hollow region according to the hollow factor and the interpolation data.

[0109] In some embodiments, the determining the interpolation hole-filling data of the hole region according to the hole factor and the interpolation data comprises:

[0110] multiplying the hole factor and the interpolation data to obtain a product as the interpolation hole-filling data of the hole region.

[0111] As an example, first, the absolute amplitude value of the original seismic data is determined, then the region with the absolute amplitude value equal to the second preset threshold is taken as the hole region, and the region with the absolute amplitude value not equal to the second preset threshold is taken as the non-hole region. Wherein, the hole factor of the hole region is represented by 1, and the hole factor of the non-hole region is represented by 0.

[0112] Further, the product obtained by multiplying the hole factor and the interpolation data is taken as the interpolation hole-filling data of the hole region.

[0113] Optionally, the second preset threshold can be set to 0.

[0114] Step 150, determining the complete interpolation hole-filling data according to the original seismic data and the interpolation hole-filling data.

[0115] In some embodiments, the determining the complete interpolation hole-filling data according to the original seismic data and the interpolation hole-filling data comprises:

[0116] adding the original seismic data and the interpolation hole-filling data to obtain the complete interpolation hole-filling data.

[0117] As an example, the original data is added to the interpolation hole-filling data of the hole region to obtain the complete interpolation hole-filling data.

[0118] Step 160, performing denoising processing on the complete interpolation hole-filling data to obtain the final adaptive interpolation hole-filling data.

[0119] In some embodiments, the denoising processing on the complete interpolation hole-filling data comprises:

[0120] performing matching pursuit Fourier interpolation denoising processing on the complete interpolation hole-filling data.

[0121] As an example, the complete interpolation hole-filling data is subjected to matching pursuit Fourier interpolation denoising processing to obtain the final adaptive interpolation hole-filling data.

[0122] Further, the disclosed adaptive partial angle stack data body interpolation hole-filling method is applied to the interpolation processing of the partial angle stack data body of a certain work area, and the effect can be referred to Figure 2 and Figure 3 . Wherein, Figure 2The part a in the figure is a data body before interpolation and hole filling of the adaptive partial angle stack data body, and there is an obvious hard boundary at the hole; Figure 2 The part b in the figure is a data body after interpolation and hole filling, and it can be seen from the figure that the interpolated data at the hole has good correlation with the surrounding data, and there is no obvious signal-to-noise ratio and energy difference; Figure 2 The part c in the figure is difference data before and after interpolation and hole filling, and it can be seen from the difference data that the change of the other effective signal data which is not a hole is small, and there is basically no residual between before and after processing. Figure 3 The part a in the figure, Figure 3 The part b in the figure is a data body before and after hole filling of the adaptive partial angle stack data body, Figure 3 The part c in the figure, Figure 3 The part d in the figure is data after explanatory processing before and after hole filling, and it can be seen from the figure that the hole filling processing can effectively reduce the error problem caused by the boundary effect when other methods are applied subsequently.

[0123] The adaptive partial angle stack data body interpolation and hole filling method provided by the embodiment can adaptively monitor the holes in the data, and combines the matched pursuit Fourier interpolation method to perform hole filling processing, so that the data below the first non-zero sample value of the seismic data is complete, thereby reducing the boundary effect of the data in seismic inversion and increasing the stability and integrity of the inversion result. Specifically, the method comprises the following steps: step 110, determining a seismic trace with a hole, and deleting the data of the seismic trace with the hole from the original seismic data to obtain an input sample space for interpolation processing; step 120, determining interpolation data according to the sample space and the original seismic data; step 130, determining the absolute amplitude value of the original seismic data, determining a hole region with an absolute amplitude value equal to a second preset threshold, and generating a hole factor for identifying whether it is a hole region; step 140, determining interpolation and hole filling data of the hole region according to the hole factor and the interpolation data; step 150, determining complete interpolation and hole filling data according to the original seismic data and the interpolation and hole filling data; and step 160, performing denoising processing on the complete interpolation and hole filling data to obtain final adaptive interpolation and hole filling data. According to the adaptive interpolation and hole filling processing method disclosed in the application, the holes in the data can be adaptively monitored, and the matched pursuit Fourier interpolation method is combined to perform hole filling processing, so that the data below the first non-zero sample value of the seismic data is complete, thereby reducing the boundary effect of the data in seismic inversion and increasing the stability and integrity of the inversion result. The situation that the inversion effect is not ideal due to the boundary effect and incomplete data caused by the holes in the data during seismic inversion calculation and explanatory processing of the inversion basis data is overcome.

[0124] Embodiment three

[0125] The embodiment of the present application provides a device for adaptive partial angle stack data body interpolation hole filling. Figure 5 The structural diagram of the device provided by the embodiment of the present application is shown in the figure. Figure 5 The device 500 disclosed by the embodiment of the present application comprises:

[0126] A sample space determination module 501 is configured to determine a seismic trace with a hole, delete data of the seismic trace with the hole from original seismic data, and obtain an input sample space for interpolation processing.

[0127] An interpolation data determination module 502 is configured to determine interpolation data according to the sample space and the original seismic data.

[0128] A hole region determination module 503 is configured to determine an absolute amplitude value of the original seismic data, determine a region with an absolute amplitude value equal to a second preset threshold as a hole region, and generate a hole factor for identifying whether the region is a hole region.

[0129] An interpolation hole filling data determination module 504 is configured to determine interpolation hole filling data of the hole region according to the hole factor and the interpolation data.

[0130] A complete interpolation hole filling data determination module 505 is configured to determine complete interpolation hole filling data according to the original seismic data and the interpolation hole filling data.

[0131] A denoising processing module 506 is configured to perform denoising processing on the complete interpolation hole filling data, and obtain final adaptive interpolation hole filling data.

[0132] In some embodiments, the original seismic data comprises:

[0133] Partial angle stack seismic data.

[0134] In some embodiments, the obtaining module 501 comprises a first determination unit and a second determination unit.

[0135] The first determination unit is configured to determine a number of non-zero sample point values of a seismic trace within a preset time length.

[0136] The second determination unit is configured to determine the seismic trace as the seismic trace with the hole in a case where the number of non-zero sample point values is not less than a first preset threshold.

[0137] In some embodiments, the interpolation data determination module 502 comprises an expected output space determination unit and an interpolation processing unit.

[0138] an expected output space determination unit configured to determine an expected output space for interpolation processing according to position information of the original seismic data;

[0139] an interpolation processing unit configured to perform matched pursuit Fourier interpolation processing according to the sample space and the expected output space to obtain interpolation data.

[0140] In some embodiments, the interpolation and hole-filling data determination module 504 is configured to determine a product of the hole factor and the interpolation data as the interpolation and hole-filling data of the hole region.

[0141] In some embodiments, the complete interpolation and hole-filling data determination module 505 is configured to determine a sum of the original seismic data and the interpolation and hole-filling data as the complete interpolation and hole-filling data.

[0142] In some embodiments, the de-noising processing module 506 is configured to perform matched pursuit Fourier interpolation de-noising processing on the complete interpolation and hole-filling data.

[0143] Those skilled in the art can understand that the structure shown in the above embodiments is not a limitation of the device of the present application, and can include more or less modules / cells than the structure shown, or combine certain modules / cells, or arrange different modules / cells. Figure 5 The structure shown in the above embodiments is not a limitation of the device of the present application, and can include more or less modules / cells than the structure shown, or combine certain modules / cells, or arrange different modules / cells.

[0144] Those skilled in the art should understand that the modules or steps of the present application described above can be realized by a general computing device, which can be concentrated on a single computing device or distributed on a network composed of multiple computing devices. Alternatively, they can be realized by program codes executable by a computing device, so that they can be stored in a storage device and executed by a computing device, and in some cases, the steps shown or described can be executed in an order different from here, or they can be manufactured into individual integrated circuit modules, or multiple modules or steps among them can be manufactured into a single integrated circuit module.

[0145] Those skilled in the art can clearly understand that, for the convenience and brevity of description, the specific working process of each module of the adaptive partial angle superposition data body interpolation and hole-filling device can refer to the corresponding process in the foregoing method embodiments, and the present embodiment will not be repeated here.

[0146] Meanwhile, in order to facilitate the understanding of the technical solutions of the present application, reference can also be made to Figure 4 .

[0147] The device provided in the embodiment comprises: a sample space determination module 501 configured to determine a seismic trace with a hole and delete data of the seismic trace with a hole from original seismic data to obtain an input sample space for interpolation processing; an interpolation data determination module 502 configured to determine interpolation data according to the sample space and the original seismic data; a hole region determination module 503 configured to determine an absolute amplitude value of the original seismic data, determine a region with an absolute amplitude value equal to a second preset threshold as a hole region, and generate a hole factor used to identify whether it is a hole region; an interpolation hole filling data determination module 504 configured to determine interpolation hole filling data of the hole region according to the hole factor and the interpolation data; a complete interpolation hole filling data determination module 505 configured to determine complete interpolation hole filling data according to the original seismic data and the interpolation hole filling data; and a denoising processing module 506 configured to perform denoising processing on the complete interpolation hole filling data to obtain final adaptive interpolation hole filling data. The holes in the data can be adaptively monitored, and the matching pursuit Fourier interpolation method is combined to perform hole filling processing, so that the data below the first non-zero sample point value of the seismic data is complete, the boundary effect of the data during seismic inversion is reduced, and the stability and integrity of the inversion result are increased. The situation that the inversion effect is not ideal due to the boundary effect and the incomplete data caused by the holes in the data during the calculation of seismic inversion and the interpretative processing of the inversion basis data is overcome.

[0148] Embodiment four

[0149] The embodiment provides a computer readable storage medium. The computer readable storage medium stores a computer program. The computer program is executed by a processor to implement all or part of the steps of the method in the foregoing method embodiments:

[0150] Determine a seismic trace with a hole and delete data of the seismic trace with a hole from original seismic data to obtain an input sample space for interpolation processing;

[0151] Determine interpolation data according to the sample space and the original seismic data;

[0152] Determine an absolute amplitude value of the original seismic data, determine a region with an absolute amplitude value equal to a second preset threshold as a hole region, and generate a hole factor used to identify whether it is a hole region;

[0153] Determine interpolation hole filling data of the hole region according to the hole factor and the interpolation data;

[0154] Determine complete interpolation hole filling data according to the original seismic data and the interpolation hole filling data;

[0155] The complete interpolation hole filling data is denoised to obtain final adaptive interpolation hole filling data.

[0156] In some embodiments, the original seismic data comprises:

[0157] Partial angle stack seismic data.

[0158] In the embodiment, the partial angle stack data body is intended to fill the hole part, and to provide a regular and complete basis data for eliminating the hole hard boundary and obtaining good inversion effect and subsequent inversion and other explanatory processing method application.

[0159] In some embodiments, the determining of the seismic trace with the hole comprises:

[0160] Determining the number of non-zero sample values of the seismic trace within a preset time length.

[0161] In a case where the number of non-zero sample values is not less than a first preset threshold, the seismic trace is determined as the seismic trace with the hole.

[0162] As an example, the number of non-zero sample values of the seismic trace within a given time window is counted, and a first preset threshold is set. When the number of non-zero sample values within the given time window is greater than the set first preset threshold, the seismic trace is defined as the seismic trace with the hole, i.e., the trace that needs to be interpolated.

[0163] It should be noted that the first preset threshold can be set according to actual requirements and / or experience.

[0164] In some embodiments, the determining of the interpolation data according to the sample space and the original seismic data comprises:

[0165] The position information of the original seismic data is taken as an expected output space of the interpolation processing.

[0166] The matching pursuit Fourier interpolation processing is performed according to the sample space and the expected output space to obtain the interpolation data.

[0167] As an example, the position information of the original data is taken as an expected output space of the interpolation processing. The matching pursuit Fourier interpolation processing is performed on the input sample data and the expected output space, and the interpolation data is obtained after calculation and processing.

[0168] In some embodiments, the determining of the interpolation hole filling data of the hole region according to the hole factor and the interpolation data comprises:

[0169] The product of the hole factor and the interpolation data is taken as the interpolation hole filling data of the hole region.

[0170] As an example, first, absolute amplitude values of the original seismic data are determined, and then regions with absolute amplitude values equal to a second preset threshold are taken as the cavity regions, and regions with absolute amplitude values not equal to the second preset threshold are taken as the non-cavity regions. Wherein, a cavity factor of the cavity regions is represented by 1, and a cavity factor of the non-cavity regions is represented by 0.

[0171] Further, a product obtained by multiplying the cavity factor with the interpolation data is taken as interpolation cavity filling data of the cavity regions.

[0172] Optionally, the second preset threshold can be set as 0.

[0173] In some embodiments, the determining of the complete interpolation cavity filling data according to the original seismic data and the interpolation cavity filling data comprises:

[0174] Summing the original seismic data and the interpolation cavity filling data is taken as the complete interpolation cavity filling data.

[0175] As an example, the original data and the interpolation cavity filling data of the cavity regions are added to obtain the complete interpolation cavity filling data.

[0176] In some embodiments, the de-noising processing of the complete interpolation cavity filling data comprises:

[0177] The complete interpolation cavity filling data is subjected to matching pursuit Fourier interpolation de-noising processing.

[0178] As an example, the complete interpolation cavity filling data is subjected to matching pursuit Fourier interpolation de-noising processing to obtain final adaptive interpolation cavity filling data.

[0179] The computer readable storage medium can also include, or be, alone, a computer program, a data file, a data structure, etc., or a combination thereof. The computer readable storage medium or the computer program can be specifically designed and understood by those skilled in the computer software field, or can be known and available to those skilled in the computer software field. Examples of the computer readable storage medium include: magnetic media, such as a hard disk, a floppy disk, and a magnetic tape; optical media, such as a CD ROM disk and a DVD; a magneto-optical medium, such as an optical disk; and a hardware device specifically configured to store and execute a computer program, such as a read-only memory (ROM), a random access memory (RAM), a flash memory, or a server, an app application store, etc. Examples of the computer program include machine code (e.g., code generated by a compiler) and a file containing high-level code that can be executed by a computer by using an interpreter. The described hardware device can be configured to function as one or more software modules to perform the above-described operations and methods, and vice versa. In addition, the computer readable storage medium can be distributed in a networked computer system, and the program code or computer program can be stored and executed in a distributed manner.

[0180] Embodiment five

[0181] The embodiment provides a computer program product. The computer program product includes a computer program or instructions, which, when executed by a processor, implement all or part of the steps of the method in the foregoing method embodiments:

[0182] Determine a seismic trace with a hole, and delete data of the seismic trace with a hole from original seismic data to obtain an input sample space for interpolation processing;

[0183] Determine interpolation data according to the sample space and the original seismic data;

[0184] Determine an absolute amplitude value of the original seismic data, determine a region with an absolute amplitude value equal to a second preset threshold as a hole region, and generate a hole factor for identifying whether it is a hole region;

[0185] Determine interpolation hole-filling data of the hole region according to the hole factor and the interpolation data;

[0186] Determine complete interpolation hole-filling data according to the original seismic data and the interpolation hole-filling data;

[0187] Perform denoising processing on the complete interpolation hole-filling data to obtain final adaptive interpolation hole-filling data.

[0188] In some embodiments, the original seismic data includes:

[0189] Partial angle stack seismic data.

[0190] In the embodiment, the purpose is to fill the empty part in the partial angle stack data body, to provide a regular and complete basis data for eliminating the empty hard boundary, obtaining good inversion effect and applying subsequent inversion and other explanatory processing methods.

[0191] In some embodiments, the determining of the seismic trace with the empty part includes:

[0192] Determining the number of non-zero sample values of the seismic trace within a preset time length;

[0193] In a case where the number of non-zero sample values is not less than a first preset threshold, the seismic trace is determined as the seismic trace with the empty part.

[0194] As an example, the number of non-zero sample values of the seismic trace within a given time window is counted, and a first preset threshold is set. When the number of non-zero sample values within the given time window is greater than the set first preset threshold, the seismic trace is defined as the seismic trace with the empty part, i.e., the trace that needs to be interpolated.

[0195] It should be noted that the first preset threshold can be set according to actual needs and / or experience.

[0196] In some embodiments, the determining of the interpolation data according to the sample space and the original seismic data includes:

[0197] Taking the position information of the original seismic data as an expected output space of the interpolation processing;

[0198] Performing a matching pursuit Fourier interpolation processing according to the sample space and the expected output space to obtain the interpolation data.

[0199] As an example, the position information of the original data is taken as the expected output space of the interpolation processing, and a matching pursuit Fourier interpolation processing is performed on the input sample data and the expected output space to obtain the interpolation data after calculation and processing.

[0200] In some embodiments, the determining of the interpolation empty hole data of the empty hole region according to the empty hole factor and the interpolation data includes:

[0201] Taking the product of the empty hole factor and the interpolation data as the interpolation empty hole data of the empty hole region.

[0202] As an example, the absolute amplitude value of the original seismic data is first determined, and then the region with the absolute amplitude value equal to a second preset threshold is taken as the empty hole region, and the region with the absolute amplitude value not equal to the second preset threshold is taken as the non-empty hole region. Wherein, the empty hole factor of the empty hole region is represented by 1, and the empty hole factor of the non-empty hole region is represented by 0.

[0203] Further, the product of multiplying the hole factor and the interpolation data is taken as the interpolation hole-filling data of the hole region.

[0204] Optionally, the second preset threshold can be set as 0.

[0205] In some embodiments, the determining the complete interpolation hole-filling data according to the original seismic data and the interpolation hole-filling data comprises:

[0206] Taking the sum of the original seismic data and the interpolation hole-filling data as the complete interpolation hole-filling data.

[0207] As an example, the original data is added to the interpolation hole-filling data of the hole region to obtain the complete interpolation hole-filling data.

[0208] In some embodiments, the de-noising processing of the complete interpolation hole-filling data comprises:

[0209] The complete interpolation hole-filling data is subjected to matched pursuit Fourier interpolation de-noising processing.

[0210] As an example, the complete interpolation hole-filling data is subjected to matched pursuit Fourier interpolation de-noising processing to obtain the final adaptive interpolation hole-filling data.

[0211] Further, the computer program product can include one or more computer-executable components which, when executed by a computerized device, cause the device to carry out the steps of any of the methods disclosed herein; the computer program product can also include a computer-readable medium having stored computer-readable instructions which, when executed by a computerized device, cause the device to carry out the steps of any of the methods disclosed herein. In such embodiments, the computer program product can be downloaded and installed by a communication component from a network, and / or installed from a removable medium.

[0212] Embodiment six

[0213] The embodiment provides an electronic device. Figure 6 A connection block diagram of an electronic device provided by the embodiment of the application is shown in Figure 6 The electronic device 600 can include one or more processors 601, a memory 602, a multimedia component 603, an input / output (I / O) interface 604, and a communication component 605.

[0214] The one or more processors 601 are configured to execute all or part of the steps of the methods as previously described. The memory 602 is configured to store various types of data, which can include, for example, instructions for any of the applications or methods in the electronic device, and application-related data.

[0215] The memory 602 can be implemented by any type of volatile or nonvolatile storage devices or a combination thereof, such as static random access memory (SRAM), electrically erasable programmable read-only memory (EEPROM), erasable programmable read-only memory (EPROM), programmable read-only memory (PROM), read-only memory (ROM), magnetic storage, flash memory, magnetic or optical disk.

[0216] The one or more processors 601 can be implemented by an application specific integrated circuit (ASIC), a digital signal processor (DSP), a digital signal processing device (DSPD), a programmable logic device (PLD), a field programmable gate array (FPGA), a controller, a microcontroller, a microprocessor, or other electronic elements for performing all or part of the steps in the foregoing method embodiments:

[0217] Determine a seismic trace with a hole, and delete data of the seismic trace with a hole from the original seismic data to obtain an input sample space for interpolation processing;

[0218] Determine interpolation data according to the sample space and the original seismic data;

[0219] Determine an absolute amplitude value of the original seismic data, determine a region with an absolute amplitude value equal to a second preset threshold as a hole region, and generate a hole factor for identifying whether it is a hole region;

[0220] Determine interpolation hole filling data of the hole region according to the hole factor and the interpolation data;

[0221] Determine complete interpolation hole filling data according to the original seismic data and the interpolation hole filling data;

[0222] The complete interpolation hole filling data is denoised to obtain final adaptive interpolation hole filling data.

[0223] The multimedia component 603 can include a screen, which can be a touch screen, and an audio component for outputting and / or inputting audio signals. For example, the audio component can include a microphone for receiving external audio signals. The received audio signals can be further stored in the memory or transmitted through the communication component. The audio component also includes at least one speaker for outputting audio signals.

[0224] The I / O interface 604 provides an interface between the one or more processors 601 and other interface modules, which can be a keyboard, a mouse, a button, and the like. These buttons can be virtual buttons or physical buttons.

[0225] The communication component 605 is used for wired or wireless communication between the electronic device 600 and other devices. The wired communication includes communication through a network port, a serial port, and the like; the wireless communication includes Wi-Fi, Bluetooth, Near Field Communication (NFC), 2G, 3G, 4G, 5G, or a combination of one or more of them.

[0226] In summary, the application provides a self-adaptive partial angle superposition data body interpolation hole filling method, device, computer readable storage medium and electronic equipment. The self-adaptive partial angle superposition data body interpolation hole filling method can adaptively monitor the holes in the data, and combines the matched pursuit Fourier interpolation method to perform hole filling processing, so that the data below the first non-zero sample value of the seismic data is complete, thereby reducing the boundary effect of the data in seismic inversion, and increasing the stability and integrity of the inversion result. Specifically, the method comprises the following steps: step 110, determining the seismic trace with holes, and deleting the data of the seismic trace with holes from the original seismic data to obtain an input sample space for interpolation processing; step 120, determining interpolation data according to the sample space and the original seismic data; step 130, determining the absolute amplitude value of the original seismic data, determining the region with an absolute amplitude value equal to a second preset threshold as a hole region, and generating a hole factor for identifying whether it is a hole region; step 140, determining interpolation hole filling data of the hole region according to the hole factor and the interpolation data; step 150, determining complete interpolation hole filling data according to the original seismic data and the interpolation hole filling data; step 160, performing denoising processing on the complete interpolation hole filling data to obtain final self-adaptive interpolation hole filling data. According to the self-adaptive interpolation hole filling processing method disclosed in the application, the holes in the data can be adaptively monitored, and the matched pursuit Fourier interpolation method is combined to perform hole filling processing, so that the data below the first non-zero sample value of the seismic data is complete, thereby reducing the boundary effect of the data in seismic inversion, and further increasing the stability and integrity of the inversion result. The application overcomes the situation that the boundary effect and incomplete data caused by the holes in the data lead to unsatisfactory inversion effect during seismic inversion calculation and inversion basic data interpretation processing.

[0227] It should also be understood that, wherever implemented, the methods and / or apparatuses described herein can be implemented by other means than those shown. For example, the methods and apparatuses described herein can be implemented by one or more of: a processor executing at least one set of instructions; an application specific integrated circuit; a field programmable gate array; etc. The one or more sets of instructions can be stored on a non-transitory computer readable medium, such as a magnetic or optical disk. Further, the algorithms described herein can be implemented by a combination of hardware and software. Also, it should be understood that any logical or

[0228] In this application, the terms "include," "including," or "includes" are intended to be inclusive, in that they mean "including, but not limited to," to the extent that there is an express reference to an item. The term "in accordance with" is intended to be inclusive, in that it means "in accordance with, but not limited to," to the extent that there is an express reference to an item. The term "consisting of" is intended to be exclusive, in that it means "consisting of, but not including, unless otherwise specified. The term "consisting essentially of" is intended to be exclusive, in that it means "consisting essentially of, but not including, unless otherwise specified. The term "comprising" is intended to be inclusive, in that it means "comprising, but not limited to," to the extent that there is an express reference to an item. The terms "first," "second," "third," etc. are used to identify different elements, and are not intended to be indicative of relative importance of the elements, unless otherwise specified. The term "a number of" is intended to mean at least two, unless otherwise specified.

[0229] It should be noted that if there is a description to the server, the server can be a stand-alone physical server or terminal, or a server cluster composed of multiple physical servers, or a cloud server capable of providing cloud server, cloud database, cloud storage and CDN and other basic cloud computing services; in the present application, if there is a description to the intelligent terminal or mobile device, it should be noted that the intelligent terminal or mobile device can be a mobile phone, a tablet computer, a smart watch, a netbook, a wearable electronic device, a personal digital assistant (PDA), an augmented reality technology device (AR), a virtual reality device (VR), a smart television, a smart sound, a personal computer (PC) and the like, but is not limited thereto, and the specific form of the intelligent terminal or mobile device is not specially limited in the present application.

[0230] Finally, it should be noted that in the description of the present application, the description of the terms "one embodiment", "some embodiments", "example", "one example" or "some examples" means that the specific features, structures, materials or characteristics described in connection with the embodiment or example are included in at least one embodiment or example of the present application. In the present application, the illustrative description of the above terms does not necessarily refer to the same embodiment or example. Moreover, the specific features, structures, materials or characteristics described can be combined in any one or more embodiments or examples in a suitable manner.

[0231] Although the embodiments of the present application have been shown and described above, it should be understood that the above embodiments are exemplary, and the content described is only an implementation adopted for the purpose of facilitating the understanding of the present application, and is not intended to limit the present application. Any person skilled in the art can make any modification and change in the implementation form and details without departing from the spirit and scope of the present application, but the protection scope of the present application shall be subject to the scope defined by the appended claims.

Claims

1. An adaptive method for filling holes by interpolation of partial angle superimposed data volumes, characterized in that: The method comprises: Determine the seismic traces with voids, and delete the data of the seismic traces with voids from the original seismic data to obtain an input sample space for interpolation processing; determining interpolation data according to the sample space and the original seismic data; determining an absolute amplitude value of the original seismic data, determining a region where the absolute amplitude value is equal to a second preset threshold as a cavity region, and generating a cavity factor for identifying whether the region is a cavity region; Determining interpolation hole filling data for the hole area according to the hole factor and the interpolation data; Determining complete interpolation void-filling data based on the original seismic data and the interpolation void-filling data; The complete interpolation hole filling data is subjected to denoising processing to obtain final adaptive interpolation hole filling data.

2. The adaptive partial angle superposition data volume interpolation method for filling holes according to claim 1 is characterized in that: The original seismic data includes: Seismic data with partial angle stacking.

3. The adaptive partial angle superposition data volume interpolation method for filling holes according to claim 1 is characterized in that: The determining of the seismic traces where cavities exist includes: Determine the number of non-zero sample values ​​of the seismic trace within a preset time length; When the number of non-zero sample values ​​is not less than a first preset threshold, the seismic trace is determined to be a seismic trace with a cavity.

4. The adaptive partial angle superposition data volume interpolation method for filling holes according to claim 1 is characterized in that: The determining of interpolation data according to the sample space and the original seismic data includes: Using the position information of the original seismic data as the desired output space of the interpolation process; Matching pursuit Fourier interpolation processing is performed according to the sample space and the expected output space to obtain interpolation data.

5. The method for filling holes by interpolation of partial angle superimposed data volume according to claim 1, characterized in that: The determining, according to the hole factor and the interpolation data, the interpolation hole filling data of the hole area includes: The product of the hole factor and the interpolation data is used as the interpolation hole filling data of the hole area.

6. The adaptive partial angle superposition data volume interpolation method for filling holes according to claim 1, characterized in that: The determining of complete interpolation void-filling data based on the original seismic data and the interpolation void-filling data includes: The sum of the original seismic data and the interpolated void-filling data is used as the complete interpolated void-filling data.

7. The adaptive partial angle superposition data volume interpolation method for filling holes according to claim 1, characterized in that: The denoising process of the complete interpolation hole filling data includes: The complete interpolation hole filling data is subjected to matching pursuit Fourier interpolation denoising processing.

8. An adaptive partial angle superposition data volume interpolation hole filling device, characterized in that: include: A sample space determination module is used to determine the seismic traces with voids and delete the data of the seismic traces with voids from the original seismic data to obtain an input sample space for interpolation processing; An interpolation data determination module, configured to determine interpolation data based on the sample space and the original seismic data; a cavity region determination module, configured to determine an absolute amplitude value of the original seismic data, determine a region where the absolute amplitude value is equal to a second preset threshold as a cavity region, and generate a cavity factor for identifying whether the region is a cavity region; an interpolation hole filling data determination module, configured to determine the interpolation hole filling data of the hole area according to the hole factor and the interpolation data; A complete interpolation void filling data determination module is used to determine complete interpolation void filling data based on the original seismic data and the interpolation void filling data; The denoising processing module is used to perform denoising processing on the complete interpolation hole filling data to obtain the final adaptive interpolation hole filling data.

9. A computer-readable storage medium, characterized in that The computer program stored in the computer-readable storage medium, when executed by one or more processors, implements the adaptive partial angle stacking data volume interpolation hole filling method according to any one of claims 1 to 7.

10. An electronic device, characterized in that: The method comprises 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 for filling holes by interpolation of an adaptive partial angle stacked data volume as claimed in any one of claims 1 to 7 is implemented.

Images