Method for selecting optimal frequency band of deep reflection seismic data, electronic device and storage medium

By sorting and filtering deep-reflected seismic data, the intermediate frequency of the initial advantageous frequency band is expanded to obtain the optimal frequency band, which solves the problem of inability to effectively suppress high-frequency energy and improve the signal-to-noise ratio in the prior art, and effectively highlighting and observability of deep-reflected information is achieved.

CN114722861BActive Publication Date: 2025-04-04GUANGDONG LABORATORY OF SOUTHERN OCEAN SCIENCE AND ENGINEERING (GUANGZHOU) +1
View PDF 2 Cites 0 Cited by

Patent Information

Application Number
CN202210255001.5
Authority / Receiving Office
CN · China
Patent Type
Patents(China)
Current Assignee / Owner
Filing Date
2022-03-15
Publication Date
2025-04-04
Estimated Expiration
2042-03-15

AI Technical Summary

Technical Problem

In the existing deep reflection seismic data processing, high-frequency energy cannot be effectively suppressed, resulting in the inability to highlight the low-frequency information below the substrate, and the signal-to-noise ratio is low, making it difficult to determine the location of the target layer.

Method used

By reordering the deep-reflected seismic data, converting it into a common detection point channel data, forming a single-channel profile, and setting up 6 filter bands to divide the single-channel profile, selecting the initial advantageous frequency band, expanding the intermediate frequency, obtaining the low-frequency and high-frequency dominant frequency bands, and finally finding the optimal frequency band.

Benefits of technology

It realizes the optimal frequency band for accurately obtaining deep-reflection seismic data, highlights deep-reflection information, improves signal-to-noise ratio, and enhances observability of the layers below the substrate.

✦ Generated by Eureka AI based on patent content.

Smart Images

  • Figure CN114722861B_ABST
    Figure CN114722861B_ABST
Patent Text Reader

Abstract

The present application is about a method, electronic device and storage medium for selecting the optimal frequency band for deep reflection seismic data. The method comprises: performing frequency division filtering on a single-channel profile by setting a filter frequency band to obtain a filtered profile, and selecting an initial dominant frequency band according to the signal-to-noise ratio and contrast of the filtered profile. Extending the intermediate frequency of the initial dominant frequency band to obtain an extended dominant frequency band. Obtaining a low-frequency dominant frequency band according to the lowest extended frequency and the lowest frequency of the initial dominant frequency band, and obtaining a high-frequency dominant frequency band according to the highest extended frequency and the highest frequency of the initial dominant frequency band. Obtaining the optimal frequency band according to the low-frequency dominant frequency band and the high-frequency dominant frequency band. The low-frequency dominant frequency band contains more deep reflection information of the target layer, but the image signal-to-noise ratio is low. The high-frequency dominant frequency band contains less deep reflection information of the target layer, but the image signal-to-noise ratio is high. Combining the low-frequency dominant frequency band and the high-frequency dominant frequency band, the reliability of deep reflection seismic data and the observability of the image are taken into account.
Need to check novelty before this filing date? Find Prior Art

Description

Technical Field

[0001] The present application relates to the field of deep reflection seismic data processing, and in particular to a method for selecting an optimal frequency band for deep reflection seismic data, an electronic device and a storage medium. Background Art

[0002] Imaging the crust and mantle structure of the lithosphere in ocean basins is an important basis for studying regional tectonic activities. At present, there is a common problem of difficulty in imaging the Moho surface using offshore seismic data. According to the principle of seismic wave propagation, when seismic waves propagate from soft sedimentary layers to hard igneous rock basements, due to the huge differences in the speed and density properties of the two media, strong reflections will occur, and only a small amount of low-frequency energy can pass through the interface and propagate deep into the crust.

[0003] The current deep reflection seismic data processing still uses full-band processing, which does not fully suppress high-frequency energy and cannot highlight the low-frequency information below the basement. Deep reflection seismic data processing usually cannot determine the location of the target layer, and there is a lot of noise below the basement, and the signal-to-noise ratio is much lower than conventional seismic data.

[0004] Therefore, there is an urgent need for a method for selecting the optimal frequency band of deep reflection seismic data. Summary of the invention

[0005] In order to overcome the problems existing in the related art, the present application provides a method for selecting the optimal frequency band of deep reflection seismic data, an electronic device and a storage medium. The method can accurately obtain the optimal frequency band of deep reflection seismic data, thereby highlighting the deep reflection information.

