Seismic data processing method, device, equipment, storage medium and program product

Abstract

The present application discloses a method, apparatus, device, storage medium, and program product for processing seismic data, relating to the field of geological exploration. The method comprises: acquiring first seismic data and second seismic data of a target area; wherein the first seismic data and the second seismic data are acquired in different ways, and the resolution of the first seismic data is higher than the resolution of the second seismic data; determining a set of matching factors for the first seismic data and the second seismic data in a target frequency domain based on the first seismic data and the second seismic data; and determining time-frequency domain shaping data for the second seismic data based on the second seismic data and the matching factor set, wherein the time-frequency domain shaping data is seismic data with a higher resolution than the second seismic data. The present application effectively broadens the frequency bandwidth of seismic data and significantly improves the resolution of seismic data.

Description

Technical Field

[0001] The present application relates to the field of geological exploration, and in particular to a method, apparatus, device, storage medium and program product for processing seismic data. Background Art

[0002] When studying underground strata, the use of high-resolution seismic data can reasonably restore the high-frequency and low-frequency information of seismic records, which is of great significance for studying more detailed issues such as distinguishing thin reservoirs, identifying small fault blocks, and identifying geological boundaries.

[0003] In related technologies, surface seismic data is generally used directly to study underground strata.

[0004] However, due to the low resolution of ground seismic data, the clarity of the geological structure reflected by the data is low. Summary of the Invention

[0005] The present invention provides a method, apparatus, device, storage medium, and program product for processing seismic data. The technical solution is as follows:

[0006] According to one aspect of an embodiment of the present application, a method for processing seismic data is provided, the method comprising:

[0007] Acquire first seismic data and second seismic data of a target area; wherein the first seismic data and the second seismic data are acquired in different ways, and a resolution of the first seismic data is higher than a resolution of the second seismic data;

[0008] determining a matching factor set for the first seismic data and the second seismic data in a target frequency domain based on the first seismic data and the second seismic data;

[0009] Time-frequency domain shaping data of the second seismic data is determined according to the second seismic data and the matching factor set, where the time-frequency domain shaping data is seismic data with a higher resolution than the second seismic data.

[0010] According to one aspect of an embodiment of the present application, a device for processing seismic data is provided, the device comprising:

[0011] a first acquisition module, configured to acquire first seismic data and second seismic data of a target area; wherein the first seismic data and the second seismic data are acquired in different ways, and the resolution of the first seismic data is higher than that of the second seismic data;

[0012] a first determining module, configured to determine a matching factor set of the first seismic data and the second seismic data in a target frequency domain based on the first seismic data and the second seismic data;

[0013] The second determination module is configured to determine time-frequency domain shaped data of the second seismic data based on the second seismic data and the matching factor set, wherein the time-frequency domain shaped data is seismic data with a higher resolution than the second seismic data.

[0014] According to one aspect of an embodiment of the present application, a computer device is provided, comprising a processor and a memory, wherein a computer program is stored in the memory, and the computer program is loaded and executed by the processor to implement the above-mentioned method for processing seismic data.

[0015] According to one aspect of an embodiment of the present application, a computer-readable storage medium is provided, in which a computer program is stored. The computer program is loaded and executed by a processor to implement the above-mentioned method for processing seismic data.

[0016] According to one aspect of an embodiment of the present application, a computer program product is provided, wherein the computer program product includes a computer program, and the computer program is loaded and executed by a processor to implement the above-mentioned method for processing seismic data.

[0017] The technical solutions provided in the embodiments of the present application can bring the following beneficial effects:

[0018] The technical solution provided in the present application calculates the matching factor between the first seismic data and the second seismic data, and uses the first seismic data with higher resolution to adjust the second seismic data with lower resolution, thereby obtaining adjusted time-frequency domain shaped data. Compared with the second seismic data before shaping, the shaped seismic data effectively broadens the frequency bandwidth of the seismic data and significantly improves the resolution of the seismic data. BRIEF DESCRIPTION OF THE DRAWINGS

[0019] Figure 1 This is a schematic diagram of an implementation environment for a solution provided by an embodiment of the present application;

[0020] Figure 2 is a flow chart of a method for processing seismic data provided by one embodiment of the present application;

[0021] Figure 3 This is a comparative cross-sectional diagram of VSP seismic data and surface seismic data provided by one embodiment of the present application;

