Methods, apparatus, equipment and storage media for correcting shear wave splitting of pure shear waves

Abstract

This disclosure provides a method, apparatus, device, and storage medium for shear wave splitting correction of pure shear waves, belonging to the field of seismic exploration technology. The method includes: acquiring fracture direction information of the strata in the area to be measured; and acquiring first pure shear wave data received by a target geophone in the area; when the fracture direction information is used to indicate a fixed fracture direction, determining the fast and slow shear waves corresponding to the first pure shear wave data based on the fracture direction information; determining the time difference between the fast and slow shear waves; correcting the slow shear wave based on the time difference to obtain a corrected slow shear wave; and determining the first corrected shear wave field corresponding to the first pure shear wave data based on the fracture direction information, the fast shear wave, and the corrected slow shear wave. Because the slow shear wave is corrected using the time difference, a shear wave field after removing azimuth anisotropy is obtained, thereby eliminating the influence of anisotropy and improving the accuracy of the determined first corrected shear wave field.

Description

Technical Field

[0001] This application relates to the field of seismic exploration technology, and in particular to a method, apparatus, equipment and storage medium for correcting shear wave splitting of pure shear waves. Background Technology

[0002] Currently, pure shear wave splitting technology is increasingly widely used in seismic exploration. For example, it is used to achieve structural imaging, oil and gas prediction, and fracture detection in the target seismic area. However, when using pure shear wave splitting to predict oil and gas in the target seismic area, if the formation medium is anisotropic, splitting correction is sometimes necessary to determine the shear wave field corresponding to the pure shear wave. Summary of the Invention

[0003] This application provides a method, apparatus, device, and storage medium for correcting shear wave splitting in pure shear waves, which can improve the accuracy of the shear wave field obtained after shear wave splitting correction. The technical solution is as follows:

[0004] On the one hand, this application provides a method for correcting shear wave splitting of pure shear waves, the method comprising:

[0005] Acquire the fracture direction information of the strata in the work area to be tested, and acquire the first pure shear wave data received by the target geophone in the work area;

[0006] When the crack direction information is used to indicate that the crack direction is a fixed direction, the fast shear wave and slow shear wave corresponding to the first pure shear wave data are determined according to the crack direction information. The fast shear wave is a shear wave parallel to the crack direction, and the slow shear wave is a shear wave perpendicular to the crack direction.

[0007] Determine the time difference between the fast shear wave and the slow shear wave, and correct the slow shear wave according to the time difference to obtain a corrected slow shear wave;

[0008] Based on the crack direction information, the fast shear wave, and the corrected slow shear wave, the first corrected shear wave field corresponding to the first pure shear wave data is determined.

[0009] In one possible implementation, determining the fast and slow shear waves corresponding to the first pure shear wave data based on the crack direction information includes:

[0010] Obtain the direction of the shot-detector line between the target seismic source and the target geophone, determine the first angle difference between the crack direction and the shot-detector line direction, and the target seismic source is the seismic source in the work area that generates the first pure shear wave data;

[0011] Determine the pure shear wave of the first target and the pure shear wave of the second target;

[0012] Based on the first angle difference, the first target pure shear wave and the second target pure shear wave are rotated in two components to obtain the fast shear wave and slow shear wave corresponding to the first pure shear wave data.

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

[0014] The determination of the first target pure shear wave and the second target pure shear wave includes:

[0015] Obtain the survey line directions of multiple detectors within the work area, and determine the second angle difference between the shot-detector connection direction and the survey line direction;

[0016] The first component wave and the second component wave corresponding to the first pure shear wave are determined, as well as the third component wave and the fourth component wave corresponding to the second pure shear wave. The first component wave is the component wave of the first pure shear wave along the direction of the measuring line, the second component wave is the component wave of the first pure shear wave in the direction perpendicular to the direction of the measuring line, the third component wave is the component wave of the second pure shear wave along the direction of the measuring line, and the fourth component wave is the component wave of the second pure shear wave in the direction perpendicular to the direction of the measuring line.

[0017] Based on the second angle difference, the first component wave, the second component wave, the third component wave, and the fourth component wave, the first target pure shear wave and the second target pure shear wave are determined.

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

[0019] Determine the fast shear wave profile corresponding to the fast shear wave and the slow shear wave profile corresponding to the slow shear wave. The fast shear wave profile includes the layer identifiers and layer times of multiple first layer waves, and the slow shear wave profile includes the layer identifiers and layer times of multiple second layer waves.

[0020] Based on the stratigraphic identifiers and stratigraphic times corresponding to the plurality of first stratigraphic waves and the stratigraphic identifiers and stratigraphic times corresponding to the plurality of second stratigraphic waves, the time difference between each first stratigraphic wave and the second stratigraphic wave corresponding to the first stratigraphic wave is determined, thereby obtaining the time difference between the fast shear wave and the slow shear wave.

[0021] Wherein, the depth of the second stratum wave corresponding to the first stratum wave is the same as the depth of the first stratum wave.

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

[0023] Obtain the measurement line directions of multiple geophones within the work area, and determine the third angle difference between the crack direction and the measurement line direction;

[0024] The first component wave and the second component wave corresponding to the first pure shear wave are determined, as well as the third component wave and the fourth component wave corresponding to the second pure shear wave. The first component wave is the component wave of the first pure shear wave along the direction of the measuring line, the second component wave is the component wave of the first pure shear wave in the direction perpendicular to the direction of the measuring line, the third component wave is the component wave of the second pure shear wave along the direction of the measuring line, and the fourth component wave is the component wave of the second pure shear wave in the direction perpendicular to the direction of the measuring line.

[0025] Based on the third angle difference, the first component wave, the second component wave, the third component wave, and the fourth component wave, determine the fast shear wave profile corresponding to the fast shear wave and the slow shear wave profile corresponding to the slow shear wave.

[0026] In one possible implementation, the first corrected shear wave field includes the SH shear wave field and the SV shear wave field; determining the first corrected shear wave field corresponding to the first pure shear wave data based on the crack direction information, the fast shear wave, and the corrected slow shear wave includes:

[0027] For any target seismic source within the work area, the direction of the shot-detector connection between the target seismic source and the target geophone is obtained, and the first angle difference between the crack direction and the shot-detector connection direction is determined. The target seismic source is the seismic source within the work area that generates the first pure shear wave data.

[0028] Based on the first angle difference, the fast shear wave and the corrected slow shear wave are reverse-rotated to obtain the SH shear wave field corresponding to the first pure shear wave data and the SV shear wave field corresponding to the first pure shear wave data.

[0029] In one possible implementation, the method further includes:

[0030] When the crack direction information is used to represent the change of crack direction with depth, multiple crack time windows corresponding to the crack direction information are determined.

[0031] Acquire the second pure shear wave data corresponding to the multiple crack time windows, as well as the target time difference and target crack direction corresponding to each crack time window;

[0032] Based on the target time difference and target crack direction corresponding to each crack time window, the second pure shear wave data is subjected to shear wave splitting correction to obtain the second corrected shear wave field corresponding to the second pure shear wave data.

[0033] On the other hand, this application provides a device for determining the splitting wave information of a transverse wave, the device comprising:

[0034] The acquisition module is used to acquire information on the direction of fractures in the strata of the work area to be tested, and to acquire the first pure shear wave data received by the target geophone in the work area.

