Land longitudinal wave and converted transverse wave river channel sand body comparison interpretation method and system, electronic equipment and storage medium

Abstract

The invention provides a land longitudinal wave and converted transverse wave river channel sand body comparison interpretation method and system, electronic equipment and a storage medium. The method comprises the steps that longitudinal wave AVO attributes and converted wave channel integral attributes are comprehensively utilized to jointly solve river channel sand body recognition and depiction, on one hand, the effectiveness of medium impedance sand body recognition through converted transverse waves is utilized, so that the problem that medium impedance sand body recognition through longitudinal waves is insufficient is solved, comprehensive recognition of high, medium and low impedance sand bodies is achieved, and the sand body recognition range is expanded; secondly, the sand body period is determined on the basis that the transverse range is determined by utilizing the advantage of relatively high longitudinal resolution of longitudinal wave data; and 3, the wave crest horizon of the channel integral represents the sand body, and the wave crest of the longitudinal wave profile represents the wave impedance interface, so that the three-dimensional engraving of the sand body is facilitated, the joint interpretation precision of multi-wave seismic data is improved, and a foundation is laid for subsequent effective solving of geologic body description, reservoir dessert prediction and fluid detection.

Description

Technical Field

[0001] The present invention belongs to the technical field of geophysical exploration, and particularly relates to a method, a system, an electronic device and a storage medium for comparing and interpreting land longitudinal waves and converted shear waves in channel sand bodies. Background Technique

[0002] Tight sandstone gas (referred to as tight gas) belongs to the category of unconventional natural gas. Currently, it has the largest scale in the development of unconventional natural gas. In recent years, rapid development has been achieved in both the theory and technology of tight gas exploration and development.

[0003] In order to support the goal of increasing reserves and production in this field, targeted technical research and application have been carried out using longitudinal wave seismic data, and effective breakthroughs have been made in aspects such as seismic data amplitude and fidelity processing, channel identification and characterization. However, with the deepening of application and actual drilling, it is found that due to the large variation in the velocity of tight sand bodies, there are deficiencies in depicting the sand body boundary using longitudinal wave seismic data. Especially when the sandstone is a medium-impedance sand body and its longitudinal wave velocity is close to that of the surrounding rock mudstone, the longitudinal wave seismic response shows weak reflection, which is easily regarded as a non-sand body of "dark spot" and ignored, affecting the identification effect of the sand body.

[0004] Land multi-wave exploration technology is an important technology for comprehensively using longitudinal waves and converted shear waves to carry out oil and gas exploration. The important difference from conventional single longitudinal wave exploration is that it uses a single source for excitation and a three-component geophone for reception, and can collect the three-component wave field information of underground geological target bodies, providing the possibility for accurately describing the structural form of underground geological bodies, accurately predicting fine fractures, and finely depicting the spatial distribution of oil and gas reservoirs.

[0005] This technology has experienced more than forty years of development and application. It has certain advantages in aspects such as lithology and reservoir identification, and is one of the effective means to improve the exploration and development accuracy of complex oil and gas reservoirs such as tight and low-permeability reservoirs. Its technical advantages mainly include the following aspects:

[0006] (1) For areas such as gas cloud areas, tight gas reservoirs in clastic rocks, carbonate reservoirs, and shale oil and gas reservoirs, the longitudinal wave data has characteristics such as low signal-to-noise ratio and weak energy, while using the converted shear wave information can achieve accurate imaging and reservoir characterization.

[0007] (2) The wave field of land multi-wave seismic exploration data can simultaneously obtain longitudinal wave and converted shear wave field information. In addition to accurately extracting the longitudinal wave velocity parameters, it can also more accurately obtain the shear wave velocity parameters, and further extract accurate reservoir elastic parameters such as Poisson's ratio, longitudinal-to-shear wave velocity ratio, Young's modulus, and brittleness index, which are used to carry out lithology prediction and reservoir characteristic description work.

[0008] (3) The converted shear wave is essentially a shear wave and is more sensitive to reservoir fractures, which can better solve the difficult problems of quantitatively predicting the development characteristics of reservoir fractures and fluid detection by using the shear wave splitting characteristics.

[0009] Therefore, in order to meet the high-precision exploration and development of channel-type tight lithologic gas reservoirs in the Shaximiao Formation, it is proposed to comprehensively utilize the advantages of P-waves and converted shear waves to finely depict the spatial distribution of channel sand bodies, further supporting new breakthroughs in the exploration and development of tight gas in the Shaximiao Formation.

[0010] In land multi-wave seismic exploration, the joint contrast interpretation of P-waves and converted shear waves is a bottleneck problem restricting the development of this technology. Especially in the characterization of tight gas channel sand bodies, due to the different kinematic and dynamic principles of the two, mainly including: the P-wave velocity is higher than the converted shear wave velocity, and the propagation time is less than that of the converted shear wave; the propagation of P-waves is easily affected by fluids, while the converted shear wave only propagates in the rock skeleton and is less affected by fluids; the main frequency and frequency band width of P-waves are higher than those of the converted wave. Therefore, there are differences in the geophysical responses of the two wave fields to the amplitude strength, phase change, and vertical resolution of channel sand bodies. How to establish an interpretation bridge between the two has become the key to the application of this technology.

[0011] Currently, the commonly used methods in the interpretation of tight gas channels are mainly: analyzing the logging response characteristics of the target horizon and target body based on logging curve data, and then combining the synthetic seismogram calibration with P-wave and converted wave seismic profiles. Then, horizon interpretation and sand body interpretation are carried out respectively according to the calibration results. Finally, planar results or three-dimensional solid results of sand body characterization are formed. In practical applications, the defect of this method is that due to the different propagation times, seismic responses, and resolutions of the two wave fields, the contrast interpretation is relatively independent, resulting in certain errors, low accuracy in the characterization of channel sand bodies, and low efficiency in repeated contrast interpretation.

[0012] The above technical problems need to be solved urgently. Summary of the Invention

[0013] To solve the above technical problems, the present invention proposes a method, system, electronic device, and storage medium for contrasting and interpreting channel sand bodies of land P-waves and converted shear waves to solve the above technical problems.

[0014] The first aspect of the present invention discloses a method for contrasting and interpreting channel sand bodies of land P-waves and converted shear waves, the method comprising:

[0015] Step S1, performing prestack time migration processing on P-wave and converted shear wave seismic data to obtain a prestack time migration gather of P-waves, a P-wave profile, and a converted shear wave profile;

[0016] Step S2, loading longitudinal wave profile data, converting shear wave profile data and loading logging curves; loading geological stratification data of target formations and sand bodies to obtain geological stratification;

[0017] Step S3, using the logging curve, carry out the calibration of the P-wave, converted wave and shear wave synthesis record; clarify the time position of the geological layer at the logging well point used on the corresponding seismic section;

[0018] Step S4, clarifying the seismic response characteristics of different target layers and sand bodies on the P-wave profile and the converted S-wave profile according to the synthetic record calibration and the time position;

[0019] Step S5, according to the seismic response characteristics, performing target layer tracking interpretation in the entire area on the P-wave section and the converted S-wave section;

[0020] Step S6, performing AV0 attribute analysis on the P-wave prestack time migration gather data, extracting intercept and gradient attributes, and obtaining a PG attribute data volume by calculation;

[0021] Step S7, performing channel integral attribute calculations on the PG attribute data volume and the converted shear wave data volume respectively to obtain a PG channel integral data volume and a converted shear wave channel integral data volume;

[0022] Step S8, according to the result of the target layer tracking interpretation, based on the converted shear wave channel integral data volume, interpret the sand body layer, and save the converted shear wave sand body layer picking result;

[0023] Step S9, projecting the converted shear wave sand body layer position picking result onto the PG channel integral data volume, performing correction and supplement on the PG channel integral data volume, and saving the PG channel integral data volume sand body layer position picking result;

[0024] Step S10, projecting the sand body layer picking result of the PG channel integral data volume onto the P-wave section, picking up the top and bottom interfaces of the sand body along the crests and troughs on the P-wave section upward and downward according to the sand body layer of the PG channel integral data volume, and carrying out the interpretation and characterization of the top and bottom interfaces of the sand body.

[0025] According to the method of the first aspect of the present invention, in step S1, the method of performing prestack time migration processing on P-wave and converted S-wave seismic data to obtain P-wave prestack time migration gathers, P-wave profiles and converted S-wave profiles comprises:

[0026] Static correction, noise attenuation, amplitude compensation and consistency, horizontal stacking and residual static correction, prestack time migration and migration stacking are performed on P-wave and converted S-wave seismic data to obtain P-wave prestack time migration gathers and converted S-wave prestack time migration gathers.

[0027] Perform post-stack high-resolution and display optimization processing on the longitudinal wave pre-stack time migration gather and converted shear wave pre-stack time migration gather to obtain a longitudinal wave profile and a converted shear wave profile.

[0028] According to the method of the first aspect of the present invention, in the step S2, the method of using well logging curves to carry out calibration of synthetic records of longitudinal waves and converted shear waves; and clarifying the time positions of geological horizons at the well logging well points adopted on the corresponding seismic profiles includes:

[0029] Calculate the reflection times of each layer of longitudinal waves and the reflection times of each layer of shear waves according to the longitudinal wave velocity, shear wave velocity and density of the well logging curves.

[0030] Carry out spectral analysis of the target horizons in the longitudinal wave profile data and the converted shear wave profile data to determine the main frequency of the longitudinal wave data and the main frequency of the converted shear wave data.

[0031] Respectively set the values of the main frequencies of the longitudinal wave data and the converted shear wave data, respectively set two wavelets of the main frequencies of the longitudinal wave data and the converted shear wave data, and combine the longitudinal wave velocity, shear wave velocity and density values in the well logging curves to make theoretical synthetic records.

[0032] Use the synthetic records to compare and calibrate with the longitudinal wave profile and the converted shear wave profile, and project the geological stratification onto the seismic profile, so as to clarify the time positions of the geological horizons at the well logging well points adopted on the corresponding seismic profiles.

[0033] According to the method of the first aspect of the present invention, in the step S4, the method of clarifying the seismic response characteristics of different target horizons and sand bodies on the longitudinal wave profile and the converted shear wave profile according to the synthetic record calibration and the time positions includes:

[0034] Check and match the geological horizons with the gamma curve, porosity curve, shale content and water saturation in the well logging curves.

[0035] According to the inspection and matching results, form the seismic response characteristics of the longitudinal wave and the converted shear wave to the sand body, including: the longitudinal wave is a wave crest or a wave trough, and the converted shear wave is a wave crest or a wave trough; the amplitude of the longitudinal wave is strong or weak, and the amplitude of the converted shear wave is strong or weak.

[0036] According to the method of the first aspect of the present invention, in the step S5, the method of performing full-area target horizon tracking and interpretation on the longitudinal wave profile and the converted shear wave profile according to the seismic response characteristics includes:

[0037] According to the seismic response characteristics, carry out horizon interpretation at all well point positions in the work area.

[0038] According to the horizon interpretation at the well point positions, carry out horizon interpretation of the cross-well seismic data in the work area.

[0039] According to the horizon interpretation of the cross-well seismic data, horizon picking is carried out for grids of 80×80, 40×40, and 20×20 in the work area, and finally, it is encrypted and interpolated into horizon picking of 1×1 grid.

[0040] According to the method of the first aspect of the present invention, in the step S6, the method for performing AVO attribute analysis on the prestack time migration gather data of the P-wave, extracting the intercept and gradient attributes, and obtaining the P-G attribute data volume by calculation includes:

[0041]

[0042] Wherein, P represents the intercept; G represents the gradient attribute; PG represents the P-G attribute data volume; v s represents the shear wave velocity; Δv S represents the shear wave velocity difference between the upper and lower strata; ρ represents the density; Δρ represents the density difference between the upper and lower strata.

[0043] According to the method of the first aspect of the present invention, in the step S7, the method for performing trace integral attribute calculation on the P-G attribute data volume and the converted shear wave data volume respectively to obtain the P-G trace integral data volume and the converted shear wave trace integral data volume includes:

[0044]

[0045] Wherein, Spg represents the P-G trace integral data volume; PG(t) represents the P-G attribute data volume; wpg(t) is the wavelet of the P-G attribute data volume, and k is the reflection coefficient of the top interface;

[0046]

[0047] Wherein, Sps represents the converted shear wave trace integral data volume; PS(t) represents the converted shear wave data volume; wps(t) is the wavelet of the converted shear wave data volume.

[0048] The second aspect of the present invention discloses a land P-wave and converted shear wave channel sand body contrast interpretation system, and the system includes:

[0049] A first processing module, configured to perform prestack time migration processing on P-wave and converted shear wave seismic data to obtain a prestack time migration gather of the P-wave, a P-wave profile, and a converted shear wave profile;

[0050] A second processing module, configured to load P-wave profile data, converted shear wave profile data, and load logging curves; load geological stratification data of the target formation and sand bodies to obtain geological stratification;

[0051] The third processing module is configured to use well logging curves to carry out calibration of the synthetic records of P-waves and converted S-waves; and determine the time positions of the geological horizons at the well logging points adopted on the corresponding seismic profiles.

[0052] The fourth processing module is configured to determine the seismic response characteristics of different target horizons and sand bodies on the P-wave profile and the converted S-wave profile according to the synthetic record calibration and the time positions.

[0053] The fifth processing module is configured to perform full-area target horizon tracing and interpretation on the P-wave profile and the converted S-wave profile according to the seismic response characteristics.

[0054] The sixth processing module is configured to perform AVO attribute analysis on the pre-stack time migration gather data of P-waves, extract intercept and gradient attributes, and obtain a P-G attribute data volume through calculation.

[0055] The seventh processing module is configured to perform trace integration attribute calculation on the P-G attribute data volume and the converted S-wave data volume respectively to obtain a P-G trace integration data volume and a converted S-wave trace integration data volume.

[0056] The eighth processing module is configured to interpret the sand body horizons based on the converted S-wave trace integration data volume according to the results of the target horizon tracing and interpretation, and save the picking results of the converted S-wave sand body horizons.

[0057] The ninth processing module is configured to project the picking results of the converted S-wave sand body horizons onto the P-G trace integration data volume, make corrections and supplements on the P-G trace integration data volume, and save the picking results of the sand body horizons of the P-G trace integration data volume.

[0058] The tenth processing module is configured to project the picking results of the sand body horizons of the P-G trace integration data volume onto the P-wave profile, pick the top and bottom interfaces of the sand body along the wave peaks and wave valleys on the P-wave profile upward and downward according to the sand body horizons of the P-G trace integration data volume, and carry out interpretation and characterization of the top and bottom interfaces of the sand body.

[0059] A third aspect of the present invention discloses an electronic device. The electronic device includes a memory and a processor. The memory stores a computer program. When the processor executes the computer program, the steps in any one of the methods for comparing and interpreting land P-wave and converted S-wave channel sand bodies in the first aspect of the present disclosure are implemented.

[0060] A fourth aspect of the present invention discloses a computer-readable storage medium. A computer program is stored on the computer-readable storage medium. When the computer program is executed by a processor, the steps in any one of the methods for comparing and interpreting land P-wave and converted S-wave channel sand bodies in the first aspect of the present disclosure are implemented.

[0061] In summary, the solution proposed by the present invention is based on the theoretical basis that the P-G attribute in the AVO analysis of longitudinal wave seismic data has a correspondence with the converted shear wave, and the trace integration is to integrate the seismic trace data to quickly obtain the relative wave impedance profile. It is a direct inversion method and an effective formation lithology interpretation technique for tracing sand bodies, and it is easier to directly distinguish the sand body range compared with the original profile. Thus, a new idea for jointly and integrally depicting channel sand bodies with multi-wave data on land is established. It can effectively establish the idea of jointly identifying and interpreting sand bodies with longitudinal waves and converted shear waves, and fully utilize the information of the two wave fields: on the one hand, utilize the effectiveness of the converted shear wave in identifying medium-impedance sand bodies to solve the problem of insufficient identification of medium-impedance sand bodies by longitudinal waves, so as to comprehensively identify high-, medium-, and low-impedance sand bodies and expand the sand body identification range. On the other hand, utilize the advantage of higher vertical resolution of longitudinal wave data to determine the sand body stage on the basis of the determined lateral range. Compared with traditional methods and techniques, the implementation mode of the present invention is more effective, the implementation mode is convenient, the method basis is reliable, the application effect is obvious, it can effectively solve the problem of the accuracy of depicting tight gas channel sand bodies in land multi-wave seismic exploration, and improve the promotion and application value of the present invention. BRIEF DESCRIPTION OF THE DRAWINGS

[0062] In order to more clearly illustrate the specific embodiments of the present invention or the technical solutions in the prior art, the following will briefly introduce the drawings required for the description of the specific embodiments or the prior art. Obviously, the drawings in the following description are 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.

[0063] Figure 1 FIG. is a flow chart of a method for comparative interpretation of channel sand bodies with longitudinal waves and converted shear waves on land according to an embodiment of the present invention;

[0064] Figure 2 FIG. is the depiction and characterization of channel sand bodies with longitudinal wave, converted shear wave trace integration data according to an embodiment of the present invention;

[0065] Figure 3 FIG. is the calibration of channel sand bodies with longitudinal waves, converted shear waves and trace integration according to an embodiment of the present invention;

[0066] Figure 4 FIG. is the effect of sand body identification and depiction according to an embodiment of the present invention;

[0067] Figure 5 FIG. is the effect of channel sand body depiction according to an embodiment of the present invention;

[0068] Figure 6 FIG. is a structural diagram of a system for comparative interpretation of channel sand bodies with longitudinal waves and converted shear waves on land according to an embodiment of the present invention;

[0069] Figure 7 Structural diagram of an electronic device according to an embodiment of the present invention. Detailed implementation manners

[0070] To make the objectives, technical solutions and advantages of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention. Apparently, the described embodiments are some but not all of the embodiments of the present invention. All other embodiments obtained by those of ordinary skill in the art based on the embodiments of the present invention without making creative efforts shall fall within the protection scope of the present invention.

[0071] The first aspect of the present invention discloses a method for comparing and interpreting land P-wave and converted S-wave channel sand bodies. Figure 1 Flowchart of a method for comparing and interpreting land P-wave and converted S-wave channel sand bodies according to an embodiment of the present invention, as Figures 1 to 5 shown, the method includes:

[0072] Step S1: Perform prestack time migration processing on P-wave and converted S-wave seismic data to obtain a prestack time migration gather of P-wave, a P-wave profile, and a converted S-wave profile;

[0073] Step S2: Load P-wave profile data, converted S-wave profile data, and logging curves; load geological stratification data of the target formation and sand bodies to obtain geological stratification;

[0074] Step S3: Use logging curves to carry out calibration of P-wave and converted S-wave synthetic seismograms; clarify the time positions of the geological horizons at the logging well points adopted on the corresponding seismic profiles;