[0006] The first aspect of the present application provides a method for selecting an optimal frequency band for deep reflection seismic data, comprising:

[0007] Rearrange the deep reflection seismic data, convert the common shot point gather data into common receiver point gather data, and form a single-channel profile;

[0008] Performing frequency division filtering on the single-channel profile by setting a filter frequency band to obtain a filtered profile, and selecting an initial dominant frequency band according to a signal-to-noise ratio and a contrast ratio of the filtered profile;

[0009] The intermediate frequency of the initial dominant frequency band is expanded to obtain an extended dominant frequency band, wherein the range of the extended dominant frequency band is from the lowest extended frequency to the highest extended frequency; a low-frequency dominant frequency band is obtained according to the lowest extended frequency and the lowest frequency of the initial dominant frequency band, and a high-frequency dominant frequency band is obtained according to the highest extended frequency and the highest frequency of the initial dominant frequency band; an optimal frequency band is obtained according to the low-frequency dominant frequency band and the high-frequency dominant frequency band.

[0010] The single-channel profile is subjected to frequency division filtering by setting a filter frequency band, comprising:

[0011] Setting 6 filter frequency bands, wherein the filter frequency bands are respectively set to 1-6 Hz, 6-12 Hz, 12-24 Hz, 24-36 Hz, 36-60 Hz and above 60 Hz;

[0012] The single-channel profile is subjected to frequency division filtering using the six filter frequency bands that are set respectively.

[0013] The step of extending the intermediate frequency of the initial dominant frequency band comprises:

[0014] Calculating the middle frequency of the initial dominant frequency band, and then extending the middle frequency towards the low frequency and high frequency directions respectively;

[0015] The intermediate frequency is an average of the lowest frequency of the initial dominant frequency band and the highest frequency of the initial dominant frequency band.

[0016] The step of extending the intermediate frequency towards the low frequency and the high frequency respectively comprises:

[0017] Subtracting the correction frequency from the intermediate frequency to obtain the lowest extension frequency;

[0018] The intermediate frequency is increased by the correction frequency to obtain a maximum extended frequency;

[0019] Using the lowest extended frequency as the lowest frequency of the extended dominant frequency band, and using the highest extended frequency as the highest frequency of the extended dominant frequency band;

[0020] The correction frequency is half of the intermediate frequency.

[0021] The step of obtaining a low-frequency dominant frequency band according to the lowest extended frequency and the lowest frequency of the initial dominant frequency band comprises:

[0022] The lowest extended frequency is used as the highest frequency of the low-frequency dominant frequency band, and the lowest frequency of the initial dominant frequency band is used as the lowest frequency of the low-frequency dominant frequency band.

[0023] The step of obtaining a high frequency dominant frequency band according to the highest extended frequency and the highest frequency of the initial dominant frequency band comprises:

[0024] The highest extended frequency is used as the lowest frequency of the high-frequency dominant frequency band, and the highest frequency of the initial dominant frequency band is used as the highest frequency of the high-frequency dominant frequency band.

[0025] The obtaining of the optimal frequency band according to the low-frequency dominant frequency band and the high-frequency dominant frequency band comprises:

[0026] Using the middle frequency of the low-frequency dominant frequency band as the lowest frequency of the optimal frequency band, and using the middle frequency of the high-frequency dominant frequency band as the highest frequency of the optimal frequency band;

[0027] The middle frequency of the low frequency dominant frequency band is the average of the lowest frequency of the low frequency dominant frequency band and the highest frequency of the low frequency dominant frequency band;

[0028] The middle frequency of the high frequency dominant frequency band is an average value of the lowest frequency of the high frequency dominant frequency band and the highest frequency of the high frequency dominant frequency band.

[0029] The selecting an initial dominant frequency band according to the signal-to-noise ratio and contrast of the filtering profile comprises:

[0030] Performing frequency division filtering on the single-channel profile to obtain a filtered profile;

[0031] An initial dominant frequency band is selected, where the initial dominant frequency band is the filter frequency band corresponding to the filter profile having the highest signal-to-noise ratio and contrast.

[0032] A second aspect of the present application provides an electronic device, including:

[0033] Processor; and

[0034] The memory stores executable codes thereon, and when the executable codes are executed by the processor, the processor is caused to execute the method as described above.

