Method, device and equipment for effectively suppressing near-offset residual multiple of gather

Abstract

The application discloses a method, device and equipment for effectively suppressing near-offset residual multiple of a gather, and the method comprises the following steps: performing prestack migration processing on preprocessed seismic data to obtain CRP gathers; sorting the CRP gathers in ascending order of offset distance; randomly sorting the data contained in each CRP gather; attenuating the multiple in the randomly sorted data; obtaining filtered CRP gather data based on the attenuated data for each CRP gather; and sorting the filtered CRP gather data in ascending order of offset distance. The method for effectively suppressing near-offset residual multiple of a gather can suppress near-offset multiple under the premise of fidelity, restore the true AVO effect, improve imaging accuracy, improve migration velocity spectrum, improve velocity picking accuracy, and improve the imaging quality of migration stacking records, thereby providing guarantee for high-precision seismic imaging and fine interpretation.

Description

TECHNICAL FIELD

[0001] The present application belongs to the technical field of geophysical exploration seismic data processing, and particularly relates to a method and device for effectively suppressing near-offset residual multiple waves of a gather, a medium and an electronic equipment. BACKGROUND

[0002] Multiple wave suppression has always been one of the focuses and difficulties of seismic data processing work. There are many processing methods for multiple wave suppression. The main evaluation object of the current multiple wave suppression method is the migrated stack record, and little attention is paid to the quality of the migrated gather. However, these multiple wave suppression methods all have a certain residual multiple wave, which is particularly obvious in the near offset.

[0003] With the continuous deepening of oil and gas exploration and development, the types of oil and gas reservoirs gradually change from structure to lithology, so the prestack reservoir prediction work is particularly important. A prerequisite for prestack reservoir prediction work is to have a gather with relatively amplitude preservation. However, the existence of near-offset residual multiple waves destroys the amplitude preservation and needs to be suppressed.

[0004] Currently, there are not many methods and technologies in the industry that are specifically aimed at near-offset residual multiple wave suppression, and the effect is not ideal. Therefore, effective suppression of near-offset residual multiple waves is a major difficulty in current research.

[0005] Therefore, there is an urgent need for a method for effectively suppressing near-offset residual multiple waves. SUMMARY

[0006] The purpose of the present application is to provide a method for effectively suppressing near-offset residual multiple waves.

[0007] In a first aspect, the present application provides a method for effectively suppressing near-offset residual multiple waves of a gather, comprising: performing prestack migration processing on preprocessed seismic data to obtain a CRP gather; sorting the CRP gather in order of increasing offset; randomly sorting the data contained in each CRP gather; attenuating the multiple waves in the randomly sorted data; for each CRP gather, obtaining filtered CRP gather data based on the attenuated data; and sorting the filtered CRP gather data in order of increasing offset.

[0008] Optionally, the obtaining of the filtered CRP gather data based on the attenuated data comprises: cutting off the data within the near offset, and performing stack processing on the data after the near-offset cut-off to obtain the stack record of the CRP gather; and performing matched filtering processing on the stack record of the CRP gather and the near-offset data corresponding to the CRP gather to obtain the filtered CRP gather data.

[0009] Optionally, the least square adaptive matched filter method is used to perform matched filter processing on the CRP trace stack record and the near offset data corresponding to the CRP trace.

[0010] Optionally, the time-frequency domain multi-channel statistics plus interpolation method is used to attenuate the multiple waves in the randomly sorted data.

[0011] In a second aspect, the present application further provides an electronic device, comprising: a memory storing executable instructions; and a processor running the executable instructions in the memory to implement the method for effectively suppressing near offset residual multiple waves of a gather.

[0012] In a third aspect, the present application further provides a computer readable storage medium storing a computer program, which is executed by a processor to implement the method for effectively suppressing near offset residual multiple waves of a gather.

[0013] In a fourth aspect, the present application further provides a device for effectively suppressing near offset residual multiple waves of a gather, comprising: a CRP gather acquisition module, which performs pre-stack migration processing on preprocessed seismic data to obtain a CRP gather; a first gather sorting module, which sorts the CRP gather in ascending order of offset distance; a channel data sorting module, which randomly sorts the data contained in each CRP trace; an attenuation module, which attenuates the multiple waves in the randomly sorted data; a filter module, which obtains filtered CRP trace data based on the attenuated data for each CRP trace; and a second gather sorting module, which sorts the filtered CRP trace data in ascending order of offset distance.

[0014] Optionally, the obtaining of the filtered CRP trace data based on the attenuated data comprises: cutting off the data within the near offset distance, and performing stack processing on the data after the near offset distance cutting to obtain a stack record of the CRP trace; and performing matched filter processing on the stack record of the CRP trace and the near offset data corresponding to the CRP trace to obtain the filtered CRP trace data.

[0015] Optionally, the least square adaptive matched filter method is used to perform matched filter processing on the CRP trace stack record and the near offset data corresponding to the CRP trace.

[0016] Optionally, the time-frequency domain multi-channel statistics plus interpolation method is used to attenuate the multiple waves in the randomly sorted data.

[0017] The method for effectively suppressing near-offset residual multiple of a gather according to the application can suppress near-offset multiple well, restore true AVO effect, improve imaging precision, improve migration velocity spectrum, improve velocity picking precision, improve imaging quality of migration stack record, and provide guarantee for high-precision seismic imaging and fine interpretation under the premise of fidelity.

[0018] The application has other characteristics and advantages, which will be apparent or will be described in detail in the drawings and subsequent specific embodiments incorporated herein, which are collectively used to explain the specific principles of the application. BRIEF DESCRIPTION OF DRAWINGS

[0019] The above and other objects, features and advantages of the present application will become more apparent from the following detailed description when taken in conjunction with the accompanying drawings, in which like reference characters refer to like parts throughout the figures, and wherein:

[0020] Figure 1 A flow chart of a method for effectively suppressing near-offset residual multiple of a gather according to one embodiment of the application is shown.

[0021] Figure 2a A pre-migration post-CRP gather of a method for effectively suppressing near-offset residual multiple of a gather according to one embodiment of the application is shown.

[0022] Figure 2b A post-offset random sorting CRP gather of a method for effectively suppressing near-offset residual multiple of a gather according to one embodiment of the application is shown.

[0023] Figure 3a A pre-residual near-offset multiple suppression CRP gather is shown.

[0024] Figure 3b A post-suppression CRP gather of a method for effectively suppressing near-offset residual multiple of a gather according to one embodiment of the application is shown.

[0025] Figure 4a A pre-residual near-offset multiple suppression velocity spectrum is shown.

[0026] Figure 4b A post-suppression velocity spectrum of a method for effectively suppressing near-offset residual multiple of a gather according to one embodiment of the application is shown.

[0027] Figure 5aA method for effectively suppressing residual multiple waves at near-offset in a gather according to an embodiment of the present invention is shown, including pre-migration stacking of residual near-offset multiple waves suppression.

[0028] Figure 5b A method for effectively suppressing residual multiple waves at close offsets in a gather according to an embodiment of the present invention is shown, which includes post-migration stacking.

[0029] Figure 6 The present invention shows a structural block diagram of a device for effectively suppressing residual multiple waves at close offsets in a gather according to an embodiment of the present invention.

[0030] Description of Reference Numerals

[0031] 102, CRP gather acquisition module; 104, first gather sorting module; 106, intra-channel data sorting module; 108, attenuation module; 110, filtering module; 112, second gather sorting module. DETAILED DESCRIPTION

[0032] The preferred embodiments of the present invention will be described in more detail below. Although the preferred embodiments of the present invention are described below, it should be understood that the present invention can be implemented in various forms and should not be limited to the embodiments set forth herein.

[0033] The present invention provides a method for effectively suppressing near-offset residual multiple waves of a gather, comprising: performing pre-stack migration processing on pre-processed seismic data to obtain a CRP gather; sorting the CRP gather in ascending order of offset; randomly sorting data contained in each CRP trace; attenuating multiple waves in the randomly sorted data; obtaining filtered CRP trace data for each CRP trace based on the attenuated data; and sorting each filtered CRP trace data in ascending order of offset.

[0034] Specifically, the data after the conventional preprocessing is subjected to pre-stack migration processing to obtain a CRP (common reflection point) gather, the CRP gather is indexed by CRP number as the first index word and by offset as the second index word, and the CRP gather is sorted in ascending order of offset. On the CRP gather, the data in each CRP trace is reordered by offset, the coherence of the near-offset residual multiple is broken by random reordering to make it behave as random noise, and the reordered data in each CRP trace is subjected to multiple suppression by time-frequency domain multi-channel statistics and interpolation. On this basis, for each CRP trace, the near-offset data after the multiple suppression is subjected to internal cutting and stacking, and then subjected to matched filtering to obtain the filtered CRP trace data.

[0035] According to the example embodiment, the method for effectively suppressing near-offset residual multiple of a gather can well suppress the near-offset multiple under the premise of fidelity, restore the true AVO effect, improve imaging accuracy, improve migration velocity spectrum, improve velocity picking accuracy, and improve the imaging quality of migration stacking records, thereby providing guarantee for high-precision seismic imaging and fine interpretation.

[0036] As an optional solution, obtaining the filtered CRP trace data based on the decayed data includes: performing near-offset internal cutting on the decayed data, stacking the data after the near-offset internal cutting to obtain stacking records of the CRP traces, and performing matched filtering on the stacking records of the CRP traces and the near-offset data corresponding to the CRP traces to obtain the filtered CRP trace data.

[0037] Specifically, the near-offset internal cutting and stacking are performed on the decayed CRP gather to obtain stacking records, the least square adaptive matched filtering method is used to perform matched filtering on the near-offset data and the stacking records, and the filtered results are sorted in ascending order of offset with CRP number as the first index word and offset as the second index word.

[0038] As an optional solution, the least square adaptive matched filtering method is used to perform matched filtering on the stacking records of the CRP traces and the near-offset data corresponding to the CRP traces.

[0039] Specifically, the least square adaptive matched filtering method is used to perform matched filtering on the stacking records of the CRP traces and the near-offset data corresponding to the CRP traces.

[0040] As an optional solution, the time-frequency domain multi-channel statistics and interpolation method is used to attenuate the multiple in the reordered data.

[0041] Specifically, the multiple waves in the randomly sorted data are attenuated in a time-frequency domain multi-channel statistical interpolation manner.

[0042] In a second aspect, the present application also provides an electronic device, which comprises a memory storing executable instructions, and a processor running the executable instructions in the memory to implement the method for effectively suppressing near-offset residual multiple waves of a gather.

[0043] In a third aspect, the present application also provides a computer readable storage medium storing a computer program, which is executed by a processor to implement the method for effectively suppressing near-offset residual multiple waves of a gather.

[0044] In a fourth aspect, the present application also provides a device for effectively suppressing near-offset residual multiple waves of a gather, which comprises a CRP gather obtaining module for obtaining a CRP gather by pre-stack migration processing on preprocessed seismic data, a first gather sorting module for sorting the CRP gather in ascending order of offset distance, an in-gather data sorting module for randomly sorting the data contained in each CRP gather, an attenuation module for attenuating the multiple waves in the randomly sorted data, a filtering module for obtaining filtered CRP gather data based on the attenuated data for each CRP gather, and a second gather sorting module for sorting the filtered CRP gather data in ascending order of offset distance.

[0045] Embodiment One

[0046] Figure 1 A flow chart of a method for effectively suppressing near-offset residual multiple waves of a gather according to an embodiment of the present application is shown. Figure 2a A post-pre-stack migration CRP gather of a method for effectively suppressing near-offset residual multiple waves of a gather according to an embodiment of the present application is shown. Figure 2b A post-offset distance random sorting CRP gather of a method for effectively suppressing near-offset residual multiple waves of a gather according to an embodiment of the present application is shown. Figure 3a A pre-suppression near-offset residual multiple wave CRP gather is shown. Figure 3b A post-suppression CRP gather of a method for effectively suppressing near-offset residual multiple waves of a gather according to an embodiment of the present application is shown. Figure 4a A pre-suppression near-offset residual multiple wave velocity spectrum is shown. Figure 4b A post-suppression velocity spectrum of a method for effectively suppressing near-offset residual multiple waves of a gather according to an embodiment of the present application is shown. Figure 5aA residual near-offset moveout stack before the residual near-offset multiple suppression of a method for effectively suppressing residual near-offset multiple of a gather is shown. Figure 5b A residual near-offset moveout stack after the residual near-offset multiple suppression of a method for effectively suppressing residual near-offset multiple of a gather is shown.

[0047] In combination Figure 1 , Figure 2a , Figure 2b , Figure 3a , Figure 3b , Figure 4a , Figure 4b , Figure 5a and Figure 5b The method for effectively suppressing residual near-offset multiple of a gather includes:

[0048] Step 1: pre-stack migration processing is performed on the pre-processed seismic data to obtain CRP gathers;

[0049] Step 2: the CRP gathers are sorted in ascending order of offset;

[0050] Step 3: the data contained in each CRP gather is randomly sorted;

[0051] Step 4: the multiple in the randomly sorted data is attenuated;

[0052] Step 5: for each CRP gather, filtered CRP gather data is obtained based on the attenuated data;

[0053] Step 6: each filtered CRP gather data is sorted in ascending order of offset.

[0054] Wherein, the filtered CRP gather data is obtained based on the attenuated data, including: the data within the near-offset is cut off, and the data after the near-offset cut-off is stacked to obtain the stack record of the CRP gather; the stack record of the CRP gather and the near-offset data corresponding to the CRP gather are matched and filtered to obtain the filtered CRP gather data.

[0055] Wherein, the least square adaptive matched filtering method is used to match and filter the stack record of the CRP gather and the near-offset data corresponding to the CRP gather.

[0056] Wherein, the multiple in the randomly sorted data is attenuated by using the time-frequency domain multi-channel statistical interpolation method.

[0057] Embodiment two

[0058] Figure 6A structural block diagram of an apparatus for effectively suppressing near-offset residual multiples of a gather is shown.

[0059] As shown in the figure, the apparatus for effectively suppressing near-offset residual multiples of a gather comprises: Figure 6

[0060] The CRP gather obtaining module 102 performs pre-stack migration processing on the preprocessed seismic data to obtain a CRP gather;

[0061] The first gather sorting module 104 sorts the CRP gather in ascending order of offset distance;

[0062] The in-gather data sorting module 106 randomly sorts the data contained in each CRP gather;

[0063] The attenuation module 108 attenuates the multiples in the randomly sorted data;

[0064] The filtering module 110 obtains filtered CRP gather data based on the attenuated data for each CRP gather;

[0065] The second gather sorting module 112 sorts the CRP gather in ascending order of offset distance and sorts the filtered CRP gather data in ascending order of offset distance.

[0066] The filtered CRP gather data is obtained based on the attenuated data, including: cutting off the data within the near offset distance from the attenuated data, and performing stacking processing on the data cut off within the near offset distance to obtain stacking records of the CRP gather; and performing matched filtering processing on the stacking records of the CRP gather and the near offset distance data corresponding to the CRP gather to obtain the filtered CRP gather data.

[0067] The matched filtering processing on the stacking records of the CRP gather and the near offset distance data corresponding to the CRP gather is performed by using a least square adaptive matched filtering method.

[0068] The multiples in the randomly sorted data are attenuated by using a time-frequency domain multi-channel statistical interpolation method.

[0069] Figure 2a For the CRP gather after pre-stack migration, it can be seen that there is obvious residual multiple at the near offset distance; Figure 2b For the CRP gather after random sorting of offset distance, it can be seen that the multiples become more random.

[0070] Figure 3a And Figure 3b ​For the CRP gather comparison chart before and after residual near offset multiple wave suppression by using the method of the application, it can be obviously seen that the residual near offset multiple wave energy is well suppressed, and the real amplitude relationship of the gather is restored.

[0071] Figure 4a and Figure 4b For the velocity spectrum comparison chart before and after residual near offset multiple wave suppression by using the method of the application, it can be obviously seen that the residual near offset multiple wave energy is well suppressed, the primary wave energy is enhanced, and the velocity focusing and explainability are better.

[0072] Figure 5a and Figure 5b For the migration stack comparison chart before and after residual near offset multiple wave suppression by using the method of the application, it can be seen that the real response of the formation is restored, the continuity is increased, and the structural false image is eliminated after the application of the technology.

[0073] Embodiment three

[0074] The present disclosure provides an electronic device, which includes a memory storing executable instructions, and a processor running the executable instructions in the memory to implement the method for effectively suppressing near offset residual multiple waves of a gather.

[0075] The electronic device according to the embodiments of the present disclosure includes a memory and a processor.

[0076] The memory is configured to store non-transitory computer-readable instructions. Specifically, the memory can include one or more computer program products, which can include various forms of computer-readable storage media, such as volatile memory and / or non-volatile memory. The volatile memory may, for example, include random access memory (RAM), cache, and / or the like. The non-volatile memory may, for example, include read-only memory (ROM), hard disk, flash memory, and / or the like.

[0077] The processor can be a central processing unit (CPU) or other forms of processing units having data processing and / or instruction execution capabilities, and can control other components in the electronic device to perform desired functions. In an embodiment of the present disclosure, the processor is configured to run the computer-readable instructions stored in the memory.

[0078] Those skilled in the art should understand that, in order to solve the technical problem of how to obtain a good user experience effect, the present embodiment can also include well-known structures such as a communication bus, an interface, and the like, which should also be included in the protection scope of the present disclosure.

[0079] Detailed description of the present embodiment can refer to the corresponding description in the foregoing embodiments, which will not be repeated here.

[0080] Embodiment Four

[0081] The present disclosure provides a computer readable storage medium, which stores a computer program, and the computer program is executed by a processor to implement the method for effectively suppressing near-offset residual multiple of a gather.

[0082] The computer readable storage medium according to the embodiments of the present disclosure stores non-transitory computer readable instructions. When the non-transitory computer readable instructions are run by a processor, all or part of the steps of the method of the embodiments of the present disclosure are executed.

[0083] The computer readable storage medium includes, but is not limited to, optical storage media (for example, CD-ROM and DVD), magneto-optical storage media (for example, MO), magnetic storage media (for example, magnetic tape or mobile hard disk), media with built-in rewritable non-volatile memory (for example, memory card), and media with built-in ROM (for example, ROM cartridge).

[0084] The above has described the embodiments of the present disclosure, and the above description is exemplary, not exhaustive, and is not limited to the disclosed embodiments. Many modifications and changes are obvious to those skilled in the art without departing from the scope and spirit of the described embodiments.

Claims

1. A method for effectively suppressing residual multiple waves at close offsets in a gather, characterized in that: include: Perform pre-stack migration on the pre-processed seismic data to obtain CRP gathers; sorting the CRP gathers in ascending order of offset; Randomly sort the data contained in each CRP channel; Attenuate multiple waves in randomly sorted data; For each CRP channel, based on the attenuated data, filtered CRP channel data are obtained; Sort each filtered CRP trace data in ascending order of offset.

2. The method for effectively suppressing near-offset residual multiple waves in a gather according to claim 1, characterized in that: The obtaining of filtered CRP channel data based on the attenuated data includes: Performing near-offset excision on the attenuated data, and performing stacking processing on the data after near-offset excision to obtain the stacked record of the CRP trace; Matched filtering is performed on the superimposed record of the CRP trace and the close-offset data corresponding to the CRP trace to obtain filtered CRP trace data.

3. The method for effectively suppressing near-offset residual multiple waves in a gather according to claim 2, characterized in that: The stacked records of the CRP trace and the close-offset data corresponding to the CRP trace are processed by matching filter using the least squares adaptive matched filter method.

4. The method for effectively suppressing near-offset residual multiple waves in a gather according to claim 1, characterized in that: The multiple waves in the randomly sorted data are attenuated by using the multi-channel statistics and interpolation method in the time-frequency domain.

5. An electronic device, characterized in that: The electronic device comprises: a memory storing executable instructions; A processor is configured to execute the executable instructions in the memory to implement the method for effectively suppressing near-offset residual multiple waves of a gather according to any one of claims 1 to 4.

6. A computer-readable storage medium, characterized in that The computer-readable storage medium stores a computer program, which, when executed by a processor, implements the method for effectively suppressing near-offset residual multiple waves of a gather according to any one of claims 1 to 4.

7. A device for effectively suppressing residual multiple waves at close offsets in a gather, characterized in that: include: CRP gather acquisition module, which performs prestack migration on preprocessed seismic data to obtain CRP gathers; A first gather sorting module sorts the CRP gathers in ascending order of offset; The intra-channel data sorting module randomly sorts the data contained in each CRP channel; Attenuation module, which attenuates multiple waves in the randomly sorted data; A filtering module, for each CRP channel, obtains filtered CRP channel data based on the attenuated data; The second gather sorting module sorts the CRP gathers in the order of offset from small to large, and sorts each filtered CRP trace data in the order of offset from small to large.

8. The method for effectively suppressing near-offset residual multiple waves of a gather according to claim 7, characterized in that: The obtaining of filtered CRP channel data based on the attenuated data includes: Performing near-offset excision on the attenuated data, and performing stacking processing on the data after near-offset excision to obtain the stacked record of the CRP trace; Matched filtering is performed on the superimposed record of the CRP trace and the close-offset data corresponding to the CRP trace to obtain filtered CRP trace data.

9. The device for effectively suppressing near-offset residual multiple waves of a gather according to claim 8, characterized in that: The stacked records of the CRP trace and the close-offset data corresponding to the CRP trace are processed by matching filter using the least squares adaptive matched filter method.

10. The device for effectively suppressing near-offset residual multiple waves of a gather according to claim 7, characterized in that: The multiple waves in the randomly sorted data are attenuated by using the multi-channel statistics and interpolation method in the time-frequency domain.

Images