Quadrilateral Hybrid Location Method for Microseismic Events in Wells, Electronic Devices and Media

Abstract

The present application discloses a quadrilateral hybrid positioning method for in-well microseismic events, an electronic device and a medium. The method may include: Step 1: respectively establish two P-wave non-linear positioning equations and two S-wave non-linear positioning equations; Step 2: calculate the initial P-wave travel time positioning result and the initial S-wave travel time positioning result of the microseismic event; Step 3: construct a quadrilateral according to the initial P-wave travel time positioning result and the initial S-wave travel time positioning result, and calculate the center point of the quadrilateral; Step 4: calculate the new initial P-wave travel time positioning result and the new initial S-wave travel time positioning result according to the center point coordinates, and repeat Steps 3-4 until the center point coordinates converge to determine the final positioning result. The present invention realizes stable and rapid positioning of microseismic events through continuous iteration of simulating the center of the quadrilateral.

Description

Technical Field

[0001] The present invention relates to the field of borehole microseismic signal processing, and more particularly, to a quadrilateral hybrid positioning method for borehole microseismic events, an electronic device, and a medium. Background Art

[0002] Borehole microseismic monitoring is one of the microseismic observation methods. Its characteristic is that downhole three-component geophones receive the full wavefield signals of microseisms. Compared with surface microseismic monitoring, the data received in the borehole has a higher signal-to-noise ratio, and the number and types of microseismic events are richer. However, due to the limited number of borehole microseismic geophones (generally 12 - 32 levels of three-component borehole geophones), the monitoring range is small, and 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 geophones of the observation well is 200 - 800 meters internationally, but the actual distance is often relatively large, which usually leads to unstable and inaccurate microseismic positioning phenomena.

[0003] In addition, the positioning accuracy of microseismic events is also related to the positioning method. The main borehole microseismic positioning methods are as follows: One is based on the forward modeling of the travel times of P-wave and S-wave events. Representative algorithms include the network search method, the simulated annealing method, the geiger method, etc. The advantage is that it is easy to implement, and the disadvantage is that it is difficult to accurately pick up the travel times of P-wave and S-wave of microseismic events, which affects the positioning result. The second is based on wave equation convolution. Representative algorithms include the interference method, the reverse time migration method, and the passive source imaging method. The advantage is that it does not require picking up the first arrival of events, and the disadvantages are that it has high requirements for the signal-to-noise ratio of data, the velocity model, and the number of geophones, and the calculation cost is high. The third is the difference between the calculation of travel times in anisotropic and isotropic media. In anisotropic media, the calculation error of using isotropic travel times is relatively large, and the corresponding positioning error is also relatively large.

[0004] In recent years, with the continuous increase in the world's demand for energy and the progress of exploration technologies, the exploration and development of oil and gas resources have been continuously developed in depth. As the most effective stimulation measure for such reservoirs, fracturing has also received more and more extensive attention at home and abroad. Borehole microseismic monitoring is often a single-well monitoring method. Also, due to the small number of geophones, the positioning aperture angle is small, resulting in unstable borehole microseismic positioning results. Therefore, for the unstable single microseismic positioning result, the linear or non-linear positioning of P-wave and S-wave is usually used to mutually constrain, so that the positioning process is both stable and a high-precision positioning result can be obtained.

[0005] Therefore, it is necessary to develop a quadrilateral hybrid positioning method for borehole microseismic events, an electronic device, and a medium that simulate the continuous iteration of the quadrilateral center.

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

[0007] The present invention provides a quadrilateral hybrid positioning method for in-well microseismic events, an electronic device, and a medium. It can construct a quadrilateral through two sets of P-wave travel-time non-linear positioning and two sets of S-wave travel-time non-linear positioning, calculate the position of the center point of the quadrilateral, and then construct a new quadrilateral with the center point and the four positioning results, calculate the position of the new center point, compare the two center points, if not convergent, continue to construct new quadrilaterals and calculate the center points until convergent, and output the final microseismic event positioning result, so as to achieve stable and rapid positioning of microseismic events.

[0008] In a first aspect, an embodiment of the present disclosure provides a quadrilateral hybrid positioning method for in-well microseismic events, including:

[0009] Step 1: Establish two P-wave non-linear positioning equations and two S-wave non-linear positioning equations respectively;

[0010] Step 2: Calculate the initial P-wave travel-time positioning result and the initial S-wave travel-time positioning result of the microseismic event;

[0011] Step 3: Construct a quadrilateral according to the initial P-wave travel-time positioning result and the initial S-wave travel-time positioning result, and calculate the center point of the quadrilateral;

[0012] Step 4: Calculate the new initial P-wave travel-time positioning result and the new initial S-wave travel-time positioning result according to the center point coordinates, and repeat Steps 3-4 until the center point coordinates converge, and determine the final positioning result.

[0013] Preferably, the P-wave non-linear positioning equation is:

[0014]

[0015]

[0016] where f(P1) is the first P-wave non-linear positioning equation, f(P2) is the second P-wave non-linear positioning equation, and ΔTp is the P-wave travel-time difference.

[0017] Preferably, the S-wave non-linear positioning equation is:

[0018]

[0019]

[0020] Among them, f(S1) is the first S-wave nonlinear positioning equation, f(S2) is the second S-wave nonlinear positioning equation, and ΔTs is the S-wave travel time difference.

[0021] Preferably, through the grid search method, when the P-wave time difference exponential equations are respectively minimized, the corresponding grid points are the P-wave travel time positioning results.

[0022] Preferably, through the grid search method, when the S-wave time difference exponential equations are respectively minimized, the corresponding grid points are the S-wave travel time positioning results.

[0023] Preferably, the initial P-wave travel time positioning results are (x P1 , y P1 ), (x P2 , y P2 ), the initial S-wave travel time positioning results are (x S1 , y S1 ), (x S2 , y S2 ), the center point coordinates are (x0, y0), then the new initial P-wave travel time positioning results are The new initial S-wave travel time positioning results are

[0024] Preferably, it is judged whether to converge through the error formula. If the error is less than or equal to 0.1, it converges.

[0025] Preferably, the error formula is:

[0026] Δ = |x * 0 - x0| + |y * 0 - y0| (5)

[0027] Among them, (x0, y0) are the center point coordinates, and (x * 0, y * 0) are the new center point coordinates.

[0028] As a specific implementation manner of the embodiments of the present disclosure,

[0029] In a second aspect, the embodiments of the present disclosure further provide an electronic device, and the electronic device includes:

[0030] A memory storing executable instructions;

[0031] A processor, and the processor runs the executable instructions in the memory to implement the wellbore microseismic event quadrilateral hybrid positioning method.

[0032] In a third aspect, embodiments of the present disclosure further provide a computer-readable storage medium storing a computer program, which when executed by a processor implements the wellbore microseismic event quadrilateral hybrid positioning method described above.

[0033] The method and apparatus of the present invention have other characteristics and advantages, which will be apparent from the accompanying drawings incorporated herein and the subsequent detailed description, or will be described in detail in the accompanying drawings incorporated herein and the subsequent detailed description, and these accompanying drawings and detailed description are used together to explain the specific principles of the present invention. BRIEF DESCRIPTION OF THE DRAWINGS

[0034] By describing the exemplary embodiments of the present invention in more detail in conjunction with the accompanying drawings, the above and other objects, features, and advantages of the present invention will become more apparent, wherein in the exemplary embodiments of the present invention, the same reference numerals generally represent the same components.

[0035] Figure 1 A flowchart showing the steps of a wellbore microseismic event quadrilateral hybrid positioning method according to an embodiment of the present invention is shown.

[0036] Figure 2 A schematic diagram of a wellbore microseismic monitoring model according to an embodiment of the present invention is shown.

[0037] Figure 3 A schematic diagram of the P-wave travel time of a microseismic event picked based on a noise-free model according to an embodiment of the present invention is shown.

[0038] Figure 4 A schematic diagram of the S-wave travel time of a microseismic event picked based on a noise-free model according to an embodiment of the present invention is shown.

[0039] Figure 5 A schematic diagram of the positioning distribution of a microseismic event of the present invention picked based on a noise-free model according to an embodiment of the present invention is shown.

[0040] Figure 6 A schematic diagram of the P-wave travel time of a microseismic event picked based on a noise interference model according to an embodiment of the present invention is shown.

[0041] Figure 7 A schematic diagram of the S-wave travel time of a microseismic event picked based on a noise interference model according to an embodiment of the present invention is shown.

[0042] Figure 8 A schematic diagram of the positioning distribution of a microseismic event of the present invention picked based on a noise interference model according to an embodiment of the present invention is shown. DETAILED DESCRIPTION

[0043] 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 by the embodiments set forth herein.

[0044] The present invention provides a quadrilateral hybrid positioning method for microseismic events in a well, including:

[0045] Step 1: Establish two P-wave nonlinear positioning equations and two S-wave nonlinear positioning equations respectively;

[0046] Step 2: Calculate the initial P-wave travel time positioning result and the initial S-wave travel time positioning result of the microseismic event;

[0047] Step 3: Construct a quadrilateral based on the initial P-wave travel time positioning result and the initial S-wave travel time positioning result, and calculate the center point of the quadrilateral;

[0048] Step 4: Calculate the new initial P-wave travel time positioning result and the new initial S-wave travel time positioning result according to the center point coordinates, and repeat steps 3-4 until the center point coordinates converge to determine the final positioning result.

[0049] In one example, the P-wave nonlinear positioning equation is:

[0050]

[0051]

[0052] where f(P1) is the first P-wave nonlinear positioning equation, f(P2) is the second P-wave nonlinear positioning equation, and ΔTp is the P-wave travel time difference.

[0053] In one example, the S-wave nonlinear positioning equation is:

[0054]

[0055]

[0056] where f(S1) is the first S-wave nonlinear positioning equation, f(S2) is the second S-wave nonlinear positioning equation, and ΔTs is the S-wave travel time difference.

[0057] In one example, by the grid search method, when the P-wave time difference exponential equations are respectively the smallest, the corresponding grid points are the P-wave travel time positioning results.

[0058] In one example, by the grid search method, when the S-wave time difference exponential equations are respectively the smallest, the corresponding grid points are the S-wave travel time positioning results.

[0059] In one example, the initial P-wave travel time positioning result is (xP1 , y P1 ), (x P2 , y P2 ), the initial S-wave travel time location result is (x S1 , y S1 ), (x S2 , y S2 ), the center point coordinates are (x0, y0), then the new initial P-wave travel time location result is The new initial S-wave travel time location result is

[0060] In one example, it is judged whether to converge through the error formula. If the error is less than or equal to 0.1, it converges.

[0061] In one example, the error formula is:

[0062] Δ = |x * 0 - x0| + |y * 0 - y0| (5)

[0063] where, (x0, y0) are the center point coordinates, and (x * 0, y * 0) are the new center point coordinates.

[0064] Specifically, pick up the P-wave travel time T Pi 、S-wave travel time T Si , and define two independent variables:

[0065] The P-wave travel time difference is:

[0066] ΔTp = ∑|T Pi - T P0i | (6)

[0067] The S-wave travel time difference is:

[0068] ΔTs = ∑|T Si - T S0i | (7)

[0069] where, T Pi 、T Si 、T P0i 、T S0i respectively represent the P-wave travel time actually picked up by the i-th detector, the S-wave travel time actually picked up by the i-th detector, the P-wave travel time of the i-th detector's theoretical forward modeling, and the S-wave travel time of the i-th detector's theoretical forward modeling.

[0070] Taking the P-wave travel-time difference ΔTp and the S-wave travel-time difference ΔTs as independent variables respectively, construct their corresponding positioning equations. The two groups of P-wave travel-time non-linear positioning equations are constructed as Formulas (1) and (2) respectively, and the two groups of S-wave travel-time non-linear positioning equations are constructed as Formulas (3) and (4) respectively

[0071] Then, according to the known picked travel times and logging data, solve the positioning equations (1) to (4).

[0072] Theoretically, solve the above positioning equations by calculating the extreme values