[0035] A third aspect of the present application provides a non-transitory machine-readable storage medium having executable code stored thereon. When the executable code is executed by a processor of an electronic device, the processor is caused to execute the method as described above.

[0036] The technical solution provided by this application may have the following beneficial effects:

[0037] The technical solution of the present application performs frequency division filtering on a single-channel profile by setting a filter band, selects an initial dominant band according to the signal-to-noise ratio and contrast of the filtered profile, and narrows the range of the filter band. Then, the middle frequency of the initial dominant band is extended to obtain an extended dominant band. A low-frequency dominant band is obtained according to the lowest extended frequency and the lowest frequency of the initial dominant band. The low-frequency dominant band contains more deep reflection information of the target layer, but the image signal-to-noise ratio is low. A high-frequency dominant band is obtained according to the highest extended frequency and the highest frequency of the initial dominant band. The high-frequency dominant band contains less deep reflection information of the target layer, but the image signal-to-noise ratio is high. Then, the optimal band is obtained according to the low-frequency dominant band and the high-frequency dominant band, and the reliability of deep reflection seismic data and the observability of the image are taken into account by combining the advantages of the low-frequency dominant band containing more deep reflection information of the target layer and the high-frequency dominant band having a higher image signal-to-noise ratio.

[0038] It should be understood that the foregoing general description and the following detailed description are exemplary and explanatory only and are not restrictive of the present application. BRIEF DESCRIPTION OF THE DRAWINGS

[0039] The above and other objects, features and advantages of the present application will become more apparent through a more detailed description of exemplary embodiments of the present application in conjunction with the accompanying drawings, wherein the same reference numerals generally represent the same components in the exemplary embodiments of the present application.

[0040] Figure 1 It is a flow chart of a method for selecting an optimal frequency band for deep reflection seismic data shown in an embodiment of the present application;

[0041] Figure 2 is a schematic diagram of a filter cross section shown in an embodiment of the present application;

[0042] Figure 3 is another schematic diagram of a filter cross section shown in an embodiment of the present application;

[0043] Figure 4 is a schematic diagram of a single-channel cross section after iterative filter banding shown in an embodiment of the present application;

[0044] Figure 5 is a schematic diagram of an original superimposed cross section shown in an embodiment of the present application;

[0045] Figure 6 is a schematic diagram of a superimposed cross section of a dominant frequency band shown in an embodiment of the present application;

[0046] Figure 7 It is a schematic diagram of the structure of an electronic device shown in an embodiment of the present application. DETAILED DESCRIPTION

[0047] The preferred embodiments of the present application will be described in more detail below with reference to the accompanying drawings. Although the preferred embodiments of the present application are shown in the accompanying drawings, it should be understood that the present application can be implemented in various forms and should not be limited by the embodiments described herein. On the contrary, these embodiments are provided to make the present application more thorough and complete, and to fully convey the scope of the present application to those skilled in the art.

[0048] The terms used in this application are for the purpose of describing specific embodiments only and are not intended to limit this application. The singular forms of "a", "said" and "the" used in this application and the appended claims are also intended to include plural forms unless the context clearly indicates other meanings. It should also be understood that the term "and / or" used herein refers to and includes any or all possible combinations of one or more associated listed items.

[0049] It should be understood that although the terms "first", "second", "third", etc. may be used in this application to describe various information, this information should not be limited to these terms. These terms are only used to distinguish the same type of information from each other. For example, without departing from the scope of this application, the first information may also be referred to as the second information, and similarly, the second information may also be referred to as the first information. Thus, the features defined as "first" and "second" may explicitly or implicitly include one or more of the features. In the description of this application, the meaning of "multiple" is two or more, unless otherwise clearly and specifically defined.

[0050] Embodiment 1

[0051] At present, deep reflection seismic data processing still uses full-band processing, which does not fully suppress high-frequency energy and cannot highlight the low-frequency information below the basement.

[0052] In view of the above problems, an embodiment of the present application provides a method for selecting an optimal frequency band for deep reflection seismic data. The technical solution of the embodiment of the present application is described in detail below in conjunction with the accompanying drawings.

[0053] Figure 1 This is a flow chart of the method for selecting the optimal frequency band for deep reflection seismic data shown in the embodiment of the present application. Figure 1 , the optimal frequency band selection method for deep reflection seismic data is as follows:

