A horizontal well fracturing microseismic monitoring s-wave segment energy compensation method and system

Abstract

The application provides a horizontal well fracturing microseismic monitoring S-wave segmented energy compensation method and system, and belongs to the microseismic signal processing field.The method firstly constructs a microseismic event S-wave first arrival picking function to obtain S-wave signal energy, then uses microseismic event positions and detector positions to segmentally construct a linear energy compensation function and a nonlinear energy compensation function of the microseismic event S-wave, finally inputs the spatial position of any microseismic event, obtains the S-wave signal energy compensation factor of the microseismic event according to the horizontal monitoring distance, and realizes the horizontal well fracturing microseismic monitoring S-wave signal energy compensation by using the S-wave signal energy compensation factor.The application overcomes the problem of weak signals caused by the observation distance, and corrects the S-wave signal energy to a reasonable observation distance according to the S-wave signal energy law, so that the microseismic signals received at different positions by the same level of seismic source point are basically consistent.

Description

Technical Field

[0001] The present invention belongs to the field of microseismic signal processing, and in particular relates to an S-wave segmented energy compensation method and system for horizontal well fracturing microseismic monitoring. Background Art

[0002] In actual fracturing monitoring, the primary source of earthquakes generated by rock compression is often shear waves. Microseismic events generate both P-wave and S-wave signals, but S-wave signals are far more powerful than P-waves, making S-waves a more reliable source for microseismic event research. Typically, microseismic monitoring data refers to the signals received by a geophone as a single P-wave, a single S-wave, or a combination of P and S waves during the artificial fractures created by hydraulic fracturing.

[0003] Microseismic monitoring generally involves both surface and borehole microseismic monitoring. Downhole three-component geophones are placed in the observation well section to receive full-wavefield microseismic signals. Compared to surface microseismic monitoring, borehole data has a higher signal-to-noise ratio and a richer number and types of microseismic events. Compared to single-component surface microseismic monitoring, three-component borehole microseismic monitoring is more sensitive to S-waves, meaning more and stronger S-wave signals are detected in the well.

[0004] However, due to the limited number of microseismic geophones in the well (generally 12 to 32 three-component geophones are set), the monitoring range is small. Moreover, due to the different monitoring distances, the S waves of microseismic events received by the geophones vary greatly, causing the well microseismic positioning method to easily produce unstable microseismic interpretation results.

[0005] Energy compensation for weak microseismic signals can reduce the impact of microseismic monitoring distance in wells. Current methods for signal enhancement include: removing random noise and linear interference to improve the signal-to-noise ratio (SNR), thereby enhancing the effective signal; and performing energy compensation based on surface or formation absorption factors.

[0006] According to the principle of conservation of energy, when S-waves are transmitted and reflected at an interface, the smaller the transmission angle, the stronger the transmitted wave energy and the weaker the reflected wave energy. Conversely, the larger the transmission angle, the weaker the transmitted wave energy and the stronger the reflected wave energy. In microseismic observations in horizontal fracturing wells, the geophone primarily receives transmitted waves. The closer the monitoring distance and the smaller the transmission angle, the stronger the received S-wave signal energy. Conversely, the farther the monitoring distance and the larger the transmission angle, the weaker the received S-wave signal energy. This phenomenon results in significant differences in the energy of the S-wave signals received by the geophone when the same energy level source is excited at different monitoring distances, leading to erroneous interpretations of microseismic events. Summary of the Invention

[0007] The present application aims to solve the problems existing in the prior art, and provides a horizontal well fracturing microseismic monitoring S-wave segmented energy compensation method and system, which eliminates or partially eliminates the influence of monitoring distance on S-wave signal energy, so that the S-wave signal energy received by the geophone from the same energy level seismic source excited at different monitoring distances is equivalent or close.

[0008] The present application is realized by the following technical solutions:

[0009] In a first aspect, the present application provides a horizontal well fracturing microseismic monitoring S-wave segmented energy compensation method, which first constructs a microseismic event S-wave first arrival picking function to obtain S-wave signal energy, then uses the microseismic event position and the geophone position to segmentally construct a linear energy compensation function and a nonlinear energy compensation function of the microseismic event S-wave, finally inputs the spatial position of any microseismic event, obtains the S-wave signal energy compensation factor of the microseismic event according to the horizontal monitoring distance, and uses the S-wave signal energy compensation factor to realize the S-wave signal energy compensation of the horizontal well fracturing microseismic monitoring.

[0010] Preferably, the method comprises:

[0011] (1) inputting the original microseismic event S-wave signal, picking up the S-wave first arrival time and the S-wave signal energy;

[0012] (2) segmentally constructing the S-wave linear energy compensation function and the S-wave nonlinear energy compensation function;

[0013] (3) using the S-wave linear energy compensation function or the S-wave nonlinear energy compensation function to compensate the energy of the S-wave of any microseismic event.

[0014] Preferably, the operation of (1) comprises:

[0015] (11) performing data processing on the microseismic data to obtain the microseismic event S-wave data and spatial position at N monitoring distances from near to far;

[0016] (12) constructing the S-wave first arrival picking function to pick up the S-wave first arrival time and the S-wave signal energy.

[0017] Preferably, the operation of (12) comprises:

[0018] (121) defining a long time window L1 and a short time window L2, the time window length of L1 is greater than that of L2, and constructing an S-wave first arrival picking function based on the S-wave amplitude weight coefficient:

[0019]

[0020]

[0021] where i represents the i-th microseismic sampling event, j represents the j-th time sample, N is the total amount of microseismic sampling data, M is the total amount of time samples, AS i,j , ER Si,j , fS i (t j ) are the S-wave amplitude value, the long-short time window energy ratio, and the S-wave first arrival picking function of the i-th microseismic sampling event at the j-th time sample, respectively;

[0022] (122), the S-wave first arrival time t * is obtained by searching the maximum value of the S-wave first arrival picking function;

[0023] (123), a time window W S is defined, and the microseismic event S-wave signal energy SEN i is obtained by using the following formula:

[0024]

[0025] Preferably, the size of the time window W S is the length of time sample corresponding to the length of a wavelet.

[0026] Preferably, the operation of (2) comprises:

[0027] (21), the horizontal monitoring distance and the receiver transmission angle of the N microseismic events are calculated by using the following formula:

[0028] L Si = |Sx i -Rx|, (i = 1, …, N) (4)

[0029]

[0030] where (Sx i , Sy i ) is the spatial position of the microseismic event, (Rx, Ry) is the coordinate of the receiver, L Si , θ Si are the horizontal monitoring distance and the receiver transmission angle of the microseismic event, respectively;

[0031] (22), the independent variable function g(θ Si , L Si ) combined by the receiver transmission angle and the horizontal monitoring distance is defined:

[0032]

[0033] The values of the function g(θ Si , L Si ) corresponding to the N microseismic events are calculated;

[0034] (23), constructing the S-wave linear energy compensation function fSEN1(θ Si , L Si ) of the microseismic event with the horizontal monitoring distance greater than 1000 meters:

[0035] fSEN1(θ Si , L Si ) = C1*g(θ Si , L Si ) + C2, (i = N * +1, ……N) (7)

[0036] wherein, when N * +1≤i≤N, L Si >1000, C1, C2 are linear coefficients of the compensation function fSEN1(θ Si , L Si );

[0037] (24), using the S-wave signal energy SEN * of the microseismic event with the serial number from (N i +1) to N corresponding to the monitoring distance from near to far, forming an over-determined equation group:

[0038]

[0039] (25), solving the over-determined equation group to obtain the values of the coefficients C1, C2;

[0040] (26), constructing the S-wave nonlinear energy compensation function fSEN2(θ Si , L Si ) of the microseismic event with the horizontal monitoring distance less than or equal to 1000 meters:

[0041] fSEN2(θ Si , L Si ) = D1*(g(θ Si , L Si )) 2 +D2*g(θ Si , L Si )+D3, (i = 1, ……N * ) (9)

[0042] wherein, when i≤N * , L Si ≤1000, D1, D2, D3 are nonlinear coefficients of the compensation function fSEN2(θ Si , L Si );

[0043] (27), using the S-wave signal energy of the microseismic event with the serial number from 1 to N *The S-wave signal energy SEN of the microseismic event i , forming an overdetermined system of equations:

[0044]

[0045] (28), solve the overdetermined system of equations and obtain the values ​​of coefficients D1, D2, and D3.

[0046] Preferably, in (25) and (28), the singular value decomposition (SVD) method is used to solve the overdetermined system of equations.

[0047] Preferably, the operation of (3) includes:

[0048] (31), assuming that the spatial location of the microseismic event is (event X , event Y ), and calculate the monitoring distance L of the microseismic event respectively. event , detector transmission angle θ event With the independent variable function g(θ event , L event );

