Weak bead reflection extraction method, system and equipment based on seismic section and medium

Abstract

The invention discloses a weak bead reflection extraction method, system and equipment based on a seismic section, and a medium, and belongs to the technical field of oil-gas exploration. The extraction method comprises the following steps: selecting a target work area, and obtaining original data of the target work area and a stratum containing bead reflection in the target work area; performing leveling treatment along the top layer of the stratum containing beaded reflection, and obtaining and stripping the leveled stratum to obtain an extraction area; analyzing the extraction area by adopting a principal component analysis method, and removing influence to obtain principal components in the extraction area; smoothing the principal component to eliminate the random component and obtain the smoothed principal component; carrying out anti-leveling treatment on the smoothed principal component to obtain a principal component subjected to the anti-leveling treatment; and obtaining secondary components in the extraction area based on the primary components after the anti-flattening processing and the original data of the target work area, highlighting bead reflection, and completing extraction.

Description

Technical Field

[0001] The present invention belongs to the technical field of oil and gas exploration, and particularly relates to a method, system, device and medium for extracting weak bead reflections based on seismic profiles. Background Art

[0002] Beads are the main types of oil and gas reservoirs in carbonate rock oilfields in the western region. The aggregate of beads and bead fractures is collectively called the fracture-vug body. Improving the prediction and identification accuracy of small-scale reservoirs such as beads and fracture-vug bodies is the key to further tapping potential production and reserves in western oilfields. In seismic exploration, the seismic records artificially excited and received by geophones contain both reflected wave information of continuous geological interfaces and rich diffracted wave information generated by rough interfaces and discontinuous interfaces. The former mainly reflects the characteristics of macroscopic geological bodies, and the latter mainly reflects the characteristics of small-scale geological bodies.

[0003] Beads are the reflection and characteristics of diffracted waves on seismic profiles. Their shapes and intensities often also reflect the properties and sizes of reservoirs, which is also the focus of carbonate rock oil and gas exploration.

[0004] How to better identify bead characteristics has also become a research hotspot in seismic exploration. Researchers usually use methods such as waveform decomposition to identify beads. However, such methods not only require parallel computing, are slow (3 - 7 days), but also have poor fidelity and poor effects, and the application effects are not good. And the weak bead reflections are masked by the reflected waves with strong energy on the seismic profile, and sometimes it is difficult to judge, which often brings troubles to seismic exploration. Summary of the Invention

[0005] In order to overcome the above-mentioned disadvantages of the prior art, the purpose of the present invention is to provide a method, system, device and medium for extracting weak bead reflections based on seismic profiles, extract bead information from seismic profiles, remove continuous strong in-phase axes in seismic data through principal component analysis, and then highlight bead reflections, so as to achieve the purpose of quickly identifying bead reflections and solve the technical problem of difficult identification of weak bead reflections.

[0006] In order to achieve the above purpose, the present invention adopts the following technical solutions:

[0007] The present invention provides a method for extracting weak bead reflections based on seismic profiles, including the following steps:

[0008] S1: Select a target work area, obtain the original data of the target work area and the formation containing bead reflections in the target work area;

[0009] S2: Perform flattening processing along the top layer of the formation containing bead reflections, obtain and strip the flattened formation to obtain an extraction area;

[0010] S3: Use the principal component analysis method to analyze the extraction area, remove the influence, and obtain the principal components in the extraction area;

[0011] S4: Smooth the principal components to eliminate the random components and obtain the smoothed principal components;

[0012] S5: Perform an inverse flattening process on the smoothed principal components to obtain the principal components after the inverse flattening process;

[0013] S6: Based on the principal components after the inverse flattening process and the original data of the target work area, obtain the secondary components in the extraction area, highlight the bead reflections, and complete the extraction.

[0014] In the specific implementation process, the original data of the target work area is a seismic profile.

[0015] In the specific implementation process, the principal components are parallel normal formations.

[0016] In the specific implementation process, the secondary components are beads.

[0017] In the specific implementation process, the removal of the influence is to remove the influence of the strong amplitude continuous event in the original data of the target work area on the bead-shaped reflection reservoir.

[0018] In the specific implementation process, the influence of the strong amplitude continuous event in the original data of the target work area on the bead-shaped reflection reservoir includes the suppression influence and the submergence influence of the strong amplitude continuous event on the bead-shaped reflection reservoir.

[0019] In the specific implementation process, the process of obtaining the secondary components in the extraction area is as follows:

[0020] Subtract the principal components after the inverse flattening process from the original data of the target work area adaptively to obtain the secondary components in the extraction area.

[0021] The present invention also provides a weak bead reflection extraction system based on a seismic profile, including:

[0022] An acquisition module, configured to select a target work area, obtain the original data of the target work area and the formation containing bead reflections in the target work area;

[0023] An extraction area obtaining module, configured to perform a flattening process along the top layer of the formation containing bead reflections, obtain and strip the flattened formation, and obtain an extraction area;

[0024] An analysis module, configured to use the principal component analysis method to analyze the extraction area, remove the influence, and obtain the principal components in the extraction area;

[0025] A smoothing processing module, configured to smooth the principal components to eliminate the random components and obtain the smoothed principal components;

[0026] An inverse flattening processing module for performing inverse flattening processing on the smoothed principal components to obtain the principal components after inverse flattening processing;

[0027] An extraction module for obtaining sub-components in the extraction area based on the principal components after inverse flattening processing and the original data of the target work area, highlighting bead reflections, and completing the extraction.

[0028] The present invention also provides an electronic device, including a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the computer program, the steps of any one of the weak bead reflection extraction methods based on seismic profiles are implemented.

[0029] The present invention also provides a computer-readable storage medium storing a computer program, and when the computer program is executed by a processor, the steps of any one of the weak bead reflection extraction methods based on seismic profiles are implemented.

[0030] Compared with the prior art, the present invention has the following beneficial effects:

[0031] The present invention provides a method for extracting weak bead reflections from a seismic profile, extracting bead information from the seismic profile, flattening the original seismic data, stripping it along the top and bottom horizons, applying the principal component analysis method to extract the bead emission characteristics, and performing adaptive subtraction with the original data after inverse flattening to obtain the bead reflections. It can successfully remove the suppression and inundation effects of strong-amplitude continuous coherent reflectors on the reservoir, and lay a solid data foundation for reservoir inversion. BRIEF DESCRIPTION OF THE DRAWINGS

[0032] Figure 1 It is a flowchart of the weak bead reflection extraction method based on seismic profiles according to an embodiment of the present invention;

[0033] Figure 2 It is a schematic diagram of the original data, principal components, and sub-components according to an embodiment of the present invention;

[0034] Figure 3 It is the original profile diagram of the target work area according to an embodiment of the present invention;

[0035] Figure 4 It is the flattened profile diagram of the top layer TO of Ordovician limestone according to an embodiment of the present invention; 3t Flattened profile display diagram;

[0036] Figure 5 It is the display diagram after stripping in the extraction area according to an embodiment of the present invention;

[0037] Figure 6 It is the display diagram after moderately smoothing the principal components according to an embodiment of the present invention;

[0038] Figure 7 The figure shows the flattened result of the smoothed principal components in the embodiments of the present invention;

[0039] Figure 8 The figure shows the beaded cross-section highlighting reflections in the embodiments of the present invention;

[0040] Figure 9 The figure shows the structural block diagram of the weak beaded reflection extraction system based on seismic profiles in the embodiments of the present invention. Detailed implementation manners

[0041] In order to enable those skilled in the art to better understand the solutions of the present invention, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention. Obviously, the described embodiments are only a part of the embodiments of the present invention, rather than all the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those of ordinary skill in the art without creative efforts shall fall within the protection scope of the present invention.

[0042] It should be noted that the terms "first", "second", etc. in the description and claims of the present invention and the above-mentioned drawings are used to distinguish similar objects, and do not necessarily need to be used to describe a specific order or sequence. It should be understood that such data can be interchanged under appropriate circumstances so that the embodiments of the present invention described herein can be implemented in an order other than those illustrated or described herein. In addition, the terms "comprising" and "having" and any variations thereof are intended to cover non-exclusive inclusion. For example, a process, method, system, product, or device that includes a series of steps or units does not necessarily have to be limited to those steps or units clearly listed, but may include other steps or units not clearly listed or inherent to these processes, methods, products, or devices.

[0043] The present invention will be further described in detail below with reference to the accompanying drawings:

[0044] The present invention provides a method, system, device, and medium for extracting weak beaded reflections from seismic profiles. The method extracts beaded information from seismic profiles, removes continuous strong in-phase axes in seismic data through principal component analysis, and then highlights beaded reflections, achieving the purpose of quickly identifying beaded reflections and solving the technical problem of difficult identification of weak beaded reflections.

[0045] The above method for extracting weak beaded reflections is a brand-new method for separating reflected waves and diffracted wave fields, namely principal component analysis. Principal component analysis is a dimensionality reduction analysis method, which is a statistical analysis method for extracting a few main comprehensive variables from the original multiple variables. That is, through orthogonal transformation, a set of variables that may be correlated are converted into a set of linearly uncorrelated variables, and this set of variables after transformation is called the principal components.

[0046] Principal component analysis was first introduced by K. Pearson for non-random variables, and later H. Hotelling extended this method to the case of random vectors. Therefore, principal component analysis is more suitable for analyzing multi-variable problems where there is a certain correlation between variables. By means of dimensionality reduction, multiple correlated variables can be reduced to a few variables to find out the most important comprehensive variables affecting the research problem. When there are thousands of sample data and thousands of dimensions, principal component analysis is a very effective dimensionality reduction method, and it can play two roles: one is to save storage space. When there is too much data, a lot of space can be saved by reducing a few dimensions; the other is to provide the calculation speed. After the data is reduced in dimension, whether it is during the training of sample data or when responding to new data, the speed will be greatly improved.

[0047] Principal component analysis is applicable to data with strong correlations between variables. If the original data has weak correlations, it cannot play a good role in dimensionality reduction. After dimensionality reduction, there is a small amount of information loss and it is impossible to contain 100% of the original data. The so-called dimensionality reduction is to reduce the number of correlated variables and replace the original variables with fewer variables. If the original variables are orthogonal, that is, there is no correlation, then principal component analysis has no effect.

[0048] The above-mentioned principal component analysis method for removing the influence of continuous event axes on beaded reflections in seismic data is based on the following understanding:

[0049] Seismic data consists of principal components (nearly parallel normal strata) and secondary components (such as beads). After removing the principal components, the remaining beaded part can be highlighted. Moreover, the better the parallelism of the strata, the stronger the correlation and the better the effect. The practical application of this technology in the Tarim Oilfield has significantly improved the detection ability and recognition accuracy of fracture-cavity bodies.

[0050] The first aspect of the present invention provides a method for extracting weak beaded reflections based on a seismic profile, including the following steps:

[0051] S1: Select a target work area, obtain the original data of the target work area and the strata containing beaded reflections in the target work area;

[0052] S2: Perform flattening processing along the top layer of the strata containing beaded reflections, obtain and strip the flattened strata to obtain an extraction area;

[0053] S3: Use principal component analysis to analyze the extraction area, remove the influence, and obtain the principal components in the extraction area;

[0054] S4: Smooth the principal components to eliminate random components and obtain the smoothed principal components;

[0055] S5: Perform an inverse flattening process on the smoothed principal components to obtain the principal components after the inverse flattening process;

[0056] S6: Based on the principal components after the inverse flattening process and the original data of the target work area, obtain the secondary components in the extraction area, highlight the bead-like reflections, and complete the extraction.

[0057] In the specific implementation process, the original data of the target work area is a seismic profile; the principal components are parallel normal formations; the secondary components are beads.

[0058] The above-mentioned removal of influence is to remove the influence of strong-amplitude continuous coherent horizons in the original data of the target work area on the bead-like reflection reservoir.

[0059] More specifically, the influence of strong-amplitude continuous coherent horizons in the original data of the target work area on the bead-like reflection reservoir includes the suppression influence and the submergence influence of the strong-amplitude continuous coherent horizons on the bead-like reflection reservoir.

[0060] More specifically, the method for extracting weak bead-like reflections based on a seismic profile is as follows:

[0061] Step 1: Select the target work area, obtain the seismic profile of the target work area and the formation containing bead-like reflections in the target work area;

[0062] Step 2: Perform a flattening process along the top layer of the formation containing bead-like reflections, obtain and strip the flattened formation, and obtain the extraction area;

[0063] Step 3: Use the principal component analysis method to analyze the extraction area, remove the suppression influence and the submergence influence of strong-amplitude continuous coherent horizons on the bead-like reflection reservoir, and obtain the principal components in the extraction area;

[0064] Step 4: Smooth the principal components to eliminate random components, and obtain the smoothed principal components, i.e., parallel normal formations;

[0065] Step 5: Perform an inverse flattening process on the smoothed principal components to obtain the principal components after the inverse flattening process;

[0066] Step 6: Based on the principal components after the inverse flattening process and the original data of the target work area, perform an adaptive subtraction between the principal components after the inverse flattening process and the original data of the target work area, obtain the secondary components in the extraction area, i.e., beads, highlight the bead-like reflections, and complete the extraction.

[0067] The second aspect of the present invention provides a system for extracting weak bead-like reflections based on a seismic profile, including an acquisition module, an extraction area acquisition module, an analysis module, a smoothing processing module, an inverse flattening processing module, and an extraction module.

[0068] Among them, the acquisition module is used to select a target work area, obtain the original data of the target work area, and the formation containing bead reflections in the target work area;

[0069] The acquisition area module is used to perform flattening processing along the top layer of the formation containing bead reflections, obtain and strip the flattened formation, and obtain the extraction area;

[0070] The analysis module is used to analyze the extraction area by using the principal component analysis method, remove the influence, and obtain the principal components in the extraction area;

[0071] The smoothing processing module is used to smooth the principal components to eliminate the random components and obtain the smoothed principal components;

[0072] The unflattening processing module is used to perform unflattening processing on the smoothed principal components to obtain the unflattened principal components;

[0073] The extraction module is based on the unflattened principal components and the original data of the target work area to obtain the secondary components in the extraction area, highlight the bead reflections, and complete the extraction.

[0074] The third aspect of the present invention provides an electronic device, including a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the computer program, the steps of any one of the weak bead reflection extraction methods based on seismic profiles are implemented.

[0075] The fourth aspect of the present invention provides a computer-readable storage medium. The computer-readable storage medium stores a computer program. When the computer program is executed by a processor, the steps of any one of the weak bead reflection extraction methods based on seismic profiles are implemented.

[0076] Embodiment

[0077] To prove the effectiveness and adaptability of the present invention, in the embodiment, a work area of carbonate rocks in the Tarim Basin is selected. Figure 3 The seismic profile of the work area, that is, the original data, shows that there are many bead reflections in the nearly parallel formations of the Ordovician interior on the profile.

[0078] This embodiment provides a weak bead reflection extraction method based on seismic profiles, as Figure 1 shown. The specific process is as follows:

[0079] Step 1: Select the target work area and find the formation containing bead information, that is, the target work area.

[0080] The bead reservoirs in the carbonate rocks of the Tarim Basin are mainly concentrated in the Ordovician limestone, and the Ordovician interior formations in the northern Tarim area are relatively parallel, which is suitable for the application of this technology. For the top of the Ordovician limestone formation (TO 3t)、 Bottom (TO 1y ) is explained (as Figure 3 shown);

[0081] Step 2: Flatten along the top layer position TO of the formation containing beaded information 3t and it can be seen that the formations are approximately parallel, indicating good data correlation (as Figure 4 shown). And strip the data in the target area along the top and bottom layer positions of the formation containing beaded information (as Figure 5 shown) to obtain the extraction area.

[0082] Step 3: Apply the principal component analysis method to extract the principal components of the data in the extraction area. Here, it is necessary to test the key parameter - the number of eigenvectors. The more eigenvectors, the closer it is to the original data. The initial parameter is selected starting from 10 - that is, 10 eigenvectors are used to represent the principal components of the data. And the accurate representation of the original data may require hundreds of eigenvectors. This process is the dimensionality reduction.

[0083] Step 4: Moderately smooth the extracted principal components to eliminate the random components (the part of the energy that may contain beads) (as Figure 6 shown).

[0084] Step 5: Unflatten the moderately smoothed principal components (as Figure 7 shown) to obtain the unflattened principal components.

[0085] Step 6: Subtract the unflattened principal components from the original data adaptively to obtain the data highlighting the beads (as Figure 8 shown).

[0086] Refer to Figure 9 , in another embodiment of the present invention, a weak bead reflection extraction system based on a seismic profile is provided, which can be used to implement the above-mentioned weak bead reflection extraction method based on a seismic profile. Specifically, the system includes an acquisition module, an extraction area acquisition module, an analysis module, a smoothing processing module, an unflattening processing module, and an extraction module.

[0087] Among them, the acquisition module is used to select a target work area, obtain the original data of the target work area and the formation containing bead reflections in the target work area;

[0088] The extraction area acquisition module is used to perform flattening processing along the top layer of the formation containing bead reflections, obtain and strip the flattened formation, and obtain the extraction area;

[0089] The analysis module is used to analyze the extraction area by using the principal component analysis method, remove the influence, and obtain the principal components in the extraction area;

[0090] The smoothing processing module is used to smooth the principal components to eliminate the random components and obtain the smoothed principal components;

[0091] An inverse flattening processing module, configured to perform inverse flattening processing on the smoothed principal components to obtain the principal components after inverse flattening processing;

[0092] An extraction module, based on the principal components after inverse flattening processing and the original data of the target work area, obtains the secondary components in the extraction area, highlights the bead reflections, and completes the extraction.

[0093] All relevant contents of each step involved in the embodiment of the foregoing weak bead reflection extraction method based on seismic profiles can be cited in the function description of the corresponding functional modules of the weak bead reflection extraction system based on seismic profiles in the embodiments of the present invention, and will not be elaborated herein.

[0094] The division of modules in the embodiments of the present invention is illustrative, and is only a logical function division. In actual implementation, there may be other division methods. In addition, in each embodiment of the present invention, the functional modules can be integrated in a processor, or can exist separately physically, or two or more modules can be integrated in one module. The above integrated modules can be implemented in the form of hardware or in the form of software functional modules.

[0095] In another embodiment of the present invention, a computer device is provided. The computer device includes a processor and a memory. The memory is used to store a computer program. The computer program includes program instructions. The processor is used to execute the program instructions stored in the computer storage medium. The processor may be a Central Processing Unit (CPU), or may also be other general-purpose processors, Digital Signal Processors (DSPs), Application Specific Integrated Circuits (ASICs), Field-Programmable Gate Arrays (FPGAs), or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. It is the computing core and control core of the terminal, and is suitable for implementing one or more instructions. Specifically, it is suitable for loading and executing one or more instructions in the computer storage medium to implement the corresponding method flow or corresponding function; the processor described in the embodiments of the present invention can be used for the operation of the weak bead reflection extraction method based on seismic profiles.

[0096] In another embodiment of the present invention, the present invention further provides a storage medium, specifically a computer-readable storage medium (Memory). The computer-readable storage medium is a memory device in a computer device, used to store programs and data. It can be understood that the computer-readable storage medium here can include both the built-in storage medium in the computer device and, of course, the extended storage medium supported by the computer device. The computer-readable storage medium provides a storage space, and this storage space stores the operating system of the terminal. And, in this storage space, there is also stored one or more instructions suitable for being loaded and executed by the processor. These instructions can be one or more computer programs (including program codes). It should be noted that the computer-readable storage medium here can be a high-speed RAM memory or a non-volatile memory, such as at least one disk memory. One or more instructions stored in the computer-readable storage medium can be loaded and executed by the processor to implement the corresponding steps of the method for extracting weak bead reflections based on seismic profiles in the above embodiments.

[0097] Those skilled in the art should understand that the embodiments of the present invention can be provided as a method, a system, or a computer program product. Therefore, the present invention can take the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware aspects. Moreover, the present invention can take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk memories, CD-ROMs, optical memories, etc.) containing computer-usable program codes.

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

