Hybrid microseismic event location method and device in well, electronic equipment and medium

Abstract

The application discloses a well microseismic event mixed positioning method and device, electronic equipment and medium. The method can comprise: calculating P-wave travel time time difference, S-wave travel time time difference and PS-wave travel time time difference; constructing P-wave time difference exponential equation, S-wave time difference exponential equation and PS-wave time difference exponential equation respectively; calculating P-wave travel time positioning result, S-wave travel time positioning result and PS-wave time difference positioning result of the microseismic event respectively; and determining final microseismic event spatial positioning result according to the P-wave, S-wave and PS-wave time difference positioning result. The application constructs a nonlinear exponential positioning equation with stable calculation process through P-wave travel time, S-wave travel time and PS-wave time difference, simulates a spatial triangle through the three positioning results, and then calculates the perpendicular center position of the triangle, that is, the final microseismic event positioning result, so that stable and rapid positioning of the microseismic event is realized.

Description

Technical Field

[0001] This invention relates to the field of borehole microseismic signal processing, and more specifically, to a hybrid localization method, apparatus, electronic device, and medium for borehole microseismic events. Background Technology

[0002] Downhole microseismic monitoring is one of the methods of microseismic observation. Its characteristic is that downhole three-component geophones receive the full wavefield signal of microseismic events. Compared to surface microseismic monitoring, downhole data has a higher signal-to-noise ratio and a richer variety of microseismic events and types. However, due to the limited number of downhole microseismic geophones (generally 12–32 magnitude three-component geophones), the monitoring range is relatively small. Furthermore, there are certain requirements for the selection of monitoring and observation wells. For example, internationally, the optimal radial distance from the fracturing section of the monitoring well to the geophone in the observation well is 200–800 meters, but the actual distance is often much larger. These factors typically lead to instability and low accuracy in microseismic location.

[0003] Furthermore, the accuracy of microseismic event location is also related to the location method. The main microseismic location methods in wells include: First, forward modeling based on P-wave and S-wave event travel times, with representative algorithms such as network search, simulated annealing, and Geiger's method. The advantage is ease of implementation, but the disadvantage is the difficulty in accurately picking up the P-wave and S-wave travel times of microseismic events, affecting the location results. Second, wave equation convolution, with representative algorithms such as interferometry, reverse time migration, and passive source imaging. The advantage is that it does not require picking up the first arrival of events, but the disadvantage is high requirements for data signal-to-noise ratio, velocity model, and a large number of detectors, resulting in high computational costs. Third, the difference between anisotropic and isotropic travel time calculations. In anisotropic media, isotropic travel time calculations result in larger errors, leading to correspondingly larger location errors.

[0004] In recent years, with the continuous increase in global energy demand and the advancement of exploration technology, oil and gas resource exploration and development have continued to deepen. Fracturing, as the most effective production enhancement measure for such reservoirs, has also received increasing attention both domestically and internationally. In-well microseismic monitoring is often a single-well monitoring method, and due to the limited number of geophones, the positioning angle is small, leading to unstable in-well microseismic positioning results.

[0005] Therefore, it is necessary to develop a hybrid method, device, electronic equipment, and medium for locating microseismic events in wells.

[0006] The information disclosed in the background section of this invention is intended only to enhance the understanding of the general background of this invention, and should not be construed as an admission or in any way implying that such information constitutes prior art known to those skilled in the art. Summary of the Invention

[0007] This invention proposes a hybrid location method, device, electronic equipment, and medium for microseismic events in wells. It can construct a relatively stable nonlinear exponential location equation by using P-wave travel time, S-wave travel time, and PS-wave time difference. The three location results obtained are used to simulate a spatial triangle, and the orthocenter of the triangle is calculated to obtain the final location result of the microseismic event, thereby achieving stable and rapid location of microseismic events.

[0008] In a first aspect, embodiments of this disclosure provide a hybrid localization method for microseismic events in wells, including:

[0009] Calculate the travel time difference of P-wave, S-wave, and PS-wave;

[0010] Construct the P-wave time difference exponential equation, the S-wave time difference exponential equation, and the PS-wave time difference exponential equation respectively;

[0011] Calculate the P-wave travel time location results, S-wave travel time location results, and PS-wave time difference location results for microseismic events, respectively.

[0012] Based on the time difference positioning results of P-wave, S-wave, and PS-wave, the final spatial positioning results of the microseismic event were determined.

[0013] Preferably, the P-wave time difference exponential equation is:

[0014]

[0015] Where ΔTp is the travel time difference of the P-wave.

[0016] Preferably, the S-wave time difference exponential equation is:

[0017]

[0018] Where ΔTs is the travel time difference of the S-wave.

[0019] Preferably, the PS wave time difference exponential equation is:

[0020]

[0021] Where ΔTps is the travel time difference of the PS wave.

[0022] Preferably, by using the grid search method, when the P-wave time difference exponential equation is minimized, the corresponding grid point is the P-wave travel time positioning result;

[0023] When the S-wave time difference exponential equation is minimized, the corresponding grid point is the S-wave travel time positioning result;

[0024] When the PS wave time difference exponential equation is minimized, the corresponding grid point is the PS wave travel time positioning result.

[0025] Preferably, a spatial triangle is established based on the time difference positioning results of P-wave, S-wave, and PS-wave, and the orthocenter position of the spatial triangle is calculated, which is the final spatial positioning result of the microseismic event.

[0026] Preferably, the orthocenter position of the spatial triangle is calculated using formula (4):

[0027]

[0028] Among them, (x P ,y P (x) represents the P-wave time difference localization result. S ,y S (x) represents the S-wave time difference positioning result. PS ,y PS (x) represents the PS wave time difference localization result. * ,y * () represents the final spatial location result of the microseismic event.

[0029] As one specific implementation of this disclosure,

[0030] Secondly, embodiments of this disclosure also provide a hybrid location device for microseismic events in wells, comprising:

[0031] The time difference calculation module calculates the time difference of P-wave, S-wave, and PS-wave.

[0032] The module is used to construct the P-wave time difference index equation, the S-wave time difference index equation, and the PS-wave time difference index equation, respectively.

[0033] The location result calculation module calculates the P-wave travel time location results, S-wave travel time location results, and PS-wave time difference location results for microseismic events, respectively.

[0034] The positioning module determines the final spatial positioning result of the microseismic event based on the time difference positioning results of P-wave, S-wave, and PS-wave.

[0035] Preferably, the P-wave time difference exponential equation is:

[0036]

[0037] Where ΔTp is the travel time difference of the P-wave.

[0038] Preferably, the S-wave time difference exponential equation is:

[0039]

[0040] Where ΔTs is the travel time difference of the S-wave.

[0041] Preferably, the PS wave time difference exponential equation is:

[0042]

[0043] Where ΔTps is the travel time difference of the PS wave.

[0044] Preferably, by using the grid search method, when the P-wave time difference exponential equation is minimized, the corresponding grid point is the P-wave travel time positioning result;

[0045] When the S-wave time difference exponential equation is minimized, the corresponding grid point is the S-wave travel time positioning result;

[0046] When the PS wave time difference exponential equation is minimized, the corresponding grid point is the PS wave travel time positioning result.

[0047] Preferably, a spatial triangle is established based on the time difference positioning results of P-wave, S-wave, and PS-wave, and the orthocenter position of the spatial triangle is calculated, which is the final spatial positioning result of the microseismic event.

[0048] Preferably, the orthocenter position of the spatial triangle is calculated using formula (4):

[0049]

[0050] Among them, (x P ,y P (x) represents the P-wave time difference localization result. S ,y S (x) represents the S-wave time difference positioning result. PS ,y PS (x) represents the PS wave time difference localization result. * ,y * () represents the final spatial location result of the microseismic event.

[0051] Thirdly, embodiments of this disclosure also provide an electronic device, the electronic device comprising:

[0052] Memory, which stores executable instructions;

[0053] A processor that executes the executable instructions in the memory to implement the hybrid location method for microseismic events in wells.

[0054] Fourthly, embodiments of this disclosure also provide a computer-readable storage medium storing a computer program that, when executed by a processor, implements the hybrid location method for microseismic events in wells.

[0055] The methods and apparatus of the present invention have other features and advantages that will be apparent from or will be set forth in detail in the accompanying drawings and following detailed description, which together serve to explain the particular principles of the invention. Attached Figure Description

[0056] The above and other objects, features and advantages of the present invention will become more apparent from the more detailed description of exemplary embodiments of the invention in conjunction with the accompanying drawings, wherein the same reference numerals generally represent the same parts.

[0057] Figure 1 A flowchart illustrating the steps of a hybrid location method for microseismic events in a well according to an embodiment of the present invention is shown.

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

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

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

[0061] Figure 5 A schematic diagram of the location distribution of microseismic events based on a noise-free model according to an embodiment of the present invention is shown.

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

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

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

[0065] Figure 9 A block diagram of a hybrid well microseismic event location device according to an embodiment of the present invention is shown.

[0066] Explanation of reference numerals in the attached figures:

[0067] 201. Travel time difference calculation module; 202. Construction module; 203. Positioning result calculation module; 204. Positioning module. Detailed Implementation

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

[0069] This invention provides a hybrid method for locating microseismic events in a well, comprising:

[0070] Calculate the travel time difference of P-wave, S-wave, and PS-wave;

[0071] Construct the P-wave time difference exponential equation, the S-wave time difference exponential equation, and the PS-wave time difference exponential equation respectively;

[0072] Calculate the P-wave travel time location results, S-wave travel time location results, and PS-wave time difference location results for microseismic events, respectively.

[0073] Based on the time difference positioning results of P-wave, S-wave, and PS-wave, the final spatial positioning results of the microseismic event were determined.

[0074] In one example, the P-wave time difference exponential equation is:

[0075]

[0076] Where ΔTp is the travel time difference of the P-wave.

[0077] In one example, the S-wave time difference exponential equation is:

[0078]

[0079] Where ΔTs is the travel time difference of the S-wave.

[0080] In one example, the PS wave time difference exponential equation is:

[0081]

[0082] Where ΔTps is the travel time difference of the PS wave.

[0083] In one example, using the grid search method, when the P-wave time difference exponential equation is minimized, the corresponding grid point is the P-wave travel time positioning result;

[0084] When the S-wave time difference exponential equation is minimized, the corresponding grid point is the S-wave travel time positioning result;

[0085] When the PS wave time difference exponential equation is minimized, the corresponding grid point is the PS wave travel time positioning result.

[0086] In one example, a spatial triangle is established based on the time difference positioning results of P-wave, S-wave, and PS-wave, and the orthocenter position of the spatial triangle is calculated, which is the final spatial positioning result of the microseismic event.

[0087] In one example, the orthocenter of the spatial triangle is calculated using formula (4):

[0088]

[0089] Among them, (x P ,y P (x) represents the P-wave time difference localization result. S ,y S (x) represents the S-wave time difference positioning result. PS ,y PS (x) represents the PS wave time difference localization result. * ,y * () represents the final spatial location result of the microseismic event.

[0090] Specifically, the relevant location index equations for P-waves and S-waves are constructed: the travel time T of the P-wave of the microseismic event is picked. Pi S-wave travel time T Si Define three independent variables:

[0091] The travel time difference of the P wave is:

[0092] ΔTp=∑|T Pi -T P0i | (5)

[0093] S-wave travel time difference is

[0094] ΔTs=∑|T Si -T S0i | (6)

[0095] PS wave travel time difference is

[0096] ΔTps=∑|(T Pi -T Si )-(T P0i -T S0i (7)

[0097] Among them, T Pi T Si T P0i T S0i Let represent the actual P-wave travel time picked up by the i-th detector, the actual S-wave travel time picked up by the i-th detector, the theoretical forward modeling P-wave travel time of the i-th detector, and the theoretical forward modeling S-wave travel time of the i-th detector, respectively.

[0098] Using the P-wave travel time difference ΔTp, S-wave travel time difference ΔTs, and PS-wave travel time difference ΔTps as independent variables, respectively, construct their corresponding positioning index equations:

[0099] P-wave travel time difference exponential equation f(e) ΔTp The equation for the travel time difference of the S-wave is given by formula (1), f(e) ΔTs ) is Equation (2), the PS wave time difference exponential equation f(e) ΔTps ) is formula (3).

[0100] Then, based on the known pick-up travel time and well logging data, solve the location equations, i.e., formulas (1)-(3).

[0101] Theoretically, the above exponential equation can be solved by calculating the extreme values:

[0102]

[0103]

[0104]

[0105] In practice, a grid search method is usually used. The first step is to define the grid (x... j ,y k The range (j = 1, 2, ..., M, k = 1, 2, ..., N) is used to ensure that the grid covers the possible spatial distribution of microseismic events as much as possible.

[0106] The second step is to establish P-wave and S-wave velocity models based on the known logging data, and then use the ray tracing method to calculate the theoretical P-wave travel time T at each grid point to reach the detector. P0i S-wave travel time T S0i Combined with the actual P-wave travel time T Pi S-wave travel time T Si Calculate the travel time difference of P-wave, S-wave and PS-wave, and substitute it into the corresponding positioning index equation.

[0107] The third step is to calculate and compare the exponential equation f(e) for all grid points. ΔTp ), and the exponential equation f(e ΔTp The grid point corresponding to the minimum value is denoted as (x). P ,y P This is the P-wave travel time location result. Using the same operation, the exponential equation f(e) can be searched for. ΔTs The grid point corresponding to the minimum value is the S-wave travel time localization result (x). S ,y S ), and can also search for the exponential equation f(e ΔTps The grid point corresponding to the minimum value is the PS wave time difference localization result (x).PS ,y PS ).

[0108] Finally, based on the calculated three location results of the microseismic events (x P ,y P ), (x S ,y S ), (x PS ,y PS Construct a triangle. Using the property that the product of the slopes of two perpendicular lines is -1, establish the orthocenter (x0). * ,y * The system of equations, namely formula (4), yields the final microseismic event location result (x). * ,y * This enables the final positioning process.

[0109] The present invention also provides a hybrid location device for microseismic events in wells, comprising:

[0110] The time difference calculation module calculates the time difference of P-wave, S-wave, and PS-wave.

[0111] The module is used to construct the P-wave time difference index equation, the S-wave time difference index equation, and the PS-wave time difference index equation, respectively.

[0112] The location result calculation module calculates the P-wave travel time location results, S-wave travel time location results, and PS-wave time difference location results for microseismic events, respectively.

[0113] The positioning module determines the final spatial positioning result of the microseismic event based on the time difference positioning results of P-wave, S-wave, and PS-wave.

[0114] In one example, the P-wave time difference exponential equation is:

[0115]

[0116] Where ΔTp is the travel time difference of the P-wave.

[0117] In one example, the S-wave time difference exponential equation is:

[0118]

[0119] Where ΔTs is the travel time difference of the S-wave.

[0120] In one example, the PS wave time difference exponential equation is:

[0121]

[0122] Where ΔTps is the travel time difference of the PS wave.

[0123] In one example, using the grid search method, when the P-wave time difference exponential equation is minimized, the corresponding grid point is the P-wave travel time positioning result;

[0124] When the S-wave time difference exponential equation is minimized, the corresponding grid point is the S-wave travel time positioning result;

[0125] When the PS wave time difference exponential equation is minimized, the corresponding grid point is the PS wave travel time positioning result.

[0126] In one example, a spatial triangle is established based on the time difference positioning results of P-wave, S-wave, and PS-wave, and the orthocenter position of the spatial triangle is calculated, which is the final spatial positioning result of the microseismic event.

[0127] In one example, the orthocenter of the spatial triangle is calculated using formula (4):

[0128]

[0129] Among them, (x P ,y P (x) represents the P-wave time difference localization result. S ,y S (x) represents the S-wave time difference positioning result. PS ,y PS (x) represents the PS wave time difference localization result. * ,y * () represents the final spatial location result of the microseismic event.

[0130] Specifically, the relevant location index equations for P-waves and S-waves are constructed: the travel time T of the P-wave of the microseismic event is picked. Pi S-wave travel time T Si Define three independent variables:

[0131] The travel time difference of the P wave is:

[0132] ΔTp=∑|T Pi -T P0i | (5)

[0133] S-wave travel time difference is

[0134] ΔTs=∑|T Si -T S0i | (6)

[0135] PS wave travel time difference is

[0136] ΔTps=∑|(T Pi -T Si )-(T P0i -T S0i (7)

[0137] Among them, T Pi T Si T P0i T S0i Let represent the actual P-wave travel time picked up by the i-th detector, the actual S-wave travel time picked up by the i-th detector, the theoretical forward modeling P-wave travel time of the i-th detector, and the theoretical forward modeling S-wave travel time of the i-th detector, respectively.

[0138] Using the P-wave travel time difference ΔTp, S-wave travel time difference ΔTs, and PS-wave travel time difference ΔTps as independent variables, respectively, construct their corresponding positioning index equations:

[0139] P-wave travel time difference exponential equation f(e) ΔTp The equation for the travel time difference of the S-wave is given by formula (1), f(e) ΔTs ) is Equation (2), the PS wave time difference exponential equation f(e) ΔTps ) is formula (3).

[0140] Then, based on the known pick-up travel time and well logging data, solve the location equations, i.e., formulas (1)-(3).

[0141] Theoretically, the above exponential equation can be solved by calculating the extreme values:

[0142]

[0143]

[0144]

[0145] In practice, a grid search method is usually used. The first step is to define the grid (x... j ,y k The range (j = 1, 2, ..., M, k = 1, 2, ..., N) is used to ensure that the grid covers the possible spatial distribution of microseismic events as much as possible.

[0146] The second step is to establish P-wave and S-wave velocity models based on the known logging data, and then use the ray tracing method to calculate the theoretical P-wave travel time T at each grid point to reach the detector. P0i S-wave travel time T S0i Combined with the actual P-wave travel time T Pi S-wave travel time T Si Calculate the travel time difference of P-wave, S-wave and PS-wave, and substitute it into the corresponding positioning index equation.

[0147] The third step is to calculate and compare the exponential equation f(e) for all grid points. ΔTp ), and the exponential equation f(e ΔTp The grid point corresponding to the minimum value is denoted as (x). P,y P This is the P-wave travel time location result. Using the same operation, the exponential equation f(e) can be searched for. ΔTs The grid point corresponding to the minimum value is the S-wave travel time localization result (x). S ,y S ), and can also search for the exponential equation f(e ΔTps The grid point corresponding to the minimum value is the PS wave time difference localization result (x). PS ,y PS ).

[0148] Finally, based on the calculated three location results of the microseismic events (x P ,y P ), (x S ,y S ), (x PS ,y PS Construct a triangle. Using the property that the product of the slopes of two perpendicular lines is -1, establish the orthocenter (x0). * ,y * The system of equations, namely formula (4), yields the final microseismic event location result (x). * ,y * This enables the final positioning process.

[0149] The present invention also provides an electronic device, comprising: a memory storing executable instructions; and a processor executing the executable instructions in the memory to implement the above-described hybrid location method for microseismic events in wells.

[0150] The present invention also provides a computer-readable storage medium storing a computer program that, when executed by a processor, implements the above-described hybrid location method for microseismic events in wells.

[0151] 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 these examples are merely for the purpose of understanding the present invention, and any specific details therein are not intended to limit the present invention in any way.

[0152] Example 1

[0153] Figure 1 A flowchart illustrating the steps of a hybrid location method for microseismic events in a well according to an embodiment of the present invention is shown.

[0154] like Figure 1As shown, the hybrid location method for microseismic events in the well includes: Step 101, calculating the travel time difference of P-waves, S-waves, and PS-waves; Step 102, constructing the exponential equations for P-wave, S-wave, and PS-wave travel time differences, respectively; Step 103, calculating the location results for P-waves, S-waves, and PS-waves based on the location results for the microseismic events; Step 104, determining the final spatial location result for the microseismic events based on the location results for P-waves, S-waves, and PS-waves.

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

[0156] like Figure 2 As shown, a 16-level downhole three-component geophone was deployed in the monitoring well (the geophone coordinates are known), and 16 event signals were generated. The spatial coordinates are shown in Table 1, in meters. The formation is divided into three layers.

[0157] Table 1

[0158] Event X Real Value The true value of event Y Event X Error Event Y Error 800 1400 1.8 -8.5 850 1405 0.6 -1.1 900 1410 1.5 -6 950 1415 1 -6.5 1000 1420 0.8 -7.6 1050 1425 -0.1 0.6 1100 1430 -0.6 2.5 1150 1435 -0.1 1.3 1200 1440 0 0.5 1250 1445 0.1 1.1 1300 1450 -0.3 1 950 1330 0.5 -4.5 820 1450 0.6 -1.8 1250 1350 0.3 -3.3 1230 1500 -1.1 8.1 1030 1380 0 0

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

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

[0161] Using existing high-precision ray tracing algorithms, based on Figure 2 The observation method is to simulate the P-wave travel time of each event arriving at the detector, as shown in the figure. Figure 3 As shown, the S-wave travel time is as follows Figure 4 As shown, these are known actual observations.

[0162] First, the P-wave and S-wave travel times are picked up as inputs, and independent variables are defined (formulas (5) to (7)). The P-wave positioning index equation, S-wave positioning index equation, and PS-wave positioning index equation (formulas (1) to (3)) are constructed.

[0163] Then, the three positioning equations (formulas (8) to (10)) are solved by the grid search method. That is, a grid range is defined—600 to 1600 horizontally and 1200 to 1600 vertically with a minimum grid unit of 0.1*0.1. Based on the known velocity model, the P-wave positioning index equation, S-wave positioning index equation, and PS-wave positioning index equation (formulas (1) to (3)) of all grid points are calculated using the ray tracing method. The grid points corresponding to the minimum positioning index equation are compared and searched, which are the corresponding P-wave positioning results (x). P,y P S-wave localization results (x) S ,y S PS wave localization results (x) PS ,y PS ).

[0164] Figure 5 A schematic diagram of the location distribution of microseismic events based on a noise-free model according to an embodiment of the present invention is shown.

[0165] Finally, construct a triangle from the three positioning results, calculate and output the orthocenter (x) of the triangle. * ,y * (Formula (4)) is the final spatial location result of the microseismic event. Figure 5 This invention provides a spatial distribution for a noise-free model.

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

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

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

[0169] However, actual data may contain noise interference, and the positioning effect of this invention needs to be tested under noisy conditions. Input P-wave travel time and S-wave travel time (e.g., with significant noise interference) Figure 6 , Figure 7 Following the same procedure, the final microseismic event location results are as follows: Figure 8 As shown.

[0170] The statistics of microseismic event location errors in the noisy model are shown in Table 2, in meters.

[0171] Table 2

[0172] Event X Real Value The true value of event Y Event X Error Event Y Error 800 1400 4.6 -31.5 850 1405 0.1 -30.1 900 1410 -3 10.1 950 1415 5.6 -35.8 1000 1420 -1.1 4.3 1050 1425 1.6 -8.9 1100 1430 -4.6 27.3 1150 1435 1.5 -6.3 1200 1440 2.3 -10.3 1250 1445 -0.5 1.75 1300 1450 2 -3 950 1330 -1.5 3.5 820 1450 3.1 -13.5 1250 1350 -1.3 3.3 1230 1500 -3.5 15 1030 1380 -1.6 -1.1

[0173] According to the error statistics in Tables 1 and 2, the overall error of the positioning results of the present invention is controllable, which verifies that the present invention has high microseismic positioning accuracy and has certain promotion and application value.

[0174] Example 2

[0175] Figure 9A block diagram of a hybrid well microseismic event location device according to an embodiment of the present invention is shown.

[0176] like Figure 9 As shown, the well-drilled microseismic event hybrid location device includes:

[0177] The time difference calculation module 201 calculates the time difference of P-wave, S-wave, and PS-wave.

[0178] Module 202 is used to construct the P-wave time difference index equation, the S-wave time difference index equation, and the PS-wave time difference index equation, respectively.

[0179] The positioning result calculation module 203 calculates the P-wave travel time positioning result, S-wave travel time positioning result, and PS-wave time difference positioning result of the microseismic event, respectively.

[0180] The positioning module 204 determines the final spatial positioning result of the microseismic event based on the positioning results of the time difference of P-wave, S-wave, and PS-wave.

[0181] As an alternative, the equation for the P-wave time difference index is:

[0182]

[0183] Where ΔTp is the travel time difference of the P-wave.

[0184] As an alternative, the S-wave time difference exponential equation is:

[0185]

[0186] Where ΔTs is the travel time difference of the S-wave.

[0187] As an alternative, the equation for the PS wave time difference index is:

[0188]

[0189] Where ΔTps is the travel time difference of the PS wave.

[0190] As an alternative, the grid search method can be used to find the P-wave travel time location result when the P-wave time difference exponential equation is minimized.

[0191] When the S-wave time difference exponential equation is minimized, the corresponding grid point is the S-wave travel time positioning result;

[0192] When the PS wave time difference exponential equation is minimized, the corresponding grid point is the PS wave travel time positioning result.

[0193] As an optional approach, a spatial triangle is established based on the time difference positioning results of P-wave, S-wave, and PS-wave, and the orthocenter position of the spatial triangle is calculated, which is the final spatial positioning result of the microseismic event.

[0194] As an alternative, the position of the orthocenter of the spatial triangle can be calculated using formula (4):

[0195]

[0196] Among them, (x P ,y P (x) represents the P-wave time difference localization result. S ,y S (x) represents the S-wave time difference positioning result. PS ,y PS (x) represents the PS wave time difference localization result. * ,y * () represents the final spatial location result of the microseismic event.

[0197] Example 3

[0198] This disclosure provides an electronic device comprising: a memory storing executable instructions; and a processor executing the executable instructions in the memory to implement the aforementioned hybrid location method for microseismic events in wells.

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

[0200] This memory is used to store non-transitory 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. The non-volatile memory may, for example, include read-only memory (ROM), hard disk, flash memory, etc.

[0201] The processor may be a central processing unit (CPU) or other form of processing unit with 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 this disclosure, the processor is used to execute computer-readable instructions stored in the memory.

[0202] Those skilled in the art will understand that, in order to solve the technical problem of how to achieve a good user experience, this embodiment may also include well-known structures such as communication buses and interfaces, and these well-known structures should also be included within the protection scope of this disclosure.

[0203] For a detailed description of this embodiment, please refer to the corresponding descriptions in the foregoing embodiments, which will not be repeated here.

[0204] Example 4

[0205] This disclosure provides a computer-readable storage medium storing a computer program that, when executed by a processor, implements the hybrid location method for microseismic events in wells.

[0206] A computer-readable storage medium according to embodiments of the present disclosure stores non-transitory computer-readable instructions. When these non-transitory computer-readable instructions are executed by a processor, all or part of the steps of the methods described in the foregoing embodiments of the present disclosure are performed.

[0207] The aforementioned 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 portable hard drive), media with built-in rewritable non-volatile memory (e.g., memory card), and media with built-in ROM (e.g., ROM cartridge).

[0208] Those skilled in the art should understand that the above description of the embodiments of the present invention is only intended to 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 of the examples given.

[0209] The various embodiments of the present invention have been described above. These descriptions are exemplary and not exhaustive, nor are they limited to the disclosed embodiments. Many modifications and variations will be apparent to those skilled in the art without departing from the scope and spirit of the described embodiments.

Claims

1. A method of mixed location of microseismic events in a well, characterized in that, The method comprises the following steps: calculating P-wave travel time difference, S-wave travel time difference and PS-wave travel time difference; constructing P-wave travel time difference index equation, S-wave travel time difference index equation and PS-wave travel time difference index equation respectively; calculating P-wave travel time positioning result, S-wave travel time positioning result and PS-wave travel time difference positioning result of the microseismic event respectively; determining the final microseismic event spatial positioning result according to the P-wave, S-wave and PS-wave travel time difference positioning results; wherein the P-wave travel time difference index equation is: （1） wherein, is the P-wave travel time time difference; wherein the S-wave travel time difference index equation is: （2） wherein, is the S-wave travel time time difference; wherein the PS-wave travel time difference index equation is: （3） wherein, is the PS wave travel time time difference; wherein a spatial triangle is established according to the P-wave, S-wave and PS-wave travel time difference positioning results, and the orthocenter position of the spatial triangle is calculated, which is the final microseismic event spatial positioning result; wherein the orthocenter position of the spatial triangle is calculated by formula (4): （4） wherein, is a P-wave traveltime positioning result, is an S-wave traveltime positioning result, is a PS-wave traveltime positioning result, is a final microseismic event spatial positioning result.

2. The hybrid microseismic event location method in a well according to claim 1, wherein, when the P-wave travel time difference index equation is the smallest, the corresponding grid point is the P-wave travel time positioning result through the grid search method; when the S-wave travel time difference index equation is the smallest, the corresponding grid point is the S-wave travel time positioning result; when the PS-wave travel time difference index equation is the smallest, the corresponding grid point is the PS-wave travel time difference positioning result.

3. A device for mixed positioning of microseismic events in a well, using the method for mixed positioning of microseismic events in a well according to claim 1 or 2, characterized in that The method comprises the following steps: a travel time difference calculation module for calculating P-wave travel time difference, S-wave travel time difference and PS-wave travel time difference; a construction module for constructing P-wave travel time difference index equation, S-wave travel time difference index equation and PS-wave travel time difference index equation respectively; a positioning result calculation module for calculating P-wave travel time positioning result, S-wave travel time positioning result and PS-wave travel time difference positioning result of the microseismic event respectively; a positioning module for determining the final microseismic event spatial positioning result according to the P-wave, S-wave and PS-wave travel time difference positioning results.

4. An electronic device, comprising: The electronic device comprises: a memory storing executable instructions; a processor running the executable instructions in the memory to implement the well microseismic event hybrid positioning method of claim 1 or 2.

5. A computer readable storage medium, characterized in that, The computer readable storage medium stores a computer program which is executed by the processor to implement the well microseismic event hybrid positioning method of claim 1 or 2.

Images