[0049] (32), when L event When the distance is greater than 1000 meters, the S-wave signal energy compensation factor SQ of the microseismic event is obtained using the following formula: * :

[0050]

[0051] When L event When the distance is less than 1000 meters, the energy compensation factor SQ of the microseismic event S wave signal is obtained using the following formula: * :

[0052]

[0053] (33) The root mean square of the S-wave signal energy compensation factor is multiplied by the S-wave amplitude value of the microseismic event to obtain the compensated S-wave amplitude.

[0054] A second aspect of the present invention provides an S-wave segmented energy compensation system for microseismic monitoring of horizontal well fracturing, the system comprising:

[0055] An input and pickup unit, used for inputting the original microseismic event S-wave signal and picking up the S-wave first arrival time and S-wave signal energy;

[0056] A segmented construction unit is connected to the input and pickup unit and is used to segmentally construct an S-wave linear energy compensation function and an S-wave nonlinear energy compensation function;

[0057] The segmented compensation unit is connected to the segmented construction unit and is used to perform energy compensation on the S wave of any microseismic event by using an S wave linear energy compensation function or an S wave nonlinear energy compensation function.

[0058] The third aspect of the present invention provides a computer-readable storage medium, which stores at least one computer-executable program. When the at least one program is executed by the computer, the computer executes the steps of the S-wave segmented energy compensation method for horizontal well fracturing microseismic monitoring of the present invention.

[0059] Compared with the existing technology, the beneficial effects of the present invention are: the present invention overcomes the problem of weak signal caused by observation distance, and corrects the S-wave signal energy to a reasonable observation distance in segments according to the law of S-wave signal energy, so that the microseismic signals received at different locations of the same level of source point are roughly consistent, providing reliable guarantee for subsequent microseismic event magnitude statistics. BRIEF DESCRIPTION OF THE DRAWINGS

[0060] Figure 1 It is a block diagram of the steps of the method of the present invention;

[0061] Figure 2 It is a theoretical model of microseismicity, which includes known velocity layers, event locations, and geophone locations;

[0062] Figure 3 It is the S-wave data profile of the original microseismic event after forward modeling and signal processing;

[0063] Figure 4 It is the intersection diagram of the energy of the microseismic S-wave linear energy compensation function and the independent variable;

[0064] Figure 5 It is the intersection diagram of the energy of the microseismic S-wave nonlinear energy compensation function and the independent variable;

[0065] Figure 6 yes Figure 3 S-wave data profile of a microseismic event after energy compensation processing. DETAILED DESCRIPTION

[0066] The present invention is further described in detail below with reference to the accompanying drawings:

[0067] The application provides a horizontal well fracturing microseismic monitoring S-wave segmented energy compensation method, which comprises the following three steps: firstly, constructing a microseismic event S-wave first arrival picking function to obtain S-wave signal energy; secondly, using known microseismic event positions and detector positions to segmentally construct a linear energy compensation function and a nonlinear energy compensation function of the microseismic event S-wave; and finally, inputting the spatial position of any microseismic event, obtaining the S-wave signal energy compensation factor of the microseismic event according to the horizontal monitoring distance, and using the S-wave signal energy compensation factor to realize the S-wave signal energy compensation of the horizontal well fracturing microseismic monitoring.

[0068] As shown in the figure, the application provides a horizontal well fracturing microseismic monitoring S-wave segmented energy compensation method, which comprises the following three steps: Figure 1

[0069] (1), inputting an original microseismic event S-wave signal, picking up S-wave first arrival time and S-wave signal energy;

[0070] (2), segmentally constructing an S-wave linear energy compensation function and an S-wave nonlinear energy compensation function;

[0071] (3), using the S-wave linear energy compensation function or the S-wave nonlinear energy compensation function to perform energy compensation on the S-wave of any microseismic event.

[0072] The implementation of the method of the application is as follows:

[0073] Example 1

[0074] The operation of (1) comprises the following steps:

[0075] (11), using an existing method to perform data processing on microseismic data to obtain microseismic event S-wave data and spatial positions at N monitoring distances from near to far, and taking these as input data;

[0076] (12), constructing an S-wave first arrival picking function to pick up S-wave first arrival time and S-wave signal energy SEN i .

[0077] The operation of (12) comprises the following steps:

[0078] (121), defining a long time window L1 and a short time window L2, the time window length of L1 is greater than that of L2, and constructing an S-wave first arrival picking function based on an S-wave amplitude weight coefficient:

[0079]

[0080]

[0081] wherein i represents the i-th microseismic sampling event, j represents the j-th time sample, N is the total amount of microseismic sampling data, M is the total amount of time samples, and AS​i,j , ER Si,j , fS i (t j ) are the amplitude value of S-wave, long-short time window energy ratio, S-wave first arrival picking function of the jth time sample of the ith microseismic sampling event, respectively;

[0082] (122), the S-wave first arrival time is obtained by searching the maximum value of the S-wave first arrival picking function: when t j = t * , the value of fS i (t j ) is the maximum, and t * is the S-wave first arrival time of the ith microseismic sampling event.

[0083] (123), a time window W S is defined (generally with a length of a subwave length time sample), and the S-wave signal energy SEN i of the microseismic event is obtained by using the following formula:

[0084]

[0085] Embodiment two:

[0086] The operation of (2) comprises:

[0087] The spatial positions of the N microseismic events are divided into two groups according to the order of the horizontal monitoring distance being greater than 1000 meters and being less than or equal to 1000 meters, and then the S-wave linear energy compensation function fSEN1(θ Si , L Si ) and the S-wave nonlinear energy compensation function fSEN2(θ Si , L Si ) are respectively constructed with the detector transmission angle θ Si and the monitoring distance L Si as independent variables. Specifically, it comprises:

[0088] (21), the horizontal monitoring distance and the detector transmission angle of the N microseismic events are calculated by using the following formula:

[0089] L Si = |Sx i -Rx|, (i = 1, ……N) (4)

[0090]

[0091] Wherein, (Sx i , Sy i ) is the spatial position of the known microseismic event, (Rx, Ry) is the coordinate of the known detector, L Si , θ Sirespectively, are the horizontal monitoring distance and the detector transmission angle of a microseismic event.

[0092] (22), defining the independent variable function g(θ Si ,L Si ) of the combination of the detector transmission angle and the horizontal monitoring distance.

[0093]

[0094] The values of the function g(θ Si ,L Si ) corresponding to N microseismic events are calculated.

[0095] (23), constructing the S-wave linear energy compensation function fSEN1(θ Si ,L Si ) of the microseismic event with horizontal monitoring distance greater than 1000 meters:

[0096] fSEN1(θ Si ,L Si ) = C1*g(θ Si ,L Si ) + C2, (i = N * +1, ……N) (7)

[0097] Wherein, when N * +1≤i≤N, L Si >1000, C1, C2 are linear coefficients of the compensation function fSEN1(θ Si ,L Si ).

[0098] (24), using the S-wave signal energy SEN * of the microseismic event with serial number from (N i +1) to N corresponding to the monitoring distance from near to far, forming an overdetermined equation group:

[0099]

[0100] (25), using the existing singular value decomposition SVD method to solve the above overdetermined equation group (process omitted), calculating the coefficients C1, C2, and realizing the construction of the S-wave linear energy compensation function fSEN1(θ Si ,L Si ) of the microseismic event.

[0101] (26), constructing the S-wave nonlinear energy compensation function fSEN2(θ Si ,L Si ) of the microseismic event with horizontal monitoring distance less than or equal to 1000 meters:

[0102] fSEN2(θ SiL Si ) = D1 * (g (θ Si ,L Si )) 2 + D2 * g (θ Si ,L Si ) + D3, (i = 1, ……, N * ) (9)

[0103] Wherein, when i≤N * , L Si ≤1000, D1, D2, D3 are nonlinear coefficients of compensation function fSEN2(θ Si ,L Si )。

[0104] (27), using monitoring distance from near to far corresponding to the microseismic event S-wave signal energy SEN * of serial number 1 to N i , form an overdetermined equation group:

[0105]

[0106] (28), using existing singular value decomposition SVD method to solve the overdetermined equation group (process is omitted), the coefficients D1, D2, D3 are calculated, that is, the construction of microseismic event S-wave nonlinear energy compensation function fSEN2(θ Si ,L Si ) is realized.

[0107] Embodiment three:

[0108] The operation of (3) comprises:

[0109] Through long-term seismic data processing, it is found that the S-wave signal energy varies with the monitoring distance and the transmission angle and presents different rules, when the monitoring distance is greater than 1000 meters, it presents a linear rule, otherwise, when the monitoring distance is less than or equal to 1000 meters, it presents a nonlinear rule, therefore, the compensation function is processed in sections, that is, the compensation function is constructed in two parts of linear and nonlinear, so that the S-wave energy compensation of different monitoring distances can be more targeted.

[0110] Specifically, input any one microseismic event, calculate the geophone transmission angle Monitoring distance and compare with 1000 meters, then substitute into the corresponding S-wave linear or nonlinear energy compensation function, construct the microseismic event S-wave signal energy compensation factor SQ * , and the event S-wave amplitude value is compared with Point multiplication is performed to finally obtain the energy-compensated microseismic event S-wave signal, and the specific steps include:

[0111] (31), assuming that the spatial position of the microseismic event is (event X , event Y ), the monitoring distance L event , the receiver transmission angle θ event and the argument function g(θ event , L event ) of the microseismic event are respectively calculated by using the above formulas (4)-(6):

[0112] L event = |event X -Rx|

[0113]

[0114]

[0115] (32), when L event >1000 meters, the S-wave signal energy compensation factor SQ * of the microseismic event is obtained by using the following formula:

[0116]

[0117] When L event ≤1000 meters, the S-wave signal energy compensation factor SQ * of the microseismic event is obtained by using the following formula:

[0118]

[0119] (33), the root mean square of the S-wave signal energy compensation factor (i.e. ) is multiplied by the S-wave amplitude value of the microseismic event to obtain the compensated S-wave amplitude, and thus the S-wave signal energy compensation of the microseismic event is finally realized.

[0120] Embodiment four:

[0121] The application also provides a horizontal well fracturing microseismic monitoring S-wave segmented energy compensation system, which comprises:

[0122] An input and pickup unit is used for inputting an original microseismic event S-wave signal and picking up S-wave first arrival time and S-wave signal energy;

[0123] A segmented construction unit is connected with the input and pickup unit and is used for segmented construction of an S-wave linear energy compensation function and an S-wave nonlinear energy compensation function;

[0124] The segmented compensation unit is connected to the segmented construction unit and is used to perform energy compensation on the S wave of any microseismic event by using an S wave linear energy compensation function or an S wave nonlinear energy compensation function.

[0125] The method of the present invention is tested using theoretical model data, and the energy compensation effect of the present invention on microseismic S-wave signals is analyzed.

[0126] First, a microseismic monitoring model with different monitoring distances in the well is designed (e.g. Figure 2 ), and set 1 detector position ( Figure 2 white triangle) and 150 event locations ( Figure 2 The horizontal distances from microseismic events to geophones range from 150 to 3130 meters, with the intervals between events being 20 meters.

[0127] Using the differential wave equation, we simulate the three-component microseismic data in the well. After processing the three-component vector signal of the microseismic data, we can obtain 150 microseismic S-wave data, such as Figure 3 As shown in Figure 2, the energy of the S-wave signal of the microseismic event is calculated using formulas (1) to (3).

[0128] Then, the microseismic S-wave segmented energy compensation function of the theoretical model is constructed, including linear compensation function and nonlinear compensation function: According to formulas (4) to (6), the transmission angle, horizontal monitoring distance and independent variable function of the microseismic event detector of 150 known monitoring distances from near to far are calculated. Then, according to whether the horizontal monitoring distance is greater than 1000 meters, the formulas (7) and (9) are respectively substituted to establish two sets of overdetermined matrices (such as formula (8) and formula (10)). The singular value decomposition SVD is used to solve the corresponding linear coefficients C1, C2 and nonlinear coefficients D1, D2, D3, and then the microseismic event S-wave signal energy compensation function is successfully constructed. Figure 4 、 Figure 5 As shown in the figure, the gray dots are the intersection diagram of the energy compensation function of the microseismic S-wave signal and the independent variable, and the corresponding black straight lines are the S-wave linear and nonlinear energy compensation function curves.

[0129] Finally, according to the S-wave signal energy compensation function calculated in the previous step, the energy compensation is performed on the S-wave data of 150 original microseismic events, as shown in the following example: Figure 6 Compare Figure 3 From the original S data in the image, it can be seen that the S-wave signal energy is well compensated. The farther the distance, the better the compensation effect, especially when the horizontal distance is above 1000 meters. Even when the horizontal distance reaches about 3000 meters, the S-wave signal becomes stronger from weak, which verifies that the present invention has a good compensation effect on the microseismic S-wave signal energy.

[0130] In the description of the present application, unless otherwise specified, the terms "upper", "lower", "left", "right", "inner", "outer" and the like indicate the orientation or positional relationship based on the orientation or positional relationship shown in the drawings, and are only for the convenience of describing the present application and simplifying the description, and do not indicate or imply that the devices or elements referred to must have a particular orientation, be constructed and operated in a particular orientation, and therefore cannot be understood as a limitation on the present application.

[0131] The technical solution described above is only one embodiment of the present application. For those skilled in the art, on the basis of the principles disclosed in the present application, various types of improvements or modifications can be easily made, and are not limited to the technical solution described in the above specific embodiments of the present application. Therefore, the foregoing description is only preferred, and is not limiting.

Claims

1. A method for S-wave segmented energy compensation in microseismic monitoring of horizontal well fracturing, characterized by: The method first constructs a microseismic event S-wave first arrival picking function to obtain S-wave signal energy, then uses the microseismic event position and the geophone position to segmentally construct a linear energy compensation function and a nonlinear energy compensation function for the microseismic event S-wave. Finally, the spatial position of any microseismic event is input, and the S-wave signal energy compensation factor of the microseismic event is obtained according to the horizontal monitoring distance. The S-wave signal energy compensation factor is used to realize S-wave signal energy compensation in horizontal well fracturing microseismic monitoring. The method comprises: (1) Input the original microseismic event S-wave signal, pick up the S-wave first arrival time and S-wave signal energy; (2) S-wave linear energy compensation function and S-wave nonlinear energy compensation function are constructed in segments; (3) Using the S-wave linear energy compensation function or the S-wave nonlinear energy compensation function to compensate the S-wave energy of any microseismic event; The operation of step (2) includes: (21) The horizontal monitoring distance and the geophone transmission angle of N microseismic events are calculated using the following formula: L Si =|Sx i -Rx|,i=1,……N (4) Among them, (Sx i , Sy i ) is the spatial position of the microseismic event, (Rx, Ry) is the coordinate of the detector, L Si ,θ Si are the horizontal monitoring distance and the geophone transmission angle of microseismic events, respectively; (22), define the independent variable function g(θ) of the combination of the detector transmission angle and the horizontal monitoring distance Si ,L Si ): Calculate the function g(θ corresponding to N microseismic events Si ,L Si ) (23) Construct the S-wave linear energy compensation function fSEN1(θ Si ,L Si ): fSEN1(θ Si ,L Si )=C1*g(θ Si ,L Si )+C2,i=N*+1,……N(7) Among them, when N * +1≤i≤N, L Si >1000, C1, C2 are the compensation function fSEN1(θ Si ,L Si )’s linear coefficient; (24), using the monitoring distance from near to far corresponding to the sequence number N * +1 to N microseismic events S wave signal energy SEN i , forming an overdetermined system of equations: (25), solve the overdetermined system of equations and obtain the values ​​of coefficients C1 and C2; (26) Construct the S-wave nonlinear energy compensation function fSEN2(θ Si ,L Si ): fSEN2(θ Si ,L Si )=D1*(g(θ Si ,L Si )) 2 +D2*g(θ Si ,L Si )+D3,i=1,……N*(9) Among them, when i≤N * When L Si ≤1000, D1, D2, D3 are the compensation function fSEN2(θ Si ,L Si )’s nonlinear coefficient; (27), using the monitoring distance from near to far corresponding to the sequence number 1 to N * The S-wave signal energy SEN of the microseismic event i , forming an overdetermined system of equations: (28), solve the overdetermined system of equations and obtain the values ​​of coefficients D1, D2, and D3.

2. The S-wave segmented energy compensation method for horizontal well fracturing microseismic monitoring according to claim 1 is characterized by: The operations of (1) include: (11) Processing the microseismic data to obtain the S-wave data and spatial positions of N microseismic events with monitoring distances from near to far; (12) , construct the S-wave first arrival picking function, pick up the S-wave first arrival time and S-wave signal energy.

3. The S-wave segmented energy compensation method for horizontal well fracturing microseismic monitoring according to claim 2 is characterized by: The operation of (12) includes: (121), define a long time window L1 and a short time window L2, where the time window length of L1 is greater than the time window length of L2, and construct an S-wave first arrival picking function based on the S-wave amplitude weight coefficient: Where i represents the i-th microseismic sampling event, j represents the j-th time sampling point, N is the total amount of microseismic sampling data, M is the total amount of time sampling points, and AS i,j , ER Si,j 、fS i (t j ) are the S-wave amplitude, the energy ratio of the long-short time window, and the S-wave first arrival picking function corresponding to the j-th time sample point of the i-th microseismic sampling event; (122), the S-wave first arrival time t is obtained by searching for the maximum value of the S-wave first arrival picking function * ; (123), define a time window W S , the S-wave signal energy SEN of the microseismic event is obtained using the following formula: i :

4. The S-wave segmented energy compensation method for horizontal well fracturing microseismic monitoring according to claim 3 is characterized by: Time window W S The size of is the time sample length corresponding to the length of a wavelet.

5. The S-wave segmented energy compensation method for horizontal well fracturing microseismic monitoring according to claim 1 is characterized by: In steps (25) and (28), the singular value decomposition (SVD) method is used to solve the overdetermined system of equations.

6. The S-wave segmented energy compensation method for horizontal well fracturing microseismic monitoring according to claim 1 is characterized by: The operations of (3) include: (31), assuming that the spatial location of the microseismic event is (event X , event Y ), and calculate the monitoring distance L of the microseismic event respectively. event , detector transmission angle θ event With the independent variable function g(θ event , L event ); (32), when L event When the distance is greater than 1000 meters, the S-wave signal energy compensation factor SQ of the microseismic event is obtained using the following formula: * : When L event When the distance is less than 1000 meters, the energy compensation factor SQ of the microseismic event S wave signal is obtained using the following formula: * : (33) The root mean square of the S-wave signal energy compensation factor is multiplied by the S-wave amplitude value of the microseismic event to obtain the compensated S-wave amplitude.

7. A horizontal well fracturing microseismic monitoring S-wave segmented energy compensation system, characterized by: The system comprises: An input and pickup unit, used for inputting the original microseismic event S-wave signal and picking up the S-wave first arrival time and S-wave signal energy; A segmented construction unit is connected to the input and pickup unit and is used to segmentally construct an S-wave linear energy compensation function and an S-wave nonlinear energy compensation function; A segmented compensation unit is connected to the segmented construction unit and is used to perform energy compensation on the S wave of any microseismic event using an S wave linear energy compensation function or an S wave nonlinear energy compensation function; The segment construction unit performs the following operations: (21) The horizontal monitoring distance and the geophone transmission angle of N microseismic events are calculated using the following formula: L Si =|Sx i -Rx|,i=1,……N (4) Among them, (Sx i , Sy i ) is the spatial position of the microseismic event, (Rx, Ry) is the coordinate of the detector, L Si ,θ Si are the horizontal monitoring distance and the geophone transmission angle of microseismic events, respectively; (22), define the independent variable function g(θ) of the combination of the detector transmission angle and the horizontal monitoring distance Si ,L Si ): Calculate the function g(θ corresponding to N microseismic events Si ,L Si ) (23) Construct the S-wave linear energy compensation function fSEN1(θ Si ,L Si ): fSEN1(θ Si ,L Si )=C1*g(θ Si ,L Si )+C2,i=N*+1,……N(7) Among them, when N * +1≤i≤N, L Si >1000, C1, C2 are the compensation function fSEN1(θ Si ,L Si )’s linear coefficient; (24), using the monitoring distance from near to far corresponding to the sequence number N * +1 to N microseismic events S wave signal energy SEN i , forming an overdetermined system of equations: (25), solve the overdetermined system of equations and obtain the values ​​of coefficients C1 and C2; (26) Construct the S-wave nonlinear energy compensation function fSEN2(θ Si ,L Si ): fSEN2(θ Si ,L Si )=D1*(g(θ Si ,L Si )) 2 +D2*g(θ Si ,L Si )+D3,i=1,……N*(9) Among them, when i≤N * When L Si ≤1000, D1, D2, D3 are the compensation function fSEN2(θ Si ,L Si )’s nonlinear coefficient; (27), using the monitoring distance from near to far corresponding to the sequence number 1 to N * The S-wave signal energy SEN of the microseismic event i , forming an overdetermined system of equations: (28), solve the overdetermined system of equations and obtain the values ​​of coefficients D1, D2, and D3.

8. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores at least one computer-executable program, and when the at least one program is executed by the computer, the computer executes the steps in the S-wave segmented energy compensation method for microseismic monitoring of horizontal well fracturing according to any one of claims 1 to 6.

Images