Self-adaptive partial angle superposition data volume interpolation void filling method and device

Abstract

The invention relates to the field of seismic signal processing, and discloses a self-adaptive partial angle superposition data volume interpolation void filling method and a self-adaptive partial angle superposition data volume interpolation void filling device. The method comprises the following steps: determining a seismic trace with a hole, deleting data of the seismic trace with the hole from original seismic data, and obtaining an input sample space; determining interpolation data according to the sample space and the original seismic data; determining an absolute amplitude value of the original seismic data, determining an area of which the absolute amplitude value is equal to a second preset threshold value as a cavity area, and generating a cavity factor for identifying whether the area is the cavity area; according to the void factor and the interpolation data, determining interpolation void filling data of the void area; determining complete interpolation void filling data according to the original seismic data and the interpolation void filling data; and de-noising processing is carried out on the complete interpolation void filling data to obtain final adaptive interpolation void filling data. And the conditions of boundary effect and data incompleteness caused by holes and unsatisfactory inversion effect caused by the holes are overcome.

Description

Technical Field

[0001] The present application relates to the technical field of seismic signal processing, and in particular to an adaptive partial angle stacking data volume interpolation hole filling method, device, storage medium and electronic device. Background Art

[0002] The purpose of the background technology description provided here is to give an overall background of the present application. The statements in this section merely provide background related to the present application and do not necessarily constitute prior art.

[0003] Prestack inversion combines the advantages of AVO (Amplitude variation with offset, AVO for short) inversion and poststack impedance inversion. The elastic impedance inversion method using multiple partial angle stacking data volumes has become a hot topic of research for geophysicists.

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

[0005] First, multiple partial angle stacking data volumes still contain AVO information, which has more information than conventional stacking data and has the potential to be used to obtain subsurface elastic parameters through inversion methods;

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

[0007] 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, when the target layer is deep, the voids in the data will seriously affect the inversion results of the seismic data. The existing pre-stack inversion interpretation processing does not consider the impact of voids in the data, especially when there are no voids in the target layer data, the void problem is basically ignored. However, in actual inversion, even if there are no voids in the target layer, the boundary effect will affect the application effect of the inversion method and the integrity of the inversion results.

[0008] Therefore, there is an urgent need for an adaptive partial angle stacking data volume interpolation method to fill holes in order to solve the above-mentioned defects. Summary of the invention

[0009] In response to the above problems, the present application proposes an adaptive partial angle stacking data volume interpolation hole filling method, device, storage medium and electronic device. The adaptive interpolation hole filling processing method based on the data itself proposed in the present application can adaptively monitor the holes in the data, and combine the matching pursuit Fourier interpolation method to fill the holes, 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 and increase the stability and integrity of the inversion results.

[0010] In a first aspect of the present application, an adaptive method for filling holes by interpolation of partial angle stacked data volumes is provided, the method comprising:

[0011] 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;

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

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

[0014] Determine interpolation hole filling data for the hole area according to the hole factor and the interpolation data;

[0015] Determining complete interpolation gap filling data according to the original seismic data and the interpolation gap filling data;

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

[0017] Furthermore, the original seismic data includes:

[0018] Seismic data with partial angle stacking.

[0019] Further, the determining of the seismic traces with voids includes:

[0020] Determine the number of non-zero sample point values ​​of the seismic trace within a preset time length;

[0021] 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.

[0022] Further, determining the interpolation data according to the sample space and the original seismic data includes:

[0023] Using the position information of the original seismic data as the desired output space of the interpolation process;

[0024] Matching pursuit Fourier interpolation processing is performed according to the sample space and the expected output space to obtain interpolation data.

[0025] Further, the determining the interpolation hole filling data of the hole area according to the hole factor and the interpolation data includes:

[0026] The product of the hole factor and the interpolation data is used as the interpolation hole filling data of the hole area.

[0027] Further, the determining of complete interpolation gap filling data based on the original seismic data and the interpolation gap filling data includes:

[0028] The sum of the original seismic data and the interpolated void-filling data is used as the complete interpolated void-filling data.

[0029] Furthermore, the denoising process of the complete interpolation hole filling data includes:

[0030] The complete interpolation hole filling data is subjected to matching pursuit Fourier interpolation denoising processing.

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

[0032] 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;

[0033] An interpolation data determination module, used to determine interpolation data according to the sample space and the original seismic data;

[0034] a hole region determination module, configured to determine the 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 hole region, and generate a hole factor for identifying whether the region is a hole region;

[0035] An interpolation hole filling data determination module, used to determine the interpolation hole filling data of the hole area according to the hole factor and the interpolation data;

[0036] A complete interpolation gap filling data determination module, used to determine complete interpolation gap filling data according to the original seismic data and the interpolation gap filling data;

[0037] 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.

[0038] According to a third aspect of the present application, a computer-readable storage medium is provided, wherein a computer program stored in the computer-readable storage medium can be executed by one or more processors to implement the steps of the method described above.

[0039] The fourth aspect of the present application provides an electronic device, including a memory and one or more processors, wherein the memory stores a computer program, the memory and the one or more processors are communicatively connected to each other, and when the computer program is executed by the one or more processors, the steps of the method described above are implemented.

[0040] Compared with the prior art, the advantages or beneficial effects of the technical solution of the present application include:

[0041] According to the method for processing holes by adaptive interpolation based on the data itself disclosed in this 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 and increase the stability and integrity of the inversion results. This 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 the seismic inversion calculation and the interpretive processing of the inversion basic data. BRIEF DESCRIPTION OF THE DRAWINGS

[0042] In order to more clearly illustrate the embodiments of the present application or the technical solutions in the prior art, the drawings required for use in the embodiments or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are merely embodiments of the present application. For ordinary technicians in the relevant field, other drawings can be obtained based on the provided drawings without paying any creative work.

[0043] It should also be noted that, for the convenience of description, only the parts related to the present disclosure are shown in the drawings. The drawings constituting a part of the specification of this application are used to provide a further understanding of the present application. The illustrative embodiments and their descriptions in this 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 stacking data volume interpolation method for filling holes provided in an embodiment of the present application;

[0045] Figure 2 A schematic diagram of the effect of an adaptive partial angle stacking data volume interpolation hole filling method provided in an embodiment of the present application;

[0046] Figure 3A schematic diagram of the explanatory processing effects before and after the interpolation and hole filling of an adaptive partial angle stacking data volume provided in an embodiment of the present application;

[0047] Figure 4 A flowchart of another adaptive partial angle stacking data volume interpolation hole filling method provided in an embodiment of the present application;

[0048] Figure 5 A schematic diagram of the structure of an adaptive partial angle stacking data volume interpolation hole filling device provided in an embodiment of the present application;

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

[0050] The following will describe the implementation methods of the present application in detail with reference to the accompanying drawings and embodiments, so that the implementation process of how the present application applies technical means to solve technical problems and achieve corresponding technical effects can be fully understood and implemented accordingly. The embodiments of the present application and the various features in the embodiments can be combined with each other without conflict, and the technical solutions formed 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, rather than all of the embodiments. Based on the embodiments in the present application, all other embodiments obtained by technicians in the relevant field without making creative work are within the scope of protection of the present application.

[0052] From the background technology, we can know that pre-stack inversion combines the advantages of AVO (Amplitude variation with offset, referred to as AVO) inversion and post-stack impedance inversion. The elastic wave impedance inversion method using multiple partial angle stacking data volumes has become a hot topic for geophysicists.

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

[0054] First, multiple partial angle stacking data volumes still contain AVO information, which has more information than conventional stacking data and has the potential to be used to obtain subsurface elastic parameters through inversion methods;

[0055] Secondly, since multiple partial angle stacking data volumes are the result of superimposing partial data within a certain angle range, they have a higher signal-to-noise ratio than conventional pre-stack gather data, which helps 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, when the target layer is deep, the voids in the data will seriously affect the inversion results of the seismic data. The existing pre-stack inversion interpretation processing does not consider the impact of voids in the data, especially when there are no voids in the target layer data, the void problem is basically ignored. However, in actual inversion, even if there are no voids in the target layer, the boundary effect will affect the application effect of the inversion method and the integrity of the inversion results.

[0057] In view of this, the present application proposes an adaptive method for interpolating and filling holes in partial angle stacked data volumes. 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 the seismic inversion calculation and the interpretive processing of the inversion basic data.