[0022] Figure 4 This is a comparative cross-sectional diagram of time-frequency domain shaping data (time domain data volume) and VSP seismic data provided by one embodiment of the present application;

[0023] Figure 5 This is a comparative cross-sectional diagram of time-frequency domain shaping data (time domain data volume) and ground seismic data provided by an embodiment of the present application;

[0024] Figure 6 This is a spectrum characteristic diagram of VSP seismic data provided by an embodiment of the present application;

[0025] Figure 7 This is a frequency spectrum characteristic diagram of ground seismic data provided by an embodiment of the present application;

[0026] Figure 8 This is a spectrum feature diagram of time-frequency domain shaping data (time domain data volume) provided by an embodiment of the present application;

[0027] Figure 9 is a flow chart of a method for processing seismic data provided by one embodiment of the present application;

[0028] Figure 10 is a block diagram of a seismic data processing device provided by one embodiment of the present application;

[0029] Figure 11 This is a structural block diagram of a computer device provided in one embodiment of the present application. DETAILED DESCRIPTION

[0030] In order to make the objectives, technical solutions and advantages of this application clearer, the implementation methods of this application will be further described in detail below with reference to the accompanying drawings.

[0031] Before introducing the embodiments of the present application, in order to facilitate understanding of the present solution, the following explanations are given for the terms used in the present solution:

[0032] 1. Seismic data: refers to seismic signals received and recorded through seismic exploration projects. Seismic waves are artificially excited at the surface. As they propagate underground, they are reflected and refracted at interfaces between rock layers with different properties. Detectors are used to receive these reflected and refracted waves at the surface or in wells. The received seismic wave signals are related to the characteristics of the earthquake source, the location of the detection points, and the properties and structure of the underground rock formations through which the waves pass. By processing and interpreting these seismic wave signals, the properties and morphology of the underground rock formations can be inferred.

[0033] Vertical Seismic Profile (VSP) seismic data refers to the seismic waves generated by a surface source and received by geophones located deep within a well. This refers to the observed signal data obtained by observing a one-dimensional artificial field in the vertical direction. Surface seismic data refers to the seismic waves generated by a surface source and received by geophones located on the surface.

[0034] 2. Resolution of seismic data: It is the ability to distinguish various geological bodies and stratum details, which can reflect the clarity of various geological structures, including two aspects: vertical resolution and lateral resolution. Vertical resolution is also called vertical resolution or time resolution, which refers to the thickness of the thinnest stratum that can be distinguished by seismic data in the vertical direction. Generally, there are two meanings: one is to be able to correctly identify the reflection waves of the top and bottom interfaces of the stratum from the seismic data; the other is to be able to determine the reflection waves of thin strata from the stratum data, thereby determining the existence of thin underground layers. Lateral resolution is also called horizontal resolution or spatial resolution, which refers to the width of the smallest geological body that can be distinguished by seismic records in the horizontal direction. The resolution of seismic data described in this application refers to the longitudinal resolution of seismic data.

[0035] Please refer to Figure 1 , which shows a schematic diagram of an implementation environment provided by an exemplary embodiment of the present application. The implementation environment includes a computer device 110 and a data acquisition system 120.

[0036] The computer device 110 is an electronic device used to process, analyze and interpret seismic data. The computer device may be a device such as a PC (Personal Computer), a tablet computer, a smart phone, a wearable device, an intelligent robot, or the like.

[0037] Data acquisition system 120 includes an in-well receiving system 121 and a surface receiving system 122. In-well receiving system 121 includes an in-well geophone, and surface receiving system 122 includes a surface excitation device and a surface geophone. There can be one or more surface geophones, and there can be multiple in-well geophones, each at a different depth.

[0038] For example, a surface excitation device, a surface geophone, and multiple borehole geophones are installed in the target area. The surface excitation device excites seismic waves that propagate underground. After reflection and refraction at the interfaces of different rock formations, the seismic waves are received by the borehole geophones in the well and the surface geophones on the surface. The seismic data received by the borehole geophones is VSP seismic data, and the seismic data received by the surface geophones is surface seismic data. The received seismic data is sent to a computer 110, which processes and analyzes the VSP seismic data and surface seismic data of the target area to obtain shaped surface seismic data.

[0039] Please refer to Figure 2 , which shows a flowchart of a method for processing seismic data provided by an exemplary embodiment of the present application. In this embodiment, the method is mainly applied to Figure 1The method may include at least one of the following steps 210 to 230:

[0040] Step 210 , acquiring first seismic data and second seismic data of the target area; wherein the first seismic data and the second seismic data are acquired in different ways, and the resolution of the first seismic data is higher than the resolution of the second seismic data.

[0041] pass Figure 1 The data acquisition system 120 shown can acquire first seismic data and second seismic data of a target area, wherein the two types of seismic data are acquired using different acquisition methods. Furthermore, due to the different acquisition methods, the two types of seismic data have different resolutions, with the resolution of the first seismic data being higher than the resolution of the second seismic data.

[0042] For example, the first seismic data may be VSP seismic data, and the second seismic data may be surface seismic data. VSP seismic data is seismic data obtained by placing detectors at different depths in a well, while surface seismic data is seismic data obtained by placing detectors on the surface. Seismic data received by detectors placed in a well is more stable and accurate, while seismic data received by detectors placed on the surface is easily affected by other surface signals and other interferences, resulting in instability and disorder in the surface seismic data. Therefore, the resolution of VSP seismic data is higher than that of surface seismic data, but the method of obtaining VSP seismic data is more difficult and more expensive than the method of obtaining surface seismic data. Therefore, by using higher-resolution VSP seismic data to reshape lower-resolution surface seismic data, higher-resolution surface seismic data can be obtained, which can also be applied to other geological exploration areas.

[0043] refer to Figure 3 The VSP seismic data and the ground seismic data are compared in the cross-section. The horizontal axis is the ground distance and the vertical axis is the formation depth. The formation depth can be calculated based on the time of receiving the seismic data and the speed of seismic wave transmission in the formation. Therefore, the vertical axis can also be regarded as the target time domain of the seismic data. Figure 3 In the figure, VSP seismic data 301 and surface seismic data 302 are shown. It can be seen that the clarity of VSP seismic data 301 is higher than that of surface seismic data 302, that is, the resolution of VSP seismic data is higher than that of surface seismic data. Both VSP seismic data 301 and surface seismic data 302 are displayed as lines of varying depths. The depth of the lines reflects the signal energy received by the detectors. The variations in the lines can also reveal the geological conditions of the subsurface, such as the presence of rocks or altered layers.

[0044] Step 220 : Determine a set of matching factors between the first seismic data and the second seismic data in a target frequency domain based on the first seismic data and the second seismic data.

[0045] The target frequency domain refers to the frequency domain range of the first and second seismic data. For example, the matching factor set includes multiple matching factors, each of which corresponds to a frequency within the target frequency domain. A frequency corresponds to a first and a second seismic data. The shaped second seismic data can be obtained by using the first seismic data corresponding to that frequency and the matching factor.

[0046] In some embodiments, the convolution of the first seismic data and the reflection coefficient of the target area is calculated; the convolution of the second seismic data and the reflection coefficient is calculated; the time-frequency domain data of the first seismic data and the reflection coefficient and the time-frequency domain data of the second seismic data are obtained by performing wavelet transform on the convolution of the first seismic data and the reflection coefficient, and the convolution of the second seismic data and the reflection coefficient, respectively; and the matching factor set of the first seismic data and the second seismic data in the target frequency domain is determined based on the time-frequency domain data of the first seismic data and the time-frequency domain data of the second seismic data.

[0047] When seismic waves excited by the ground enter the medium vertically and reach the interface between rock layers, some of the energy is reflected, while the remaining energy passes through the interface. The ratio of the amplitude of the reflected wave to the amplitude of the incident wave is called the reflection coefficient. The magnitude of the reflection coefficient is determined by the impedance difference between the rock layers above and below the interface; a greater impedance difference results in a greater reflection coefficient. Since both the first and second seismic data were acquired from the target area, the reflection coefficients of the first and second seismic data can be considered equal.

[0048] According to seismic data, it is the convolution of seismic wavelet and reflection coefficient, that is:

[0049] X h (t) = W h (t)*R(t)

[0050] X l (t) = W l (t)*R(t)

[0051] Among them, W h (t) is the first seismic wavelet, W l (t) is the second seismic wavelet, R(t) is the reflection coefficient of the target area, X h (t) is the first seismic data, which is the convolution of the first seismic wavelet and the reflection coefficient, X l (t) is the second seismic data, which is the convolution of the second seismic wavelet and the reflection coefficient.

[0052] X h (t) and X l (t) Perform wavelet transform to obtain the time-frequency domain data of the first seismic data and the time-frequency domain data of the second seismic data, namely:

[0053] X h (f, t) = W h (f, t) × R(f, t)

[0054] X l (f, t) = W l (f, t) × R(f, t)

[0055] Among them, the value range of f is the target frequency domain, the value range of t is the target time domain, and X h (f, t) is the time-frequency domain data of the first seismic data, X l (f, t) is the time-frequency domain data of the second seismic data.

[0056] According to the obtained X h (f, t) and X l (f, t), the matching factor set of the first seismic data and the second seismic data in the target frequency domain can be derived, namely:

[0057] H(f, t) = X h (f, t) / X l (f, t) = W h (f, t) / W l (f, t)

[0058] Among them, H(f, t) is the matching factor set.

[0059] Through the above steps, a matching factor set between the first seismic data and the second seismic data can be obtained. The matching factors in the matching factor set correspond to a frequency in the target frequency domain respectively. Through the one-to-one corresponding matching factors, the one-to-one shaped second seismic data can be obtained.

[0060] Step 230 : Determine time-frequency domain shaped data of the second seismic data based on the second seismic data and the matching factor set. The time-frequency domain shaped data has a higher resolution than the second seismic data.

[0061] Exemplarily, based on the target frequency in the target frequency domain, the target matching factor corresponding to the target frequency and the target time-frequency domain data of the second seismic data are determined; based on the target matching factor and the target time-frequency domain data, the target time-frequency domain shaping data of the second seismic data are determined.

[0062] Applying the matching factor to the time-frequency domain data of the second seismic data can obtain the time-frequency domain shaped data of the second seismic data, namely:

[0063] X l,h (f, t) = X l (f, t) × H(f, t)

[0064] Among them, X l,h (f, t) is the time-frequency domain shaping data of the second seismic data.

[0065] By respectively determining the matching factor and the time-frequency domain shaping data corresponding to the frequency in the target frequency domain, and according to the above formula, the time-frequency domain shaping data of the second seismic data can be obtained. Figure 4 The comparison profile of the time-frequency domain shaping data and VSP seismic data is shown, and Figure 5 The comparison profile of the time-frequency domain shaping data and the ground seismic data shown in the figure shows that the clarity of the time-frequency domain shaping data is higher than that of the ground seismic data, that is, the resolution of the time-frequency domain shaping data is higher than that of the ground seismic data. Figure 5 The comparison shows that the continuity of the data lines after shaping is better, and the contact of the curved parts of the data lines is more obvious and clear, which can reflect that the continuity of the phase axis of the data after shaping has been improved, and the wave resistance characteristics have also been strengthened. Figure 4 and Figure 5 From the overall effect, the overall quality of the shaped ground seismic data has been greatly improved and enhanced.

[0066] For example, refer to Figure 6 、 Figure 7 and Figure 8 , the horizontal axis represents the frequency, and the vertical axis represents the energy ratio of the seismic data to the energy of the original excited seismic wave. Figure 6 The spectral characteristics of the VSP seismic data shown are compared with Figure 7 The spectrum characteristic diagram of the ground seismic data shown in FIG, VSP seismic data has a wider frequency band, richer frequency components, and is closer to zero phase. The spectrum characteristic diagram of the seismic data after shaping is as follows Figure 8 As shown, it is very close to zero phase. Since geological exploration is affected by many factors, it is almost impossible to make the seismic data completely zero phase. Computer processing can be used to make the seismic data tend to zero phase. In addition, Figure 8 The high-frequency and low-frequency information of the seismic data after the shaping process are very rich, and the bandwidth is also wider than Figure 7 The surface seismic data shown is wider. From the overall effect, the shaped seismic data effectively broadens the frequency bandwidth of the seismic data, improves the resolution of the seismic data, and enriches the frequency components of the seismic data in the target frequency domain, which is helpful for subsequent research using the shaped seismic data.

[0067] The technical solution provided in this application calculates the matching factor between first and second seismic data, using the higher-resolution first seismic data to adjust the lower-resolution second seismic data, thereby generating adjusted time-frequency domain shaped data. Compared to the unshaped second seismic data, the shaped seismic data effectively broadens the frequency bandwidth and significantly improves the resolution of the seismic data. Furthermore, from the perspective of the overall data profile, the continuity of the seismic data events is improved, the wave resistance characteristics are strengthened, and the overall quality of the profile is significantly improved.

