CN115576008A

CN115576008A - Shear wave splitting correction method, device, equipment and storage medium for pure shear waves

Info

Publication number: CN115576008A
Application number: CN202110757354.0A
Authority: CN
Inventors: 岳媛媛; 张少华; 钱忠平; 杨军; 邹振; 张文波
Original assignee: China National Petroleum Corp; BGP Inc
Current assignee: China National Petroleum Corp; BGP Inc
Priority date: 2021-07-05
Filing date: 2021-07-05
Publication date: 2023-01-06

Abstract

The disclosure provides a shear wave splitting correction method, a device, equipment and a storage medium of pure shear waves, and belongs to the technical field of seismic exploration. The method comprises the following steps: acquiring crack direction information of a stratum of a work area to be detected, and acquiring first pure transverse wave data received by a target detector in the work area; when the crack direction information is used for indicating that the crack direction is a fixed direction, determining a fast transverse wave and a slow transverse wave corresponding to the first pure transverse wave data according to the crack direction information; determining the time difference between the fast transverse wave and the slow transverse wave, and correcting the slow transverse wave according to the time difference to obtain a corrected slow transverse wave; and determining a first corrected transverse wave field corresponding to the first pure transverse wave data according to the crack direction information, the fast transverse wave and the corrected slow transverse wave. The slow transverse wave is corrected through the time difference, and the transverse wave field with the azimuth anisotropy removed is obtained, so that the influence of the anisotropy is eliminated, and the accuracy of the determined first corrected transverse wave field is improved.

Description

Shear wave splitting correction method, device, equipment and storage medium for pure shear waves

Technical Field

The application relates to the technical field of seismic exploration, in particular to a shear wave splitting correction method, a device, equipment and a storage medium for pure shear waves.

Background

At present, the pure shear wave splitting technology is more and more widely applied to the field of seismic exploration. For example, the structural imaging, oil gas prediction, crack detection and the like of the work area to be detected are realized through a pure shear wave splitting technology. In the process of predicting oil gas of a work area to be measured by a pure shear wave splitting technology, when a stratum medium of the work area to be measured is an anisotropic medium, splitting correction needs to be performed on pure shear waves sometimes in order to determine a shear wave field corresponding to the pure shear waves.

Disclosure of Invention

The embodiment of the application provides a shear wave splitting correction method, a shear wave splitting correction device, shear wave splitting correction equipment and a storage medium for pure shear waves, and can improve the accuracy of a shear wave field obtained after shear wave splitting correction. The technical scheme is as follows:

in one aspect, the present application provides a shear wave splitting correction method for pure shear waves, including:

acquiring crack direction information of a stratum of a work area to be detected, and acquiring first pure transverse wave data received by a target detector in the work area;

when the crack direction information is used for indicating that the crack direction is a fixed direction, determining a fast transverse wave and a slow transverse wave corresponding to the first pure transverse wave data according to the crack direction information, wherein the fast transverse wave is a transverse wave parallel to the crack direction, and the slow transverse wave is a transverse wave perpendicular to the crack direction;

determining the time difference between the fast transverse wave and the slow transverse wave, and correcting the slow transverse wave according to the time difference to obtain a corrected slow transverse wave;

and determining a first corrected shear wave field corresponding to the first pure shear wave data according to the crack direction information, the fast shear wave and the corrected slow shear wave.

In a possible implementation manner, the determining, according to the fracture direction information, a fast shear wave and a slow shear wave corresponding to the first pure shear wave data includes:

acquiring a shot-and-test connecting line direction between a target seismic source and the target detector, and determining a first angle difference between the crack direction and the shot-and-test connecting line direction, wherein the target seismic source is a seismic source for exciting the first pure transverse wave data in the work area;

determining a first target pure transverse wave and a second target pure transverse wave;

and according to the first angle difference, performing two-component rotation on the first target pure transverse wave and the second target pure transverse wave to obtain a fast transverse wave and a slow transverse wave corresponding to the first pure transverse wave data.

In one possible implementation, the first pure shear wave data includes a first pure shear wave and a second pure shear wave;

the determining the first target pure shear wave and the second target pure shear wave comprises:

acquiring the line measuring directions of a plurality of detectors in the work area, and determining a second angle difference between the shot detection connecting line direction and the line measuring direction;

determining a first component wave and a second component wave corresponding to the first pure transverse wave, and a third component wave and a fourth component wave corresponding to the second pure transverse wave, wherein the first component wave is a component wave of the first pure transverse wave along the line measuring direction, the second component wave is a component wave of the first pure transverse wave along the direction perpendicular to the line measuring direction, the third component wave is a component wave of the second pure transverse wave along the line measuring direction, and the fourth component wave along the direction perpendicular to the line measuring direction;

and determining a first target pure shear wave and a second target pure shear wave according to the second angle difference, the first component wave, the second component wave, the third component wave and the fourth component wave.

In one possible implementation, the determining a time difference between the fast shear wave and the slow shear wave includes:

determining a fast transverse wave profile corresponding to the fast transverse waves and a slow transverse wave profile corresponding to the slow transverse waves, wherein the fast transverse wave profile comprises a plurality of first horizon identifications and horizon times of first horizon waves, and the slow transverse wave profile comprises a plurality of second horizon identifications and horizon times of second horizon waves;

determining a time difference between each first level wave and a second level wave corresponding to the first level wave according to the level identification and the level time corresponding to the first level waves and the level identification and the level time corresponding to the second level waves to obtain a time difference between the fast transverse wave and the slow transverse wave;

and the depth of the second level wave corresponding to the first level wave is the same as the depth of the first level wave.

In one possible implementation, the first pure shear wave data includes a first pure shear wave and a second pure shear wave; the determining a fast shear wave profile corresponding to the fast shear wave and a slow shear wave profile corresponding to the slow shear wave includes:

obtaining the line measuring directions of a plurality of detectors in the work area, and determining a third angle difference between the crack direction and the line measuring direction;

and determining a fast transverse wave section corresponding to the fast transverse wave and a slow transverse wave section corresponding to the slow transverse wave according to the third angle difference, the first component wave, the second component wave, the third component wave and the fourth component wave.

In one possible implementation, the first correction shear wave field comprises an SH shear wave field and an SV shear wave field; determining a first corrected shear wave field corresponding to the first pure shear wave data according to the fracture direction information, the fast shear wave and the corrected slow shear wave, including:

for any target seismic source in the work area, acquiring a shot-and-test connecting line direction between the target seismic source and the target detector, and determining a first angle difference between the crack direction and the shot-and-test connecting line direction, wherein the target seismic source is a seismic source which excites the first pure shear wave data in the work area;

and according to the first angle difference, performing reverse rotation on the fast transverse wave and the corrected slow transverse wave to obtain an SH transverse wave field corresponding to the first pure transverse wave data and an SV transverse wave field corresponding to the first pure transverse wave data.

In one possible implementation, the method further includes:

when the crack direction information is used for representing that the crack direction changes along with the depth, determining a plurality of crack time windows corresponding to the crack direction information;

acquiring second pure shear wave data corresponding to the plurality of fracture time windows, and a target time difference and a target fracture direction corresponding to each fracture time window;

and performing shear wave splitting correction on the second pure shear wave data according to the target time difference and the target fracture direction corresponding to each fracture time window to obtain a second corrected shear wave field corresponding to the second pure shear wave data.

In another aspect, the present application provides an apparatus for determining splitting wave information of a shear wave, the apparatus including:

the acquisition module is used for acquiring the crack direction information of the stratum of a work area to be detected and acquiring first pure transverse wave data received by a target detector in the work area;

a first determining module, configured to determine, according to the fracture direction information, a fast transverse wave and a slow transverse wave corresponding to the first pure transverse wave data when the fracture direction information is used to indicate that the fracture direction is a fixed direction, where the fast transverse wave is a transverse wave parallel to the fracture direction, and the slow transverse wave is a transverse wave perpendicular to the fracture direction;

the second determining module is used for determining the time difference between the fast transverse wave and the slow transverse wave, and correcting the slow transverse wave according to the time difference to obtain a corrected slow transverse wave;

and the third determining module is used for determining a first corrected shear wave field corresponding to the first pure shear wave data according to the crack direction information, the fast shear wave and the corrected slow shear wave.

In one possible implementation manner, the first determining module includes:

the acquiring unit is used for acquiring a shot-geophone link direction between a target seismic source and the target detector, and determining a first angle difference between the crack direction and the shot-geophone link direction, wherein the target seismic source is a seismic source for exciting the first pure transverse wave data in the work area;

the first determining unit is used for determining a first target pure transverse wave and a second target pure transverse wave;

and the rotating unit is used for performing two-component rotation on the first target pure transverse wave and the second target pure transverse wave according to the first angle difference to obtain a fast transverse wave and a slow transverse wave obtained by splitting the first pure transverse wave data along the crack.

In another possible implementation, the first pure shear wave data includes a first pure shear wave and a second pure shear wave;

the first determining unit is used for acquiring the line measuring directions of a plurality of detectors in the work area and determining a second angle difference between the shot detection connecting line direction and the line measuring direction; determining a first component wave and a second component wave corresponding to the first pure transverse wave, and a third component wave and a fourth component wave corresponding to the second pure transverse wave, where the first component wave is a component wave of the first pure transverse wave along the direction of the measuring line, the second component wave is a component wave of the first pure transverse wave along a direction perpendicular to the direction of the measuring line, the third component wave is a component wave of the second pure transverse wave along the direction of the measuring line, and the fourth component wave is a component wave of the second pure transverse wave along a direction perpendicular to the direction of the measuring line; and determining a first target pure transverse wave and a second target pure transverse wave according to the second angle difference, the first component wave, the second component wave, the third component wave and the fourth component wave.

In another possible implementation manner, the second determining module includes:

a second determining unit, configured to determine a fast shear wave profile corresponding to the fast shear wave and a slow shear wave profile corresponding to the slow shear wave, where the fast shear wave profile includes a plurality of first horizon identifiers and horizon times of the first horizon waves, and the slow shear wave profile includes horizon identifiers and horizon times of the plurality of second horizon waves;

a third determining unit, configured to determine a time difference between each first level wave and a second level wave corresponding to the first level wave according to the level identifier and the level time corresponding to the plurality of first level waves and the level identifier and the level time corresponding to the plurality of second level waves, so as to obtain a time difference between the fast transverse wave and the slow transverse wave;

the second determining unit is used for acquiring the line measuring directions of a plurality of detectors in the work area and determining a third angle difference between the crack direction and the line measuring directions; determining a first component wave and a second component wave corresponding to the first pure transverse wave, and a third component wave and a fourth component wave corresponding to the second pure transverse wave, where the first component wave is a component wave of the first pure transverse wave along the direction of the measuring line, the second component wave is a component wave of the first pure transverse wave along a direction perpendicular to the direction of the measuring line, the third component wave is a component wave of the second pure transverse wave along the direction of the measuring line, and the fourth component wave is a component wave of the second pure transverse wave along a direction perpendicular to the direction of the measuring line; and determining a fast transverse wave section corresponding to the fast transverse wave and a slow transverse wave section corresponding to the slow transverse wave according to the third angle difference, the first component wave, the second component wave, the third component wave and the fourth component wave.

In one possible implementation, the first correction shear wave wavefield includes an SH shear wave wavefield and an SV shear wave wavefield; the third determining module is configured to acquire, for any target seismic source in the work area, a shot-geophone link direction between the target seismic source and the target geophone, and determine a first angle difference between the fracture direction and the shot-geophone link direction, where the target seismic source is a seismic source that excites the first pure shear wave data in the work area; and according to the first angle difference, performing reverse rotation on the fast shear wave and the corrected slow shear wave to obtain an SH shear wave field corresponding to the first pure shear wave data and an SV shear wave field corresponding to the first pure shear wave data.

In one possible implementation, the apparatus further includes:

the fourth determining module is used for determining a plurality of crack time windows corresponding to the crack direction information when the crack direction information is used for representing that the crack direction changes along with the depth; acquiring second pure shear wave data corresponding to the plurality of fracture time windows, and a target time difference and a target fracture direction corresponding to each fracture time window; and performing shear wave splitting correction on the second pure shear wave data according to the target time difference and the target fracture direction corresponding to each fracture time window to obtain a second corrected shear wave field corresponding to the second pure shear wave data.

In another aspect, an embodiment of the present application provides a computer device, where the computer device includes: a processor and a memory, the memory having stored therein at least one program code, the at least one program code being loaded into and executed by the processor to implement the operations performed in the shear-wave-only shear wave splitting correction method of any of the possible implementations described above.

In another aspect, an embodiment of the present application provides a computer-readable storage medium, where at least one program code is stored, and the at least one program code is loaded and executed by a processor to implement the operations performed in the shear wave splitting correction method for shear wave only according to any one of the foregoing possible implementation manners.

The technical scheme provided by the embodiment of the application has the beneficial effects that at least:

the embodiment of the application provides a shear wave splitting correction method for pure shear waves, and due to the fact that for the stratum of an anisotropic medium, slow shear waves are corrected through the time difference between fast shear waves and slow shear waves, the influence of azimuth anisotropy is eliminated, and then through carrying out reverse rotation on the fast shear waves and the corrected slow shear waves, the accuracy of an obtained shear wave field is high.

Drawings

In order to more clearly illustrate the technical solutions in the embodiments of the present application, the drawings needed to be used in the description of the embodiments are briefly introduced below, and it is obvious that the drawings in the following description are only some embodiments of the present application, and it is obvious for those skilled in the art to obtain other drawings based on these drawings without creative efforts.

FIG. 1 is a flow chart illustrating a method of shear wave splitting correction of pure shear waves in accordance with an exemplary embodiment;

FIG. 2 is a schematic diagram illustrating four-component data of a pure shear wave in accordance with an exemplary embodiment;

FIG. 3 is a schematic illustration of a rotated quarter-wave transverse data according to an exemplary embodiment;

FIG. 4 is a schematic diagram illustrating a fast shear wave profile and a slow shear wave profile in accordance with an exemplary embodiment;

FIG. 5 is a schematic diagram of an SH and SV shear wave fields in accordance with an exemplary embodiment;

FIG. 6 is a block diagram illustrating a shear wave splitting correction apparatus for pure shear waves in accordance with an exemplary embodiment;

FIG. 7 is a block diagram illustrating a shear wave splitting correction apparatus for pure shear waves in accordance with an exemplary embodiment;

FIG. 8 is a block diagram illustrating a configuration of a computer device, according to an example embodiment.

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.

FIG. 1 is a flow chart illustrating a method of pure shear wave splitting correction in accordance with an exemplary embodiment.

Referring to fig. 1, the method includes:

101. the computer equipment obtains the crack direction information of the stratum of the work area to be detected and obtains first pure transverse wave data received by a target detector in the work area.

In the present embodiment, fracture direction information of the formation has been obtained. Optionally, the crack direction information is determined by pure shear wave splitting analysis methods such as alford rotation.

In one possible implementation, the direction of the fracture may be fixed with depth, with the fracture direction being the same at different depths in the formation. At this time, the computer device may perform shear wave splitting correction of pure shear waves through steps 102 to 104.

In another possible implementation, the fracture direction is different at different depths of the formation, and the fracture direction changes with the depth. At this time, the computer device may perform shear wave splitting correction of pure shear waves through steps 105 to 106.

It should be noted that a plurality of seismic sources and a plurality of detectors are distributed in the work area to be measured. The target detector is a detector for receiving first pure transverse wave data excited by the target seismic source. Wherein the target source may be any one of a plurality of sources. Optionally, the target seismic source is a controllable seismic source excited in the horizontal direction, and the seismic wave signal excited by the target seismic source is a pure transverse wave mainly including a transverse wave, and is excited once in the parallel line direction and the vertical line direction respectively. That is, the target detector receives two pure transverse waves. Correspondingly, the first pure transverse wave data comprises a first pure transverse wave and a second pure transverse wave, wherein the first pure transverse wave is a pure transverse wave which is received by the target detector and is excited by the target seismic source along a parallel measuring line direction, and the second pure transverse wave is a pure transverse wave which is excited by the target seismic source along a vertical measuring line direction.

In one possible implementation, the target detector is a three-component detector having two components in the horizontal direction, namely an x-component and a y-component, wherein the x-component and the y-component are perpendicular to each other. There is a component in the vertical direction, the z component, which is perpendicular to both the x and y components.

The two horizontal components of the three-component detector are respectively parallel to the line measuring direction and perpendicular to the line measuring direction, so that a component wave of a first pure transverse wave parallel to the line measuring direction, a component wave perpendicular to the line measuring direction, a component wave of a second pure transverse wave parallel to the line measuring direction, and a component wave perpendicular to the line measuring direction are received, namely, initial four-component data corresponding to the first pure transverse wave data.

102. When the crack direction information is used for indicating that the crack direction is a fixed direction along with the depth, the computer equipment determines a fast transverse wave and a slow transverse wave corresponding to the first pure transverse wave data according to the crack direction information, wherein the fast transverse wave is a transverse wave parallel to the crack direction, and the slow transverse wave is a transverse wave perpendicular to the crack direction.

In one possible implementation, the present step may include the following steps (1) to (3):

(1) The computer equipment acquires the shot and inspection connecting line direction between the target seismic source and the target detector and determines a first angle difference between the crack direction and the shot and inspection connecting line direction.

The target seismic source is a seismic source for exciting first pure transverse wave data in a work area; the target detector is a detector for receiving first pure transverse wave data excited by the target seismic source. Wherein the crack direction may be represented by the letter θ _fracture Represents; the shot-examination connecting line direction can be represented by the letter theta _azimuth Represents; the first angle difference between the crack direction and the shot-checking connecting line direction is as follows: theta _fracture -θ _azimuth . Can be used forAlternatively, the first angular difference may be represented by the letter γ; that is, γ = θ _fracture -θ _azimuth 。

In this step, the shot link direction is the direction in which the seismic source position points to the detector position. Alternatively, the shot-to-shot link direction may be referred to as the radial direction, generally denoted by the letter r; the vertical direction of the shot-to-shot link may be referred to as the tangential direction, generally indicated by the letter t.

(2) The computer device determines a first target pure shear wave and a second target pure shear wave.

In one possible implementation, the first pure shear wave data includes a first pure shear wave and a second pure shear wave; accordingly, the step of the computer device determining the first target pure shear wave and the second target pure shear wave may include the steps (21) to (23) of:

(21) And the computer equipment acquires the line measuring directions of a plurality of detectors in the work area and determines a second angle difference between the shot detection line connecting direction and the line measuring direction.

Wherein, the line measuring direction can be represented by the letter theta _inline And showing that a second angle difference between the shot-check connecting line direction and the survey line direction is as follows: theta _azimuth -θ _inline . Alternatively, the second angular difference may be represented by the letter β; that is, β = θ _azimuth -θ _inline 。

(22) The computer device determines a first component wave and a second component wave corresponding to the first pure shear wave, and a third component wave and a fourth component wave corresponding to the second pure shear wave.

It should be noted that the target detector is a three-component detector, and the three-component detector has two components in the horizontal direction, namely an x component and a y component, wherein the x component and the y component are perpendicular to each other. There is a component in the vertical direction, the z component, which is perpendicular to both the x and y components.

In one possible implementation, the x-component is parallel to the line direction and the y-component is perpendicular to the line direction. The target seismic source respectively excites the pure transverse wave once along the direction of the survey line and the direction vertical to the survey line. The first component wave is a component wave of a first pure transverse wave received by the target detector along the line measuring direction, and the second component wave is a component wave of a first pure transverse wave received by the target detector along the direction perpendicular to the line measuring direction. The third component wave is a component wave of a second pure transverse wave received by the target detector along the line measuring direction, and the fourth component wave is a component wave of the second pure transverse wave received by the target detector along the direction perpendicular to the line measuring direction.

In one possible implementation, the first component wave may be S _x R _x Indicating that the second component wave can be S _x R _y And (4) showing. The third component wave may be S _y R _x The fourth component wave can be represented by S _y R _y And (4) showing. Optionally, referring to fig. 2, the initial quarter data corresponding to the first pure shear wave data includes a first component wave, a second component wave, a third component wave, and a fourth component wave.

(23) And the computer equipment determines a first target pure transverse wave and a second target pure transverse wave according to the second angle difference, the first component wave, the second component wave, the third component wave and the fourth component wave.

In one possible implementation, the first target pure shear wave comprises a first radial pure shear wave and a first tangential pure shear wave; the second target pure shear wave comprises a second radial pure shear wave and a second tangential pure shear wave; and the computer equipment performs four-component rotation on the initial four-component data of the pure transverse wave according to a second angle difference between the shot-detection connecting line direction and the measuring line direction to obtain the rotated four-component transverse wave data. Correspondingly, the method comprises the following steps: the computer equipment determines a first radial pure transverse wave, a first tangential pure transverse wave, a second radial pure transverse wave and a second tangential pure transverse wave according to the second angle difference, the first component wave, the second component wave, the third component wave and the fourth component wave through a formula I;

the formula I is as follows:

wherein S is _x R _x Representing the first component wave, S _x R _y To representSecond component wave, S _y R _x Represents the third component wave, S _y R _y Representing a fourth component wave, beta representing a second angular difference, S _t R _r Representing the first radial pure transverse wave, S _t R _t Representing the first tangential pure transverse wave, S _t R _r Representing the second radial pure transverse wave, S _t R _t Representing a second tangential pure shear wave.

Optionally, referring to fig. 3, the rotated quarter-component shear wave data includes a first radial pure shear wave S _t R _r The first tangential pure transverse wave is S _t R _t Second radial pure transverse wave S _r R _r And a second tangential pure shear wave S _r R _t 。

(3) And the computer equipment performs two-component rotation on the first target pure shear wave and the second target pure shear wave according to the first angle difference to obtain a fast shear wave and a slow shear wave corresponding to the first pure shear wave data.

The first target pure shear wave comprises a first radial pure shear wave and a first tangential pure shear wave; the second target pure shear wave comprises a second radial pure shear wave and a second tangential pure shear wave. Correspondingly, the method comprises the following steps: the computer equipment performs two-component rotation on the first radial pure shear wave and the first tangential pure shear wave to obtain a first fast shear wave and a first slow shear wave corresponding to the first target pure shear wave; and the computer equipment performs two-component rotation on the second radial pure transverse wave and the second tangential pure transverse wave to obtain a second fast transverse wave and a second slow transverse wave corresponding to the second target pure transverse wave.

In one possible implementation, the two-component rotation of the first radial pure shear wave and the first tangential pure shear wave is achieved by rotating the receiving direction of the detector to be parallel to the crack direction and perpendicular to the crack direction. Correspondingly, the step of performing two-component rotation on the first radial pure shear wave and the first tangential pure shear wave by the computer device is as follows: according to the first angle difference, the computer equipment performs two-component rotation on the first radial pure transverse wave and the first tangential pure transverse wave through the following formula II to obtain a first fast transverse wave and a first slow transverse wave corresponding to the first target pure transverse wave;

the formula II is as follows:

wherein S is _t R _s1 Representing the first fast transverse wave, S _t R _s2 Representing the first slow transverse wave, S _t R _r Representing the first radial pure transverse wave, S _t R _t Representing a first tangential pure shear wave and gamma a first angular difference.

The step of performing two-component rotation on the second radial pure transverse wave and the second tangential pure transverse wave by the computer equipment comprises the following steps: according to the first angle difference and through the following formula III, the computer equipment performs two-component rotation on the second radial pure transverse wave and the second tangential pure transverse wave to obtain a second fast transverse wave and a second slow transverse wave corresponding to the second target pure transverse wave;

the formula III is as follows:

wherein S is _r R _s1 Representing the second fast transverse wave, S _r R _s2 Representing the second slow transverse wave, S _r R _r Representing the second radial pure transverse wave, S _r R _t Representing a second tangential pure shear wave and gamma representing a first angular difference.

103. And the computer equipment determines the time difference between the fast transverse wave and the slow transverse wave, and corrects the slow transverse wave according to the time difference to obtain the corrected slow transverse wave.

In one possible implementation, the computer device determining the time difference of the fast shear wave and the slow shear wave may be implemented by the following steps (1) to (2):

(1) The computer equipment determines a fast transverse wave section corresponding to the fast transverse wave and a slow transverse wave section corresponding to the slow transverse wave, wherein the fast transverse wave section comprises a plurality of first horizon marks and horizon time of the first horizon waves corresponding to the fast transverse wave in a preset time period, and the slow transverse wave section comprises a plurality of second horizon marks and horizon time of the second horizon waves corresponding to the slow transverse wave in the preset time period.

In one possible implementation, the first pure shear wave data includes a first pure shear wave and a second pure shear wave. Correspondingly, the method comprises the following steps: the computer equipment acquires the direction of the survey line of the work area and determines a third angle difference between the direction of the crack and the direction of the survey line; determining a first component wave and a second component wave corresponding to the first pure transverse wave, and a third component wave and a fourth component wave corresponding to the second pure transverse wave, wherein the first component wave is a component wave of the first pure transverse wave along the measuring line direction, the second component wave is a component wave of the first pure transverse wave along the perpendicular direction of the measuring line direction, the third component wave is a component wave of the second pure transverse wave along the measuring line direction, and the fourth component wave along the perpendicular direction of the second pure transverse wave along the measuring line direction; and determining a fast transverse wave section corresponding to the fast transverse wave and a slow transverse wave section corresponding to the slow transverse wave according to the third angle difference, the first component wave, the second component wave, the third component wave and the fourth component wave.

In a possible implementation manner, the first component wave, the second component wave, the third component wave and the fourth component wave are initial four-component data corresponding to the first pure shear wave data, and the computer device performs fast-slow wave separation on the four-component data through a rotation formula. Correspondingly, the step of determining, by the computer device, a fast transverse wave profile corresponding to the fast transverse wave and a slow transverse wave profile corresponding to the slow transverse wave according to the third angle difference, the first component wave, the second component wave, the third component wave and the fourth component wave is as follows: the computer equipment determines a fast transverse wave section corresponding to the fast transverse wave and a slow transverse wave section corresponding to the slow transverse wave according to the third angle difference, the first component wave, the second component wave, the third component wave and the fourth component wave through the following formula IV;

the formula four is as follows:

wherein S is _x R _x Representing the first component wave, S _x R _y Representing the second component wave, S _y R _x Represents the third component wave, S _y R _y A fourth component wave is shown which is,

represents a third angle difference, S _S1 R _S1 Showing the fast transverse wave profile, S _S2 R _S2 Showing a slow shear wave profile.

Wherein the crack direction may be represented by the letter θ _fracture The direction of the measuring line can be indicated by the letter theta _inline A third angular difference between the crack direction and the gage line direction is shown as: theta.theta. _fracture -θ _inline . Alternatively, the third angular difference may be in the form of letters

Represents; that is to say that the first and second electrodes,

the fast shear wave section corresponds to the section information of the fast shear wave, and the slow shear wave section corresponds to the section information of the slow shear wave. Optionally, referring to fig. 4, the computer device performs fast-slow wave separation on the shear wave initial four-component data through a rotation formula, and the obtained profile information of the fast shear wave and the profile information of the slow shear wave are shown in fig. 4, where the fast shear wave profile is a leftmost shear wave profile in fig. 4, and the slow shear wave profile is a rightmost shear wave profile in fig. 4.

(2) And the computer equipment determines the time difference between each first layer wave and the second layer wave corresponding to the first layer wave according to the layer identification and the layer time corresponding to the first layer waves of the fast transverse wave profile and the layer identification and the layer time corresponding to the second layer waves of the slow transverse wave profile, so as to obtain the time difference between the fast transverse wave and the slow transverse wave.

The level marks of a plurality of first level waves of the fast transverse wave profile and the level marks of a plurality of second level waves of the slow transverse wave profile correspond to each other from shallow to deep one by one. For example, the horizon of a plurality of first horizon waves of a fast shear wave profile is identified as h ₁ s ₁ 、h ₂ s ₁ 、h ₃ s ₁ The level of a plurality of second level waves of the slow transverse wave profile is marked as h ₁ s ₂ 、h ₂ s ₂ 、h ₃ s ₂ (ii) a Wherein h is ₁ s ₁ 、h ₂ s ₁ 、h ₃ s ₁ And h ₁ s ₂ 、h ₂ s ₂ 、h ₃ s ₂ One-to-one correspondence from shallow to deep, i.e., h ₁ s ₁ Corresponds to h ₁ s ₂ ，h ₂ s ₁ Corresponds to h ₂ s ₂ ，h ₃ s ₁ Corresponds to h ₃ s ₂ 。

In one possible implementation, the method includes the following steps: according to the sequence from shallow to deep, the computer equipment determines the time difference between a first level wave and a first second level wave according to the level time of a first level wave of a fast transverse wave section and the level time of a first second level wave of a slow transverse wave section, determines the time difference between a second first level wave and a second level wave according to the level time of a second first level wave of the fast transverse wave section and the level time of a second level wave of the slow transverse wave section, and determines the time difference between each first level wave and the second level wave corresponding to the first level wave from shallow to deep in sequence to obtain the time difference between the fast transverse wave and the slow transverse wave.

In step 103, the computer device corrects the slow transverse wave according to the time difference, and the step of obtaining the corrected slow transverse wave is as follows: the computer equipment corrects a first slow transverse wave corresponding to the first target pure transverse wave according to the time difference and a formula V to obtain a first corrected slow transverse wave; according to the time difference, correcting a second slow transverse wave corresponding to a second target pure transverse wave through a formula six to obtain a second corrected slow transverse wave;

the formula five is as follows: s _t R _s2up (t)＝S _t R _s2 (t+Δt)

Formula six: s _r R _s2up (t)＝S _r R _s2 (t+Δt)

Wherein S is _t R _s2up (t) denotes the first corrective slow transverse wave, S, corresponding to time t _r R _s2up (t) represents a second corrected slow transverse wave corresponding to time t, Δ t represents a time difference, S _t R _s2 (t + Δ t) represents a first slow transverse wave corresponding to the time t + Δ t, S _r R _s2 (t + Δ t) represents a second slow transverse wave corresponding to the time t + Δ t. Alternatively, t may represent a time point or a time period.

It should be noted that the time difference may vary with the depth, that is, the time difference between the fast and slow waves corresponding to different levels may be different, and at this time, for the pure transverse wave between two levels, the time difference between the fast and slow waves may be the interpolation result of the time difference between the previous level and the time difference between the next level.

It should be noted that the slow transverse wave profile is an amplitude value corresponding to the slow transverse wave, where the amplitude value is obtained based on first pure transverse wave data periodically collected by the target detector. In a possible implementation manner, in the process of performing time difference correction on the slow wave, when t + Δ t has no corresponding amplitude value, the computer device performs interpolation according to the sinc function to determine the amplitude value corresponding to t + Δ t.

In the embodiment of the application, when the time after the slow transverse wave time difference correction does not have a corresponding amplitude value, the computer equipment can determine the amplitude value corresponding to the time through sinc function interpolation, so that the accuracy of the time difference correction is improved.

104. And the computer equipment determines a first corrected transverse wave field corresponding to the first pure transverse wave data according to the crack direction information, the fast transverse wave and the corrected slow transverse wave.

In one possible implementation, the first correction shear wave wavefield includes an SH shear wave wavefield and an SV shear wave wavefield. Correspondingly, the method comprises the following steps: the computer equipment acquires the shot-examination connecting line direction between a target seismic source and a target detector for any target seismic source in a work area, determines a first angle difference between the crack direction and the shot-examination connecting line direction, and the target seismic source is a seismic source which excites the first pure transverse wave data in the work area; and according to the first angle difference, performing reverse rotation on the fast transverse wave and the correction slow transverse wave to obtain an SH transverse wave field corresponding to the first pure transverse wave data and an SV transverse wave field corresponding to the first pure transverse wave data.

In one possible implementation, the fast shear waves include a first fast shear wave corresponding to a first target pure shear wave and a second fast shear wave corresponding to a second target pure shear wave; the corrected slow shear waves comprise first corrected slow shear waves corresponding to the first target pure shear waves and second corrected slow shear waves corresponding to the second target pure shear waves. Correspondingly, the step that the computer equipment reversely rotates the fast transverse wave and the correction slow transverse wave according to the first angle difference to obtain the SH transverse wave field corresponding to the first pure transverse wave data and the SV transverse wave field corresponding to the first pure transverse wave data is as follows: the computer equipment reversely rotates a first fast transverse wave and a first correction slow transverse wave corresponding to the first target pure transverse wave according to the first angle difference and the following formula seven to obtain an SH transverse wave field corresponding to the first pure transverse wave data; according to the first angle difference, carrying out reverse rotation on a second fast transverse wave and a second correction slow transverse wave corresponding to a second target pure transverse wave through the following formula eight to obtain an SV transverse wave field corresponding to the first pure transverse wave data;

the formula seven:

the formula eight:

wherein S is _t R _s1 Representing the first fast transverse wave, S _t R _s2up Representing the first corrected slow transverse wave, S _r R _s1 Representing the second fast transverse wave, S _r R _s2up Representing a second corrected slow transverse wave, S _t R _t ' denotes the SH shear wave field, S, corresponding to the first pure shear wave data _r R _r ' denotes the SV shear wave field to which the first pure shear wave data corresponds.

Note that S is _t R _r ' and S _t R _t ' correcting the post-shear wave field for the first target pure shear wave splitting, S _r R _r ' and S _r R _t ' denotes the second target pure shear wave splitting corrected shear wave field. Wherein S is _t R _t ' is the corrected SH wave field, S _t R _r ' effective energy will become 0, i.e., S _t R _t ' to eliminate the anisotropic SH shear wave field. Wherein S is _r R _r ' is the corrected SV wave field, S _r R _t ' the effective energy will become 0, i.e., S _r R _r ' to eliminate the anisotropic SV shear wave field.

In one possible implementation, see FIG. 5, the shear wave splitting corrected SV wave is shown in the shear wave wavefield on the left side of FIG. 5, i.e., S _r R _r '; the shear wave splitting corrected SH wave is shown in the shear wave field on the right side of FIG. 5, i.e., S _r R _t ', the middle two shear wave fields of FIG. 5 are S _r R _t ' and S _t R _r ', indicates that the energy has become 0.

It should be noted that, with continued reference to fig. 1, when the crack direction information is used to indicate that the crack direction changes with the depth, the computer device completes step 101, directly executes steps 105 to 107, and does not execute steps 102 to 104.

105. When the crack direction information is used to indicate that the crack direction changes with depth, the computer device determines a plurality of crack time windows corresponding to the crack direction information.

In a possible implementation manner, when the crack direction changes with the depth, the crack directions corresponding to different depths of the crack are different, and accordingly, the step of determining, by the computer device, a plurality of crack time windows corresponding to the crack direction information is as follows: the computer equipment determines a plurality of crack directions corresponding to the cracks under different depths, and determines a plurality of crack time windows according to the crack directions, wherein one crack direction corresponds to one crack time window.

106. And the computer equipment acquires second pure shear wave data corresponding to the plurality of fracture time windows, and a target time difference and a target fracture direction corresponding to each fracture time window.

For each fracture time window, the target time difference is the time difference between the split fast shear wave and the split slow shear wave when the second pure shear wave data passes through the fracture time window. The target crack direction is the crack direction below the crack time window.

In the embodiment of the application, the target time difference and the crack direction can be obtained by a worker through shear wave splitting analysis. Correspondingly, the method comprises the following steps: and the staff transmits the target time difference and the crack direction to the computer equipment, and the computer equipment acquires second pure shear wave data corresponding to the plurality of crack time windows and the target time difference and the target crack direction corresponding to each crack time window.

107. And the computer equipment performs shear wave splitting correction on the second pure shear wave data according to the target time difference and the target fracture direction corresponding to each fracture time window to obtain a second corrected shear wave field corresponding to the second pure shear wave data.

In one possible implementation, the second pure shear wave data includes a third pure shear wave and a fourth pure shear wave; the second corrected shear wave field comprises an SV wave field parallel to the shot-to-shot link direction and an SH wave field perpendicular to the shot-to-shot link direction. Correspondingly, the method comprises the following steps: according to the target time difference and the target crack direction corresponding to each crack time window, the computer equipment carries out shear wave splitting correction on the third pure shear wave through the following formula nine to obtain an SV wave field of the second pure shear wave data parallel to the shot-inspection connecting line direction, and according to the target time difference and the target crack direction corresponding to each crack time window, carries out shear wave splitting correction on the fourth pure shear wave through the following formula ten to obtain an SH wave field of the second pure shear wave data perpendicular to the shot-inspection connecting line direction;

the formula is nine:

U _0sr (ω)＝R ₁ ^-1 D ₁ …R _n-1 ^-1 D _n-1 R _n ^-1 D _n R _n D _n-1 R _n-1 …D ₁ R ₁ U _sr (ω)

formula ten:

U _0st (ω)＝R ₁ ^-1 D ₁ …R _n-1 ^-1 D _n-1 R _n ^-1 D _n R _n D _n-1 R _n-1 …D ₁ R ₁ U _st (ω)

where n represents the number of multiple fracture time windows, U _sr (ω) represents the third pure transverse wave, U _st (ω) represents the fourth pure transverse wave, U _0sr (omega) is the SV shear wave field, U, after shear wave splitting correction _0st (ω) represents the SH shear wave field, Δ t, after shear wave splitting correction _n Represents the target time difference, theta, corresponding to the nth crack window _n And representing the target crack direction corresponding to the nth crack time window, wherein alpha represents the shot-check connecting line direction.

Optionally, U _0sr (ω) is the SV shear wave field with the anisotropy removed. U shape _0st (ω) is the SH shear wave field with the anisotropy removed.

The embodiment of the application provides a shear wave splitting correction method for pure shear waves, and as for the stratum of an anisotropic medium, slow shear waves are corrected through the time difference between fast shear waves and slow shear waves, the influence of azimuth anisotropy is eliminated, and then the fast shear waves and the corrected slow shear waves are subjected to reverse rotation, so that the accuracy of an obtained shear wave field is high.

Fig. 6 is a block diagram illustrating a shear wave splitting correction apparatus for pure shear waves according to an exemplary embodiment. Referring to fig. 6, the apparatus includes:

the acquisition module 601 is configured to acquire crack direction information of a stratum of a work area to be detected, and acquire first pure transverse wave data received by a target detector in the work area;

a first determining module 602, configured to determine, according to the fracture direction information, a fast transverse wave and a slow transverse wave corresponding to the first pure transverse wave data when the fracture direction information is used to indicate that the fracture direction is a fixed direction, where the fast transverse wave is a transverse wave parallel to the fracture direction, and the slow transverse wave is a transverse wave perpendicular to the fracture direction;

a second determining module 603, configured to determine a time difference between the fast transverse wave and the slow transverse wave, and correct the slow transverse wave according to the time difference to obtain a corrected slow transverse wave;

and a third determining module 604, configured to determine a first corrected shear wave field corresponding to the first pure shear wave data according to the fracture direction information, the fast shear wave, and the corrected slow shear wave.

In one possible implementation, referring to fig. 7, the first determining module 602 includes:

the acquiring unit 6021 is configured to acquire a shot-geophone link direction between a target seismic source and a target geophone, and determine a first angle difference between a crack direction and the shot-geophone link direction, where the target seismic source is a seismic source for exciting first pure shear wave data in a work area;

a first determination unit 6022 for determining a first target pure shear wave and a second target pure shear wave;

and the rotating unit 6023 is configured to perform two-component rotation on the first target pure transverse wave and the second target pure transverse wave according to the first angle difference to obtain a fast transverse wave and a slow transverse wave corresponding to the first pure transverse wave data.

the first determining unit 6021 is configured to obtain line measurement directions of the plurality of detectors in the work area, and determine a second angle difference between the shot detection line connection direction and the line measurement direction; determining a first component wave and a second component wave corresponding to the first pure transverse wave, and a third component wave and a fourth component wave corresponding to the second pure transverse wave, wherein the first component wave is a component wave of the first pure transverse wave along the measuring line direction, the second component wave is a component wave of the first pure transverse wave along the perpendicular direction of the measuring line direction, the third component wave is a component wave of the second pure transverse wave along the measuring line direction, and the fourth component wave along the perpendicular direction of the second pure transverse wave along the measuring line direction; and determining a first target pure transverse wave and a second target pure transverse wave according to the second angle difference, the first component wave, the second component wave, the third component wave and the fourth component wave.

In another possible implementation, with continued reference to fig. 7, the second determining module 603 includes:

a second determining unit 6031 configured to determine a fast transverse wave profile corresponding to a fast transverse wave and a slow transverse wave profile corresponding to a slow transverse wave, where the fast transverse wave profile includes a plurality of level identifiers and level times of the first level waves, and the slow transverse wave profile includes a plurality of level identifiers and level times of the second level waves;

a third determining unit 6032, configured to determine, according to the level identifiers and the level times corresponding to the plurality of first level waves and the level identifiers and the level times corresponding to the plurality of second level waves, a time difference between each first level wave and a second level wave corresponding to the first level wave, so as to obtain a time difference between a fast transverse wave and a slow transverse wave;

a second determining unit 6031, configured to obtain line measurement directions of multiple detectors in the work area, and determine a third angle difference between the crack direction and the line measurement directions; determining a first component wave and a second component wave corresponding to the first pure transverse wave, and a third component wave and a fourth component wave corresponding to the second pure transverse wave, wherein the first component wave is a component wave of the first pure transverse wave along the measuring line direction, the second component wave is a component wave of the first pure transverse wave along the perpendicular direction of the measuring line direction, the third component wave is a component wave of the second pure transverse wave along the measuring line direction, and the fourth component wave along the perpendicular direction of the second pure transverse wave along the measuring line direction; and determining a fast transverse wave section corresponding to the fast transverse wave and a slow transverse wave section corresponding to the slow transverse wave according to the third angle difference, the first component wave, the second component wave, the third component wave and the fourth component wave.

In one possible implementation, the first correction shear wave wavefield includes an SH shear wave wavefield and an SV shear wave wavefield; a third determining module 604, configured to, for any target seismic source in the work area, obtain a shot-detection connection line direction between the target seismic source and a target detector, and determine a first angle difference between a crack direction and the shot-detection connection line direction, where the target seismic source is a seismic source that excites first pure shear wave data in the work area; and according to the first angle difference, performing reverse rotation on the fast transverse wave and the correction slow transverse wave to obtain an SH transverse wave field corresponding to the first pure transverse wave data and an SV transverse wave field corresponding to the first pure transverse wave data.

In one possible implementation, with continued reference to fig. 7, the apparatus further comprises:

a fourth determining module 605, configured to determine a plurality of fracture time windows corresponding to the fracture direction information when the fracture direction information is used to indicate that the fracture direction changes with the depth; acquiring second pure shear wave data corresponding to a plurality of fracture time windows, and a target time difference and a target fracture direction corresponding to each fracture time window; and performing shear wave splitting correction on the second pure shear wave data according to the target time difference and the target fracture direction corresponding to each fracture time window to obtain a second corrected shear wave field corresponding to the second pure shear wave data.

The embodiment of the application provides a shear wave splitting correcting unit of pure shear wave, because to the stratum of anisotropic medium, through the time difference of fast shear wave and slow shear wave, rectifies slow shear wave, has eliminated the anisotropic influence in position, and then through carrying out the derotation with fast shear wave and correction slow shear wave, the accuracy of the shear wave field that obtains is high.

Fig. 8 is a block diagram illustrating a terminal 800 according to an exemplary embodiment of the present invention. The terminal 800 may be: a smart phone, a tablet computer, an MP3 player (Moving Picture Experts Group Audio Layer III, motion Picture Experts compression standard Audio Layer 3), an MP4 player (Moving Picture Experts Group Audio Layer IV, motion Picture Experts compression standard Audio Layer 4), a notebook computer, or a desktop computer. The terminal 800 may also be referred to by other names such as user equipment, portable terminal, laptop terminal, desktop terminal, etc.

In general, the terminal 800 includes: a processor 801 and a memory 802.

The processor 801 may include one or more processing cores, such as a 4-core processor, an 8-core processor, and so forth. The processor 801 may be implemented in at least one hardware form of a DSP (Digital Signal Processing), an FPGA (Field-Programmable Gate Array), and a PLA (Programmable Logic Array). The processor 801 may also include a main processor and a coprocessor, where the main processor is a processor for Processing data in an awake state, and is also called a Central Processing Unit (CPU); a coprocessor is a low power processor for processing data in a standby state. In some embodiments, the processor 801 may be integrated with a GPU (Graphics Processing Unit), which is responsible for rendering and drawing the content required to be displayed on the display screen. In some embodiments, the processor 801 may further include an AI (Artificial Intelligence) processor for processing computing operations related to machine learning.

Memory 802 may include one or more computer-readable storage media, which may be non-transitory. Memory 802 may also include high speed random access memory, as well as non-volatile memory, such as one or more magnetic disk storage devices, flash memory storage devices. In some embodiments, a non-transitory computer-readable storage medium in memory 802 is used to store at least one instruction for execution by processor 801 to implement a shear-wave-only shear-wave splitting correction method provided by method embodiments herein.

In some embodiments, the terminal 800 may further include: a peripheral interface 803 and at least one peripheral. The processor 801, memory 802 and peripheral interface 803 may be connected by bus or signal lines. Various peripheral devices may be connected to the peripheral interface 803 by a bus, signal line, or circuit board. Specifically, the peripheral device includes: at least one of a radio frequency circuit 804, a display 805, a camera 806, an audio circuit 807, a positioning component 808, and a power supply 809.

The peripheral interface 803 may be used to connect at least one peripheral related to I/O (Input/Output) to the processor 801 and the memory 802. In some embodiments, the processor 801, memory 802, and peripheral interface 803 are integrated on the same chip or circuit board; in some other embodiments, any one or two of the processor 801, the memory 802, and the peripheral interface 803 may be implemented on separate chips or circuit boards, which are not limited by this embodiment.

The Radio Frequency circuit 804 is used for receiving and transmitting RF (Radio Frequency) signals, also called electromagnetic signals. The radio frequency circuitry 804 communicates with communication networks and other communication devices via electromagnetic signals. The rf circuit 804 converts an electrical signal into an electromagnetic signal to be transmitted, or converts a received electromagnetic signal into an electrical signal. Optionally, the radio frequency circuit 804 includes: an antenna system, an RF transceiver, one or more amplifiers, a tuner, an oscillator, a digital signal processor, a codec chipset, a subscriber identity module card, and so forth. The radio frequency circuit 804 may communicate with other terminals via at least one wireless communication protocol. The wireless communication protocols include, but are not limited to: metropolitan area networks, various generation mobile communication networks (2G, 3G, 4G, and 5G), wireless local area networks, and/or WiFi (Wireless Fidelity) networks. In some embodiments, the radio frequency circuit 804 may further include NFC (Near Field Communication) related circuits, which are not limited in this application.

The display screen 805 is used to display a UI (User Interface). The UI may include graphics, text, icons, video, and any combination thereof. When the display 805 is a touch display, the display 805 also has the ability to capture touch signals on or above the surface of the display 805. The touch signal may be input to the processor 801 as a control signal for processing. At this point, the display 805 may also be used to provide virtual buttons and/or a virtual keyboard, also referred to as soft buttons and/or a soft keyboard. In some embodiments, the display 805 may be one, providing the front panel of the terminal 800; in other embodiments, the display 805 may be at least two, respectively disposed on different surfaces of the terminal 800 or in a folded design; in still other embodiments, the display 805 may be a flexible display disposed on a curved surface or a folded surface of the terminal 800. Even further, the display 805 may be arranged in a non-rectangular irregular pattern, i.e., a shaped screen. The Display 805 can be made of LCD (Liquid Crystal Display), OLED (Organic Light-Emitting Diode), and other materials.

The camera assembly 806 is used to capture images or video. Optionally, camera assembly 806 includes a front camera and a rear camera. Generally, a front camera is disposed at a front panel of the terminal, and a rear camera is disposed at a rear surface of the terminal. In some embodiments, the number of the rear cameras is at least two, and each rear camera is any one of a main camera, a depth-of-field camera, a wide-angle camera and a telephoto camera, so that the main camera and the depth-of-field camera are fused to realize a background blurring function, and the main camera and the wide-angle camera are fused to realize panoramic shooting and VR (Virtual Reality) shooting functions or other fusion shooting functions. In some embodiments, camera assembly 806 may also include a flash. The flash lamp can be a monochrome temperature flash lamp or a bicolor temperature flash lamp. The double-color-temperature flash lamp is a combination of a warm-light flash lamp and a cold-light flash lamp, and can be used for light compensation at different color temperatures.

The audio circuit 807 may include a microphone and a speaker. The microphone is used for collecting sound waves of a user and the environment, converting the sound waves into electric signals, and inputting the electric signals to the processor 801 for processing or inputting the electric signals to the radio frequency circuit 804 to realize voice communication. For the purpose of stereo sound collection or noise reduction, a plurality of microphones may be provided at different portions of the terminal 800. The microphone may also be an array microphone or an omni-directional pick-up microphone. The speaker is used to convert electrical signals from the processor 801 or the radio frequency circuit 804 into sound waves. The loudspeaker can be a traditional film loudspeaker or a piezoelectric ceramic loudspeaker. When the speaker is a piezoelectric ceramic speaker, the speaker can be used for purposes such as converting an electric signal into a sound wave audible to a human being, or converting an electric signal into a sound wave inaudible to a human being to measure a distance. In some embodiments, the audio circuitry 807 may also include a headphone jack.

The positioning component 808 is used to locate the current geographic Location of the terminal 800 for navigation or LBS (Location Based Service). The Positioning component 808 may be a Positioning component based on the GPS (Global Positioning System) in the united states, the beidou System in china, the graves System in russia, or the galileo System in the european union.

