Time matching method and device for surface seismic P-S converted wave data

Abstract

The application relates to a time matching method and device for ground seismic P-wave converted wave data, which comprises the following steps: obtaining vertical seismic profile data, ground seismic P-wave data and ground seismic converted wave data; according to the vertical seismic profile data, solving a P-wave interval velocity curve and a S-wave interval velocity curve at a well point position; processing the vertical seismic profile data to separate out vertical seismic P-wave data and vertical seismic converted wave data; performing seismic inversion, analysis and calculation on the P-wave and S-wave interval velocity curves, the vertical seismic P-wave data and the vertical seismic converted wave data to obtain a relationship formula of P-wave velocity and S-wave velocity; combining the P-wave interval velocity curve at the well point position and ground seismic P-wave velocity data obtained by inverting the ground seismic P-wave data to further calculate ground seismic S-wave velocity; and finally converting the converted wave to P-wave two-way travel time according to the ground seismic P-wave velocity, the ground seismic S-wave velocity and the ground seismic converted wave data.

Description

Technical Field

[0001] This invention relates to the field of seismic exploration technology, and in particular to a time matching method and apparatus for ground seismic P-wave converted wave data. Background Technology

[0002] Multiwave seismic exploration can obtain wavefield information of various types, including P-waves, S-waves, and converted waves, and has unique advantages such as multiwave joint imaging, multiwave joint calibration, multiwave joint inversion, and multi-attribute fusion. It shows great promise in analyzing reservoir lithology (referring to the rock type, mineral composition, structure, physical properties, etc. of the reservoir) and predicting hydrocarbon content. However, because P-waves and S-waves have different propagation speeds (S-waves take longer to propagate than their corresponding P-waves), and in addition to the difference in propagation time, the reflection coefficients and polarities of P-wave and S-wave reflections also differ, resulting in different dynamic characteristics of the P-wave and S-wave wavefields. Therefore, to fully utilize the advantages of multiwave seismic exploration technology, the matching problem of multiwave seismic data must first be solved. Time matching of P-wave converted wave data is a crucial step in the interpretation of multiwave seismic exploration data.

[0003] Currently, time matching methods for ground seismic P-wave converted wave data are mainly based on stratigraphic matching. This involves manually interpreting P-wave and converted wave seismic exploration data separately, tracing the corresponding stratigraphic levels, and then using these interpreted stratigraphic levels to control the P-wave and converted wave data, generating time-shifted volumes for P-wave and converted wave matching. These time-shifted volumes are then applied to the P-wave and converted wave data to achieve time matching, i.e., unifying them to the two-way travel time of either the P-wave or the converted wave. Generally, the converted wave data is transformed from its time domain to the corresponding P-wave two-way travel time, laying the foundation for further joint interpretation of P-wave and converted wave data. This matching method heavily relies on manual interpretation experience, lacks sufficient control between stratigraphic levels, and is prone to waveform distortion, exhibiting significant uncertainty. Especially with data where stratigraphic reflection characteristics are unclear, high-precision matching cannot be achieved.

[0004] Therefore, how to propose a method to improve the temporal matching accuracy of multi-wave seismic data has become an urgent technical problem to be solved. Summary of the Invention

[0005] In view of this, the present invention aims to solve the problem of low time matching accuracy between P-wave and converted wave data, thereby reducing the difficulty of time matching between P-wave and converted wave and improving the time matching accuracy and efficiency between P-wave and converted wave.

[0006] Specifically, the present invention is achieved through the following technical solution:

[0007] According to a first aspect of the present invention, a time-matching method for ground seismic P-wave converted wave data is provided, comprising: acquiring vertical seismic profile data, ground seismic P-wave data, and ground seismic converted wave data; obtaining P-wave layer velocity curves and S-wave layer velocity curves at well locations based on the vertical seismic profile data; performing data processing on the vertical seismic profile data to obtain vertical seismic P-wave data and vertical seismic converted wave data; performing seismic inversion on the P-wave layer velocity curves, S-wave layer velocity curves, vertical seismic P-wave data, and vertical seismic converted wave data to obtain vertical seismic P-wave velocity and vertical seismic S-wave velocity; and performing time matching on the vertical seismic P-wave data. The relationship between P-wave velocity and vertical seismic S-wave velocity is obtained by analyzing the P-wave velocity and vertical seismic S-wave velocity. The ground seismic P-wave velocity is obtained by inversion from the P-wave layer velocity curve and ground seismic P-wave data. The ground seismic P-wave velocity is converted to the ground seismic S-wave velocity according to the relationship. Based on the ground seismic P-wave velocity and the ground seismic S-wave velocity, a well-to-surface fusion spatiotemporally varying P-wave and S-wave velocity ratio field is constructed. Based on the well-to-surface fusion spatiotemporally varying P-wave and S-wave velocity ratio field and the two-way travel time of the ground seismic converted wave data, the corresponding two-way travel time of the ground seismic P-wave data is determined, and the ground seismic converted wave data transformed to the P-wave two-way travel time is obtained.

[0008] Optionally, the steps of obtaining the P-wave velocity curves and S-wave velocity curves at the wellpoint location based on the vertical seismic profile data include: performing wavefield analysis on the vertical seismic profile data, and performing P-wave first arrival picking and S-wave first arrival picking to calculate the P-wave velocity and S-wave velocity at the wellpoint location; and obtaining the P-wave velocity curves and S-wave velocity curves based on the P-wave velocity and S-wave velocity.