[0073]

[0074]

[0075]

[0076]

[0077] In fact, the grid search method is usually used to solve the equations. First step, define the grid (x j , y k ) (j = 1, 2,..., M, k = 1, 2,..., N) range, so that the grid can cover the possible distribution space of the microseismic events as much as possible; Second step, according to the known logging data, establish the P-wave and S-wave velocity models, and use the ray tracing method to calculate the theoretical P-wave travel time T P0i and the S-wave travel time T S0i of each grid point reaching the geophone, and further calculate the travel-time differences of the P-wave, S-wave and PS-wave, and substitute them into the corresponding positioning equations; Third step, calculate and compare the numerical values of the positioning equations of all grid points, and record the grid point corresponding to the minimum value as the corresponding travel-time positioning result. According to this operation, two groups of P-wave travel-time positioning results (x P1 , y P1 ) and (x P2 , y P2 ) can be searched out. Also, two groups of S-wave travel-time positioning results (x S1 , y S1 ) and (x S2 , y S2 ) can be searched out. At the same time, construct a quadrilateral with the four positioning results, and calculate the intersection point of the two center lines of the quadrilateral, that is, the center point position (x0, y0) is

[0078]

[0079] According to the two groups of P-wave positioning results (x P , y P ) and (x S , yS ), The travel time positioning results of two groups of S-waves (x S1 , y S1 ), (x S2 , y S2 ), and the center point position (x0, y0) are used to construct a new quadrilateral Calculate the new center point position (x * 0, y * 0) as follows:

[0080]

[0081] Compare the center point (x0, y0) of the previous step with the new center point (x * 0, y * 0). Define the error function as formula (5). If the error is less than or equal to 0.1, it converges; otherwise, it does not converge.

[0082] When the error does not meet the convergence condition, repeat the above steps to continue constructing a quadrilateral and calculate the center point of the quadrilateral until the error between the current center point and the center point of the previous step meets the convergence condition. Take the position of the current latest quadrilateral center point as the final microseismic event positioning result (x * , y * ).

[0083] The present invention also provides an electronic device, which includes: a memory storing executable instructions; a processor that runs the executable instructions in the memory to implement the above-mentioned hybrid quadrilateral positioning method for in-well microseismic events.

[0084] The present invention also provides a computer-readable storage medium, which stores a computer program. When the computer program is executed by a processor, it implements the above-mentioned hybrid quadrilateral positioning method for in-well microseismic events.

[0085] To facilitate the understanding of the solutions and effects of the embodiments of the present invention, the following gives three specific application examples. Those skilled in the art should understand that this example is only for facilitating the understanding of the present invention, and any specific details are not intended to limit the present invention in any way.

[0086] Example 1

[0087] Figure 1 Shows a flowchart of the steps of the hybrid quadrilateral positioning method for in-well microseismic events according to an embodiment of the present invention.

[0088] As Figure 1As shown in the figure, the quadrilateral hybrid location method for microseismic events in a well includes: Step 1: Establish two P-wave non-linear location equations and two S-wave non-linear location equations respectively; Step 2: Calculate the initial P-wave travel time location result and the initial S-wave travel time location result of the microseismic event; Step 3: Construct a quadrilateral based on the initial P-wave travel time location result and the initial S-wave travel time location result, and calculate the center point of the quadrilateral; Step 4: Calculate the new initial P-wave travel time location result and the new initial S-wave travel time location result according to the center point coordinates, and repeat Steps 3-4 until the center point coordinates converge to determine the final location result.

[0089] Figure 2 The schematic diagram of the downhole microseismic monitoring model according to an embodiment of the present invention is shown, with geophones Δ, events ◆, and stratification interfaces -.

[0090] As Figure 2 shown, 16-level downhole three-component geophones (geophone coordinates are known) are lowered into the monitoring well, 16 event signals, and the spatial coordinates are shown in Table 1, unit: meters, and the formation is divided into three layers.

[0091] Table 1

[0092]

[0093]

[0094] Figure 3 The schematic diagram of the P-wave travel time of the microseismic event picked up based on the noise-free model according to an embodiment of the present invention is shown.

[0095] Figure 4 The schematic diagram of the S-wave travel time of the microseismic event picked up based on the noise-free model according to an embodiment of the present invention is shown.

[0096] Before giving examples of the present invention, using the existing high-precision ray tracing algorithm, according to Figure 2 the observation method, the P-wave travel time ( Figure 3 ) and the S-wave travel time ( Figure 4 ) of each event reaching the geophone are forward modeled as known true observation values.

[0097] First, the picked P-wave and S-wave travel times (such as Figure 4 , Figure 5 ) are used as inputs, the independent variables ΔTp and ΔTs are defined, and two sets of P-wave non-linear location equations and two sets of S-wave non-linear location equations are constructed.

[0098] Solve the positioning equation by the grid search method, that is: define a grid range - 600 to 1600 horizontally, 1200 to 1600 vertically, and the minimum grid unit is 0.1*0.1. According to the known velocity model, use the ray tracing method to calculate the two groups of P-wave nonlinear positioning equations and the two groups of S-wave nonlinear positioning equations for all grid points, and compare and search for the grid point corresponding to the minimum value of the positioning equation, which is the corresponding P-wave positioning result (x P1 , y P1 ), (x P2 , y P2 ), S-wave positioning result (x S1 , y S1 ), (x S2 , y S2 ). At the same time, construct a quadrilateral with these four positioning results, calculate and output the center point position (x0, y0) of the quadrilateral.

[0099] According to the two groups of P-wave positioning results (x P , y P ), (x S , y S ) of the microseismic event calculated, the two groups of S-wave travel time positioning results (x S1 , y S1 ), (x S2 , y S2 ) and the center point position (x0, y0), construct a new quadrilateral Calculate the new center point position (x * 0, y * 0).

[0100] Compare the center point (x0, y0) in the previous step with the new center point (x * 0, y * 0), calculate the error through formula (5), and determine whether it converges. When the error does not converge, continue to construct a quadrilateral with the four points of the new quadrilateral and the new center point, and continue to calculate the center point of the quadrilateral until the error between the current center point and the center point in the previous step meets the convergence condition. Take the position of the center point of the current latest quadrilateral as the final microseismic event positioning result (x * , y * ).

[0101] Figure 5 Fig. shows a schematic diagram of the microseismic event positioning distribution based on the noise-free model according to an embodiment of the present invention, with very small positioning error and very high accuracy.

[0102] Figure 6 Fig. shows a schematic diagram of the P-wave travel time of the microseismic event picked up based on the noise interference model according to an embodiment of the present invention.

[0103] Figure 7 Shows a schematic diagram of the S-wave travel time of microseismic events picked based on a noise interference model according to an embodiment of the present invention.

[0104] Figure 8 Shows a schematic diagram of the location distribution of microseismic events of the present invention based on a noise interference model according to an embodiment of the present invention.

[0105] Given the theoretical P-wave travel time and S-wave travel time (such as Figure 6 , Figure 7 ), adding a certain interference signal to simulate the possible noise interference in actual data. Through the operation process of the present invention, the final microseismic event location result is as shown in Figure 8 . The statistical results of the location error of microseismic events in the noisy model are shown in Table 2, unit: meter.

[0106] Table 2

[0107]

[0108]

[0109] Through the error statistical tables 1 and 2, whether it is a noise-free model or adding a large amount of noise interference, the location result error of the present invention is small and overall controllable. Especially, the lateral error is smaller than the depth error, overcoming the shortcoming of unstable lateral location of traditional borehole microseismic, verifying that the present invention has a high microseismic location accuracy and has certain popularization and application value.

[0110] Example 2

[0111] The present disclosure provides an electronic device, which includes: a memory storing executable instructions; a processor that runs the executable instructions in the memory to implement the above-mentioned quadrilateral hybrid location method for borehole microseismic events.

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

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

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

[0115] Those skilled in the art should understand that, in order to solve the technical problem of how to obtain good user experience effects, known structures such as communication buses and interfaces may also be included in this embodiment, and these known structures should also be included in the protection scope of the present disclosure.

[0116] For a detailed description of this embodiment, reference may be made to the corresponding descriptions in the foregoing embodiments, and details will not be repeated here.

[0117] Example 3

[0118] An embodiment of the present disclosure provides a computer-readable storage medium, which stores a computer program. When the computer program is executed by a processor, the wellbore microseismic event quadrilateral hybrid positioning method described above is implemented.

[0119] According to the computer-readable storage medium of the embodiment of the present disclosure, non-temporary computer-readable instructions are stored thereon. When the non-temporary computer-readable instructions are run by a processor, all or part of the steps of the methods of the foregoing embodiments of the present disclosure are executed.

[0120] The above computer-readable storage media include, but are not limited to: optical storage media (such as CD-ROMs and DVDs), magneto-optical storage media (such as MOs), magnetic storage media (such as tapes or external hard drives), media with built-in rewritable non-volatile memories (such as memory cards), and media with built-in ROMs (such as ROM cartridges).

[0121] 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 example given.

[0122] The above has described the embodiments of the present invention. The above description is exemplary and not exhaustive, and is not limited to the disclosed embodiments. Many modifications and variations are obvious to those of ordinary skill in the art in the technical field without departing from the scope and spirit of the described embodiments.

Claims

1. A quadrilateral hybrid location method for in-well microseismic events, characterized in that Including: Step 1: Establish two P-wave nonlinear positioning equations and two S-wave nonlinear positioning equations respectively; Step 2: Calculate the initial P-wave travel time positioning result and the initial S-wave travel time positioning result of the microseismic event; Step 3: Construct a quadrilateral based on the initial P-wave travel time positioning result and the initial S-wave travel time positioning result, and calculate the center point of the quadrilateral; Step 4: Calculate the new initial P-wave travel time positioning result and the new initial S-wave travel time positioning result according to the center point coordinates, and repeat steps 3 - 4 until the center point coordinates converge to determine the final positioning result; Among them, the initial P-wave travel time positioning results are (x P1 , y P1 ), (x P2 , y P2 ), the initial S-wave travel time positioning results are (x S1 , y S1 ), (x S2 , y S2 ), the center point coordinates are (x0, y0), then the new initial P-wave travel time positioning result is The new initial S-wave travel time positioning result is 2. The quadrilateral hybrid location method for microseismic events in a well according to claim 1, wherein, The P-wave nonlinear positioning equation is: where f(P1) is the first P-wave nonlinear positioning equation, f(P2) is the second P-wave nonlinear positioning equation, and ΔTp is the P-wave travel time difference.

3. The quadrilateral hybrid location method for in-well microseismic events according to claim 1, wherein The S-wave nonlinear positioning equation is: where f(S1) is the first S-wave nonlinear positioning equation, f(S2) is the second S-wave nonlinear positioning equation, and ΔTs is the S-wave travel time difference.

4. The quadrilateral hybrid location method for in-well microseismic events according to claim 1, wherein, Through the grid search method, when the P-wave time difference exponential equations are respectively the smallest, the corresponding grid points are the P-wave travel time positioning results.

5. The method for quadrilateral hybrid positioning of microseismic events in a well according to claim 1, wherein, Through the grid search method, when the S-wave time difference exponential equations are respectively the smallest, the corresponding grid points are the S-wave travel time positioning results.

6. The method for quadrilateral hybrid location of in-well microseismic events according to claim 1, wherein, Judge whether to converge through the error formula. If the error is less than or equal to 0.1, it converges.

7. The method for quadrilateral hybrid positioning of microseismic events in a well according to claim 6, wherein, The error formula is: Δ = |x * 0 - x0| + |y * 0 - y0| (5) Among them, (x0, y0) is the center point coordinate, (x * 0, y * 0) is the new center point coordinate.

8. An electronic device, characterized in that, The electronic device includes: A memory storing executable instructions; A processor that runs the executable instructions in the memory to implement the in-well microseismic event quadrilateral hybrid positioning method according to any one of claims 1 - 7.

9. 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, it implements the in-well microseismic event quadrilateral hybrid positioning method according to any one of claims 1 - 7.

Images