[0068] On the other hand, in the related art, when shaping ground seismic data, it is necessary to extract and deconvolve the seismic wavelets of the ground seismic data. There are many uncertain factors in this process, which makes the processing process more complicated and the quality control more difficult. In addition, due to the many factors affecting the wavelet morphology, it is also difficult to accurately extract the seismic wavelets from the seismic data. Moreover, VSP seismic data is only used as prior information to compare and optimize ground seismic data, and does not participate in the processing and calculation of seismic wavelets. The present application applies VSP seismic data to the processing of ground seismic data. While obtaining seismic data with a wider bandwidth and higher resolution, it also optimizes the processing of ground seismic data, omits the step of extracting seismic wavelets, and makes the implementation of the solution simpler.

[0069] In some embodiments, time domain shaping data corresponding to the time-frequency domain shaping data is obtained by performing inverse transformation on the time-frequency domain shaping data; wherein the time domain shaping data refers to seismic data of the time-frequency domain shaping data in the target time domain.

[0070] By the above X l,h By performing an inverse transform on (f, t), we can obtain the corresponding time-domain shaped data, namely, the range of values ​​for the target time domain t. By acquiring the target time-domain data and combining it with the propagation velocity of the ground excitation wave in the subsurface rock formation, we can determine the stratigraphic depth of all seismic data received by the ground detectors. The waveform morphology, spectral characteristics, amplitude energy, and phase of the received seismic data facilitate the study of the geological morphology and stratigraphic conditions of different strata. The shaped seismic data facilitates more detailed and accurate study of the subsurface rock formations.

[0071] For example, Figure 9As shown in the flowchart, the technical solution provided by this application first acquires high-resolution VSP seismic data and lower-resolution surface seismic data within the target area. Data analysis is then performed on the VSP seismic data and the surface seismic data to obtain the convolution of the seismic data and the reflection coefficient. This convolution is then subjected to a wavelet transform to obtain time-frequency domain data for the VSP seismic data and the surface seismic data, respectively, within the target frequency domain. Based on the time-frequency domain data of the VSP seismic data and the surface seismic data, a set of matching factors for the VSP seismic data and the surface seismic data within the target frequency domain is then obtained, where each matching factor corresponds to a frequency within the target frequency domain.

[0072] Applying the matching factor set within the target frequency domain to the time-frequency domain data of lower-resolution ground seismic data yields time-frequency domain shaped data with higher resolution within the target frequency domain. This shaped data broadens the frequency bandwidth of the seismic data and improves its resolution. By performing an inverse transform on the time-frequency domain shaped data, the corresponding time-domain shaped data can be obtained. Analysis of this high-resolution time-domain shaped data facilitates the study of geological morphology and stratigraphic characteristics at different depths.

[0073] The following are device embodiments of the present application, which can be used to implement the method embodiments of the present application. For details not disclosed in the device embodiments of the present application, please refer to the method embodiments of the present application.

[0074] Please refer to Figure 10 , which shows a block diagram of a seismic data processing device provided by an embodiment of the present application. The device has the function of implementing the above-mentioned seismic data processing method, and the function can be implemented by hardware or by hardware executing corresponding software. The device can be the computer device described above, or it can be set in a computer device. Figure 10 As shown, the apparatus 1000 may include: a first acquisition module 1010 , a first determination module 1020 , and a second determination module 1030 .

[0075] The first acquisition module 1010 is used to acquire first seismic data and second seismic data of a target area; wherein the first seismic data and the second seismic data are acquired in different ways, and the resolution of the first seismic data is higher than that of the second seismic data.

[0076] The first determination module 1020 is configured to determine a matching factor set of the first seismic data and the second seismic data in a target frequency domain based on the first seismic data and the second seismic data.

[0077] The second determination module 1030 is configured to determine time-frequency domain shaped data of the second seismic data according to the second seismic data and the matching factor set, wherein the time-frequency domain shaped data has a higher resolution than the second seismic data.

[0078] In some embodiments, the first determining module 1020 is configured to:

[0079] Obtaining time-frequency domain data of the first seismic data and time-frequency domain data of the second seismic data by performing wavelet transform on the first seismic data and the second seismic data respectively;

[0080] A matching factor set of the first seismic data and the second seismic data in the target frequency domain is determined based on the time-frequency domain data of the first seismic data and the time-frequency domain data of the second seismic data.

[0081] In some embodiments, the matching factor set includes multiple matching factors, each matching factor corresponding to a frequency in the target frequency domain.

[0082] In some embodiments, the second determining module 1030 is configured to:

[0083] determining, according to a target frequency in the target frequency domain, a target matching factor corresponding to the target frequency and target time-frequency domain data of the second seismic data;

[0084] Target time-frequency domain shaping data of the second seismic data is determined according to the target matching factor and the target time-frequency domain data.

[0085] In some embodiments, the apparatus 1000 further includes:

[0086] The second acquisition module 1040 is used to obtain time domain shaping data corresponding to the time-frequency domain shaping data by performing inverse transformation on the time-frequency domain shaping data; wherein the time domain shaping data refers to seismic data of the time-frequency domain shaping data in the target time domain.

[0087] In some embodiments, the first seismic data is vertical seismic profile (VSP) seismic data, and the second seismic data is surface seismic data.

[0088] The technical solution provided in this application calculates the matching factor between first and second seismic data, using the higher-resolution first seismic data to adjust the lower-resolution second seismic data, thereby generating adjusted time-frequency domain shaped data. Compared to the unshaped second seismic data, the shaped seismic data effectively broadens the frequency bandwidth and significantly improves the resolution of the seismic data. Furthermore, from the perspective of the overall data profile, the continuity of the seismic data events is improved, the wave resistance characteristics are strengthened, and the overall quality of the profile is significantly improved.

[0089] It should be noted that the apparatus provided in the above embodiments, when implementing its functions, is only illustrated by the division of the above functional modules. In actual applications, the above functions can be assigned to different functional modules as needed, that is, the content structure of the device can be divided into different functional modules to complete all or part of the functions described above. In addition, the apparatus and method embodiments provided in the above embodiments are based on the same concept. The specific implementation process is detailed in the method embodiment and will not be repeated here.

[0090] Please refer to Figure 11 , which shows a block diagram of a terminal device 1100 provided in one embodiment of the present application. The terminal device 1100 can be any electronic device with data calculation, processing, and storage functions. The terminal device 1100 can be used to implement the seismic data processing method provided in the above embodiment.

[0091] Typically, the terminal device 1100 includes a processor 1101 and a memory 1102 .

[0092] The processor 1101 may include one or more processing cores, such as a 4-core processor, an 8-core processor, and the like. The processor 1101 may be implemented in at least one hardware form of DSP (Digital Signal Processing), FPGA (Field Programmable Gate Array), or PLA (Programmable Logic Array). The processor 1101 may also include a main processor and a coprocessor. The main processor is a processor for processing data in an awake state, also known as a CPU (Central Processing Unit); the coprocessor is a low-power processor for processing data in a standby state. In some embodiments, the processor 1101 may be integrated with a GPU (Graphics Processing Unit), which is responsible for rendering and drawing the content to be displayed on the display screen. In some embodiments, the processor 1101 may also include an AI processor for processing computing operations related to machine learning.

[0093] Memory 1102 may include one or more computer-readable storage media, which may be non-transitory. Memory 1102 may also include high-speed random access memory and non-volatile memory, such as one or more magnetic disk storage devices or flash memory storage devices. In some embodiments, the non-transitory computer-readable storage media in memory 1102 is used to store a computer program, which is configured to be executed by one or more processors to implement the above-described seismic data processing method.

[0094] Those skilled in the art will understand that Figure 11 The structure shown in the figure does not constitute a limitation on the terminal device 1100, and may include more or fewer components than shown in the figure, or combine certain components, or adopt a different component arrangement.

[0095] In an illustrative embodiment, a computer-readable storage medium is also provided, storing a computer program that, when executed by a processor of a terminal device, implements the aforementioned seismic data processing method. Optionally, the computer-readable storage medium may be a ROM (Read-Only Memory), a RAM (Random Access Memory), a CD-ROM (Compact Disc Read-Only Memory), a magnetic tape, a floppy disk, an optical data storage device, or the like.

[0096] In an exemplary embodiment, a computer program product is also provided, the computer program product including a computer program stored in a computer-readable storage medium. A processor of a terminal device reads the computer program from the computer-readable storage medium and executes the computer program, causing the terminal device to perform the above-described method for processing seismic data.

