Cracking crack network statistical method and device based on minimum bounding rectangle envelope

Abstract

The invention relates to the technical field of petroleum geophysics, and particularly discloses a fracturing fracture network statistical method and device based on a minimum enclosing rectangle envelope, and the method comprises the steps: obtaining the height of a fracturing fracture network and the enclosing rectangle area of each angle in an azimuth range based on the position of a microseismic point set, the position of a perforation point and the azimuth range; obtaining an angle with the minimum rectangular area based on the circumscribed rectangular area of each angle in the azimuth range; and obtaining the length and width of the fracture net based on the angle with the minimum rectangular area. According to the fracturing fracture network statistical method based on the minimum enclosing rectangle envelope, the minimum enclosing rectangle envelope is obtained through the position information of the microseism event point set, the geometrical information of the fracturing fracture network is obtained, the fracturing fracture network statistical method based on the minimum enclosing rectangle envelope is formed, and automatic statistics of the geometrical information of the fracturing fracture network is achieved.

Description

Technical Field

[0001] The invention relates to the technical field of petroleum geophysics, and in particular to a method and device for calculating the statistics of a hydraulic fracture network based on a minimum circumscribed rectangular envelope. Background Art

[0002] In recent years, the development of unconventional oil and gas reservoirs, especially shale gas reservoirs, has been included in the national energy strategy. Hydraulic fracturing is a necessary way to commercially develop shale oil and gas. It connects natural fractures through fracturing to form artificial fracture networks. Microseismic monitoring technology is an important way to monitor hydraulic fracturing. It monitors the tiny earthquake events generated by hydraulic fracturing and then analyzes the fracture of underground rocks.

[0003] Microseismic monitoring results are generally based on spatial position, energy intensity, length, width, height and orientation of the fracture network, and further calculated by reservoir transformation volume. In terms of the statistics of fracture network length, width, height and orientation, it is currently mainly carried out through manual measurement. On the one hand, manual measurement introduces human errors, and on the other hand, there are also deficiencies in automation and efficiency.

[0004] Based on this technical background, the present invention studies a method and device for calculating the statistics of a fracture network based on a minimum circumscribed rectangular envelope. Summary of the invention

[0005] In view of the shortcomings of the prior art, the present invention provides a method and device for statistics of a hydraulic fracturing network based on a minimum circumscribed rectangular envelope. The method obtains the minimum circumscribed rectangular envelope of a microseismic event point set through its position information, obtains the geometric information of the hydraulic fracturing network, and forms a method for statistics of a hydraulic fracturing network based on the minimum circumscribed rectangular envelope, thereby realizing automatic statistics of the geometric information of the hydraulic fracturing network.

[0006] In order to achieve the above object, a first aspect of the present invention provides a method for calculating the statistics of a hydraulic fracture network based on a minimum circumscribed rectangular envelope, comprising:

[0007] Based on the microseismic point set position, perforation point position and azimuth range, the height of the fracture network and the circumscribed rectangular area of ​​each angle within the azimuth range are obtained respectively;

[0008] Based on the area of ​​the circumscribed rectangle of each angle within the azimuth range, the angle with the smallest rectangular area is obtained;

[0009] The length and width of the fracture network are obtained based on the angle with the smallest rectangular area.

[0010] A second aspect of the present invention provides a device for calculating the statistics of a hydraulic fracture network based on a minimum circumscribed rectangular envelope, comprising:

[0011] The area calculation module is used to obtain the circumscribed rectangular area of ​​each angle within the height and azimuth range of the fracture network based on the microseismic point set position, the perforation point position, and the azimuth range;

[0012] An angle obtaining module, used for obtaining the angle with the smallest rectangular area based on the circumscribed rectangular area of ​​each angle within the azimuth range;

[0013] The length and width obtaining module is used to obtain the length and width of the fracture network based on the angle with the smallest rectangular area.

[0014] A third aspect of the present invention provides an electronic device, the electronic device comprising:

[0015] A memory storing executable instructions;

