Method, device, electronic equipment and medium for positioning microseismic longitudinal and transverse waves in inclined wells

Abstract

The present application discloses a method, device, electronic device and medium for positioning microseismic longitudinal and transverse waves in an inclined well. The method may include: calculating the propagation direction of the longitudinal wave signal in the three-component horizontal plane of the first monitoring well, and determining the ray equation of the detector passing through the first monitoring well; calculating the propagation direction of the longitudinal wave signal in the three-component horizontal plane of the second monitoring well, and determining the ray equation of the detector passing through the second monitoring well; establishing an overdetermined set of equations based on the ray equation of the detector passing through the first monitoring well and the ray equation of the detector passing through the second monitoring well; solving the overdetermined set of equations, and calculating the positioning result of the microseism in the horizontal coordinate system; establishing a microseismic depth grid range, calculating the position of the depth grid point, and obtaining the spatial positioning result of the longitudinal and transverse waves of the microseismic in the inclined well. The present invention adopts a simultaneous monitoring method of two inclined wells with arbitrary well trajectories. In horizontal positioning, all the detector azimuths are used to constrain each other, and the longitudinal and transverse wave travel time is used as the target equation, which can enhance the positioning stability.

Description

Technical Field

[0001] The present invention relates to the field of microseismic signal processing in wells, and more specifically, to a method, device, electronic equipment and medium for positioning microseismic longitudinal and transverse waves in inclined wells. Background Art

[0002] Well microseismic monitoring is one of the microseismic observation methods. Its characteristic is that the three-component geophones in the well receive the full wavefield signal of microseismic. Compared with the ground microseismic monitoring, the data received in the well has a higher signal-to-noise ratio and a richer number and type of microseismic events. However, due to the limited number of microseismic geophones in the well (generally 12 to 32 three-component geophones in the well), the monitoring range is small. At the same time, there are certain requirements for the selection of monitoring wells and observation wells. For example, the optimal radial distance from the fracturing section of the monitoring well to the geophone of the observation well is 200 to 800 meters internationally, but the actual distance is often larger, which usually leads to unstable and low-precision microseismic positioning phenomena.

[0003] At present, the main methods for microseismic positioning in wells are: one is based on the forward modeling of the travel time of P-wave and S-wave events. Representative algorithms include network search method, simulated annealing method, Geiger method, etc. The advantage is that it is easy to implement, and the disadvantage is that the travel time of P-wave and S-wave of microseismic events is difficult to pick up accurately, which affects the positioning results; the second is based on wave equation convolution. Representative algorithms include interference method, reverse time migration method, and passive source imaging method. The advantage is that it does not need to pick up the first arrival of the event. The disadvantage is that it has high requirements on the data signal-to-noise ratio and velocity model, requires a large number of detectors, and has high calculation cost; the third is the difference between anisotropic and isotropic travel time calculations. In anisotropic media, the error of isotropic travel time calculation is large, and the corresponding positioning error is also large.

[0004] In recent years, with the continuous increase in the world's demand for energy and the advancement of exploration technology, the exploration and development of oil and gas resources has continued to develop in depth. As the most effective measure to increase the production of such reservoirs, fracturing has also received more and more attention at home and abroad. Microseismic monitoring in wells is often a single well monitoring method. Due to the small number of geophones, the positioning angle is small, resulting in unstable horizontal positioning results of microseismic in wells.

[0005] Therefore, it is necessary to develop a method, device, electronic equipment and medium for the joint rapid positioning of microseismic longitudinal and transverse waves for simultaneous in-well monitoring of two inclined wells.

[0006] The information disclosed in the background technology section of the present invention is only intended to deepen the understanding of the general background technology of the present invention, and should not be regarded as acknowledging or suggesting in any form that the information constitutes the prior art already known to those skilled in the art. Summary of the invention

[0007] The present invention proposes a method, device, electronic equipment and medium for positioning microseismic longitudinal and transverse waves in an inclined well. The method adopts a simultaneous monitoring mode of two inclined wells with arbitrary well trajectories. In horizontal positioning, the azimuth angles of all detectors are mutually constrained, and the longitudinal and transverse wave travel times are used as target equations, which can enhance positioning stability.

[0008] In a first aspect, the present disclosure provides a method for microseismic longitudinal and transverse wave positioning in a deviated well, comprising:

[0009] Calculate the propagation direction of the longitudinal wave signal of the microseismic event in the three-component horizontal plane of each geophone of the first monitoring well, and determine the ray equation of the geophone through the first monitoring well;

[0010] Calculate the propagation direction of the longitudinal wave signal of the microseismic event in the three-component horizontal plane of each detector in the second monitoring well, and determine the ray equation passing through the detector in the second monitoring well;

[0011] Establishing an overdetermined set of equations according to the ray equation passing through the first monitoring well detector and the ray equation passing through the second monitoring well detector;

[0012] Solving the overdetermined equations, calculating the vertical and horizontal coordinates, which are the positioning results of the microseismic movement in the horizontal coordinate system;

[0013] The microseismic depth grid range is established, and the depth grid point position is calculated according to the longitudinal and transverse coordinates to determine the microseismic depth direction positioning result, thereby obtaining the inclined well microseismic longitudinal and transverse wave spatial positioning result.

[0014] Preferably, the ray equation passing through the first monitoring well geophone is:

[0015] yy 1,i =tan(θ 1,i )*(xx 1,i ) (1)

[0016] Wherein, M is the number of detectors in the first monitoring well.

[0017] Preferably, the ray equation passing through the second monitoring well geophone is:

[0018] yy 2,j =tan(θ 2,j )*(xx 2,j ) (2)

[0019] Wherein, N is the number of detectors in the second monitoring well.

[0020] Preferably, the overdetermined system of equations is:

[0021]

[0022] Among them, x 1,i ,y1,i ,θ 1,i are the horizontal coordinate, vertical coordinate and propagation direction angle of the longitudinal wave signal of the microseismic event in the horizontal plane coordinate system of the i-th detector in the first monitoring well, respectively. There are M detectors in total, x 1,j ,y 1,j ,θ 1,j are the horizontal coordinate, vertical coordinate and propagation direction angle of the longitudinal wave signal of the microseismic event in the horizontal plane coordinate system of the j-th geophone in the second monitoring well, respectively. There are N geophones in total.

[0023] Preferably, calculating the depth grid point positions and determining the microseismic depth direction positioning results includes:

[0024] Calculate the arrival time of microseismic P-wave and S-wave signals from each grid point to the first monitoring well and the second monitoring well geophone, compare with the actual microseismic P-wave and S-wave travel time, and calculate the target error of the P-wave and S-wave combined travel time;

[0025] The depth grid point position corresponding to the minimum target error is the microseismic depth positioning result.

[0026] Preferably, the combined travel time target error of longitudinal wave and shear wave is calculated by formula (4):

[0027]

[0028] Where, ΔT PS,k is the target error of the combined travel time of longitudinal wave and shear wave, T 1,P,k is the arrival time of the microseismic P-wave signal from the grid point to the first monitoring well geophone, T 1,S,k is the arrival time of the microseismic shear wave signal from the grid point to the first monitoring well geophone, T 2,P,k is the arrival time of the microseismic P-wave signal from the grid point to the second monitoring well geophone, T 2,S,k is the arrival time of the microseismic shear wave signal from the grid point to the second monitoring well geophone, T 0 1,P is the known actual picked up microseismic P-wave travel time of the first monitoring well, T 0 1,S is the known actual picked up shear wave travel time of the first monitoring well, T 0 2,P is the known actual microseismic P-wave travel time of the second monitoring well, T 0 2,S It is the travel time of the shear wave actually picked up in the second monitoring well.

[0029] Preferably, the depth grid point position is calculated by formula (5):

[0030] z k =z0+(k-1)*Δz (5)

[0031] Among them, z k is the depth value corresponding to the grid point k, Δz is the unit grid size, z1 and z2 are the upper and lower limits of the search range of the known fracturing stage perforation depth z0, respectively, k ∈(z0-z1,z0+z2).

[0032] As a specific implementation method of the embodiment of the present disclosure,

[0033] In a second aspect, the present disclosure also provides a device for positioning microseismic longitudinal and transverse waves of a deviated well, comprising:

[0034] A first monitoring well geophone ray equation construction module is used to calculate the propagation direction of the microseismic event longitudinal wave signal in the three-component horizontal plane of each geophone in the first monitoring well, and determine the geophone ray equation through the first monitoring well;

[0035] A second monitoring well detector ray equation construction module is used to calculate the propagation direction of the microseismic event longitudinal wave signal in the three-component horizontal plane of each detector in the second monitoring well, and determine the second monitoring well detector ray equation;

[0036] An overdetermined equation group building module is used to build an overdetermined equation group according to the ray equation passing through the first monitoring well detector and the ray equation passing through the second monitoring well detector;

[0037] A horizontal positioning result calculation module solves the overdetermined equations and calculates the vertical and horizontal coordinates, which are the positioning results of the microseismic movement in the horizontal coordinate system;

[0038] The spatial positioning result calculation module establishes the microseismic depth grid range, calculates the depth grid point position according to the longitudinal and transverse coordinates, determines the microseismic depth direction positioning result, and then obtains the inclined well microseismic longitudinal and transverse wave spatial positioning result.

[0039] Preferably, the ray equation passing through the first monitoring well geophone is:

[0040] yy 1,i =tan(θ 1,i )*(xx 1,i ) (1)

[0041] Wherein, M is the number of detectors in the first monitoring well.

[0042] Preferably, the ray equation passing through the second monitoring well geophone is:

[0043] yy 2,j =tan(θ 2,j )*(xx 2,j ) (2)

[0044] Wherein, N is the number of detectors in the second monitoring well.

[0045] Preferably, the overdetermined system of equations is:

[0046]

[0047] Among them, x 1,i ,y 1,i ,θ 1,i are the horizontal coordinate, vertical coordinate and propagation direction angle of the longitudinal wave signal of the microseismic event in the horizontal plane coordinate system of the i-th detector in the first monitoring well, respectively. There are M detectors in total, x 1,j ,y 1,j ,θ 1,j are the horizontal coordinate, vertical coordinate and propagation direction angle of the longitudinal wave signal of the microseismic event in the horizontal plane coordinate system of the j-th geophone in the second monitoring well, respectively. There are N geophones in total.

[0048] Preferably, calculating the depth grid point positions and determining the microseismic depth direction positioning results includes:

[0049] Calculate the arrival time of microseismic P-wave and S-wave signals from each grid point to the first monitoring well and the second monitoring well geophone, compare with the actual microseismic P-wave and S-wave travel time, and calculate the target error of the P-wave and S-wave combined travel time;

[0050] The depth grid point position corresponding to the minimum target error is the microseismic depth positioning result.

[0051] Preferably, the combined travel time target error of longitudinal wave and shear wave is calculated by formula (4):

[0052]

[0053] Where, ΔT PS,k is the target error of the combined travel time of longitudinal wave and shear wave, T 1,P,k is the arrival time of the microseismic P-wave signal from the grid point to the first monitoring well geophone, T 1,S,k is the arrival time of the microseismic shear wave signal from the grid point to the first monitoring well geophone, T 2,P,k is the arrival time of the microseismic P-wave signal from the grid point to the second monitoring well geophone, T 2,S,k is the arrival time of the microseismic shear wave signal from the grid point to the second monitoring well geophone, T 0 1,P is the known actual picked up microseismic P-wave travel time of the first monitoring well, T 0 1,S is the known actual picked up shear wave travel time of the first monitoring well, T 0 2,P is the known actual microseismic P-wave travel time of the second monitoring well, T 0 2,SIt is the travel time of the shear wave actually picked up in the second monitoring well.

[0054] Preferably, the depth grid point position is calculated by formula (5):

[0055] z k =z0+(k-1)*Δz (5)

[0056] Among them, z k is the depth value corresponding to the grid point k, Δz is the unit grid size, z1 and z2 are the upper and lower limits of the search range of the known fracturing stage perforation depth z0, respectively, k ∈(z0-z1,z0+z2).

[0057] In a third aspect, an embodiment of the present disclosure further provides an electronic device, the electronic device comprising:

[0058] A memory storing executable instructions;

[0059] A processor runs the executable instructions in the memory to implement the inclined well microseismic longitudinal and transverse wave positioning method.

[0060] In a fourth aspect, an embodiment of the present disclosure further provides a computer-readable storage medium, which stores a computer program, and when the computer program is executed by a processor, the method for locating microseismic longitudinal and transverse waves in inclined wells is implemented.

[0061] Its beneficial effects are:

[0062] The methods and apparatus of the present invention have other features and advantages that will be apparent from, or will be described in detail in, the accompanying drawings and subsequent detailed descriptions incorporated herein, which together serve to explain the specific principles of the present invention. BRIEF DESCRIPTION OF THE DRAWINGS