[0009] Optionally, the steps of processing the vertical seismic profile data to obtain vertical seismic P-wave data and vertical seismic converted wave data include: performing wavefield separation and imaging processing on the vertical seismic profile data to obtain vertical seismic P-wave data and vertical seismic converted wave data.

[0010] Optionally, the steps of performing seismic inversion on the P-wave velocity curve, the S-wave velocity curve, the vertical seismic P-wave data, and the vertical seismic converted wave data to obtain the vertical seismic P-wave velocity and the vertical seismic S-wave velocity include: performing synthetic record calibration on the vertical seismic P-wave data and the vertical seismic converted wave data according to the P-wave velocity curve and the S-wave velocity curve, and performing layer interpretation processing, and obtaining the vertical seismic P-wave velocity and the vertical seismic S-wave velocity through seismic inversion.

[0011] Optionally, the relationship between the P-wave velocity and the S-wave velocity is as follows:

[0012] V S =0.58V P -149;

[0013] Among them, VS Expressed as transverse wave velocity, in units of m / s, V p It is expressed as longitudinal wave velocity, with the unit being m / s.

[0014] Optionally, the steps for obtaining the ground seismic P-wave velocity from the P-wave layer velocity curve and the ground seismic P-wave data include: performing synthetic recording and calibration of the ground seismic P-wave data based on the P-wave layer velocity curve, and layer interpretation processing, and obtaining the ground seismic P-wave velocity through seismic inversion.

[0015] Optionally, the step of converting the ground seismic P-wave velocity to the ground seismic S-wave velocity according to the relationship includes: performing amplitude processing on the ground seismic P-wave velocity according to the relationship between P-wave velocity and S-wave velocity to obtain the ground seismic S-wave velocity.

[0016] Optionally, the steps for constructing a well-to-ground integrated spatiotemporally varying P-wave and S-wave velocity ratio field based on the ground seismic P-wave velocity and the ground seismic S-wave velocity include: dividing the ground seismic P-wave velocity and the ground seismic S-wave velocity to obtain the velocity ratio of the ground seismic P-wave and S-wave; and performing spatial and temporal processing on the velocity ratio of the ground seismic P-wave and S-wave to construct a well-to-ground integrated spatiotemporally varying P-wave and S-wave velocity ratio field.

[0017] Optionally, the formula for calculating the two-way travel time of the corresponding ground seismic P-wave data is as follows:

[0018]

[0019] Among them, T pp T represents the two-way travel time of ground seismic P-wave data, in milliseconds (ms), where k represents the P-wave to S-wave velocity ratio. ps The two-way travel time for ground seismic converted wave data is expressed in milliseconds (ms).

[0020] According to a second aspect of the present invention, a time-matching device for ground seismic P-wave converted wave data is provided, comprising: an acquisition module for acquiring vertical seismic profile data, ground seismic P-wave data, and ground seismic converted wave data; a calculation module for calculating P-wave layer velocity curves and S-wave layer velocity curves at well locations based on the vertical seismic profile data; a processing module for processing the vertical seismic profile data to obtain vertical seismic P-wave data and vertical seismic converted wave data; an inversion module for performing seismic inversion on the P-wave layer velocity curves, S-wave layer velocity curves, vertical seismic P-wave data, and vertical seismic converted wave data to obtain vertical seismic P-wave velocity and vertical seismic S-wave velocity; and an analysis module for analyzing the vertical seismic profile data. The system analyzes the P-wave velocity and vertical S-wave velocity to obtain the relationship between them. The inversion module is used to invert the ground P-wave velocity based on the P-wave layer velocity curve and ground P-wave data. The conversion module converts the ground P-wave velocity to the ground S-wave velocity based on the relationship. The construction module constructs a well-to-surface fused spatiotemporally varying P-wave / S-wave velocity ratio field based on the ground P-wave and S-wave velocities. The determination module determines the corresponding two-way travel time of the ground P-wave data based on the well-to-surface fused spatiotemporally varying P-wave / S-wave velocity ratio field and the two-way travel time of the ground converted wave data, thus obtaining the ground converted wave data transformed to the P-wave two-way travel time.

[0021] According to a third aspect of the invention, a storage medium is provided having a computer program stored thereon, which, when executed by a processor, implements the steps of the time-matching method for ground seismic P-wave converted wave data in the first aspect of the invention or any possible implementation thereof.

[0022] According to a fourth aspect of the present invention, an electronic device is provided, including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the computer program to implement the steps of the time matching method for ground seismic P-wave converted wave data in the first aspect of the present invention or any possible implementation thereof.

[0023] The technical solution provided by this invention brings at least the following beneficial effects:

[0024] This invention performs multiple data processing steps on vertical seismic profile data, surface seismic P-wave data, and surface seismic converted wave data to obtain the velocity ratio fields of P-waves and S-waves. Then, by combining the velocity ratio fields of P-waves and S-waves with the two-way travel time of the surface seismic converted wave data, time-matched surface seismic converted wave data is obtained. This effectively achieves time matching between surface seismic P-wave and converted wave data, thereby improving the accuracy of time matching and further enhancing the accuracy of the joint interpretation results of surface seismic P-wave and converted wave data. Furthermore, the method of this invention has broad application prospects in multi-wave seismic exploration and is highly practical and efficient in production. Attached Figure Description

[0025] The accompanying drawings, which are incorporated in and form part of this specification, illustrate embodiments consistent with the invention and, together with the description, serve to explain the principles of the invention.

[0026] To more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the drawings used in the description of the embodiments or related technologies will be briefly introduced below. Obviously, those skilled in the art can obtain other drawings based on these drawings without creative effort.

