Horizontal well fracturing micro-seismic monitoring P-wave nonlinear energy compensation method and system

Abstract

The invention provides a P-wave nonlinear energy compensation method and system for horizontal well fracturing microseism monitoring, and belongs to the field of microseism signal processing. According to the method, micro-seismic signals monitored in the middle and long distance are corrected to be monitored in the short distance, so that weak signal energy is enhanced. According to the method, a nonlinear rule between energy of each perforation and a transmission angle and a nonlinear rule between energy of each perforation and a horizontal distance in the horizontal direction are searched by utilizing forward micro-seismic P-wave signal energy, a known perforation position and a detector position, a micro-seismic P-wave nonlinear energy compensation function taking the transmission angle as an independent variable is constructed, and a micro-seismic P-wave nonlinear energy compensation function can be obtained only by inputting a micro-seismic event space position subsequently. The P-wave energy compensation factor of the event can be calculated, and then P-wave energy compensation for horizontal well fracturing microseism monitoring is achieved. After the method is used, the P-wave energy is well compensated, the farther the distance is, the better the compensation effect is, the calculation process is efficient, and the method has very high practicability.

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 nonlinear energy of P waves in microseismic monitoring of horizontal well fracturing. Background Art

[0002] During the process of artificial fractures caused by hydraulic fracturing and breaking of underground rocks, the detector will receive a single P wave or a single S wave or a P&S wave combination signal, which becomes microseismic monitoring data.

[0003] Microseismicity is mainly divided into ground microseismicity and borehole microseismicity. The characteristics of ground microseismicity are: conventional ground geophones are used to collect microseismic signals, and its collection method is similar to borehole seismic VSP. There are many ground microseismic geophones with various arrangements and distributions, which can fully collect microseismic signals. However, due to the long distance between the underground source and the geophone, the received microseismic signals are weak, so the ground microseismic positioning is relatively stable but not accurate enough. The characteristics of borehole microseismicity are: the three-component geophone downhole is placed in the observation well section to receive the full wavefield signal of the microseismicity. 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. However, due to the limited number of microseismic geophones in the well (generally 12 to 32 three-component borehole geophones), the monitoring range is small, so the borehole microseismic positioning method is prone to unstable microseismic positioning results.

[0004] In the process of microseismic monitoring in wells such as hydraulic fracturing, carbon dioxide fracturing, and gas storage injection, the types of microseismic events are complex, and the signal energy is unequal, and even the magnitude is relatively low. A large number of weak events will appear, and these events are the true reflection of the development of cracks during the fracturing process. At the same time, it is affected by formation factors, the environment around the well, and various noises, making it more difficult to pick up the first arrival time of the weak signal of the microseismic event. Therefore, weak signal processing in microseismic monitoring data records becomes the key, which is conducive to event identification and positioning.

[0005] In the processing of microseismic weak signals, there are many means to enhance signals, such as: various denoising methods, which improve the signal-to-noise ratio of data by removing random noise and linear interference, thereby enhancing the effective signal; various energy compensation methods, such as the spherical diffusion energy compensation method based on the propagation path, and energy compensation based on the surface or formation absorption factor.

[0006] However, due to the monitoring distance, the same source is excited at different horizontal well monitoring distances, and the signal energy received by the geophone is different, which affects the energy analysis and interpretation of microseismic events. This is a defect of microseismic monitoring in wells. The current method only compensates for the energy weakening caused by the signal in the propagation path, and cannot solve the disadvantage of systematic energy weakening in well monitoring. 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 nonlinear energy compensation of P waves in microseismic monitoring of horizontal well fracturing. By forward modeling data, a microseismic P-wave nonlinear energy compensation function is constructed. Subsequently, the microseismic P-wave signal can be enhanced by simply inputting the spatial position of the microseismic event, thereby compensating for the defect of weak microseismic signal caused by the long-distance observation method.

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

[0009] A first aspect of the present invention provides a method for nonlinear energy compensation of P waves in microseismic monitoring of horizontal well fracturing, which corrects microseismic signals monitored at medium and long distances to those monitored at short distances, thereby enhancing the energy of weak signals.

[0010] Preferably, the method comprises:

[0011] The first step is to obtain microseismic P-wave forward modeling data based on known acoustic logging data, perforation locations, and geophone locations;

[0012] The second step is to obtain the first arrival time of the microseismic perforation P wave and the signal energy of the microseismic P wave forward modeling data;

[0013] The third step is to construct the nonlinear energy compensation function of microseismic P waves;

[0014] The fourth step is to perform microseismic P-wave energy compensation for any microseismic event.

[0015] Preferably, the first step includes:

[0016] First, the interval velocity model is established based on the known sonic logging data;

[0017] Secondly, set the coordinates of one geophone position (Rx, Ry) and the coordinates of N perforations (Sx i , Sy i ), i = 1, ... N, i represents the i-th perforation, N is the number of perforations, from the first perforation to the N-th perforation, the distance between the perforation and the detector is from near to far;

[0018] Finally, the microseismic wave field forward modeling was carried out to obtain N microseismic P wave forward modeling data APi,j , j=1,...M, j represents the jth time sample point, and M is the number of time sample points.

[0019] Preferably, the second step includes:

[0020] The first arrival time of microseismic perforation P wave is obtained by using the energy ratio of long and short time windows;

[0021] Get the microseismic P wave forward data signal energy EN i .

[0022] Preferably, the operation of obtaining the first arrival time of the microseismic perforation P wave by using the energy ratio of the long and short time windows includes:

[0023] First, define the long time window L1 and the short time window L2, and calculate the energy ratio curve ER of the microseismic P wave forward modeling data. Pi,j :

[0024]

[0025] Among them, ER Pi,j is the energy ratio curve corresponding to the time sample point j of the i-th perforation in the forward modeling;

[0026] Secondly, define the threshold value K P , perform the following operations on the energy ratio curve of each microseismic P-wave forward modeling data in turn:

[0027] Search ER Pi,j The first time ER appeared Pi,j ≥K P The time sample T corresponding to KP ;

[0028] Define the first arrival search window size W P , search for the maximum peak value F in the time window P : F P =max{ER Pi,j} and j∈[T KP ,T KP +W P ];

[0029] Maximum peak F P The corresponding time sample t P That is, the first arrival time of the microseismic perforation P wave.

[0030] Preferably, the microseismic P wave forward modeling data signal energy EN is obtained i The operations include:

[0031] The signal energy EN of N microseismic P-wave forward modeling data is calculated using the following formula: i :

[0032]

[0033] Preferably, the operation of the third step includes:

[0034] First, by setting the known perforation position coordinates (Sx i , Sy i ) and a detector position coordinate (Rx, Ry), calculate the perforation relative to the detector transmission angle θ i :

[0035]

[0036] Then, construct the transmission angle θ i The microseismic P-wave nonlinear energy compensation function fPEN(θ i ):

[0037]

[0038] Among them, L i =|Sx i -Rx| is the horizontal distance of the ith perforation relative to the detector;

[0039] Finally, establish and solve the overdetermined matrix to obtain coefficients C1, C2, and C3:

[0040]

[0041] Preferably, the fourth step includes:

[0042] First, assume that the spatial coordinates of any event are (event X , event Y ), and use the following formula to calculate the transmission angle of the event relative to the detector:

[0043]

[0044] Secondly, substitute the transmission angle of the event relative to the detector into the microseismic P wave nonlinear energy compensation function to obtain fPEN(θ * );

[0045] Again, construct the event energy compensation factor:

[0046]

[0047] Finally, multiply the P wave amplitude of the event by The compensated microseismic P-wave signal value is obtained.

[0048] A second aspect of the present invention provides a P-wave nonlinear energy compensation system for microseismic monitoring of horizontal well fracturing, the system comprising:

[0049] A forward modeling unit is used to obtain microseismic P-wave forward modeling data based on known acoustic logging data, perforation locations, and geophone locations;

[0050] The first arrival time and energy acquisition unit is connected to the forward modeling unit and is used to obtain the first arrival time of the microseismic perforation P wave and the energy of the microseismic P wave forward modeling data signal;

[0051] A compensation function construction unit, connected with the first arrival time and energy acquisition unit, is used to construct a microseismic P-wave nonlinear energy compensation function;

[0052] The energy compensation unit is connected to the compensation function construction unit and is used to perform microseismic P wave energy compensation for any given microseismic event.

[0053] 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 nonlinear energy compensation of P waves in microseismic monitoring of horizontal well fracturing according to the present invention.

[0054] Compared with the prior art, the beneficial effects of the present invention are as follows: in the weak signal processing of microseismic events, the present invention uses the forward modeled microseismic P-wave signal energy, known perforation positions, and detector positions to find the nonlinear law between the energy of each perforation in the horizontal direction and the transmission angle and horizontal distance, and constructs a microseismic P-wave nonlinear energy compensation function with the transmission angle as the independent variable. Subsequently, the P-wave energy compensation factor of the event can be calculated by simply inputting the spatial position of the microseismic event, thereby realizing P-wave energy compensation for microseismic monitoring of horizontal well fracturing.

[0055] Through data testing, the P wave energy is well compensated by the present invention, and the longer the distance, the better the compensation effect. In particular, when the horizontal distance is above 1,500 meters, the P wave signal changes from weak to strong, proving that the present invention has a good compensation effect on microseismic P wave energy. Moreover, the calculation process is efficient and has strong practicality. BRIEF DESCRIPTION OF THE DRAWINGS

[0056] Figure 1 It is a flowchart of the steps of the method of the present invention;

[0057] Figure 2 It is the microseismic forward velocity model and the perforation and geophone positions;

[0058] Figure 3 It is the data profile of the original microseismic P-wave forward modeling data;

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

[0060] Figure 5 yes Figure 3 Microseismic P-wave data profile after nonlinear energy compensation processing. DETAILED DESCRIPTION

[0061] The present invention is further described in detail below in conjunction with the accompanying drawings:

[0062] Since the observation method leads to weak signals, the present invention proposes a method for compensating microseismic signal energy for microseismic events in horizontal well monitoring. By counting the signal energies received by detectors at different monitoring distances, a compensation function related to the monitoring distance is constructed to carry out microseismic energy compensation, correct the energy to the signal enhancement performed by reasonable observation, and minimize the influence of the monitoring distance on the weakening of the energy of microseismic events, so that the microseismic signals received at different locations of the same level of source points are roughly consistent, providing reliable guarantee for subsequent magnitude statistics of microseismic events.

[0063] Different from the conventional energy compensation method, the present invention uses energy compensation to correct the microseismic signals of medium and long-distance monitoring to short-distance monitoring, so that the energy of some weak signals is enhanced.

[0064] Embodiments of the present invention are as follows:

[0065] Embodiment 1:

[0066] like Figure 1 As shown, the present invention provides a method for P-wave nonlinear energy compensation for microseismic monitoring of horizontal well fracturing, the method comprising:

[0067] The first step is to obtain the elastic wave forward modeling perforation P-wave signal monitored by the geophone, i.e., the microseismic P-wave forward modeling data, based on the known acoustic logging data, perforation location, and geophone location, as follows:

[0068] Establishing interval velocity model based on known sonic logging data;

[0069] Set the position coordinates (Rx, Ry) of a geophone and the position coordinates (Sx) of N horizontal wells with known perforations from near to far from the geophone (i.e., among the N perforations on the horizontal well, from the first perforation to the Nth perforation, the distances between the perforations and the geophone are from near to far). i , Sy i ), i=1,...N, i represents the ith perforation, and these N+1 coordinates are all known.

[0070] The existing common elastic wave equation is used to carry out microseismic wave field forward modeling, and N microseismic P wave forward modeling data AP are obtained. i,j (i.e., P-wave amplitude value), i=1,...N, j=1,...M, where i represents the i-th perforation, j represents the j-th time sample point, N is the number of perforations, M is the number of time sample points, and each perforation P-wave data has a total of M time sample points.

[0071] The second step is to obtain the first arrival time of the microseismic perforation P wave and the signal energy EN of the microseismic P wave forward modeling data. i , as follows:

[0072] The first arrival time of the microseismic perforation P wave is obtained by using the energy ratio of the long and short time windows. A certain time window is set to further obtain the microseismic P wave forward modeling data signal energy EN i :

[0073] First, define the long time window L1 and the short time window L2, where L1 is the time window length of the long time window, L2 is the time window length of the short time window, and L1>L2. Calculate the P wave data AP obtained in the first step i The long and short time window energy ratio curve:

[0074]

[0075] Among them, AP i,j , ER Pi,j They are respectively the P-wave amplitude value and the long-short time window energy ratio curve corresponding to the time sample point j of the i-th perforation in the forward modeling.

[0076] Define the threshold value K P , respectively perform signal recognition on N forward modeled microseismic P wave data: search for the long-short time window energy ratio curve ER Pi,j The first time ER appeared Pi,j ≥K P The time sample T corresponding to KP , that is, the following conditions are met: j = T KP And ER Pi,j ≥K P . Compare the energy ratio data with the threshold value. If ER appears for the first time Pi,j ≥K P When , it indicates that a P wave signal of a microseismic event occurs near that time point.

[0077] At the same time, define the first arrival search window size W P (usually the length of the wavelet), find the maximum peak value F in the time window P : F P =max{ER Pi,j} and j∈[T KP ,T KP+W P ] The maximum peak value F P The corresponding time sample t P Output, time sample t P That is, the first arrival time of the microseismic perforation P wave.

[0078] On this basis, we further calculate the signal energy EN of N microseismic P-wave forward modeling data. i :

[0079]

[0080] The third step is to calculate the transmission angle θ of the perforation relative to the detector using the known perforation position and detector position. i , combined with the second step to obtain the microseismic P wave forward data signal energy EN i , construct the microseismic P wave nonlinear energy compensation function fPEN(θ i ), as follows:

[0081] Construct the microseismic P wave energy compensation function fPEN(θ i ):

[0082] First, by setting the known perforation position coordinates (Sx i , Sy i ) and a detector position coordinate (Rx, Ry), calculate the perforation relative to the detector transmission angle θ i :

[0083]

[0084] Then, construct the transmission angle θ i The microseismic P-wave nonlinear energy compensation function fPEN(θ i ):

[0085]

[0086] Among them, L i =|Sx i -Rx| is the horizontal distance of the ith perforation relative to the detector.

[0087] By combining the microseismic P-wave forward data signal energy EN calculated in formula (2) i , forming N equations and establishing an overdetermined matrix:

[0088]

[0089] The present invention uses the existing singular value decomposition (SVD) method to solve the matrix of the above formula (5) and calculate the coefficients C1, C2, and C3, thereby realizing the P wave energy compensation function fPEN (θi ) construction.

[0090] The fourth step is to implement microseismic P-wave energy compensation for any microseismic event:

[0091] Given any microseismic event, according to the event position (microseismic events include perforation, that is, perforation is a microseismic event with known design coordinates. Other microseismic events need to be positioned to obtain the event position, but the positioning process can be done by existing methods and is not within the scope of protection of the present invention), calculate the event relative to the geophone transmission angle and substitute it into the third step microseismic P wave nonlinear energy compensation function fPEN (θ i ), construct the event energy compensation factor Q * , multiply the P wave amplitude value of the event by Finally, the compensated P wave signal value is obtained, as follows:

[0092] Assume that the spatial coordinates of any event are (event X , event Y ), calculate the transmission angle of the event relative to the detector:

[0093]

[0094] Substituting into formula (4) we can obtain the P wave energy compensation function fPEN(θ * );

[0095] Construct the energy compensation factor for this event:

[0096]

[0097] Multiply the P wave amplitude of the event by The compensated microseismic P-wave signal value can be obtained.

[0098] The method of the present invention takes N perforation P-wave energies as known samples and constructs a P-wave energy compensation function fPEN(θ i ), then input any microseismic event (including perforation), calculate the corresponding P-wave energy compensation function, and compare it with the strongest P-wave energy EN1 at the nearest monitoring distance Q * , multiplied by the event P wave amplitude value That is, the energy of microseismic events at any monitoring distance is corrected to the size of EN1, thereby achieving energy compensation).

[0099] The present invention utilizes the forward modeled microseismic P-wave signal energy, known perforation locations, and detector locations to find the nonlinear law between each perforation energy, the transmission angle, and the horizontal distance in the horizontal direction, and constructs a microseismic P-wave nonlinear energy compensation function with the transmission angle as the independent variable. Subsequently, the P-wave energy compensation factor of the event can be calculated by simply inputting the spatial position of the microseismic event, thereby realizing P-wave energy compensation for microseismic monitoring of horizontal well fracturing.

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

[0101] Embodiment 2:

[0102] First, according to Figure 1 The process shown in the figure performs the following operations: a velocity model is established from known velocity logs, such as Figure 2 As shown, and set 1 detector position ( Figure 2 white triangles) and 111 perforation locations ( Figure 2 The horizontal distance from the horizontal well perforation to the geophone is simulated to be 150 to 2350 meters, and the distance between perforations is 20 meters. Then, the microseismic P wave elastic wave forward modeling is carried out using the differential wave equation to obtain the original microseismic P wave data, such as Figure 3 By combining the microseismic P wave forward modeling data signal energy calculated in formula (2), 111 equations are formed, an overdetermined matrix is ​​established, and the coefficients in formula (4) of the present invention are calculated using singular value decomposition SVD.

[0103] Secondly, the previous calculation process uses the energy ratio of the long and short time windows (such as formula (1)) and sets a wave field size time window to obtain the microseismic P wave energy (such as formula (2)). The microseismic P wave nonlinear energy compensation function of the present invention is constructed (such as formula (4)), and C1, C2, and C3 are calculated, thereby realizing the construction of the microseismic P wave nonlinear energy compensation function. Figure 4 As shown, the dots in the figure are the intersection diagram of the energy and independent variables of the microseismic P-wave energy compensation function, and the black line is the microseismic P-wave nonlinear energy compensation function curve.

[0104] Finally, according to the microseismic P-wave nonlinear energy compensation function calculated in the previous step, the energy compensation of 111 perforation data is performed, as follows: Figure 5 Compare Figure 3 From the original P data, it can be seen that the perforation P wave energy is well compensated. The farther the distance, the better the compensation effect. In particular, when the horizontal distance is above 1500 meters, the P wave signal changes from weak to strong, which verifies that the present invention has a good compensation effect on microseismic P wave energy.

[0105] The present invention also provides a P-wave nonlinear energy compensation system for horizontal well fracturing microseismic monitoring, and the embodiments of the system are as follows:

[0106] Embodiment three:

[0107] The system comprises:

[0108] A forward modeling unit is used to obtain microseismic P-wave forward modeling data based on known acoustic logging data, perforation locations, and geophone locations;

[0109] The first arrival time and energy acquisition unit is connected to the forward modeling unit and is used to obtain the first arrival time of the microseismic perforation P wave and the energy of the microseismic P wave forward modeling data signal;

[0110] A compensation function construction unit, connected with the first arrival time and energy acquisition unit, is used to construct a microseismic P-wave nonlinear energy compensation function;

[0111] The energy compensation unit is connected to the compensation function construction unit and is used to perform microseismic P wave energy compensation for any given microseismic event.

[0112] The present invention also 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 a P-wave nonlinear energy compensation method for horizontal well fracturing microseismic monitoring of the present invention.

[0113] 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 compensating nonlinear energy of P waves in microseismic monitoring of horizontal well fracturing, characterized by: The method corrects the microseismic signals monitored at medium and long distances to those monitored at short distances, thereby enhancing the energy of weak signals.

2. The method for P-wave nonlinear energy compensation in microseismic monitoring of horizontal well fracturing according to claim 1, characterized in that: The method comprises: The first step is to obtain microseismic P-wave forward modeling data based on known acoustic logging data, perforation locations, and geophone locations; The second step is to obtain the first arrival time of the microseismic perforation P wave and the signal energy of the microseismic P wave forward modeling data; The third step is to construct the nonlinear energy compensation function of microseismic P waves; The fourth step is to perform microseismic P-wave energy compensation for any microseismic event.

3. The method for P-wave nonlinear energy compensation in microseismic monitoring of horizontal well fracturing according to claim 2 is characterized by: The first step includes: First, the interval velocity model is established based on the known sonic logging data; Secondly, set the coordinates of one geophone position (Rx, Ry) and the coordinates of N perforations (Sx i , Sy i ), i = 1, ... N, i represents the i-th perforation, N is the number of perforations, from the first perforation to the N-th perforation, the distance between the perforation and the detector is from near to far; Finally, the microseismic wave field forward modeling was carried out to obtain N microseismic P wave forward modeling data AP i,j , j=1,...M, j represents the jth time sample point, and M is the number of time sample points.

4. The method for P-wave nonlinear energy compensation in microseismic monitoring of horizontal well fracturing according to claim 2, characterized in that: The second step includes: The first arrival time of microseismic perforation P wave is obtained by using the energy ratio of long and short time windows; Get the microseismic P wave forward data signal energy EN i .

5. The method for P-wave nonlinear energy compensation in microseismic monitoring of horizontal well fracturing according to claim 4 is characterized by: The operation of obtaining the first arrival time of the microseismic perforation P wave by using the energy ratio of the long and short time windows includes: First, define the long time window L1 and the short time window L2, and calculate the energy ratio curve ER of the microseismic P wave forward modeling data. Pi,j : Among them, ER Pi,j is the energy ratio curve corresponding to the time sample point j of the i-th perforation in the forward modeling; Next, define the threshold value K P , perform the following operations on the energy ratio curve of each microseismic P-wave forward modeling data in turn: Search ER Pi,j The first time ER appeared Pi,j ≥K P The time sample T corresponding to KP ; Define the first arrival search window size W P , search for the maximum peak value F in the time window P : F P =max{ER Pi,j } and j∈[T KP ,T KP +W P ]; Maximum peak F P The corresponding time sample t P That is, the first arrival time of the microseismic perforation P wave.

6. The method for P-wave nonlinear energy compensation in microseismic monitoring of horizontal well fracturing according to claim 5, characterized in that: Get the microseismic P wave forward data signal energy EN i The operations include: The signal energy EN of N microseismic P-wave forward modeling data is calculated using the following formula: i :

7. The method for P-wave nonlinear energy compensation in microseismic monitoring of horizontal well fracturing according to claim 2, characterized in that: The third step includes: First, by setting the known perforation position coordinates (Sx i , Sy i ) and a detector position coordinate (Rx, Ry), calculate the perforation relative to the detector transmission angle θ i : Then, construct the transmission angle θ i The microseismic P-wave nonlinear energy compensation function fPEN(θ i ): Among them, L i =|Sx i -Rx| is the horizontal distance of the ith perforation relative to the detector; Finally, establish and solve the overdetermined matrix to obtain coefficients C1, C2, and C3:

8. The method for P-wave nonlinear energy compensation in microseismic monitoring of horizontal well fracturing according to claim 2, characterized in that: The fourth step includes: First, assume that the spatial coordinates of any event are (event X , event Y ), and use the following formula to calculate the transmission angle of the event relative to the detector: Secondly, substitute the transmission angle of the event relative to the detector into the microseismic P wave nonlinear energy compensation function to obtain fPEN(θ * ); Again, construct the event energy compensation factor: Finally, multiply the P wave amplitude of the event by The compensated microseismic P-wave signal value is obtained.

9. A P-wave nonlinear energy compensation system for microseismic monitoring of horizontal well fracturing, characterized by: The system comprises: A forward modeling unit is used to obtain microseismic P-wave forward modeling data based on known acoustic logging data, perforation locations, and geophone locations; The first arrival time and energy acquisition unit is connected to the forward modeling unit and is used to obtain the first arrival time of the microseismic perforation P wave and the energy of the microseismic P wave forward modeling data signal; A compensation function construction unit, connected with the first arrival time and energy acquisition unit, is used to construct a microseismic P-wave nonlinear energy compensation function; The energy compensation unit is connected to the compensation function construction unit and is used to perform microseismic P wave energy compensation for any given microseismic event.

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 P-wave nonlinear energy compensation method for microseismic monitoring of horizontal well fracturing as described in any one of claims 1-8.

Images