[0016] A processor runs the executable instructions in the memory to implement the method for calculating the statistics of the fracture network based on the minimum circumscribed rectangular envelope described in the first aspect.

[0017] A fourth aspect of the present invention provides a computer-readable storage medium storing a computer program, which, when executed by a processor, implements the method for calculating the statistics of a fracture network based on a minimum circumscribed rectangular envelope as described in the first aspect.

[0018] The beneficial effects of the present invention include:

[0019] (1) The present invention proposes a method for statistically analyzing a hydraulic fracture network based on a minimum circumscribed rectangular envelope, which obtains the minimum circumscribed rectangular envelope of a microseismic event point set through its position information, and obtains the geometric information of the hydraulic fracture network, thereby forming a method for statistically analyzing a hydraulic fracture network based on a minimum circumscribed rectangular envelope, and realizing automatic statistics of the geometric information of the hydraulic fracture network.

[0020] (2) The present invention proposes a method for calculating the statistics of a hydraulic fracturing network based on the minimum circumscribed rectangle envelope. The method obtains the minimum circumscribed rectangle by traversing. The method has good stability and practicability, and the method has a small amount of calculation, which is suitable for real-time statistics of hydraulic fracturing network information.

[0021] Other features and advantages of the present invention will be described in detail in the following detailed description. BRIEF DESCRIPTION OF THE DRAWINGS

[0022] The above and other objects, features and advantages of the present invention will become more apparent through a more detailed description of exemplary embodiments of the present invention with reference to the accompanying drawings.

[0023] Figure 1 This is a schematic flow chart of the method for calculating the statistics of a hydraulic fracture network based on the minimum circumscribed rectangular envelope proposed in the present invention.

[0024] Figure 2A schematic flow chart of a specific implementation of the method for calculating the statistics of a fracture network based on the minimum circumscribed rectangular envelope proposed in the present invention.

[0025] Figure 3 The microseismic point set and perforation point positions in a specific implementation of the fracture network statistical method based on the minimum circumscribed rectangular envelope proposed by the present invention.

[0026] Figure 4 This is a schematic diagram of fracture network information statistics in a specific implementation of the fracture network statistics method based on the minimum circumscribed rectangular envelope proposed in the present invention.

[0027] Figure 5 This is a microseismic point set and perforation point positions in another specific implementation of the fracture network statistics method based on the minimum circumscribed rectangular envelope proposed by the present invention.

[0028] Figure 6 This is a schematic diagram of fracture network information statistics in another specific implementation of the fracture network statistics method based on the minimum circumscribed rectangular envelope proposed in the present invention. DETAILED DESCRIPTION

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

[0030] The present invention provides a method for calculating the statistics of a hydraulic fracture network based on a minimum circumscribed rectangular envelope, such as Figure 1 As shown, including:

[0031] Based on the microseismic point set position, perforation point position and azimuth range, the height of the fracture network and the circumscribed rectangular area of ​​each angle within the azimuth range are obtained respectively;

[0032] Based on the area of ​​the circumscribed rectangle of each angle within the azimuth range, the angle with the smallest rectangular area is obtained;

[0033] The length and width of the fracture network are obtained based on the angle with the smallest rectangular area.

[0034] In the present invention, the minimum circumscribed rectangular envelope of a microseismic event point set is obtained through its position information, and the geometric information of the fracture network is obtained, thereby forming a fracture network statistical method based on the minimum circumscribed rectangular envelope, and realizing automatic statistics of the geometric information of the fracture network.

[0035] According to the present invention, the calculation formula of the height of the fracture network is:

[0036] H=maxz i -minz i ;

[0037] Among them, zi is the height coordinate of the i-th point in the new microseismic point set formed by combining the microseismic point set position and the perforation point position, and H is the height of the fracture network.

[0038] According to the present invention, the calculation formula for the area of ​​the circumscribed rectangle of each angle within the azimuth range is:

[0039] A j =(maxL i -minL i ) × (maxW i -minW i );

[0040] Among them, Li is the position projection of the i-th point in the new microseismic point set on the axial vector, Wi is the position projection of the i-th point in the new microseismic point set on the lateral vector, and Aj is the area of ​​the circumscribed rectangle of each angle within the azimuth range.

[0041] According to the present invention, the calculation formula for the position projection of the i-th point in the new microseismic point set on the axial vector is:

[0042]

[0043] The calculation formula for the position projection of the i-th point in the new microseismic point set on the lateral vector is:

[0044]

[0045] Among them, x i ,y i are the horizontal and vertical coordinates of the i-th point in the new microseismic point set, l j is the axial vector at the jth azimuth angle within the azimuth range, w j is the lateral vector at the jth azimuth angle within the azimuth range.

[0046] Preferably, the expression of the axial vector at the jth azimuth angle within the azimuth range is:

[0047] l i =(cosθ j , sinθ j );

[0048] The expression of the lateral vector at the jth azimuth angle within the azimuth range is:

[0049] w j =(sinθ j , cosθ j );

[0050] Among them, θ jIn the azimuth range (θ 1 ,θ 2 ) at any angle.

[0051] According to the present invention, the calculation formula for the angle with the smallest rectangular area is:

[0052]

[0053] The angles with the smallest rectangular area based on the circumscribed rectangular area of ​​each angle within the azimuth range include: 1 ,θ 2 ) traverses to find the area of ​​the rectangle and obtains the angle θ corresponding to the area of ​​the minimum circumscribed rectangle.

[0054] Preferably, obtaining the length and width of the fracture network based on the angle with the smallest rectangular area includes:

[0055] Calculate the axial vector and lateral vector at the angle that minimizes the area of ​​the rectangle;

[0056] The position projection of each point in the new microseismic point set on the axial vector and the position projection on the lateral vector at the angle with the smallest rectangular area are obtained respectively;

[0057] The calculation formula for the fracture network length is:

[0058] L=maxL k -minL k ;

[0059] The calculation formula of the fracture network width is:

[0060] W=maxW k -min W k ;

[0061] Among them, L k is the position projection of the kth point in the new microseismic point set on the axial vector at the angle with the smallest rectangular area, W k is the position projection of the kth point in the new microseismic point set on the lateral vector at the angle with the smallest rectangular area, L is the length of the fracture network, and W is the width of the fracture network.

[0062] In the present invention, the minimum circumscribed rectangle is obtained by traversal, and the method has good stability and practicability, and the method has a small amount of calculation, and is suitable for real-time statistics of fracture network information.

[0063] The present invention will be described in more detail below by way of examples.

[0064] Embodiment 1:

[0065] In this embodiment, 200 microseismic event points are randomly generated in a cube with a length of 200 meters, a width of 100 meters, and a height of 100 meters at 60° north-east of a certain earthquake zone to perform parameter statistics of the hydraulic fracture network, wherein the azimuth range of the hydraulic fracture network is -90° to 90°, and the angle interval Δθ is 1°;

[0066] like Figure 2 As shown, this embodiment provides a method for calculating the statistics of a hydraulic fracture network based on a minimum circumscribed rectangular envelope to perform fracture network statistics on the above-mentioned earthquake event points. The specific steps are as follows:

[0067] (1) Input the microseismic point set location, perforation point location, and azimuth range:

[0068] Input microseismic point set location (x i ,y i , z i ), perforation point position (x o ,y o , z o ), azimuth range (θ 1 ,θ 2 ), the perforation point position and the microseismic point set position are combined to form a new microseismic point set position (x i ,y i , z i );

[0069] (2) obtaining the height of the hydraulic fracture network;

[0070] Through the microseismic point set position (x i ,y i , z i ) in the depth information z i , the height H of the fracture network is obtained by the following formula:

[0071] H=maxz i -minz i ;

[0072] (3) Obtain the area of ​​the circumscribed rectangle of each angle within the azimuth range;

[0073] The azimuth range (θ 1 ,θ 2 ) within any angle θ j , find the axial vector l at this angle j :

[0074] l j =(coSθ j , sinθ j );

[0075] Find the lateral vector w at this angle j :