[0058] Embodiment 1

[0059] This embodiment provides an adaptive partial angle superposition data volume interpolation method for filling holes. Figure 1 A flowchart of an adaptive partial angle stacking data volume interpolation method for filling holes provided in an embodiment of the present application is shown in FIG. Figure 1 As shown, the method disclosed in this embodiment includes the following steps:

[0060] Step 110, determining the seismic traces with voids, and deleting 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 the absolute amplitude value of the original seismic data, determining the area where the absolute amplitude value is equal to the second preset threshold as a hole area, and generating a hole factor for identifying whether it is a hole area;

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

[0064] Step 150: determining complete interpolation gap filling data according to the original seismic data and the interpolation gap filling data;

[0065] Step 160: De-noise the complete interpolation hole filling data to obtain final adaptive interpolation hole filling data.

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

[0067] Seismic data with partial angle stacking.

[0068] In this embodiment, the purpose is to fill the holes in some angle stacked data volumes, so as to provide regular and complete basic data with hard boundaries of holes eliminated for obtaining good inversion results and the application of subsequent inversion and other interpretive processing methods.

[0069] In some embodiments, the step of determining a seismic trace in which a cavity exists comprises:

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

[0071] 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.

[0072] As an example, the number of non-zero sample values ​​of a 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 a seismic trace with holes, that is, a 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, determining interpolation data according to the sample space and the original seismic data includes:

[0075] Using the position information of the original seismic data as the desired output space of the interpolation process;

[0076] Matching pursuit Fourier interpolation processing is performed according to the sample space and the expected output space to obtain interpolation data.

[0077] As an example, the position information of the original data is used as the expected output space of the interpolation processing, and 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 processing.

[0078] In some embodiments, determining the interpolation hole filling data of the hole area according to the hole factor and the interpolation data includes:

[0079] The product of the hole factor and the interpolation data is used as the interpolation hole filling data of the hole area.

[0080] As an example, the absolute amplitude value of the original seismic data is first determined, and then the area where the absolute amplitude value is equal to the second preset threshold is regarded as the hole area, and the area where the absolute amplitude value is not equal to the second preset threshold is regarded as the non-hole area. The hole factor of the hole area is represented by 1, and the hole factor of the non-hole area is represented by 0.

[0081] Furthermore, the product obtained by multiplying the hole factor and the interpolation data is used as the interpolation hole filling data of the hole area.

[0082] Optionally, the second preset threshold may be set to 0.

[0083] In some embodiments, determining complete interpolation gap filling data according to the original seismic data and the interpolation gap filling data includes:

[0084] The sum of the original seismic data and the interpolated void-filling data is used as the complete interpolated void-filling data.

[0085] As an example, the original data is added to the interpolated hole-filling data of the hole area to obtain the complete interpolated hole-filling data.

[0086] In some embodiments, the denoising process of the complete interpolation hole filling data includes:

[0087] The complete interpolation hole filling data is subjected to matching pursuit Fourier interpolation denoising processing.

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

[0089] The adaptive partial angle stacking data volume interpolation hole filling method provided in this embodiment can adaptively monitor the 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, so as to reduce the boundary effect of the data during seismic inversion and increase 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 according to 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 the second preset threshold as a hole area, and generating a hole factor for identifying whether it is a hole area; Step 140, determining the interpolation hole-filling data of the hole area 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, denoising the complete interpolation hole-filling data to obtain final adaptive interpolation hole-filling data. According to the method for processing holes by adaptive interpolation 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. This 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 the seismic inversion calculation and the interpretive processing of the inversion basic data.

[0090] Embodiment 2

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

[0092] As an example, see Figure 1 , Figure 1 A flowchart of an adaptive partial angle stacking data volume interpolation method for filling holes provided in an embodiment of the present application is shown in FIG. 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 this embodiment, the purpose is to fill the holes in some angle stacked data volumes, so as to provide regular and complete basic data with hard boundaries of holes eliminated for obtaining good inversion results and the application of subsequent inversion and other interpretive processing methods.

[0097] In some embodiments, the step of determining a seismic trace in which a cavity exists comprises:

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

[0099] 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.

