A multi-well joint monitoring microseismic positioning method and system

Abstract

The application provides a multi-well joint monitoring microseismic positioning method and system, and belongs to the field of microseismic data processing. The method uses microseismic data of multiple wells to perform multi-well joint microseismic positioning, and obtains a positioning result of a microseismic event. The method comprises the following steps: (1) inputting microseismic data of multiple wells, a P-S wave velocity model and initialization parameters; (2) multi-well microseismic event identification and matching; (3) P wave and S wave travel time picking of multi-well microseismic events; (4) obtaining a P wave polarization angle of a multi-well microseismic event at each geophone; (5) multi-well joint microseismic positioning to obtain a positioning result of the microseismic event; and (6) outputting the positioning result of the microseismic event. The multi-well joint positioning increases travel time and polarization information of multi-well data, and is beneficial to improving the accuracy of microseismic positioning; and the multi-well joint positioning improves the reliability of identifying microseismic events and the stability of the microseismic positioning method.

Description

TECHNICAL FIELD

[0001] The present application belongs to the field of microseismic data processing, and particularly relates to a microseismic positioning method and system for multi-well joint monitoring. BACKGROUND

[0002] Microseismic fracture monitoring technology has now become a commonly used monitoring technology in the development of tight reservoir oil and gas fields, and its most important application is the monitoring of hydraulic fracture distribution, and it is also used for monitoring the fractures generated by secondary exploitation such as oil and gas exploitation or water injection and gas injection. Microseismic technology has been a commonly used method in the fields of mine disaster monitoring and geothermal development since the 1970s and 1980s. Since the 21st century, with the rapid development of unconventional oil and gas, especially shale gas development, microseismic technology plays an important role in optimizing fracturing schemes and developing well network deployment, which enables the microseismic monitoring technology in oil and gas field development to develop rapidly.

[0003] Microseismic monitoring includes wellbore microseismic monitoring and surface microseismic monitoring. Surface microseismic monitoring and wellbore microseismic monitoring each have advantages and disadvantages. The surface monitoring method sets up several receiving points or measuring lines on the ground above the fracturing well section for microseismic monitoring. Compared with wellbore monitoring, the surface monitoring data signal energy is weak, and the signal-to-noise ratio is low due to the reasons such as more environmental noise interference, larger stratum absorption attenuation, and complex propagation path. However, the surface detector is easy to set up, the setting range is wide, and it does not need an observation well, avoiding the restriction of factors such as no observation well. The wellbore monitoring method sets up a detector string in the wellbore around the monitoring target area for microseismic monitoring. The wellbore detector is located in the wellbore, and the microseismic signal attenuation is small, the received microseismic signal is strong, and the environmental noise interference is small, so the signal-to-noise ratio is high. However, the horizontal direction positioning accuracy of single-well monitoring microseismic events is greatly affected by the monitoring distance, and the horizontal direction error increases with the increase of the distance between the event and the monitoring well. With the gradual popularization of ultra-long horizontal well fracturing, the monitoring distance of single-well microseismic monitoring is greatly increased, and the accuracy of fracturing crack characterization cannot effectively meet the engineering requirements.

[0004] Chinese patent publication CN103105622A discloses a same type wave travel time positioning method based on database technology, which relates to a geophysical exploration microseismic monitoring technology, and more particularly relates to a same type wave travel time positioning method combining database technology, which greatly improves the calculation efficiency while ensuring the calculation accuracy. The method mainly includes the following steps: (1) establishment and grid division of the microseismic occurrence space region; (2) using the ray tracing algorithm to write the P-wave or S-wave forward result of each grid into the database; (3) combining the actual microseismic occurrence space azimuth to determine whether the azimuth of the i-th grid Volume i (x, y, z) is within the range determined by the microseismic event azimuth (θ m -θ′, θm +θ′), if at a given azimuth angle (θ) m -θ′,θ m Within the range of +θ′), it is determined as a spatial point where microseismic events may occur. i (x, y, z, θ) m (4) Using the formula, in the spatial region (θ′) of the microseismic occurrence volume, determine its azimuth angle within (θ′). m -θ′,θ m Volume of all grid points within +θ′) i (x, y, z, θ) m The optimal seismic source location is searched using ±θ′).

[0005] However, the methods disclosed in the aforementioned patents have the following shortcomings:

[0006] (1) Utilize single microseismic data from wells or the surface;

[0007] (2) The positioning accuracy of single-well microseismic monitoring is low.

[0008] Considering the limitations of single-well microseismic monitoring, multi-well combined microseismic monitoring has emerged in recent years to overcome the application limitations of single-well microseismic monitoring in long horizontal well fracturing monitoring. Multi-well combined microseismic monitoring has its unique advantages: it improves the detection capability of microseismic events; in addition, it increases the observation angle, improving the location accuracy of microseismic events. Summary of the Invention

[0009] The purpose of this invention is to solve the problems existing in the prior art and provide a microseismic location method and system for multi-well joint monitoring. This method comprehensively utilizes the microseismic data information from multiple wells to perform multi-well joint microseismic location, improves the location accuracy of microseismic events, and further improves the overall accuracy of the spatial distribution of hydraulic fracturing fractures to meet engineering requirements.

[0010] This invention is achieved through the following technical solution:

[0011] In a first aspect, the present invention provides a microseismic location method based on multi-well joint monitoring, wherein the method utilizes microseismic data from multiple wells to perform multi-well joint microseismic location and obtain the location results of microseismic events.

[0012] A further improvement of the present invention is that:

[0013] The method includes:

[0014] (1) Input microseismic data, P-wave and S-wave velocity models and initialization parameters from multiple wells;

[0015] (2) Identification and matching of microseismic events in multiple wells;

[0016] (3) P-wave and S-wave travel time picking of multi-well microseismic events;

[0017] (4) Obtain the P-wave polarization angle of multi-well microseismic events in each detector;

[0018] (5) Multi-well joint microseismic location to obtain the location results of microseismic events;

[0019] (6) Output the location results of microseismic events.

[0020] A further improvement of the present invention is that:

[0021] The initialization parameters input in step (1) include: the range of the earthquake source location space, the size of the spatial discretization interval, and the weighting coefficient λ.

[0022] A further improvement of the present invention is that:

[0023] The operation of step (2) includes:

[0024] (21) Dynamic correction and stacking processing were performed on the microseismic data of each well to obtain the stacked data trace;

[0025] (22) The microseismic excitation time and magnitude were obtained from the superimposed data channels of each well;

[0026] (23) Search out all microseismic events in all wells that have the same excitation time and magnitude, and use them as the identified microseismic events.

[0027] A further improvement of the present invention is that:

[0028] The operation of step (3) includes:

[0029] P-wave and S-wave travel times were picked up for the identified microseismic events obtained in step (2).

[0030] A further improvement of the present invention is that:

[0031] The operation of step (4) includes:

[0032] Polarization analysis was performed on the identified microseismic events obtained in step (2) to obtain the P-wave polarization angle of each microseismic event in each detector.

[0033] A further improvement of the present invention is that:

[0034] The operation of step (5) includes:

[0035] A multi-well joint microseismic location objective function is established, and the location results of microseismic events are obtained using the multi-well joint microseismic location objective function.

[0036] A further improvement of the present invention is that:

[0037] The objective function for multi-well joint microseismic location established in step (5) is:

[0038] F = F t +λ*F θ (1)

[0039] Among them, F t F θ Let these represent the travel time objective function and the angle objective function, respectively.

[0040]

[0041] in, These represent the theoretical P-wave travel time and the theoretical S-wave travel time from the earthquake source location to detector i, respectively. The actual travel times of the P-wave and S-wave of the microseismic event are respectively represented, that is, the P-wave travel times and S-wave travel times obtained in step (3); This is the average of the actual P-wave travel times of all detectors; <t p > is the average of the theoretical P-wave travel times of all detectors; N r This refers to the total number of geophones in all wells during multi-well monitoring.

[0042]

[0043] Where, θ i This is the theoretical polarization angle at detector i; The actual horizontal polarization angle at detector i is the P-wave polarization angle obtained in step (4).

[0044] A further improvement of the present invention is that:

[0045] The weighting coefficient λ is a number greater than 0;

[0046] If neither time nor angle is given preference, the weighting coefficient is set to 1.

[0047] If a bias towards time travel is required, the weighting coefficient should be less than 1.

[0048] If an emphasis is needed on the angle, the weighting coefficient should be greater than 1.

[0049] A further improvement of the present invention is that:

[0050] The operation in step (5) of obtaining the location results of microseismic events using the multi-well joint microseismic location objective function includes:

[0051] (51) Input the P-wave travel time, S-wave travel time, P-wave polarization angle, P-wave and S-wave velocity model and initialization parameters of the microseismic events obtained in steps (3) and (4);

[0052] (52) The spatial discretization of the seismic source location is performed to obtain discrete seismic source points;

[0053] (53) Calculate the theoretical P-wave travel time, S-wave travel time and theoretical polarization angle from the discrete source point to each detector;

[0054] (54) Calculate the objective function value of multi-well joint microseismic location;

[0055] (55) Repeat steps (53)-(54) to obtain the objective function values ​​corresponding to all discrete source points;

[0056] (56) Sort the objective function values, and the coordinates of the discrete source point corresponding to the smallest objective function value are the location results of the microseismic event.

[0057] A second aspect of the present invention provides a microseismic location system for multi-well joint monitoring, the system comprising:

[0058] Input unit: Used to input microseismic data, P-wave and S-wave velocity models, and initialization parameters from multiple wells;

[0059] Identification unit: connected to the input unit, used for identifying and matching microseismic events from multiple wells;

[0060] Travel time acquisition unit: connected to the identification unit, used for P-wave and S-wave travel time acquisition of multi-well microseismic events;

[0061] Polarization angle acquisition unit: connected to the identification unit, used to obtain the P-wave polarization angle of multi-well microseismic events in each detector;

[0062] The positioning unit is connected to the input unit, the travel time acquisition unit, and the polarization angle acquisition unit, respectively, and is used to perform multi-well joint microseismic positioning to obtain the positioning results of microseismic events;

[0063] Output unit: Connected to the positioning unit, used to output the positioning results of microseismic events.

[0064] A third aspect of the present invention provides a computer-readable storage medium storing at least one computer-executable program, which, when executed by the computer, causes the computer to perform the steps in the above-described microseismic location method for multi-well joint monitoring.

[0065] Compared with the prior art, the beneficial effects of the present invention are:

[0066] (1) Multi-well joint positioning increases the travel time and polarization information of multiple well data, which is beneficial to improving the accuracy of microseismic positioning;

[0067] (2) Multi-well joint location improves the reliability of microseismic event identification and the stability of microseismic location methods. Attached Figure Description

[0068] Figure 1 A flowchart illustrating the steps of the method of this invention;

[0069] Figure 2 A flowchart of step (5) in the method of the present invention;

[0070] Figure 3 Single well location result diagram;

[0071] Figure 4 The positioning result diagram of the method of the present invention. Detailed Implementation

[0072] The present invention will now be described in further detail with reference to the accompanying drawings:

[0073] This invention provides a microseismic location method based on multi-well joint monitoring. The method utilizes microseismic data from multiple wells to perform multi-well joint microseismic location and obtain the location results of microseismic events.

[0074] like Figure 1 As shown, embodiments of the method of the present invention are as follows:

[0075] Example 1

[0076] The method includes:

[0077] (1) Input microseismic data, P-wave and S-wave velocity models and initialization parameters from multiple wells;

[0078] The initialization parameters input in step (1) include: the range of the seismic source location space, the size of the spatial discretization interval, and the weighting coefficient λ.

[0079] The spatial range for earthquake source location includes: Xmin, Ymin, Zmin, Xmax, Ymax, Zmax, which represent the minimum and maximum values ​​in the X, Y, and Z directions of the earthquake source location space, respectively.

[0080] The spatial discretization intervals include dx, dy, and dz, which represent the size of the discrete grid in the X, Y, and Z directions, respectively.

[0081] (2) Multi-well microseismic event identification and matching, as detailed below:

[0082] (21) Dynamic correction and stacking processing are performed on the microseismic data of each well to obtain the stacked data trace. This step can be implemented using existing technology and will not be described in detail here.

[0083] (22) The microseismic excitation time and magnitude are obtained from the stacked data channels of each well. This step can be achieved using existing technology and will not be described in detail here.

[0084] (23) Based on the criterion that the same microseismic event has the same excitation time and magnitude, the same microseismic event is identified and matched in the data of multiple wells. That is, the microseismic events with the same excitation time and magnitude in all wells are searched out and used as the identified microseismic events.

[0085] The following steps (3) to (6) are all for processing the microseismic events with the same time and magnitude identified in step (2).

[0086] (3) P-wave and S-wave travel time picking of multi-well microseismic events: P-wave and S-wave travel time (also known as travel time) are picked for the microseismic events identified in step (2); this step can be implemented using existing technology and will not be described in detail here.

[0087] (4) Perform polarization analysis on the identified microseismic events obtained in step (2) to obtain the P-wave polarization angle of each microseismic event in each detector. This step can be implemented using existing technology and will not be described in detail here.

[0088] (5) Multi-well joint microseismic location, as detailed below:

[0089] A multi-well joint microseismic location objective function is established, and the location result of the microseismic event is obtained using the multi-well joint microseismic location objective function:

[0090] Data on microseismic events were recorded in multiple wells. The P-wave and S-wave travel times and P-wave polarization angles obtained in steps (3) and (4) were used to establish a multi-well joint microseismic location objective function.

[0091] The established multi-well joint microseismic location objective function can be expressed as:

[0092] F = F t +λ*F θ (1)

[0093] Among them, F t F θ Let these represent the travel time objective function and the angle objective function, respectively.

[0094] λ is a weighting coefficient used to adjust the weights of the objective function of the angle, as follows:

[0095] The weighting coefficient λ is a number greater than 0;

[0096] If neither time nor angle is given preference, the weighting coefficient can be set to 1.

[0097] If a bias towards time travel is required, the weighting coefficient should be less than 1.

[0098] If an emphasis is needed on the angle, the weighting coefficient should be greater than 1.

[0099] The specific values ​​of the weighting coefficients can be determined after multiple tests based on the accuracy of the travel time and polarization angle of the microseismic data.

[0100]

[0101] in, These represent the theoretical P-wave travel time and the theoretical S-wave travel time from the source location obtained from theoretical calculations to detector i, respectively. The actual travel times of the P-wave and S-wave of the microseismic event are respectively represented as the P-wave and S-wave travel times obtained in step (3); This is the average (arithmetic mean) of the actual P-wave travel times of all detectors; <t p > is the average of the theoretical P-wave travel times of all detectors; N r This refers to the total number of geophones in all wells during multi-well monitoring.

[0102]

[0103] Where, θ i The horizontal polarization angle from the theoretically calculated source location to detector i is the theoretical polarization angle. The horizontal polarization angle at the actual detector i is the P-wave polarization angle obtained in step (4).

[0104] like Figure 2 As shown, the specific steps of step 5 are as follows:

[0105] (51) Input the P and S wave travel times and polarization angles of the microseismic events obtained in steps (3) and (4), the P and S wave velocity models and initialization parameters;

[0106] (52) Spatial discretization of the seismic source location is performed to obtain discrete seismic source points s. j j = 1, 2, ..., n, s j Let n be a discrete seismic source point, and n represent the total number of discrete seismic source points.

[0107]

[0108] (53), calculate discrete source point s jThe theoretical travel time (i.e., P-wave and S-wave theoretical travel time) and theoretical polarization angle of each detector i are calculated using existing ray tracing methods to obtain the discrete source point s. j The theoretical time (i.e., theoretical time of P-wave and S-wave) and theoretical polarization angle of each detector i;

[0109] (54) Calculate the objective function value of multi-well joint microseismic location;

[0110] (55) Repeat steps (53)-(54) to obtain the objective function values ​​corresponding to all discrete source points;

[0111] (56) Sort the objective function values, and the coordinates of the discrete source point corresponding to the smallest objective function value are the location results of the microseismic event.

[0112] When the objective function value is minimized, the travel time and polarization angle from the discrete source point to the detector are closest to the actual travel time and polarization angle. Theoretically, the difference between the two should be zero, but in reality, the difference can only be minimized. Therefore, the coordinates of the discrete source point corresponding to the minimum objective function value are the positioning results.

[0113] Step 6: Output the location results of the microseismic events.

[0114] The application examples of this invention are as follows:

[0115]

Example 2

[0116] According to the multi-well joint microseismic location method of the present invention, the method is verified using theoretical data. The actual location of the microseismic event (0m, 0m, 2400m) is taken, and random errors (±2ms) are added to the P-wave and S-wave travel times, and random errors (±1 degree) are added to the polarization angle. Then, the multi-well joint microseismic location method of the present invention is repeated 50 times, with different travel time and polarization errors each time. That is to say, it is equivalent to locating 50 microseismic events with the same theoretical location.

[0117] To test the method, this embodiment incorporates random errors into the theoretically calculated travel time and polarization angle of microseismic events, using these as the actual travel time and polarization angle (the specific random values ​​added each time are different, but the range of random values ​​remains consistent). Without random errors, all 50 positioning attempts should pinpoint the same true location. However, with random errors in travel time and polarization angle, the positioning should be within a certain range near the theoretical location; the better the positioning method, the smaller the corresponding range.

[0118] Figure 3 This is a top view showing the location results using only single-well data. Figure 4 The localization results obtained using this method are compared. Figure 3 andFigure 4 It can be observed that, Figure 4 The range is smaller, therefore theoretical data tests show that the positioning method of the present invention has a smaller positioning error, which is beneficial to improving the positioning accuracy of microseismic locations.

[0119] The present invention also provides a microseismic location system for multi-well joint monitoring, and an embodiment of the system is as follows:

[0120]

Example 3

[0121] The system includes:

[0122] Input unit: Used to input microseismic data, P-wave and S-wave velocity models, and initialization parameters from multiple wells;

[0123] Identification unit: connected to the input unit, used for identifying and matching microseismic events from multiple wells;

[0124] Travel time acquisition unit: connected to the identification unit, used for P-wave and S-wave travel time acquisition of multi-well microseismic events;

[0125] Polarization angle acquisition unit: connected to the identification unit, used to obtain the P-wave polarization angle of multi-well microseismic events in each detector;

[0126] The positioning unit is connected to the input unit, the travel time acquisition unit, and the polarization angle acquisition unit, respectively, and is used to perform multi-well joint microseismic positioning to obtain the positioning results of microseismic events;

[0127] Output unit: Connected to the positioning unit, used to output the positioning results of microseismic events.

[0128] Specifically, the initialization parameters input into the input unit include: the range of the seismic source location space, the size of the spatial discretization interval, and the weighting coefficient λ.

[0129] The spatial range for earthquake source location includes: Xmin, Ymin, Zmin, Xmax, Ymax, Zmax, which represent the minimum and maximum values ​​in the X, Y, and Z directions of the earthquake source location space, respectively.

[0130] The spatial discretization intervals include dx, dy, and dz, which represent the size of the discrete grid in the X, Y, and Z directions, respectively.

[0131] The identification unit performs the following processing:

[0132] (21) Dynamic correction and stacking processing are performed on the microseismic data of each well to obtain the stacked data trace. This step can be implemented using existing technology and will not be described in detail here.

[0133] (22) The microseismic excitation time and magnitude are obtained from the stacked data channels of each well. This step can be achieved using existing technology and will not be described in detail here.

[0134] (23) Based on the criterion that the same microseismic event has the same excitation time and magnitude, the same microseismic event is identified and matched in the data of multiple wells. That is, the microseismic events with the same excitation time and magnitude in all wells are searched out and used as the identified microseismic events.

[0135] The time-travel pickup unit performs the following processing:

[0136] P-wave and S-wave travel times were picked up from the microseismic events identified by the identification unit.

[0137] The polarization angle acquisition unit performs the following processing:

[0138] Polarization analysis was performed on the microseismic events identified by the identification unit to obtain the P-wave polarization angle of each microseismic event in each detector.

[0139] The positioning unit performs the following processing:

[0140] A multi-well joint microseismic positioning objective function is established using the P-wave and S-wave travel times and P-wave polarization angles obtained from the travel time picking unit and polarization angle acquisition unit.

[0141] Discrete source points are obtained by spatial discretization of the earthquake source location.

[0142] Calculate the theoretical travel time and theoretical polarization angle from the discrete seismic source point to each detector;

[0143] Obtain the objective function values ​​corresponding to all discrete seismic source points;

[0144] The objective function values ​​are sorted, and the coordinates of the discrete source point corresponding to the smallest objective function value are the location results of the microseismic event.

[0145] This invention utilizes the P-wave and S-wave travel time information of microseismic events from multi-well data to achieve high-precision localization of microseismic events.

[0146] Finally, it should be noted that the above technical solution is only one embodiment of the present invention. For those skilled in the art, based on the application methods and principles disclosed in the present invention, it is easy to make various types of improvements or modifications, and not limited to the methods described in the above specific embodiments of the present invention. Therefore, the methods described above are only preferred and have no limiting significance.

Claims

1. A microseismic location method based on multi-well joint monitoring, characterized in that: The method utilizes microseismic data from multiple wells to perform multi-well joint microseismic localization, thereby obtaining the localization results of microseismic events. The method includes: (1) Input microseismic data, P-wave and S-wave velocity models and initialization parameters from multiple wells; (2) Identification and matching of microseismic events in multiple wells; (3) P-wave and S-wave travel time picking of multi-well microseismic events; (4) Obtain the P-wave polarization angle of multi-well microseismic events at each detector; (5) Multi-well joint microseismic location to obtain the location results of microseismic events; (6) Output the location results of microseismic events; Step (5) includes the following operations: Establish a multi-well joint microseismic location objective function; （1） in, Indicates the weighting coefficient. Let these represent the travel time objective function and the angle objective function, respectively. （2） in, , These represent the distance from the seismic source to the detector. P-wave theory timekeeping and S-wave theory timekeeping; , The actual travel times of the P-wave and S-wave of the microseismic event are respectively represented as the P-wave travel times and S-wave travel times obtained in step (3); This is the average of the actual P-wave travel times of all detectors; The average of the theoretical P-wave travel times of all detectors; This refers to the total number of geophones in all wells during multi-well monitoring. （3） in, For detector The theoretical polarization angle at that location; For detector The actual horizontal polarization angle at that point, i.e., the P-wave polarization angle obtained in step (4); The location results of microseismic events are obtained by using a multi-well combined microseismic location objective function, including: (51) Input the P-wave travel time, S-wave travel time, P-wave polarization angle, P-wave and S-wave velocity model and initialization parameters of the microseismic event obtained in steps (3) and (4); (52) Spatial discretization of the seismic source location is performed to obtain discrete seismic source points; (53) Calculate the theoretical travel time of P-wave, theoretical travel time of S-wave, and theoretical polarization angle from the discrete source point to each detector; (54) Calculate the objective function value for multi-well joint microseismic location; (55) Repeat steps (53)-(54) to obtain the objective function values ​​corresponding to all discrete seismic source points; (56) Sort the objective function values, and the coordinates of the discrete source point corresponding to the smallest objective function value are the location results of the microseismic event.

2. The microseismic location method based on multi-well joint monitoring according to claim 1, characterized in that: The initialization parameters input in step (1) include: the range of the seismic source location space, the size of the spatial discretization interval, and the weighting coefficient. .

3. The microseismic location method based on multi-well joint monitoring according to claim 1, characterized in that: The operation of step (2) includes: (21) Dynamic correction and stacking processing were performed on the microseismic data of each well to obtain the stacked data trace; (22) The microseismic excitation time and magnitude were obtained from the stacked data channels of each well; (23) Search out all microseismic events in all wells that have the same excitation time and magnitude, and use them as the identified microseismic events.

4. The microseismic location method based on multi-well joint monitoring according to claim 3, characterized in that: P-wave and S-wave travel times were picked up for the identified microseismic events obtained in step (2).

5. The microseismic location method based on multi-well joint monitoring according to claim 1, characterized in that: The operation of step (4) includes: Polarization analysis was performed on the identified microseismic events obtained in step (2) to obtain the P-wave polarization angle of each microseismic event in each detector.

6. The microseismic location method based on multi-well joint monitoring according to claim 1, characterized in that: The weighting coefficients Numbers greater than 0; If neither time nor angle is given preference, the weighting coefficient is set to 1. If a bias towards time travel is required, the weighting coefficient should be less than 1. If an emphasis is needed on the angle, the weighting coefficient should be greater than 1.

7. A microseismic location system based on multi-well joint monitoring, wherein the microseismic location method based on multi-well joint monitoring is described in any one of claims 1-6, characterized in that: The system includes: Input unit: Used to input microseismic data, P-wave and S-wave velocity models, and initialization parameters from multiple wells; Identification unit: connected to the input unit, used for identifying and matching microseismic events from multiple wells; Travel time acquisition unit: connected to the identification unit, used for P-wave and S-wave travel time acquisition of multi-well microseismic events; Polarization angle acquisition unit: connected to the identification unit, used to obtain the P-wave polarization angle of multi-well microseismic events in each detector; The positioning unit is connected to the input unit, the travel time acquisition unit, and the polarization angle acquisition unit, respectively, and is used to perform multi-well joint microseismic positioning to obtain the positioning results of microseismic events; Output unit: Connected to the positioning unit, used to output the positioning results of microseismic events.

8. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores at least one computer-executable program, which, when executed by the computer, causes the computer to perform the steps in the microseismic location method for multi-well joint monitoring as described in any one of claims 1-6.

Images