[0054] A1. Data sorting to obtain single-channel profiles.

[0055] The deep reflection seismic data are reordered and converted from common shot point gather data to common receiver point gather data to form a single-channel profile.

[0056] A2. Set the filter band to perform frequency division filtering on the single-channel profile.

[0057] Set the filter frequency band: divide the filter frequency band into 6 different filter frequency bands, namely, 1-6 Hz, 6-12 Hz, 12-24 Hz, 24-36 Hz, 36-60 Hz and above 60 Hz. Use 6 different filter frequency bands to perform frequency division filtering on the single-channel profile to obtain the filtered profile.

[0058] A3. Select the initial dominant frequency band based on the signal-to-noise ratio and contrast of the filtered profile.

[0059] The filter band corresponding to the filter profile with the highest signal-to-noise ratio and contrast is selected as the initial dominant frequency band.

[0060] For example, after comparison, it was found that the filtering profiles obtained using 1-6 Hz, 6-12 Hz, 12-24 Hz and 24-36 Hz frequency bands had better effects, and potential deep reflection signals were vaguely visible in the layers below the basement. The layers below the basement were taken as the target layers, and 1-36 Hz was selected as the initial dominant frequency band.

[0061] Preferably, energy balance adjustment is performed on multiple filter profiles to adjust the color weight of the filter profiles and highlight the target layer.

[0062] Figure 2 It is a schematic diagram of the filter profile shown in the embodiment of the present application, including schematic diagrams of the filter profiles corresponding to 1-6 Hz, 6-12 Hz and 12-24 Hz respectively. Figure 3 It is another schematic diagram of the filter profile shown in the embodiment of the present application, including schematic diagrams of the filter profiles corresponding to 24-36 Hz, 36-60 Hz and above 60 Hz respectively.

[0063] A4. Expand the middle frequency of the initial dominant frequency band to obtain an expanded dominant frequency band.

[0064] The average of the lowest frequency of the initial dominant frequency band and the highest frequency of the initial dominant frequency band is calculated to obtain an intermediate frequency. The correction frequency is subtracted from the intermediate frequency to obtain a minimum extended frequency. The correction frequency is added to the intermediate frequency to obtain a maximum extended frequency. The correction frequency is half of the intermediate frequency.

[0065] The lowest extension frequency is taken as the lowest frequency of the extended dominant frequency band, and the highest extension frequency is taken as the highest frequency of the extended dominant frequency band.

[0066] For example, the initial dominant frequency band is 1-36 Hz, the intermediate frequency is 18.5 Hz, and the correction frequency is 9.25 Hz. The lowest extended frequency is 18.5-9.25=9.25 Hz, and the highest extended frequency corresponding to the lowest extended frequency is 18.5+9.25=27.75 Hz. The extended dominant frequency band is 9.25-27.75 Hz.

[0067] A5. Obtain the low-frequency dominant frequency band and the high-frequency dominant frequency band.

[0068] The low-frequency dominant frequency band is a frequency band containing more target layer deep reflection information but having a lower image signal-to-noise ratio.

[0069] The high-frequency dominant frequency band is a frequency band that contains less deep reflection information of the target layer but has a higher image signal-to-noise ratio.

[0070] The lowest extended frequency is used as the highest frequency of the low-frequency dominant frequency band, and the lowest frequency of the initial dominant frequency band is used as the lowest frequency of the low-frequency dominant frequency band.

[0071] The highest extended frequency is used as the lowest frequency of the high-frequency dominant frequency band, and the highest frequency of the initial dominant frequency band is used as the highest frequency of the high-frequency dominant frequency band.

[0072] Exemplarily, the highest frequency of the low frequency dominant band is 9.25 Hz, the lowest frequency of the low frequency dominant band is 1 Hz, the low frequency dominant band is 1-9.25 Hz, the lowest frequency of the high frequency dominant band is 27.75 Hz, the highest frequency of the high frequency dominant band is 36 Hz, and the high frequency dominant band is 27.75-36 Hz.

[0073] A6. Obtain the optimal frequency band based on the low-frequency dominant frequency band and the high-frequency dominant frequency band.

[0074] The optimal frequency band is a frequency band that contains more deep reflection information of the target layer and has a higher signal-to-noise ratio. The optimal frequency band has the advantages of both low-frequency advantage band and high-frequency advantage band.

[0075] The average value of the lowest frequency of the low-frequency dominant frequency band and the highest frequency of the low-frequency dominant frequency band is calculated, and the average value is used as the intermediate frequency of the low-frequency dominant frequency band.

[0076] The average value of the lowest frequency of the high frequency dominant frequency band and the highest frequency of the high frequency dominant frequency band is calculated, and the average value is used as the intermediate frequency of the high frequency dominant frequency band.

[0077] The middle frequency of the low-frequency dominant frequency band is used as the lowest frequency of the optimal frequency band, and the middle frequency of the high-frequency dominant frequency band is used as the highest frequency of the optimal frequency band.

[0078] For example, the middle frequency of the low-frequency dominant frequency band is (1+9.25) / 2=5.125 Hz, the middle frequency of the high-frequency dominant frequency band is (27.75+36) / 2=31.875 Hz, and the optimal frequency band is 5.125-31.875 Hz.

[0079] Figure 4 is a schematic diagram of a single channel cross section after the iterative filter band shown in the embodiment of the present application, please refer to Figure 4 , Figure 4 The left is a schematic diagram of the filter profile obtained by filtering using the initial dominant frequency band. Figure 4 The figure in the figure is a schematic diagram of the filter profile obtained by filtering using the extended dominant frequency band. Figure 4 The right is a schematic diagram of the filter profile obtained by filtering using the optimal frequency band.

[0080] A7. Convert the filtering profile corresponding to the optimal frequency band into a superposition profile.

[0081] The stacked profile is obtained by breaking up the horizontal sequence of the single-channel profile while keeping the vertical time series of the single-channel profile relatively unchanged. The broken time series are then reorganized and superimposed to obtain the profile.

[0082] The single-channel profile is filtered using the optimal frequency band to obtain the filtered profile corresponding to the optimal frequency band. The horizontal sequence of the filtered profile is broken up, and the vertical time series of the filtered profile remains relatively unchanged. The broken up time series are then reorganized and superimposed to obtain the dominant frequency band superimposed profile.

[0083] The filtered profile is a single-channel profile after filtering according to the filter band.

[0084] Figure 5 is a schematic diagram of the original superimposed cross section shown in the embodiment of the present application, Figure 6 It is a schematic diagram of the superposition cross section of the advantageous frequency bands shown in the embodiment of the present application.

[0085] By comparison Figure 5 and Figure 6 It can be seen that the dominant frequency band stacked profile has more low-frequency information below the basement than the original stacked profile.

[0086] The embodiment of the present application performs frequency division filtering on a single-channel profile by setting a filter band, selects an initial dominant band according to the signal-to-noise ratio and contrast of the filtered profile, and narrows the range of the filter band. Then, the intermediate frequency of the initial dominant band is extended to obtain an extended dominant band. A low-frequency dominant band is obtained according to the lowest extended frequency and the lowest frequency of the initial dominant band. The low-frequency dominant band contains more deep reflection information of the target layer, but the image signal-to-noise ratio is low. A high-frequency dominant band is obtained according to the highest extended frequency and the highest frequency of the initial dominant band. The high-frequency dominant band contains less deep reflection information of the target layer, but the image signal-to-noise ratio is high. Then, the optimal band is obtained according to the low-frequency dominant band and the high-frequency dominant band, and the advantages of the low-frequency dominant band containing more deep reflection information of the target layer and the high-frequency dominant band image signal-to-noise ratio are combined to take into account the reliability of deep reflection seismic data and the observability of the image.

[0087] Embodiment 2

[0088] Corresponding to the aforementioned application function implementation method embodiment, the present application also provides an electronic device and corresponding embodiments.

[0089] Figure 7 It is a schematic diagram of the structure of an electronic device shown in an embodiment of the present application.

[0090] See also Figure 7 , the electronic device 1000 includes a memory 1010 and a processor 1020 .

[0091] The processor 1020 may be a central processing unit (CPU), or other general-purpose processors, digital signal processors (DSP), application-specific integrated circuits (ASIC), field-programmable gate arrays (FPGA), or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. A general-purpose processor may be a microprocessor or any conventional processor, etc.

[0092] The memory 1010 may include various types of storage units, such as system memory, read-only memory (ROM), and permanent storage devices. Among them, ROM can store static data or instructions required by the processor 1020 or other modules of the computer. The permanent storage device may be a readable and writable storage device. The permanent storage device may be a non-volatile storage device that does not lose the stored instructions and data even after the computer is powered off. In some embodiments, the permanent storage device uses a large-capacity storage device (such as a magnetic or optical disk, flash memory) as a permanent storage device. In some other embodiments, the permanent storage device may be a removable storage device (such as a floppy disk, optical drive). The system memory may be a readable and writable storage device or a volatile readable and writable storage device, such as a dynamic random access memory. The system memory may store some or all instructions and data required by the processor at runtime. In addition, the memory 1010 may include any combination of computer-readable storage media, including various types of semiconductor memory chips (DRAM, SRAM, SDRAM, flash memory, programmable read-only memory), and disks and / or optical disks may also be used. In some embodiments, the memory 1010 may include a readable and / or writable removable storage device, such as a laser disc (CD), a read-only digital versatile disc (such as a DVD-ROM, a double-layer DVD-ROM), a read-only Blu-ray disc, an ultra-density optical disc, a flash memory card (such as an SD card, a mini SD card, a Micro-SD card, etc.), a magnetic floppy disk, etc. The computer-readable storage medium does not include carrier waves and transient electronic signals transmitted wirelessly or wired.

[0093] The memory 1010 stores executable codes, and when the executable codes are processed by the processor 1020 , the processor 1020 can execute part or all of the methods described above.

[0094] The scheme of the present application has been described in detail above with reference to the accompanying drawings. In the above embodiments, the description of each embodiment has its own emphasis. For parts that are not described in detail in a certain embodiment, reference can be made to the relevant descriptions of other embodiments. Those skilled in the art should also be aware that the actions and modules involved in the description are not necessarily required for the present application. In addition, it can be understood that the steps in the method of the embodiment of the present application can be adjusted in order, merged and deleted according to actual needs, and the modules in the device of the embodiment of the present application can be merged, divided and deleted according to actual needs.

[0095] In addition, the method according to the present application may also be implemented as a computer program or a computer program product, which includes computer program code instructions for executing some or all of the steps in the above method of the present application.

[0096] Alternatively, the present application can also be implemented as a non-temporary machine-readable storage medium (or computer-readable storage medium, or machine-readable storage medium) on which executable code (or computer program, or computer instruction code) is stored. When the executable code (or computer program, or computer instruction code) is executed by a processor of an electronic device (or electronic device, server, etc.), the processor executes part or all of the steps of the above-mentioned method according to the present application.

[0097] Those skilled in the art will further appreciate that the various illustrative logical blocks, modules, circuits, and algorithm steps described in connection with the application herein may be implemented as electronic hardware, computer software, or combinations of both.

[0098] The flow chart and block diagram in the accompanying drawings show the possible architecture, function and operation of the system and method according to multiple embodiments of the present application. In this regard, each square box in the flow chart or block diagram can represent a part of a module, a program segment or a code, and the part of the module, the program segment or the code contains one or more executable instructions for realizing the specified logical function. It should also be noted that in some alternative implementations, the functions marked in the square box can also occur in a sequence different from that marked in the accompanying drawings. For example, two continuous square boxes can actually be executed substantially in parallel, and they can sometimes be executed in the opposite order, depending on the functions involved. It should also be noted that each square box in the block diagram and / or the flow chart, and the combination of the square boxes in the block diagram and / or the flow chart can be implemented with a dedicated hardware-based system that performs the specified function or operation, or can be implemented with a combination of dedicated hardware and computer instructions.

[0099] The embodiments of the present application have been described above, and the above description is exemplary, not exhaustive, and is not limited to the disclosed embodiments. Many modifications and changes will be apparent to those of ordinary skill in the art without departing from the scope and spirit of the described embodiments. The selection of terms used herein is intended to best explain the principles of the embodiments, practical applications, or improvements to the technology in the market, or to enable other persons of ordinary skill in the art to understand the embodiments disclosed herein.

Claims

