CN117991345A - Method, device, equipment, storage medium and program product for processing seismic data - Google Patents
Method, device, equipment, storage medium and program product for processing seismic data Download PDFInfo
- Publication number
- CN117991345A CN117991345A CN202211385799.1A CN202211385799A CN117991345A CN 117991345 A CN117991345 A CN 117991345A CN 202211385799 A CN202211385799 A CN 202211385799A CN 117991345 A CN117991345 A CN 117991345A
- Authority
- CN
- China
- Prior art keywords
- seismic data
- data
- frequency domain
- seismic
- time
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 47
- 238000012545 processing Methods 0.000 title claims abstract description 38
- 238000007493 shaping process Methods 0.000 claims abstract description 54
- 238000004590 computer program Methods 0.000 claims description 16
- 230000009466 transformation Effects 0.000 claims description 8
- 230000001131 transforming effect Effects 0.000 claims description 2
- 230000015572 biosynthetic process Effects 0.000 description 12
- 238000005755 formation reaction Methods 0.000 description 12
- 238000010586 diagram Methods 0.000 description 10
- 230000008569 process Effects 0.000 description 8
- 230000009286 beneficial effect Effects 0.000 description 5
- 230000005284 excitation Effects 0.000 description 5
- 230000006870 function Effects 0.000 description 5
- 239000011435 rock Substances 0.000 description 5
- 238000001228 spectrum Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 4
- 238000003672 processing method Methods 0.000 description 3
- 230000003595 spectral effect Effects 0.000 description 3
- 230000000052 comparative effect Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 230000005055 memory storage Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 238000003908 quality control method Methods 0.000 description 1
- 238000009877 rendering Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
The application discloses a method, a device, equipment, a storage medium and a program product for processing seismic data, and relates to the field of geological exploration. The method comprises the following steps: acquiring first seismic data and second seismic data of a target area; the first seismic data and the second seismic data are acquired in different modes, and the resolution of the first seismic data is higher than that of the second seismic data; determining a matching factor set of the first seismic data and the second seismic data in a target frequency domain according to the first seismic data and the second seismic data; and determining time-frequency domain shaping data of the second seismic data according to the second seismic data and the matching factor set, wherein the time-frequency domain shaping data is seismic data with higher resolution than the second seismic data. The application effectively widens the frequency bandwidth of the seismic data and obviously improves the resolution of the seismic data.
Description
Technical Field
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
When the underground stratum is researched, high-frequency and low-frequency information of the seismic records can be reasonably recovered by using high-resolution seismic data, and the method has great significance in researching the finer problems of distinguishing thin reservoirs, identifying small fault blocks, distinguishing geological boundaries and the like.
In the related art, the subsurface formation is generally studied directly by using the surface seismic data.
However, due to the lower resolution of the surface seismic data, the clarity of the geologic structure reflected by the data is lower.
Disclosure of Invention
The embodiment of the application provides a method, a device, equipment, a storage medium and a program product for processing seismic data. The technical scheme is as follows:
according to an aspect of an embodiment of the present application, there is provided a method for processing seismic data, the method including:
Acquiring first seismic data and second seismic data of a target area; the first seismic data and the second seismic data are acquired in different manners, and the resolution of the first seismic data is higher than that of the second seismic data;
determining a matching factor set of the first seismic data and the second seismic data in a target frequency domain according to the first seismic data and the second seismic data;
And determining time-frequency domain shaping data of the second seismic data according to the second seismic data and the matching factor set, wherein the time-frequency domain shaping data is seismic data with higher resolution than the second seismic data.
According to an aspect of an embodiment of the present application, there is provided a seismic data processing apparatus, the apparatus including:
the first acquisition module is used for acquiring first seismic data and second seismic data of the target area; the first seismic data and the second seismic data are acquired in different manners, and the resolution of the first seismic data is higher than that of the second seismic data;
The first determining module is used for determining a matching factor set of the first seismic data and the second seismic data in a target frequency domain according to the first seismic data and the second seismic data;
And the second determining module is used for determining time-frequency domain shaping data of the second seismic data according to the second seismic data and the matching factor set, wherein the time-frequency domain shaping data is seismic data with higher resolution than the second seismic data.
According to an aspect of an embodiment of the present application, there is provided a computer apparatus including a processor and a memory, the memory storing a computer program loaded and executed by the processor to implement the method of processing seismic data described above.
According to an aspect of an embodiment of the present application, there is provided a computer-readable storage medium having stored therein a computer program loaded and executed by a processor to implement the above-described seismic data processing method.
According to an aspect of an embodiment of the present application, there is provided a computer program product comprising a computer program loaded and executed by a processor to implement the method of processing seismic data described above.
The technical scheme provided by the embodiment of the application can bring the following beneficial effects:
According to the technical scheme provided by the application, the matching factors of the first seismic data and the second seismic data are calculated, the first seismic data with higher resolution is used for adjusting the second seismic data with lower resolution, so that the adjusted time-frequency domain shaping data are obtained, and compared with the second seismic data before shaping, the shaped seismic data effectively widen the frequency bandwidth of the seismic data, and the resolution of the seismic data is obviously improved.
Drawings
FIG. 1 is a schematic illustration of an implementation environment for an embodiment of the present application;
FIG. 2 is a flow chart of a method of processing seismic data provided by one embodiment of the application;
FIG. 3 is a cross-sectional view of VSP seismic data and surface seismic data as provided by one embodiment of the application;
FIG. 4 is a comparative cross-sectional view of time-frequency domain shaping data (time domain data volume) and VSP seismic data provided by one embodiment of the application;
FIG. 5 is a cross-sectional view of time-frequency domain shaping data (time domain data volume) and surface seismic data, as provided by one embodiment of the present application;
FIG. 6 is a spectral signature of VSP seismic data provided by one embodiment of the application;
FIG. 7 is a spectral signature of surface seismic data provided by one embodiment of the application;
FIG. 8 is a diagram of spectral characteristics of time-frequency domain shaping data (time domain data volume) provided by one embodiment of the present application;
FIG. 9 is a flow chart of a method of processing seismic data provided by one embodiment of the application;
FIG. 10 is a block diagram of a seismic data processing device provided in one embodiment of the application;
FIG. 11 is a block diagram of a computer device according to one embodiment of the present application.
Detailed Description
For the purpose of making the objects, technical solutions and advantages of the present application more apparent, the embodiments of the present application will be described in further detail with reference to the accompanying drawings.
Before describing embodiments of the present application, in order to facilitate understanding of the present solution, the following explanation is made on terms appearing in the present solution:
1. Seismic data: refers to seismic signals received and recorded by seismic exploration engineering. When the earthquake wave is excited by artificial method on the earth surface, the earthquake wave can be reflected and refracted by the rock stratum interfaces with different medium properties when propagating to the underground, and the reflected and refracted earthquake wave can be received in the earth surface or the well by using the wave detector. The received seismic signals are related to the nature of the source, the location of the detector points, the nature and structure of the subsurface formation through which the seismic waves pass, and by processing and interpreting the seismic signals, the nature and morphology of the subsurface formation can be inferred.
VSP (VERTICAL SEISMIC Profile ) seismic data refers to seismic waves received by a set of detectors in a well with a source excitation at the surface. Namely, the one-dimensional artificial field is observed in the vertical direction, and observed signal data is obtained. The ground seismic data refers to seismic waves which are excited by a seismic source arranged on the ground surface and received by a wave detector arranged on the ground surface.
2. Resolution of seismic data: the method is used for distinguishing various geological bodies and stratum details, and can reflect the definition of various geological structures, including longitudinal resolution and transverse resolution. Longitudinal resolution, also referred to as vertical resolution or time resolution, refers to the thickness of the thinnest formation that the seismic data can resolve in the vertical direction. There are generally two meanings: one is a reflected wave that can accurately identify the formation top and bottom interfaces from the seismic data; another is that reflected waves of the thin formations can be determined from the formation data to determine the presence of the thin layers in the subsurface. The lateral resolution is also referred to as horizontal resolution or spatial resolution, which refers to the width of the smallest volume of land that a seismic record can resolve in the horizontal direction. The resolution of the seismic data depicted in the present application refers to the longitudinal resolution of the seismic data.
Referring to FIG. 1, a schematic diagram of an implementation environment provided by an exemplary embodiment of the present application is shown. Included in this implementation environment are computer devices 110 and data acquisition systems 120.
The computer device 110 is an electronic device for processing, analyzing, and interpreting seismic data, which may be a device such as a PC (Personal Computer ), tablet, smart phone, wearable device, smart robot, or the like.
The data acquisition system 120 includes an in-well receiving system 121 and a surface receiving system 122, the in-well receiving system 121 including an in-well detector, and the surface receiving system 122 including a surface excitation device and a surface detector. The surface detectors can be arranged in one or more than one, the downhole detectors can be arranged in a plurality of the downhole detectors, and the depths of the downhole detectors are different.
Illustratively, a surface excitation device, a surface detector, and a plurality of borehole detectors are disposed in the target area. The ground excitation device excites the seismic waves to spread underground, the seismic waves can be received by an in-well detector arranged in a well and a ground detector arranged on the ground after being reflected and refracted by rock stratum interfaces with different properties of the underground, the seismic data received by the in-well detector is VSP seismic data, and the seismic data received by the ground detector is ground seismic data. The received seismic data is sent to the computer device 110, and the computer device 110 processes and analyzes the VSP seismic data and the ground seismic data of the target area to obtain shaped ground seismic data.
Referring to fig. 2, a flowchart of a method for processing seismic data according to an exemplary embodiment of the application is shown. In the present embodiment, the method is mainly applied to the computer device shown in fig. 1 for illustration. The method may include at least one of the following steps 210-230:
Step 210, acquiring first seismic data and second seismic data of a target area; the first seismic data and the second seismic data are acquired in different manners, and the resolution of the first seismic data is higher than that of the second seismic data.
Through the data acquisition system 120 shown in fig. 1, the first seismic data and the second seismic data of the target area may be acquired, where the acquisition modes for acquiring the two seismic data are different. And because the two types of seismic data are acquired in different modes, the two types of seismic data have different resolutions, and the resolution of the first seismic data is higher than that of the second seismic data.
Illustratively, 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 measured at different depths with receivers disposed in the well, while surface seismic data is seismic data measured with receivers disposed at the surface. The seismic data received by the wave detector in the well is more stable and accurate, the seismic data received by the wave detector in the ground surface is easily influenced by other signals and other interferents on the ground surface, so that the resolution of the VSP seismic data is higher than that of the ground seismic data, but the acquisition mode of the VSP seismic data is more difficult and higher in cost than that of the ground seismic data. Therefore, the ground seismic data with higher resolution can be obtained by shaping the ground seismic data with lower resolution by using VSP seismic data with higher resolution, and the method can be also suitable for other geological exploration areas.
Referring to the comparative cross-sectional view of VSP seismic data and ground seismic data shown in fig. 3, the abscissa is the ground distance, and the ordinate is the formation depth, which can be calculated from the time of receiving the seismic data and the speed of the seismic wave transmitted in the formation, so that the ordinate can also be regarded as the target time domain of the seismic data. From the VSP seismic data 301 and the ground seismic data 302 in FIG. 3, it can be seen that the sharpness of the VSP seismic data 301 is higher than the sharpness of the ground seismic data 302, i.e., the resolution of the VSP seismic data is higher than the resolution of the ground seismic data. Both VSP seismic data 301 and surface seismic data 302 are shown as a number of lines of varying depths, the depths of which may reflect the amount of signal energy received by the detectors. Whereas the change in the line may reflect the geological condition of the subsurface, such as whether rock is present, the layer is changed, etc.
Step 220, determining a matching factor set of the first seismic data and the second seismic data in the target frequency domain according to the first seismic data and the second seismic data.
The target frequency domain refers to a frequency domain range of the first seismic data and the second seismic data. Illustratively, the set of matching factors includes a plurality of matching factors, each matching factor corresponding to a frequency in the target frequency domain. A frequency corresponds to a first seismic data and a second seismic data, from which the shaped second seismic data can be derived.
In some embodiments, convolutions of the first seismic data and the reflection coefficients are calculated from the reflection coefficients of the first seismic data and the target region; calculating convolution of the second seismic data and the reflection coefficient according to the second seismic data and the reflection coefficient; the time-frequency domain data of the first seismic data and the time-frequency domain data of the second seismic data are obtained by respectively carrying out wavelet transformation on the convolution of the first seismic data and the reflection coefficient and the convolution of the second seismic data and the reflection coefficient; and determining a matching factor set of the first seismic data and the second seismic data in the target frequency domain according to the time-frequency domain data of the first seismic data and the time-frequency domain data of the second seismic data.
The surface-excited seismic waves travel perpendicularly from the medium into the formation interface, with some energy being reflected and some energy being transmitted through the interface, where the ratio of the reflected wave to the incident wave may be referred to as the reflection coefficient. The magnitude of the reflection coefficient is determined by the difference in wave impedance between the rock formations above and below the interface, and the larger the difference in wave impedance, the larger the reflection coefficient. Since both the first seismic data and the second seismic data are acquired from the target region, the reflection coefficients of the first seismic data and the second seismic data can be regarded as equal.
From the seismic data is a convolution of the seismic wavelet and the reflection coefficient, namely:
Xh(t)=Wh(t)*R(t)
Xl(t)=Wl(t)*R(t)
Wherein W h (t) is a first seismic wavelet, W l (t) is a second seismic wavelet, R (t) is a reflection coefficient of a target area, X h (t) is first seismic data, X l (t) is a convolution of the first seismic wavelet and the reflection coefficient, and X l (t) is second seismic data, which is a convolution of the second seismic wavelet and the reflection coefficient.
Performing wavelet transformation on X h (t) and X l (t) respectively to obtain time-frequency domain data of the first seismic data and time-frequency domain data of the second seismic data, namely:
Xh(f,t)=Wh(f,t)×R(f,t)
Xl(f,t)=Wl(f,t)×R(f,t)
Wherein, the value range of f is the target frequency domain, the value range of t is the target time domain, X h (f, t) is the time-frequency domain data of the first seismic data, and X l (f, t) is the time-frequency domain data of the second seismic data.
From the obtained X h (f, t) and X l (f, t), a set of matching factors for the first and second seismic data in the target frequency domain can be derived, namely:
H(f,t)=Xh(f,t)/Xl(f,t)=Wh(f,t)/Wl(f,t)
Wherein H (f, t) is a set of matching factors.
Through the steps, the 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 one frequency in the target frequency domain respectively, and the one-to-one corresponding matching factors can be used for obtaining the one-to-one corresponding shaped second seismic data.
Step 230, determining time-frequency domain shaping data of the second seismic data according to the second seismic data and the matching factor set, wherein the time-frequency domain shaping data is the seismic data with higher resolution than the second seismic data.
Illustratively, determining a target matching factor corresponding to the target frequency and target time-frequency domain data of the second seismic data according to the target frequency in the target frequency domain; and determining target time-frequency domain shaping data of the second seismic data according to the target matching factor and the target time-frequency domain data.
Applying the matching factor to the time-frequency domain data of the second seismic data can obtain time-frequency domain shaping data of the second seismic data, namely:
Xl,h(f,t)=Xl(f,t)×H(f,t)
Wherein X l,h (f, t) is the time-frequency domain shaping data of the second seismic data.
And respectively determining a matching factor corresponding to the frequency in the target frequency domain and the time-frequency domain shaping data, and obtaining the time-frequency domain shaping data of the second seismic data according to the formula. Illustratively, referring to the contrasted cross-sectional view of the time-frequency domain shaped data and the VSP seismic data shown in fig. 4, and the contrasted cross-sectional view of the time-frequency domain shaped data and the ground seismic data shown in fig. 5, it can be seen that the sharpness of the time-frequency domain shaped data is higher than the sharpness of the ground seismic data, and that the resolution of the instantaneous frequency domain shaped data is higher than the resolution of the ground seismic data. The comparison of fig. 5 shows that the continuity of the shaped data lines is better, the contact of the bent parts of the data lines is more obvious and clear, the continuity of the same phase axes of the shaped data can be reflected to be improved, the wave drag characteristic is also enhanced, and the overall quality of the shaped ground seismic data is greatly improved and promoted from the overall effect of fig. 4 and 5.
For example, referring to fig. 6, 7 and 8, the horizontal axis represents frequency and the vertical axis represents energy ratio of the seismic data compared to the originally excited seismic wave. The frequency spectrum characteristic diagram of the VSP seismic data shown in fig. 6 is wider in frequency band, richer in frequency components, and closer to zero phase than the frequency spectrum characteristic diagram of the ground seismic data shown in fig. 7. The spectrum characteristic diagram of the shaped seismic data is shown in fig. 8, and is very close to zero phase, and because the spectrum characteristic diagram is influenced by a plurality of factors in the geological exploration process, the obtained seismic data is almost impossible to completely zero phase, and the seismic data can be made to approach to zero phase through computer processing. In addition, the high-frequency information and the low-frequency information of the shaped seismic data in fig. 8 are rich, and the frequency bandwidth is wider than that of the ground seismic data shown in fig. 7. From the overall effect, the shaped seismic data effectively widens the frequency bandwidth of the seismic data, improves the resolution of the seismic data, ensures that the frequency components of the seismic data in the target frequency domain are more abundant, and is beneficial to the subsequent research using the shaped seismic data.
According to the technical scheme provided by the application, the matching factors of the first seismic data and the second seismic data are calculated, the first seismic data with higher resolution is used for adjusting the second seismic data with lower resolution, so that the adjusted time-frequency domain shaping data are obtained, and compared with the second seismic data before shaping, the shaped seismic data effectively widen the frequency bandwidth of the seismic data, and the resolution of the seismic data is obviously improved. And from the effect of the whole data section, the continuity of the same phase axis of the seismic data is improved, the wave drag characteristic is enhanced, and the whole quality of the section is greatly improved.
On the other hand, in the related art, when the ground seismic data is shaped, the seismic wavelets of the ground seismic data need to be extracted and deconvoluted, and the uncertain factors in the processing process are more, so that the processing process is more complex, and the quality control difficulty is higher. And because of the more influencing factors of wavelet morphology, it is difficult to accurately extract the seismic wavelet from the seismic data. And VSP seismic data is only used as priori information to compare and optimize ground seismic data, and is not involved in the processing operation process of seismic wavelets. According to the application, VSP seismic data is applied to the processing process of ground seismic data, the processing process of the ground seismic data is optimized while the seismic data with wider frequency bandwidth and higher resolution is obtained, the step of extracting seismic wavelets is omitted, and the scheme implementation is simpler.
In some embodiments, the time domain shaping data corresponding to the time domain shaping data is obtained by performing inverse transformation on the time domain shaping data; the time domain shaping data refers to seismic data of which the time domain shaping data is in a target time domain.
By performing inverse transformation on the X l,h (f, t), corresponding time domain shaping data, that is, a value range of the target time domain t, can be obtained. The formation depth of all seismic data received by the surface detectors can be determined by acquiring data of a target time domain and combining the propagation velocity of the surface excitation wave in the underground rock stratum. And then according to the waveform form, frequency spectrum characteristic, amplitude energy, phase and other characteristic aspects of the received seismic data, the method is beneficial to researching the geological form and the stratum condition of different strata, and the shaped seismic data is beneficial to researching the underground strata more finely and accurately.
Illustratively, as shown in the flowchart of fig. 9, the technical solution provided by the present application first obtains VSP seismic data with higher resolution and ground seismic data with lower resolution in the target area. And respectively carrying out data analysis on the VSP seismic data and the ground seismic data to obtain convolution of the seismic data and the reflection coefficient, and then carrying out wavelet transformation on the convolution to respectively obtain time-frequency domain data of the VSP seismic data and the ground seismic data in a target frequency domain. Thus, a set of matching factors of the VSP seismic data and the ground seismic data in the target frequency domain can be obtained according to the time-frequency domain data of the VSP seismic data and the time-frequency domain data of the ground seismic data, wherein each matching factor corresponds to one frequency in the target frequency domain.
The matching factor set in the target frequency domain is applied to the time-frequency domain data of the ground seismic data with lower resolution, so that the time-frequency domain shaping data with higher resolution in the target frequency domain can be obtained, the frequency bandwidth of the seismic data is widened by the time-frequency domain shaping data, and the resolution of the seismic data is improved. By inversely transforming the time-frequency domain shaping data, time-domain shaping data corresponding to the time-frequency domain shaping data can be obtained, and the analysis of the high-resolution time-domain shaping data is beneficial to researching the geological morphology and stratum characteristics of different stratum depths.
The following are examples of the apparatus of the present application that may be used to perform the method embodiments of the present application. For details not disclosed in the embodiments of the apparatus of the present application, please refer to the embodiments of the method of the present application.
Referring to fig. 10, a block diagram of a seismic data processing apparatus according to an embodiment of the application is shown. The device has the function of realizing the seismic data processing method, and the function can be realized by hardware or corresponding software executed by hardware. The apparatus may be the computer device described above or may be provided in a computer device. As shown in fig. 10, the apparatus 1000 may include: a first acquisition module 1010, a first determination module 1020, and a second determination module 1030.
A first acquisition module 1010 configured to acquire first seismic data and second seismic data of a target area; the first seismic data and the second seismic data are acquired in different manners, and the resolution of the first seismic data is higher than that of the second seismic data.
A first determining module 1020, configured to determine a set of matching factors of the first seismic data and the second seismic data in a target frequency domain according to the first seismic data and the second seismic data.
A second determining module 1030 is configured to determine time-frequency domain shaping data of the second seismic data according to the second seismic data and the set of matching factors, where the time-frequency domain shaping data is seismic data with a resolution higher than that of the second seismic data.
In some embodiments, the first determining module 1020 is configured to:
The time-frequency domain data of the first seismic data and the time-frequency domain data of the second seismic data are obtained by respectively carrying out wavelet transformation on the first seismic data and the second seismic data;
and determining a matching factor set of the first seismic data and the second seismic data in the target frequency domain according to the time-frequency domain data of the first seismic data and the time-frequency domain data of the second seismic data.
In some embodiments, the set of matching factors includes a plurality of matching factors, each matching factor corresponding to one frequency in the target frequency domain.
In some embodiments, the second determining module 1030 is configured to:
determining a target matching factor corresponding to the target frequency and target time-frequency domain data of the second seismic data according to the target frequency in the target frequency domain;
And determining target time-frequency domain shaping data of the second seismic data according to the target matching factor and the target time-frequency domain data.
In some embodiments, the apparatus 1000 further comprises:
The second obtaining module 1040 is configured to obtain time domain shaping data corresponding to the time domain shaping data by performing inverse transformation on the time domain shaping data; the time domain shaping data refers to seismic data of the time domain shaping data in a target time domain.
In some embodiments, the first seismic data is vertical seismic profile VSP seismic data and the second seismic data is ground seismic data.
According to the technical scheme provided by the application, the matching factors of the first seismic data and the second seismic data are calculated, the first seismic data with higher resolution is used for adjusting the second seismic data with lower resolution, so that the adjusted time-frequency domain shaping data are obtained, and compared with the second seismic data before shaping, the shaped seismic data effectively widen the frequency bandwidth of the seismic data, and the resolution of the seismic data is obviously improved. And from the effect of the whole data section, the continuity of the same phase axis of the seismic data is improved, the wave drag characteristic is enhanced, and the whole quality of the section is greatly improved.
It should be noted that, in the apparatus provided in the foregoing embodiment, when implementing the functions thereof, only the division of the foregoing functional modules is used as an example, in practical application, the foregoing functional allocation may be implemented by different functional modules, that is, the content structure of the device is divided into different functional modules, so as to implement all or part of the functions described above. In addition, the apparatus and the method embodiments provided in the foregoing embodiments belong to the same concept, and specific implementation processes of the apparatus and the method embodiments are detailed in the method embodiments and are not repeated herein.
Referring to fig. 11, a block diagram of a terminal device 1100 according to an embodiment of the present application is shown. The terminal device 1100 may be any electronic device having data computing, processing, and storage functions. The terminal apparatus 1100 may be used to implement the method of processing seismic data provided in the above-described embodiments.
In general, the terminal apparatus 1100 includes: a processor 1101 and a memory 1102.
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 ), PLA (Programmable Logic Array, programmable logic array). The processor 1101 may also include a main processor, which is a processor for processing data in an awake state, also called a CPU (Central Processing Unit ), and a coprocessor; a coprocessor is a low-power processor for processing data in a standby state. In some embodiments, the processor 1101 may integrate a GPU (Graphics Processing Unit, image processor) for rendering and drawing of content required to be displayed by the display screen. In some embodiments, the processor 1101 may also include an AI processor for processing computing operations related to machine learning.
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, as well as non-volatile memory, such as one or more magnetic disk storage devices, flash memory storage devices. In some embodiments, a non-transitory computer readable storage medium in memory 1102 is used to store a computer program configured to be executed by one or more processors to implement the above-described methods of processing seismic data.
It will be appreciated by those skilled in the art that the structure shown in fig. 11 is not limiting and that terminal device 1100 may include more or fewer components than shown, or may combine certain components, or may employ a different arrangement of components.
In an exemplary embodiment, a computer readable storage medium is also provided, in which a computer program is stored, which computer program, when being executed by a processor of a terminal device, implements the above-mentioned method of processing seismic data. Alternatively, the above-mentioned 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.
In an exemplary embodiment, a computer program product is also provided, the computer program product comprising a computer program stored in a computer readable storage medium. The processor of the terminal device reads the computer program from the computer-readable storage medium, and the processor executes the computer program so that the terminal device executes the above-described processing method of the seismic data.
It should be understood that references herein to "a plurality" are to two or more. "and/or", describes an association relationship of an association object, and indicates that there may be three relationships, for example, a and/or B, and may indicate: a exists alone, A and B exist together, and B exists alone. The character "/" generally indicates that the context-dependent object is an "or" relationship. In addition, the step numbers described herein are merely exemplary of one possible execution sequence among steps, and in some other embodiments, the steps may be executed out of the order of numbers, such as two differently numbered steps being executed simultaneously, or two differently numbered steps being executed in an order opposite to that shown, which is not limiting.
The foregoing description of the exemplary embodiments of the application is not intended to limit the application to the particular embodiments disclosed, but on the contrary, the intention is to cover all modifications, equivalents, and alternatives falling within the spirit and scope of the application.
Claims (10)
1. A method of processing seismic data, the method comprising:
Acquiring first seismic data and second seismic data of a target area; the first seismic data and the second seismic data are acquired in different manners, and the resolution of the first seismic data is higher than that of the second seismic data;
determining a matching factor set of the first seismic data and the second seismic data in a target frequency domain according to the first seismic data and the second seismic data;
And determining time-frequency domain shaping data of the second seismic data according to the second seismic data and the matching factor set, wherein the time-frequency domain shaping data is seismic data with higher resolution than the second seismic data.
2. The method of claim 1, wherein the determining a set of matching factors for the first seismic data and the second seismic data in a target frequency domain from the first seismic data and the second seismic data comprises:
The time-frequency domain data of the first seismic data and the time-frequency domain data of the second seismic data are obtained by respectively carrying out wavelet transformation on the first seismic data and the second seismic data;
and determining a matching factor set of the first seismic data and the second seismic data in the target frequency domain according to the time-frequency domain data of the first seismic data and the time-frequency domain data of the second seismic data.
3. The method of claim 1, wherein the set of matching factors includes a plurality of matching factors, each matching factor corresponding to a frequency in the target frequency domain.
4. The method of claim 1, wherein the determining time-frequency domain shaping data of the second seismic data from the second seismic data and the set of matching factors comprises:
determining a target matching factor corresponding to the target frequency and target time-frequency domain data of the second seismic data according to the target frequency in the target frequency domain;
And determining target time-frequency domain shaping data of the second seismic data according to the target matching factor and the target time-frequency domain data.
5. The method according to claim 1, wherein the method further comprises:
Obtaining time domain shaping data corresponding to the time domain shaping data by inversely transforming the time domain shaping data;
the time domain shaping data refers to seismic data of the time domain shaping data in a target time domain.
6. The method of claim 1, wherein the first seismic data is vertical seismic profile VSP seismic data and the second seismic data is ground seismic data.
7. A seismic data processing apparatus, the apparatus comprising:
the first acquisition module is used for acquiring first seismic data and second seismic data of the target area; the first seismic data and the second seismic data are acquired in different manners, and the resolution of the first seismic data is higher than that of the second seismic data;
The first determining module is used for determining a matching factor set of the first seismic data and the second seismic data in a target frequency domain according to the first seismic data and the second seismic data;
And the second determining module is used for determining time-frequency domain shaping data of the second seismic data according to the second seismic data and the matching factor set, wherein the time-frequency domain shaping data is seismic data with higher resolution than the second seismic data.
8. A computer device comprising a processor and a memory, the memory having stored therein a computer program that is loaded and executed by the processor to implement a method of processing seismic data as claimed in any one of claims 1 to 6.
9. A computer readable storage medium having stored therein a computer program that is loaded and executed by a processor to implement a method of processing seismic data as claimed in any one of claims 1 to 6.
10. A computer program product, characterized in that it comprises a computer program that is loaded and executed by a processor to implement a method of processing seismic data according to any of claims 1to 6.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211385799.1A CN117991345A (en) | 2022-11-07 | 2022-11-07 | Method, device, equipment, storage medium and program product for processing seismic data |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211385799.1A CN117991345A (en) | 2022-11-07 | 2022-11-07 | Method, device, equipment, storage medium and program product for processing seismic data |
Publications (1)
Publication Number | Publication Date |
---|---|
CN117991345A true CN117991345A (en) | 2024-05-07 |
Family
ID=90898098
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211385799.1A Pending CN117991345A (en) | 2022-11-07 | 2022-11-07 | Method, device, equipment, storage medium and program product for processing seismic data |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117991345A (en) |
-
2022
- 2022-11-07 CN CN202211385799.1A patent/CN117991345A/en active Pending
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11009617B2 (en) | Method for fast calculation of seismic attributes using artificial intelligence | |
US5673191A (en) | Method and apparatus for identifying geological structures using wavelet analysis of potential fields | |
CN107065013B (en) | A kind of interval velocity under earthquake scale determines method and device | |
US11561312B2 (en) | Mapping near-surface heterogeneities in a subterranean formation | |
US10983235B2 (en) | Characterizing depositional features by geologic-based seismic classification | |
CN109188520A (en) | Thin reservoir thickness prediction method and device | |
EA030770B1 (en) | System and method for seismic adaptive optics | |
CN109471162A (en) | Interbed multiple processing method, system, electronic equipment and readable medium | |
Lin et al. | Time-frequency mixed domain multi-trace simultaneous inversion method | |
CA2693247C (en) | Seismic attributes for structural analysis | |
CN111123361B (en) | Method and device for regularly reconstructing seismic data of vertical seismic profile and storage medium | |
US11740374B2 (en) | System and method for randomness measurement in sesimic image data using vectorized disorder algorithm | |
CN117991345A (en) | Method, device, equipment, storage medium and program product for processing seismic data | |
EP4384855A1 (en) | Method and system for seismic denoising using omnifocal reformation | |
US11867857B2 (en) | Method and system for updating a seismic velocity model | |
CN114755740A (en) | Rock distribution determination method, device, equipment and medium | |
CN113806674A (en) | Method and device for quantifying longitudinal dimension of ancient river channel, electronic equipment and storage medium | |
CN114002741B (en) | Pre-stack depth migration method and device and computer readable storage medium | |
CN113625338B (en) | Differential-based seismic weak signal extraction method | |
Devi et al. | High-resolution seismic signals from band-limited data using scaling laws of wavelet transforms | |
US20230266493A1 (en) | Seismic denoising by weighted multiplanar reconstruction | |
CN112083491B (en) | Wellbore wave pressing method and device based on signal attribute characteristics | |
CN113466936B (en) | Method, device, equipment and storage medium for acquiring CRP gather of fault shadow area | |
US20240210581A1 (en) | Implicit finite-difference method for wave equation based on recursive deconvolution | |
WO2023178632A1 (en) | Rock facies identification method based on seismic attribute classification using a machine learning network |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination |