Shear wave vector inversion modeling method and device

Abstract

The present application provides a shear wave vector inversion modeling method and device, which belongs to the field of shear wave static correction technology. This technical solution obtains a primary inversion shear wave velocity model by adopting primary constrained inversion. According to the model, shear wave velocity control points can be selected at different positions, so that the shear wave model can be effectively constrained. Based on the shear wave velocity control points, the corresponding common center point domain data is obtained, and secondary constrained inversion is performed to obtain a secondary inversion velocity model, which has high accuracy. The present application is applicable to shear wave seismic exploration projects. By picking up the first arrival of shear waves on StRt component data and applying the secondary constrained inversion modeling method, it better adapts to the rapid changes in shear wave velocity, improves the accuracy of the shear wave model, eliminates the low-amplitude anomaly caused by the surface layer, and lays a solid model foundation for the solution of shear wave static correction and pre-stack depth migration imaging.

Description

Technical Field

[0001] The present application relates to the technical field of shear wave static correction, and in particular to a shear wave vector inversion modeling method and device. Background Art

[0002] In natural gas extraction technology, it is necessary to process the seismic data of the formation to obtain the structure of the gas-bearing area in the formation. The seismic data includes longitudinal wave seismic data and shear wave seismic data. Due to the gas content in the formation, the imaging accuracy of longitudinal wave seismic data is low, and the changes in the surface shear wave structure in the shear wave seismic data will cause anomalies. The anomaly is mixed with the low-amplitude structure and is difficult to distinguish. Therefore, research on shear wave vector inversion modeling has been carried out.

[0003] The commonly implemented shear-wave vector inversion modeling method currently uses three types of source excitation: P waves, SH waves, and SV waves. A three-component geophone receives nine-component three-dimensional vector seismic data to perform shear-wave vector inversion modeling. Affected by the three-dimensional observation position and the direction of shear-wave propagation, the shear-wave energy received by the three-component geophone constantly changes, making it difficult to pick up the first arrival of shear waves. At the same time, the surface shear wave velocity changes more dramatically than that of longitudinal waves, making it difficult to accurately obtain constraint information and control the shear wave variation characteristics, which affects the accuracy of the first arrival tomographic inversion. Summary of the Invention

[0004] The present invention provides a method and apparatus for shear wave vector inversion modeling, which can improve the accuracy of shear wave models and eliminate low-amplitude anomalies caused by surface layers, laying a solid model foundation for solving shear wave static correction and prestack depth migration imaging. The technical solution is as follows:

[0005] In one aspect, a shear wave vector inversion modeling method is provided, the method comprising:

[0006] The initial SH wave component seismic data in the area are rotated toward the direction of the shot-detection line;

[0007] Based on the rotationally transformed SH wave component seismic data, the StRt component seismic data perpendicular to the shot-detection line is obtained;

[0008] Based on the StRt component seismic data, the first arrival data of the StRt component is obtained;

[0009] Based on the first arrival data of the StRt component and the first constraint model, the first inversion shear wave velocity model is obtained through the first constraint inversion;

[0010] Based on the initial inversion shear wave velocity model, a corresponding number of shear wave velocity control points are obtained at different positions;

[0011] Based on the StRt component seismic data of the shear wave velocity control point, the shot point information of the corresponding position is extracted;

[0012] Based on the shot point information, corresponding common center point domain data is obtained;

[0013] Perform refraction layer interpretation in the common center point domain of the shot point to obtain the velocity and thickness of the shear wave velocity control point;

[0014] Based on the first arrival data of the StRt component, the primary constraint model, and the velocity and thickness of the shear wave velocity control points, the secondary inversion velocity model is obtained through secondary constraint inversion.

[0015] In one possible implementation, the initial SH wave component seismic data has three components.

[0016] In one possible implementation, based on the initial inversion of the shear wave velocity model, a corresponding number of shear wave velocity control points are obtained at different positions, including:

[0017] Based on the initial inversion of the shear wave velocity model, the velocity interface fluctuations and velocity variation characteristics are obtained;

[0018] Based on the characteristics of velocity interface fluctuation and velocity variation, a corresponding number of shear wave velocity control points are obtained at different positions.

[0019] In one possible implementation, based on the velocity interface fluctuation and velocity variation characteristics, a corresponding number of shear wave velocity control points are obtained at different positions, including:

[0020] Based on the characteristics of velocity interface fluctuation and velocity change, the corresponding area is extracted at different positions;

[0021] Based on the area of the region, a corresponding number of shear wave velocity control points are obtained.

[0022] In one possible implementation, obtaining corresponding common center point domain data based on the shot point information includes:

[0023] Based on the shot point information, the corresponding near-arrangement first arrival data is obtained;

[0024] Get the common center point domain data corresponding to the nearly arranged first arrival data.

[0025] In one aspect, a shear wave vector inversion modeling apparatus is provided, the apparatus comprising:

[0026] The rotation module is used to rotate the initial SH wave component seismic data in the area toward the direction of the shot-detection line;

[0027] A component acquisition module is used to acquire StRt component seismic data perpendicular to the shot-detection line based on the SH wave component seismic data after rotation transformation;

[0028] A first arrival acquisition module is used to acquire first arrival data of the StRt component based on the StRt component seismic data;

[0029] The primary inversion module is used to obtain the primary inversion shear wave velocity model based on the first arrival data of the StRt component and the primary constraint model through the primary constraint inversion;

[0030] A control point acquisition module is used to obtain a corresponding number of shear wave velocity control points at different positions based on the initial inversion shear wave velocity model;

[0031] The shot point acquisition module is used to extract the shot point information of the corresponding position based on the StRt component seismic data of the shear wave velocity control point;

[0032] A constraint condition module is used to obtain corresponding common center point domain data based on the shot point information;

[0033] The constraint condition module is also used to perform refraction layer interpretation in the common center point domain of the shot point to obtain the velocity and thickness of the shear wave velocity control point;

[0034] The secondary inversion module is used to obtain the secondary inversion velocity model through secondary constrained inversion based on the first arrival data of the StRt component, the primary constraint model, and the velocity and thickness of the shear wave velocity control point.

[0035] In one possible implementation, the initial SH wave component seismic data has three components.

[0036] In one possible implementation, the control point acquisition module is configured to:

[0037] Based on the initial inversion of the shear wave velocity model, the velocity interface fluctuations and velocity variation characteristics are obtained;

[0038] Based on the characteristics of velocity interface fluctuation and velocity variation, a corresponding number of shear wave velocity control points are obtained at different positions.

[0039] In one possible implementation, the control point acquisition module is configured to:

[0040] Based on the characteristics of velocity interface fluctuation and velocity change, the corresponding area is extracted at different positions;

[0041] Based on the area of the region, a corresponding number of shear wave velocity control points are obtained.

[0042] In one possible implementation, the constraint condition module is used to:

[0043] Based on the shot point information, the corresponding near-arrangement first arrival data is obtained;

[0044] Get the common center point domain data corresponding to the nearly arranged first arrival data.

[0045] The technical solution provided in the embodiment of the present application obtains a primary inversion shear wave velocity model by adopting primary constrained inversion. According to the model, a corresponding number of shear wave velocity control points can be selected at different positions. Compared with blind selection of shear wave velocity control points, the selection of the positions of the above-mentioned shear wave velocity control points can effectively constrain the shear wave model. Based on the shear wave velocity control points, the common center point domain data of the corresponding shot point is selected to perform secondary constrained inversion, thereby obtaining a secondary inversion velocity model with high accuracy. The present application is applicable to shear wave seismic exploration projects. By picking up the first arrival of shear waves on StRt component data and applying the secondary constrained inversion modeling method, it better adapts to the rapid changes in shear wave velocity, improves the accuracy of the shear wave model, eliminates the low-amplitude anomaly caused by the surface layer, and lays a solid model foundation for solving shear wave static correction and pre-stack depth migration imaging. BRIEF DESCRIPTION OF THE DRAWINGS

[0046] In order to more clearly illustrate the technical solutions in the embodiments of the present application, the following briefly introduces the drawings required for use in the description of the embodiments. Obviously, the drawings described below are only some embodiments of the present application. For ordinary technicians in this field, other drawings can be obtained based on these drawings without any creative work.