[0099] These computer program instructions can also be stored in a computer-readable memory that can direct a computer or other programmable data processing device to work in a specific manner, so that the instructions stored in the computer-readable memory generate a manufactured article including an instruction device, and this instruction device implements the functions in Figure 1 one or more flows and / or blocksFigure 1 The functions specified in one or more boxes.

[0100] These computer program instructions can also be loaded onto a computer or other programmable data processing device, so that a series of operation steps are executed on the computer or other programmable device to generate a computer-implemented process. Thus, the instructions executed on the computer or other programmable device provide for implementing the steps of the functions specified in one Figure 1 process or multiple processes and / or boxes Figure 1 or more boxes.

[0101] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention and are not intended to limit them. Although the present invention has been described in detail with reference to the above embodiments, those of ordinary skill in the art should understand that: the specific implementation manners of the present invention can still be modified or equivalently replaced, and any modification or equivalent replacement that does not depart from the spirit and scope of the present invention shall be covered by the protection scope of the claims of the present invention.

Claims

1. A method for extracting weak bead reflections based on seismic profiles, characterized in that, It includes the following steps: S1: Select a target work area, obtain the original data of the target work area and the formation containing bead reflections in the target work area; S2: Flatten along the top layer of the formation containing bead reflections, obtain and strip the flattened formation to obtain an extraction area; S3: Analyze the extraction area using the principal component analysis method, remove the influence, and obtain the principal components in the extraction area; S4: Smooth the principal components to eliminate random components and obtain the smoothed principal components; S5: Perform an inverse flattening process on the smoothed principal components to obtain the principal components after the inverse flattening process; S6: Based on the principal components after the inverse flattening process and the original data of the target work area, obtain the secondary components in the extraction area, highlight the bead reflections, and complete the extraction.

2. The method for extracting weak bead reflections based on seismic profiles according to claim 1, wherein, The original data of the target work area is a seismic profile.

3. The method for extracting weak bead reflections based on seismic profiles according to claim 1, characterized in that The principal components are parallel normal formations.

4. The method for extracting weak bead reflections based on seismic profiles according to claim 1, wherein The secondary components are beads.

5. The method for extracting weak bead reflections based on seismic profiles according to claim 1, characterized in that The removal of the influence is to remove the influence of strong-amplitude continuous coherent events in the original data of the target work area on the bead-shaped reflection reservoir.

6. The method for extracting weak string of beads reflections based on seismic profiles according to claim 5, characterized in that, The influence of the strong-amplitude continuous coherent events in the original data of the target work area on the bead-shaped reflection reservoir includes the suppression influence and the inundation influence of the strong-amplitude continuous coherent events on the bead-shaped reflection reservoir.

7. The method for extracting weak bead reflections based on seismic profiles according to claim 1, characterized in that, The process of obtaining the secondary components in the extraction area is as follows: The principal components after the inverse flattening process and the original data of the target work area are adaptively subtracted to obtain the secondary components in the extraction area.

8. A weak bead reflection extraction system based on seismic profiles, characterized in that, It includes: An acquisition module, configured to select a target work area, obtain the original data of the target work area and the formation containing bead reflections in the target work area; An extraction area acquisition module, configured to flatten along the top layer of the formation containing bead reflections, obtain and strip the flattened formation to obtain an extraction area; An analysis module, configured to analyze the extraction area using the principal component analysis method, remove the influence, and obtain the principal components in the extraction area; A smoothing processing module, configured to smooth the principal components to eliminate random components and obtain the smoothed principal components; An inverse flattening processing module, configured to perform an inverse flattening process on the smoothed principal components to obtain the principal components after the inverse flattening process; An extraction module, based on the principal components after the inverse flattening process and the original data of the target work area, obtains the secondary components in the extraction area, highlights the bead reflections, and completes the extraction.

9. An electronic device, including a memory, a processor, and a computer program stored in the memory and executable on the processor, where when the processor executes the computer program, it implements the steps of the weak bead reflection extraction method based on a seismic profile according to any one of claims 1-7.

10. A computer-readable storage medium, where the computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, it implements the steps of the weak bead reflection extraction method based on a seismic profile according to any one of claims 1-7.

Images