[0100] As an example, the number of non-zero sample values ​​of a 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 a seismic trace with holes, that is, a trace that needs to be interpolated.

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

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

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

[0104] Using the position information of the original seismic data as the desired output space of the interpolation process;

[0105] Matching pursuit Fourier interpolation processing is performed according to the sample space and the expected output space to obtain interpolation data.

[0106] As an example, the position information of the original data is used as the expected output space of the interpolation processing, and 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 processing.

[0107] Step 130: determine the absolute amplitude value of the original seismic data, determine the area where the absolute amplitude value is equal to the second preset threshold as a hole area, and generate a hole factor for identifying whether it is a hole area.

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

[0109] In some embodiments, determining the interpolation hole filling data of the hole area according to the hole factor and the interpolation data includes:

[0110] The product of the hole factor and the interpolation data is used as the interpolation hole filling data of the hole area.

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

[0112] Furthermore, the product obtained by multiplying the hole factor and the interpolation data is used as the interpolation hole filling data of the hole area.

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

[0114] Step 150: Determine complete interpolation gap filling data based on the original seismic data and the interpolation gap filling data.

[0115] In some embodiments, determining complete interpolation gap filling data according to the original seismic data and the interpolation gap filling data includes:

[0116] The sum of the original seismic data and the interpolated void-filling data is used as the complete interpolated void-filling data.

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

[0118] Step 160: De-noise the complete interpolation hole filling data to obtain final adaptive interpolation hole filling data.

[0119] In some embodiments, the denoising process of the complete interpolation hole filling data includes:

[0120] The complete interpolation hole filling data is subjected to matching pursuit Fourier interpolation denoising processing.

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

[0122] Furthermore, the disclosed adaptive partial angle stacking data volume interpolation method for filling holes is applied to the partial angle stacking data volume interpolation processing of a certain work area, and the effect can be referred to Figure 2 and Figure 3 .in, Figure 2Part a in the figure is the data volume before the hole is filled by interpolation of the adaptive partial angle superposition data volume, and there is an obvious hard boundary at the hole; Figure 2 Part b in the figure is the data volume after interpolation and hole filling. It can be seen that the interpolated data at the hole has a good correlation with the surrounding data, and there is no obvious difference in signal-to-noise ratio and energy. Figure 2 Part c in the figure is the difference data before and after the interpolation to fill the hole. From the difference data, it can be seen that the present invention has little effect on other valid signal data that are not holes, and there is basically no residual before and after the processing. Figure 3 Part a of Figure 3 Part b in the figure is the data before and after the adaptive partial angle superposition data volume filling. Figure 3 Part c of Figure 3 Part d in the figure shows the data before and after the explanatory processing of the hole filling. It can be seen that the hole filling process can effectively reduce the error problems caused by boundary effects when other subsequent methods are applied.

[0123] The adaptive partial angle stacking data volume interpolation hole filling method provided in this embodiment can adaptively monitor the 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, so as to reduce the boundary effect of the data during seismic inversion and increase 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 according to 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 the second preset threshold as a hole area, and generating a hole factor for identifying whether it is a hole area; Step 140, determining the interpolation hole-filling data of the hole area 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, denoising the complete interpolation hole-filling data to obtain final adaptive interpolation hole-filling data. According to the method for processing holes by adaptive interpolation 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. This 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 the seismic inversion calculation and the interpretive processing of the inversion basic data.

[0124] Embodiment 3

[0125] This embodiment provides an adaptive partial angle stacking data volume interpolation hole filling device. This device embodiment can be used to execute the method embodiment of this application. For details not disclosed in this device embodiment, please refer to the method embodiment of this application. Figure 5 A schematic diagram of a structure of a device provided in an embodiment of the present application is shown in FIG. Figure 5 As shown, the device 500 disclosed in this embodiment includes:

[0126] The sample space determination module 501 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 the input sample space for interpolation processing;

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

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

[0129] An interpolation hole filling data determination module 504 is used to determine the interpolation hole filling data of the hole area according to the hole factor and the interpolation data;

[0130] A complete interpolation gap filling data determination module 505 is used to determine complete interpolation gap filling data according to the original seismic data and the interpolation gap filling data;

[0131] The denoising processing module 506 is used to perform denoising processing on the complete interpolation hole filling data to obtain final adaptive interpolation hole filling data.

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