[0076] w j =(sinθ j , cosθ j );

[0077] The location of each microseismic point (x i ,y i , z i ) in the axial vector l i The position projection on is:

[0078]

[0079] The location of each microseismic point (x i ,y i , z i ) in the lateral vector w j The position projection on is:

[0080]

[0081] Angle θ j The area of ​​the circumscribed rectangle below is:

[0082] A j =(maxL i -minL i )×(maxW i -minW i );

[0083] (4) Obtaining the angle θ at which the rectangular area is minimized;

[0084] In the azimuth range (θ 1 ,θ 2 ) Traverse to find the area of ​​the rectangle and obtain the angle θ corresponding to the minimum circumscribed rectangle area:

[0085]

[0086] (5) Obtain the length and width of the rectangle corresponding to the angle θ;

[0087] Obtain the angle θ, axial vector l and lateral vector w, and calculate the position of each microseismic point (x i ,y i , z i ) is the position projection L on the axial vector l i , the position projection W on the axial vector w i , and then the length of the fracture network is obtained:

[0088] L=maxL i -minL i ;

[0089] Calculate the width of the fracture network:

[0090] W=maxW i -minW i ;

[0091] (6) Output the length, width, height and orientation of the fracture network;

[0092] Output the length L, width W, height H, and orientation θ of the crack network.

[0093] Figure 3 is the microseismic point set and perforation point location, Figure 4 The statistical results of the hydraulic fracture network information of the method of the present invention are as follows; the test results show that the minimum circumscribed rectangular envelope and the statistical results of the hydraulic fracture network of the method of the present invention are consistent with the simulated fracture network parameters.

[0094] Embodiment 2:

[0095] The present embodiment is different from the first embodiment in that the parameters of the fracture network are statistically analyzed for the microseismic event points monitored at a certain fracture section in a horizontal well, wherein the azimuth range of the fracture network is -90° to 90°, and the angle interval Δθ is 1°.

[0096] Figure 5 is the microseismic point set and perforation point location, Figure 6 The statistical results of the hydraulic fracture network information of the method of the present invention are shown in the test results. The statistical results of the hydraulic fracture network of the present invention are consistent with the fracture morphology.

[0097] Embodiment three:

[0098] This embodiment provides a method for calculating the statistics of a hydraulic fracture network based on the minimum circumscribed rectangular envelope. Figure 1 As shown, including:

[0099] Based on the microseismic point set position, perforation point position and azimuth range, the height of the fracture network and the circumscribed rectangular area of ​​each angle within the azimuth range are obtained respectively;

[0100] Based on the area of ​​the circumscribed rectangle of each angle within the azimuth range, the angle with the smallest rectangular area is obtained;

[0101] The length and width of the fracture network are obtained based on the angle with the smallest rectangular area;

[0102] The calculation formula for the height of the fracture network is:

[0103] H=maxz i -minz i ;

[0104] Among them, zi is the height coordinate of the i-th point in the new microseismic point set formed by combining the microseismic point set position and the perforation point position, and H is the height of the fracture network;

[0105] The calculation formula for the area of ​​the circumscribed rectangle of each angle within the azimuth range is:

[0106] A j =(maxL i -minL i ) × (maxW i -minW i );

[0107] Among them, Li is the position projection of the i-th point in the new microseismic point set on the axial vector, Wi is the position projection of the i-th point in the new microseismic point set on the lateral vector, and Aj is the area of ​​the circumscribed rectangle of each angle within the azimuth range;

[0108] The calculation formula for the position projection of the i-th point in the new microseismic point set on the axial vector is:

[0109]

[0110] The calculation formula for the position projection of the i-th point in the new microseismic point set on the lateral vector is:

[0111]

[0112] Among them, x i ,y i are the horizontal and vertical coordinates of the i-th point in the new microseismic point set, l j is the axial vector at the jth azimuth angle within the azimuth range, Wj is the lateral vector at the jth azimuth angle within the azimuth range;

[0113] The expression of the axial vector at the jth azimuth angle within the azimuth range is:

[0114] l j =(cosθ j , sinθ j );

[0115] The expression of the lateral vector at the jth azimuth angle within the azimuth range is:

[0116] w j =(sinθ j , cosθ j );

[0117] Among them, θ j In the azimuth range (θ 1 ,θ2 ) Any angle;

[0118] The formula for calculating the angle with the smallest area of ​​a rectangle is:

[0119]

[0120] The angles with the smallest rectangular area based on the circumscribed rectangular area of ​​each angle within the azimuth range include: 1 ,θ 2 ) Traverse and find the area of ​​the rectangle, and obtain the angle θ corresponding to the area of ​​the minimum circumscribed rectangle;

[0121] The length and width of the fracture network are obtained based on the angle with the smallest rectangular area:

[0122] Calculate the axial vector and lateral vector at the angle that minimizes the area of ​​the rectangle;

[0123] The position projection of each point in the new microseismic point set on the axial vector and the position projection on the lateral vector at the angle with the smallest rectangular area are obtained respectively;

[0124] The calculation formula for the fracture network length is:

[0125] L=maxL k -minL k ;

[0126] The calculation formula of the fracture network width is:

[0127] W=maxW k -minW k ;

[0128] Among them, Lk is the position projection of the kth point in the new microseismic point set on the axial vector at the angle with the smallest rectangular area, Wk is the position projection of the kth point in the new microseismic point set on the lateral vector at the angle with the smallest rectangular area, L is the length of the fracture network, and W is the width of the fracture network.

[0129] Embodiment 4:

[0130] This embodiment provides a device for calculating the statistics of a hydraulic fracture network based on a minimum circumscribed rectangular envelope, including:

[0131] The area calculation module is used to obtain the circumscribed rectangular area of ​​each angle within the height and azimuth range of the fracture network based on the microseismic point set position, the perforation point position, and the azimuth range;

[0132] An angle obtaining module, used for obtaining the angle with the smallest rectangular area based on the circumscribed rectangular area of ​​each angle within the azimuth range;

[0133] The length and width obtaining module is used to obtain the length and width of the fracture network based on the angle with the smallest rectangular area;

[0134] The calculation formula for the height of the fracture network is:

[0135] H=maxz i -minz i ;

[0136] Among them, zi is the height coordinate of the i-th point in the new microseismic point set formed by combining the microseismic point set position and the perforation point position, and H is the height of the fracture network;

[0137] The calculation formula for the area of ​​the circumscribed rectangle of each angle within the azimuth range is:

[0138] A j =(maxL i -minL i )×(maxW i -minW i );

[0139] Among them, Li is the position projection of the i-th point in the new microseismic point set on the axial vector, Wi is the position projection of the i-th point in the new microseismic point set on the lateral vector, and Aj is the area of ​​the circumscribed rectangle of each angle within the azimuth range;

[0140] The calculation formula for the position projection of the i-th point in the new microseismic point set on the axial vector is:

[0141]

[0142] The calculation formula for the position projection of the i-th point in the new microseismic point set on the lateral vector is:

[0143]

[0144] Among them, x i ,y i are the horizontal and vertical coordinates of the i-th point in the new microseismic point set, l j is the axial vector at the jth azimuth angle within the azimuth range, W j is the lateral vector at the jth azimuth angle within the azimuth range;

[0145] The expression of the axial vector at the jth azimuth angle within the azimuth range is:

[0146] l j =(cosθ j , sinθ j );

[0147] The expression of the lateral vector at the jth azimuth angle within the azimuth range is:

[0148] w j =(sinθ j , cosθ j );

[0149] Among them, θ j In the azimuth range (θ 1 ,θ 2 ) Any angle;

[0150] The formula for calculating the angle with the smallest area of ​​a rectangle is:

[0151]

[0152] The angles with the smallest rectangular area based on the circumscribed rectangular area of ​​each angle within the azimuth range include: 1 ,θ 2 ) Traverse and find the area of ​​the rectangle, and obtain the angle θ corresponding to the area of ​​the minimum circumscribed rectangle;