[0027] Figure 1 A flowchart illustrating a time-matching method for ground seismic P-wave converted wave data provided in an embodiment of the present invention;

[0028] Figure 2 This is a schematic diagram of zero well source distance in the Z component provided in an embodiment of the present invention;

[0029] Figure 3 This is a schematic diagram of zero well source distance in the X component provided in an embodiment of the present invention;

[0030] Figure 4 This is a schematic diagram of zero well source distance in the Y component provided in an embodiment of the present invention;

[0031] Figure 5 A schematic diagram of ground seismic P-wave data provided in an embodiment of the present invention;

[0032] Figure 6 A schematic diagram of ground seismic converted wave data provided in an embodiment of the present invention;

[0033] Figure 7 This is a schematic diagram of two-dimensional vertical seismic data in the Z component provided in an embodiment of the present invention;

[0034] Figure 8 This is a schematic diagram of two-dimensional vertical seismic data in the X component provided in an embodiment of the present invention;

[0035] Figure 9 This is a schematic diagram of two-dimensional vertical seismic data in the Y component provided in an embodiment of the present invention;

[0036] Figure 10 A schematic diagram of vertical seismic P-wave data provided in an embodiment of the present invention;

[0037] Figure 11 A schematic diagram of vertical seismic converted wave data provided in an embodiment of the present invention;

[0038] Figure 12 This is a schematic diagram of vertical seismic P-wave velocity data provided in an embodiment of the present invention;

[0039] Figure 13 This is a schematic diagram of vertical seismic shear wave velocity data provided in an embodiment of the present invention;

[0040] Figure 14 A schematic diagram illustrating the analysis results of vertical seismic P-wave velocity and vertical seismic S-wave velocity provided in an embodiment of the present invention;

[0041] Figure 15 This is a schematic diagram of ground seismic P-wave velocity data provided in an embodiment of the present invention;

[0042] Figure 16 A schematic diagram of the spatiotemporally variable P-wave and S-wave velocity ratio field provided in an embodiment of the present invention;

[0043] Figure 17 A schematic diagram of ground seismic converted wave data before time matching provided in an embodiment of the present invention;

[0044] Figure 18 A schematic diagram of time-matched ground seismic converted wave data provided in an embodiment of the present invention;

[0045] Figure 19 A flowchart illustrating another time-matching method for ground seismic P-wave converted wave data provided in an embodiment of the present invention;

[0046] Figure 20 A block diagram of a time matching device for ground seismic P-wave converted wave data provided in an embodiment of the present invention;

[0047] Figure 21 This is a block diagram of an electronic device provided in an embodiment of the present invention. Detailed Implementation

[0048] To make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer, the technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0049] like Figure 1 As shown, this embodiment of the invention provides a time matching method for ground seismic P-wave converted wave data, including:

[0050] S101: Acquire vertical seismic profile data, surface seismic P-wave data, and surface seismic converted wave data;

[0051] S102: The P-wave velocity curves and S-wave velocity curves at the well locations are obtained based on the vertical seismic profile data.

[0052] S103: Process the vertical seismic profile data to obtain vertical seismic P-wave data and vertical seismic converted wave data;

[0053] S104: Perform seismic inversion on the P-wave velocity curve, S-wave velocity curve, vertical seismic P-wave data and vertical seismic converted wave data to obtain the vertical seismic P-wave velocity and vertical seismic S-wave velocity.

[0054] S105: Analyze the vertical seismic P-wave velocity and vertical seismic S-wave velocity to obtain the relationship between P-wave velocity and S-wave velocity;

[0055] S106: The ground seismic P-wave velocity is obtained by inversion from the P-wave layer velocity curve and ground seismic P-wave data;

[0056] S107: The ground seismic P-wave velocity is converted to the ground seismic S-wave velocity according to the formula;

[0057] S108: Construct a well-ground integrated spatiotemporally variable P-wave and S-wave velocity ratio field based on ground seismic P-wave velocity and ground seismic S-wave velocity;

[0058] S109: Determine the corresponding two-way travel time of the ground seismic P-wave data based on the spatiotemporally variable P-wave velocity ratio field of the well-ground fusion and the two-way travel time of the ground seismic converted wave data, and obtain the ground seismic converted wave data transformed into the P-wave two-way travel time.

[0059] The method in this embodiment is particularly suitable for cases with zero well source spacing, such as... Figure 2 , Figure 3 and Figure 4As shown, this embodiment illustrates the data at zero well-source distance, i.e., when the distance between the seismic source and the observation point is zero. First, vertical seismic profile data, surface P-wave data, and surface converted wave data are acquired, as shown below. Figure 5 and Figure 6 As shown, schematic diagrams of ground seismic P-wave data and ground seismic converted wave data provided in this embodiment are illustrated respectively. These seismic data are raw data collected by seismic instruments. Vertical seismic profile data refers to seismic wave data located underground, while ground seismic P-wave data and ground seismic converted wave data are seismic wave data received at the ground. Figure 7 , Figure 8 and Figure 9 The diagram shows the two-dimensional vertical seismic profile data provided in this embodiment. Based on the vertical seismic profile data, the P-wave velocity curves and S-wave velocity curves at the wellpoint locations are then determined. Further data processing of the vertical seismic profile data yields vertical seismic P-wave data and vertical seismic converted wave data, as shown below. Figure 10 As shown, a schematic diagram of the vertical seismic P-wave data provided in this embodiment is illustrated. Figure 11 The diagram illustrates the vertical seismic converted wave data provided in this embodiment. Further seismic inversion is performed based on the P-wave velocity curves, S-wave velocity curves, vertical seismic P-wave data, and vertical seismic converted wave data to obtain the vertical seismic P-wave velocity and vertical seismic S-wave velocity, as shown below. Figure 12 As shown, a schematic diagram of the vertical seismic P-wave velocity data provided in this embodiment is illustrated. Figure 13 The diagram shows a schematic representation of the vertical seismic shear wave velocity provided in this embodiment. Next, the vertical seismic P-wave velocity and vertical seismic shear wave velocity are analyzed to obtain the relationship between P-wave velocity and shear wave velocity, as shown below. Figure 14 The diagram illustrates the analysis results of vertical seismic P-wave velocity and vertical seismic S-wave velocity provided in this embodiment. Simultaneously, the ground seismic P-wave velocity is obtained by inverting the P-wave layer velocity curve and ground seismic P-wave data, as shown below. Figure 15 The diagram illustrates the data of ground seismic P-wave velocity provided in this embodiment. Using the aforementioned relationship, the ground seismic S-wave velocity is obtained from the P-wave velocity, thereby calculating the ground seismic P-wave velocity and the ground seismic S-wave velocity. This further determines the velocity ratio field between the P-wave and S-wave waves. Figure 16 The diagram illustrates the spatiotemporally varying P-wave and S-wave velocity ratio field provided in this embodiment. Finally, the corresponding two-way travel time of the P-wave data is calculated using the velocity ratio field of the surface P-waves and S-waves and the two-way travel time of the surface seismic converted wave data. This ensures that the P-wave and converted wave data are on the same two-way travel time, resulting in surface seismic converted wave data transformed to the P-wave two-way travel time. Figure 17As shown, this diagram illustrates the ground seismic converted wave data before time matching provided in this embodiment. Figure 18 The diagram illustrates the time-matched ground seismic converted wave data provided in this embodiment. This embodiment can improve the efficiency and accuracy of P-wave conversion and time matching, thereby improving the accuracy of interpretation results during multi-wave seismic exploration.

[0060] Understandably, multiwave seismic exploration is an exploration method used to analyze reservoir lithology and invert reservoir lithology and hydrocarbon content. Utilizing seismic waves to understand the subsurface medium can effectively reduce the potential for multiple interpretations during seismic interpretation, leading to more accurate results. However, in multiwave seismic exploration, different waves have different propagation speeds and reflection characteristics, making it difficult for technicians to perform joint interpretations based on multiwave seismic exploration data, thus reducing the accuracy of the interpretation results. The time matching method for P-wave and converted wave data provided in this embodiment is applicable to zero-source-source distance situations, i.e., when the distance between the seismic source and the observation point is zero. It can perform time matching of P-wave and converted wave data to obtain surface seismic converted wave data transformed into P-wave two-way travel time, solving the time matching problem of surface seismic P-wave converted wave data, improving matching accuracy, and further enhancing the accuracy of the joint interpretation results of surface seismic P-wave and converted wave data.

[0061] In addition, vertical seismic surveys involve placing the seismic source in a well and the geophone on the surface (or vice versa). These surveys primarily include zero-source-distance vertical seismic surveys, non-zero-source-distance vertical seismic surveys, two-dimensional vertical seismic surveys, and three-dimensional vertical seismic surveys. Generally, in zero-source-distance vertical seismic surveys, the horizontal distance between the seismic source and the wellhead is less than or equal to 150 meters, while in non-zero-source-distance vertical seismic surveys, the horizontal distance is greater than 150 meters. The vertical seismic profile data obtained in this embodiment can be either zero-source-distance or non-zero-source-distance data. However, for the matching method in this embodiment, the effect is particularly significant when the horizontal distance between the seismic source and the wellhead is less than or equal to 150 meters. Furthermore, the obtained vertical seismic profile data can be either two-dimensional or three-dimensional. The time matching method proposed in this invention has strong applicability, enabling time matching of two-dimensional or three-dimensional ground seismic P-wave converted wave data using vertical seismic data, thereby improving the accuracy of data interpretation results.

[0062] Optionally, the steps of obtaining the P-wave velocity curves and S-wave velocity curves at the wellpoint location based on the vertical seismic profile data include: performing wavefield analysis on the vertical seismic profile data, and performing P-wave first arrival picking and S-wave first arrival picking to calculate the P-wave velocity and S-wave velocity at the wellpoint location; and obtaining the P-wave velocity curves and S-wave velocity curves based on the P-wave velocity and S-wave velocity.

[0063] In this embodiment, the process of obtaining the P-wave velocity curve and S-wave velocity curve at the well point location specifically includes: performing wavefield analysis on the vertical seismic profile data. Wavefield analysis can remove noise from the data and improve the accuracy of data processing. Then, the first arrivals of the P-wave and S-wave are picked up, that is, the first arrival data of the seismic P-wave and S-wave are obtained. The P-wave velocity and S-wave velocity at the well point location are further obtained. Then, the P-wave velocity curve and S-wave velocity curve are plotted based on the P-wave velocity and S-wave velocity at the well point location.

[0064] Optionally, the steps of processing the vertical seismic profile data to obtain vertical seismic P-wave data and vertical seismic converted wave data include: performing wavefield separation and imaging processing on the vertical seismic profile data to obtain vertical seismic P-wave data and vertical seismic converted wave data.

[0065] In this embodiment, the process of obtaining vertical seismic P-wave data and vertical seismic converted wave data specifically includes: performing wavefield separation and imaging processing on the vertical seismic profile data, that is, performing wavefield separation on the raw vertical seismic data, separating the P-wave data and converted wave data, and performing data imaging processing to obtain vertical seismic P-wave data and vertical seismic converted wave data.