[0063] The above and other objects, features and advantages of the present invention will become more apparent through a more detailed description of exemplary embodiments of the present invention in conjunction with the accompanying drawings, wherein like reference numerals generally represent like components throughout the exemplary embodiments of the present invention.

[0064] Figure 1 A flow chart showing the steps of a method for locating microseismic longitudinal and transverse waves in a deviated well according to an embodiment of the present invention.

[0065] Figure 2a and Figure 2b Schematic diagrams showing a side view and a top view of a system for simultaneous microseismic monitoring of two inclined wells according to an embodiment of the present invention are shown respectively.

[0066] Figure 3 A schematic diagram showing a top view of the ray equation path of the same microseismic event passing through two monitoring wells according to an embodiment of the present invention.

[0067] Figure 4a and Figure 4b Schematic diagrams showing a side view and a top view of a three-dimensional spatial positioning result of a microseismic event according to an embodiment of the present invention are shown respectively.

[0068] Figure 5 A block diagram of a device for locating microseismic longitudinal and transverse waves of a deviated well according to an embodiment of the present invention is shown.

[0069] Description of reference numerals:

[0070] 201. A module for constructing the ray equation of the first monitoring well detector; 202. A module for constructing the ray equation of the second monitoring well detector; 203. A module for constructing an overdetermined set of equations; 204. A module for calculating the horizontal positioning result; 205. A module for calculating the spatial positioning result. DETAILED DESCRIPTION

[0071] The preferred embodiments of the present invention will be described in more detail below. Although the preferred embodiments of the present invention are described below, it should be understood that the present invention can be implemented in various forms and should not be limited to the embodiments set forth herein.

[0072] Compared with ground microseismic monitoring, the acquisition cost of microseismic data in wells is lower, and it is less affected by weather and surrounding environmental noise. The monitoring distance is closer than that of ground microseismic data, and it can identify more microseismic events and enrich the subsequent microseismic fracturing interpretation. However, due to the small number of detectors, the monitoring distance exceeds the ideal range, and the insufficient depth below, the horizontal azimuth positioning results of microseismic data in a single well are unstable and inaccurate. In order to meet the production market's demand for high-precision microseismic positioning in wells, two wells are monitored simultaneously and mutually constrained to improve the horizontal positioning accuracy. The present invention utilizes the polarization principle of P-wave signals to respectively calculate the detector azimuths of the same microseismic event in two monitoring wells, removes the maximum and minimum azimuths, and the remaining detector azimuths and detector coordinates form a series of ray equations to construct an overdetermined equation group. The vertical and horizontal coordinates calculated by solving the binary linear overdetermined equations are the horizontal positioning results of the microseismic events monitored by the two inclined wells. The grid search method is then used to further search for the grid points with the smallest travel time errors with the longitudinal and transverse waves of the microseismic events of the two wells within the depth range, which are the depth positioning results of the microseismic events, thereby realizing rapid positioning of the microseismic events simultaneously monitored by the two inclined wells.

[0073] The present invention provides a method for positioning microseismic longitudinal and transverse waves in an inclined well, comprising:

[0074] The propagation direction of the longitudinal wave signal of the microseismic event in the three-component horizontal plane of each geophone of the first monitoring well is calculated, and the ray equation passing through the geophone of the first monitoring well is determined; in one example, the ray equation passing through the geophone of the first monitoring well is:

[0075] yy 1,i =tan(θ 1,i )*(xx 1,i ) (1)

[0076] Wherein, M is the number of detectors in the first monitoring well.

[0077] Specifically, according to the principle of longitudinal wave signal polarization: the propagation direction of the longitudinal wave is consistent with the polarization direction of the longitudinal wave in the horizontal plane of the three-component detector, and the polarization direction is represented by the horizontal azimuth of the three-component detector. The present invention adopts the existing algorithm - energy histogram analysis method, assuming that the original instantaneous amplitudes of the longitudinal wave signal of the microseismic event received by the X and Y components of the three-component detector in the first monitoring well are A and X,i , A Y,i , by rotating at an angle α, we obtain two new components - the horizontal radial component R and the horizontal tangential component T:

[0078] A R,i =A X,i cos(α)+A Y,i sin(α) (6)

[0079] A T,i =-A X,i sin(α)+A Y,i cos(α) (7)

[0080] Among them, i is the time sample point, R i , T i They are the instantaneous amplitudes of the horizontal radial component R and the horizontal tangential component T after rotation, and the angle α ranges from 0 to 360 degrees.

[0081] According to formulas (6) and (7), the energy difference ΔE between the horizontal radial component R and the horizontal tangential component T after rotation is established:

[0082] ΔE=∑|A R,i | 2 -∑|A T,i | 2 (8)

[0083] Statistically analyze the energy difference ΔE of angle α from 0 degrees to 360 degrees. When ΔE is the largest, the corresponding angle α is the horizontal azimuth of one of the three-component detectors in the first monitoring well that is consistent with the propagation direction of the longitudinal wave signal of the microseismic event. The same operation is performed to calculate the horizontal azimuth of all detectors in the first monitoring well, and the maximum and minimum horizontal azimuth of the detector can be found. The remaining detector azimuths θ 1,i The horizontal coordinate of its detector (x 1,i ,y 1,i ) forms a series of ray equations L 1,i is formula (1).

[0084] The propagation direction of the longitudinal wave signal of the microseismic event in the three-component horizontal plane of each detector in the second monitoring well is calculated to determine the ray equation of the detector through the second monitoring well; in one example, the ray equation of the detector through the second monitoring well is:

[0085] yy 2,j =tan(θ 2,j )*(xx 2,j ) (2)

[0086] Wherein, N is the number of detectors in the second monitoring well.

[0087] Specifically, the same operation process is used to calculate the azimuths of all the geophones in the second monitoring well for the same microseismic event P-wave signal by energy histogram analysis. The maximum and minimum horizontal azimuths of the geophones are removed, and the remaining geophone azimuths θ 2,j The horizontal coordinate of its detector (x 2,j ,y 2,j ), forming a series of ray equations L for the second monitoring well 2,j is formula (2).

[0088] An overdetermined set of equations is established based on the ray equations passing through the first monitoring well detector and the ray equations passing through the second monitoring well detector. In one example, the overdetermined set of equations is:

[0089]

[0090] Among them, x 1,i ,y 1,i ,θ 1,i are the horizontal coordinate, vertical coordinate and propagation direction angle of the longitudinal wave signal of the microseismic event in the horizontal plane coordinate system of the i-th detector in the first monitoring well, respectively. There are M detectors in total, x 1,j ,y 1,j ,θ 1,j are the horizontal coordinate, vertical coordinate and propagation direction angle of the longitudinal wave signal of the microseismic event in the horizontal plane coordinate system of the j-th geophone in the second monitoring well, respectively. There are N geophones in total.

[0091] Specifically, the horizontal positioning results of the microseismic events monitored simultaneously in two inclined wells are obtained by using the ray equations of the two wells monitoring the detectors. The ray equations of the first monitoring well from the same microseismic event are obtained by using the formulas (1) and (2) in the previous step: 1,i , the ray equation group L through the first monitoring well geophone 2,j , forming an overdetermined system of equations, which is formula (3).

[0092] Solve the overdetermined equations and calculate the vertical and horizontal coordinates, which are the positioning results of the micro-earthquake in the horizontal coordinate system;

[0093] Specifically, the two vertical and horizontal coordinate values ​​(x * ,y * ), which is the positioning result of the microseismic event in the horizontal coordinate system.

[0094] Establish a microseismic depth grid range, calculate the depth grid point position according to the vertical and horizontal coordinates, determine the microseismic depth direction positioning result, and then obtain the inclined well microseismic longitudinal and transverse wave spatial positioning result. In one example, calculating the depth grid point position and determining the microseismic depth direction positioning result includes: calculating the arrival time of the microseismic longitudinal and transverse wave signals from each grid point to the first monitoring well and the second monitoring well geophone, comparing with the actual microseismic longitudinal and transverse wave travel time, and calculating the target error of the longitudinal and transverse wave combined travel time; the corresponding depth grid point position when the target error is minimized is the microseismic depth direction positioning result.

[0095] In one example, the combined travel time target error of longitudinal wave and shear wave is calculated by formula (4):

[0096]

[0097] Where, ΔT PS,k is the target error of the combined travel time of longitudinal wave and shear wave, T 1,P,k is the arrival time of the microseismic P-wave signal from the grid point to the first monitoring well geophone, T 1,S,k is the arrival time of the microseismic shear wave signal from the grid point to the first monitoring well geophone, T 2,P,k is the arrival time of the microseismic P-wave signal from the grid point to the second monitoring well geophone, T 2,S,k is the arrival time of the microseismic shear wave signal from the grid point to the second monitoring well geophone, T 0 1,P is the known actual picked up microseismic P-wave travel time of the first monitoring well, T 0 1,S is the known actual picked up shear wave travel time of the first monitoring well, T 0 2,Pis the known actual microseismic P-wave travel time of the second monitoring well, T 0 2,S It is the travel time of the shear wave actually picked up in the second monitoring well.

[0098] In one example, the depth grid point position is calculated by formula (5):

[0099] z k =z0+(k-1)*Δz (5)

[0100] Among them, z k is the depth value corresponding to the grid point k, Δz is the unit grid size, z1 and z2 are the upper and lower limits of the search range of the known fracturing stage perforation depth z0, respectively, k ∈(z0-z1,z0+z2).

[0101] Specifically, through the grid search method, the depth positioning results of microseismic events monitored simultaneously in two inclined wells are obtained, and finally the rapid spatial positioning of microseismic events is achieved.

[0102] With the known perforation depth z0 of the fracturing stage as the center, the microseismic event depth grid is established by formula (5). Assume that each depth grid point z k Place a seismic source and combine the calculated microseismic event horizontal location results (x * ,y * ), based on the known acoustic logging and the spatial positions of the geophones in the two inclined monitoring wells, ray tracing forward modeling was performed using Snell's law to calculate the arrival travel time T of the microseismic P-wave signal from each grid point to the geophone of the first monitoring well. 1,P,k , Transverse wave signal arrival time T 1,S,k The arrival time of the microseismic P-wave signal from the second monitoring well geophone is T 2,P,k , Transverse wave signal arrival time T 2,S,k , and the known actual picked-up microseismic P-wave travel time T of the first monitoring well 0 1,P , transverse wave travel time T 0 1,S The microseismic P-wave travel time T0 of the second monitoring well 2,P , transverse wave travel time T 0 2,S , establish the target error formula of longitudinal and transverse wave travel time ΔT PS,k is formula (4).

[0103] Using the search method, the travel time error ΔT of formula (4) is searched within the grid range centered on the known perforation depth of the fracturing stage. PS,k The depth grid point z corresponding to the minimum * , which is the depth positioning result of the microseismic event, combined with the calculated horizontal positioning result (x *,y * ), and finally obtain the spatial positioning results of microseismic events monitored simultaneously in two inclined wells (x * ,y * ,z * ).

[0104] The present invention also provides a device for positioning microseismic longitudinal and transverse waves of inclined wells, comprising:

[0105] The first monitoring well detector ray equation construction module calculates the propagation direction of the microseismic event longitudinal wave signal in the three-component horizontal plane of each detector in the first monitoring well, and determines the detector ray equation through the first monitoring well; in one example, the detector ray equation through the first monitoring well is:

[0106] yy 1,i =tan(θ 1,i )*(xx 1,i ) (1)

[0107] Wherein, M is the number of detectors in the first monitoring well.

[0108] Specifically, according to the principle of longitudinal wave signal polarization: the propagation direction of the longitudinal wave is consistent with the polarization direction of the longitudinal wave in the horizontal plane of the three-component detector, and the polarization direction is represented by the horizontal azimuth of the three-component detector. The present invention adopts the existing algorithm - energy histogram analysis method, assuming that the original instantaneous amplitudes of the longitudinal wave signal of the microseismic event received by the X and Y components of the three-component detector in the first monitoring well are A and X,i , A Y,i , by rotating at an angle α, we obtain two new components: the horizontal radial component R and the horizontal tangential component T, which are formulas (6) and (7) respectively.

[0109] According to formulas (6) and (7), the energy difference ΔE between the horizontal radial component R and the horizontal tangential component T after rotation is established as formula (8).

[0110] Statistically analyze the energy difference ΔE of angle α from 0 degrees to 360 degrees. When ΔE is the largest, the corresponding angle α is the horizontal azimuth of one of the three-component detectors in the first monitoring well that is consistent with the propagation direction of the longitudinal wave signal of the microseismic event. The same operation is performed to calculate the horizontal azimuth of all detectors in the first monitoring well, and the maximum and minimum horizontal azimuth of the detector can be found. The remaining detector azimuths θ 1,i The horizontal coordinate of its detector (x 1,i ,y 1,i ) forms a series of ray equations L 1,i is formula (1).