Power supply 809 is used to provide power to various components in terminal 800. The power supply 809 can be ac, dc, disposable or rechargeable. When the power source 809 comprises a rechargeable battery, the rechargeable battery may support wired or wireless charging. The rechargeable battery may also be used to support fast charge technology.

In some embodiments, terminal 800 also includes one or more sensors 810. The one or more sensors 810 include, but are not limited to: acceleration sensor 811, gyro sensor 812, pressure sensor 813, fingerprint sensor 814, optical sensor 815 and proximity sensor 816.

The acceleration sensor 811 may detect the magnitude of acceleration in three coordinate axes of the coordinate system established with the terminal 800. For example, the acceleration sensor 811 may be used to detect the components of the gravitational acceleration in three coordinate axes. The processor 801 may control the display 805 to display the user interface in a landscape view or a portrait view according to the gravitational acceleration signal collected by the acceleration sensor 811. The acceleration sensor 811 may also be used for acquisition of motion data of a game or a user.

The gyro sensor 812 may detect a body direction and a rotation angle of the terminal 800, and the gyro sensor 812 may cooperate with the acceleration sensor 811 to acquire a 3D motion of the user with respect to the terminal 800. From the data collected by the gyro sensor 812, the processor 801 may implement the following functions: motion sensing (such as changing the UI according to a user's tilting operation), image stabilization at the time of photographing, game control, and inertial navigation.

Pressure sensors 813 may be disposed on the side frames of terminal 800 and/or underneath display 805. When the pressure sensor 813 is disposed on the side frame of the terminal 800, the holding signal of the user to the terminal 800 can be detected, and the processor 801 performs left-right hand recognition or shortcut operation according to the holding signal collected by the pressure sensor 813. When the pressure sensor 813 is disposed at a lower layer of the display screen 805, the processor 801 controls the operability control on the UI interface according to the pressure operation of the user on the display screen 805. The operability control comprises at least one of a button control, a scroll bar control, an icon control and a menu control.

The fingerprint sensor 814 is used for collecting a fingerprint of the user, and the processor 801 identifies the identity of the user according to the fingerprint collected by the fingerprint sensor 814, or the fingerprint sensor 814 identifies the identity of the user according to the collected fingerprint. Upon identifying that the user's identity is a trusted identity, the processor 801 authorizes the user to perform relevant sensitive operations including unlocking a screen, viewing encrypted information, downloading software, paying for and changing settings, etc. Fingerprint sensor 814 may be disposed on the front, back, or side of terminal 800. When a physical button or a vendor Logo is provided on the terminal 800, the fingerprint sensor 814 may be integrated with the physical button or the vendor Logo.

The optical sensor 815 is used to collect ambient light intensity. In one embodiment, processor 801 may control the display brightness of display 805 based on the ambient light intensity collected by optical sensor 815. Specifically, when the ambient light intensity is high, the display brightness of the display screen 805 is increased; when the ambient light intensity is low, the display brightness of the display 805 is reduced. In another embodiment, the processor 801 may also dynamically adjust the shooting parameters of the camera assembly 806 based on the ambient light intensity collected by the optical sensor 815.

A proximity sensor 816, also known as a distance sensor, is typically disposed on a front panel of the terminal 800. The proximity sensor 816 is used to collect the distance between the user and the front surface of the terminal 800. In one embodiment, when the proximity sensor 816 detects that the distance between the user and the front surface of the terminal 800 gradually decreases, the processor 801 controls the display 805 to switch from the bright screen state to the dark screen state; when the proximity sensor 816 detects that the distance between the user and the front surface of the terminal 800 becomes gradually larger, the display 805 is controlled by the processor 801 to switch from the breath-screen state to the bright-screen state.

Those skilled in the art will appreciate that the configuration shown in fig. 8 is not intended to be limiting of terminal 800 and may include more or fewer components than those shown, or some components may be combined, or a different arrangement of components may be used.

In an exemplary embodiment, a storage medium comprising program code, such as the memory 804 comprising program code, executable by the processor 820 of the apparatus 800 to perform the above-described method is also provided. Alternatively, the storage medium may be a non-transitory computer readable storage medium, which may be, for example, a ROM (Read-Only Memory), a RAM (Random Access Memory), a CD-ROM (Compact Disc Read-Only Memory), a magnetic tape, a floppy disk, an optical data storage device, and the like.

Other embodiments of the disclosure will be apparent to those skilled in the art from consideration of the specification and practice of the disclosure disclosed herein. This application is intended to cover any variations, uses, or adaptations of the disclosure following, in general, the principles of the disclosure and including such departures from the present disclosure as come within known or customary practice within the art to which the disclosure pertains. It is intended that the specification and examples be considered as exemplary only, with a true scope and spirit of the disclosure being indicated by the following claims.

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

1. A method for shear wave splitting correction of pure shear waves, the method comprising:

and determining a first correction transverse wave field corresponding to the first pure transverse wave data according to the crack direction information, the fast transverse wave and the correction slow transverse wave.

2. The method of claim 1, wherein the determining fast shear waves and slow shear waves corresponding to the first pure shear wave data according to the fracture direction information comprises:

3. The method of claim 2, wherein the first pure shear wave data comprises a first pure shear wave and a second pure shear wave;

the determining a first target pure shear wave and a second target pure shear wave includes:

and determining a first target pure transverse wave and a second target pure transverse wave according to the second angle difference, the first component wave, the second component wave, the third component wave and the fourth component wave.

4. The method of claim 1, wherein said determining a time difference between said fast shear wave and said slow shear wave comprises:

5. The method of claim 4, wherein the first pure shear wave data comprises a first pure shear wave and a second pure shear wave;

the determining of the fast shear wave profile corresponding to the fast shear wave and the slow shear wave profile corresponding to the slow shear wave includes:

6. The method of claim 1, wherein the first correction shear wave wavefield comprises an SH shear wave wavefield and an SV shear wave wavefield;

determining a first corrected shear wave field corresponding to the first pure shear wave data according to the fracture direction information, the fast shear wave and the corrected slow shear wave, including:

for any target seismic source in the work area, acquiring a shot-detection connecting line direction between the target seismic source and the target detector, and determining a first angle difference between the crack direction and the shot-detection connecting line direction, wherein the target seismic source is a seismic source which excites the first pure shear wave data in the work area;

and according to the first angle difference, performing reverse rotation on the fast shear wave and the corrected slow shear wave to obtain an SH shear wave field corresponding to the first pure shear wave data and an SV shear wave field corresponding to the first pure shear wave data.

7. The method of claim 1, further comprising:

8. A shear wave splitting correction device for pure shear waves, the device comprising:

and the third determining module is used for determining a first corrected transverse wave field corresponding to the first pure transverse wave data according to the crack direction information, the fast transverse wave and the corrected slow transverse wave.

9. A computer device, characterized in that the computer device comprises:

a processor and a memory, the memory having stored therein at least one program code, the at least one program code being loaded and executed by the processor to perform operations performed in the shear-wave-only shear wave splitting correction method of any one of claims 1 to 7.

10. A computer-readable storage medium having at least one program code stored therein, the at least one program code being loaded into and executed by a processor to implement the operations performed in the shear wave-splitting correction method for shear waves only according to any one of claims 1 to 7.