Horizontal well fracturing microseism monitoring S-wave segmented energy compensation method and horizontal well fracturing microseism monitoring S-wave segmented energy compensation system

Abstract

The invention provides a horizontal well fracturing micro-seismic monitoring S-wave segmented energy compensation method and system, and belongs to the field of micro-seismic signal processing. The method comprises the steps of firstly constructing a microseismic event S-wave first arrival pickup function, obtaining S-wave signal energy, then constructing a linear energy compensation function and a non-linear energy compensation function of microseismic event S-waves in a segmented manner by using the position of a microseismic event and the position of a detector, and finally inputting the spatial position of any microseismic event. According to the horizontal monitoring distance, an S-wave signal energy compensation factor of the micro-seismic event is obtained, and horizontal well fracturing micro-seismic monitoring S-wave signal energy compensation is achieved through the S-wave signal energy compensation factor. According to the method, the problem of weak signals caused by the observation distance is solved, the S-wave signal energy is corrected to the reasonable observation distance in a segmented mode according to the law of the S-wave signal energy, and micro-seismic signals received by the same-level seismic source point at different positions are roughly consistent.

Description

Technical Field

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

[0002] In actual fracturing monitoring, the source of the earthquake generated by rock compression is often dominated by shear waves. Microseismic events will generate both P-wave and S-wave signals, but the energy of S-wave signals is much stronger than that of P-waves. Therefore, S-waves are more reliable in the study of the source of microseismic events. Usually, microseismic monitoring data refers to the single P-wave or single S-wave or P&S wave combination signal received by the detector during the process of artificial fractures caused by hydraulic fracturing and underground rock crushing.

[0003] Microseismic monitoring generally includes: ground microseismic monitoring and well microseismic monitoring. The three-component geophone is placed in the observation well section to receive the full wavefield signal of microseismic monitoring. Compared with ground microseismic monitoring, the data received in the well has a higher signal-to-noise ratio and a richer number and type of microseismic events. Compared with single-component monitoring of ground microseismic monitoring, three-component monitoring of well microseismic monitoring is more sensitive to S waves, that is, 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 wellbore geophones are set), the monitoring range is small, and the S waves of microseismic events received by the geophones vary greatly due to different monitoring distances, which makes the well microseismic positioning method prone to unstable microseismic interpretation results.

[0005] In the energy compensation processing of microseismic weak signals, the influence of microseismic monitoring distance in the well can be reduced through energy compensation. The current means of signal enhancement include: improving the signal-to-noise ratio of data by removing random noise and linear interference, 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 on the interface, if the transmission angle is smaller, the energy of the transmitted wave is stronger and the energy of the reflected wave is weaker. Conversely, if the transmission angle is larger, the energy of the transmitted wave is weaker and the energy of the reflected wave is stronger. For the microseismic observation method in horizontal well fracturing wells, the detector mainly receives the transmitted wave. The closer the monitoring distance and the smaller the transmission angle, the stronger the energy of the received S-wave signal. Conversely, the farther the monitoring distance and the larger the transmission angle, the weaker the energy of the corresponding received S-wave signal. Therefore, this phenomenon leads to the same energy level source being excited at different monitoring distances, and the energy of the S-wave signal received by the detector is quite different, which causes the microseismic misinterpretation. Summary of the invention

[0007] The purpose of the present invention is to solve the above-mentioned difficulties existing in the prior art and to provide a method and system for S-wave segmented energy compensation for microseismic monitoring of horizontal well fracturing, which can eliminate or partially eliminate the influence of monitoring distance on the energy of S-wave signals, so that the energy of S-wave signals excited by the same energy level source at different monitoring distances received by the detector is equivalent or close.

[0008] The present invention is achieved through the following technical solutions:

[0009] The first aspect of the present invention provides a method for segmented energy compensation of S waves in microseismic monitoring of horizontal well fracturing. The method first constructs a microseismic event S wave first arrival picking function to obtain the S wave signal energy, and then uses the microseismic event position and the detector position to segmentally construct a linear energy compensation function and a nonlinear energy compensation function of 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 microseismic monitoring of horizontal well fracturing.

[0010] Preferably, the method comprises:

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

[0012] (2) S-wave linear energy compensation function and S-wave nonlinear energy compensation function are constructed in sections;

[0013] (3) Energy compensation of the S wave of any microseismic event is performed using an S wave linear energy compensation function or an S wave nonlinear energy compensation function.

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

[0015] (11) Processing the microseismic data to obtain S-wave data and spatial positions of N microseismic events with monitoring distances from near to far;

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

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

[0018] (121), define a long time window L1 and a short time window L2, the time window length of L1 is greater than the time window length of L2, and construct the 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 point, N is the total amount of microseismic sampling data, M is the total amount of time sample points, and AS i,j , ER Si,j ,fS i (t j ) are respectively the S-wave amplitude value, 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;

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

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

[0024]

[0025] Preferably, the time window W S The size of is the time sample length corresponding to the length of a sub-wave.

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

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

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

[0029]

[0030] 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;

[0031] (22), define the independent variable function g(θ) of the combination of the detector transmission angle and the horizontal monitoring distance Si ,L Si ):

[0032]

[0033] Calculate the function g(θ corresponding to N microseismic events Si ,L Si )

[0034] (23) Construct the S-wave linear energy compensation function fSEN1(θ Si ,L Si ):

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

[0036] Among them, when N * +1≤i≤N, L Si >1000, C1, C2 are the compensation function fSEN1(θ Si ,L Si )’s linear coefficient;

[0037] (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:

[0038]

[0039] (25), solve the overdetermined system of equations and obtain the values ​​of coefficients C1 and C2;

[0040] (26) to construct the S-wave nonlinear energy compensation function fSEN2(θ Si ,L Si ):

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

[0042] Among them, when i≤N * When L Si ≤1000, D1, D2, D3 are the compensation function fSEN2(θ Si ,L Si )’s nonlinear coefficient;

[0043] (27), using the monitoring distance from near to far corresponding to the serial number 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 >1000 m, 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 ≤1000 m, the energy compensation factor SQ of the S-wave signal of the microseismic event is obtained by 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, picking up the S wave first arrival time and S wave signal energy;

[0056] A segmented construction unit, connected to the input and pickup unit, for segmentally constructing 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] According to a third aspect of the present invention, a computer-readable storage medium is provided, wherein 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 of a method for S-wave segmented energy compensation for microseismic monitoring of horizontal well fracturing according to the present invention.

[0059] Compared with the prior art, 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 points are generally consistent, providing a reliable guarantee for the subsequent magnitude statistics of microseismic events. BRIEF DESCRIPTION OF THE DRAWINGS

[0060] Figure 1 It is a flowchart 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 and independent variable of the microseismic S-wave nonlinear energy compensation function;

[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 in conjunction with the accompanying drawings:

[0067] The present invention provides a method for segmented energy compensation of S waves in microseismic monitoring of horizontal well fracturing. The method first constructs a microseismic event S wave first arrival picking function to obtain S wave signal energy, then uses known microseismic event positions and detector positions to segmentally construct linear energy compensation functions and nonlinear energy compensation functions of the microseismic event S waves, 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 achieve S wave signal energy compensation for microseismic monitoring of horizontal well fracturing.

[0068] like Figure 1 As shown, the S-wave segmented energy compensation method for horizontal well fracturing microseismic monitoring of the present invention comprises the following three steps:

[0069] (1) Input the original microseismic event S-wave signal, pick up the S-wave first arrival time and S-wave signal energy;

[0070] (2) S-wave linear energy compensation function and S-wave nonlinear energy compensation function are constructed in sections;

[0071] (3) Energy compensation of the S wave of any microseismic event is performed using an S wave linear energy compensation function or an S wave nonlinear energy compensation function.

[0072] The embodiments of the method of the present invention are as follows:

[0073] Embodiment 1:

[0074] The operation of (1) includes:

[0075] (11) Using existing methods to process microseismic data, obtain S-wave data and spatial positions of N microseismic events with monitoring distances from near to far, and use these as input data;

[0076] (12) Construct the S wave first arrival picking function, pick the S wave first arrival time and S wave signal energy SEN i .

[0077] The operation of (12) includes:

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

[0079]

[0080]

[0081] Where i represents the i-th microseismic sampling event, j represents the j-th time sample point, N is the total amount of microseismic sampling data, M is the total amount of time sample points, and ASi,j , ER Si,j ,fS i (t j ) are respectively the S-wave amplitude value, 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;

[0082] (122), the S-wave first arrival time is obtained by searching for the maximum value of the S-wave first arrival picking function: when t j =t * When fS i (t j ) has the largest value, at this time t * is the S-wave first arrival time of the i-th microseismic sampling event.

[0083] (123), define a time window W S (The size is generally the length of a sub-wavelength time sample point), and the S-wave signal energy SEN of the microseismic event is obtained using the following formula i :

[0084]

[0085] Embodiment 2:

[0086] (2) The operation includes:

[0087] The spatial positions of N microseismic events are divided into two groups in sequence: one with a horizontal monitoring distance greater than 1000 meters and one with a horizontal monitoring distance less than or equal to 1000 meters. Si Monitoring distance L Si As independent variables, construct the S wave linear energy compensation function fSEN1(θ Si ,L Si ), S-wave nonlinear energy compensation function fSEN2(θ Si ,L Si ). Specifically include:

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

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

[0090]

[0091] Among them, (Sx i , Sy i ) is the known spatial position of the microseismic event, (Rx, Ry) is the known coordinates of the geophone, L Si ,θ Siare the horizontal monitoring distance of microseismic events and the transmission angle of the geophone, respectively.

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

[0093]

[0094] Calculate the function g(θ corresponding to N microseismic events Si ,L Si ) is the value of .

[0095] (23) Construct the S-wave linear energy compensation function fSEN1(θ Si ,L Si ):

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

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

[0098] (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:

[0099]

[0100] (25), the existing singular value decomposition (SVD) method is used to solve the above overdetermined equations (the process is omitted), and the coefficients C1 and C2 are calculated to realize the microseismic event S-wave linear energy compensation function fSEN1(θ Si ,L Si ) construction.

[0101] (26) to construct the S-wave nonlinear energy compensation function fSEN2(θ Si ,L Si ):

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

[0103] Among them, when i≤N * When L Si ≤1000, D1, D2, D3 are the compensation function fSEN2(θ Si ,L Si ) is the nonlinear coefficient of .

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

[0105]

[0106] (28) The existing singular value decomposition (SVD) method is used to solve the above overdetermined equations (the process is omitted) and the coefficients D1, D2, and D3 are calculated, that is, the nonlinear energy compensation function fSEN2(θ Si ,L Si ) construction.

[0107] Embodiment three:

[0108] The operation of (3) includes:

[0109] After long-term seismic data processing, it is found that the energy of the S-wave signal shows different rules with the change of monitoring distance and transmission angle. When the monitoring distance is greater than 1000 meters, it shows a linear rule. Conversely, when the monitoring distance is less than or equal to 1000 meters, it shows a nonlinear rule. Therefore, the present invention performs segmented processing on the compensation function, that is, divides it into linear and nonlinear parts to construct compensation functions respectively, so that the energy compensation of S waves at different monitoring distances can be more targeted.

[0110] Specifically, input any microseismic event and calculate the detector transmission angle of the event Monitoring distance And compare With a size of 1000 m, the corresponding S-wave linear or nonlinear energy compensation function is then substituted to construct the microseismic event S-wave signal energy compensation factor SQ * , and the S wave amplitude of the event Perform point multiplication to finally obtain the S-wave signal of the microseismic event after energy compensation, which specifically includes:

[0111] (31), assuming that the spatial location of the microseismic event is (event X , event Y ), using the above formulas (4) to (6), the monitoring distance L of the microseismic event is calculated respectively. event , detector transmission angle θ event With the independent variable function g(θ event , L event ):

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

[0113]

[0114]

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

[0116]

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

[0118]

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

[0120] Embodiment 4:

[0121] The present invention also provides a horizontal well fracturing microseismic monitoring S-wave segmented energy compensation system, the system comprising:

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

[0123] A segmented construction unit, connected to the input and pickup unit, for segmentally constructing 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 by 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 the geophones range from 150 to 3130 meters, with 20 meters between events.

[0127] Using the differential wave equation, we simulated the microseismic three-component data in the well. After processing the microseismic three-component vector signal, we can obtain 150 microseismic S-wave data, such as Figure 3 As shown in Figure 2, the S-wave signal energy 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 detectors of 150 known monitoring distances from near to far are calculated. Then, according to whether the horizontal monitoring distance is greater than 1000 meters, they are respectively substituted into formulas (7) and (9) to establish two sets of overdetermined matrices (such as formulas (8) and (10)), and the singular value decomposition SVD is used to solve them to obtain the corresponding linear coefficients C1, C2 and nonlinear coefficients D1, D2, D3, and then the energy compensation function of the microseismic event S-wave signal 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 150 original microseismic event S-wave data, as shown in Figure 6 Compare Figure 3 From the original S data in the figure, 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 invention, unless otherwise specified, the terms "upper", "lower", "left", "right", "inside", "outside", etc. indicate directions or positional relationships based on the directions or positional relationships shown in the accompanying drawings. They are only for the convenience of describing the present invention and simplifying the description, and do not indicate or imply that the device or element referred to must have a specific direction, be constructed and operated in a specific direction. Therefore, they cannot be understood as limitations on the present invention.

[0131] The above technical solution is only one implementation mode of the present invention. For those skilled in the art, it is easy to make various types of improvements or modifications based on the principles disclosed in the present invention, and it is not limited to the technical solution described in the above specific embodiments of the present invention. Therefore, the above description is only preferred and does not have a restrictive meaning.

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 detector position to construct the linear energy compensation function and the nonlinear energy compensation function of the microseismic event S-wave in sections, 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 S-wave signal energy compensation in horizontal well fracturing microseismic monitoring.

2. The S-wave segmented energy compensation method for horizontal well fracturing microseismic monitoring according to claim 1 is characterized by: 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 sections; (3) Energy compensation of the S wave of any microseismic event is performed using an S wave linear energy compensation function or an S wave nonlinear energy compensation function.

3. The S-wave segmented energy compensation method for horizontal well fracturing microseismic monitoring according to claim 2 is characterized by: The operation of (1) includes: (11) Processing the microseismic data to obtain 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 to pick up the S-wave first arrival time and S-wave signal energy.

4. The S-wave segmented energy compensation method for horizontal well fracturing microseismic monitoring according to claim 3 is characterized by: The operation of (12) includes: (121), define a long time window L1 and a short time window L2, the time window length of L1 is greater than the time window length of L2, and construct the 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 sample point, N is the total amount of microseismic sampling data, M is the total amount of time sample points, and AS i,j , ER Si,j ,fS i (t j ) are respectively the S-wave amplitude value, 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 :

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

6. The S-wave segmented energy compensation method for horizontal well fracturing microseismic monitoring according to claim 2 is characterized by: (2) The operation 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) to construct the S-wave nonlinear energy compensation function fSEN2(θ Si ,L Si ): <h2 style=";text-align:left;direction:ltr">fSEN2(θ<h2 style=";text-align:left;direction:ltr"> Si <h2 style=";text-align:left;direction:ltr"> ,L<h2 style=";text-align:left;direction:ltr"> Si <h2 style=";text-align:left;direction:ltr"> )=D1*(g(θ<h2 style=";text-align:left;direction:ltr"> Si <h2 style=";text-align:left;direction:ltr"> ,L<h2 style=";text-align:left;direction:ltr"> Si <h2 style=";text-align:left;direction:ltr"> ))<h2 style=";text-align:left;direction:ltr"> 2 <h2 style=";text-align:left;direction:ltr"> +D2*g(θ<h2 style=";text-align:left;direction:ltr"> Si <h2 style=";text-align:left;direction:ltr"> ,L<h2 style=";text-align:left;direction:ltr"> Si <h2 style=";text-align:left;direction:ltr"> )+D3,(i=1,……N<h2 style=";text-align:left;direction:ltr"> * <h2 style=";text-align:left;direction:ltr"> ) (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 serial 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.

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

8. The S-wave segmented energy compensation method for horizontal well fracturing microseismic monitoring according to claim 2 is characterized by: The operation of (3) includes: (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 >1000 m, 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 ≤1000 m, the energy compensation factor SQ of the S-wave signal of the microseismic event is obtained by 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.

9. 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, picking up the S wave first arrival time and S wave signal energy; A segmented construction unit, connected to the input and pickup unit, for segmentally constructing an S-wave linear energy compensation function and an S-wave nonlinear energy compensation function; 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.

10. 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 as described in any one of claims 1-8.

Images