[0047] Figure 1 This is a flow chart of a shear wave vector inversion modeling method provided in an embodiment of the present application;

[0048] Figure 2 This is a flow chart of a shear wave vector inversion modeling method provided in an embodiment of the present application;

[0049] Figure 3 This is a schematic diagram of a rotation principle provided by an embodiment of the present application;

[0050] Figure 4 This is a schematic diagram of an initial inversion result analysis provided in an embodiment of the present application;

[0051] Figure 5 This is a schematic diagram of a shot point and its near arrangement provided in an embodiment of the present application;

[0052] Figure 6 This is a schematic diagram of a near arrangement provided in an embodiment of the present application;

[0053] Figure 7 This is a first arrival time distance diagram provided by an embodiment of the present application;

[0054] Figure 8 Schematic diagram of the structure of a shear wave vector inversion modeling device provided in an embodiment of the present application;

[0055] Figure 9 It is a structural diagram of a computer device provided in an embodiment of the present application. DETAILED DESCRIPTION

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

[0057] In this application, unless otherwise specified or limited, terms such as "mounted," "connected," "connect," and "fixed" should be understood in a broad sense. For example, they may refer to fixed, detachable, or integral connections; mechanical or electrical connections; direct or indirect connections through an intermediate medium; and internal communication between two components. Those skilled in the art will understand the specific meanings of these terms in this application based on specific circumstances.

[0058] Figure 1 This is a flow chart of a shear wave vector inversion modeling method provided in an embodiment of the present application. Figure 1 , the method can be applied to a computer device, the method comprising:

[0059] 101. The initial SH wave component seismic data in the area are rotated in the direction of the shot-detection line.

[0060] 102. Based on the rotationally transformed SH wave component seismic data, obtain the StRt component seismic data perpendicular to the shot-detection line.

[0061] 103. Based on the StRt component seismic data, obtain the first arrival data of the StRt component.

[0062] 104. Based on the first arrival data of the StRt component and the first constraint model, the first inversion shear wave velocity model is obtained through the first constraint inversion.

[0063] 105. Based on the initial inversion of the shear wave velocity model, obtain a corresponding number of shear wave velocity control points at different locations.

[0064] 106. Based on the StRt component seismic data of the shear wave velocity control point, the shot point information of the corresponding position is extracted.

[0065] 107. Based on the shot point information, corresponding common center point domain data is obtained.

[0066] 108. Perform refraction layering interpretation in the common center point domain of the shot point to obtain the velocity and thickness of the shear wave velocity control point.

[0067] 109. Based on the first arrival data of the StRt component, the primary constraint model, and the velocity and thickness of the shear wave velocity control point, the secondary inversion velocity model is obtained through secondary constraint inversion.

[0068] The method provided in the embodiment of the present application obtains a primary inversion shear wave velocity model by adopting primary constrained inversion. According to the model, a corresponding number of shear wave velocity control points can be selected at different positions. Compared with blind selection of shear wave velocity control points, the selection of the positions of the above-mentioned shear wave velocity control points can effectively constrain the shear wave model. Based on the shear wave velocity control points, the common center point domain data of the corresponding shot point is selected and secondary constrained inversion is performed to obtain a secondary inversion velocity model with high accuracy. The present application is applicable to shear wave seismic exploration projects. By picking the first arrival of shear waves on StRt component data and applying the secondary constrained inversion modeling method, it better adapts to the rapid changes in shear wave velocity, improves the accuracy of the shear wave model, eliminates the low-amplitude anomaly caused by the surface layer, and lays a solid model foundation for solving shear wave static correction and pre-stack depth migration imaging.

[0069] In one possible implementation, the initial SH wave component seismic data is acquired through a three-component geophone.

[0070] In one possible implementation, based on the initial inversion of the shear wave velocity model, a corresponding number of shear wave velocity control points are obtained at different positions, including:

[0071] Based on the initial inversion of the shear wave velocity model, the velocity interface fluctuations and velocity variation characteristics are obtained;

[0072] Based on the characteristics of velocity interface fluctuation and velocity variation, a corresponding number of shear wave velocity control points are obtained at different positions.

[0073] In one possible implementation, based on the velocity interface fluctuation and velocity variation characteristics, a corresponding number of shear wave velocity control points are obtained at different positions, including:

[0074] Based on the characteristics of velocity interface fluctuation and velocity change, the corresponding area is extracted at different positions;

[0075] Based on the area of the region, a corresponding number of shear wave velocity control points are obtained.

[0076] In one possible implementation, obtaining corresponding common center point domain data based on the shot point information includes:

[0077] Based on the shot point information, the corresponding near-arrangement first arrival data is obtained;

[0078] Get the common center point domain data corresponding to the nearly arranged first arrival data.

[0079] Figure 2 This is a flow chart of a shear wave vector inversion modeling method provided in an embodiment of the present application. Figure 2 , the method can be applied to a computer device, the method comprising:

[0080] 201. Rotate the initial SH wave component seismic data in the area in the direction of the shot-detection line.

[0081] In seismic data, P waves represent longitudinal waves, SH waves represent waves in which particle vibrations occur in a plane parallel to the wave propagation plane, and SV waves represent waves in which particle vibrations occur in a plane perpendicular to the wave propagation plane.

[0082] In one possible implementation, the initial SH wave component seismic data is acquired using a three-component geophone. In this step, because shear waves in seismic data exhibit directional propagation, the geophone receives maximum shear wave energy only when the shear wave's vibration direction is parallel to the geophone's horizontal component (X or Y) reception direction. To enhance the discernibility of the shear wave's first arrival, a rotational transformation can be performed on the geophone's received signal, converting the X and Y components of the SH wave signal to the Y component.

[0083] The principle and process of this rotation step are as follows:

[0084] In three-dimensional SH wave data acquisition, for seismic data received by a three-component geophone, the X component is parallel to the inline direction, the Y component is parallel to the crossline direction, and the Z component is vertically downward. The SH wave with maximum energy will be received only when the SH wave vibration direction is parallel to the receiving direction.

[0085] The three-dimensional seismic observation system arranges shot lines and receiver lines. The position of each shot point varies relative to the position of the receiver point. There is no way to arrange a single component of the receiver at a fixed angle to all shot points. In other words, with the existing three-component receiver, there is no way to make one component capture the full-energy SH wave. This is because the other component will inevitably capture a portion of the SH wave energy. Moreover, it is difficult to change the receiving direction of all the receivers at each excitation. Therefore, during data acquisition, the three-component receivers can only be arranged along a certain survey line direction. After the data is acquired, it is processed to eliminate the directional influence of the SH wave, that is, the rotation change.

[0086] Three-dimensional SH waves collect seismic signals tangentially along a radial circle centered on the earthquake source, projected onto components in the X and Y directions. The receiving angles for each shot and each receiver are different, requiring a certain degree of rotation. The resulting local coordinates change relative to the global coordinates (x, y, z). Given that the Z component contains no SH wave energy, it does not need to change during the rotation. Therefore, the two coordinate systems are converted using the following formula 1:

[0087]

[0088]

[0089] Among them, x'=x'(t), y'=y'(t), z'=z'(t), represent the rotated SH wave seismic time signal or time series; x=x(t), y=y(t), z=z(t), represent the original observed SH wave seismic time signal or time series; φ represents the rotation angle, that is, the angle between the line connecting the shot point and the receiver point and the survey line, which can be calculated from the position coordinates of the source and the receiver point. The rotation process can be seen in Figure 3 .

[0090] 202. Based on the rotationally transformed SH wave component seismic data, obtain the StRt component seismic data perpendicular to the shot-detection line.

[0091] In this step, for the rotationally transformed SH wave component seismic data, a component perpendicular to the shot-detection line, namely the StRt component seismic data, is extracted for use in subsequent steps.

[0092] 203. Based on the StRt component seismic data, obtain the first arrival data of the StRt component.

[0093] Among them, the first arrival data refers to the time when the seismic wave first arrives.

[0094] 204. Based on the first arrival data of the StRt component and the first constraint model, the first inversion shear wave velocity model is obtained through the first constraint inversion.

[0095] The initial constraint model can be generated using gradient or other constraint information. For a certain depth range from the surface to the underground (for example, 500-1000 meters below the surface), this range is gridded. Appropriate parameters can be selected based on actual needs. Each layer of grid is assigned an initial velocity value, such as 500 m / s for the first layer, 510 m / s for the second layer, and 520 m / s for the third layer. The entire range is filled according to a certain pattern. This pattern can vary in a gradient, ultimately forming velocity data that varies with depth from top (surface) to bottom (underground).

[0096] The first constraint model is a nonlinear model inversion technique that uses the first arrival wave information of seismic records to invert the velocity model of the near-surface structure. It reduces the complex surface geological model to microelement, assumes that the medium in the microelement is stable, and uses the network method for ray forward modeling. The travel time of the seismic wave is described as the line integral of the medium slowness function along the ray path, as shown in the following formula 2:

[0097]

[0098] Where: S(x,z) is the slowness function of the underground medium; dl is the differential of the ray path; T is the travel time of the seismic wave from the source point s to the receiving point r.

[0099] Tomographic inversion is a method of inverting the slowness function S(x,z) of the near-surface velocity field based on the known wave travel time matrix T.

[0100] The initial inversion shear wave velocity model can represent new velocity data that changes with depth from top (surface) to bottom (underground). Compared with the initial constraint model, the initial inversion shear wave velocity model uses the first arrival to participate in the calculation, which to a certain extent reflects the changes from the surface to the underground and will be more consistent with geological laws.

[0101] 205. Based on the initial inversion of the shear wave velocity model, the velocity interface fluctuation and velocity variation characteristics are obtained.

[0102] In this step, the fluctuation of velocity interface and the characteristics of velocity change are clarified, which can make clear the general change law of shear wave model, that is, to establish the approximate position where constraint control is required, especially the reasonable position control is required in complex areas. Figure 4 As shown in Figure 2, the initial inversion shows that there is a velocity anomaly area within the range L on the surface.

[0103] 206. Based on the characteristics of velocity interface fluctuation and velocity variation, the corresponding area is extracted at different locations.

[0104] 207. Based on the area of the region, obtain the corresponding number of shear wave velocity control points.

[0105] like Figure 4 As shown in the figure, according to the position of the abnormal area, speed control is required to establish the three positions ABC, so as to clarify the position for the secondary more precise constraint control.

[0106] 208. Based on the StRt component seismic data of the shear wave velocity control point, the shot point information of the corresponding position is extracted.

[0107] In this step, the shot points at the corresponding positions are obtained according to the positions of the shear wave velocity control points.

[0108] 209. Based on the shot point information, obtain the corresponding near-arrangement first arrival data.

[0109] Select the StRt component 3D seismic data of the shear wave velocity control point (such as A, B, C), and extract the nearly arranged first arrival data of the 3D StRt shot point at the position, such as Figure 5 shown.

[0110] Selecting the first arrival of a near arrangement includes: 3D observation is a shot excitation with many arrangement line receptions, such as Figure 6 As shown, the five-pointed star is the gun point, and you can choose the nearest arrangement on the right (or left) (such as Figure 6 Output the first arrival of this shot and the selected arrangement. The picking software can output the first arrival.

[0111] The above extraction refers to the selection of the arrangement of StRt component shot point data. For example, if the three-dimensional shot point arrangement is 24-line reception, the StRt component initial arrival data of 1 line on the left and right of the shot point can be selected ( Figure 6 ), since shear waves change more dramatically than longitudinal waves, the purpose is to improve the recognition of the inflection point of refraction stratification.

[0112] 210. Obtain the common center point domain data corresponding to the nearly arranged first arrival data.

[0113] Among them, the Common Middle Point (CMP) refers to extracting the tracks with common middle points from different shot sets to form a new set.

[0114] 211. Perform refraction layering interpretation in the common center point domain of the shot point to obtain the velocity and thickness of the shear wave velocity control point.

[0115] Refraction layering is based on the positional relationship between the shot point and the receiver point, and the refraction interpretation is performed to obtain the StRt component velocity thickness parameters, which form the StRt component shear wave velocity control points. The information of the shear wave velocity control points is used as a constraint for the StRt component primary to secondary inversion to establish a more refined shear wave surface model. Figure 7 The first arrival time-distance diagram shown is used for refraction interpretation. Refraction interpretation is the usual time-distance diagram interpretation. The horizontal axis of the time-distance diagram is the offset distance and the vertical axis is the first arrival time. The diagram shows three groups of first arrival information: CMP1 first arrival, CMP2 first arrival and CMP3 first arrival. In this way, the thickness and velocity information can be interpreted based on the refraction theory.

[0116] Based on the first arrival data of the StRt component, the primary constraint model, and the velocity and thickness of the shear wave velocity control points, a secondary inversion velocity model is obtained through secondary constraint inversion.

[0117] In this step, the information of the shear wave velocity control points is used as a constraint for the primary to secondary inversion of the StRt component to establish a more refined shear wave surface model.

[0118] The method provided in the embodiment of the present application obtains a primary inversion shear wave velocity model by adopting primary constrained inversion. According to the model, a corresponding number of shear wave velocity control points can be selected at different positions. Compared with blind selection of shear wave velocity control points, the selection of the positions of the above-mentioned shear wave velocity control points is conducive to improving the accuracy of the model and can reflect... Based on the shear wave velocity control points, the common center point domain data of the corresponding shot point is selected to perform secondary constrained inversion, thereby obtaining a secondary inversion velocity model with high accuracy. The present application is applicable to shear wave seismic exploration projects. By picking up the first arrival of shear waves on StRt component data and applying a secondary constrained inversion modeling device, it better adapts to the rapid changes in shear wave velocity, improves the accuracy of the shear wave model, eliminates the low-amplitude anomaly caused by the surface layer, and lays a solid model foundation for solving shear wave static correction and pre-stack depth migration imaging.

[0119] This solution solves the problem of low shear wave velocity model inversion accuracy caused by the blindness of shear wave micro-logging constraints on shear wave velocity control points in shear wave tomography inversion, and adapts to the drastic changes in shear wave surface velocity. This solution is simple and effective, does not rely on shear wave surface survey data, avoids costly shear wave micro-logging surveys, improves the accuracy of shear wave surface inversion modeling, and lays a foundation for solving low-amplitude anomalies caused by drastic changes in surface shear wave velocity and improving the quality of shear wave seismic profiles.

[0120] This application explores and applies nine-component P- and S-wave 3D seismic exploration in the Sanhu area of the Qaidam Basin. The Sanhu area is a natural gas-rich region of the Qaidam Basin, and joint exploration using SH and SV waves was employed. To improve the processing quality of SH and SV shear-wave seismic data, a quadratic constrained inversion modeling method was applied to better adapt to rapid changes in shear wave velocity, improve the accuracy of the shear wave model, eliminate low-amplitude anomalies caused by surface layers, and lay a solid model foundation for solving shear wave static corrections and prestack depth migration imaging.

[0121] All of the above optional technical solutions can be combined in any way to form optional embodiments of the present application, and will not be described in detail here.

[0122] Figure 8 This is a schematic diagram of the structure of a shear wave vector inversion modeling device provided in an embodiment of the present application. Figure 8 , the device comprises:

[0123] The rotation module 801 is used to rotate the initial SH wave component seismic data in the area toward the direction of the shot-detection line;

[0124] Component acquisition module 802, for acquiring StRt component seismic data perpendicular to the shot-detection line direction based on the SH wave component seismic data after rotation transformation;

[0125] A first arrival acquisition module 803 is used to acquire first arrival data of the StRt component based on the StRt component seismic data;

[0126] A primary inversion module 804 is configured to obtain a primary inverted shear wave velocity model through primary constrained inversion based on the first arrival data of the StRt component and the primary constrained model;

[0127] The control point acquisition module 805 is used to acquire a corresponding number of shear wave velocity control points at different positions based on the initial inversion shear wave velocity model;

[0128] The shot point acquisition module 806 is used to extract the shot point information of the corresponding position based on the StRt component seismic data of the shear wave velocity control point;

[0129] The constraint condition module 807 is used to obtain corresponding common center point domain data based on the shot point information;

[0130] The constraint condition module 807 is also used to perform refraction layer interpretation in the common center point domain of the shot point to obtain the velocity and thickness of the shear wave velocity control point;

[0131] The secondary inversion module 808 is used to obtain a secondary inversion velocity model through secondary constrained inversion based on the first arrival data of the StRt component, the primary constrained model, and the velocity and thickness of the shear wave velocity control point.

[0132] In one possible implementation, the initial SH wave component seismic data is acquired through a three-component geophone.

[0133] In one possible implementation, the control point acquisition module 805 is configured to:

[0134] Based on the initial inversion of the shear wave velocity model, the velocity interface fluctuations and velocity variation characteristics are obtained;

[0135] Based on the characteristics of velocity interface fluctuation and velocity variation, a corresponding number of shear wave velocity control points are obtained at different positions.

[0136] In one possible implementation, the control point acquisition module 805 is configured to:

[0137] Based on the characteristics of velocity interface fluctuation and velocity change, the corresponding area is extracted at different positions;

[0138] Based on the area of the region, a corresponding number of shear wave velocity control points are obtained.

[0139] In one possible implementation, the constraint condition module 807 is configured to:

[0140] Based on the shot point information, the corresponding near-arrangement first arrival data is obtained;

[0141] Get the common center point domain data corresponding to the nearly arranged first arrival data.

[0142] The device provided in the embodiment of the present application obtains a primary inversion shear wave velocity model by adopting primary constrained inversion. According to the model, a corresponding number of shear wave velocity control points can be selected at different positions. Compared with blind selection of shear wave velocity control points, the selection of the positions of the above-mentioned shear wave velocity control points can effectively constrain the shear wave model. Based on the shear wave velocity control points, the common center point domain data of the corresponding shot point is selected to perform secondary constrained inversion, thereby obtaining a secondary inversion velocity model with high accuracy. The present application is applicable to shear wave seismic exploration projects. By picking up the first arrival of shear waves on StRt component data and applying the secondary constrained inversion modeling device, it better adapts to the rapid changes in shear wave velocity, improves the accuracy of the shear wave model, eliminates the low-amplitude anomaly caused by the surface layer, and lays a solid model foundation for solving shear wave static correction and pre-stack depth migration imaging.

[0143] Figure 9 : This is a schematic diagram of the structure of a computer device provided in an embodiment of the present application. The computer device 900 may have relatively large differences due to different configurations or performances, and may include one or more processors (central processing units, CPU) 901 and one or more memories 902, wherein the memory 902 stores at least one program code, and the at least one program code is loaded and executed by the processor 901 to implement the shear wave vector inversion modeling method provided in each of the above method embodiments. Of course, the computer device may also have components such as a wired or wireless network interface, a keyboard, and an input and output interface for input and output. The computer device may also include other components for implementing device functions, which will not be described in detail here.

[0144] In some embodiments, the computer program involved in the embodiments of the present application can be deployed and executed on a computer device, or on multiple computer devices located at one location, or on multiple computer devices distributed at multiple locations and interconnected through a communication network. Multiple computer devices distributed at multiple locations and interconnected through a communication network can constitute a blockchain system.

[0145] In an exemplary embodiment, a computer-readable storage medium is also provided, such as a memory including program code. The program code can be executed by a processor in a computer device to perform the shear wave vector inversion modeling method in the above embodiment. For example, the computer-readable storage medium can be a read-only memory (ROM), a random access memory (RAM), a compact disc read-only memory (CD-ROM), a magnetic tape, a floppy disk, an optical data storage device, or the like.

[0146] Those skilled in the art will understand that all or part of the steps for implementing the above embodiments may be accomplished by hardware, or may be accomplished by instructing the relevant hardware through a program, and the above program may be stored in a computer-readable storage medium, which may be a read-only memory, a disk, or an optical disk, etc.

[0147] The above are only optional embodiments of the present application and are not intended to limit the present application. Any modifications, equivalent replacements, improvements, etc. made within the spirit and principles of the present application should be included in the scope of protection of the present application.

Claims

1. A shear wave vector inversion modeling method, characterized in that: The method comprises: The initial SH wave component seismic data in the area are rotated toward the direction of the shot-detection line; Based on the rotationally transformed SH wave component seismic data, the StRt component seismic data perpendicular to the shot-detection line is obtained; Based on the StRt component seismic data, the first arrival data of the StRt component is obtained; Based on the first arrival data of the StRt component and the first constraint model, the first inversion shear wave velocity model is obtained through the first constraint inversion; Based on the initial inversion shear wave velocity model, a corresponding number of shear wave velocity control points are obtained at different positions; Based on the StRt component seismic data of the shear wave velocity control point, the shot point information of the corresponding position is extracted; Based on the shot point information, corresponding common center point domain data is acquired; Perform refraction layer interpretation in the common center point domain of the shot point to obtain the velocity and thickness of the shear wave velocity control point; Based on the first arrival data of the StRt component, the primary constraint model, and the velocity and thickness of the shear wave velocity control points, the secondary inversion velocity model is obtained through secondary constraint inversion.

2. The method according to claim 1, characterized in that The initial SH wave component seismic data has three components.

3. The method according to claim 1, characterized in that The method of obtaining a corresponding number of shear wave velocity control points at different locations based on the initial inversion shear wave velocity model includes: Based on the initial inversion of the shear wave velocity model, the velocity interface fluctuations and velocity variation characteristics are obtained; Based on the characteristics of velocity interface fluctuation and velocity variation, a corresponding number of shear wave velocity control points are obtained at different positions.

4. The method according to claim 3, characterized in that The method of obtaining a corresponding number of shear wave velocity control points at different positions based on the velocity interface fluctuation and velocity variation characteristics includes: Based on the characteristics of velocity interface fluctuation and velocity change, the corresponding area is extracted at different positions; Based on the area of the region, a corresponding number of shear wave velocity control points are obtained.

5. The method according to claim 1, wherein The acquiring corresponding common center point domain data based on the shot point information includes: Based on the shot point information, obtaining corresponding near-arrangement first arrival data; Get the common center point domain data corresponding to the nearly arranged first arrival data.

6. A shear wave vector inversion modeling device, characterized in that: The device comprises: The rotation module is used to rotate the initial SH wave component seismic data in the area toward the direction of the shot-detection line; A component acquisition module is used to acquire StRt component seismic data perpendicular to the shot-detection line based on the SH wave component seismic data after rotation transformation; A first arrival acquisition module is used to acquire first arrival data of the StRt component based on the StRt component seismic data; The primary inversion module is used to obtain the primary inversion shear wave velocity model based on the first arrival data of the StRt component and the primary constraint model through the primary constraint inversion; A control point acquisition module is used to obtain a corresponding number of shear wave velocity control points at different positions based on the initial inversion shear wave velocity model; The shot point acquisition module is used to extract the shot point information of the corresponding position based on the StRt component seismic data of the shear wave velocity control point; A constraint condition module, configured to obtain corresponding common center point domain data based on the shot point information; The constraint condition module is also used to perform refraction layer interpretation in the common center point domain of the shot point to obtain the velocity and thickness of the shear wave velocity control point; The secondary inversion module is used to obtain the secondary inversion velocity model through secondary constrained inversion based on the first arrival data of the StRt component, the primary constraint model, and the velocity and thickness of the shear wave velocity control point.

7. The device according to claim 6, characterized in that The initial SH wave component seismic data has three components.

8. The device according to claim 6, characterized in that The control point acquisition module is used to: Based on the initial inversion of the shear wave velocity model, the velocity interface fluctuations and velocity variation characteristics are obtained; Based on the characteristics of velocity interface fluctuation and velocity variation, a corresponding number of shear wave velocity control points are obtained at different positions.

9. The device according to claim 8, characterized in that The control point acquisition module is used to: Based on the characteristics of velocity interface fluctuation and velocity change, the corresponding area is extracted at different positions; Based on the area of the region, a corresponding number of shear wave velocity control points are obtained.

10. The device according to claim 6, characterized in that The constraint condition module is used to: Based on the shot point information, obtaining corresponding near-arrangement first arrival data; Get the common center point domain data corresponding to the nearly arranged first arrival data.

Images