CN115774288A - Shear wave vector inversion modeling method and device - Google Patents
Shear wave vector inversion modeling method and device Download PDFInfo
- Publication number
- CN115774288A CN115774288A CN202111052075.0A CN202111052075A CN115774288A CN 115774288 A CN115774288 A CN 115774288A CN 202111052075 A CN202111052075 A CN 202111052075A CN 115774288 A CN115774288 A CN 115774288A
- Authority
- CN
- China
- Prior art keywords
- inversion
- model
- acquiring
- shear wave
- component
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Images
Abstract
The application provides a shear wave vector inversion modeling method and device, and belongs to the technical field of shear wave static correction. According to the technical scheme, a primary inversion transverse wave velocity model is obtained by adopting primary constraint inversion, transverse wave velocity control points can be selected at different positions according to the model, so that effective constraint can be carried out on the transverse wave model, corresponding common center point domain data is obtained based on the transverse wave velocity control points, secondary constraint inversion is carried out, and a secondary inversion velocity model is obtained and has high precision. The method is suitable for the transverse wave seismic exploration project, transverse wave first arrivals are picked up on StRt component data, and a quadratic constraint inversion modeling method is applied, so that the method is better suitable for the rapid change of the transverse wave speed, the transverse wave model precision is improved, the low-amplitude abnormal phenomenon caused by the surface layer is eliminated, and a solid model foundation is laid for the solution of transverse wave static correction and prestack depth migration imaging.
Description
Technical Field
The application relates to the technical field of shear wave static correction, in particular to a shear wave vector inversion modeling method and device.
Background
In the natural gas exploitation technology, seismic data of a stratum need to be processed to obtain a gas-containing zone structure in the stratum, the seismic data comprise longitudinal wave seismic data and transverse wave seismic data, and due to the fact that the stratum contains gas, imaging precision of the longitudinal wave seismic data is low, surface layer transverse wave structural change in the transverse wave seismic data can cause abnormity, the abnormity and a low-amplitude structure are mixed and are difficult to distinguish, and therefore research on transverse wave vector inversion modeling is conducted.
The currently commonly implemented shear wave vector inversion modeling method adopts three seismic sources of P wave, SH wave and SV wave for excitation, receives nine-component three-dimensional vector seismic data through a three-component detector, thereby performing shear wave vector inversion modeling, and is influenced by a three-dimensional observation position and the propagation direction of shear waves, and the energy of the shear waves received by the three-component detector is constantly changed, so that the shear waves are difficult to pick up at the first arrival; meanwhile, the surface shear wave velocity change is more severe than the longitudinal wave, the constraint information is difficult to obtain accurately, the shear wave change characteristic is difficult to control, and the first arrival chromatography inversion precision is influenced.
Disclosure of Invention
The embodiment of the application provides a shear wave vector inversion modeling method and device, which can improve the precision of a shear wave model, eliminate low-amplitude abnormal phenomena caused by a surface layer and lay a solid model foundation for solving shear wave static correction and prestack depth migration imaging. The technical scheme is as follows:
in one aspect, a shear wave vector inversion modeling method is provided, and the method includes:
carrying out rotation transformation on the initial SH wave component seismic data in the area towards the direction of a shot-detection connecting line;
acquiring StRt component seismic data vertical to the shot-inspection connecting line direction based on the SH wave component seismic data after the rotation transformation;
acquiring first arrival data of the StRt component based on the StRt component seismic data;
based on the first arrival data of the StRt component and the primary constraint model, obtaining a primary inversion shear wave velocity model through primary constraint inversion;
acquiring a corresponding number of shear wave velocity control points at different positions based on the primary inversion shear wave velocity model;
extracting shot point information of a corresponding position based on StRt component seismic data of the transverse wave velocity control point;
acquiring corresponding common center point domain data based on the shot point information;
performing refraction and layered interpretation in a common center point region of the shot point to obtain the speed and the thickness of a transverse wave speed control point;
and obtaining a quadratic inversion velocity model through quadratic constraint inversion based on the first arrival data of the StRt component, the primary constraint model and the velocity and the thickness of the shear wave velocity control point.
In one possible implementation, the initial SH wave component seismic data has three components.
In one possible implementation manner, the obtaining, at different positions, a corresponding number of shear wave velocity control points based on the primary inversion shear wave velocity model includes:
acquiring the fluctuation and speed change characteristics of a speed interface based on a primary inversion shear wave speed model;
and acquiring corresponding quantity of transverse wave speed control points at different positions based on the characteristics of fluctuation and speed change of the speed interface.
In one possible implementation manner, the acquiring, at different positions, a corresponding number of shear wave velocity control points based on the characteristics of the velocity interface fluctuation and the velocity variation includes:
extracting regions with corresponding areas at different positions based on the characteristics of fluctuation and speed change of a speed interface;
and acquiring a corresponding number of shear wave speed control points based on the area of the region.
In one possible implementation, the obtaining, based on the shot point information, corresponding common-center-point domain data includes:
acquiring corresponding near-ranking first-arrival data based on the shot point information;
and acquiring the co-center point domain data corresponding to the near-arrangement first-motion data.
In one aspect, a shear wave vector inversion modeling apparatus is provided, the apparatus including:
the rotation module is used for performing rotation transformation on the initial SH wave component seismic data in the region in the direction of the shot-detection connecting line;
the component acquisition module is used for acquiring StRt component seismic data which are vertical to the shot-detection connecting line direction based on the SH wave component seismic data after the rotation transformation;
the first arrival acquisition module is used for acquiring first arrival data of the StRt component based on the StRt component seismic data;
the primary inversion module is used for obtaining a primary inversion shear wave velocity model through primary constraint inversion based on the primary arrival data of the StRt component and a primary constraint model;
the control point acquisition module is used for acquiring the corresponding number of shear wave speed control points at different positions based on the primary inversion shear wave speed model;
the shot point acquisition module is used for extracting shot point information of a corresponding position based on StRt component seismic data of the transverse wave velocity control point;
the constraint condition module is used for acquiring corresponding common center point domain data based on the shot point information;
the constraint condition module is also used for carrying out refraction layering explanation on the common center point region of the shot point to obtain the speed and the thickness of the transverse wave speed control point;
and the quadratic inversion module is used for obtaining a quadratic inversion velocity model through quadratic constraint inversion based on the first arrival data of the StRt component, the primary constraint model and the velocity and the thickness of the shear wave velocity control point.
In one possible implementation, the initial SH wave component seismic data has three components.
In one possible implementation manner, the control point obtaining module is configured to:
acquiring the fluctuation and speed change characteristics of a speed interface based on a primary inversion shear wave speed model;
and acquiring corresponding quantity of transverse wave speed control points at different positions based on the fluctuation and speed change characteristics of the speed interface.
In one possible implementation manner, the control point obtaining module is configured to:
extracting regions with corresponding areas at different positions based on the characteristics of fluctuation and speed change of a speed interface;
and acquiring a corresponding number of shear wave speed control points based on the area of the region.
In one possible implementation, the constraint module is configured to:
acquiring corresponding near-ranking first-arrival data based on the shot point information;
and acquiring the co-center point domain data corresponding to the near-arrangement first-motion data.
The technical scheme that this application embodiment provided, through adopting first constraint inversion, obtain first inversion shear wave velocity model, according to this model, can select the shear wave velocity control point of corresponding quantity in the position of difference, compare in blind selection shear wave velocity control point, the selection of above-mentioned shear wave velocity control point position can carry out effectual constraint to the shear wave model, based on this shear wave velocity control point, select the common center point domain data that corresponds the shot point, carry out secondary constraint inversion, thereby obtain secondary inversion velocity model, this model has higher precision. The method is suitable for the transverse wave seismic exploration project, transverse wave first arrivals are picked up on StRt component data, and a quadratic constraint inversion modeling method is applied, so that the method is better suitable for the rapid change of the transverse wave speed, the transverse wave model precision is improved, the low-amplitude abnormal phenomenon caused by the surface layer is eliminated, and a solid model foundation is laid for the solution of transverse wave static correction and prestack depth migration imaging.
Drawings
In order to more clearly illustrate the technical solutions in the embodiments of the present application, the drawings required to be used in the description of the embodiments are briefly introduced below, and it is obvious that the drawings in the description below are only some embodiments of the present application, and it is obvious for those skilled in the art to obtain other drawings without creative efforts.
Fig. 1 is a flowchart of a shear wave vector inversion modeling method according to an embodiment of the present application;
fig. 2 is a flowchart of a shear wave vector inversion modeling method according to an embodiment of the present application;
FIG. 3 is a schematic diagram of a rotation principle provided by an embodiment of the present application;
FIG. 4 is a schematic diagram illustrating analysis of primary inversion results according to an embodiment of the present disclosure;
FIG. 5 is a shot and near arrangement diagram provided by an embodiment of the present application;
FIG. 6 is a close-up schematic view provided by an embodiment of the present application;
FIG. 7 is a first-arrival time-distance diagram provided by an embodiment of the present application;
fig. 8 is a schematic structural diagram of a shear wave vector inversion modeling apparatus provided in an embodiment of the present application;
fig. 9 is a schematic structural diagram of a computer device according to an embodiment of the present application.
Detailed Description
To make the objects, technical solutions and advantages of the present application more clear, embodiments of the present application will be described in further detail below with reference to the accompanying drawings.
In this application, unless expressly stated or limited otherwise, the terms "mounted," "connected," "secured," and the like are to be construed broadly and encompass, for example, being fixedly connected, releasably connected, or integrally connected; can be mechanically or electrically connected; they may be connected directly or indirectly through intervening media, or they may be interconnected between two elements. The specific meaning of the above terms in the present application can be understood by those of ordinary skill in the art as appropriate.
Fig. 1 is a flowchart of a shear wave vector inversion modeling method provided in an embodiment of the present application, please refer to fig. 1, which can be applied to a computer device, and the method includes:
101. and carrying out rotation transformation on the initial SH wave component seismic data in the region in the direction of the shot-survey connecting line.
102. And acquiring StRt component seismic data vertical to the shot-survey connecting line direction based on the SH wave component seismic data after the rotation transformation.
103. Based on the StRt component seismic data, first arrival data of the StRt component is acquired.
104. And obtaining a primary inversion shear wave velocity model through primary constraint inversion based on the primary motion data of the StRt component and the primary constraint model.
105. And acquiring corresponding quantity of shear wave velocity control points at different positions based on the primary inversion shear wave velocity model.
106. And extracting shot point information of a corresponding position based on the StRt component seismic data of the transverse wave velocity control point.
107. And acquiring corresponding common center point domain data based on the shot point information.
108. And performing refraction layered interpretation in the common center point region of the shot point to obtain the speed and the thickness of the shear wave speed control point.
109. And obtaining a quadratic inversion velocity model through quadratic constraint inversion based on the first arrival data of the StRt component, the primary constraint model and the velocity and the thickness of the shear wave velocity control point.
According to the method provided by the embodiment of the application, a primary inversion transverse wave velocity model is obtained by adopting primary constraint inversion, corresponding quantity of transverse wave velocity control points can be selected at different positions according to the model, compared with blind selection of transverse wave velocity control points, the selection of the positions of the transverse wave velocity control points can effectively constrain the transverse wave model, common center point domain data corresponding to shot points are selected based on the transverse wave velocity control points, secondary constraint inversion is carried out, and therefore a secondary inversion velocity model is obtained and has high precision. The method is suitable for the transverse wave seismic exploration project, transverse wave first arrivals are picked up on StRt component data, and a quadratic constraint inversion modeling method is applied, so that the method is better suitable for the rapid change of the transverse wave speed, the transverse wave model precision is improved, the low-amplitude abnormal phenomenon caused by the surface layer is eliminated, and a solid model foundation is laid for the solution of transverse wave static correction and prestack depth migration imaging.
In one possible implementation, the initial SH wave component seismic data is acquired by a three-component detector.
In one possible implementation manner, the obtaining, at different positions, a corresponding number of shear wave velocity control points based on the primary inversion shear wave velocity model includes:
acquiring the fluctuation and speed change characteristics of a speed interface based on a primary inversion shear wave speed model;
and acquiring corresponding quantity of transverse wave speed control points at different positions based on the characteristics of fluctuation and speed change of the speed interface.
In one possible implementation manner, the acquiring, at different positions, a corresponding number of shear wave velocity control points based on the characteristics of the velocity interface fluctuation and the velocity variation includes:
extracting regions with corresponding areas at different positions based on the characteristics of fluctuation and speed change of a speed interface;
and acquiring a corresponding number of shear wave speed control points based on the area of the region.
In one possible implementation, the obtaining, based on the shot point information, corresponding common-center-point domain data includes:
acquiring corresponding near-ranking first-arrival data based on the shot point information;
and acquiring the co-center point domain data corresponding to the near-arrangement first-motion data.
Fig. 2 is a flowchart of a shear wave vector inversion modeling method provided in an embodiment of the present application, please refer to fig. 2, which can be applied to a computer device, and the method includes:
201. and carrying out rotation transformation on the initial SH wave component seismic data in the region towards the shot-detection connecting line direction.
In seismic data, a P wave represents a longitudinal wave, an SH wave represents a wave in which particle vibration occurs in a plane parallel to the propagation plane of the wave, and an SV wave represents a wave in which particle vibration occurs in a plane perpendicular to the propagation plane of the wave.
In one possible implementation, the initial SH wave component seismic data is acquired by a three-component detector. In this step, since the transverse waves in the seismic data have directivity in propagation, the maximum transverse wave energy is received by the detector only when the vibration direction of the transverse waves is parallel to the reception direction of the horizontal component (X or Y) of the detector. In order to increase the first arrival discrimination of the transverse wave, the signal received by the detector may be subjected to a rotation transformation to convert the SH wave signals of the X and Y components into the Y component.
The principle and process of this rotation step is as follows:
in the data acquisition of three-dimensional SH waves, for seismic data received by a three-component detector, an X component is parallel to an inline direction, a Y component is parallel to a crossline direction, a Z component is vertically downward, and the SH waves with maximum energy can be received only when the SH wave vibration direction is parallel to the receiving direction.
The three-dimensional seismic observation system is provided with a shot line and a wave detection line, the position of each shot point is changed relative to the position of the wave detection point, and there is no way to arrange a certain component of a wave detector to form a fixed angle with all the shot points, that is, the existing three-component wave detector has no way to enable one component to obtain the full-energy SH wave, because the other component must obtain a part of SH wave energy, and the receiving directions of all the wave detectors are difficult to change during each excitation, therefore, the three-component wave detector can only be arranged according to a certain measuring line direction during the acquisition work, certain processing is carried out after the data are obtained, and the directional influence of the SH wave, namely the rotation change, is eliminated.
The three-dimensional SH wave is acquired by projecting seismic signals along the tangential direction on a radiation circle with a seismic source as the center onto components in the X direction and the Y direction, the receiving angle of each geophone point of each geophone is different, rotation of a certain angle is required, and the position of a local coordinate obtained after rotation is changed relative to a global coordinate (X, Y, z). Considering that the Z component has no SH wave energy, the Z component need not change during rotation. Thus, the two coordinate systems are converted by the following equation 1:
wherein x '= x' (t), y '= y' (t), z '= z' (t), which represents the rotated SH-wave seismic time signal or time series; x = x (t), y = y (t), z = z (t), representing the original observed SH wave seismic time signal or time series; phi represents a rotation angle, namely an included angle between a shot point-demodulator probe connecting line and a survey line, and can be calculated by the position coordinates of the seismic source and the position coordinates of the demodulator probe, and the rotation process can be seen in fig. 3.
202. And acquiring StRt component seismic data vertical to the shot-detection connecting line direction based on the SH wave component seismic data after the rotation transformation.
In this step, for the SH wave component seismic data after rotation transformation, a component perpendicular to the shot-to-survey line direction, that is, stRt component seismic data, is extracted for the subsequent step.
203. Based on the StRt component seismic data, first arrival data of the StRt component is acquired.
The first arrival data refers to the time of the first arrival of seismic waves.
204. And obtaining a primary inversion shear wave velocity model through primary constraint inversion based on the primary motion data of the StRt component and the primary constraint model.
The primary constraint model may be generated by gradient or other constraint information, and for a certain depth range from the earth surface to the underground (for example, 500-1000m below the earth surface), the range is gridded, and suitable parameters may be selected according to actual needs, each layer of grid is assigned with an initial velocity value, for example, the first layer is 500m/s, the second layer is 510m/s, the third layer is 520m/s, and the whole range is filled according to a certain rule. This law may be gradient, resulting in a velocity data over depth from top (surface) to bottom (subsurface).
The primary constraint model is a nonlinear model inversion technology, velocity model inversion is carried out on a near-surface structure by using the information of the primary waves of seismic records, a complex surface geological model is subjected to infinitesimal modeling, a network method is used for carrying out ray forward modeling on the assumption that a medium in the infinitesimal is stable and unchangeable, and the travel time of seismic waves is described as the line integral of a medium slowness function along a ray path, and the following formula 2 is adopted:
in the formula: s (x, z) is a slowness function of the subsurface 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.
Tomographic inversion is a method for inverting the near-surface velocity field slowness function S (x, z) according to the travel time matrix T of the known wave.
The primary inversion shear wave velocity model can represent new velocity data changing with depth from top (earth surface) to bottom (underground), and compared with a primary constraint model, the primary inversion shear wave velocity model utilizes the first arrival to participate in calculation, reflects the change from the earth surface to the underground to a certain extent, and can better accord with geological rules.
205. And acquiring the characteristics of speed interface fluctuation and speed change based on the primary inversion shear wave speed model.
In the step, the fluctuation and speed change characteristics of the speed interface are clarified, and the general change rule of the transverse wave model can be clarified, namely the general position needing to be restrained and controlled is determined, and particularly the control of the reasonable position in a complex area is needed. As shown in fig. 4, a velocity abnormality region exists in the range of the table L by the first inversion.
206. And extracting regions with corresponding areas at different positions based on the fluctuation and speed change characteristics of the speed interface.
207. And acquiring a corresponding number of shear wave speed control points based on the area of the region.
As shown in fig. 4, establishing three ABC positions requires velocity control based on the abnormal area position, and the positions are cleared for secondary finer constraint control.
208. And extracting shot point information of a corresponding position based on the StRt component seismic data of the transverse wave velocity control point.
In this step, the shot point at the corresponding position is obtained according to the position of the shear wave velocity control point.
209. And acquiring corresponding near-ranking first-arrival data based on the shot point information.
The three-dimensional seismic data of the StRt components at the shear wave velocity control point (for example, A, B, C) are selected, and near alignment first arrival data around the three-dimensional StRt shot at the position is extracted, as shown in fig. 5.
Selecting a near alignment first arrival comprises: the three-dimensional observation is that a shot firing has many lines of alignment to receive, as shown in fig. 6, and the five-pointed star is the shot point, and the right (or left) nearest line of alignment (as shown by the dotted line in fig. 6), or left and right, can be selected. The cannon and the selected arrangement are output corresponding to first arrivals, and the picking software can output the first arrivals.
The above extraction is to select the arrangement of the StRt component shot data, and if the three-dimensional shot arrangement is 24-line reception, the StRt component first arrival data of 1 line of each of the left and right sides of the shot (fig. 6) can be selected, and the transverse wave changes more sharply than the longitudinal wave, so as to improve the refraction stratification inflection point discrimination.
210. And acquiring the co-center point domain data corresponding to the near-arrangement first-motion data.
Wherein, common Middle Point (CMP) means that the channels having Common Middle points in different shot sets are extracted to form a new set.
211. And performing refraction layered interpretation in the common center point region of the shot point to obtain the speed and the thickness of the shear wave speed control point.
The refraction layering is to perform refraction explanation according to the position relation of a shot point and a receiving point to obtain a StRt component speed thickness parameter and form a StRt component transverse wave speed control point. And (3) using the information of the transverse wave speed control point as a constraint condition for StRt component first-arrival quadratic inversion to establish a more refined transverse wave surface model. Referring to the first arrival time interval diagram shown in fig. 7, refraction explanation is a normal time interval diagram explanation, the horizontal axis of the time interval diagram is offset, the vertical axis is first arrival time, and three sets of first arrival information, namely CMP1 first arrival, CMP2 first arrival and CMP3 first arrival, are shown in the diagram, so that thickness and speed information can be interpreted according to refraction theory.
212. And obtaining a quadratic inversion velocity model through quadratic constraint inversion based on the first arrival data of the StRt component, the primary constraint model and the velocity and the thickness of the shear wave velocity control point.
The step uses the information of the transverse wave speed control point as a constraint condition for StRt component first-arrival quadratic inversion, and establishes a more refined transverse wave surface model.
According to the method provided by the embodiment of the application, a primary inversion transverse wave velocity model is obtained by adopting primary constraint inversion, corresponding quantity of transverse wave velocity control points can be selected at different positions according to the model, compared with blind selection of the transverse wave velocity control points, the selection of the position of the transverse wave velocity control point is beneficial to improving the precision of the model and can reflect, and secondary constraint inversion is carried out by selecting common-center-point domain data of corresponding shot points based on the transverse wave velocity control points, so that a secondary inversion velocity model is obtained, and the model has higher precision. The method is suitable for the transverse wave seismic exploration project, transverse wave first arrivals are picked up on StRt component data, and the application of the secondary constraint inversion modeling device better adapts to the rapid change of the transverse wave speed, improves the transverse wave model precision, eliminates the low-amplitude abnormal phenomenon caused by the surface layer, and lays a solid model foundation for the solution of transverse wave static correction and prestack depth migration imaging.
The scheme solves the problem of low inversion accuracy of the transverse wave velocity model caused by the blindness of the transverse wave micro-logging constraint transverse wave velocity control point in transverse wave chromatography inversion, adapts to the severe change of the transverse wave surface velocity, is simple and effective, does not depend on transverse wave surface survey data, avoids the transverse wave micro-logging survey with high cost, improves the transverse wave surface inversion modeling accuracy, and lays a foundation for solving the low-amplitude abnormal phenomenon caused by the severe change of the surface transverse wave velocity and improving the transverse wave seismic profile quality.
The method is practically explored and applied to three-dimensional seismic exploration of nine components of longitudinal and transverse waves in three lakes of the Chadamu basin. The three-lake region is a natural gas enrichment region of a firewood basin and is explored by adopting SH and SV waves in a combined mode. In order to improve the processing quality of SH and SV transverse wave seismic data, the method better adapts to the rapid change of transverse wave speed through the application of a quadratic constraint inversion modeling method, improves the precision of a transverse wave model, eliminates the low-amplitude abnormal phenomenon caused by a surface layer, and lays a solid model foundation for the solution of transverse wave static correction and prestack depth migration imaging.
All the above optional technical solutions may be combined arbitrarily to form optional embodiments of the present application, and are not described in detail herein.
Fig. 8 is a schematic structural diagram of a shear wave vector inversion modeling apparatus provided in an embodiment of the present application, please refer to fig. 8, the apparatus includes:
a rotation module 801, configured to perform rotation transformation on the initial SH wave component seismic data in the region in the shot-survey connecting line direction;
a component obtaining module 802, configured to obtain, based on the SH wave component seismic data after the rotation transformation, stRt component seismic data perpendicular to the shot-survey line direction;
a first arrival acquisition module 803, configured to acquire first arrival data of the StRt component based on the StRt component seismic data;
the primary inversion module 804 is used for obtaining a primary inversion shear wave velocity model through primary constraint inversion based on the primary arrival data of the StRt component and the primary constraint model;
a control point obtaining module 805, configured to obtain, at different positions, a corresponding number of shear wave velocity control points based on the primary inversion shear wave velocity model;
a shot point obtaining module 806, configured to extract shot point information at a corresponding position based on the StRt component seismic data of the shear wave velocity control point;
a constraint condition module 807 for acquiring corresponding common center point domain data based on the shot point information;
the constraint condition module 807 is further configured to perform refraction and layering interpretation in the common-center-point region of the shot point to obtain the speed and thickness of the transverse wave speed control point;
and the quadratic inversion module 808 is configured to obtain a quadratic inversion velocity model through quadratic constraint 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.
In one possible implementation, the initial SH wave component seismic data is acquired by a three-component detector.
In one possible implementation, the control point obtaining module 805 is configured to:
acquiring the fluctuation and speed change characteristics of a speed interface based on a primary inversion shear wave speed model;
and acquiring corresponding quantity of transverse wave speed control points at different positions based on the characteristics of fluctuation and speed change of the speed interface.
In one possible implementation, the control point obtaining module 805 is configured to:
extracting regions with corresponding areas at different positions based on the characteristics of fluctuation and speed change of a speed interface;
and acquiring a corresponding number of shear wave speed control points based on the area of the region.
In one possible implementation, the constraint module 807 is configured to:
acquiring corresponding near-ranking first-arrival data based on the shot point information;
and acquiring the co-center point domain data corresponding to the near-arrangement first-motion data.
The device that this application embodiment provided, through adopting first constraint inversion, obtain first inversion shear wave velocity model, according to this model, can select the shear wave velocity control point of corresponding quantity in the position of difference, compare in blind selection shear wave velocity control point, the selection of above-mentioned shear wave velocity control point position can carry out effectual constraint to the shear wave model, based on this shear wave velocity control point, select the common midpoint territory data that corresponds the shot point, carry out secondary constraint inversion, thereby obtain secondary inversion velocity model, this model has higher precision. The method is suitable for the transverse wave seismic exploration project, transverse wave first arrivals are picked up on StRt component data, and the application of the secondary constraint inversion modeling device better adapts to the rapid change of the transverse wave speed, improves the transverse wave model precision, eliminates the low-amplitude abnormal phenomenon caused by the surface layer, and lays a solid model foundation for the solution of transverse wave static correction and prestack depth migration imaging.
Fig. 9 is a schematic structural diagram of a computer device provided in an embodiment of the present application, where the computer device 900 may generate relatively large differences due to different configurations or performances, and may include one or more processors (CPUs) 901 and one or more memories 902, where the memory 902 stores at least one program code, and the at least one program code is loaded and executed by the processors 901 to implement the shear wave vector inversion modeling method provided in the foregoing method embodiments. Certainly, the computer device may further have components such as a wired or wireless network interface, a keyboard, and an input/output interface, so as to perform input and output, and the computer device may further include other components for implementing the functions of the device, which is not described herein again.
In some embodiments, the computer program according to the embodiments of the present application may be deployed to be executed on one computer device or on multiple computer devices located at one site, or may be executed on multiple computer devices distributed at multiple sites and interconnected by a communication network, and the multiple computer devices distributed at the multiple sites and interconnected by the communication network may constitute a block chain system.
In an exemplary embodiment, a computer-readable storage medium, such as a memory, including program code, which is executable by a processor in a computer device to perform the shear wave vector inversion modeling method in the above embodiments is also provided. For example, the computer-readable storage medium may 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, and the like.
It will be understood by those skilled in the art that all or part of the steps for implementing the above embodiments may be implemented by hardware, or may be implemented by a program instructing relevant hardware, and the program may be stored in a computer-readable storage medium, and the storage medium may be a read-only memory, a magnetic disk or an optical disk.
The above description is only exemplary of the present application and should not be taken as limiting, as any modification, equivalent replacement, or improvement made within the spirit and principle of the present application should be included in the protection scope of the present application.
Claims (10)
1. A shear wave vector inversion modeling method, the method comprising:
carrying out rotation transformation on the initial SH wave component seismic data in the area towards the direction of a shot-detection connecting line;
acquiring StRt component seismic data vertical to the shot-inspection connecting line direction based on the SH wave component seismic data after the rotation transformation;
acquiring first-arrival data of the StRt component based on the seismic data of the StRt component;
based on the first arrival data of the StRt component and the primary constraint model, obtaining a primary inversion shear wave velocity model through primary constraint inversion;
acquiring a corresponding number of shear wave velocity control points at different positions based on the primary inversion shear wave velocity model;
extracting shot point information of a corresponding position based on StRt component seismic data of the transverse wave velocity control point;
acquiring corresponding common center point domain data based on the shot point information;
performing refraction and layered interpretation in a common center point region of the shot point to obtain the speed and the thickness of a transverse wave speed control point;
and obtaining a quadratic inversion velocity model through quadratic constraint inversion based on the first arrival data of the StRt component, the primary constraint model and the velocity and the thickness of the shear wave velocity control point.
2. The method of claim 1 wherein the initial SH wave component seismic data has three components.
3. The method of claim 1, wherein obtaining a corresponding number of shear wave velocity control points at different locations based on the primary inversion shear wave velocity model comprises:
acquiring the fluctuation and speed change characteristics of a speed interface based on a primary inversion shear wave speed model;
and acquiring corresponding quantity of transverse wave speed control points at different positions based on the characteristics of fluctuation and speed change of the speed interface.
4. The method of claim 3, wherein the obtaining a corresponding number of shear wave velocity control points at different locations based on the velocity interface undulation and velocity variation characteristics comprises:
extracting regions with corresponding areas at different positions based on the characteristics of fluctuation and speed change of a speed interface;
and acquiring a corresponding number of shear wave speed control points based on the area of the region.
5. The method of claim 1, wherein said obtaining corresponding isocenter domain data based on the shot point information comprises:
acquiring corresponding near-ranking first-arrival data based on the shot point information;
and acquiring common center point domain data corresponding to the near-arrangement first-arrival data.
6. A shear wave vector inversion modeling apparatus, the apparatus comprising:
the rotation module is used for performing rotation transformation on the initial SH wave component seismic data in the region in the direction of the shot-detection connecting line;
the component acquisition module is used for acquiring StRt component seismic data which are vertical to the shot-detection connecting line direction based on the SH wave component seismic data after the rotation transformation;
the first arrival acquisition module is used for acquiring first arrival data of the StRt component based on the StRt component seismic data;
the primary inversion module is used for obtaining a primary inversion shear wave velocity model through primary constraint inversion based on the primary arrival data of the StRt component and a primary constraint model;
the control point acquisition module is used for acquiring the corresponding number of shear wave speed control points at different positions based on the primary inversion shear wave speed model;
the shot point acquisition module is used for extracting shot point information of a corresponding position based on StRt component seismic data of the transverse wave velocity control point;
the constraint condition module is used for acquiring corresponding common center point domain data based on the shot point information;
the constraint condition module is also used for carrying out refraction layering explanation on the common center point region of the shot point to obtain the speed and the thickness of the transverse wave speed control point;
and the quadratic inversion module is used for obtaining a quadratic inversion velocity model through quadratic constraint inversion based on the first arrival data of the StRt component, the primary constraint model and the velocity and the thickness of the shear wave velocity control point.
7. The apparatus of claim 6 wherein the initial SH wave component seismic data has three components.
8. The apparatus of claim 6, wherein the control point obtaining module is configured to:
acquiring the fluctuation and speed change characteristics of a speed interface based on a primary inversion shear wave speed model;
and acquiring corresponding quantity of transverse wave speed control points at different positions based on the characteristics of fluctuation and speed change of the speed interface.
9. The apparatus of claim 8, wherein the control point obtaining module is configured to:
extracting regions with corresponding areas at different positions based on the fluctuation and speed change characteristics of a speed interface;
and acquiring a corresponding number of shear wave speed control points based on the area of the region.
10. The apparatus of claim 6, wherein the constraint module is configured to:
acquiring corresponding near-ranking first-arrival data based on the shot point information;
and acquiring common center point domain data corresponding to the near-arrangement first-arrival data.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111052075.0A CN115774288A (en) | 2021-09-08 | 2021-09-08 | Shear wave vector inversion modeling method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111052075.0A CN115774288A (en) | 2021-09-08 | 2021-09-08 | Shear wave vector inversion modeling method and device |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN115774288A true CN115774288A (en) | 2023-03-10 |
Family
ID=85388216
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111052075.0A Pending CN115774288A (en) | 2021-09-08 | 2021-09-08 | Shear wave vector inversion modeling method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115774288A (en) |
-
2021
- 2021-09-08 CN CN202111052075.0A patent/CN115774288A/en active Pending
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107526101B (en) | A kind of acquisition and processing method obtaining earthquake reflected wave | |
| RU2319982C2 (en) | Definition of characteristics of inclined cracks with the aid of exchanged reflected waves in combination with the process of consistent exclusion of influence of the superposed layers | |
| US8908471B2 (en) | Method for building velocity models for imaging in multi-azimuth marine seismic surveys | |
| AU2014227491B2 (en) | Construction and application of angle gathers from three-dimensional imaging of multiples wavefields | |
| AU742887B2 (en) | Hydrocarbon edge detection using seismic amplitude | |
| Rentsch et al. | Fast location of seismicity: A migration-type approach with application to hydraulic-fracturing data | |
| US20080130411A1 (en) | Seismic imaging with natural green's functions derived from vsp data | |
| EP3341759B1 (en) | Optimal survey design | |
| RU2339056C2 (en) | Generic three-dimension prediction of multiple events from surface | |
| CN109765615A (en) | A kind of inversion method for stratigraphic quality factor and device | |
| Brodic et al. | Three-component seismic land streamer study of an esker architecture through S-and surface-wave imaging | |
| AU2011380936B2 (en) | Seismic imaging systems and methods employing correlation-based stacking | |
| EP3341757B1 (en) | Nodal hybrid gather | |
| CN106443791B (en) | The method for asking for tilted stratum or anisotropic formation shear wave Value of residual static correction | |
| Li et al. | Gaussian beam imaging of fractures near the wellbore using sonic logging tools after removing dispersive borehole waves | |
| CN111077577B (en) | Well-ground combined reservoir description method and device | |
| Robinson | Velocity structure of the Wellington region, New Zealand, from local earthquake data and its implications for subduction tectonics | |
| EA025952B1 (en) | Extracting sv shear data from p-wave seismic data | |
| CN107422375B (en) | The determination method and apparatus of the CRP gather of subsea node | |
| CN115774288A (en) | Shear wave vector inversion modeling method and device | |
| Cosma et al. | VSP in crystalline rocks—From downhole velocity profiling to 3-D fracture mapping | |
| US5629905A (en) | Method of datuming seismic data and method of processing seismic data | |
| CN110869813A (en) | Method for acquiring a seismic dataset over a region of interest and related system | |
| CN111474580B (en) | Azimuth angle gather extraction method and system based on offset vector piece | |
| US4682307A (en) | Underwater seismic testing |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination |