CN116009069A - Quadrilateral hybrid positioning method, electronic equipment and medium for microseismic events in wells - Google Patents

Abstract

The application discloses a quadrilateral mixed positioning method for microseismic events in wells, electronic equipment and media. The method may include: Step 1: Establishing two P-wave nonlinear positioning equations and two S-wave nonlinear positioning equations respectively; Step 2: Calculating the initial P-wave traveltime positioning results and the initial S-wave traveltime positioning results of the microseismic event; Step 3: Construct a quadrilateral according to the initial P-wave travel-time positioning results and the initial S-wave travel-time positioning results, and calculate the center point of the quadrilateral; Step 4: Calculate the new initial P-wave travel-time positioning results and the new initial S-wave travel-time positioning based on the coordinates of the center point As a result, repeat steps 3‑4 until the coordinates of the center point converge to determine the final positioning result. The invention achieves stable and fast positioning of microseismic events by simulating the center of the quadrilateral to iterate continuously.

Description

Quadrilateral hybrid positioning method, electronic equipment and medium for microseismic events in wells

Technical Field

The present invention relates to the field of microseismic signal processing in a well, and more specifically, to a quadrilateral hybrid positioning method, electronic equipment and medium for microseismic events in a well.

Background Art

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.

In addition, the accuracy of microseismic event positioning is also related to the positioning method. The main methods for microseismic positioning in wells are: First, it is based on the forward modeling of P-wave and S-wave event travel time. Representative algorithms include network search method, simulated annealing method, Geiger method, etc. The advantage is that it is easy to implement, but the disadvantage is that the P-wave and S-wave travel time of microseismic events are difficult to accurately pick up, which affects the positioning results; second, it 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 require the first arrival of the event to be picked up, but the disadvantage is that it has high requirements for data signal-to-noise ratio and velocity model, requires a large number of geophones, and has high calculation cost; third, 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.

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. Fracturing, as the most effective measure to increase production in such reservoirs, 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 microseismic positioning results in the well. Therefore, the positioning results of a single microseismic positioning are unstable. Usually, P-wave and S-wave linear or nonlinear positioning are used to constrain each other, so that the positioning process is both stable and can obtain high-precision positioning results.

Therefore, it is necessary to develop a quadrilateral hybrid positioning method, electronic equipment and medium for microseismic events in wells that simulates the continuous iteration of the quadrilateral center.

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

The present invention proposes a method, electronic device and medium for hybrid positioning of microseismic events in a well using quadrilaterals, which can construct a quadrilateral through nonlinear positioning of two groups of P-wave travel times and two groups of S-wave travel times, calculate the position of its center point, and then construct a new quadrilateral using the center point and four positioning results, calculate the position of the new center point, compare the two center points, and if there is no convergence, continue to construct a new quadrilateral and calculate the center point until convergence, and output the final microseismic event positioning result, thereby realizing stable and rapid positioning of microseismic events.

In a first aspect, an embodiment of the present disclosure provides a method for quadrilateral hybrid positioning of microseismic events in a well, comprising:

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 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, repeat steps 3-4 until the center point coordinates converge, and determine the final positioning result.

Preferably, the P-wave nonlinear positioning equation is:

Wherein, f(P ₁ ) is the first P wave nonlinear positioning equation, f(P ₂ ) is the second P wave nonlinear positioning equation, and ΔTp is the P wave travel time difference.

Preferably, the S-wave nonlinear positioning equation is:

Wherein, f(S ₁ ) is the first S-wave nonlinear positioning equation, f(S ₂ ) is the second S-wave nonlinear positioning equation, and ΔTs is the S-wave travel time difference.

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.

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.

新的初始S波走时定位结果为

Preferably, the initial P-wave travel time positioning result is (x _P1 , y _P1 ), (x _P2 , y _P2 ), the initial S-wave travel time positioning result is (x _S1 , y _S1 ), (x _S2 , y _S2 ), and the center point coordinates are (x ₀ , y ₀ ), then the new initial P-wave travel time positioning result is

The new initial S wave travel time positioning result is

Preferably, whether convergence is achieved is determined by an error formula, and convergence is achieved if the error is less than or equal to 0.1.

Preferably, the error formula is:

Δ＝|x ^* ₀ -x ₀ |+|y ^* ₀ -y ₀ | (5)

Among them, (x ₀ ,y ₀ ) are the coordinates of the center point, and (x ^* ₀ ,y ^* ₀ ) are the new coordinates of the center point.

As a specific implementation method of the embodiment of the present disclosure,

In a second aspect, an embodiment of the present disclosure further provides an electronic device, the electronic device comprising:

A memory storing executable instructions;

A processor runs the executable instructions in the memory to implement the quadrilateral hybrid positioning method for microseismic events in a well.

In a third aspect, an embodiment of the present disclosure further provides a computer-readable storage medium storing a computer program, which, when executed by a processor, implements the method for quadrilateral hybrid positioning of microseismic events in a well.

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

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.

FIG1 is a flow chart showing the steps of a method for quadrilateral hybrid positioning of microseismic events in a well according to an embodiment of the present invention.

FIG. 2 shows a schematic diagram of a microseismic monitoring model in a well according to an embodiment of the present invention.

FIG3 is a schematic diagram showing the P-wave travel time of a microseismic event picked up based on a noise-free model according to an embodiment of the present invention.

FIG4 is a schematic diagram showing the S-wave travel time of a microseismic event picked up based on a noise-free model according to an embodiment of the present invention.

FIG5 shows a schematic diagram of the microseismic event location distribution of the present invention based on the noise-free model according to an embodiment of the present invention.

FIG6 is a schematic diagram showing the P-wave travel time of a microseismic event picked up based on a noise interference model according to an embodiment of the present invention.

FIG. 7 is a schematic diagram showing the S-wave travel time of a microseismic event picked up based on a noise interference model according to an embodiment of the present invention.

FIG8 shows a schematic diagram of the microseismic event location distribution of the present invention based on the noise interference model according to an embodiment of the present invention.

DETAILED DESCRIPTION

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.

The present invention provides a method for quadrilateral hybrid positioning of microseismic events in a well, comprising:

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, repeat steps 3-4 until the center point coordinates converge, and determine the final positioning result.

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

Wherein, f(P ₁ ) is the first P wave nonlinear positioning equation, f(P ₂ ) is the second P wave nonlinear positioning equation, and ΔTp is the P wave travel time difference.

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

Wherein, f(S ₁ ) is the first S-wave nonlinear positioning equation, f(S ₂ ) is the second S-wave nonlinear positioning equation, and ΔTs is the S-wave travel time difference.

In one example, 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.

In one example, 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.

新的初始S波走时定位结果为

In an example, 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 ), and the center point coordinates are (x ₀ , y ₀ ). Then the new initial P-wave travel time positioning result is

The new initial S wave travel time positioning result is

In one example, whether convergence is achieved is determined by using an error formula, and if the error is less than or equal to 0.1, convergence is achieved.

In one example, the error formula is:

Δ＝|x ^* ₀ -x ₀ |+|y ^* ₀ -y ₀ | (5)

Among them, (x ₀ ,y ₀ ) are the coordinates of the center point, and (x ^* ₀ ,y ^* ₀ ) are the new coordinates of the center point.

Specifically, the P-wave travel time T _Pi and S-wave travel time T _Si of the microseismic event are picked up, and two independent variables are defined:

The P wave travel time difference is:

ΔTp＝∑|T _Pi -T _P0i | (6)

The S wave travel time difference is:

ΔTs＝∑|T _Si -T _S0i | (7)

Among them, T _Pi , T _Si , T _P0i and T _S0i represent the P-wave travel time actually picked up by the ith detector, the S-wave travel time actually picked up by the ith detector, the P-wave travel time theoretically forward modeled by the ith detector and the S-wave travel time theoretically forward modeled by the ith detector, respectively.

The P-wave travel time difference ΔTp and the S-wave travel time difference ΔTs are taken as independent variables to construct their corresponding positioning equations. The two sets of P-wave travel time nonlinear positioning equations are formulas (1) and (2), respectively. The two sets of S-wave travel time nonlinear positioning equations are formulas (3) and (4), respectively.

Then, based on the known picked travel time and logging data, the positioning equations (1) to (4) are solved.

In theory, the above positioning equation is solved by calculating the extreme value:

In fact, the grid search method is usually used to solve the equation. The first step is to define the range of the grid (x _j , y _k ) (j = 1, 2, ..., M, k = 1, 2, ..., N) so that the grid covers the possible distribution space of microseismic events as much as possible; the second step is to establish the P-wave and S-wave velocity models based on the known logging data, and use the ray tracing method to calculate the theoretical P-wave travel time T _P0i and S-wave travel time T _S0i of each grid point reaching the detector, and further calculate the travel time difference of the P-wave, S-wave and PS-wave, and substitute them into the corresponding positioning equation; the third step is to calculate and compare the values of the positioning equation 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, and two groups of S-wave travel time positioning results (x _S1 , y _S1 ) and (x _S2 , y _S2 ) can also be searched out. At the same time, the four positioning results are used to construct a quadrilateral, and the intersection point of the two center lines of the quadrilateral, that is, the center point position (x ₀ ,y ₀ ) is calculated as:

计算新的中心点位置(x^* ₀,y^* ₀)为：Based on the two sets of P-wave location results (x _P ,y _P ) and (x _S ,y _S ) of the microseismic events, the two sets of S-wave travel time location results (x _S1 ,y _S1 ) and (x _S2 ,y _S2 ) and the center point position (x ₀ ,y ₀ ), a new quadrilateral is constructed.

Calculate the new center point position (x ^* ₀ , y ^* ₀ ) as:

Compare the center point of the previous step (x ₀ , y ₀ ) with the new center point (x ^* ₀ , y ^* ₀ ), and 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.

When the error does not meet the convergence condition, repeat the above steps, continue to build a quadrilateral, and continue to 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. The current latest quadrilateral center point position is used as the final microseismic event location result (x ^* , y ^* ).

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

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 quadrilateral hybrid positioning method for microseismic events in a well is implemented.

To facilitate understanding of the solutions and effects of the embodiments of the present invention, three 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.

Example 1

FIG1 is a flow chart showing the steps of a method for quadrilateral hybrid positioning of microseismic events in a well according to an embodiment of the present invention.

As shown in Figure 1, the quadrilateral hybrid positioning method for microseismic events in the well includes: step 1: establishing two P-wave nonlinear positioning equations and two S-wave nonlinear positioning equations respectively; step 2: calculating the initial P-wave travel time positioning result and the initial S-wave travel time positioning result of the microseismic event; step 3: constructing a quadrilateral according to the initial P-wave travel time positioning result and the initial S-wave travel time positioning result, and calculating the center point of the quadrilateral; step 4: calculating the new initial P-wave travel time positioning result and the new initial S-wave travel time positioning result according to the coordinates of the center point, repeating steps 3-4 until the coordinates of the center point converge, and determining the final positioning result.

FIG2 shows a schematic diagram of a microseismic monitoring model in a well according to an embodiment of the present invention, including a detector Δ, an event ◆, and a layered interface -.

As shown in Figure 2, 16 downhole three-component geophones are placed in the monitoring well (the coordinates of the geophones are known), and the spatial coordinates of the 16 event signals are shown in Table 1 (unit: meter). The stratum is divided into three layers.

Table 1

FIG3 is a schematic diagram showing the P-wave travel time of a microseismic event picked up based on a noise-free model according to an embodiment of the present invention.

FIG4 is a schematic diagram showing the S-wave travel time of a microseismic event picked up based on a noise-free model according to an embodiment of the present invention.

Before the present invention is given as an example, the existing high-precision ray tracing algorithm is used to perform the P-wave travel time (FIG. 3) and S-wave travel time (FIG. 4) of each event arriving at the detector according to the observation method of FIG2 as known true observation values.

First, the picked P-wave and S-wave travel times (as shown in Figures 4 and 5) are used as input, and the independent variables ΔTp and ΔTs are defined to construct two sets of P-wave nonlinear positioning equations and two sets of S-wave nonlinear positioning equations.

The positioning equation is solved by the grid search method, that is, a grid range is defined - 600-1600 horizontally, 1200-1600 vertically and the minimum grid unit is 0.1*0.1. According to the known velocity model, the ray tracing method is used to calculate two sets of P-wave nonlinear positioning equations and two sets of S-wave nonlinear positioning equations for all grid points, and the grid points corresponding to the minimum positioning equation values are compared and searched, which are the corresponding P-wave positioning results (x _P1 ,y _P1 ), (x _P2 ,y _P2 ), S-wave positioning results (x _S1 ,y _S1 ), (x _S2 ,y _S2 ). At the same time, these four positioning results are used to construct a quadrilateral, and the center point position (x ₀ ,y ₀ ) of the quadrilateral is calculated and output.

计算新的中心点位置(x^* ₀,y^* ₀)。Based on the two sets of P-wave location results (x _P ,y _P ) and (x _S ,y _S ) of the microseismic events, the two sets of S-wave travel time location results (x _S1 ,y _S1 ) and (x _S2 ,y _S2 ) and the center point position (x ₀ ,y ₀ ), a new quadrilateral is constructed.

Calculate the new center point position (x ^* ₀ ,y ^* ₀ ).

Compare the center point of the previous step (x ₀ , y ₀ ) with the new center point (x ^* ₀ , y ^* ₀ ), calculate the error using formula (5), and determine whether it converges. If the error does not converge, continue to construct a quadrilateral by adding the four points of the new quadrilateral to 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 of the previous step meets the convergence condition. The current latest quadrilateral center point position is used as the final microseismic event location result (x ^* , y ^* ).