[0035] The first determining module is used to determine the fast shear wave and slow shear wave corresponding to the first pure shear wave data according to the crack direction information when the crack direction information is used to indicate that the crack direction is a fixed direction. The fast shear wave is a shear wave parallel to the crack direction, and the slow shear wave is a shear wave perpendicular to the crack direction.

[0036] The second determining module is used to determine the time difference between the fast shear wave and the slow shear wave, and to correct the slow shear wave according to the time difference to obtain a corrected slow shear wave.

[0037] The third determining module is used to determine the first corrected shear wave field corresponding to the first pure shear wave data based on the crack direction information, the fast shear wave, and the corrected slow shear wave.

[0038] In one possible implementation, the first determining module includes:

[0039] The acquisition unit is used to acquire the direction of the shot-detector connection between the target seismic source and the target geophone, and to determine the first angle difference between the crack direction and the shot-detector connection direction. The target seismic source is the seismic source in the work area that generates the first pure shear wave data.

[0040] The first determining unit is used to determine the pure shear wave of the first target and the pure shear wave of the second target.

[0041] The rotation unit is used to perform two-component rotation on the first target pure shear wave and the second target pure shear wave according to the first angle difference, so as to obtain the fast shear wave and slow shear wave obtained by splitting the first pure shear wave data along the crack.

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

[0043] The first determining unit is used to acquire the survey line directions of multiple detectors within the work area, determine the second angle difference between the shot-detector connection direction and the survey line direction; determine the first component wave and the second component wave corresponding to the first pure shear wave, and the third component wave and the fourth component wave corresponding to the second pure shear wave, wherein the first component wave is the component wave of the first pure shear wave along the survey line direction, the second component wave is the component wave of the first pure shear wave along the direction perpendicular to the survey line direction, the third component wave is the component wave of the second pure shear wave along the survey line direction, and the fourth component wave is the component wave of the second pure shear wave along the direction perpendicular to the survey line direction; and determine the first target pure shear wave and the second target pure shear wave based on the second angle difference, the first component wave, the second component wave, the third component wave, and the fourth component wave.

[0044] In another possible implementation, the second determining module includes:

[0045] The second determining unit is used to determine the fast shear wave profile corresponding to the fast shear wave and the slow shear wave profile corresponding to the slow shear wave. The fast shear wave profile includes the layer identifiers and layer times of a plurality of first layer waves, and the slow shear wave profile includes the layer identifiers and layer times of a plurality of second layer waves.

[0046] The third determining unit is used to determine the time difference between each first-layer wave and the second-layer wave corresponding to the first-layer wave based on the layer identifier and layer time corresponding to the plurality of first-layer waves and the layer identifier and layer time corresponding to the plurality of second-layer waves, thereby obtaining the time difference between the fast shear wave and the slow shear wave.

[0047] Wherein, the depth of the second stratum wave corresponding to the first stratum wave is the same as the depth of the first stratum wave.

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

[0049] The second determining unit is used to acquire the measurement line directions of multiple detectors within the work area, determine the third angle difference between the crack direction and the measurement line direction; determine the first component wave and the second component wave corresponding to the first pure shear wave, and the third component wave and the fourth component wave corresponding to the second pure shear wave, wherein the first component wave is the component wave of the first pure shear wave along the measurement line direction, the second component wave is the component wave of the first pure shear wave along the direction perpendicular to the measurement line direction, the third component wave is the component wave of the second pure shear wave along the measurement line direction, and the fourth component wave is the component wave of the second pure shear wave along the direction perpendicular to the measurement line direction; and determine the fast shear wave profile corresponding to the fast shear wave and the slow shear wave profile corresponding to the slow shear wave based on the third angle difference, the first component wave, the second component wave, the third component wave, and the fourth component wave.

[0050] In one possible implementation, the first corrected shear wave field includes an SH shear wave field and an SV shear wave field; the third determining module is used to, for any target seismic source within the work area, obtain the direction of the shot-detector line between the target seismic source and the target geophone, determine the first angle difference between the crack direction and the shot-detector line direction, wherein the target seismic source is the seismic source within the work area that generates the first pure shear wave data; and based on the first angle difference, perform reverse rotation on the fast shear wave and the corrected slow shear wave to obtain the SH shear wave field corresponding to the first pure shear wave data and the SV shear wave field corresponding to the first pure shear wave data.

[0051] In one possible implementation, the device further includes:

[0052] The fourth determining module is used to determine multiple crack time windows corresponding to the crack direction information when the crack direction information is used to represent the change of crack direction with depth; to acquire second pure shear wave data corresponding to the multiple crack time windows and target time difference and target crack direction corresponding to each crack time window; and to perform shear wave splitting correction on the second pure shear wave data according to the target time difference and target crack direction corresponding to each crack time window to obtain the second corrected shear wave field corresponding to the second pure shear wave data.

[0053] On the other hand, embodiments of this application provide a computer device, the computer device including: a processor and a memory, the memory storing at least one piece of program code, the at least one piece of program code being loaded and executed by the processor to implement the operations performed in the pure shear wave splitting correction method described in any of the above possible implementations.

[0054] On the other hand, embodiments of this application provide a computer-readable storage medium storing at least one piece of program code, which is loaded and executed by a processor to implement the operations performed in the pure shear wave splitting correction method described in any of the above possible implementations.

[0055] The beneficial effects of the technical solutions provided in this application include at least the following:

[0056] This application provides a method for correcting shear wave splitting of pure shear waves. For anisotropic media, the method corrects the slow shear wave by using the time difference between the fast and slow shear waves, thus eliminating the influence of azimuth anisotropy. Furthermore, by reversing the fast shear wave and the corrected slow shear wave, the accuracy of the obtained shear wave field is high. Attached Figure Description

[0057] To more clearly illustrate the technical solutions in the embodiments of this application, the accompanying drawings used in the description of the embodiments will be briefly introduced below. Obviously, the accompanying drawings described below are only some embodiments of this application. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0058] Figure 1 This is a flowchart illustrating a shear wave splitting correction method for a pure shear wave according to an exemplary embodiment;

[0059] Figure 2 This is a schematic diagram illustrating four-component data of a pure shear wave according to an exemplary embodiment;

[0060] Figure 3 This is a schematic diagram illustrating a rotated four-component shear wave data according to an exemplary embodiment;

[0061] Figure 4 This is a schematic diagram illustrating a fast shear wave profile and a slow shear wave profile according to an exemplary embodiment;

[0062] Figure 5 This is a schematic diagram illustrating an SH shear wave field and an SV shear wave field according to an exemplary embodiment;

[0063] Figure 6 This is a block diagram illustrating a shear wave splitting correction device for pure shear waves according to an exemplary embodiment;

[0064] Figure 7 This is a block diagram illustrating a shear wave splitting correction device for pure shear waves according to an exemplary embodiment;

[0065] Figure 8This is a structural block diagram of a computer device according to an exemplary embodiment. Detailed Implementation

[0066] To make the objectives, technical solutions, and advantages of this application clearer, the embodiments of this application will be described in further detail below with reference to the accompanying drawings.

[0067] Figure 1 This is a flowchart illustrating a method for correcting the splitting of a pure shear wave according to an exemplary embodiment.

[0068] See Figure 1 The method includes:

[0069] 101. The computer equipment acquires information on the fracture direction of the strata in the work area to be tested, and acquires the first pure shear wave data received by the target geophone in the work area.

[0070] In this embodiment of the application, the fracture orientation information of the formation has been obtained. Optionally, this fracture orientation information is determined using pure shear wave splitting analysis methods such as Alford rotation.

[0071] In one possible implementation, the direction of the fracture can be fixed with depth, and the fractures at different depths in the formation have the same direction. In this case, the computer equipment can perform shear wave splitting correction for pure shear waves through steps 102 to 104.

[0072] In another possible implementation, the fractures at different depths in the formation have different orientations, with the fracture orientation varying with depth. In this case, the computer equipment can perform shear wave splitting correction for pure shear waves through steps 105 to 106.

[0073] It should be noted that multiple seismic sources and multiple geophones are deployed within the area to be measured. The target geophone is the one that receives the first pure shear wave data excited by the target seismic source. The target seismic source can be any one of multiple seismic sources. Optionally, the target seismic source is a horizontally excited, controllable source, and the seismic wave signal excited by the target seismic source is a pure shear wave dominated by shear waves, excited once along the direction parallel to the survey line and once perpendicular to the survey line. That is, the target geophone receives two pure shear waves. Correspondingly, the first pure shear wave data includes a first pure shear wave and a second pure shear wave, wherein the first pure shear wave is the pure shear wave excited by the target seismic source along the direction parallel to the survey line, and the second pure shear wave is the pure shear wave excited by the target seismic source along the direction perpendicular to the survey line, as received by the target geophone.

[0074] In one possible implementation, the target detector is a three-component detector. The three-component detector has two components in the horizontal direction, namely the x-component and the y-component, which are perpendicular to each other. It also has one component in the vertical direction, the z-component, which is perpendicular to both the x-component and the y-component.

[0075] Among them, the two horizontal components of the three-component detector are parallel to the measurement line direction and perpendicular to the measurement line direction, respectively, so as to receive the component wave of the first pure shear wave parallel to the measurement line direction and the component wave of the second pure shear wave parallel to the measurement line direction and the component wave of the second pure shear wave perpendicular to the measurement line direction, that is, the initial four-component data corresponding to the first pure shear wave data.

[0076] 102. When the crack direction information is used to indicate that the crack direction is fixed with depth, the computer device determines the fast shear wave and slow shear wave corresponding to the first pure shear wave data based on the crack direction information. The fast shear wave is a shear wave parallel to the crack direction, and the slow shear wave is a shear wave perpendicular to the crack direction.

[0077] In one possible implementation, this step may include the following steps (1) to (3):

[0078] (1) The computer equipment obtains the direction of the shot-detector line between the target seismic source and the target detector, and determines the first angle difference between the crack direction and the shot-detector line direction.

[0079] The target seismic source is the source within the work area that generates the first pure shear wave data; the target geophone is the geophone that receives the first pure shear wave data generated by the target seismic source. The crack direction can be represented by the letter θ. fracture The direction of the shot-receiver connection can be indicated by the letter θ. azimuth The first angular difference between the crack direction and the line connecting the shot and receiver is: θ fracture -θ azimuth Optionally, the first angular difference can be represented by the letter γ; that is, γ = θ fracture -θ azimuth .

[0080] In this step, the direction of the shot-detector connection is the direction from the source position to the detector position. Optionally, the direction of the shot-detector connection can be called the radial direction, generally represented by the letter r; the perpendicular direction of the shot-detector connection can be called the transverse direction, generally represented by the letter t.

[0081] (2) The computer equipment determines the first target pure transverse wave and the second target pure transverse wave.

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

[0083] (21) The computer equipment acquires the measurement line direction of multiple detectors in the work area and determines the second angle difference between the shot-detector connection direction and the measurement line direction.

[0084] The direction of the survey line can be represented by the letter θ. inline This indicates that the second angular difference between the direction of the shot-receiver line and the direction of the survey line is: θ azimuth -θ inline Optionally, the second angular difference can be represented by the letter β; that is, β = θ. azimuth -θ inline .

[0085] (22) The computer device determines the first component wave and the second component wave corresponding to the first pure transverse wave, and the third component wave and the fourth component wave corresponding to the second pure transverse wave.

[0086] It should be noted that the target detector is a three-component detector. A three-component detector has two components in the horizontal direction, namely the x-component and the y-component, which are perpendicular to each other. In the vertical direction, it has one component, the z-component, which is perpendicular to both the x-component and the y-component.

[0087] In one possible implementation, the x-component is parallel to the survey line direction, and the y-component is perpendicular to the survey line direction. The target source generates one pure shear wave each along the survey line direction and in the direction perpendicular to the survey line direction. Specifically, the first component wave is a component of the first pure shear wave received by the target detector along the survey line direction; the second component wave is a component of the first pure shear wave received by the target detector in the direction perpendicular to the survey line direction; the third component wave is a component of the second pure shear wave received by the target detector along the survey line direction; and the fourth component wave is a component of the second pure shear wave received by the target detector in the direction perpendicular to the survey line direction.

[0088] In one possible implementation, the first component wave can be represented by S. x R x This means that the second component wave can be represented by S. x R y The third component wave can be represented by S. y R x This means that the fourth component wave can be represented by S. y R y Indicated. Optional, see [link to other documentation]. Figure 2 The initial four-component data corresponding to the first pure shear wave data include the first component wave, the second component wave, the third component wave, and the fourth component wave.

[0089] (23) The computer equipment determines the pure transverse wave of the first target and the pure transverse wave of the second target based on the second angle difference, the first component wave, the second component wave, the third component wave and the fourth component wave.

[0090] In one possible implementation, the first target pure shear wave includes a first radial pure shear wave and a first tangential pure shear wave; the second target pure shear wave includes a second radial pure shear wave and a second tangential pure shear wave; the computer equipment performs a four-component rotation on the initial four-component data of the pure shear wave based on the second angular difference between the shot-receiver connection direction and the survey line direction to obtain the rotated four-component shear wave data. Accordingly, this step is as follows: the computer equipment determines the first radial pure shear wave, the first tangential pure shear wave, the second radial pure shear wave, and the second tangential pure shear wave based on the second angular difference, the first component wave, the second component wave, the third component wave, and the fourth component wave using the following formula:

[0091] Formula 1:

[0092] Among them, S x R x S represents the first component wave. x R y S represents the second component wave. y R x S represents the third component wave. y R y This represents the fourth component wave, β represents the second angular difference, and S... t R r S represents the first radial pure transverse wave. t R t S represents the first tangential pure transverse wave. t R r S represents the second radial pure transverse wave. t R t This indicates the second tangential pure transverse wave.

[0093] Optional, see Figure 3 The rotated four-component shear wave data includes the 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 Second tangential pure transverse wave S r R t .

[0094] (3) 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, and obtains the fast shear wave and slow shear wave corresponding to the first pure shear wave data.

[0095] The first target pure shear wave includes a first radial pure shear wave and a first tangential pure shear wave; the second target pure shear wave includes a second radial pure shear wave and a second tangential pure shear wave. Accordingly, this step includes: the computer device performing a 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 device performing a two-component rotation on the second radial pure shear wave and the second tangential pure shear wave to obtain a second fast shear wave and a second slow shear wave corresponding to the second target pure shear wave.

[0096] In one possible implementation, the first radial pure shear wave and the first tangential pure shear wave are rotated in two components by rotating the receiving direction of the detector to be parallel to and perpendicular to the crack direction. Correspondingly, the computer device performs the two-component rotation of the first radial pure shear wave and the first tangential pure shear wave as follows: based on the first angular difference, the computer device performs the two-component rotation of the first radial pure shear wave and the first tangential pure shear wave using the following formula two, to obtain the first fast shear wave and the first slow shear wave corresponding to the first target pure shear wave.