[0066] Optionally, the steps of performing seismic inversion on the P-wave velocity curve, the S-wave velocity curve, the vertical seismic P-wave data, and the vertical seismic converted wave data to obtain the vertical seismic P-wave velocity and the vertical seismic S-wave velocity include: performing synthetic record calibration on the vertical seismic P-wave data and the vertical seismic converted wave data according to the P-wave velocity curve and the S-wave velocity curve, and performing layer interpretation processing, and obtaining the vertical seismic P-wave velocity and the vertical seismic S-wave velocity through seismic inversion.

[0067] In this embodiment, the process of obtaining the vertical seismic P-wave velocity and vertical seismic S-wave velocity through seismic inversion includes: synthesizing and calibrating the vertical seismic P-wave data and vertical seismic converted wave data according to the P-wave layer velocity curve and the S-wave layer velocity curve; recording and calibrating the vertical seismic P-wave data and vertical seismic converted wave data in layers with different velocities to improve the accuracy and reliability of data processing, and also to provide a basis for subsequent data analysis and processing; then performing layer interpretation processing, which can be understood as tracing the distribution and orientation of each layer on the seismic profile and drawing the layer contour map; combining geological theory and actual geological conditions to interpret and infer the geological structure; and further obtaining the vertical seismic P-wave velocity and vertical seismic S-wave velocity through seismic inversion, thereby improving the accuracy of the calculation of the vertical seismic P-wave and S-wave velocities.

[0068] Optionally, the relationship between the P-wave velocity and the S-wave velocity is as follows:

[0069] V S =0.58V P -149;

[0070] Among them, V S Expressed as transverse wave velocity, in units of m / s, V p It is expressed as longitudinal wave velocity, with the unit being m / s.

[0071] In this embodiment, V S Expressed as transverse wave velocity, V p The longitudinal wave velocity is expressed as the longitudinal wave velocity, and both are in m / s. The transverse wave velocity can be further calculated using the above relationship and the longitudinal wave velocity that has already been obtained.

[0072] Optionally, the steps for obtaining the ground seismic P-wave velocity from the P-wave layer velocity curve and the ground seismic P-wave data include: performing synthetic recording and calibration of the ground seismic P-wave data based on the P-wave layer velocity curve, and layer interpretation processing, and obtaining the ground seismic P-wave velocity through seismic inversion.

[0073] In this embodiment, the process of inverting to obtain the ground seismic P-wave velocity specifically includes: performing synthetic recording and calibration of the ground seismic P-wave data based on the P-wave layer velocity curve, as well as layer interpretation processing, and using seismic inversion processing of the ground seismic P-wave data to obtain the ground seismic P-wave velocity, thereby improving the calculation accuracy of the ground seismic P-wave velocity.

[0074] Optionally, the step of converting the ground seismic P-wave velocity to the ground seismic S-wave velocity according to the relationship includes: performing amplitude processing on the ground seismic P-wave velocity according to the relationship between P-wave velocity and S-wave velocity to obtain the ground seismic S-wave velocity.

[0075] In this embodiment, the ground seismic P-wave velocity is processed by amplitude adjustment according to the relationship between P-wave velocity and S-wave velocity to obtain the ground seismic S-wave velocity. Amplitude adjustment involves calculating the ground seismic S-wave velocity using the relationship between P-wave velocity and S-wave velocity.

[0076] Optionally, the steps for constructing a well-to-ground integrated spatiotemporally varying P-wave and S-wave velocity ratio field based on the ground seismic P-wave velocity and the ground seismic S-wave velocity include: dividing the ground seismic P-wave velocity and the ground seismic S-wave velocity to obtain the velocity ratio of the ground seismic P-wave and S-wave; and performing spatial and temporal processing on the velocity ratio of the ground seismic P-wave and S-wave to construct a well-to-ground integrated spatiotemporally varying P-wave and S-wave velocity ratio field.

[0077] In this embodiment, after obtaining the surface seismic P-wave velocity and surface seismic S-wave velocity, the P-wave velocity and S-wave velocity are divided to calculate the velocity ratio of the P-wave to the S-wave. Continuous P-wave to S-wave velocity ratios form a velocity ratio field. Then, spatial and temporal processing is performed on the surface seismic P-wave and S-wave velocity ratios to construct a well-ground integrated spatiotemporally variable P-wave and S-wave velocity ratio field. That is, spatial and temporal processing is performed on the velocity ratio field to take into account the spatiotemporal variations of surface seismic P-waves and S-waves. The constructed well-ground integrated spatiotemporally variable P-wave and S-wave velocity ratio field can improve the accuracy of seismic exploration.

[0078] Optionally, the formula for calculating the two-way travel time of the corresponding ground seismic P-wave data is as follows:

[0079]

[0080] Among them, T pp T represents the two-way travel time of ground seismic P-wave data, in milliseconds (ms), where k represents the P-wave to S-wave velocity ratio. ps The two-way travel time for ground seismic converted wave data is expressed in milliseconds (ms).

[0081] In this embodiment, the two-way travel time of the ground seismic P-wave data is calculated according to the above calculation formula. The P-wave / S-wave velocity ratio has been obtained in the previous steps. The ground seismic converted wave data contains two-way travel time data. That is, the two-way travel time of the ground seismic P-wave data can be calculated based on the two-way travel time of the ground seismic converted wave data, thereby obtaining the ground seismic converted wave data transformed into the P-wave two-way travel time.

[0082] In a specific application, such as Figure 19 As shown, this embodiment of the invention provides another time-matching method for ground seismic P-wave converted wave data, including:

[0083] S201: Acquire vertical seismic profile data, surface seismic P-wave data, and surface seismic converted wave data;

[0084] S202: Perform wavefield analysis on vertical seismic profile data, and pick up the first arrival of P-waves and the first arrival of S-waves to calculate the P-wave and S-wave velocities at the well point location.

[0085] S203: Based on the P-wave velocity and S-wave velocity, obtain the P-wave velocity curve and the S-wave velocity curve;

[0086] S204: Perform wavefield separation and imaging processing on vertical seismic profile data to obtain vertical seismic P-wave data and vertical seismic converted wave data;

[0087] S205: Based on the P-wave layer velocity curve and the S-wave layer velocity curve, the vertical seismic P-wave data and vertical seismic converted wave data are synthesized, calibrated, and interpreted by layer. The vertical seismic P-wave velocity and vertical seismic S-wave velocity are obtained through seismic inversion.

[0088] S206: Analyze the vertical seismic P-wave velocity and vertical seismic S-wave velocity to obtain the relationship between P-wave velocity and S-wave velocity;

[0089] S207: Based on the P-wave layer velocity curve, the ground seismic P-wave data is synthesized, recorded, calibrated, and interpreted by layer. The ground seismic P-wave velocity is obtained through seismic inversion.

[0090] S208: The ground earthquake P-wave velocity is obtained by amplitude processing based on the relationship between P-wave velocity and S-wave velocity;

[0091] S209: Calculate the velocity ratio of ground earthquake P-waves to ground earthquake S-waves by dividing the ground earthquake P-wave velocity and the ground earthquake S-wave velocity;

[0092] S210: Spatial and temporal processing of the velocity ratio of P-waves and S-waves in ground earthquakes to construct a well-ground integrated spatiotemporally variable P-wave and S-wave velocity ratio field;

[0093] S211: Determine the corresponding two-way travel time of the ground seismic P-wave data based on the spatiotemporally variable P-wave velocity ratio field of the well-ground fusion and the two-way travel time of the ground seismic converted wave data, and obtain the ground seismic converted wave data transformed into the P-wave two-way travel time.

[0094] The time matching method for ground seismic P-wave converted wave data proposed in this embodiment is characterized by strong operability, obvious effect, and high practical value in production. It effectively solves the problem of time matching between ground seismic P-wave and converted wave data, which not only improves the matching accuracy but also contributes to the quality and accuracy of the joint interpretation results of ground seismic P-wave and converted wave data.

[0095] like Figure 20As shown, according to an embodiment of the second aspect of the present invention, a time-matching device 10 for ground seismic P-wave converted wave data is provided, comprising: an acquisition module 11 for acquiring vertical seismic profile data, ground seismic P-wave data, and ground seismic converted wave data; a calculation module 12 for calculating the P-wave layer velocity curve and the S-wave layer velocity curve at the well point location based on the vertical seismic profile data; a processing module 13 for processing the vertical seismic profile data to obtain vertical seismic P-wave data and vertical seismic converted wave data; an inversion module 14 for performing seismic inversion on the P-wave layer velocity curve, the S-wave layer velocity curve, the vertical seismic P-wave data, and the vertical seismic converted wave data to obtain the vertical seismic P-wave velocity and the vertical seismic S-wave velocity; and an analysis module 15. The system is used to analyze vertical seismic P-wave velocity and vertical seismic S-wave velocity to obtain the relationship between P-wave velocity and S-wave velocity; the inversion module 14 is also used to invert the ground seismic P-wave velocity based on the P-wave layer velocity curve and the ground seismic P-wave data; the conversion module 16 is used to convert the ground seismic P-wave velocity to the ground seismic S-wave velocity based on the relationship; the construction module 17 is used to construct the well-to-surface fused spatiotemporally varying P-wave and S-wave velocity ratio field based on the ground seismic P-wave velocity and the ground seismic S-wave velocity; the determination module 18 is used to determine the corresponding two-way travel time of the ground seismic P-wave data based on the well-to-surface fused spatiotemporally varying P-wave and S-wave velocity ratio field and the two-way travel time of the ground seismic converted wave data, and obtain the ground seismic converted wave data transformed to the P-wave two-way travel time.

[0096] In this embodiment, the time-matching device 10 for ground seismic P-wave converted wave data first acquires vertical seismic profile data, ground seismic P-wave data, and ground seismic converted wave data through the acquisition module 11. These seismic data are raw data collected by seismic instruments. The vertical seismic profile data refers to seismic wave data received underground, while the ground seismic P-wave data and ground seismic converted wave data are seismic wave data received at the ground. Then, the calculation module 12 determines the P-wave velocity curve and S-wave velocity curve at the well point location based on the vertical seismic profile data. Next, the processing module 13 processes the vertical seismic profile data to obtain the vertical seismic P-wave data and vertical seismic converted wave data. Finally, the inversion module 14 performs seismic inversion based on the P-wave velocity curve, S-wave velocity curve, vertical seismic P-wave data, and vertical seismic converted wave data to obtain the vertical seismic P-wave velocity and vertical seismic S-wave velocity. Next, the vertical seismic P-wave velocity and vertical seismic S-wave velocity are analyzed by the analysis module 15 to obtain the relationship between P-wave velocity and S-wave velocity. Simultaneously, the surface seismic P-wave velocity is obtained by inverting the P-wave layer velocity curve and surface seismic P-wave data using the inversion module 14. The conversion module 16 then uses the above relationship to obtain the surface seismic S-wave velocity based on the surface seismic P-wave velocity, thereby calculating the surface seismic P-wave velocity and surface seismic S-wave velocity. Further, the velocity ratio field of the surface P-wave and S-wave is determined by the construction module 17. Finally, the determination module 18 calculates the corresponding two-way travel time of the surface seismic P-wave data based on the velocity ratio field of the surface P-wave and S-wave and the two-way travel time of the surface seismic converted wave data. This ensures that the surface seismic P-wave and converted wave data are on the same two-way travel time, i.e., the surface seismic converted wave data transformed to the P-wave two-way travel time. This embodiment can improve the efficiency and accuracy of P-wave conversion and wave conversion time matching, thereby improving the accuracy of interpretation results during multi-wave seismic exploration.