[0075] Step S4: According to the synthetic seismogram calibration and the time positions, clarify the seismic response characteristics of different target horizons and sand bodies on the P-wave profile and the converted S-wave profile;

[0076] Step S5: According to the seismic response characteristics, perform full-area target horizon tracking and interpretation on the P-wave profile and the converted S-wave profile;

[0077] Step S6: Perform AV0 attribute analysis on the prestack time migration gather data of P-wave, extract intercept and gradient attributes, and obtain a P-G attribute data volume through calculation;

[0078] Step S7: Perform trace integration attribute calculation on the P-G attribute data volume and the converted S-wave data volume respectively to obtain a P-G trace integration data volume and a converted S-wave trace integration data volume;

[0079] Step S8: According to the results of the target horizon tracking and interpretation, based on the converted S-wave trace integration data volume, interpret the sand body horizons and save the picking results of the converted S-wave sand body horizons;

[0080] Step S9. Project the converted shear wave sand body horizon picking result onto the P-G trace integral data volume, make corrections and supplements on the P-G trace integral data volume, and save the picking result of the sand body horizon of the P-G trace integral data volume.

[0081] Step S10. Project the picking result of the sand body horizon of the P-G trace integral data volume onto the P-wave section, pick the top and bottom interfaces of the sand body along the wave peaks and troughs on the P-wave section upward and downward according to the picking result of the sand body horizon of the P-G trace integral data volume, and carry out the interpretation and characterization of the top and bottom interfaces of the sand body.

[0082] In step S1, perform prestack time migration processing on the P-wave and converted shear wave seismic data to obtain the P-wave prestack time migration gather, the P-wave section, and the converted shear wave section.

[0083] In some embodiments, in the step S1, the method for performing prestack time migration processing on the P-wave and converted shear wave seismic data to obtain the P-wave prestack time migration gather, the P-wave section, and the converted shear wave section includes:

[0084] Perform static correction processing, noise attenuation processing, amplitude compensation and consistency processing, horizontal stacking and residual static correction processing, and prestack time migration processing and migration stacking on the P-wave and converted shear wave seismic data to obtain the P-wave prestack time migration gather and the converted shear wave prestack time migration gather;

[0085] Perform post-stack high-resolution and display optimization processing on the P-wave prestack time migration gather and the converted shear wave prestack time migration gather to obtain the P-wave section and the converted shear wave section.

[0086] In step S2, load the P-wave section data, the converted shear wave section data, and load the logging curves; load the geological stratification data of the target formation and sand body to obtain the geological stratification.

[0087] Specifically, loading the logging curves mainly includes the P-wave velocity (Vp), the shear wave velocity (Vs), the density (ρ), the gamma curve γ, and the porosity curve the shale content (Vsh) and the water saturation (Sw).

[0088] In step S3, use the logging curves to carry out the calibration of the P-wave and converted shear wave synthetic seismograms; clarify the time positions of the geological horizons at the logging well points adopted on the corresponding seismic sections.

[0089] In some embodiments, in the step S3, the method for using the logging curves to carry out the calibration of the P-wave and converted shear wave synthetic seismograms; clarify the time positions of the geological horizons at the logging well points adopted on the corresponding seismic sections includes:

[0090] Calculate the reflection time of each layer of the longitudinal wave and the reflection time of each layer of the shear wave based on the longitudinal wave velocity, shear wave velocity, and density of the logging curve;

[0091] Carry out spectral analysis of the target layer in the longitudinal wave profile data and the converted shear wave profile data to determine the main frequency of the longitudinal wave data and the main frequency of the converted shear wave data;

[0092] Respectively set the value of the main frequency of the longitudinal wave data and the value of the main frequency of the converted shear wave data, respectively set two wavelets of the main frequency of the longitudinal wave data and the main frequency of the converted shear wave data, and combine the longitudinal wave velocity, shear wave velocity, and density values in the logging curve to produce a theoretical synthetic record;

[0093] Use the synthetic record to compare and calibrate with the longitudinal wave profile and the converted shear wave profile, and project the geological stratification onto the seismic profile, so as to clarify the time position of the geological horizons at the logging well points adopted on the corresponding seismic profile.

[0094] Specifically, the method for calculating the reflection time of each layer of the longitudinal wave and the reflection time of each layer of the shear wave based on the longitudinal wave velocity, shear wave velocity, and density of the logging curve includes:

[0095]

[0096] Among them, H p，i is the formation depth of the i-th layer of the longitudinal wave, H p，i+1 is the formation depth of the (i + 1)-th layer of the longitudinal wave, T p，i is the reflection time of the i-th layer of the longitudinal wave, V p，i+1 is the longitudinal wave layer velocity of the (i + 1)-th layer of the longitudinal wave, V p,i is the longitudinal wave layer velocity of the i-th layer;

[0097]

[0098] Among them, H s，i is the formation depth of the i-th layer of the shear wave, H s，i+1 is the formation depth of the (i + 1)-th layer of the shear wave, T s，i is the reflection time of the i-th layer of the shear wave, V s，i+1 is the shear wave layer velocity of the (i + 1)-th layer of the shear wave, V s,i is the shear wave layer velocity of the i-th layer.

[0099] In step S4, based on the synthetic record calibration and the time position, clarify the seismic response characteristics of different target horizons and sand bodies on the longitudinal wave profile and the converted shear wave profile.

[0100] In some embodiments, in the step S4, the method for clarifying the seismic response characteristics of different target horizons and sand bodies on the longitudinal wave profile and the converted shear wave profile based on the synthetic record calibration and the time position includes:

[0101] Check and match the geological horizons with gamma ray curves, porosity curves, shale content, and water saturation in well logging curves;

[0102] According to the check and match results, form the seismic response characteristics of the P-wave and converted S-wave to sand bodies, including: whether the P-wave is a peak or a trough, and whether the converted S-wave is a peak or a trough; whether the P-wave amplitude is strong or weak, and whether the converted S-wave amplitude is strong or weak.

[0103] In step S5, according to the seismic response characteristics, perform full-area target horizon tracing and interpretation on the P-wave profile and the converted S-wave profile.

[0104] In some embodiments, in step S5, the method of performing full-area target horizon tracing and interpretation on the P-wave profile and the converted S-wave profile according to the seismic response characteristics includes:

[0105] According to the seismic response characteristics, carry out horizon interpretation at all well point positions in the work area;

[0106] According to the horizon interpretation at the well point positions, carry out horizon interpretation on the cross-well seismic data in the work area;

[0107] According to the horizon interpretation of the cross-well seismic data, carry out horizon picking for the 80×80, 40×40, and 20×20 grids in the work area, and finally interpolate and encrypt it into horizon picking for the 1×1 grid.

[0108] In step S6, perform AVO attribute analysis on the P-wave pre-stack time migration gather data, extract the intercept and gradient attributes, and obtain the P-G attribute data volume through calculation.

[0109] In some embodiments, in step S6, the method of performing AVO attribute analysis on the P-wave pre-stack time migration gather data, extracting the intercept and gradient attributes, and obtaining the P-G attribute data volume through calculation includes:

[0110]

[0111] Where P represents the intercept; G represents the gradient attribute; PG represents the P-G attribute data volume; v s represents the shear wave velocity; Δv S represents the shear wave velocity difference between the upper and lower strata; ρ represents the density; Δρ represents the density difference between the upper and lower strata.

[0112] Specifically, perform angle domain conversion processing on the P-wave pre-stack time migration gather Gather_P to form the P-wave angle domain gather data Angel_P. Use the following formula:

[0113]

[0114] Among them, V0 and k in the formula can be obtained by least squares fitting.

[0115]

[0116]

[0117] In the formula: x is the shot-receiver offset; Z is the depth of the target layer, n is the number of sample points, v i ,z i are the velocity and depth at each sample point respectively;

[0118] According to the P-wave angle domain gather data, intercept (P attribute) and gradient attribute (G attribute) in AVO attribute analysis are extracted within the target layer range:

[0119] R(α) = P + G sin 2 α

[0120] P is the intercept of this linear equation; G is the slope or gradient of this equation. Under certain assumptions, that is, Δρ, Δv P , Δv S are relatively small compared to ρ, v P , v S respectively, and v p / v s = 2, thus obtaining:

[0121]

[0122]

[0123] Calculate to obtain the P-G attribute data volume:

[0124]

[0125] In step S7, perform trace integration attribute calculations on the P-G attribute data volume and the converted shear wave data volume respectively to obtain the P-G trace integration data volume and the converted shear wave trace integration data volume.

[0126] In some embodiments, in the step S7, the method of performing trace integration attribute calculations on the P-G attribute data volume and the converted shear wave data volume respectively to obtain the P-G trace integration data volume and the converted shear wave trace integration data volume includes:

[0127]

[0128] Among them, Spg represents the P-G trace integration data volume; PG(t) represents the P-G attribute data volume; wpg(t) is the wavelet of the P-G attribute data volume, and k is the reflection coefficient of the top interface;

[0129]

[0130] Among them, Sps represents the converted shear wave channel integral data volume; PS(t) represents the converted shear wave data volume; and wps(t) is the wavelet of the converted shear wave data volume.

[0131] In step S10, project the sand body horizon picking result of the P-G channel integral data volume onto the longitudinal wave section, pick the top and bottom interfaces of the sand body along the wave peaks and troughs on the longitudinal wave section upward and downward according to the sand body horizon of the P-G channel integral data volume, and carry out the interpretation and characterization of the top and bottom interfaces of the sand body.

[0132] Specifically, save the top and bottom interfaces of the sand body; use the attribute interpretation and display function to complete the three-dimensional characterization of the sand body.

[0133] In summary, the solution proposed by the present invention can adopt the theoretical basis that the P-G attribute in the AVO analysis of longitudinal wave seismic data has a correspondence with the converted shear wave, and the trace integration is to integrate the seismic trace data to quickly obtain the relative wave impedance profile, which is a direct inversion method and an effective formation lithology interpretation technique for tracing sand bodies, and it is easier to directly distinguish the sand body range compared with the original profile. Thus, a new idea for the joint integrated characterization of land multi-wave data of river channel sand bodies is established. It can effectively establish the joint sand body identification and interpretation idea of longitudinal wave and converted shear wave, and fully apply the information of the two wave fields: on the one hand, utilize the effectiveness of the converted shear wave in identifying medium impedance sand bodies, so as to solve the problem of insufficient identification of medium impedance sand bodies by longitudinal waves, and achieve the comprehensive identification of high, medium, and low impedance sand bodies and expand the sand body identification range. On the other hand, utilize the advantage of the higher longitudinal resolution of the longitudinal wave data to determine the sand body stages on the basis of the determined lateral range. Compared with the traditional method and technology, the implementation mode of the present invention is more effective, the implementation mode is convenient, the method basis is reliable, the application effect is obvious, it can effectively solve the problem of the characterization accuracy of tight gas river channel sand bodies in land multi-wave seismic exploration, and improve the promotion and application value of the present invention.

[0134] The second aspect of the present invention discloses a land longitudinal wave and converted shear wave river channel sand body comparison and interpretation system. Figure 6 The structural diagram of a land longitudinal wave and converted shear wave river channel sand body comparison and interpretation system according to an embodiment of the present invention; as Figure 6 shown, the system 100 includes:

[0135] The first processing module 101 is configured to perform pre-stack time migration processing on the longitudinal wave and converted shear wave seismic data to obtain a pre-stack time migration trace gather of the longitudinal wave, a longitudinal wave section, and a converted shear wave section;

[0136] The second processing module 102 is configured to load the longitudinal wave section data, the converted shear wave section data, and load the logging curves; load the geological stratification data of the target formation and sand bodies to obtain the geological stratification;

[0137] The third processing module 103 is configured to use the logging curve to carry out the calibration of the P-wave, converted wave and shear wave synthesis record; to clarify the time position of the geological layer at the logging well point used on the corresponding seismic section;

[0138] The fourth processing module 104 is configured to clarify the seismic response characteristics of different target layers and sand bodies on the longitudinal wave section and the converted shear wave section according to the synthetic record calibration and the time position;

[0139] The fifth processing module 105 is configured to perform target layer tracking interpretation of the entire area on the longitudinal wave section and the converted shear wave section according to the seismic response characteristics;

[0140] The sixth processing module 106 is configured to perform AV0 attribute analysis on the P-wave prestack time migration gather data, extract intercept and gradient attributes, and obtain a PG attribute data volume by calculation;

[0141] The seventh processing module 107 is configured to perform channel integral attribute calculation on the PG attribute data volume and the converted shear wave data volume respectively to obtain a PG channel integral data volume and a converted shear wave channel integral data volume;