[0097] Formula 2:

[0098] Among them, S t R s1 S represents the first transverse wave. t R s2 S represents the first slow transverse wave. t R r S represents the first radial pure transverse wave. t R t γ represents the first tangential pure transverse wave, and γ represents the first angular difference.

[0099] The steps for the computer equipment to perform two-component rotation on the second radial pure shear wave and the second tangential pure shear wave are as follows: The computer equipment performs two-component rotation on the second radial pure shear wave and the second tangential pure shear wave according to the first angle difference and the following formula three to obtain the second fast shear wave and the second slow shear wave corresponding to the second target pure shear wave.

[0100] Formula 3:

[0101] Among them, S r R s1 S represents the second fastest transverse wave. r R s2 Indicates the second slow transverse wave, S r R r S represents the second radial pure transverse wave. r R t This represents the second tangential pure transverse wave, and γ represents the first angular difference.

[0102] 103. The computer equipment determines the time difference between the fast and slow shear waves, and corrects the slow shear wave based on the time difference to obtain the corrected slow shear wave.

[0103] In one possible implementation, the computer device determines the time difference between the fast and slow shear waves through the following steps (1) to (2):

[0104] (1) The computer equipment determines the fast shear wave profile corresponding to the fast shear wave and the slow shear wave profile corresponding to the slow shear wave. The fast shear wave profile includes the layer identifier and layer time of multiple first layer waves corresponding to the fast shear wave within a preset time period. The slow shear wave profile includes the layer identifier and layer time of multiple second layer waves corresponding to the slow shear wave within a preset time period.

[0105] In one possible implementation, the first pure shear wave data includes a first pure shear wave and a second pure shear wave. Accordingly, this step is as follows: the computer equipment acquires the survey line direction of the work area and determines the third angle difference between the crack direction and the survey line direction; it determines the first and second component waves corresponding to the first pure shear wave, as well as the third and fourth component waves corresponding to the second pure shear wave, wherein the first component wave is the component wave of the first pure shear wave along the survey line direction, the second component wave is the component wave of the first pure shear wave perpendicular to the survey line direction, the third component wave is the component wave of the second pure shear wave along the survey line direction, and the fourth component wave is the component wave of the second pure shear wave perpendicular to the survey line direction; based on the third angle difference, the first component wave, the second component wave, the third component wave, and the fourth component wave, it determines the fast shear wave profile corresponding to the fast shear wave and the slow shear wave profile corresponding to the slow shear wave.

[0106] In one possible implementation, the first component wave, the second component wave, the third component wave, and the fourth component wave are the initial four-component data corresponding to the first pure shear wave data. The computer device separates the fast and slow waves of the four-component data using a rotation formula. Accordingly, the steps for the computer device to determine the fast shear wave profile corresponding to the fast shear wave and the slow shear wave profile corresponding to the slow shear wave based on the third angle difference, the first component wave, the second component wave, the third component wave, and the fourth component wave are as follows: The computer device determines the fast shear wave profile corresponding to the fast shear wave and the slow shear wave profile corresponding to the slow shear wave using the following formula four based on the third angle difference, the first component wave, the second component wave, the third component wave, and the fourth component wave.

[0107] Formula 4:

[0108] Among them, S x R x S represents the first component wave. x R y S represents the second component wave. y R x S represents the third component wave.y R y Indicates the fourth component wave, S represents the difference in the third angle. S1 R S1 S represents the fast shear wave profile. S2 R S2 This represents a slow transverse wave profile.

[0109] The direction of the crack can be represented by the letter θ. fracture The direction of the survey line can be indicated by the letter θ. inline This indicates that the third angle difference between the crack direction and the survey line direction is: θ fracture -θ inline Optionally, the third angle difference can be represented by letters. In other words,

[0110] It should be noted that the fast shear wave profile corresponds to the profile information of a fast shear wave, while the slow shear wave profile corresponds to the profile information of a slow shear wave. Optional, see [link to relevant documentation]. Figure 4 The computer equipment uses a rotation formula to separate the fast and slow waves from the initial four-component data of the shear wave, obtaining the profile information of the fast shear wave and the profile information of the slow shear wave, as shown below. Figure 4 As shown, the fast shear wave profile is Figure 4 The leftmost shear wave profile is the slow shear wave profile. Figure 4 The rightmost transverse wave profile.

[0111] (2) The computer equipment determines the time difference between each first-layer wave and the corresponding second-layer wave based on the layer identifier and layer time of multiple first-layer waves in the fast shear wave profile and the layer identifier and layer time of multiple second-layer waves in the slow shear wave profile, thus obtaining the time difference between the fast shear wave and the slow shear wave.

[0112] In this section, the layer markers of multiple first-layer waves in the fast shear wave profile correspond one-to-one with the layer markers of multiple second-layer waves in the slow shear wave profile, from shallowest to deepest. For example, the layer markers of multiple first-layer waves in the fast shear wave profile are h1s1, h2s1, and h3s1, and the layer markers of multiple second-layer waves in the slow shear wave profile are h1s2, h2s2, and h3s2; where h1s1, h2s1, and h3s1 correspond one-to-one with h1s2, h2s2, and h3s2 from shallowest to deepest, that is, h1s1 corresponds to h1s2, h2s1 corresponds to h2s2, and h3s1 corresponds to h3s2.

[0113] In one possible implementation, this step is as follows: The computer device determines the time difference between the first first-layer wave and the first second-layer wave in a sequence from shallow to deep, based on the layer time of the first first-layer wave in the fast shear wave profile and the layer time of the first second-layer wave in the slow shear wave profile. Based on the layer time of the second first-layer wave in the fast shear wave profile and the layer time of the second second-layer wave in the slow shear wave profile, the time difference between the second first-layer wave and the second second-layer wave is determined sequentially from shallow to deep, thus obtaining the time difference between the fast shear wave and the slow shear wave.

[0114] In step 103, the computer device corrects the slow shear wave according to the time difference to obtain the corrected slow shear wave. The steps are as follows: the computer device corrects the first slow shear wave corresponding to the pure shear wave of the first target according to the time difference using the following formula five to obtain the first corrected slow shear wave; and, according to the time difference, corrects the second slow shear wave corresponding to the pure shear wave of the second target using the following formula six to obtain the second corrected slow shear wave.

[0115] Formula 5: S t R s2up (t)=S t R s2 (t+Δt)

[0116] Formula Six: S r R s2up (t)=S r R s2 (t+Δt)

[0117] Among them, S t R s2up (t) represents the first corrected slow transverse wave at time t, S r R s2up (t) represents the second corrected slow transverse wave at time t, Δt represents the time difference, and S t R s2 (t+Δt) represents the first slow transverse wave corresponding to time t+Δt, S r R s2 (t+Δt) represents the second slow transverse wave corresponding to time t+Δt. Optionally, t can represent a point in time or a time interval.

[0118] It should be noted that this time difference can vary with depth. In other words, the fast and slow wave time differences can be different for different layers. In this case, for a pure shear wave between two layers, the fast and slow wave time difference is the interpolation result of the time difference between the upper layer and the lower layer.

[0119] Another point to note is that the slow shear wave profile represents the amplitude value corresponding to the slow shear wave, where the amplitude value is obtained based on the first pure shear wave data periodically acquired by the target detector. In one possible implementation, during the time difference correction process for the slow wave, when there is no corresponding amplitude value for t+Δt, the computer device determines the amplitude value corresponding to t+Δt based on the sinc function interpolation.

[0120] In this embodiment of the application, when there is no corresponding amplitude value for the time after the slow shear wave time difference correction, the computer device can determine the amplitude value corresponding to the time by interpolation using the sinc function, thereby improving the accuracy of time difference correction.

[0121] 104. The computer equipment determines the first corrected shear wave field corresponding to the first pure shear wave data based on the crack direction information, the fast shear wave, and the corrected slow shear wave.

[0122] In one possible implementation, the first corrected shear wave field includes the SH shear wave field and the SV shear wave field. Accordingly, this step is as follows: For any target seismic source within the work area, the computer equipment acquires the direction of the shot-detector line between the target seismic source and the target geophone, determines the first angular difference between the crack direction and the shot-detector line direction, and the target seismic source is the seismic source within the work area that generates the first pure shear wave data; based on the first angular difference, the fast shear wave and the corrected slow shear wave are reverse-rotated to obtain the SH shear wave field corresponding to the first pure shear wave data and the SV shear wave field corresponding to the first pure shear wave data.

[0123] In one possible implementation, the fast shear wave includes a first fast shear wave corresponding to the first target pure shear wave and a second fast shear wave corresponding to the second target pure shear wave; the corrected slow shear wave includes a first corrected slow shear wave corresponding to the first target pure shear wave and a second corrected slow shear wave corresponding to the second target pure shear wave. Accordingly, the computer device performs an inverse rotation on the fast shear wave and the corrected slow shear wave based on the first angle difference to obtain the SH shear wave field and the SV shear wave field corresponding to the first pure shear wave data. The steps are as follows: the computer device performs an inverse rotation on the first fast shear wave and the first corrected slow shear wave corresponding to the first target pure shear wave according to the first angle difference using the following formula (7) to obtain the SH shear wave field corresponding to the first pure shear wave data; and, according to the first angle difference, performs an inverse rotation on the second fast shear wave and the second corrected slow shear wave corresponding to the second target pure shear wave using the following formula (8) to obtain the SV shear wave field corresponding to the first pure shear wave data.

[0124] Formula 7:

[0125] Formula 8:

[0126] Among them, S t R s1S represents the first transverse wave. t R s2up Indicates the first corrected slow transverse wave, S r R s1 S represents the second fastest transverse wave. r R s2up Indicates the second corrected slow transverse wave, S t R t ' represents the SH shear wave field corresponding to the first pure shear wave data, S r R r ' indicates the SV shear wave field corresponding to the first pure shear wave data.

[0127] One point that needs to be clarified is that S t R r 'and S t R t 'For the first objective, the pure shear wave, after shear wave splitting correction, the shear wave field, S r R r 'and S r R t ' represents the shear wave field after shear wave splitting correction for the second target pure shear wave. Where S... t R t 'This is the corrected SH wave field, S t R r 'The effective energy will become 0, that is, S' t R t 'To eliminate the anisotropic SH transverse wave field. Wherein, S r R r 'This is the corrected SV wave field, S r R t 'The effective energy will become 0, that is, S' r R r 'To eliminate the anisotropic SV transverse wave field.'

[0128] In one possible implementation, see Figure 5 SV wave after shear wave splitting correction Figure 5 The transverse wave field on the left, i.e., S r R r ';SH wave after shear wave splitting correction Figure 5 The transverse wave field on the right, i.e., S r R t ', Figure 5 The two transverse wave fields in the middle are S r R t 'and S t R r ' indicates that the energy has become 0.

[0129] One point to note is that, please continue to see... Figure 1When the crack direction information is used to represent the change of crack direction with depth, the computer device executes step 101 and directly executes steps 105 to 107, without executing steps 102 to 104.

[0130] 105. When crack direction information is used to represent the change of crack direction with depth, the computer device determines multiple crack time windows corresponding to the crack direction information.

[0131] In one possible implementation, when the crack direction changes with depth, the crack direction corresponding to different depths is different. Accordingly, the computer device determines the multiple crack time windows corresponding to the crack direction information as follows: the computer device determines the multiple crack directions corresponding to the crack at different depths, and determines multiple crack time windows based on the multiple crack directions, wherein one crack direction corresponds to one crack time window.

[0132] 106. The computer equipment acquires the second pure shear wave data corresponding to multiple crack time windows, as well as the target time difference and target crack direction corresponding to each crack time window.

[0133] For each crack time window, the target time difference is the time difference between the fast and slow shear waves that split off when the second pure shear wave data passes through that crack time window. The target crack direction is the crack direction under that crack time window.

[0134] In this embodiment, the target time difference and crack direction can be obtained by the operator through shear wave splitting analysis. Accordingly, this step involves the operator uploading the target time difference and crack direction to a computer device, which then acquires the second pure shear wave data corresponding to multiple crack time windows, as well as the target time difference and target crack direction for each crack time window.

[0135] 107. The computer equipment performs shear wave splitting correction on the second pure shear wave data according to the target time difference and target crack direction corresponding to each crack time window, and obtains the second corrected shear wave field corresponding to the second pure shear wave data.

[0136] 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 includes an SV wave field parallel to the shot-receiver connection direction and an SH wave field perpendicular to the shot-receiver connection direction. Accordingly, this step is as follows: the computer equipment, based on the target time difference and target crack direction corresponding to each crack time window, performs shear wave splitting correction on the third pure shear wave using the following formula (nine) to obtain the SV wave field of the second pure shear wave data parallel to the shot-receiver connection direction; and, based on the target time difference and target crack direction corresponding to each crack time window, performs shear wave splitting correction on the fourth pure shear wave using the following formula (ten) to obtain the SH wave field of the second pure shear wave data perpendicular to the shot-receiver connection direction.

[0137] Formula Nine:

[0138] U 0sr (ω)=R1 -1 D1…R n-1 -1 D n-1 R n -1 D n R n D n-1 R n-1 …D1R1U sr (ω)

[0139] Formula 10:

[0140] U 0st (ω)=R1 -1 D1…R n-1 -1 D n-1 R n -1 D n R n D n-1 R n-1 …D1R1U st (ω)

[0141]

[0142]

[0143]

[0144]

[0145] Where n represents the number of windows when there are multiple cracks, U sr (ω) represents the third pure transverse wave, U st (ω) represents the fourth pure transverse wave, U 0sr (ω) represents the SV shear wave field after shear wave splitting correction, U 0st (ω) represents the SH shear wave field after shear wave splitting correction, Δt n θ represents the target time difference corresponding to the nth crack window. n The direction of the target crack corresponds to the nth crack window, and α represents the direction of the shot-receiver connection.

[0146] Optional, U 0sr (ω) represents the SV transverse wave field after anisotropy has been eliminated. U 0st (ω) represents the SH transverse wave field with anisotropy eliminated.

[0147] This application provides a method for correcting shear wave splitting of pure shear waves. For anisotropic media, the method corrects the slow shear wave by using the time difference between the fast and slow shear waves, thus eliminating the influence of azimuth anisotropy. Furthermore, by reversing the fast shear wave and the corrected slow shear wave, the accuracy of the obtained shear wave field is high.

[0148] Figure 6 This is a block diagram illustrating a shear wave splitting correction device for pure shear waves according to an exemplary embodiment. See also Figure 6 The device includes:

[0149] The acquisition module 601 is used to acquire the fracture direction information of the strata in the work area to be tested, and to acquire the first pure shear wave data received by the target geophone in the work area.

[0150] The first determining module 602 is used to determine the fast shear wave and slow shear wave corresponding to the first pure shear wave data according to the crack direction information when the crack direction information is used to indicate that the crack direction is a fixed direction. The fast shear wave is a shear wave parallel to the crack direction, and the slow shear wave is a shear wave perpendicular to the crack direction.

[0151] The second determining module 603 is used to determine the time difference between the fast shear wave and the slow shear wave, and to correct the slow shear wave based on the time difference to obtain the corrected slow shear wave.

[0152] The third determining module 604 is used to determine the first corrected shear wave field corresponding to the first pure shear wave data based on the crack direction information, the fast shear wave, and the corrected slow shear wave.

[0153] In one possible implementation, see Figure 7 The first determining module 602 includes:

[0154] The acquisition unit 6021 is used to acquire the direction of the shot-detector connection between the target source and the target detector, and to determine the first angle difference between the crack direction and the shot-detector connection direction. The target source is the source that generates the first pure shear wave data in the work area.

[0155] The first determining unit 6022 is used to determine the pure shear wave of the first target and the pure shear wave of the second target.

[0156] The rotation unit 6023 is used to perform two-component rotation on the first target pure shear wave and the second target pure shear wave according to the first angle difference, so as to obtain the fast shear wave and slow shear wave corresponding to the first pure shear wave data.

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

[0158] The first determining unit 6021 is used to acquire the survey line directions of multiple detectors within the work area, determine the second angle difference between the shot-detector connection direction and the survey line direction; determine the first component wave and the second component wave corresponding to the first pure shear wave, as well as the third component wave and the fourth component wave corresponding to the second pure shear wave, wherein the first component wave is the component wave of the first pure shear wave along the survey line direction, the second component wave is the component wave of the first pure shear wave perpendicular to the survey line direction, the third component wave is the component wave of the second pure shear wave along the survey line direction, and the fourth component wave is the component wave of the second pure shear wave perpendicular to the survey line direction; and determine the first target pure shear wave and the second target pure shear wave based on the second angle difference, the first component wave, the second component wave, the third component wave, and the fourth component wave.

[0159] In another possible implementation, see [link to previous section]. Figure 7 The second determining module 603 includes:

[0160] The second determining unit 6031 is used to determine the fast shear wave profile corresponding to the fast shear wave and the slow shear wave profile corresponding to the slow shear wave. The fast shear wave profile includes the layer identifiers and layer times of multiple first layer waves, and the slow shear wave profile includes the layer identifiers and layer times of multiple second layer waves.

[0161] The third determining unit 6032 is used to determine the time difference between each first-layer wave and the second-layer wave corresponding to the first-layer wave based on the layer identifier and layer time corresponding to multiple first-layer waves and the layer identifier and layer time corresponding to multiple second-layer waves, so as to obtain the time difference between fast shear waves and slow shear waves.

[0162] Among them, the depth of the second stratum wave corresponding to the first stratum wave is the same as the depth of the first stratum wave.

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

[0164] The second determining unit 6031 is used to acquire the measurement line directions of multiple detectors in the work area, determine the third angle difference between the crack direction and the measurement line direction; determine the first component wave and the second component wave corresponding to the first pure shear wave, as well as the third component wave and the fourth component wave corresponding to the second pure shear wave, wherein the first component wave is the component wave of the first pure shear wave along the measurement line direction, the second component wave is the component wave of the first pure shear wave perpendicular to the measurement line direction, the third component wave is the component wave of the second pure shear wave along the measurement line direction, and the fourth component wave is the component wave of the second pure shear wave perpendicular to the measurement line direction; and determine the fast shear wave profile corresponding to the fast shear wave and the slow shear wave profile corresponding to the slow shear wave based on the third angle difference, the first component wave, the second component wave, the third component wave and the fourth component wave.

[0165] In one possible implementation, the first corrected shear wave field includes the SH shear wave field and the SV shear wave field; the third determining module 604 is used to obtain the direction of the shot-detector connection between the target source and the target detector for any target source in the work area, and determine the first angle difference between the crack direction and the shot-detector connection direction, wherein the target source is the source in the work area that generates the first pure shear wave data; based on the first angle difference, the fast shear wave and the corrected slow shear wave are reverse-rotated to obtain the SH shear wave field corresponding to the first pure shear wave data and the SV shear wave field corresponding to the first pure shear wave data.

[0166] In one possible implementation, see [link to previous section] Figure 7 The device also includes:

[0167] The fourth determining module 605 is used to determine multiple crack time windows corresponding to the crack direction information when the crack direction information is used to represent the change of crack direction with depth; to acquire the second pure shear wave data corresponding to the multiple crack time windows and the target time difference and target crack direction corresponding to each crack time window; and to perform shear wave splitting correction on the second pure shear wave data according to the target time difference and target crack direction corresponding to each crack time window to obtain the second corrected shear wave field corresponding to the second pure shear wave data.

[0168] This application provides a shear wave splitting correction device for pure shear waves. For anisotropic media formations, the slow shear wave is corrected by the time difference between the fast and slow shear waves, thus eliminating the influence of azimuth anisotropy. Furthermore, by reversing the fast and slow shear waves, the accuracy of the obtained shear wave field is high.

[0169] Figure 8 The diagram illustrates a structural block diagram of a terminal 800 provided in an exemplary embodiment of the present invention. The terminal 800 may be a smartphone, tablet computer, MP3 player (Moving Picture Experts Group Audio Layer III), MP4 player (Moving Picture Experts Group Audio Layer IV), laptop computer, or desktop computer. The terminal 800 may also be referred to as a user device, portable terminal, laptop terminal, desktop terminal, or other names.

[0170] Typically, terminal 800 includes a processor 801 and a memory 802.

[0171] Processor 801 may include one or more processing cores, such as a quad-core processor or an octa-core processor. Processor 801 may be implemented using at least one hardware form selected from DSP (Digital Signal Processing), FPGA (Field-Programmable Gate Array), and PLA (Programmable Logic Array). Processor 801 may also include a main processor and a coprocessor. The main processor, also known as a CPU (Central Processing Unit), is used to process data in the wake-up state; the coprocessor is a low-power processor used to process data in the standby state. In some embodiments, processor 801 may integrate a GPU (Graphics Processing Unit), which is responsible for rendering and drawing the content to be displayed on the screen. In some embodiments, processor 801 may also include an AI (Artificial Intelligence) processor, which is used to handle computational operations related to machine learning.

[0172] The memory 802 may include one or more computer-readable storage media, which may be non-transitory. The memory 802 may also include high-speed random access memory and non-volatile memory, such as one or more disk storage devices or flash memory devices. In some embodiments, the non-transitory computer-readable storage media in the memory 802 is used to store at least one instruction, which is executed by the processor 801 to implement the shear wave splitting correction method for pure shear waves provided in the method embodiments of this application.

[0173] In some embodiments, the terminal 800 may also optionally include a peripheral device interface 803 and at least one peripheral device. The processor 801, memory 802, and peripheral device interface 803 can be connected via a bus or signal line. Each peripheral device can be connected to the peripheral device interface 803 via a bus, signal line, or circuit board. Specifically, the peripheral device includes at least one of the following: a radio frequency circuit 804, a display screen 805, a camera 806, an audio circuit 807, a positioning component 808, and a power supply 809.

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

[0175] The radio frequency (RF) circuit 804 is used to receive and transmit RF (Radio Frequency) signals, also known as electromagnetic signals. The RF circuit 804 communicates with communication networks and other communication devices via electromagnetic signals. The RF circuit 804 converts electrical signals into electromagnetic signals for transmission, or converts received electromagnetic signals back into electrical signals. Optionally, the RF 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 user identity module card, etc. The RF circuit 804 can communicate with other terminals through at least one wireless communication protocol. This wireless communication protocol includes, but is not limited to: metropolitan area networks (MANs), various generations of mobile communication networks (2G, 3G, 4G, and 5G), wireless local area networks (WLANs), and / or WiFi (Wireless Fidelity) networks. In some embodiments, the RF circuit 804 may also include circuitry related to NFC (Near Field Communication), which is not limited in this application.

[0176] Display screen 805 is used to display a UI (User Interface). This UI may include graphics, text, icons, videos, and any combination thereof. When display screen 805 is a touch display screen, it also has the ability to collect touch signals on or above its surface. These touch signals can be input as control signals to processor 801 for processing. In this case, display screen 805 can also be used to provide virtual buttons and / or a virtual keyboard, also known as soft buttons and / or a soft keyboard. In some embodiments, there may be one display screen 805, which serves as the front panel of terminal 800; in other embodiments, there may be at least two display screens, respectively disposed on different surfaces of terminal 800 or in a folded design; in still other embodiments, display screen 805 may be a flexible display screen, disposed on a curved or folded surface of terminal 800. Furthermore, display screen 805 may be configured as a non-rectangular irregular shape, i.e., a non-rectangular screen. Display screen 805 may be made of materials such as LCD (Liquid Crystal Display) or OLED (Organic Light-Emitting Diode).

[0177] The camera assembly 806 is used to acquire images or videos. Optionally, the camera assembly 806 includes a front-facing camera and a rear-facing camera. Typically, the front-facing camera is located on the front panel of the terminal, and the rear-facing camera is located on the back of the terminal. In some embodiments, there are at least two rear-facing cameras, which are any one of a main camera, a depth-sensing camera, a wide-angle camera, and a telephoto camera, to achieve background blurring by fusion of the main camera and the depth-sensing camera, panoramic shooting by fusion of the main camera and the wide-angle camera, VR (Virtual Reality) shooting, or other fusion shooting functions. In some embodiments, the camera assembly 806 may also include a flash. The flash can be a single-color temperature flash or a dual-color temperature flash. A dual-color temperature flash refers to a combination of a warm-light flash and a cool-light flash, which can be used for light compensation at different color temperatures.

[0178] The audio circuit 807 may include a microphone and a speaker. The microphone is used to collect sound waves from the user and the environment, converting the sound waves into electrical signals that are input to the processor 801 for processing, or input to the radio frequency circuit 804 to achieve voice communication. For stereo sound acquisition or noise reduction purposes, multiple microphones may be used, each located at a different part of the terminal 800. The microphone may also be an array microphone or an omnidirectional microphone. The speaker is used to convert the electrical signals from the processor 801 or the radio frequency circuit 804 into sound waves. The speaker may be a conventional diaphragm speaker or a piezoelectric ceramic speaker. When the speaker is a piezoelectric ceramic speaker, it can convert electrical signals not only into audible sound waves but also into inaudible sound waves for purposes such as distance measurement. In some embodiments, the audio circuit 807 may also include a headphone jack.

[0179] The positioning component 808 is used to determine the current geographic location of the terminal 800 in order to enable navigation or LBS (Location Based Service). The positioning component 808 can be a positioning component based on the US GPS (Global Positioning System), China's BeiDou system, Russia's Granas system, or the European Union's Galileo system.

[0180] Power supply 809 is used to supply power to the various components in terminal 800. Power supply 809 can be AC ​​power, DC power, a disposable battery, or a rechargeable battery. When power supply 809 includes a rechargeable battery, the rechargeable battery can support wired or wireless charging. The rechargeable battery can also be used to support fast charging technology.

[0181] In some embodiments, the terminal 800 further includes one or more sensors 810. The one or more sensors 810 include, but are not limited to: an accelerometer 811, a gyroscope 812, a pressure sensor 813, a fingerprint sensor 814, an optical sensor 815, and a proximity sensor 816.

[0182] Accelerometer 811 can detect the magnitude of acceleration on the three coordinate axes of a coordinate system established by terminal 800. For example, accelerometer 811 can be used to detect the components of gravitational acceleration on the three coordinate axes. Processor 801 can control display screen 805 to display the user interface in either a landscape or portrait view based on the gravitational acceleration signal acquired by accelerometer 811. Accelerometer 811 can also be used for games or for acquiring user motion data.

[0183] The gyroscope sensor 812 can detect the orientation and rotation angle of the terminal 800. The gyroscope sensor 812, in conjunction with the accelerometer sensor 811, can collect 3D motion data from the user on the terminal 800. Based on the data collected by the gyroscope sensor 812, the processor 801 can perform the following functions: motion sensing (e.g., changing the UI based on the user's tilt), image stabilization during shooting, game control, and inertial navigation.

[0184] The pressure sensor 813 can be disposed on the side bezel of the terminal 800 and / or on the lower layer of the display screen 805. When the pressure sensor 813 is disposed on the side bezel of the terminal 800, it can detect the user's grip signal on the terminal 800, and the processor 801 can perform left / right hand recognition or quick operation based on the grip signal collected by the pressure sensor 813. When the pressure sensor 813 is disposed on the lower layer of the display screen 805, the processor 801 can control the operable controls on the UI interface based on the user's pressure operation on the display screen 805. The operable controls include at least one of button controls, scroll bar controls, icon controls, and menu controls.

[0185] The fingerprint sensor 814 is used to collect the user's fingerprint. The processor 801 identifies the user's identity based on the fingerprint collected by the fingerprint sensor 814, or the fingerprint sensor 814 identifies the user's identity based on the collected fingerprint. When the user's identity is identified as trusted, the processor 801 authorizes the user to perform relevant sensitive operations, including unlocking the screen, viewing encrypted information, downloading software, making payments, and changing settings. The fingerprint sensor 814 can be located on the front, back, or side of the terminal 800. When the terminal 800 has physical buttons or a manufacturer's logo, the fingerprint sensor 814 can be integrated with the physical buttons or manufacturer's logo.

[0186] An optical sensor 815 is used to collect ambient light intensity. In one embodiment, the processor 801 can control the display brightness of the display screen 805 based on the ambient light intensity collected by the 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 screen 805 is decreased. In another embodiment, the processor 801 can also dynamically adjust the shooting parameters of the camera assembly 806 based on the ambient light intensity collected by the optical sensor 815.

[0187] The proximity sensor 816, also known as a distance sensor, is typically located on the front panel of the terminal 800. The proximity sensor 816 is used to detect the distance between the user and the front of the terminal 800. In one embodiment, when the proximity sensor 816 detects that the distance between the user and the front of the terminal 800 is gradually decreasing, the processor 801 controls the display screen 805 to switch from a screen-on state to a screen-off state; when the proximity sensor 816 detects that the distance between the user and the front of the terminal 800 is gradually increasing, the processor 801 controls the display screen 805 to switch from a screen-off state to a screen-on state.

[0188] Those skilled in the art will understand that Figure 8 The structure shown does not constitute a limitation on terminal 800 and may include more or fewer components than shown, or combine certain components, or use different component arrangements.

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

[0190] Other embodiments of this disclosure will readily occur to those skilled in the art upon consideration of the specification and practice of the invention disclosed herein. This application is intended to cover any variations, uses, or adaptations of this disclosure that follow the general principles of this disclosure and include common knowledge or customary techniques in the art not disclosed herein. The specification and examples are to be considered exemplary only, and the true scope and spirit of this disclosure are indicated by the following claims.

[0191] The above description is merely an optional embodiment of this application and is not intended to limit this application. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of this application should be included within the protection scope of this application.

Claims

1. A method for correcting shear wave splitting in pure shear waves, characterized in that, The method includes: Acquire the fracture direction information of the strata in the work area to be tested, and acquire the first pure shear wave data received by the target geophone in the work area; When the crack direction information is used to indicate that the crack direction is a fixed direction, the fast shear wave and slow shear wave corresponding to the first pure shear wave data are determined according to the crack direction information. The fast shear wave is a shear wave parallel to the crack direction, and the slow shear wave is a shear wave perpendicular to the crack direction. Determine the time difference between the fast shear wave and the slow shear wave, and correct the slow shear wave according to the time difference to obtain a corrected slow shear wave; Based on the crack direction information, the fast shear wave, and the corrected slow shear wave, the first corrected shear wave field corresponding to the first pure shear wave data is determined.

2. The method according to claim 1, characterized in that, The step of determining the fast and slow shear waves corresponding to the first pure shear wave data based on the crack direction information includes: Obtain the direction of the shot-detector line between the target seismic source and the target geophone, determine the first angle difference between the crack direction and the shot-detector line direction, and the target seismic source is the seismic source in the work area that generates the first pure shear wave data; Determine the pure shear wave of the first target and the pure shear wave of the second target; Based on the first angle difference, the first target pure shear wave and the second target pure shear wave are rotated in two components to obtain the fast shear wave and slow shear wave corresponding to the first pure shear wave data.

3. The method according to claim 2, characterized in that, The first pure shear wave data includes a first pure shear wave and a second pure shear wave; The determination of the first target pure shear wave and the second target pure shear wave includes: Obtain the survey line directions of multiple detectors within the work area, and determine the second angle difference between the shot-detector connection direction and the survey line direction; The first component wave and the second component wave corresponding to the first pure shear wave are determined, as well as the third component wave and the fourth component wave corresponding to the second pure shear wave. The first component wave is the component wave of the first pure shear wave along the direction of the measuring line, the second component wave is the component wave of the first pure shear wave in the direction perpendicular to the direction of the measuring line, the third component wave is the component wave of the second pure shear wave along the direction of the measuring line, and the fourth component wave is the component wave of the second pure shear wave in the direction perpendicular to the direction of the measuring line. Based on the second angle difference, the first component wave, the second component wave, the third component wave, and the fourth component wave, the first target pure shear wave and the second target pure shear wave are determined.

4. The method according to claim 1, characterized in that, Determining the time difference between the fast shear wave and the slow shear wave includes: Determine the fast shear wave profile corresponding to the fast shear wave and the slow shear wave profile corresponding to the slow shear wave. The fast shear wave profile includes the layer identifiers and layer times of multiple first layer waves, and the slow shear wave profile includes the layer identifiers and layer times of multiple second layer waves. Based on the stratigraphic identifiers and stratigraphic times corresponding to the plurality of first stratigraphic waves and the stratigraphic identifiers and stratigraphic times corresponding to the plurality of second stratigraphic waves, the time difference between each first stratigraphic wave and the second stratigraphic wave corresponding to the first stratigraphic wave is determined, thereby obtaining the time difference between the fast shear wave and the slow shear wave. Wherein, the depth of the second stratum wave corresponding to the first stratum wave is the same as the depth of the first stratum wave.

5. The method according to claim 4, characterized in that, The first pure shear wave data includes a first pure shear wave and a second pure shear wave; Determining the fast shear wave profile corresponding to the fast shear wave and the slow shear wave profile corresponding to the slow shear wave includes: Obtain the measurement line directions of multiple geophones within the work area, and determine the third angle difference between the crack direction and the measurement line direction; The first component wave and the second component wave corresponding to the first pure shear wave are determined, as well as the third component wave and the fourth component wave corresponding to the second pure shear wave. The first component wave is the component wave of the first pure shear wave along the direction of the measuring line, the second component wave is the component wave of the first pure shear wave in the direction perpendicular to the direction of the measuring line, the third component wave is the component wave of the second pure shear wave along the direction of the measuring line, and the fourth component wave is the component wave of the second pure shear wave in the direction perpendicular to the direction of the measuring line. Based on the third angle difference, the first component wave, the second component wave, the third component wave, and the fourth component wave, determine the fast shear wave profile corresponding to the fast shear wave and the slow shear wave profile corresponding to the slow shear wave.

6. The method according to claim 1, characterized in that, The first corrected shear wave field includes the SH shear wave field and the SV shear wave field; The step of determining the first corrected shear wave field corresponding to the first pure shear wave data based on the crack direction information, the fast shear wave, and the corrected slow shear wave includes: For any target seismic source within the work area, the direction of the shot-detector connection between the target seismic source and the target geophone is obtained, and the first angle difference between the crack direction and the shot-detector connection direction is determined. The target seismic source is the seismic source within the work area that generates the first pure shear wave data. Based on the first angle difference, the fast shear wave and the corrected slow shear wave are reverse-rotated to obtain the SH shear wave field corresponding to the first pure shear wave data and the SV shear wave field corresponding to the first pure shear wave data.

7. The method according to claim 1, characterized in that, The method further includes: When the crack direction information is used to represent the change of crack direction with depth, multiple crack time windows corresponding to the crack direction information are determined. Acquire the second pure shear wave data corresponding to the multiple crack time windows, as well as the target time difference and target crack direction corresponding to each crack time window; Based on the target time difference and target crack direction corresponding to each crack time window, the second pure shear wave data is subjected to shear wave splitting correction to obtain the second corrected shear wave field corresponding to the second pure shear wave data.

8. A transverse wave splitting correction device for pure transverse waves, characterized in that, The device includes: The acquisition module is used to acquire information on the direction of fractures in the strata of the work area to be tested, and to acquire the first pure shear wave data received by the target geophone in the work area. The first determining module is used to determine the fast shear wave and slow shear wave corresponding to the first pure shear wave data according to the crack direction information when the crack direction information is used to indicate that the crack direction is a fixed direction. The fast shear wave is a shear wave parallel to the crack direction, and the slow shear wave is a shear wave perpendicular to the crack direction. The second determining module is used to determine the time difference between the fast shear wave and the slow shear wave, and to correct the slow shear wave according to the time difference to obtain a corrected slow shear wave. The third determining module is used to determine the first corrected shear wave field corresponding to the first pure shear wave data based on the crack direction information, the fast shear wave, and the corrected slow shear wave.

9. A computer device, characterized in that, The computer device includes: A processor and a memory, wherein the memory stores at least one piece of program code, which is loaded and executed by the processor to perform the operations performed in the shear wave splitting correction method for pure shear waves according to any one of claims 1 to 7.

10. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores at least one piece of program code, which is loaded and executed by a processor to perform the operations performed in the shear wave splitting correction method for pure shear waves as described in any one of claims 1 to 7.

Images