[0111] The second monitoring well detector ray equation construction module calculates the propagation direction of the microseismic event longitudinal wave signal in the three-component horizontal plane of each detector in the second monitoring well, and determines the detector ray equation through the second monitoring well; in an example, the detector ray equation through the second monitoring well is:

[0112] yy 2,j =tan(θ 2,j )*(xx 2,j ) (2)

[0113] Wherein, N is the number of detectors in the second monitoring well.

[0114] Specifically, the same operation process is used to calculate the azimuths of all the geophones in the second monitoring well for the same microseismic event P-wave signal by energy histogram analysis. The maximum and minimum horizontal azimuths of the geophones are removed, and the remaining geophone azimuths θ 2,j The horizontal coordinate of its detector (x 2,j ,y 2,j ), forming a series of ray equations L for the second monitoring well 2,j is formula (2).

[0115] The overdetermined equation group building module establishes an overdetermined equation group according to the ray equation passing through the first monitoring well detector and the ray equation passing through the second monitoring well detector; in one example, the overdetermined equation group is:

[0116]

[0117] Among them, x 1,i ,y 1,i ,θ 1,i are the horizontal coordinate, vertical coordinate and propagation direction angle of the longitudinal wave signal of the microseismic event in the horizontal plane coordinate system of the i-th detector in the first monitoring well, respectively. There are M detectors in total, x 1,j ,y 1,j ,θ 1,j are the horizontal coordinate, vertical coordinate and propagation direction angle of the longitudinal wave signal of the microseismic event in the horizontal plane coordinate system of the j-th geophone in the second monitoring well, respectively. There are N geophones in total.

[0118] Specifically, the horizontal positioning results of the microseismic events monitored simultaneously in two inclined wells are obtained by using the ray equations of the two wells monitoring the detectors. The ray equations of the first monitoring well from the same microseismic event are obtained by using the formulas (1) and (2) in the previous step: 1,i , the ray equation group L through the first monitoring well geophone 2,j , forming an overdetermined system of equations, which is formula (3).

[0119] The horizontal positioning result calculation module solves the overdetermined equations and calculates the vertical and horizontal coordinates, which are the positioning results of the microseismic movement in the horizontal coordinate system.

[0120] Specifically, the two vertical and horizontal coordinate values ​​(x * ,y * ), which is the positioning result of the microseismic event in the horizontal coordinate system.

[0121] The spatial positioning result calculation module establishes the microseismic depth grid range, calculates the depth grid point position according to the vertical and horizontal coordinates, determines the microseismic depth direction positioning result, and then obtains the inclined well microseismic longitudinal and transverse wave spatial positioning result. In one example, calculating the depth grid point position and determining the microseismic depth direction positioning result includes: calculating the arrival time of the microseismic longitudinal and transverse wave signals from each grid point to the first monitoring well and the second monitoring well detector, comparing with the actual microseismic longitudinal and transverse wave travel time, and calculating the target error of the longitudinal and transverse wave combined travel time; the corresponding depth grid point position when the target error is minimized is the microseismic depth direction positioning result.

[0122] In one example, the combined travel time target error of longitudinal wave and shear wave is calculated by formula (4):

[0123]

[0124] Where, ΔT PS,k is the target error of the combined travel time of longitudinal wave and shear wave, T 1,P,k is the arrival time of the microseismic P-wave signal from the grid point to the first monitoring well geophone, T 1,S,k is the arrival time of the microseismic shear wave signal from the grid point to the first monitoring well geophone, T 2,P,k is the arrival time of the microseismic P-wave signal from the grid point to the second monitoring well geophone, T 2,S,k is the arrival time of the microseismic shear wave signal from the grid point to the second monitoring well geophone, T 0 1,P is the known actual picked up microseismic P-wave travel time of the first monitoring well, T 0 1,S is the known actual picked up shear wave travel time of the first monitoring well, T 0 2,P is the known actual microseismic P-wave travel time of the second monitoring well, T 0 2,S It is the travel time of the shear wave actually picked up in the second monitoring well.

[0125] In one example, the depth grid point position is calculated by formula (5):

[0126] z k =z0+(k-1)*Δz (5)

[0127] Among them, z k is the depth value corresponding to the grid point k, Δz is the unit grid size, z1 and z2 are the upper and lower limits of the search range of the known fracturing stage perforation depth z0, respectively, k ∈(z0-z1,z0+z2).

[0128] Specifically, through the grid search method, the depth positioning results of microseismic events monitored simultaneously in two inclined wells are obtained, and finally the rapid spatial positioning of microseismic events is achieved.

[0129] With the known perforation depth z0 of the fracturing stage as the center, the microseismic event depth grid is established by formula (5). Assume that each depth grid point z k Place a seismic source and combine the calculated microseismic event horizontal location results (x * ,y * ), based on the known acoustic logging and the spatial positions of the geophones in the two inclined monitoring wells, ray tracing forward modeling was performed using Snell's law to calculate the arrival travel time T of the microseismic P-wave signal from each grid point to the geophone of the first monitoring well. 1,P,k , Transverse wave signal arrival time T 1,S,k The arrival time of the microseismic P-wave signal from the second monitoring well geophone is T 2,P,k , Transverse wave signal arrival time T 2,S,k , and the known actual picked-up microseismic P-wave travel time T of the first monitoring well 0 1,P , transverse wave travel time T 0 1,S The microseismic P-wave travel time T of the second monitoring well 0 2,P , transverse wave travel time T 0 2,S , establish the target error formula of longitudinal and transverse wave travel time ΔT PS,k is formula (4).

[0130] Using the search method, the travel time error ΔT of formula (4) is searched within the grid range centered on the known perforation depth of the fracturing stage. PS,k The depth grid point z corresponding to the minimum * , which is the depth positioning result of the microseismic event, combined with the calculated horizontal positioning result (x * ,y * ), and finally obtain the spatial positioning results of microseismic events monitored simultaneously in two inclined wells (x * ,y * ,z * ).

[0131] The present invention also provides an electronic device, which includes: a memory storing executable instructions; and a processor, which runs the executable instructions in the memory to implement the above-mentioned inclined well microseismic longitudinal and transverse wave positioning method.

[0132] The present invention also provides a computer-readable storage medium, which stores a computer program. When the computer program is executed by a processor, the above-mentioned inclined well microseismic longitudinal and transverse wave positioning method is implemented.

[0133] To facilitate understanding of the solutions and effects of the embodiments of the present invention, four specific application examples are given below. Those skilled in the art should understand that the examples are only for facilitating understanding of the present invention, and any specific details thereof are not intended to limit the present invention in any way.

[0134] Example 1

[0135] Figure 1 A flow chart showing the steps of a method for locating microseismic longitudinal and transverse waves in a deviated well according to an embodiment of the present invention.

[0136] like Figure 1 As shown, the method for locating longitudinal and transverse waves of microseismic events in inclined wells includes: step 101, calculating the propagation direction of the longitudinal wave signal of the microseismic event in the three-component horizontal plane of each detector in the first monitoring well, and determining the ray equation of the detector passing through the first monitoring well; step 102, calculating the propagation direction of the longitudinal wave signal of the microseismic event in the three-component horizontal plane of each detector in the second monitoring well, and determining the ray equation of the detector passing through the second monitoring well; step 103, establishing an overdetermined equation group according to the ray equation of the detector passing through the first monitoring well and the ray equation of the detector passing through the second monitoring well; step 104, solving the overdetermined equation group, calculating the longitudinal and transverse coordinates, which is the positioning result of the microseism in the horizontal coordinate system; step 105, establishing the microseismic depth grid range, calculating the depth grid point position according to the longitudinal and transverse coordinates, determining the microseismic depth direction positioning result, and then obtaining the inclined well microseismic longitudinal and transverse wave spatial positioning result.

[0137] Figure 2a and Figure 2b Schematic diagrams showing a side view and a top view of a system for simultaneous microseismic monitoring of two inclined wells according to an embodiment of the present invention are shown respectively.

[0138] The first and second monitoring wells both have 13-stage three-component geophones deployed underground (the geophone coordinates are known), 1 known perforation signal and 9 event signals, and their spatial coordinates are shown in Table 1, in meters. The perforation horizontal coordinates are at the origin, and their geometric relationships are as follows: Figure 2a , Figure 2b As shown, the formation is divided into three layers. This system belongs to a typical dual-well simultaneous microseismic monitoring system. Figure 2a , Figure 2b The observation method is to perform the longitudinal wave travel time information and the shear wave travel time information of each perforation and event reaching the detector as the known real observation value. The present invention does not consider the influence of the lateral change of velocity as the known input.

[0139] Table 1

[0140] The actual X coordinate of the event The actual Y coordinate of the event The actual Z coordinate of the event -205 205 1260 -155 155 1270 -105 105 1280 -55 55 1290 -5 5 1300 45 -45 1305 95 -95 1310 145 -145 1315 195 -195 1320

[0141] Figure 3 A schematic diagram showing a top view of the ray equation path of the same microseismic event passing through two monitoring wells according to an embodiment of the present invention.

[0142] First, take the spatial coordinate (-5, 5, 1300) event as an example, according to the principle of longitudinal wave signal polarization: the longitudinal wave propagation direction is consistent with the longitudinal wave polarization direction in the horizontal plane of the three-component detector, and the polarization direction is represented by the horizontal azimuth of the three-component detector. The present invention uses the energy histogram analysis method to calculate the longitudinal wave of the event at different detector azimuths of the first monitoring well and different detector azimuths of the second monitoring well, as shown in Table 2, unit: degree. Search and remove the maximum azimuth and minimum azimuth of each of the two monitoring wells, and combine the corresponding detector horizontal coordinates to form a series of detector ray equations, such as Figure 3 shown.

[0143] Table 2

[0144]

[0145]

[0146] Then, due to the same microseismic event, Figure 3 A series of ray equations in the equation pass through the same event point theoretically. Therefore, through formulas (1) and (2), an overdetermined equation system (3) is formed. The ordinate and abscissa coordinate values ​​(x * ,y * ), which is the horizontal positioning result of microseismic events monitored simultaneously in two inclined wells.

[0147] Figure 4a and Figure 4b Schematic diagrams showing a side view and a top view of a three-dimensional spatial positioning result of a microseismic event according to an embodiment of the present invention are shown respectively.

[0148] Finally, with the known perforation depth z0 (depth of 1300 meters) of the fracturing stage as the center, a microseismic event depth grid range of 1200-1400 meters and a unit grid length of 1 meter are established, and the grid points of formula (5) are established. For each microseismic event, the grid point is taken as the depth, combined with the horizontal positioning result (x * ,y * ), assuming that a seismic source is placed at each depth grid point, using Snell's law, ray tracing forward modeling is performed to calculate the P-wave and S-wave travel time information of the microseismic detector of the first monitoring well from each grid point and the P-wave and S-wave travel time information of the microseismic detector of the second monitoring well from the point, and the P-wave and S-wave travel time information of the microseismic detector of the two inclined wells actually picked up, and the P-wave and S-wave travel time combination error formula (4) is established. Using the search method, when the P-wave and S-wave travel time combination error of the two inclined wells is the smallest, the corresponding depth grid point z * , which is the depth positioning result of the microseismic event, combined with the calculated horizontal positioning result (x * ,y * ), and finally obtain the spatial location distribution of microseismic events monitored simultaneously in two inclined wells (x * ,y * ,z * ),like Figure 4a , Figure 4b shown.

[0149] It can be seen from the positioning error table 3 that both the horizontal positioning error and the depth positioning error of the microseismic event are basically controlled at about 5 meters, which verifies that the present invention has a high microseismic positioning accuracy and provides technical support for the promotion and application of simultaneous microseismic monitoring of two inclined wells.

[0150] Table 3

[0151] Event x error (m) Event y error (m) Event z error (m) 1.562668 -0.800003 1 1.562119 0.40509 -5 -0.320854 0.949577 4 -3.137589 3.231972 0 3.317762 -0.576878 1 -0.116978 -2.033707 0 -1.303085 -4.124016 5 -0.691101 1.35376 -2 -0.753235 0.673416 2

[0152] Example 2

[0153] Figure 5 A block diagram of a device for locating microseismic longitudinal and transverse waves of a deviated well according to an embodiment of the present invention is shown.

[0154] like Figure 5 As shown, the inclined well microseismic longitudinal and transverse wave positioning device comprises:

[0155] The first monitoring well geophone ray equation construction module 201 calculates the propagation direction of the microseismic event longitudinal wave signal in the three-component horizontal plane of each geophone in the first monitoring well, and determines the geophone ray equation of the first monitoring well;

[0156] The second monitoring well detector ray equation construction module 202 calculates the propagation direction of the microseismic event longitudinal wave signal in the three-component horizontal plane of each detector in the second monitoring well, and determines the second monitoring well detector ray equation;