[0097] Optionally, the calculation module 12 is specifically used for: performing wavefield analysis on the vertical seismic profile data, picking up the first arrival of P-waves and the first arrival of S-waves, calculating the P-wave layer velocity and the S-wave layer velocity at the well point location; and obtaining the P-wave layer velocity curve and the S-wave layer velocity curve based on the P-wave layer velocity and the S-wave layer velocity.

[0098] Optionally, the processing module 13 is specifically used to: perform wavefield separation and imaging processing on the vertical seismic profile data to obtain vertical seismic P-wave data and vertical seismic converted wave data.

[0099] Optionally, the inversion module 14 is specifically used to: perform synthetic recording and calibration of vertical seismic P-wave data and vertical seismic converted wave data based on the P-wave layer velocity curve and the S-wave layer velocity curve, as well as layer interpretation processing, and obtain the vertical seismic P-wave velocity and vertical seismic S-wave velocity through seismic inversion.

[0100] Optionally, the inversion module 14 is also specifically used for: performing synthetic recording and calibration of ground seismic P-wave data based on the P-wave layer velocity curve, and layer interpretation processing, and obtaining the ground seismic P-wave velocity through seismic inversion.

[0101] Optionally, the conversion module 16 is specifically used to: perform amplitude processing on the ground seismic P-wave velocity according to the relationship between P-wave velocity and S-wave velocity to obtain the ground seismic S-wave velocity.

[0102] Optionally, the construction module 17 is specifically used to: calculate the ratio of ground earthquake P-wave velocity to ground earthquake S-wave velocity by dividing the ground earthquake P-wave velocity and the ground earthquake S-wave velocity; and to perform spatial and temporal processing on the ground earthquake P-wave and S-wave velocity ratio to construct a well-ground integrated spatiotemporally variable P-wave and S-wave velocity ratio field.

[0103] According to a third aspect of the invention, a storage medium is provided having a computer program stored thereon, which, when executed by a processor, implements the steps of the time-matching method for ground seismic P-wave converted wave data in the first aspect of the invention or any possible implementation thereof.

[0104] like Figure 21 As shown, according to a fourth aspect of the present invention, an electronic device 20 is provided, including a memory 21, a processor 22, and a computer program stored in the memory 21 and executable on the processor 22. When the processor executes the computer program, it implements the steps of the time matching method for ground seismic P-wave converted wave data in the first aspect of the present invention or any possible implementation thereof.

[0105] The specific implementation process of the functions and roles of each unit in the above device can be found in the implementation process of the corresponding steps in the above method, and will not be repeated here.

[0106] It should be noted that, in this document, relational terms such as "first" and "second" are used merely to distinguish one entity or operation from another, and do not necessarily require or imply any such actual relationship or order between these entities or operations. Furthermore, the terms "comprising," "including," or any other variations thereof are intended to cover non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements includes not only those elements but also other elements not expressly listed, or elements inherent to such a process, method, article, or apparatus. Unless otherwise specified, an element defined by the phrase "comprising one..." does not exclude the presence of other identical elements in the process, method, article, or apparatus that includes the element.

[0107] The above are merely specific embodiments of the present invention, enabling those skilled in the art to understand or implement the invention. Various modifications to these embodiments will be readily apparent to those skilled in the art, and the general principles defined herein may be implemented in other embodiments without departing from the spirit or scope of the invention. Therefore, the present invention is not to be limited to the embodiments shown herein, but is to be accorded the widest scope consistent with the principles and novel features claimed herein.

Claims

1. A time-matching method for ground seismic P-wave converted wave data, characterized in that, include: Acquire vertical seismic profile data, surface seismic P-wave data, and surface seismic converted wave data; Based on the vertical seismic profile data, the P-wave velocity curves and S-wave velocity curves at the well locations were obtained. The vertical seismic profile data is processed to obtain vertical seismic P-wave data and vertical seismic converted wave data; Seismic inversion is performed on the P-wave velocity curve, the S-wave velocity curve, the vertical seismic P-wave data, and the vertical seismic converted wave data to obtain the vertical seismic P-wave velocity and the vertical seismic S-wave velocity. The vertical seismic P-wave velocity and the vertical seismic S-wave velocity were analyzed to obtain the relationship between the P-wave velocity and the S-wave velocity. The ground seismic P-wave velocity is obtained by inversion from the P-wave layer velocity curve and the ground seismic P-wave data. The ground seismic P-wave velocity is converted into the ground seismic S-wave velocity according to the aforementioned relationship; Based on the ground seismic P-wave velocity and the ground seismic S-wave velocity, a well-ground integrated spatiotemporally variable P-wave and S-wave velocity ratio field is constructed. Based on the spatiotemporally variable P-wave and S-wave velocity ratio field of the well-ground fusion and the two-way travel time of the ground seismic converted wave data, the two-way travel time of the corresponding ground seismic P-wave data is determined, and the ground seismic converted wave data transformed into P-wave two-way travel time is obtained.