[0133] Seismic data with partial angle stacking.

[0134] In some embodiments, the acquisition module 501 includes: a first determination unit and a second determination unit; wherein,

[0135] A first determining unit, used to determine the number of non-zero sample point values ​​of a seismic trace within a preset time length;

[0136] The second determining unit is configured to determine that the seismic trace is a seismic trace with a void if the number of non-zero sample values ​​is not less than a first preset threshold.

[0137] In some embodiments, the interpolation data determination module 502 includes: a desired output space determination unit, an interpolation processing unit; wherein,

[0138] An expected output space determination unit, used for taking the position information of the original seismic data as the expected output space for interpolation processing;

[0139] The interpolation processing unit is used to perform matching pursuit Fourier interpolation processing according to the sample space and the expected output space to obtain interpolation data.

[0140] In some embodiments, the interpolation hole filling data determination module 504 is used to use the product of the hole factor and the interpolation data as the interpolation hole filling data of the hole area.

[0141] In some embodiments, the complete interpolation void filling data determination module 505 is used to take the sum of the original seismic data and the interpolation void filling data as the complete interpolation void filling data.

[0142] In some embodiments, the denoising module 506 is used to perform matching pursuit Fourier interpolation denoising on the complete interpolation hole filling data.

[0143] It can be understood by those skilled in the art that Figure 5 The structure shown in the figure does not constitute a limitation on the device of the embodiment of the present application, and may include more or fewer modules / units than shown in the figure, or a combination of certain modules / units, or different module / unit arrangements.

[0144] It should be understood by those skilled in the art that the above modules or steps of the present application can be implemented by a general-purpose computing device, and they can be concentrated on a single computing device, or distributed on a network composed of multiple computing devices. Optionally, they can be implemented by a program code executable by a computing device, so that they can be stored in a storage device and executed by the computing device, and in some cases, the steps shown or described can be executed in a different order than herein, or they can be made into individual integrated circuit modules, or multiple modules or steps therein can be made into a single integrated circuit module for implementation.

[0145] Technicians in the relevant field can clearly understand that, for the convenience and simplicity of description, the specific working process of each module of the adaptive partial angle superposition data volume interpolation hole filling device can refer to the corresponding process in the aforementioned method embodiment, and this embodiment will not be repeated here.

[0146] At the same time, in order to facilitate understanding of the technical solution of this application, you can also refer to Figure 4 .

[0147] The device provided in this embodiment includes: a sample space determination module 501, which is used to determine the seismic traces with holes, and delete the data of the seismic traces with holes from the original seismic data to obtain an input sample space for interpolation processing; an interpolation data determination module 502, which is used to determine the interpolation data according to the sample space and the original seismic data; a hole area determination module 503, which is used to determine the absolute amplitude value of the original seismic data, determine the area where the absolute amplitude value is equal to the second preset threshold as the hole area, and generate a hole factor for identifying whether it is a hole area; an interpolation hole filling data determination module 504, which is used to determine the interpolation hole filling data of the hole area according to the hole factor and the interpolation data; a complete interpolation hole filling data determination module 505, which is used to determine the complete interpolation hole filling data according to the original seismic data and the interpolation hole filling data; a denoising processing module 506, which is used to denoise the complete interpolation hole filling data to obtain the final adaptive interpolation hole filling data. It can adaptively monitor the holes in the data and fill the holes with matching pursuit Fourier interpolation method, 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 that the inversion effect is not ideal due to the boundary effect caused by holes in the data and the incomplete data during the seismic inversion calculation and the interpretive processing of the inversion basic data.

[0148] Embodiment 4

[0149] This embodiment provides a computer-readable storage medium. The computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, all or part of the steps of the method in the above method embodiment can be implemented:

[0150] 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;

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

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

[0153] Determine interpolation hole filling data for the hole area according to the hole factor and the interpolation data;

[0154] Determining complete interpolation gap filling data according to the original seismic data and the interpolation gap filling data;

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

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

[0157] Seismic data with partial angle stacking.

[0158] In this embodiment, the purpose is to fill the holes in some angle stacked data volumes, so as to provide regular and complete basic data with hard boundaries of holes eliminated for obtaining good inversion results and the application of subsequent inversion and other interpretive processing methods.

[0159] In some embodiments, the step of determining a seismic trace in which a cavity exists comprises:

[0160] Determine the number of non-zero sample point values ​​of the seismic trace within a preset time length;

[0161] 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.

[0162] As an example, the number of non-zero sample values ​​of a 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 a seismic trace with holes, that is, a trace that needs to be interpolated.

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

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

[0165] Using the position information of the original seismic data as the desired output space of the interpolation process;

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

[0167] As an example, the position information of the original data is used as the expected output space of the interpolation processing, and 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 processing.

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

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

[0170] As an example, the absolute amplitude value of the original seismic data is first determined, and then the area where the absolute amplitude value is equal to the second preset threshold is regarded as the hole area, and the area where the absolute amplitude value is not equal to the second preset threshold is regarded as the non-hole area. The hole factor of the hole area is represented by 1, and the hole factor of the non-hole area is represented by 0.

[0171] Furthermore, the product obtained by multiplying the hole factor and the interpolation data is used as the interpolation hole filling data of the hole area.

[0172] Optionally, the second preset threshold may be set to 0.

[0173] In some embodiments, determining complete interpolation gap filling data according to the original seismic data and the interpolation gap filling data includes:

[0174] The sum of the original seismic data and the interpolated void-filling data is used as the complete interpolated void-filling data.

[0175] As an example, the original data is added to the interpolated hole-filling data of the hole area to obtain the complete interpolated hole-filling data.

[0176] In some embodiments, the denoising process of the complete interpolation hole filling data includes:

[0177] The complete interpolation hole filling data is subjected to matching pursuit Fourier interpolation denoising processing.

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

[0179] Among them, the computer-readable storage medium may also include computer programs, data files, data structures, etc. alone, or include a combination thereof. The computer-readable storage medium or computer program may be specifically designed and understood by a technician in the field of computer software, or the computer-readable storage medium may be known and available to a technician in the field of computer software. Examples of computer-readable storage media include: magnetic media, such as hard disks, floppy disks, and tapes; optical media, such as CDROM disks and DVDs; magneto-optical media, such as optical disks; and hardware devices, specifically configured to store and execute computer programs, such as read-only memory (ROM), random access memory (RAM), flash memory; or servers, app application stores, etc. Examples of computer programs include machine code (e.g., code generated by a compiler) and files containing high-level code, which can be executed by a computer using an interpreter. The described hardware device may be configured to be used as one or more software modules to perform the operations and methods described above, and vice versa. In addition, the computer-readable storage medium may be distributed in a networked computer system, and program code or computer programs may be stored and executed in a decentralized manner.

[0180] Embodiment 5

[0181] This embodiment provides a computer program product. The computer program product includes a computer program or instructions, and when the computer program or instructions are executed by a processor, all or part of the steps of the method in the above method embodiment are implemented:

[0182] 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;

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

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

[0185] Determine interpolation hole filling data for the hole area according to the hole factor and the interpolation data;

[0186] Determining complete interpolation gap filling data according to the original seismic data and the interpolation gap filling data;

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

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

[0189] Seismic data with partial angle stacking.

[0190] In this embodiment, the purpose is to fill the holes in some angle stacked data volumes, so as to provide regular and complete basic data with hard boundaries of holes eliminated for obtaining good inversion results and the application of subsequent inversion and other interpretive processing methods.

[0191] In some embodiments, the step of determining a seismic trace in which a cavity exists comprises:

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

[0193] 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.

[0194] As an example, the number of non-zero sample values ​​of a 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 a seismic trace with holes, that is, a 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, determining interpolation data according to the sample space and the original seismic data includes:

[0197] Using the position information of the original seismic data as the desired output space of the interpolation process;

[0198] Matching pursuit Fourier interpolation processing is performed according to the sample space and the expected output space to obtain interpolation data.

[0199] As an example, the position information of the original data is used as the expected output space of the interpolation processing, and 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 processing.

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

[0201] The product of the hole factor and the interpolation data is used as the interpolation hole filling data of the hole area.