[0157] An overdetermined equation group building module 203 builds an overdetermined equation group according to a ray equation passing through the first monitoring well detector and a ray equation passing through the second monitoring well detector;

[0158] The horizontal positioning result calculation module 204 solves the overdetermined equations and calculates the vertical and horizontal coordinates, which are the positioning results of the microseismic movement in the horizontal coordinate system;

[0159] The spatial positioning result calculation module 205 establishes the microseismic depth grid range, calculates the depth grid point position according to the longitudinal and transverse coordinates, determines the microseismic depth direction positioning result, and then obtains the inclined well microseismic longitudinal and transverse wave spatial positioning result.

[0160] As an optional solution, the ray equation through the first monitoring well geophone is:

[0161] yy 1,i =tan(θ 1,i )*(xx 1,i ) (1)

[0162] Wherein, M is the number of detectors in the first monitoring well.

[0163] As an optional solution, the ray equation through the second monitoring well geophone is:

[0164] yy 2,j =tan(θ 2,j )*(xx 2,j ) (2)

[0165] Wherein, N is the number of detectors in the second monitoring well.

[0166] As an alternative, the overdetermined system of equations is:

[0167]

[0168] Among them, x 1,i ,y 1,i ,θ 1,i are the horizontal coordinate, vertical coordinate and propagation direction angle of the longitudinal wave signal of the microseismic event in the horizontal plane coordinate system of the i-th detector in the first monitoring well, respectively. There are M detectors in total, x 1,j ,y 1,j ,θ 1,j are the horizontal coordinate, vertical coordinate and propagation direction angle of the longitudinal wave signal of the microseismic event in the horizontal plane coordinate system of the j-th geophone in the second monitoring well, respectively. There are N geophones in total.

[0169] As an optional solution, the depth grid point positions are calculated to determine the microseismic depth direction positioning results including:

[0170] Calculate the arrival time of microseismic P-wave and S-wave signals from each grid point to the first monitoring well and the second monitoring well geophone, compare with the actual microseismic P-wave and S-wave travel time, and calculate the target error of the P-wave and S-wave combined travel time;

[0171] The depth grid point position corresponding to the minimum target error is the microseismic depth positioning result.

[0172] As an optional solution, the combined travel time target error of longitudinal wave and shear wave is calculated by formula (4):

[0173]

[0174] Where, ΔT PS,k is the target error of the combined travel time of longitudinal wave and shear wave, T 1,P,k is the arrival time of the microseismic P-wave signal from the grid point to the first monitoring well geophone, T 1,S,k is the arrival time of the microseismic shear wave signal from the grid point to the first monitoring well geophone, T 2,P,k is the arrival time of the microseismic P-wave signal from the grid point to the second monitoring well geophone, T 2,S,k is the arrival time of the microseismic shear wave signal from the grid point to the second monitoring well geophone, T 0 1,P is the known actual picked up microseismic P-wave travel time of the first monitoring well, T 0 1,S is the known actual picked up shear wave travel time of the first monitoring well, T 0 2,P is the known actual microseismic P-wave travel time of the second monitoring well, T 0 2,S It is the travel time of the shear wave actually picked up in the second monitoring well.

[0175] As an alternative, the depth grid point position is calculated by formula (5):

[0176] z k =z0+(k-1)*Δz (5)

[0177] Among them, z k is the depth value corresponding to the grid point k, Δz is the unit grid size, z1 and z2 are the upper and lower limits of the search range of the known fracturing stage perforation depth z0, respectively, k ∈(z0-z1,z0+z2).

[0178] Example 3

[0179] The present disclosure provides an electronic device, which includes: a memory storing executable instructions; a processor running the executable instructions in the memory to implement the above-mentioned inclined well microseismic longitudinal and transverse wave positioning method.

[0180] An electronic device according to an embodiment of the present disclosure includes a memory and a processor.

[0181] The memory is used to store non-temporary computer-readable instructions. Specifically, the memory may include one or more computer program products, which may include various forms of computer-readable storage media, such as volatile memory and / or non-volatile memory. The volatile memory may, for example, include random access memory (RAM) and / or cache memory (cache), etc. The non-volatile memory may, for example, include read-only memory (ROM), hard disk, flash memory, etc.

[0182] The processor may be a central processing unit (CPU) or other forms of processing units having data processing capabilities and / or instruction execution capabilities, and may control other components in the electronic device to perform desired functions. In one embodiment of the present disclosure, the processor is used to run the computer-readable instructions stored in the memory.

[0183] Those skilled in the art should be able to understand that in order to solve the technical problem of how to obtain a good user experience, the present embodiment may also include well-known structures such as a communication bus and an interface, and these well-known structures should also be included in the protection scope of the present disclosure.

[0184] For detailed description of this embodiment, reference may be made to the corresponding descriptions in the aforementioned embodiments, which will not be repeated here.

[0185] Example 4

[0186] An embodiment of the present disclosure provides a computer-readable storage medium storing a computer program, which implements the inclined well microseismic longitudinal and transverse wave positioning method when executed by a processor.

[0187] According to the computer-readable storage medium of the embodiment of the present disclosure, non-transitory computer-readable instructions are stored thereon. When the non-transitory computer-readable instructions are executed by a processor, all or part of the steps of the above-mentioned methods of each embodiment of the present disclosure are executed.

[0188] The above-mentioned computer-readable storage media include, but are not limited to: optical storage media (e.g., CD-ROM and DVD), magneto-optical storage media (e.g., MO), magnetic storage media (e.g., magnetic tape or mobile hard disk), media with built-in rewritable non-volatile memory (e.g., memory card) and media with built-in ROM (e.g., ROM box).

[0189] Those skilled in the art should understand that the purpose of the above description of the embodiments of the present invention is only to exemplarily illustrate the beneficial effects of the embodiments of the present invention, and is not intended to limit the embodiments of the present invention to any given examples.

[0190] The embodiments of the present invention have been described above, and the above description is exemplary, not exhaustive, and is not limited to the disclosed embodiments. Many modifications and changes will be apparent to those skilled in the art without departing from the scope and spirit of the described embodiments.

Claims