2. The time matching method for ground seismic P-wave converted wave data according to claim 1, characterized in that, The steps of obtaining the P-wave velocity curves and S-wave velocity curves at the wellpoint locations based on the vertical seismic profile data include: Wavefield analysis was performed on the vertical seismic profile data, and P-wave first arrival and S-wave first arrival were picked up to calculate the P-wave layer velocity and S-wave layer velocity at the well point location. Based on the P-wave velocity and the S-wave velocity, the P-wave velocity curve and the S-wave velocity curve are obtained.

3. The time matching method for ground seismic P-wave converted wave data according to claim 1, characterized in that, The step of processing the vertical seismic profile data to obtain vertical seismic P-wave data and vertical seismic converted wave data includes: Wavefield separation and imaging processing are performed on the vertical seismic profile data to obtain the vertical seismic P-wave data and the vertical seismic converted wave data.

4. The time matching method for ground seismic P-wave converted wave data according to claim 1, characterized in that, The step of performing seismic inversion on the P-wave velocity curve, the S-wave velocity curve, the vertical seismic P-wave data, and the vertical seismic converted wave data to obtain the vertical seismic P-wave velocity and the vertical seismic S-wave velocity includes: Based on the P-wave velocity curve and the S-wave velocity curve, the vertical seismic P-wave data and the vertical seismic converted wave data are synthesized, recorded, calibrated, and subjected to layer interpretation processing. The vertical seismic P-wave velocity and the vertical seismic S-wave velocity are obtained through seismic inversion.

5. The time matching method for ground seismic P-wave converted wave data according to claim 1, characterized in that, The relationship between the longitudinal wave velocity and the transverse wave velocity is as follows: V S =0.58V P -149; Among them, V S Expressed as transverse wave velocity, in units of m / s, V p It is expressed as longitudinal wave velocity, with the unit being m / s.

6. The time matching method for ground seismic P-wave converted wave data according to claim 1, characterized in that, The step of inverting the ground seismic P-wave velocity based on the P-wave layer velocity curve and the ground seismic P-wave data includes: Based on the P-wave velocity curve, the ground seismic P-wave data is synthesized, recorded, calibrated, and interpreted according to the stratigraphic layers. The ground seismic P-wave velocity is then obtained through seismic inversion.

7. The time matching method for ground seismic P-wave converted wave data according to claim 1, characterized in that, The step of converting the ground seismic P-wave velocity to the ground seismic S-wave velocity according to the aforementioned relationship includes: The ground seismic P-wave velocity is obtained by amplitude processing based on the relationship between the P-wave velocity and the S-wave velocity.

8. The time matching method for ground seismic P-wave converted wave data according to claim 1, characterized in that, The step of constructing a well-to-surface integrated spatiotemporally varying P-wave and S-wave velocity ratio field based on the ground seismic P-wave velocity and the ground seismic S-wave velocity includes: The velocity ratio of ground earthquake P-waves to S-waves is calculated by dividing the ground earthquake P-wave velocity by the ground earthquake S-wave velocity. Spatial and temporal processing is performed on the velocity ratio of the ground seismic P-waves and S-waves to construct the well-ground integrated spatiotemporally variable P-wave and S-wave velocity ratio field.

9. The time matching method for ground seismic P-wave converted wave data according to claim 1, characterized in that, The formula for calculating the two-way travel time of the corresponding ground seismic P-wave data is as follows: Among them, T pp T represents the two-way travel time of ground seismic P-wave data, k represents the P-wave / S-wave velocity ratio, and T represents the two-way travel time of ground seismic P-wave data. ps This represents the two-way travel time of ground-based converted wave data.

10. A time matching device for ground seismic P-wave converted wave data, characterized in that, include: The acquisition module is used to acquire vertical seismic profile data, ground seismic P-wave data, and ground seismic converted wave data; The calculation module is used to obtain the P-wave velocity curve and S-wave velocity curve at the well point location based on the vertical seismic profile data. The processing module is used to process the vertical seismic profile data to obtain vertical seismic P-wave data and vertical seismic converted wave data. The inversion module is used to perform seismic inversion on the P-wave velocity curve, the S-wave velocity curve, the vertical seismic P-wave data and the vertical seismic converted wave data to obtain the vertical seismic P-wave velocity and the vertical seismic S-wave velocity. The analysis module is used to analyze the vertical seismic P-wave velocity and the vertical seismic S-wave velocity to obtain the relationship between the P-wave velocity and the S-wave velocity. The inversion module is also used to invert the ground seismic P-wave velocity based on the P-wave layer velocity curve and the ground seismic P-wave data; A conversion module is used to convert the ground seismic P-wave velocity into the ground seismic S-wave velocity according to the relationship; The construction module is used to construct a well-ground fusion spatiotemporally variable P-wave and S-wave velocity ratio field based on the ground seismic P-wave velocity and the ground seismic S-wave velocity. The determination module is used to determine the corresponding two-way travel time of the ground seismic P-wave data based on the spatiotemporally variable P-wave velocity ratio field of the well-ground fusion and the two-way travel time of the ground seismic converted wave data, so as to obtain the ground seismic converted wave data transformed into the P-wave two-way travel time.

11. A storage medium having a computer program stored thereon, characterized in that, When the computer program is executed by a processor, it implements the steps of the method according to any one of claims 1 to 9.

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

Images