[0153] The length and width of the fracture network are obtained based on the angle with the smallest rectangular area:

[0154] Calculate the axial vector and lateral vector at the angle that minimizes the area of ​​the rectangle;

[0155] The position projection of each point in the new microseismic point set on the axial vector and the position projection on the lateral vector at the angle with the smallest rectangular area are obtained respectively;

[0156] The calculation formula for the fracture network length is:

[0157] L=maxL k -minL k ;

[0158] The calculation formula of the fracture network width is:

[0159] W=maxW k -minW k ;

[0160] Among them, Lk is the position projection of the kth point in the new microseismic point set on the axial vector at the angle with the smallest rectangular area, Wk is the position projection of the kth point in the new microseismic point set on the lateral vector at the angle with the smallest rectangular area, L is the length of the fracture network, and W is the width of the fracture network.

[0161] Embodiment five:

[0162] An embodiment of the present invention provides an electronic device including a memory and a processor.

[0163] A memory storing executable instructions;

[0164] The processor runs the executable instructions in the memory to implement a fracture network statistical method based on a minimum circumscribed rectangular envelope.

[0165] The memory is used to store non-temporary computer-readable instructions. Specifically, the memory may include one or more computer program products, which may include various forms of computer-readable storage media, such as volatile memory and / or non-volatile memory. The volatile memory may, for example, include random access memory (RAM) and / or cache memory (cache), etc. The non-volatile memory may, for example, include read-only memory (ROM), hard disk, flash memory, etc.

[0166] The processor may be a central processing unit (CPU) or other forms of processing units with data processing capabilities and / or instruction execution capabilities, and may control other components in the electronic device to perform desired functions. In one embodiment of the present invention, the processor is used to run the computer-readable instructions stored in the memory.

[0167] Those skilled in the art should be able to understand that in order to solve the technical problem of how to obtain a good user experience, the present embodiment may also include well-known structures such as a communication bus and an interface, and these well-known structures should also be included in the protection scope of the present invention.

[0168] For detailed description of this embodiment, reference may be made to the corresponding descriptions in the aforementioned embodiments, which will not be repeated here.

[0169] Embodiment six:

[0170] An embodiment of the present invention provides a computer-readable storage medium storing a computer program. When the computer program is executed by a processor, a method for calculating the statistics of a hydraulic fracture network based on a minimum circumscribed rectangular envelope is implemented.

[0171] The computer-readable storage medium according to the embodiment of the present invention stores non-transitory computer-readable instructions, and when the non-transitory computer-readable instructions are executed by a processor, all or part of the steps of the above-mentioned methods of the embodiments of the present invention are executed.

[0172] The above-mentioned computer-readable storage media include, but are not limited to: optical storage media (e.g., CD-ROM and DVD), magneto-optical storage media (e.g., MO), magnetic storage media (e.g., magnetic tape or mobile hard disk), media with built-in rewritable non-volatile memory (e.g., memory card) and media with built-in ROM (e.g., ROM box).

[0173] The embodiment of the present invention proposes a method for statistically analyzing a hydraulic fracture network based on a minimum circumscribed rectangular envelope, which obtains the minimum circumscribed rectangular envelope of a microseismic event point set through its position information, and obtains the geometric information of the hydraulic fracture network, thereby forming a method for statistically analyzing a hydraulic fracture network based on a minimum circumscribed rectangular envelope, and realizing automatic statistics of the geometric information of the hydraulic fracture network.

[0174] The embodiments of the present invention have been described above, and the above description is exemplary, not exhaustive, and is not limited to the disclosed embodiments. Many modifications and changes will be apparent to those skilled in the art without departing from the scope and spirit of the described embodiments.

Claims

1. A statistical method for hydraulic fracture network based on the minimum circumscribed rectangular envelope, It is characterized in that include: Based on the microseismic point set position, perforation point position and azimuth range, the height of the fracture network and the circumscribed rectangular area of ​​each angle within the azimuth range are obtained respectively; Based on the area of ​​the circumscribed rectangle of each angle within the azimuth range, the angle with the smallest rectangular area is obtained; The length and width of the fracture network are obtained based on the angle with the smallest rectangular area.