1. A method for selecting the optimal frequency band of deep reflection seismic data, characterized in that: include: Rearrange the deep reflection seismic data, convert the common shot point gather data into common receiver point gather data, and form a single-channel profile; Performing frequency division filtering on the single-channel profile by setting a filter frequency band to obtain a filtered profile, and selecting an initial dominant frequency band according to a signal-to-noise ratio and a contrast ratio of the filtered profile; Expanding the middle frequency of the initial dominant frequency band to obtain an expanded dominant frequency band; The method of extending the intermediate frequency of the initial dominant frequency band comprises: calculating the intermediate frequency of the initial dominant frequency band, and then extending the intermediate frequency in the direction of low frequency and high frequency respectively; the intermediate frequency is the average of the lowest frequency of the initial dominant frequency band and the highest frequency of the initial dominant frequency band; the method of extending the intermediate frequency in the direction of low frequency and high frequency respectively comprises: subtracting the correction frequency from the intermediate frequency to obtain the lowest extension frequency; adding the correction frequency to the intermediate frequency to obtain the highest extension frequency; taking the lowest extension frequency as the lowest frequency of the extended dominant frequency band, and taking the highest extension frequency as the highest frequency of the extended dominant frequency band; the correction frequency is half of the intermediate frequency; the range of the extended dominant frequency band is from the lowest extension frequency to the highest extension frequency; obtaining a low-frequency dominant frequency band according to the lowest extension frequency and the lowest frequency of the initial dominant frequency band, and obtaining a high-frequency dominant frequency band according to the highest extension frequency and the highest frequency of the initial dominant frequency band; and obtaining an optimal frequency band according to the low-frequency dominant frequency band and the high-frequency dominant frequency band.

2. The method for selecting the optimal frequency band of deep reflection seismic data according to claim 1, characterized in that: The single-channel profile is subjected to frequency division filtering by setting a filter frequency band, comprising: Setting 6 filter frequency bands, wherein the filter frequency bands are respectively set to 1-6 Hz, 6-12 Hz, 12-24 Hz, 24-36 Hz, 36-60 Hz and above 60 Hz; The single-channel profile is subjected to frequency division filtering using the six filter frequency bands that are set respectively.

3. The method for selecting the optimal frequency band of deep reflection seismic data according to claim 1, characterized in that: The step of obtaining a low-frequency dominant frequency band according to the lowest extended frequency and the lowest frequency of the initial dominant frequency band comprises: The lowest extended frequency is used as the highest frequency of the low-frequency dominant frequency band, and the lowest frequency of the initial dominant frequency band is used as the lowest frequency of the low-frequency dominant frequency band.

4. The method for selecting the optimal frequency band of deep reflection seismic data according to claim 1, characterized in that: The step of obtaining a high frequency dominant frequency band according to the highest extended frequency and the highest frequency of the initial dominant frequency band comprises: The highest extended frequency is used as the lowest frequency of the high-frequency dominant frequency band, and the highest frequency of the initial dominant frequency band is used as the highest frequency of the high-frequency dominant frequency band.

5. The method for selecting the optimal frequency band of deep reflection seismic data according to claim 1, characterized in that: The obtaining of the optimal frequency band according to the low-frequency dominant frequency band and the high-frequency dominant frequency band comprises: Using the middle frequency of the low-frequency dominant frequency band as the lowest frequency of the optimal frequency band, and using the middle frequency of the high-frequency dominant frequency band as the highest frequency of the optimal frequency band; The middle frequency of the low frequency dominant frequency band is the average of the lowest frequency of the low frequency dominant frequency band and the highest frequency of the low frequency dominant frequency band; The middle frequency of the high frequency dominant frequency band is an average value of the lowest frequency of the high frequency dominant frequency band and the highest frequency of the high frequency dominant frequency band.

6. The method for selecting the optimal frequency band of deep reflection seismic data according to claim 1, characterized in that: The selecting an initial dominant frequency band according to the signal-to-noise ratio and contrast of the filtering profile comprises: Performing frequency division filtering on the single-channel profile to obtain a filtered profile; An initial dominant frequency band is selected, where the initial dominant frequency band is the filter frequency band corresponding to the filter profile having the highest signal-to-noise ratio and contrast.

7. An electronic device, characterized in that: include: processor; as well as A memory having executable codes stored thereon, which, when executed by the processor, causes the processor to execute the method according to any one of claims 1 to 6.

8. A non-transitory machine-readable storage medium, characterized in that: An executable code is stored thereon, and when the executable code is executed by a processor of an electronic device, the processor is caused to execute the method as claimed in any one of claims 1 to 6.

Citation Information

Patent Citations

  • Pre-stack seismic data frequency band expanding method of combining well information and horizon information

    CN106597532A

  • Band spreading method and device, electronic equipment and computer readable storage medium

    CN110556122A