[0097] It should be understood that the "multiple" mentioned in this article refers to two or more. "And / or" describes the association relationship of associated objects, indicating that three relationships may exist. For example, A and / or B can represent three situations: A exists alone, A and B exist at the same time, and B exists alone. The character " / " generally indicates that the previous and subsequent associated objects are in an "or" relationship. In addition, the step numbers described in this article only illustrate a possible execution sequence between the steps. In some other embodiments, the above steps may not be executed in the order of the numbers, such as two steps with different numbers are executed at the same time, or two steps with different numbers are executed in the opposite order to the diagram. The embodiments of the present application do not limit this.

[0098] The above description is merely an exemplary embodiment of the present application and is not intended to limit the present application. Any modifications, equivalent substitutions, improvements, etc. made within the spirit and principles of the present application shall be included in the scope of protection of the present application.

Claims

1. A method for processing seismic data, characterized in that: The method comprises: Acquire first seismic data and second seismic data of a target area; wherein the first seismic data and the second seismic data are acquired in different ways, and a resolution of the first seismic data is higher than a resolution of the second seismic data; Obtaining time-frequency domain data of the first seismic data and time-frequency domain data of the second seismic data by performing wavelet transform on the first seismic data and the second seismic data, respectively; determining a matching factor set for the first seismic data and the second seismic data in a target frequency domain based on the time-frequency domain data of the first seismic data and the time-frequency domain data of the second seismic data; the matching factor set including a plurality of matching factors, each matching factor corresponding to a frequency in the target frequency domain; Determining time-frequency domain shaped data of the second seismic data based on the second seismic data and the matching factor set, wherein the time-frequency domain shaped data is seismic data with a higher resolution than the second seismic data; The step of determining the time-frequency domain shaped data of the second seismic data according to the second seismic data and the matching factor set includes: determining, according to a target frequency in the target frequency domain, a target matching factor corresponding to the target frequency and target time-frequency domain data of the second seismic data; Target time-frequency domain shaping data of the second seismic data is determined according to the target matching factor and the target time-frequency domain data.

2. The method according to claim 1, characterized in that The method further comprises: Obtaining time-domain shaping data corresponding to the time-frequency domain shaping data by performing inverse transformation on the time-frequency domain shaping data; The time-domain shaping data refers to the seismic data of the time-frequency domain shaping data in the target time domain.

3. The method according to claim 1, characterized in that The first seismic data is vertical seismic profile (VSP) seismic data, and the second seismic data is surface seismic data.

4. A seismic data processing device, characterized in that: The device comprises: a first acquisition module, configured to acquire first seismic data and second seismic data of a target area; wherein the first seismic data and the second seismic data are acquired in different ways, and the resolution of the first seismic data is higher than that of the second seismic data; a first determining module, configured to obtain time-frequency domain data of the first seismic data and time-frequency domain data of the second seismic data by performing wavelet transform on the first seismic data and the second seismic data, respectively; and determine a matching factor set of the first seismic data and the second seismic data in a target frequency domain based on the time-frequency domain data of the first seismic data and the time-frequency domain data of the second seismic data; the matching factor set including a plurality of matching factors, each matching factor corresponding to a frequency in the target frequency domain; a second determining module, configured to determine time-frequency domain shaped data of the second seismic data based on the second seismic data and the matching factor set, wherein the time-frequency domain shaped data is seismic data having a higher resolution than the second seismic data; The second determination module is also used to determine the target matching factor corresponding to the target frequency and the target time-frequency domain data of the second seismic data based on the target frequency in the target frequency domain; and determine the target time-frequency domain shaping data of the second seismic data based on the target matching factor and the target time-frequency domain data.

5. A computer device, characterized in that: The computer device includes a processor and a memory, wherein a computer program is stored in the memory, and the computer program is loaded and executed by the processor to implement the seismic data processing method according to any one of claims 1 to 3.

6. A computer-readable storage medium, characterized in that The computer-readable storage medium stores a computer program, which is loaded and executed by a processor to implement the method for processing seismic data according to any one of claims 1 to 3.

7. A computer program product, characterized in that The computer program product comprises a computer program, which is loaded and executed by a processor to implement the method for processing seismic data according to any one of claims 1 to 3.

Images