2. The method according to claim 1, It is characterized in that The calculation formula of the height of the hydraulic fracture network is: H=maxz i -minz i ; Among them, zi is the height coordinate of the i-th point in the new microseismic point set formed by combining the microseismic point set position and the perforation point position, and H is the height of the fracture network.

3. The method according to claim 1, It is characterized in that The calculation formula for the area of ​​the circumscribed rectangle of each angle within the azimuth range is: THE j =(maxL i -minL i )×(maxW i -minW i ); Among them, Li is the position projection of the i-th point in the new microseismic point set on the axial vector, Wi is the position projection of the i-th point in the new microseismic point set on the lateral vector, and Aj is the area of ​​the circumscribed rectangle of each angle within the azimuth range.

4. The method according to claim 3, It is characterized in that The calculation formula for the position projection of the i-th point in the new microseismic point set on the axial vector is: The calculation formula for the position projection of the i-th point in the new microseismic point set on the lateral vector is: Among them, x i ,y i are the horizontal and vertical coordinates of the i-th point in the new microseismic point set, l j is the axial vector at the jth azimuth angle within the azimuth range, w j is the lateral vector at the jth azimuth angle within the azimuth range.

5. The method according to claim 4, It is characterized in that The expression of the axial vector at the jth azimuth angle within the azimuth range is: l j =(cosθ j ,sinθ j ); The expression of the lateral vector at the jth azimuth angle within the azimuth range is: w j =(sinθ j ,cosθ j ); Among them, θ j In the azimuth range (θ 1 ,θ 2 ) at any angle.

6. The method according to claim 5, It is characterized in that The calculation formula for the angle with the smallest rectangular area is: The angles with the smallest rectangular area based on the circumscribed rectangular area of ​​each angle within the azimuth range include: 1 ,θ 2 ) traverses to find the area of ​​the rectangle and obtains the angle θ corresponding to the area of ​​the minimum circumscribed rectangle.

7. The method according to claim 6, It is characterized in that The length and width of the fracture network are obtained based on the angle with the smallest rectangular area: Calculate the axial vector and lateral vector at the angle that minimizes the area of ​​the rectangle; The position projection of each point in the new microseismic point set on the axial vector and the position projection on the lateral vector at the angle with the smallest rectangular area are obtained respectively; The calculation formula of the fracture network length is: L=maxL k -minL k ; The calculation formula of the fracture network width is: <h2 style=";text-align:left;direction:ltr">W = maxW<h2 style=";text-align:left;direction:ltr"> k <h2 style=";text-align:left;direction:ltr"> -minW<h2 style=";text-align:left;direction:ltr"> k <h2 style=";text-align:left;direction:ltr"> ; Among them, L k is the position projection of the kth point in the new microseismic point set on the axial vector at the angle with the smallest rectangular area, W k is the position projection of the kth point in the new microseismic point set on the lateral vector at the angle with the smallest rectangular area, L is the length of the fracture network, and W is the width of the fracture network.

8. A statistical device for hydraulic fracture network based on the minimum circumscribed rectangular envelope, It is characterized in that include: The area calculation module is used to obtain the circumscribed rectangular area of ​​each angle within the height and azimuth range of the fracture network based on the microseismic point set position, the perforation point position, and the azimuth range; An angle obtaining module, used for obtaining the angle with the smallest rectangular area based on the circumscribed rectangular area of ​​each angle within the azimuth range; The length and width obtaining module is used to obtain the length and width of the fracture network based on the angle with the smallest rectangular area.

9. An electronic device, It is characterized in that The electronic device comprises: A memory storing executable instructions; A processor, wherein the processor runs the executable instructions in the memory to implement the minimum circumscribed rectangular envelope-based fracture network statistics method according to any one of claims 1-7.

10. A computer-readable storage medium, It is characterized in that The computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the method for statistically analyzing a fracture network based on a minimum circumscribed rectangular envelope according to any one of claims 1 to 7 is implemented.

Images