1. A method for positioning microseismic longitudinal and transverse waves in a deviated well, characterized in that: include: Calculate the propagation direction of the longitudinal wave signal of the microseismic event in the three-component horizontal plane of each geophone of the first monitoring well, and determine the ray equation of the geophone through the first monitoring well; Calculate the propagation direction of the longitudinal wave signal of the microseismic event in the three-component horizontal plane of each detector in the second monitoring well, and determine the ray equation passing through the detector in the second monitoring well; Establishing an overdetermined set of equations according to the ray equation passing through the first monitoring well detector and the ray equation passing through the second monitoring well detector; Solving the overdetermined equations, calculating the vertical and horizontal coordinates, which are the positioning results of the microseismic movement in the horizontal coordinate system; Establishing a microseismic depth grid range, calculating the depth grid point position according to the longitudinal and transverse coordinates, determining the microseismic depth direction positioning result, and then obtaining the inclined well microseismic longitudinal and transverse wave spatial positioning result; Determining the ray equation of the first monitoring well detector includes: calculating the horizontal azimuth of all the detectors of the first monitoring well, removing the maximum horizontal azimuth of the detector and the minimum horizontal azimuth of the detector, and forming a series of ray equation groups with the remaining detector azimuths and their detector horizontal coordinates and ordinates; Determining the detector ray equation of the second monitoring well includes: calculating the horizontal azimuths of all detectors of the second monitoring well, removing the maximum horizontal azimuth of the detector and the minimum horizontal azimuth of the detector, and forming a series of ray equations for the second monitoring well with the remaining detector azimuths and their detector horizontal coordinates and vertical coordinates.

2. The method for positioning microseismic longitudinal and transverse waves in inclined wells according to claim 1, wherein: The ray equation passing through the first monitoring well geophone is: yy 1,i =tan(θ 1,i )*(xx 1,i ) (1) Among them, x 1,i ,y 1,i ,θ 1,i are the horizontal coordinate, vertical coordinate and propagation azimuth of the longitudinal wave signal of the microseismic event in the horizontal plane coordinate system of the i-th detector in the first monitoring well.

3. The method for positioning microseismic longitudinal and transverse waves in inclined wells according to claim 2, wherein: The ray equation through the second monitoring well geophone is: yy 2,j =tan(θ 2,j )*(xx 2,j ) (2) Among them, x 2,j ,y 2,j ,θ 2,j are the horizontal coordinate, vertical coordinate and propagation azimuth of the longitudinal wave signal of the microseismic event in the horizontal plane coordinate system of the j-th detector in the second monitoring well.

4. The method for positioning microseismic longitudinal and transverse waves in inclined wells according to claim 3, wherein: The overdetermined system of equations is: Among them, the first monitoring well has a total of M detectors, and the second monitoring well has a total of N detectors.

5. The method for positioning microseismic longitudinal and transverse waves in inclined wells according to claim 1, wherein: Calculate the depth grid point position and determine the microseismic depth direction positioning results including: Calculate the arrival time of microseismic P-wave and S-wave signals from each grid point to the first monitoring well and the second monitoring well geophone, compare with the actual microseismic P-wave and S-wave travel time, and calculate the target error of the P-wave and S-wave combined travel time; The depth grid point position corresponding to the minimum target error is the microseismic depth positioning result.

6. The method for positioning microseismic longitudinal and transverse waves in inclined wells according to claim 5, wherein: The target error of the combined travel time of longitudinal wave and shear wave is calculated by formula (4): Where, ΔT PS,k is the target error of the combined travel time of longitudinal wave and shear wave, T 1,P,k is the arrival time of the microseismic P-wave signal from the grid point to the first monitoring well geophone, T 1,S,k is the arrival time of the microseismic shear wave signal from the grid point to the first monitoring well geophone, T 2,P,k is the arrival time of the microseismic P-wave signal from the grid point to the second monitoring well geophone, T 2,S,k is the arrival time of the microseismic shear wave signal from the grid point to the second monitoring well geophone, T 0 1,P is the known actual picked up microseismic P-wave travel time of the first monitoring well, T 0 1,S is the known actual picked up shear wave travel time of the first monitoring well, T 0 2,P is the known actual microseismic P-wave travel time of the second monitoring well, T 0 2,S It is the travel time of the shear wave actually picked up in the second monitoring well.

7. The method for positioning microseismic longitudinal and transverse waves in inclined wells according to claim 5, wherein: The depth grid point position is calculated by formula (5): from k =z0+(k-1)*Δz (5) Among them, z k is the depth value corresponding to the grid point k, Δz is the unit grid size, z1 and z2 are the upper and lower limits of the search range of the known fracturing stage perforation depth z0, respectively, k ∈(z0-z1,z0+z2).

8. A microseismic longitudinal and transverse wave positioning device for inclined wells, characterized in that: include: A first monitoring well geophone ray equation construction module is used to calculate the propagation direction of the microseismic event longitudinal wave signal in the three-component horizontal plane of each geophone in the first monitoring well, and determine the geophone ray equation through the first monitoring well; A second monitoring well detector ray equation construction module is used to calculate the propagation direction of the microseismic event longitudinal wave signal in the three-component horizontal plane of each detector in the second monitoring well, and determine the second monitoring well detector ray equation; An overdetermined equation group building module is used to build an overdetermined equation group according to the ray equation passing through the first monitoring well detector and the ray equation passing through the second monitoring well detector; A horizontal positioning result calculation module solves the overdetermined equations and calculates the vertical and horizontal coordinates, which are the positioning results of the microseismic movement in the horizontal coordinate system; The spatial positioning result calculation module establishes the microseismic depth grid range, calculates the depth grid point position according to the longitudinal and transverse coordinates, determines the microseismic depth direction positioning result, and then obtains the inclined well microseismic longitudinal and transverse wave spatial positioning result; Determining the ray equation of the first monitoring well detector includes: calculating the horizontal azimuth of all the detectors of the first monitoring well, removing the maximum horizontal azimuth of the detector and the minimum horizontal azimuth of the detector, and forming a series of ray equation groups with the remaining detector azimuths and their detector horizontal coordinates and ordinates; Determining the detector ray equation of the second monitoring well includes: calculating the horizontal azimuths of all detectors of the second monitoring well, removing the maximum horizontal azimuth of the detector and the minimum horizontal azimuth of the detector, and forming a series of ray equations for the second monitoring well with the remaining detector azimuths and their detector horizontal coordinates and vertical coordinates.

9. An electronic device, characterized in that: The electronic device comprises: A memory storing executable instructions; A processor, wherein the processor runs the executable instructions in the memory to implement the inclined well microseismic longitudinal and transverse wave positioning method according to any one of claims 1 to 7.

10. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the inclined well microseismic longitudinal and transverse wave positioning method described in any one of claims 1 to 7 is implemented.

Images