[0202] As an example, the absolute amplitude value of the original seismic data is first determined, and then the area where the absolute amplitude value is equal to the second preset threshold is regarded as the hole area, and the area where the absolute amplitude value is not equal to the second preset threshold is regarded as the non-hole area. The hole factor of the hole area is represented by 1, and the hole factor of the non-hole area is represented by 0.

[0203] Furthermore, the product obtained by multiplying the hole factor and the interpolation data is used as the interpolation hole filling data of the hole area.

[0204] Optionally, the second preset threshold may be set to 0.

[0205] In some embodiments, determining complete interpolation gap filling data according to the original seismic data and the interpolation gap filling data includes:

[0206] The sum of the original seismic data and the interpolated void-filling data is used as the complete interpolated void-filling data.

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

[0208] In some embodiments, the denoising process of the complete interpolation hole filling data includes:

[0209] The complete interpolation hole filling data is subjected to matching pursuit Fourier interpolation denoising processing.

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

[0211] Further, the computer program product may include one or more computer executable components configured to perform the embodiments when the program is run; the computer program product may also include a computer program tangibly contained on a readable medium, the computer program including program code for performing any method in the embodiments of the present disclosure. In such an embodiment, the computer program may be downloaded and installed from a network through a communication portion, and / or installed from a removable medium.

[0212] Embodiment 6

[0213] This embodiment provides an electronic device. Figure 6 A connection block diagram of an electronic device provided in an embodiment of the present application, such as Figure 6 As shown, the electronic device 600 may 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 used to execute all or part of the steps in the above method embodiment. The memory 602 is used to store various types of data, which may include instructions of any application or method in the electronic device, and data related to the application.

[0215] The memory 602 can be implemented by any type of volatile or non-volatile storage device 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, disk or optical disk.

[0216] The one or more processors 601 may 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 components, and are used to execute all or part of the steps in the aforementioned method embodiment:

[0217] 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;

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

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

[0220] Determine interpolation hole filling data for the hole area according to the hole factor and the interpolation data;

[0221] Determining complete interpolation gap filling data according to the original seismic data and the interpolation gap filling data;

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

[0223] The multimedia component 603 may include a screen and an audio component, the screen may be a touch screen, and the audio component is used to output and / or input audio signals. For example, the audio component may include a microphone, and the microphone is used to receive external audio signals. The received audio signal may be further stored in a memory or sent via a 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 one or more processors 601 and other interface modules, which may be a keyboard, a mouse, buttons, etc. These buttons may 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. Wired communication includes communication through a network port, a serial port, etc.; 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 present application provides an adaptive partial angle stacking data volume interpolation hole filling method, device, computer readable storage medium and electronic device. The adaptive partial angle stacking data volume interpolation hole filling method can adaptively monitor the holes in the data, and combine the matching pursuit Fourier interpolation method to fill the holes, 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 and increase 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 according to 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 the second preset threshold as a hole area, and generating a hole factor for identifying whether it is a hole area; Step 140, determining the interpolation hole-filling data of the hole area 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, denoising the complete interpolation hole-filling data to obtain final adaptive interpolation hole-filling data. According to the method for processing holes by adaptive interpolation 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. This 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 the seismic inversion calculation and the interpretive processing of the inversion basic data.

[0227] In addition, it should be understood that the method or device disclosed in the embodiments provided in the present application can also be implemented in other ways. The method or device embodiments described above are merely schematic. For example, the flowcharts and block diagrams in the accompanying drawings show the possible architecture, functions and operations of the methods and devices according to the multiple embodiments of the present application. In this regard, each box in the flowchart or block diagram can represent a module, a computer program segment or a part of a computer program, and the module, computer program segment or a part of a computer program contains one or more computer programs for implementing the specified logical functions. It should also be noted that in some alternative implementations, the functions marked in the box can also occur in a different order from the order marked in the accompanying drawings, and can actually be executed in parallel. They can also be executed in the opposite order sometimes, depending on the functions involved. It should also be noted that each box in the block diagram and / or flowchart, and the combination of boxes in the block diagram and / or flowchart can be implemented with a dedicated hardware-based system that performs the specified function or action, or can be implemented with a combination of dedicated hardware and computer programs.

[0228] In the present application, the terms "comprises", "includes" or any other variants thereof are intended to cover non-exclusive inclusion, so that a process, method, article or device including a series of elements includes not only those elements, but also other elements not explicitly listed, or also includes elements inherent to such process, method, article or device. In the absence of further restrictions, the elements defined by the sentence "including one..." do not exclude the existence of other identical elements in the process, method, device or device including the elements; if there is a description of "first", "second", etc., it is only for descriptive purposes and cannot be understood as indicating or implying relative importance or implicitly indicating the number of indicated technical features or implicitly indicating the order of the indicated technical features; in the description of the present application, unless otherwise specified, the terms "multiple" and "multiple" mean at least two.

[0229] It should be noted that if a server is described, the server can be an independent physical server or terminal, or a server cluster composed of multiple physical servers, or a cloud server that can provide basic cloud computing services such as cloud servers, cloud databases, cloud storage and CDN; if a smart terminal or mobile device is described in this application, it should be noted that the smart terminal or mobile device can be a mobile phone, tablet computer, smart watch, netbook, wearable electronic device, personal digital assistant (PDA), augmented reality technology equipment (AR), virtual reality equipment (VR), smart TV, smart speaker, personal computer (PC), etc., but is not limited to this. This application does not specifically limit the specific form of the smart terminal or mobile device.

[0230] Finally, it should be noted that, in the description of this specification, the description with reference to the terms "one embodiment", "some embodiments", "example", "an example" or "some examples" means that the specific features, structures, materials or characteristics described in conjunction with the embodiment or example are included in at least one embodiment or example of the present application. In this specification, the schematic representation of the above terms does not necessarily refer to the same embodiment or example. Moreover, the specific features, structures, materials or characteristics described may 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 is understood that the above embodiments are exemplary, and the contents described are only embodiments adopted for the convenience of understanding the present application, and are not intended to limit the present application. Any technician in the technical field to which the present application belongs can make any modifications and changes in the form and details of the implementation without departing from the spirit and scope disclosed in the present application, but the scope of protection of the present application shall still be based on the scope defined in the attached claims.

Claims

1. An adaptive partial angle superposition data volume interpolation method for filling holes, 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; Determine interpolation data according to the sample space and the original seismic data; Determine the absolute amplitude value of the original seismic data, determine the area where the absolute amplitude value is equal to the second preset threshold as a hole area, and generate a hole factor for identifying whether it is a hole area; Determine interpolation hole filling data for the hole area according to the hole factor and the interpolation data; Determining complete interpolation gap filling data according to the original seismic data and the interpolation gap filling data; The complete interpolation hole filling data is denoised to obtain final adaptive interpolation hole filling data.

2. The method for filling holes by interpolation of partially angled superimposed data volumes according to claim 1, characterized in that: The original seismic data includes: Seismic data with partial angle stacking.

3. The method for filling holes by interpolation of partially angled superimposed data volumes according to claim 1, characterized in that: The method of determining the seismic traces where cavities exist includes: Determine the number of non-zero sample point 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 method for filling holes by interpolation of partially angled stacked data volumes according to claim 1, characterized in that: The step of determining interpolation data according to the sample space and the original seismic data comprises: 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 stacked data volume according to claim 1, characterized in that: The step of determining the interpolation hole filling data of the hole area according to the hole factor and the interpolation data comprises: The product of the hole factor and the interpolation data is used as the interpolation hole filling data of the hole area.

6. The method for filling holes by interpolation of partially angled stacked data volumes according to claim 1, characterized in that: The step of determining complete interpolation gap filling data according to the original seismic data and the interpolation gap filling data comprises: The sum of the original seismic data and the interpolated void-filling data is used as the complete interpolated void-filling data.

7. The method for filling holes by interpolation of partially angled stacked data volumes according to claim 1, characterized in that: The denoising process of the complete interpolation hole filling data comprises: 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, used to determine interpolation data according to the sample space and the original seismic data; a hole region determination module, configured to determine the 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 hole region, and generate a hole factor for identifying whether the region is a hole region; An interpolation hole filling data determination module, used to determine the interpolation hole filling data of the hole area according to the hole factor and the interpolation data; A complete interpolation gap filling data determination module, used to determine complete interpolation gap filling data according to the original seismic data and the interpolation gap 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 as described in 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 described in any one of claims 1 to 7 is implemented.

Images