Seismic data processing method and apparatus

Abstract

The present invention discloses a method and apparatus for seismic data processing. The method includes: acquiring shear wave data and longitudinal wave data in seismic data; performing pseudo-sampling processing on the shear wave data and the longitudinal wave data to obtain the shear wave data and the longitudinal wave data after pseudo-sampling processing; the pseudo-sampling processing is used to modify the sampling interval values in the trace headers and reel headers of the data header blocks of the shear wave data and the longitudinal wave data to preset values; interpreting the velocity spectra corresponding to the shear wave data and the longitudinal wave data after pseudo-sampling processing to form the correspondence between the time and velocity of the pseudo-sampled processing data; performing velocity analysis processing, residual static correction iterative processing, and time-domain imaging processing on the seismic data according to the correspondence between the time and velocity of the pseudo-sampled processing data. The present invention can change the data recording length through pseudo-sampling processing, can significantly improve the accuracy of seismic data velocity analysis; is beneficial to improving the effect of residual static correction iterative processing and enhancing the imaging quality of seismic data.

Description

Technical Field

[0001] The present invention relates to the technical field of seismic exploration data processing, and particularly to a seismic data processing method and apparatus. Background Art

[0002] This section aims to provide background or context for the embodiments of the present invention described in the claims. The description herein is not admitted to be prior art merely by virtue of its inclusion in this section.

[0003] Both P-waves and S-waves belong to elastic waves and can be used for seismic exploration. Since P-waves are easy to generate, have a single waveform, fast propagation speed, and are easy to interpret, the vast majority of previous seismic explorations have used P-waves. S-waves first entered the field of seismic exploration due to their low propagation speed. Passing through the same thickness of formation, S-waves take a longer time, and their ability to distinguish formations is generally higher than that of P-waves. Given that S-waves and P-waves have similar ray propagation paths and geometric seismological characteristics, the main processing steps of S-wave processing are basically the same as those of P-wave processing. Whether it is a P-wave or an S-wave, the original seismic record is obtained by discretely sampling the excited continuous seismic signal at a certain time interval. The product of the total number of samples in a seismic trace and the sampling interval is the record length; and the product of the serial number of a certain sample and the sampling interval is the time position of this sample in the seismic record, and this time position is fixed in the seismic record.

[0004] However, for the same set of underground formations, P-waves have a fast propagation speed and a short arrival time; S-waves have a slow propagation speed and a long time consumption; the time corresponding to the reflection wave group representing the formation in the P-wave data is different from the time corresponding to the reflection wave group in the S-wave data, and they form their own independent data characteristics. The current practice is to process P-wave and S-wave data in their respective spatio-temporal domains.

[0005] When processing S-wave data, the S-wave data faces processing problems such as low speed, low frequency, low signal-to-noise ratio, large static correction problems, and well-developed multiple waves, which make it difficult to accurately locate the S-wave velocity. The sensitivity of the S-wave velocity requires a higher velocity positioning accuracy; and this velocity analysis method of measuring the S-wave velocity change with a more concentrated velocity distribution range using the same velocity scale as P-waves and the spatio-temporal domain attribute characteristics of S-wave data fundamentally restrict the improvement of velocity interpretation accuracy.

[0006] When processing longitudinal wave data, various processing measures are taken around noise reduction and static correction methods of seismic data in order to obtain an accurate imaging velocity. Nowadays, most studies focus on solving the relevant factors affecting velocity quality, such as signal-to-noise ratio and static correction problems. After that, fine iteration of velocity and static correction is adopted to improve the accuracy of velocity analysis, and various velocity modeling methods are combined to finally improve the imaging quality. However, there is always an issue that cannot be ignored. Due to the relatively high velocity of longitudinal waves, the ability of longitudinal wave exploration to resolve strata is limited. In terms of velocity spectra, the sampling density in the time dimension is insufficient to depict the drastic changes in the vertical velocity of strata, which limits the potential of longitudinal wave data processing in the longitudinal wave domain.

[0007] In summary, when the current technology processes seismic data, the velocity spectra generated by longitudinal waves and shear waves in their respective spatio-temporal domains usually have objective problems such as low sampling density in the time or velocity dimension and poor recognition of velocity spectra. As a result, subsequent seismic data processing has problems such as poor effects in velocity analysis and residual static correction iterative processing and low seismic imaging quality. Summary of the Invention

[0008] An embodiment of the present invention provides a seismic data processing method for improving the accuracy of seismic data velocity analysis, improving the effect of residual static correction iterative processing, and improving the imaging quality of seismic data. The method includes:

[0009] Obtain shear wave data and longitudinal wave data in the seismic data;

[0010] Perform pseudo-sampling processing on the shear wave data and longitudinal wave data to obtain the shear wave data and longitudinal wave data after pseudo-sampling processing; the pseudo-sampling processing is used to modify the value corresponding to the sampling interval in the trace header and volume header of the data header block of the shear wave data and longitudinal wave data to a preset value; the shear wave data and longitudinal wave data after pseudo-sampling processing are used to represent the shear wave data and longitudinal wave data with the data record length changed after pseudo-sampling processing;

[0011] Interpret the velocity spectra corresponding to the shear wave data and longitudinal wave data after pseudo-sampling processing to form the corresponding relationship between time and velocity of the pseudo-sampled data;

[0012] According to the corresponding relationship between time and velocity of the pseudo-sampled data, perform velocity analysis processing, residual static correction iterative processing, and time-domain imaging processing of the seismic data.

[0013] An embodiment of the present invention also provides a seismic data processing device for improving the accuracy of seismic data velocity analysis, improving the effect of residual static correction iterative processing, and improving the imaging quality of seismic data. The device includes:

[0014] A data acquisition module for obtaining shear wave data and longitudinal wave data in the seismic data;

[0015] A pseudo-sampling processing module, configured to perform pseudo-sampling processing on shear wave data and longitudinal wave data to obtain the shear wave data and longitudinal wave data after pseudo-sampling processing; the pseudo-sampling processing is used to modify the value corresponding to the byte where the sampling interval is located in the trace header and volume header of the data header block of the shear wave data and longitudinal wave data to a preset value; the shear wave data and longitudinal wave data after pseudo-sampling processing are used to represent the shear wave data and longitudinal wave data whose data record lengths have changed after pseudo-sampling processing.

[0016] A velocity spectrum interpretation module, configured to interpret the velocity spectrum corresponding to the shear wave data and longitudinal wave data after pseudo-sampling processing to form the correspondence between the time and velocity of the pseudo-sampled data.

[0017] A data processing module, configured to perform velocity analysis processing, residual static correction iterative processing, and time-domain imaging processing on seismic data according to the correspondence between the time and velocity of the pseudo-sampled data.

[0018] An embodiment of the present invention further provides a computer device, including a memory, a processor, and a computer program stored on the memory and executable on the processor, where when the processor executes the computer program, the above-mentioned seismic data processing method is implemented.

[0019] An embodiment of the present invention further provides a computer-readable storage medium, where the computer-readable storage medium stores a computer program for executing the above-mentioned seismic data processing method.

[0020] In an embodiment of the present invention, shear wave data and longitudinal wave data in seismic data are acquired; the shear wave data and the longitudinal wave data are subjected to pseudo-sampling processing to obtain the shear wave data and the longitudinal wave data after pseudo-sampling processing; the pseudo-sampling processing is used to modify the sampling interval values in the trace headers and reel headers of the data header blocks of the shear wave data and the longitudinal wave data to preset values; the shear wave data and the longitudinal wave data after pseudo-sampling processing are used to represent the shear wave data and the longitudinal wave data with the data record length changed after pseudo-sampling processing; the velocity spectra corresponding to the shear wave data and the longitudinal wave data after pseudo-sampling processing are interpreted to form the correspondence between the time and velocity of the pseudo-sampled data; according to the correspondence between the time and velocity of the pseudo-sampled data, velocity analysis processing, residual static correction iterative processing, and time-domain imaging processing of the seismic data are performed. Compared with the prior art, the sampling interval values in the trace headers and reel headers of the data header blocks of the shear wave data and the longitudinal wave data can be modified through pseudo-sampling processing, thereby changing the data record lengths of the shear wave data and the longitudinal wave data, obtaining seismic data attribute characteristics different from the original spatio-temporal domain of the shear wave data and the longitudinal wave data, making the velocity distribution of the effective waves in the seismic data velocity spectrum generated therefrom more finely depicted and the distinguishability significantly enhanced. Therefore, the upward space for improving the accuracy of seismic data velocity analysis and interpretation can be expanded. Subsequently, through the correspondence between the time and velocity of the pseudo-sampled data, higher-efficiency velocity and residual static correction iteration of the seismic data and higher-quality time-domain imaging effects can be achieved. BRIEF DESCRIPTION OF THE DRAWINGS

[0021] In order to more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the following will briefly introduce the drawings required for use in the description of the embodiments or the prior art. Obviously, the following drawings are only some embodiments of the present invention. For those of ordinary skill in the art, other drawings can be obtained based on these drawings without creative efforts. In the drawings:

[0022] Figure 1 It is a flowchart of a seismic data processing method in an embodiment of the present invention;

[0023] Figure 2 It is a comparison diagram of the velocity spectra of the shear wave data before and after pseudo-sampling processing in an embodiment of the present invention;

[0024] Figure 3 It is a comparison diagram of the velocity spectra of the longitudinal wave data before and after pseudo-sampling processing in an embodiment of the present invention;

[0025] Figure 4 It is a comparison diagram of the imaging profiles of the shear wave data before and after pseudo-sampling processing in an embodiment of the present invention;

[0026] Figure 5 It is a comparison diagram of the imaging profiles of the longitudinal wave data before and after pseudo-sampling processing in an embodiment of the present invention;

[0027] Figure 6 It is a specific example diagram of a seismic data processing method in an embodiment of the present invention;

[0028] Figure 7 It is a schematic structural diagram of a seismic data processing device in an embodiment of the present invention. Specific implementation manners

[0029] To make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer and more understandable, the following further describes the embodiments of the present invention in detail with reference to the accompanying drawings. Herein, the illustrative embodiments of the present invention and their descriptions are used to explain the present invention, but not to limit the present invention.

[0030] Figure 1 It is a schematic flowchart of a seismic data processing method in an embodiment of the present invention. As Figure 1 shown, a seismic data processing method provided by an embodiment of the present invention may include:

[0031] Step 101: Obtain shear wave data and longitudinal wave data in seismic data;

[0032] Step 102: Perform pseudo-sampling processing on the shear wave data and the longitudinal wave data to obtain the shear wave data and the longitudinal wave data after pseudo-sampling processing; the pseudo-sampling processing is used to modify the value corresponding to the byte where the sampling interval is located in the trace header and the volume header of the data header block of the shear wave data and the longitudinal wave data to a preset value; the shear wave data and the longitudinal wave data after pseudo-sampling processing are used to represent the shear wave data and the longitudinal wave data whose data record length has changed after pseudo-sampling processing;

[0033] Step 103: Interpret the velocity spectra corresponding to the shear wave data and the longitudinal wave data after pseudo-sampling processing to form the corresponding relationship between the time and the velocity of the pseudo-sampled data;

[0034] Step 104: Perform velocity analysis processing, residual static correction iterative processing, and time-domain imaging processing on the seismic data according to the corresponding relationship between the time and the velocity of the pseudo-sampled data.

[0035] In an embodiment of the present invention, shear wave data and longitudinal wave data in seismic data are acquired; the shear wave data and the longitudinal wave data are subjected to pseudo-sampling processing to obtain the shear wave data and the longitudinal wave data after pseudo-sampling processing; the pseudo-sampling processing is used to modify the sampling interval values in the trace headers and volume headers of the data header blocks of the shear wave data and the longitudinal wave data to preset values; the shear wave data and the longitudinal wave data after pseudo-sampling processing are used to represent the shear wave data and the longitudinal wave data with the data record length changed after pseudo-sampling processing; the velocity spectra corresponding to the shear wave data and the longitudinal wave data after pseudo-sampling processing are interpreted to form the correspondence between the time and velocity of the pseudo-sampled data; according to the correspondence between the time and velocity of the pseudo-sampled data, velocity analysis processing, residual static correction iterative processing, and time-domain imaging processing of the seismic data are performed. Compared with the prior art, the sampling interval values in the trace headers and volume headers of the data header blocks of the shear wave data and the longitudinal wave data can be modified through pseudo-sampling processing, thereby changing the data record lengths of the shear wave data and the longitudinal wave data, obtaining seismic data attribute characteristics different from the original spatio-temporal domain of the shear wave data and the longitudinal wave data, making the velocity distribution of the effective waves in the seismic data velocity spectrum generated therefrom more finely depicted and the distinguishability significantly enhanced. Therefore, the upward space for improving the accuracy of seismic data velocity analysis and interpretation can be expanded. Subsequently, through the correspondence between the time and velocity of the pseudo-sampled data, higher-efficiency velocity and residual static correction iteration of the seismic data and higher-quality time-domain imaging effects can be achieved.

[0036] Specifically in implementation, first, shear wave data and longitudinal wave data in the seismic data are acquired.

[0037] In the embodiment, acquiring the shear wave data and the longitudinal wave data in the seismic data may include: acquiring the shear wave data and the longitudinal wave data in SEGY format from a seismic data processing system.

[0038] Specifically in implementation, a seismic data processing method provided by an embodiment of the present invention may further include: preprocessing the shear wave data and the longitudinal wave data in the seismic data; the preprocessing includes seismic data decompilation processing and definition processing of the seismic data acquisition system; performing pseudo-sampling processing on the shear wave data and the longitudinal wave data to obtain the shear wave data and the longitudinal wave data after pseudo-sampling processing may include: performing pseudo-sampling processing on the preprocessed shear wave data and longitudinal wave data to obtain the shear wave data and the longitudinal wave data after pseudo-sampling processing.

[0039] In the embodiment, before performing pseudo-sampling processing on the shear wave data and the longitudinal wave data, the preliminary processing of the seismic data mainly includes the preprocessing of the acquired seismic data, including:

[0040] In specific implementation, after obtaining the shear wave data and the longitudinal wave data in the seismic data, perform pseudo-sampling processing on the shear wave data and the longitudinal wave data to obtain the shear wave data and the longitudinal wave data after pseudo-sampling processing; the pseudo-sampling processing is used to modify the value corresponding to the byte where the sampling interval is located in the trace header and the volume header of the data header block of the shear wave data and the longitudinal wave data to a preset value; the shear wave data and the longitudinal wave data after pseudo-sampling processing are used to represent the shear wave data and the longitudinal wave data with the data record length changed after pseudo-sampling processing.

[0041] In the embodiment, perform pseudo-sampling processing on the shear wave data and the longitudinal wave data to obtain the shear wave data and the longitudinal wave data after pseudo-sampling processing, including: reducing or increasing the value corresponding to the byte where the sampling interval is located in the trace header and the volume header of the data header block of the shear wave data and the longitudinal wave data according to a preset multiple to obtain the shear wave data and the longitudinal wave data after pseudo-sampling processing.

[0042] In the above embodiment, by modifying the sampling interval value in the data volume header and the trace header, the pseudo-sampling processing of the shear wave data and the longitudinal wave data is realized. Applying the characteristics that the SEGY trace header and the data sample points are separable and combinable, the sampling interval value in the trace header and the volume header of the data header block of the seismic data is modified, but it will not affect the sample entity of the seismic data. Since the sampling interval value changes while the number of sample points remains unchanged, the data record length of the seismic record changes accordingly, and the corresponding compressed or stretched data record length is obtained. Therefore, the attribute characteristics of the seismic data will also change, such as the attribute characteristics of the velocity and time correspondence relationship of the seismic data.

[0043] Among them, the shear wave data and the longitudinal wave data after pseudo-sampling processing are used for but not limited to velocity analysis and residual static correction iteration to solve the problems of velocity, static correction and imaging of seismic data.

[0044] In specific implementation, perform pseudo-sampling processing on the shear wave data and the longitudinal wave data to obtain the shear wave data and the longitudinal wave data after pseudo-sampling processing, which may include: reducing or increasing the value corresponding to the byte where the sampling interval is located in the trace header and the volume header of the data header block of the shear wave data and the longitudinal wave data according to a preset multiple to obtain the shear wave data and the longitudinal wave data after pseudo-sampling processing.

[0045] In the embodiment, reducing the value corresponding to the byte where the sampling interval is located in the trace header and the volume header of the data header block of the shear wave data according to a preset multiple to obtain the shear wave data after pseudo-sampling processing may include: simultaneously modifying the value corresponding to the byte where the sampling interval is located in the trace header and the volume header of the binary header block of the shear wave data to half of the value corresponding to the byte where the original sampling interval is located to obtain the shear wave data after pseudo-sampling processing.

[0046] In an embodiment, increasing the value corresponding to the byte where the sampling interval is located in the trace header and the reel header of the longitudinal wave data header block by a preset multiple to obtain the longitudinally wave data after pseudo-sampling processing may include: simultaneously modifying the value corresponding to the byte where the sampling interval is located in the trace header and the reel header of the binary header block of the longitudinal wave data to twice the value corresponding to the byte where the original sampling interval is located, thereby obtaining the longitudinally wave data after pseudo-sampling processing.

[0047] In the above embodiment, by modifying the value corresponding to the byte where the sampling interval is located in the trace header and the reel header of the shear wave data and the longitudinal wave data header block, the modification of the seismic recording length of the shear wave data and the longitudinal wave data can be achieved. The processing advantages after the seismic data is compressed are mainly reflected in increasing the sampling density of the velocity dimension in the velocity spectrum. On the velocity spectrum based on the time-compressed data, the distribution range of the effective wave velocity can be increased by several times, the velocity difference between the primary reflection wave and the multiple interference waves can be increased, the distinguishability of the primary reflection energy cluster can be significantly enhanced, and the velocity analysis accuracy can be significantly improved. At the same time, the accurate primary wave velocity can make the velocity and residual static correction iteration converge faster, more effectively suppress the multiple waves, and improve the signal-to-noise ratio of the seismic data.

[0048] The processing advantages after the seismic data is stretched are mainly reflected in increasing the sampling density of the time dimension in the velocity spectrum. On the velocity spectrum based on the time-stretched data, the time span of a set of strata is increased by several times, effectively overcoming the disadvantage of insufficient time sampling of the longitudinal wave data in areas with drastic vertical velocity changes. The longitudinal fineness of velocity analysis and the data imaging quality are both significantly improved.

[0049] For example, Figure 2 as shown, Figure 2 is a comparison chart of the velocity spectra of the shear wave data before and after pseudo-sampling processing in the embodiment of the present invention. Figure 2 The area corresponding to the circular frame in [Figure] shows the different performances of the multiple waves on the velocity spectrum before and after pseudo-sampling processing. Figure 2 The area corresponding to the square frame in [Figure] shows the change in the energy cluster distribution characteristics of the shallow multiple waves and the primary wave before and after time compression. Figure 2 It shows that the multiple waves and the primary wave on the velocity spectrum after pseudo-sampling are significantly separated compared to before. Figure 2 The ordinate in [Figure] represents time, and the abscissa represents velocity. For example, Figure 3 as shown, Figure 3 is a comparison chart of the velocity spectra of the longitudinal wave data before and after pseudo-sampling processing in the embodiment of the present invention. Figure 3 The ordinate in [Figure] represents time, and the abscissa represents velocity. Figure 3 It shows that after the time of the longitudinal wave data is stretched, the sampling density in the time direction of the velocity spectrum increases, and the vertical change trend of the velocity is more finely characterized.

[0050] Another example, Figure 4This is a comparison diagram of the imaging profiles before and after pseudo-sampling processing of shear wave data in the embodiments of the present invention. Figure 4 The ordinate represents time, and the abscissa represents the horizontal arrangement of CMP points. From Figure 4 it can be seen that both the multiple wave problem and the static correction problem in the imaging profile of the shear wave data after being processed by the pseudo-sampling technique have been significantly improved. Figure 5 This is a comparison diagram of the imaging profiles before and after pseudo-sampling processing of compressional wave data in the embodiments of the present invention. Figure 5 The ordinate represents time, and the abscissa represents the horizontal arrangement of CMP points. Figure 5 It shows that after the compressional wave data is processed by the pseudo-sampling technique, the faults, diffractions, and interlayer contact relationships on the stacked profile are more clearly presented.

[0051] In a specific embodiment, let the arrival time of the wave at any point underground in the seismic data be t n , Δt be the sampling interval, and n be the sample point sequence number. Then there is: t n = Δt × n; it can be seen from the above formula that as long as the sampling interval value is changed, the arrival time of the wave can be changed; the sampling interval can be made larger or smaller;

[0052] For example, if the sampling interval value is reduced by 1 time, that is, let Δt' = Δt / 2, where Δt' is the changed sampling interval value, then there is: t n ' = Δt' × n = (Δt / 2) × n = (Δt × n) / 2 = t n / 2, where t n ' is the arrival time of the wave after changing the sampling interval value. That is, after the sampling interval value is reduced by 1 time, the arrival time of the wave is compressed by 1 time. Similarly, after the sampling interval value is increased by 1 time, it can be deduced that the arrival time of the wave is increased by 1 time, and the seismic record is stretched by 1 time.

[0053] Suppose the acquisition record length of the shear wave data or compressional wave data in the original seismic data (i.e., the seismic record length) is 10 s, and the sampling interval is 2 ms. Then each record contains 5000 sample points; if the sampling interval data is changed to 1 ms, then the total number of sample points remains unchanged, still 5000, and the seismic record length becomes 5 s; if the sampling interval is changed to 4 ms, then there are still these 5000 sample points, and the record length becomes 20 s.

[0054] As described above, before and after the time of the seismic record is stretched or compressed, the depth where the wave group is located remains unchanged; when the time changes, the velocity changes accordingly. The velocity, depth, and time follow the following corresponding relationship:

[0055] H n = v n × t n ; where H nis depth, t n is time, v n is velocity.

[0056] For example: when the time is compressed by a factor of 1, H n =(t n / 2)×(2v n ); when the time is stretched by a factor of 1, then H n =(2×t n )×(v n / 2).

[0057] In summary, the embodiments of the present invention can change the correspondence between the velocity and time of seismic data, which is also a means to obtain key imaging parameters. By applying the above correspondence, P-wave data can be switched to the spatio-temporal domain processing of S-waves; S-wave data can also be switched to the P-wave spatio-temporal domain processing. The purpose of doing this is to make full use of the seismic attribute advantages of P-wave and S-wave data, avoid their respective disadvantages, and thus expand the improvement space of velocity analysis accuracy and imaging quality.

[0058] Specifically in implementation, a seismic data processing method provided by the embodiments of the present invention may further include: directly performing denoising, velocity analysis, iteration of residual static correction, and time-domain imaging processing on the S-wave data and P-wave data after pseudo-sampling processing.

[0059] In the embodiment, by performing pseudo-sampling processing on the S-wave data and P-wave data, the reconstruction of the S-wave data and P-wave data is realized, and the reconstruction of the seismic record length of the S-wave data and P-wave data is realized.

[0060] In the above embodiment, the reconstructed S-wave data and P-wave data can be directly processed for imaging.

[0061] Specifically in implementation, after interpreting the velocity spectra corresponding to the S-wave data and P-wave data after pseudo-sampling processing to form the correspondence between the time and velocity of the pseudo-sampled processed data, according to the correspondence between the time and velocity of the pseudo-sampled processed data, perform velocity analysis processing, residual static correction iterative processing, and time-domain imaging processing of seismic data.

[0062] In the embodiment, according to the correspondence between the time and velocity of the pseudo-sampled processed data, performing velocity analysis processing, residual static correction iterative processing, and time-domain imaging processing of seismic data may include: performing conversion on the correspondence between the time and velocity of the pseudo-sampled processed data; according to the correspondence between the time and velocity of the pseudo-sampled processed data after conversion, perform iteration of velocity analysis and residual static correction and time-domain imaging processing on seismic data.

[0063] In the above embodiments, the conversion may include: reducing the speed in the correspondence between the time and speed of the pseudo-sampled processed data by a first preset multiple and increasing the time by a second preset multiple; increasing the speed in the correspondence between the time and speed of the pseudo-sampled processed data by a third preset multiple and reducing the time by a fourth preset multiple, etc.

[0064] In the embodiments, according to the correspondence between the time and speed of the pseudo-sampled processed data, performing velocity analysis processing, residual static correction iterative processing, and time-domain imaging processing on seismic data may include: when performing denoising stacking imaging processing and residual static correction iterative processing on shear wave data in seismic data, reducing the speed in the correspondence between the time and speed of the pseudo-sampled processed data by a first preset multiple and increasing the time by a second preset multiple to obtain the correspondence between the time and speed of the original shear wave data in the corresponding seismic data; according to the correspondence between the time and speed of the original shear wave data in the corresponding seismic data, performing velocity analysis processing, residual static correction iterative processing, and time-domain imaging processing on the original shear wave data in the seismic data.

[0065] In the embodiments, according to the correspondence between the time and speed of the pseudo-sampled processed data, performing velocity analysis processing, residual static correction iterative processing, and time-domain imaging processing on seismic data may include: when performing velocity analysis processing, residual static correction iterative processing, and time-domain imaging processing on seismic data, increasing the speed in the correspondence between the time and speed of the pseudo-sampled processed data by a third preset multiple and reducing the time by a fourth preset multiple to obtain the correspondence between the time and speed of the original compressional wave data in the corresponding seismic data; according to the correspondence between the time and speed of the original compressional wave data in the corresponding seismic data, performing velocity analysis processing, residual static correction iterative processing, and time-domain imaging processing on the seismic data.

[0066] In the above embodiments, by performing pseudo-sampling processing on shear wave data and compressional wave data, the reconstruction of shear wave data and compressional wave data is realized. By interpreting the velocity spectra of the reconstructed shear wave data and compressional wave data, the correspondence between the time and speed of the original shear wave data in the corresponding seismic data can be formed, and imaging processing can be performed on the original seismic data through the converted correspondence between the time and speed.

[0067] For a specific example, such as Figure 6 shown Figure 6 is a specific example diagram of a seismic data processing method in an embodiment of the present invention. This example may specifically include the following steps:

[0068] 1. Prepare shear wave data and compressional wave data in seismic data to be input into the processing system

[0069] Input the original seismic data into the processing system and perform preprocessing on the seismic data input into the processing system;

[0070] 2. Perform SEGY format conversion and output, and at the same time modify the corresponding trace headers and volume labels in the data header block to implement the present invention.

[0071] Convert the preprocessed seismic data input in step 1 into SEGY format data for output, and increase or decrease the value corresponding to the byte where the sampling interval is located in the binary header block of the SEGY format data by a preset multiple. This step does not perform true resampling on the seismic data, but pseudo-sampling, where the number of output samples remains unchanged;

[0072] 3. Obtain the reconstructed seismic data

[0073] Re-introduce the SEGY format data with the header block content reset into the processing system. In this way, the record lengths of the shear wave data and the compressional wave data after pseudo-sampling processing introduced have been compressed or stretched;

[0074] 4. Application of the reconstructed data

[0075] The reconstructed data can be directly used for velocity analysis or can be used for data processing such as denoising and stacking imaging after re-positioning the survey.

[0076] As described above, Figure 6 the core steps of the embodiment of the present invention in the dashed box. It should be noted that the embodiment of the present invention only modifies the sampling interval value, rather than performing resampling processing on the data.

[0077] A specific embodiment is given below to illustrate the specific application of the method of the present invention. In this embodiment, the following steps may be included:

[0078] 1. Preprocessing of shear wave data in seismic data

[0079] 2. Perform pseudo-sampling processing on the preprocessed shear wave data

[0080] Convert the data output in step 1 into SEGY format data, and define the value at the byte where the sampling interval is located in the binary header block of the SEGY format data as half of the original sampling interval. At the same time as modifying the trace header, the volume header must also be modified;

[0081] 3. Obtain the pseudo-sampling processed data of the shear wave record

[0082] Re-introduce the SEGY format data with the sampling interval reset into the processing system. In this way, the introduced data has achieved pseudo-sampling of the shear wave record. Since the recorded sample values remain unchanged, the data record length has been reduced by 1 time;

[0083] 4. Application of the pseudo-sampling data

[0084] The pseudo-sampled data is mainly used for velocity analysis and residual static correction iteration to obtain the velocity and static correction parameter required for imaging. After redefining the acquisition system for the pseudo-sampled data output in the third step and performing a velocity spectrum, a corresponding relationship file of time and velocity based on the pseudo-sampled data (i.e., the corresponding relationship between the time and velocity of the pseudo-sampled processed data) is formed through velocity interpretation. Dividing the velocity in this file by 2 and multiplying the time by 2, a velocity-time corresponding relationship file corresponding to the original data (i.e., the corresponding relationship between the time and velocity of the original shear wave data in the seismic data) is obtained, which can be used for the residual static correction iteration and imaging processing of the original shear wave data.

[0085] The above specific embodiment is derived from the acquisition experiment of 2D wide-line seismic data excited by a shear wave vibrator in a certain area of the Junggar Basin. In this example, by applying the seismic data processing method provided by the embodiment of the present invention, for the velocity spectrum made from the shear wave data after pseudo-sampling processing, the recognizability of the primary wave is significantly enhanced, and the energy clusters of the primary wave and the multiple waves become easy to distinguish (as Figure 2 shown), the accuracy of velocity analysis and the velocity iteration efficiency are significantly improved, and the signal-to-noise ratio and resolution of the resulting data are both significantly improved (as Figure 4 shown); the processing result quality of the shear wave data is directly related to the evaluation of the exploration value of the shear wave. The embodiment of the present invention can provide an effective processing means for the popularization and application of the shear wave exploration method in the field of seismic exploration.

[0086] Another specific embodiment is given below to illustrate the specific application of the method of the present invention. In this embodiment, the following steps may be included:

[0087] 1. Preprocessing of P-wave seismic data

[0088] 2. Performing pseudo-sampling processing on the preprocessed P-wave data

[0089] Convert the data output in step 1 into SEGY format data, and define the value at the byte where the sampling interval is located in the binary header block of the SEGY format data as 2 times the original sampling interval. When modifying the trace header, the volume header must also be modified;

[0090] 3. Obtaining the time-stretched data of the P-wave record

[0091] Reintroduce the SEGY format data with the reset sampling interval into the processing system. The data introduced in this way has achieved pseudo-sampling of the P-wave record. Since the recorded sample point values remain unchanged, the data record length is stretched by 1 time;

[0092] 4. Application of the pseudo-sampled data

[0093] The pseudo-sampled data is used for velocity analysis and stacking imaging. After redefining the acquisition system for the pseudo-sampled data output in the third step, a velocity spectrum is made, and through velocity interpretation, a corresponding relationship file of time and velocity based on the pseudo-sampled data is formed. This file can be directly used for the imaging processing of the P-wave stretched data. It is also possible to multiply the velocity in this file by 2 and divide the time by 2 to obtain a corresponding relationship file of velocity and time corresponding to the original data (i.e., the corresponding relationship of time and velocity for the original P-wave data in the seismic data), which can be used for the imaging of the original P-wave data.

[0094] The above specific embodiment was carried out in a certain area of the Junggar Basin. During the specific implementation, two cross-line 2D P-wave survey lines were collected for processing. The target concerned in the embodiment of the present invention is in the middle and shallow layers, where the formation structure is complex and the velocity changes violently. Applying the pseudo-sampling processing technology can stretch and switch the P-wave data to the S-wave spatio-temporal domain for processing, significantly improving the precision of velocity analysis in the time direction (as Figure 3 shown), making the stacked section obtained present various wave fields more clearly and richly (as Figure 6 shown). The idea of applying the pseudo-sampling technology to process P-wave data in the S-wave domain has significant application effects on mountain surface layer velocity modeling, extended application of P-wave micro-log data, etc.

[0095] In the embodiment of the present invention, S-wave data and P-wave data in seismic data are acquired; pseudo-sampling processing is performed on the S-wave data and P-wave data to obtain the S-wave data and P-wave data after pseudo-sampling processing; the pseudo-sampling processing is used to modify the sampling interval values in the trace header and reel header of the data header block of the S-wave data and P-wave data to a preset value; the S-wave data and P-wave data after pseudo-sampling processing are used to represent the S-wave data and P-wave data with the data record length changed after pseudo-sampling processing; the velocity spectra corresponding to the S-wave data and P-wave data after pseudo-sampling processing are interpreted to form the corresponding relationship of time and velocity of the pseudo-sampled data; according to the corresponding relationship of time and velocity of the pseudo-sampled data, velocity analysis processing, residual static correction iterative processing and time-domain imaging processing of seismic data are carried out. Compared with the prior art, by pseudo-sampling processing, the sampling interval values in the trace header and reel header of the data header block of the S-wave data and P-wave data can be modified, thereby changing the data record length of the S-wave data and P-wave data, obtaining seismic data attribute characteristics different from the original spatio-temporal domain of the S-wave data and P-wave data, making the velocity distribution of the effective wave in the seismic data velocity spectrum depicted more finely and the recognizability significantly enhanced. Therefore, the upward space for improving the accuracy of seismic data velocity analysis and interpretation can be expanded. Subsequently, through the corresponding relationship of time and velocity of the pseudo-sampled data, higher-efficiency velocity and residual static correction iteration and higher-quality time-domain imaging effects of seismic data can be achieved.

[0096] As described above, the embodiment of the present invention is a method of false transformation processing in the spatio-temporal domain, which can significantly improve the precision of velocity analysis of seismic data of P-waves and S-waves. It aims to flexibly transform the required spatio-temporal domain in the processing of P-wave and S-wave seismic data, fundamentally solve the problems such as low sampling density in the time or velocity dimension and low velocity analysis accuracy of the velocity spectra of P-wave and S-wave data, significantly increase the recognizability of effective waves in the velocity spectrum and the fineness of the characterization of effective reflection energy, and thus improve the imaging effect of seismic data. This method does not change the seismic data samples themselves and can be flexibly applied.

[0097] In an embodiment of the present invention, a seismic data processing device is also provided, as described in the following embodiments. Since the principle of the device for solving problems is similar to that of the seismic data processing method, the implementation of the device can refer to the implementation of the seismic data processing method, and the repeated parts will not be described again.

[0098] Figure 7 It is a schematic structural diagram of a seismic data processing device in an embodiment of the present invention, as Figure 7 shown, the seismic data processing device provided by the embodiment of the present invention may include:

[0099] A data acquisition module 01, configured to acquire S-wave data and P-wave data in seismic data;

[0100] A pseudo-sampling processing module 02, configured to perform pseudo-sampling processing on the S-wave data and P-wave data to obtain the S-wave data and P-wave data after pseudo-sampling processing; the pseudo-sampling processing is used to modify the value corresponding to the sampling interval in the trace header and volume header of the data header block of the S-wave data and P-wave data to a preset value; the S-wave data and P-wave data after pseudo-sampling processing are used to represent the S-wave data and P-wave data with the data record length changed after pseudo-sampling processing;

[0101] A velocity spectrum interpretation module 03, configured to interpret the velocity spectra corresponding to the S-wave data and P-wave data after pseudo-sampling processing to form the correspondence between the time and velocity of the pseudo-sampled data;

[0102] A data processing module 04, configured to perform velocity analysis processing, residual static correction iterative processing, and time-domain imaging processing on the seismic data according to the correspondence between the time and velocity of the pseudo-sampled data.

[0103] In one embodiment, the data acquisition module is specifically configured to: acquire S-wave data and P-wave data in SEGY format from a seismic data processing system.

[0104] In one embodiment, it further includes: a preprocessing module, configured to: perform preprocessing on the S-wave data and P-wave data in the seismic data; the preprocessing includes seismic data decompilation processing and definition processing of the seismic data acquisition system;

[0105] The pseudo-sampling processing module is specifically configured to: perform pseudo-sampling processing on the preprocessed shear wave data and longitudinal wave data to obtain the pseudo-sampled shear wave data and longitudinal wave data.

[0106] In one embodiment, the pseudo-sampling processing module is specifically configured to: decrease or increase the value corresponding to the byte where the sampling interval is located in the trace header and reel header of the data header block of the shear wave data and longitudinal wave data by a preset multiple to obtain the pseudo-sampled shear wave data and longitudinal wave data.

[0107] In one embodiment, the pseudo-sampling processing module is specifically configured to: simultaneously modify the value corresponding to the byte where the sampling interval is located in the trace header and reel header of the binary header block of the shear wave data to half of the value corresponding to the byte where the original sampling interval is located to obtain the pseudo-sampled shear wave data.

[0108] In one embodiment, the pseudo-sampling processing module is specifically configured to: simultaneously modify the value corresponding to the byte where the sampling interval is located in the trace header and reel header of the binary header block of the longitudinal wave data to twice the value corresponding to the byte where the original sampling interval is located to obtain the pseudo-sampled longitudinal wave data.

[0109] In one embodiment, it further includes: a processing and imaging module, configured to: directly perform denoising, velocity analysis, iteration of residual static correction, and time-domain imaging processing on the pseudo-sampled shear wave data and longitudinal wave data.

[0110] In one embodiment, the data processing module is specifically configured to:

[0111] When performing denoising, stacking, and imaging processing and iteration of residual static correction on the shear wave data in seismic data, decrease the velocity in the corresponding relationship between the time and velocity of the pseudo-sampled data by a first preset multiple and increase the time by a second preset multiple to obtain the corresponding relationship between the time and velocity of the original shear wave data in the corresponding seismic data;

[0112] Perform denoising, stacking, and imaging processing and iteration of residual static correction on the original shear wave data in the seismic data according to the corresponding relationship between the time and velocity of the original shear wave data in the corresponding seismic data.

[0113] In one embodiment, the data processing module is specifically configured to:

[0114] When performing velocity analysis processing, iteration of residual static correction, and time-domain imaging processing on the seismic data, increase the velocity in the corresponding relationship between the time and velocity of the pseudo-sampled data by a third preset multiple and decrease the time by a fourth preset multiple to obtain the corresponding relationship between the time and velocity of the original longitudinal wave data in the corresponding seismic data;

[0115] According to the corresponding relationship between the time and velocity of the original P-wave data in the seismic data, perform velocity analysis processing, residual static correction iterative processing, and time-domain imaging processing on the seismic data.

[0116] An embodiment of the present invention also provides a computer device, including a memory, a processor, and a computer program stored on the memory and executable on the processor. When the processor executes the computer program, the above seismic data processing method is implemented.

[0117] An embodiment of the present invention also provides a computer-readable storage medium, which stores a computer program for executing the above seismic data processing method.

[0118] In an embodiment of the present invention, S-wave data and P-wave data in seismic data are obtained; pseudo-sampling processing is performed on the S-wave data and P-wave data to obtain the S-wave data and P-wave data after pseudo-sampling processing; the pseudo-sampling processing is used to modify the sampling interval values in the trace headers and volume headers of the data header blocks of the S-wave data and P-wave data to a preset value; the S-wave data and P-wave data after pseudo-sampling processing are used to represent the S-wave data and P-wave data with the data record length changed after pseudo-sampling processing; the velocity spectra corresponding to the S-wave data and P-wave data after pseudo-sampling processing are interpreted to form the corresponding relationship between the time and velocity of the pseudo-sampled data; according to the corresponding relationship between the time and velocity of the pseudo-sampled data, perform velocity analysis processing, residual static correction iterative processing, and time-domain imaging processing on the seismic data. Compared with the prior art, the sampling interval values in the trace headers and volume headers of the data header blocks of the S-wave data and P-wave data can be modified through pseudo-sampling processing, thereby changing the data record lengths of the S-wave data and P-wave data, obtaining seismic data attribute characteristics different from the original spatio-temporal domain of the S-wave data and P-wave data, making the velocity distribution of the effective waves in the seismic data velocity spectrum depicted more finely and the distinguishability significantly enhanced. Therefore, the upward space for improving the accuracy of seismic data velocity analysis and interpretation can be expanded. Subsequently, through the corresponding relationship between the time and velocity of the pseudo-sampled data, higher-efficiency velocity and residual static correction iteration of seismic data and higher-quality time-domain imaging effects can be achieved.

[0119] 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 storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

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

[0121] 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, such that the instructions stored in the computer-readable memory generate a manufactured article including instruction means that implement the functions specified in one or more of the flows Figure 1 one or more of the flows and / or blocks Figure 1 or means for implementing the functions specified in one or more of the blocks.

[0122] These computer program instructions can also be loaded onto a computer or other programmable data processing device, such that a series of operational steps are performed on the computer or other programmable device to generate a computer-implemented process, and thus the instructions executed on the computer or other programmable device provide steps for implementing the functions specified in one or more of the flows Figure 1 one or more of the flows and / or blocks Figure 1 or means for implementing the functions specified in one or more of the blocks.

[0123] The specific embodiments described above further elaborate on the objectives, technical solutions, and beneficial effects of the present invention. It should be understood that the above are only specific embodiments of the present invention and are not used to limit the protection scope of the present invention. Any modifications, equivalent replacements, improvements, etc. made within the spirit and principles of the present invention shall be included in the protection scope of the present invention.

Claims

1. A method for processing seismic data, characterized in that, Including: Obtaining shear wave data and longitudinal wave data in seismic data; Performing pseudo-sampling processing on the shear wave data and longitudinal wave data to obtain the shear wave data and longitudinal wave data after pseudo-sampling processing; the pseudo-sampling processing is used to modify the value corresponding to the byte where the sampling interval is located in the trace header and reel header of the data header block of the shear wave data and longitudinal wave data to a preset value; the shear wave data and longitudinal wave data after pseudo-sampling processing are used to represent the shear wave data and longitudinal wave data with the data record length changed after pseudo-sampling processing; Interpret the velocity spectra corresponding to the shear wave data and longitudinal wave data after pseudo-sampling processing to form the corresponding relationship between the time and velocity of the pseudo-sampled data; Perform velocity analysis processing, residual static correction iterative processing, and time-domain imaging processing on the seismic data according to the corresponding relationship between the time and velocity of the pseudo-sampled data.

2. The method according to claim 1, characterized in that Obtaining shear wave data and longitudinal wave data in seismic data includes: obtaining shear wave data and longitudinal wave data in SEGY format from a seismic data processing system.

3. The method according to claim 1, wherein It also includes: Performing preprocessing on the shear wave data and longitudinal wave data in the seismic data; the preprocessing includes seismic data decompression processing and definition processing of the seismic data acquisition system; Performing pseudo-sampling processing on the shear wave data and longitudinal wave data to obtain the shear wave data and longitudinal wave data after pseudo-sampling processing, including: Performing pseudo-sampling processing on the shear wave data and longitudinal wave data after preprocessing to obtain the shear wave data and longitudinal wave data after pseudo-sampling processing.

4. The method according to claim 1, wherein Performing pseudo-sampling processing on the shear wave data and longitudinal wave data to obtain the shear wave data and longitudinal wave data after pseudo-sampling processing, including: Reducing or increasing the value corresponding to the byte where the sampling interval is located in the trace header and reel header of the data header block of the shear wave data and longitudinal wave data by a preset multiple to obtain the shear wave data and longitudinal wave data after pseudo-sampling processing.

5. The method according to claim 4, wherein Reducing the value corresponding to the byte where the sampling interval is located in the trace header and reel header of the data header block of the shear wave data by a preset multiple to obtain the shear wave data after pseudo-sampling processing, including: Simultaneously modifying the value corresponding to the byte where the sampling interval is located in the trace header and reel header of the binary header block of the shear wave data to half of the value corresponding to the byte where the original sampling interval is located to obtain the shear wave data after pseudo-sampling processing.

6. The method according to claim 4, wherein Increasing the value corresponding to the byte where the sampling interval is located in the trace header and reel header of the data header block of the longitudinal wave data by a preset multiple to obtain the longitudinal wave data after pseudo-sampling processing, including: Simultaneously modifying the value corresponding to the byte where the sampling interval is located in the trace header and reel header of the binary header block of the longitudinal wave data to twice the value corresponding to the byte where the original sampling interval is located to obtain the longitudinal wave data after pseudo-sampling processing.

7. The method according to claim 1, wherein It also includes: Directly performing denoising, velocity analysis, iteration of residual static correction, and time-domain imaging processing on the shear wave data and longitudinal wave data after pseudo-sampling processing.

8. The method according to claim 1, wherein Performing velocity analysis processing, residual static correction iterative processing, and time-domain imaging processing on the seismic data according to the corresponding relationship between the time and velocity of the pseudo-sampled data, including: When performing velocity analysis processing, residual static correction iterative processing, and time-domain imaging processing on seismic data, reduce the velocity in the correspondence between the time and velocity of the pseudo-sampled processing data by a first preset multiple and increase the time by a second preset multiple to obtain the correspondence between the time and velocity of the original shear wave data in the corresponding seismic data; Perform velocity analysis processing, residual static correction iterative processing, and time-domain imaging processing on the seismic data according to the correspondence between the time and velocity of the original shear wave data in the corresponding seismic data.

9. The method according to claim 1, wherein Performing velocity analysis processing, residual static correction iterative processing, and time-domain imaging processing on the seismic data according to the correspondence between the time and velocity of the pseudo-sampled processing data includes: When performing velocity analysis processing, residual static correction iterative processing, and time-domain imaging processing on seismic data, increase the velocity in the correspondence between the time and velocity of the pseudo-sampled processing data by a third preset multiple and decrease the time by a fourth preset multiple to obtain the correspondence between the time and velocity of the original compressional wave data in the corresponding seismic data; Perform velocity analysis processing, residual static correction iterative processing, and time-domain imaging processing on the seismic data according to the correspondence between the time and velocity of the original compressional wave data in the corresponding seismic data.

10. An earthquake data processing device, characterized in that, Including: A data acquisition module for acquiring shear wave data and compressional wave data in seismic data; A pseudo-sampling processing module for performing pseudo-sampling processing on the shear wave data and compressional wave data to obtain the pseudo-sampled shear wave data and compressional wave data; the pseudo-sampling processing is used to modify the value corresponding to the byte where the sampling interval is located in the trace header and reel header of the data header block of the shear wave data and compressional wave data to a preset value; the pseudo-sampled shear wave data and compressional wave data are used to represent the shear wave data and compressional wave data with the data record length changed after pseudo-sampling processing; A velocity spectrum interpretation module for interpreting the velocity spectra corresponding to the pseudo-sampled shear wave data and compressional wave data to form the correspondence between the time and velocity of the pseudo-sampled processing data; A data processing module for performing velocity analysis processing, residual static correction iterative processing, and time-domain imaging processing on the seismic data according to the correspondence between the time and velocity of the pseudo-sampled processing data.

11. The device according to claim 10, characterized in that, The data acquisition module is specifically used for: acquiring shear wave data and compressional wave data in SEGY format from a seismic data processing system.

12. The device according to claim 10, wherein It also includes: A preprocessing module for: preprocessing the shear wave data and compressional wave data in the seismic data; the preprocessing includes seismic data decoding processing and definition processing of the seismic data acquisition system; The pseudo-sampling processing module is specifically used for: performing pseudo-sampling processing on the preprocessed shear wave data and compressional wave data to obtain the pseudo-sampled shear wave data and compressional wave data.

13. The device according to claim 10, characterized in that, The pseudo-sampling processing module is specifically used for: reducing or increasing the value corresponding to the byte where the sampling interval is located in the trace header and reel header of the shear wave data and compressional wave data by a preset multiple to obtain the pseudo-sampled shear wave data and compressional wave data.

14. The device according to claim 13, characterized in that, The pseudo-sampling processing module is specifically configured to: simultaneously modify the values corresponding to the bytes where the sampling intervals are located in the trace header and the volume header of the binary header block of the shear wave data to half of the values corresponding to the bytes where the original sampling intervals are located, so as to obtain the shear wave data after pseudo-sampling processing.

15. The device according to claim 13, characterized in that, The pseudo-sampling processing module is specifically configured to: simultaneously modify the values corresponding to the bytes where the sampling intervals are located in the trace header and the volume header of the binary header block of the compressional wave data to twice the values corresponding to the bytes where the original sampling intervals are located, so as to obtain the compressional wave data after pseudo-sampling processing.

16. The device according to claim 10, wherein, It further includes: The processing and imaging module is used to: directly perform denoising, velocity analysis, iteration of residual static correction, and time-domain imaging processing on the shear wave data and the compressional wave data after pseudo-sampling processing.

17. The device according to claim 10, characterized in that, The data processing module is specifically configured to: When performing velocity analysis processing, residual static correction iteration processing, and time-domain imaging processing on seismic data, reduce the velocity in the corresponding relationship between the time and velocity of the pseudo-sampling processed data by a first preset multiple, and increase the time by a second preset multiple, so as to obtain the corresponding relationship between the time and velocity of the original shear wave data in the corresponding seismic data; According to the corresponding relationship between the time and velocity of the original shear wave data in the corresponding seismic data, perform velocity analysis processing, residual static correction iteration processing, and time-domain imaging processing on the seismic data.

18. The device according to claim 10, characterized in that, The data processing module is specifically configured to: When performing velocity analysis processing, residual static correction iteration processing, and time-domain imaging processing on seismic data, increase the velocity in the corresponding relationship between the time and velocity of the pseudo-sampling processed data by a third preset multiple, and reduce the time by a fourth preset multiple, so as to obtain the corresponding relationship between the time and velocity of the original compressional wave data in the corresponding seismic data; According to the corresponding relationship between the time and velocity of the original compressional wave data in the corresponding seismic data, perform velocity analysis processing, residual static correction iteration processing, and time-domain imaging processing on the seismic data.

19. A computer device, comprising a memory, a processor, and a computer program stored on the memory and executable on the processor, characterized in that, When the processor executes the computer program, it implements the method according to any one of claims 1 to 9.

20. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores a computer program, and when the computer program is executed by the processor, it implements the method according to any one of claims 1 to 9.

Images