[0142] The eighth processing module 108 is configured to interpret the sand body layer according to the result of the target layer tracking interpretation and based on the converted shear wave channel integral data volume, and save the converted shear wave sand body layer picking result;

[0143] The ninth processing module 109 is configured to project the converted shear wave sand body layer position picking result onto the PG channel integral data volume, perform correction and supplement on the PG channel integral data volume, and save the PG channel integral data volume sand body layer position picking result;

[0144] The tenth processing module 110 is configured to project the sand body layer picking results of the PG channel integral data volume onto the longitudinal wave section, pick up the top and bottom interfaces of the sand body upward and downward along the crests and troughs on the longitudinal wave section according to the sand body layer of the PG channel integral data volume, and carry out the interpretation and characterization of the top and bottom interfaces of the sand body.

[0145] According to the system of the second aspect of the present invention, the first processing module 101 is specifically configured as follows: the method of performing prestack time migration processing on P-wave and converted S-wave seismic data to obtain P-wave prestack time migration gathers, P-wave profiles and converted S-wave profiles comprises:

[0146] Static correction processing, noise attenuation processing, amplitude compensation and consistency processing, horizontal stacking and residual static correction processing, and prestack time migration processing and migration stacking are performed on the P-wave and converted S-wave seismic data to obtain the P-wave prestack time migration gather and the converted S-wave prestack time migration gather;

[0147] Post-stack high-resolution and display optimization processing are performed on the P-wave prestack time migration gather and the converted S-wave prestack time migration gather to obtain the P-wave profile and the converted S-wave profile.

[0148] According to the system of the second aspect of the present invention, the second processing module 102 is specifically configured to load well logging curves, mainly including the P-wave velocity (Vp), S-wave velocity (Vs), density (ρ), gamma curve γ, porosity curve shale content (Vsh) and water saturation (Sw).

[0149] According to the system of the second aspect of the present invention, the third processing module 103 is specifically configured to use the well logging curves to perform calibration of synthetic P-wave and converted S-wave records; the method for determining the time position of the geological horizons at the well logging well points adopted on the corresponding seismic profiles includes:

[0150] Calculate the reflection time of each P-wave layer and the reflection time of each S-wave layer according to the P-wave velocity, S-wave velocity, and density of the well logging curves;

[0151] Perform spectral analysis on the target layers in the P-wave profile data and the converted S-wave profile data to determine the dominant frequency of the P-wave data and the dominant frequency of the converted S-wave data;

[0152] Respectively set the values of the dominant frequencies of the P-wave data and the converted S-wave data, respectively set two wavelets of the dominant frequencies of the P-wave data and the converted S-wave data, and combine the P-wave velocity, S-wave velocity, and density values in the well logging curves to produce theoretical synthetic records;

[0153] Use the synthetic records to perform comparison and calibration with the P-wave profile and the converted S-wave profile, and project the geological stratification onto the seismic profile, so as to clarify the time position of the geological horizons at the well logging well points adopted on the corresponding seismic profiles.

[0154] Specifically, the method for calculating the reflection time of each P-wave layer and the reflection time of each S-wave layer according to the P-wave velocity, S-wave velocity, and density of the well logging curves includes:

[0155]

[0156] where H p，i is the formation depth of the i-th P-wave layer, H p，i+1 is the formation depth of the (i + 1)-th P-wave layer, T p，iis the reflection time of the i-th layer of the P-wave, V p，i+1 is the P-wave layer velocity of the (i + 1)-th layer of the P-wave, V p,i is the P-wave layer velocity of the i-th layer;

[0157]

[0158] wherein, H s，i is the formation depth of the i-th layer of the S-wave, H s，i+1 is the formation depth of the (i + 1)-th layer of the S-wave, T s，i is the reflection time of the i-th layer of the S-wave, V s，i+1 is the S-wave layer velocity of the (i + 1)-th layer of the S-wave, V s,i is the S-wave layer velocity of the i-th layer.

[0159] According to the system of the second aspect of the present invention, the fourth processing module 104 is specifically configured that the method for clarifying the seismic response characteristics of different target horizons and sand bodies on the P-wave profile and the converted S-wave profile according to the synthetic record calibration and the time position includes:

[0160] Check and match the geological horizons with the gamma curve, porosity curve, shale content, and water saturation in the logging curves;

[0161] According to the check and match results, form the seismic response characteristics of the P-wave and the converted wave to the sand body, including: the P-wave is a peak or a trough, and the converted S-wave is a peak or a trough; the amplitude of the P-wave is strong or weak, and the amplitude of the converted S-wave is strong or weak.

[0162] According to the system of the second aspect of the present invention, the fifth processing module 105 is specifically configured that the method for performing full-area target horizon tracking and interpretation on the P-wave profile and the converted S-wave profile according to the seismic response characteristics includes:

[0163] According to the seismic response characteristics, carry out horizon interpretation at all well point positions in the work area;

[0164] According to the horizon interpretation at the well point positions, carry out horizon interpretation of the cross-well seismic data in the work area;

[0165] According to the horizon interpretation of the cross-well seismic data, carry out horizon picking for the 80×80, 40×40, and 20×20 grids in the work area, and finally interpolate and encrypt it into horizon picking for the 1×1 grid.

[0166] According to the system of the second aspect of the present invention, the sixth processing module 106 is specifically configured that the method for performing AV0 attribute analysis on the P-wave pre-stack time migration gather data, extracting the intercept and gradient attributes, and obtaining the P-G attribute data volume by calculation includes:

[0167]

[0168] Among them, P represents the intercept; G represents the gradient attribute; PG represents the P-G attribute data volume; v s represents the shear wave velocity; Δv S represents the shear wave velocity difference between the upper and lower strata; ρ represents the density; Δρ represents the density difference between the upper and lower strata.

[0169] Specifically, perform angle-domain conversion processing on the P-wave pre-stack time migration gather Gather_P to form the P-wave angle-domain gather data Angel_P. Using the following formula:

[0170]

[0171] Among them, V0 and k in the formula can be obtained by least squares fitting,

[0172]

[0173]

[0174] In the formula: x is the shot-receiver offset; Z is the depth of the target layer, n is the number of sample points, v i ,z i are the velocity and depth at each sample point respectively;

[0175] According to the P-wave angle-domain gather data, extract the intercept (P attribute) and gradient attribute (G attribute) in the AVO attribute analysis within the range of the target layer:

[0176] R(α) = P + G sin 2 α

[0177] P is the intercept of this linear equation; G is the slope or gradient of this equation. Under certain assumptions, that is, Δρ, Δv P , Δv S are relatively small compared to ρ, v P , v S respectively, and v p / v s = 2, thus obtaining:

[0178]

[0179]

[0180] Calculate to obtain the P-G attribute data volume:

[0181]

[0182] For the system according to the second aspect of the present invention, the seventh processing module 107 is specifically configured that the method for respectively performing trace integration attribute calculation on the P-G attribute data volume and the converted shear wave data volume to obtain the P-G trace integration data volume and the converted shear wave trace integration data volume includes:

[0183]

[0184] Among them, Spg represents the P-G trace integration data volume; PG(t) represents the P-G attribute data volume; wpg(t) is the wavelet of the P-G attribute data volume, and k is the reflection coefficient of the top interface;

[0185]

[0186] Among them, Sps represents the converted shear wave trace integration data volume; PS(t) represents the converted shear wave data volume; wps(t) is the wavelet of the converted shear wave data volume.

[0187] A third aspect of the present invention discloses an electronic device. The electronic device includes a memory and a processor. The memory stores a computer program. When the processor executes the computer program, the steps in any one of the first aspects of the present invention for a method of comparing and interpreting land longitudinal waves and converted shear wave channel sand bodies are implemented.

[0188] Figure 7 As a structural diagram of an electronic device according to an embodiment of the present invention, as Figure 7 shown, the electronic device includes a processor, a memory, a communication interface, a display screen, and an input device connected through a system bus. Among them, the processor of the electronic device is used to provide computing and control capabilities. The memory of the electronic device includes a non-volatile storage medium and an internal memory. The non-volatile storage medium stores an operating system and a computer program. The internal memory provides an environment for the operation of the operating system and the computer program in the non-volatile storage medium. The communication interface of the electronic device is used to communicate with an external terminal in a wired or wireless manner. The wireless manner can be implemented through WIFI, a carrier network, near field communication (NFC), or other technologies. The display screen of the electronic device can be a liquid crystal display screen or an electronic ink display screen. The input device of the electronic device can be a touch layer covering the display screen, or a button, a trackball, or a touchpad provided on the housing of the electronic device, or an external keyboard, touchpad, or mouse, etc.

[0189] Those skilled in the art can understand that Figure 7 the structure shown in is only a structural diagram of a part related to the technical solution of the present disclosure, and does not constitute a limitation on the electronic device to which the solution of the present application is applied. The specific electronic device may include more or fewer components than those shown in the figure, or combine some components, or have a different component layout.

[0190] A fourth aspect of the present invention discloses a computer-readable storage medium. A computer program is stored on the computer-readable storage medium. When the computer program is executed by a processor, the steps in a method for comparing and interpreting land P-waves and converted S-waves in channel sand bodies according to any one of the first aspects disclosed in the present invention are implemented.

[0191] Please note that the technical features of the above embodiments can be combined arbitrarily. For the sake of concise description, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, it should be considered as the scope described in this specification. The above embodiments only represent several implementation manners of the present application. The description is relatively specific and detailed, but it should not be construed as a limitation on the scope of the invention patent. It should be noted that for those of ordinary skill in the art, without departing from the concept of the present application, several modifications and improvements can be made, and these all belong to the protection scope of the present application. Therefore, the protection scope of the patent of the present application shall be subject to the appended claims.

[0192] The above are the preferred implementation manners of the present invention. It should be noted that for those of ordinary skill in the technical field, without departing from the principle of the present invention, several improvements and refinements can be made, and these improvements and refinements should also be regarded as the protection scope of the present invention.

Claims

1. A method for comparative interpretation of land P-wave and converted S-wave channel sand bodies, characterized in that The method includes: Step S1: Perform prestack time migration processing on the P-wave and converted S-wave seismic data to obtain a prestack time migration gather of P-waves, a P-wave section, and a converted S-wave section; Step S2: Load the P-wave section data, the converted S-wave section data, and the well logging curves; load the geological stratification data of the target formation and sand bodies to obtain geological stratification; Step S3: Use the well logging curves to carry out calibration of synthetic P-wave and converted S-wave records; determine the time positions of the geological horizons at the well logging well points adopted on the corresponding seismic sections; Step S4: According to the calibration of the synthetic records and the time positions, determine the seismic response characteristics of different target horizons and sand bodies on the P-wave section and the converted S-wave section; Step S5: According to the seismic response characteristics, perform full-area tracking and interpretation of the target horizons on the P-wave section and the converted S-wave section; Step S6: Perform AVO attribute analysis on the prestack time migration gather data of P-waves, extract the intercept and gradient attributes, and obtain a P-G attribute data volume through calculation; Step S7: Perform trace integration attribute calculations on the P-G attribute data volume and the converted S-wave data volume respectively to obtain a P-G trace integration data volume and a converted S-wave trace integration data volume; Step S8: According to the results of the target horizon tracking and interpretation, based on the converted S-wave trace integration data volume, interpret the sand body horizons and save the picking results of the converted S-wave sand body horizons; Step S9: According to the picking results of the converted S-wave sand body horizons, project them onto the P-G trace integration data volume, perform correction and supplementation on the P-G trace integration data volume, and save the picking results of the sand body horizons of the P-G trace integration data volume; Step S10: Project the picking results of the sand body horizons of the P-G trace integration data volume onto the P-wave section, and pick the top and bottom interfaces of the sand bodies along the wave peaks and troughs on the P-wave section upward and downward according to the sand body horizons of the P-G trace integration data volume, and carry out the interpretation and characterization of the top and bottom interfaces of the sand bodies.

2. The method for comparing and interpreting land P-wave and converted S-wave channel sand bodies according to claim 1, wherein In the said Step S1, the method for performing prestack time migration processing on the P-wave and converted S-wave seismic data to obtain a prestack time migration gather of P-waves, a P-wave section, and a converted S-wave section includes: Perform static correction processing, noise attenuation processing, amplitude compensation and consistency processing, horizontal stacking and residual static correction processing, and prestack time migration processing and migration stacking on the P-wave and converted S-wave seismic data to obtain a prestack time migration gather of P-waves and a prestack time migration gather of converted S-waves; Perform post-stack high-resolution and display optimization processing on the prestack time migration gather of P-waves and the prestack time migration gather of converted S-waves to obtain a P-wave section and a converted S-wave section.

3. A method for comparing and interpreting land P-waves and converted S-waves in channel sand bodies according to claim 1, characterized in that In the said Step S3, the method for using the well logging curves to carry out calibration of synthetic P-wave and converted S-wave records; determine the time positions of the geological horizons at the well logging well points adopted on the corresponding seismic sections includes: Calculate the reflection times of each layer of P-waves and the reflection times of each layer of S-waves according to the P-wave velocity, S-wave velocity, and density of the well logging curves; Carry out spectral analysis of the target horizons in the P-wave section data and the converted S-wave section data to determine the main frequency of the P-wave data and the main frequency of the converted S-wave data; Set the main frequency values of the longitudinal wave data and the converted shear wave data respectively, and set two wavelets for the main frequency of the longitudinal wave data and the main frequency of the converted shear wave data respectively. Combine the longitudinal wave velocity, shear wave velocity and density values in the logging curve to make a theoretical synthetic record; Use the synthetic record to compare and calibrate with the longitudinal wave profile and the converted shear wave profile, and project the geological stratification into the seismic profile, so as to clarify the time position of the geological horizon at the logging well point adopted on the corresponding seismic profile.

4. A method for comparing and interpreting channel sand bodies of land P-waves and converted S-waves according to claim 1, characterized in that, In the step S4, the method for clarifying the seismic response characteristics of different target horizons and sand bodies on the longitudinal wave profile and the converted shear wave profile according to the synthetic record calibration and the time position includes: Check and match the geological horizon with the gamma curve, porosity curve, shale content, and water saturation in the logging curve; According to the inspection and matching results, form the seismic response characteristics of the longitudinal wave and the converted wave to the sand body, including: the longitudinal wave is a peak or a trough, and the converted shear wave is a peak or a trough; the amplitude of the longitudinal wave is strong or weak, and the amplitude of the converted shear wave is strong or weak.

5. A method for comparing and interpreting land P-waves and converted S-wave channel sand bodies according to claim 1, characterized in that In the step S5, the method for performing full-area target horizon tracking and interpretation on the longitudinal wave profile and the converted shear wave profile according to the seismic response characteristics includes: According to the seismic response characteristics, carry out horizon interpretation at all well point positions in the work area; According to the horizon interpretation at the well point positions, carry out horizon interpretation of the cross-well seismic data in the work area; According to the horizon interpretation of the cross-well seismic data, carry out horizon picking for the 80×80, 40×40 and 20×20 grids in the work area, and finally interpolate and encrypt it into horizon picking for the 1×1 grid.

6. The method for comparing and interpreting land P-wave and converted S-wave channel sand bodies according to claim 1 is characterized in that In the step S6, the method for performing AV0 attribute analysis on the pre-stack time migration gather data of the longitudinal wave, extracting the intercept and gradient attributes, and obtaining the P-G attribute data volume by calculation includes: Wherein, P represents the intercept; G represents the gradient attribute; PG represents the P-G attribute data volume; v s represents the shear wave velocity; Δv S represents the shear wave velocity difference between the upper and lower strata; ρ represents the density; Δρ represents the density difference between the upper and lower strata.

7. A method for comparing and interpreting land P-waves and converted S-wave channel sand bodies according to claim 1, characterized in that In the step S7, the method for respectively calculating the trace integration attributes of the P-G attribute data volume and the converted shear wave data volume to obtain the P-G trace integration data volume and the converted shear wave trace integration data volume includes: Wherein, Spg represents the P-G trace integration data volume; PG(t) represents the P-G attribute data volume; wpg(t) is the wavelet of the P-G attribute data volume, and k is the reflection coefficient of the top interface; Wherein, Sps represents the converted shear wave trace integration data volume; PS(t) represents the converted shear wave data volume; wps(t) is the wavelet of the converted shear wave data volume.

8. A system for comparative interpretation of land P-wave and converted S-wave channel sand bodies, characterized in that The system includes: The first processing module is configured to perform pre-stack time migration processing on the longitudinal wave and converted shear wave seismic data to obtain the pre-stack time migration gather of the longitudinal wave, the longitudinal wave profile and the converted shear wave profile; The second processing module is configured to load the longitudinal wave profile data, the converted shear wave profile data and load the logging curve; load the geological stratification data of the target formation and sand body to obtain the geological stratification; The third processing module is configured to use the logging curve to carry out calibration of the synthetic records of the longitudinal wave and the converted shear wave; clarify the time position of the geological horizon at the adopted logging well point on the corresponding seismic profile; A fourth processing module is configured to clarify the seismic response characteristics of different target layers and sand bodies on the longitudinal wave section and the converted shear wave section according to the synthetic record calibration and the time position; A fifth processing module is configured to perform target layer tracking interpretation of the entire area on the longitudinal wave section and the converted shear wave section according to the seismic response characteristics; The sixth processing module is configured to perform AV0 attribute analysis on the P-wave prestack time migration gather data, extract intercept and gradient attributes, and obtain a PG attribute data volume by calculation; A seventh processing module is configured to perform channel integral attribute calculation on the PG attribute data volume and the converted shear wave data volume respectively to obtain a PG channel integral data volume and a converted shear wave channel integral data volume; An eighth processing module is configured to interpret the sand body layer according to the result of the target layer tracking interpretation and based on the converted shear wave channel integral data volume, and save the converted shear wave sand body layer picking result; A ninth processing module is configured to project the converted shear wave sand body layer position picking result onto a PG channel integral data volume, perform correction and supplement on the PG channel integral data volume, and save the PG channel integral data volume sand body layer position picking result; The tenth processing module is configured to project the sand body layer picking results of the PG channel integral data volume onto the longitudinal wave section, pick up the top and bottom interfaces of the sand body upward and downward along the crests and troughs on the longitudinal wave section according to the sand body layer of the PG channel integral data volume, and carry out the interpretation and characterization of the top and bottom interfaces of the sand body.

9. An electronic device, characterized in that, The electronic device includes a memory and a processor, wherein the memory stores a computer program, and when the processor executes the computer program, the steps of a land longitudinal wave and converted shear wave channel sand body comparative interpretation method described in any one of claims 1 to 7 are implemented.

10. 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 a processor, the steps of the land longitudinal wave and converted shear wave channel sand body comparative interpretation method described in any one of claims 1 to 7 are implemented.

Images