FIG5 shows a schematic diagram of the microseismic event location distribution of the present invention based on the noise-free model according to an embodiment of the present invention, in which the location error is very small and the accuracy is very high.

FIG6 is a schematic diagram showing the P-wave travel time of a microseismic event picked up based on a noise interference model according to an embodiment of the present invention.

FIG. 7 is a schematic diagram showing the S-wave travel time of a microseismic event picked up based on a noise interference model according to an embodiment of the present invention.

FIG8 shows a schematic diagram of the microseismic event location distribution of the present invention based on the noise interference model according to an embodiment of the present invention.

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

Table 2

Through the error statistics Table 1 and Table 2, no matter it is a noise-free model or a model with large noise interference, the positioning result error of the present invention is small and controllable as a whole, especially the lateral error is smaller than the depth error, which overcomes the disadvantage of unstable lateral positioning of traditional microseismic in-wells, verifies that the present invention has high microseismic positioning accuracy and has certain promotion and application value.

Example 2

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 quadrilateral hybrid positioning method for microseismic events in a well.

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

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.

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.

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.

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

Example 3

An embodiment of the present disclosure provides a computer-readable storage medium storing a computer program. When the computer program is executed by a processor, the method for quadrilateral hybrid positioning of microseismic events in a well is implemented.

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.

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).

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.

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 (10)

1. A method for quadrilateral hybrid positioning of microseism events in a well, comprising the steps of:

step 1: respectively establishing two P-wave nonlinear positioning equations and two S-wave nonlinear positioning equations;

step 2: calculating an initial P wave travel time positioning result and an initial S wave travel time positioning result of the microseism event;

step 3: constructing a quadrangle according to the initial P wave travel time positioning result and the initial S wave travel time positioning result, and calculating a center point of the quadrangle;

step 4: and (3) calculating a new initial P wave travel time positioning result and a new initial S wave travel time positioning result according to the center point coordinates, repeating the steps 3-4 until the center point coordinates are converged, and determining a final positioning result.

2. The method of hybrid positioning of a quadrilateral of a microseismic event in a well of claim 1 wherein the P-wave nonlinear positioning equation is:

wherein f (P ₁ ) For the first P-wave nonlinear positioning equation, f (P ₂ ) For the second P-wave nonlinear positioning equation, ΔTp is the P-wave travel time difference.

3. The method of hybrid positioning of a quadrilateral of a microseismic event in a well of claim 1 wherein the S-wave nonlinear positioning equation is:

wherein f (S) ₁ ) For the first S-wave nonlinear positioning equation, f (S ₂ ) And as a second S-wave nonlinear positioning equation, deltaTs is the S-wave travel time difference.

4. The method for hybrid positioning of a well microseism event quadrangle according to claim 1, wherein when the P-wave time difference index equations are respectively minimum, the corresponding grid points are the P-wave travel time positioning results by a grid search method.

5. The method for hybrid positioning of a well microseism event quadrangle according to claim 1, wherein when the S-wave time difference index equations are respectively minimum, the corresponding grid points are the S-wave travel time positioning results by a grid search method.

6. The method of hybrid positioning of a quadrilateral of a microseismic event in a well of claim 1 wherein the initial P-wave travel time positioning result is (x _P1 ,y _P1 )、(x _P2 ,y _P2 ) The initial S-wave travel time positioning result is (x) _S1 ,y _S1 )、(x _S2 ,y _S2 ) The center point coordinates are (x ₀ ,y ₀ ) The new initial P-wave travel time positioning result is

The new initial S-wave travel time positioning result is

7. The method for hybrid positioning of a quadrilateral of a microseismic event in a well according to claim 1, wherein whether the quadrilateral of the microseismic event is converged is judged by an error formula, and if the error is less than or equal to 0.1, the quadrilateral of the microseismic event is converged.

8. The method of quadrilateral hybrid positioning of microseismic events in a well of claim 1 wherein the error formula is:

Δ＝|x ^* ₀ -x ₀ |+|y ^* ₀ -y ₀ | (5)

wherein, (x) ₀ ,y ₀ ) Is the center point coordinate, (x) ^* ₀ ,y ^* ₀ ) Is the new center point coordinates.

9. An electronic device, the electronic device comprising:

a memory storing executable instructions;

a processor executing the executable instructions in the memory to implement the method of well microseismic event quadrilateral hybrid positioning of any of claims 1-8.

10. A computer readable storage medium, characterized in that it stores a computer program which, when executed by a processor, implements the method for quadrilateral hybrid localization of microseismic events in wells according to any one of claims 1 to 8.

Images