A crack morphology analysis method, device, equipment and medium
By classifying microseismic data, calculating convex hulls, and analyzing normal vectors, and combining this with random fractal methods, the problem of inaccurate determination of crack type and orientation in existing technologies has been solved, achieving efficient crack morphology analysis.
Patent Information
- Authority / Receiving Office
- CN · China
- Patent Type
- Patents(China)
- Current Assignee / Owner
- YANCHANG OIL MINE ADMINISTRATION BUREAU NANNIWAN EXTRACTION CO
- Filing Date
- 2023-08-23
- Publication Date
- 2026-06-12
AI Technical Summary
Existing microseismic monitoring technologies cannot accurately determine the type and orientation of cracks, have low computational efficiency, and cannot comprehensively and accurately depict the true morphology of cracks.
After acquiring microseismic data and classifying the data, we perform convex hull and projected area calculations to determine the main fracture surface. We then use normal vector calculations to determine the fracture surface dip angle and orientation, and combine this with random fractal methods to analyze the fracture morphology.
Accurately determine the type and direction of cracks, reduce the amount of calculation, improve calculation efficiency, and achieve a comprehensive and accurate depiction of the true morphology of cracks.
Smart Images
Figure CN117075196B_ABST
Abstract
Description
Technical Field
[0001] This invention relates to the field of petroleum exploration and development, and in particular to a method, apparatus, equipment and medium for fracture morphology analysis. Background Technology
[0002] The basic principle of microseismic monitoring technology is to detect microseismic events caused by changes in rock internal stress. This involves deploying highly sensitive seismic sensors around the fractured well to continuously record microseismic activity generated during the fracturing process due to changes in the physical properties of the oil and gas reservoir. Current microseismic monitoring technologies can only monitor fractures during fracturing and cannot effectively monitor and identify microseismic events in other geological processes. Furthermore, they can only provide approximate fracture location information and cannot accurately characterize and analyze specific fracture parameters (such as length, width, and direction). Currently, there are two main methods for monitoring and identifying microseismic events: one involves extracting fracture surfaces from microseismic points, primarily through cluster analysis and probability density function estimation. These methods identify sets with similar characteristics based on the spatial distribution characteristics of microseismic point data, thus characterizing the fracture surface. The other method is a random fractal fracture network characterization method, which views the fracture morphology as a self-similar system composed of primary and secondary fractures. Based on the similarity between fracture propagation and plant growth, and utilizing the L-system fractal principle, fractal fracture morphology simulation is conducted according to fractal tree generation rules. Different fracture node coordinates are obtained according to different fracture generation rules, resulting in different fracture morphologies. The corresponding drawbacks are as follows: One method, because it fits a discrete crack network, ignores the differences between large-scale main cracks and small-scale secondary crack networks at the ends; the fitting process requires significant computational resources, resulting in low computational efficiency for large-scale microseismic point data. Another method, the random fractal crack network characterization method, defaults to projecting microseismic points onto a horizontal plane before iterative calculation, interpreting the crack distribution morphology only from a horizontal perspective and ignoring the spatial distribution characteristics of microseismic points: if the crack type is near-horizontal, the interpretation result is reasonable, and the spatial distribution of cracks does not need to be considered; if the crack type is near-vertical, the true crack morphology cannot be obtained from the horizontal projection, and a vertical projection is necessary; if the crack type is between the horizontal and vertical planes, neither horizontal nor vertical projection can obtain the true crack morphology.
[0003] As can be seen from the above, how to accurately determine the type and orientation of cracks, reduce the amount of computation of microseismic point data, improve computational efficiency, and achieve a comprehensive and accurate characterization and analysis of the true morphology of cracks are problems that need to be solved in this field. Summary of the Invention
[0004] In view of this, the purpose of this invention is to provide a method, apparatus, device, and medium for crack morphology analysis, which can accurately determine the crack type and orientation, reduce the computational workload of microseismic point data, improve computational efficiency, and achieve a comprehensive and accurate characterization and analysis of the true morphology of cracks. The specific solution is as follows:
[0005] In a first aspect, this application discloses a crack morphology analysis method, including:
[0006] Acquire microseismic data of the hydraulic fracturing fractures, and classify the microseismic data to obtain target microseismic data with fracture morphology as the main fracture.
[0007] The target microseismic data is subjected to convex hull calculation to obtain a three-dimensional convex hull. The projected area of the three-dimensional convex hull is calculated to obtain each two-dimensional convex hull and its corresponding area. The area of the two-dimensional convex hull with the largest value is taken as the target two-dimensional convex hull area.
[0008] The main crack surface is determined based on the target two-dimensional convex hull area corresponding to the target two-dimensional convex hull area. The normal vector of the main crack surface is calculated to obtain the crack surface inclination angle and crack surface orientation. The crack type is determined based on the crack surface inclination angle and crack surface orientation.
[0009] Crack morphology analysis is performed based on the microseismic data, the main crack surface, the crack type, and the crack surface orientation.
[0010] Optionally, the step of classifying the microseismic data to obtain target microseismic data with fracture morphology as the main fracture includes:
[0011] The microseismic data is classified using a preset data classification algorithm to obtain target microseismic data with fracture morphology as the main fracture; wherein, the data classification algorithm includes a clustering analysis algorithm.
[0012] Optionally, the step of performing convex hull calculation on the target microseismic data to obtain a three-dimensional convex hull, and calculating the projected area of the three-dimensional convex hull to obtain each two-dimensional convex hull and its corresponding area, includes:
[0013] The convex hull algorithm is used to calculate the convex hull of the target microseismic data to obtain a three-dimensional convex hull;
[0014] The three-dimensional convex hull is projected to obtain the projected two-dimensional convex hulls. The area of each two-dimensional convex hull is calculated to obtain the corresponding two-dimensional convex hull area.
[0015] Optionally, determining the main crack surface based on the target two-dimensional convex hull corresponding to the area of the target two-dimensional convex hull includes:
[0016] Determine the target two-dimensional convex hull corresponding to the area of the target two-dimensional convex hull;
[0017] The plane at the center of the microseismic data points in the target two-dimensional convex hull is taken as the main crack surface.
[0018] Optionally, the step of calculating the normal vector of the main crack surface to obtain the crack surface dip angle and crack surface orientation includes:
[0019] Determine the normal direction of the main crack surface;
[0020] The normal vector of the main crack surface is calculated based on the normal direction to obtain the crack surface dip angle and crack surface orientation.
[0021] Optionally, the fracture morphology analysis based on the microseismic data, the main fracture surface, the fracture type, and the fracture surface orientation includes:
[0022] All the microseismic data are projected onto the main fracture surface to obtain the projected main fracture surface;
[0023] Crack morphology analysis is achieved based on the projected main crack surface, the crack type, and the crack surface orientation.
[0024] Optionally, the crack morphology analysis based on the projected main crack surface, the crack type, and the crack surface orientation includes:
[0025] A random fractal crack network characterization method is adopted, and crack morphology analysis is achieved based on the projected main crack surface, the crack type, and the crack surface orientation.
[0026] Secondly, this application discloses a crack morphology analysis device, comprising:
[0027] The data classification module is used to acquire microseismic data of hydraulic fracturing fractures and classify the microseismic data to obtain target microseismic data with fracture morphology as the main fracture.
[0028] The projected area calculation module is used to perform convex hull calculation on the target microseismic data to obtain a three-dimensional convex hull, and to perform projected area calculation on the three-dimensional convex hull to obtain each two-dimensional convex hull and its corresponding area. The area of the two-dimensional convex hull with the largest value is taken as the target two-dimensional convex hull area.
[0029] The normal vector calculation module is used to determine the main crack surface based on the target two-dimensional convex hull corresponding to the area of the target two-dimensional convex hull, perform normal vector calculation on the main crack surface to obtain the crack surface inclination angle and crack surface orientation, and determine the crack type based on the crack surface inclination angle and crack surface orientation.
[0030] The crack morphology analysis module is used to perform crack morphology analysis based on the microseismic data, the main crack surface, the crack type, and the crack surface orientation.
[0031] Thirdly, this application discloses an electronic device, including:
[0032] Memory, used to store computer programs;
[0033] A processor is used to execute the computer program to implement the aforementioned crack morphology analysis method.
[0034] Fourthly, this application discloses a computer storage medium for storing a computer program; wherein, when the computer program is executed by a processor, it implements the steps of the aforementioned disclosed crack morphology analysis method.
[0035] As can be seen, this application provides a method for fracture morphology analysis, including acquiring microseismic data of hydraulic fracturing fractures, classifying the microseismic data to obtain target microseismic data with fracture morphology as the main fracture; performing convex hull calculation on the target microseismic data to obtain a three-dimensional convex hull, calculating the projected area of the three-dimensional convex hull to obtain each two-dimensional convex hull and its corresponding area, and taking the area of the two-dimensional convex hull with the largest value as the target two-dimensional convex hull area; determining the main fracture surface based on the target two-dimensional convex hull corresponding to the target two-dimensional convex hull area, calculating the normal vector of the main fracture surface to obtain the fracture surface dip angle and fracture surface orientation, and determining the fracture type based on the fracture surface dip angle and fracture surface orientation; and realizing fracture morphology analysis based on the microseismic data, the main fracture surface, the fracture type, and the fracture surface orientation. This application performs convex hull calculation and projected area calculation on the target microseismic data of the main fracture to determine the main fracture surface. The normal vector of the main fracture surface is then calculated to obtain the fracture surface dip angle and fracture surface orientation, thereby accurately determining the fracture type and orientation. This reduces the amount of calculation and improves the calculation efficiency, ultimately achieving a comprehensive and accurate characterization and analysis of the true morphology of the fracture. Attached Figure Description
[0036] To more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are only embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on the provided drawings without creative effort.
[0037] Figure 1 This is a flowchart of a crack morphology analysis method disclosed in this application;
[0038] Figure 2This is a flowchart illustrating a crack morphology analysis method disclosed in this application.
[0039] Figure 3 This is a flowchart of another crack morphology analysis method disclosed in this application;
[0040] Figure 4 This application discloses a cloud map showing the distribution of well microseismic points in target microseismic data.
[0041] Figure 5 This is a composite image of the main fracture in a fracturing segment disclosed in this application;
[0042] Figure 6(a) is a diagram of the main fracture surface and fracture network of the first section of the H well disclosed in this application.
[0043] Figure 6(b) is a diagram of the main fracture surface and fracture network of the second fracturing section of Well H disclosed in this application;
[0044] Figure 6(c) is a diagram of the main fracture surface and fracture network of the third section of the H well disclosed in this application.
[0045] Figure 6(d) is a diagram depicting the main fracture surface and fracture network of the fracturing section of the fourth segment of Well H disclosed in this application.
[0046] Figure 6(e) is a diagram of the main fracture surface and fracture network of the fracturing section of the fifth stage of Well H disclosed in this application.
[0047] Figure 6(f) is a diagram of the main fracture surface and fracture network of the fracturing section of the sixth segment of Well H disclosed in this application;
[0048] Figure 6(g) is a diagram of the main fracture surface and fracture network of the fracturing section of the 7th segment of Well H disclosed in this application;
[0049] Figure 6(h) is a diagram depicting the main fracture surface and fracture network of the fracturing section of the 8th segment of Well H disclosed in this application.
[0050] Figure 6(I) is a diagram of the main fracture surface and fracture network of the fracturing section of the 9th segment of Well H disclosed in this application;
[0051] Figure 7 This is a schematic diagram of a crack morphology analysis device disclosed in this application;
[0052] Figure 8 This application provides a structural diagram of an electronic device. Detailed Implementation
[0053] The technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.
[0054] The basic principle of microseismic monitoring technology is to detect microseismic events caused by changes in rock internal stress. This involves deploying highly sensitive seismic sensors around the fractured well to continuously record microseismic activity generated during the fracturing process due to changes in the physical properties of the oil and gas reservoir. Current microseismic monitoring technologies can only monitor fractures during fracturing and cannot effectively monitor and identify microseismic events in other geological processes. Furthermore, they can only provide approximate fracture location information and cannot accurately characterize and analyze specific fracture parameters (such as length, width, and direction). Currently, there are two main methods for monitoring and identifying microseismic events: one involves extracting fracture surfaces from microseismic points, primarily through cluster analysis and probability density function estimation. These methods identify sets with similar characteristics based on the spatial distribution characteristics of microseismic point data, thus characterizing the fracture surface. The other method is a random fractal fracture network characterization method, which views the fracture morphology as a self-similar system composed of primary and secondary fractures. Based on the similarity between fracture propagation and plant growth, and utilizing the L-system fractal principle, fractal fracture morphology simulation is conducted according to fractal tree generation rules. Different fracture node coordinates are obtained according to different fracture generation rules, resulting in different fracture morphologies. The corresponding drawbacks are as follows: Firstly, because the fitted crack network is a discrete crack network, it ignores the differences between large-scale main cracks and small-scale secondary crack networks at the ends; the fitting process requires significant computational resources, resulting in low computational efficiency for large-scale microseismic point data. Secondly, the random fractal crack network characterization method defaults to projecting microseismic points onto a horizontal plane before iterative calculation, interpreting the crack distribution morphology only from a horizontal perspective and ignoring the spatial distribution characteristics of microseismic points: if the crack type is near-horizontal, the interpretation result of this method is reasonable, and the spatial distribution of cracks does not need to be considered; if the crack type is near-vertical, the true crack morphology cannot be obtained from the horizontal plane projection, and should be obtained from the vertical direction; if the crack type is between the horizontal and vertical planes, the true crack morphology cannot be obtained from either horizontal or vertical projection. Therefore, accurately determining the crack type and orientation, reducing the computational load of microseismic point data, improving computational efficiency, and achieving a comprehensive and accurate characterization and analysis of the true crack morphology are problems that need to be solved in this field.
[0055] See Figure 1As shown in the figure, an embodiment of the present invention discloses a crack morphology analysis method, which may specifically include:
[0056] Step S11: Acquire microseismic data of the hydraulic fracturing fracture, and classify the microseismic data to obtain target microseismic data with fracture morphology as the main fracture.
[0057] In this embodiment, microseismic data of hydraulic fracturing fractures are acquired, and the microseismic data are classified using a preset data classification algorithm to obtain target microseismic data with fracture morphology as the main fracture; wherein, the data classification algorithm includes a clustering analysis algorithm.
[0058] Specifically, cluster analysis algorithms are used to classify the microseismic data of the main fracture, and then noise reduction processing is performed on the microseismic data of the main fracture to obtain the target microseismic data.
[0059] Step S12: Perform convex hull calculation on the target microseismic data to obtain a three-dimensional convex hull. Calculate the projected area of the three-dimensional convex hull to obtain each two-dimensional convex hull and its corresponding area. Take the area of the two-dimensional convex hull with the largest value as the target two-dimensional convex hull area.
[0060] Step S13: Determine the main crack surface based on the target two-dimensional convex hull corresponding to the area of the target two-dimensional convex hull, calculate the normal vector of the main crack surface to obtain the crack surface inclination angle and crack surface orientation, and determine the crack type based on the crack surface inclination angle and crack surface orientation.
[0061] Step S14: Based on the microseismic data, the main fracture surface, the fracture type, and the fracture surface orientation, perform fracture morphology analysis.
[0062] In this embodiment, all the microseismic data are projected onto the main fracture surface to obtain the projected main fracture surface; fracture morphology analysis is then performed based on the projected main fracture surface, the fracture type, and the fracture surface orientation. Specifically, a random fractal fracture network characterization method is used, and fracture morphology analysis is performed based on the projected main fracture surface, the fracture type, and the fracture surface orientation. In other words, the specific process for fracture morphology analysis is as follows: all microseismic data are projected onto the main fracture surface of the point cloud data, and then the random fractal fracture network characterization method is used to characterize the fracture morphology of the projected microseismic data.
[0063] The specific crack morphology analysis process is as follows: Figure 2As shown, firstly, microseismic data of the hydraulic fracturing fracture is acquired; then, the process of finding the plane of the point center is executed, namely: using a cluster analysis algorithm to determine the microseismic data of the main fracture, then performing noise reduction to obtain the target microseismic data, then performing convex hull calculation and projected area calculation to determine the main fracture surface, and projecting all the microseismic data onto the main fracture surface to obtain the projected main fracture surface; finally, based on the projected main fracture surface, the fracture type, and the fracture surface orientation, fracture morphology analysis is performed.
[0064] In this embodiment, microseismic data of hydraulic fracturing is acquired, and the microseismic data is classified to obtain target microseismic data with the main fracture morphology. Convex hull calculation is performed on the target microseismic data to obtain a three-dimensional convex hull. The projected area of the three-dimensional convex hull is calculated to obtain each two-dimensional convex hull and its corresponding area. The area of the two-dimensional convex hull with the largest value is taken as the target two-dimensional convex hull area. The main fracture surface is determined based on the target two-dimensional convex hull corresponding to the target two-dimensional convex hull area. The normal vector of the main fracture surface is calculated to obtain the fracture surface dip angle and fracture surface orientation. The fracture type is determined based on the fracture surface dip angle and fracture surface orientation. Fracture morphology analysis is achieved based on the microseismic data, the main fracture surface, the fracture type, and the fracture surface orientation. This application performs convex hull calculation and projected area calculation on the target microseismic data of the main fracture to determine the main fracture surface. The normal vector calculation of the main fracture surface yields the fracture surface dip angle and fracture surface orientation, thereby accurately determining the fracture type and orientation. This reduces computational load, improves computational efficiency, and ultimately achieves a comprehensive and accurate characterization and analysis of the true fracture morphology.
[0065] See Figure 3 As shown in the figure, an embodiment of the present invention discloses a crack morphology analysis method, which may specifically include:
[0066] Step S21: Acquire microseismic data of the hydraulic fracturing fracture, and classify the microseismic data to obtain target microseismic data with fracture morphology as the main fracture.
[0067] Step S22: Use the convex hull algorithm to calculate the convex hull of the target microseismic data to obtain a three-dimensional convex hull. Project the three-dimensional convex hull to obtain the projected two-dimensional convex hulls. Calculate the area of each two-dimensional convex hull to obtain the corresponding two-dimensional convex hull area. Take the area of the two-dimensional convex hull with the largest value as the target two-dimensional convex hull area.
[0068] Step S23: Determine the target two-dimensional convex hull corresponding to the area of the target two-dimensional convex hull, take the plane at the center of the microseismic data points in the target two-dimensional convex hull as the main fracture surface, determine the normal direction of the main fracture surface, calculate the normal vector of the main fracture surface based on the normal direction to obtain the fracture surface dip angle and fracture surface orientation, and determine the fracture type according to the fracture surface dip angle and fracture surface orientation.
[0069] Step S24: Based on the microseismic data, the main fracture surface, the fracture type, and the fracture surface orientation, perform fracture morphology analysis.
[0070] Taking the fracture morphology analysis of a horizontal well as an example, the fracturing section of horizontal well H is divided into 9 segments. Clustering is used to obtain the microseismic data points of the main fracture in each segment. After noise reduction processing, the microseismic data points of the main fracture are as follows: Figure 4 As shown, where, Figure 4 The left image is a top view of the distribution cloud map of microseismic points in the target microseismic data wells, and the right image is a side view of the distribution cloud map of microseismic points in the target microseismic data wells. A convex hull algorithm is used to calculate the 3D convex hull. The projected area of the 3D convex hull is then calculated to obtain the 2D convex hull and its corresponding area, thereby determining the main fracture surface of the microseismic data points in each fracturing segment. The composite image of the main fracture surface of the fracturing segments 1-9 in Well H is shown below. Figure 5 As shown in the figure. Based on the normal direction of the main fracture surface of each fracturing segment, the dip angle of the fracture surface is calculated. All microseismic data are projected onto the main fracture surface of the point cloud data. Then, the fracture morphology of the projected microseismic data is characterized by the random fractal fracture network characterization method, thereby realizing fracture morphology analysis. Specifically, the main fracture surface and fracture network characterization of the fracturing segments 1, 2, 3, 4, 5, 6, 7, 8, and 9 of well H are respectively shown in Figure 6(a), Figure 6(b), Figure 6(c), Figure 6(d), Figure 6(e), Figure 6(f), Figure 6(g), Figure 6(h), and Figure 6(I).
[0071] This application can efficiently and accurately locate the main fracture surface from microseismic data and pre-determine the fracture type, thereby enabling better characterization and analysis of fractures. For complex subsurface structures, it can better reflect the true state of the reservoir. In addition, the use of a fast convex hull algorithm greatly reduces the amount of computation and improves computational efficiency.
[0072] In this embodiment, microseismic data of hydraulic fracturing is acquired, and the microseismic data is classified to obtain target microseismic data with the main fracture morphology. Convex hull calculation is performed on the target microseismic data to obtain a three-dimensional convex hull. The projected area of the three-dimensional convex hull is calculated to obtain each two-dimensional convex hull and its corresponding area. The area of the two-dimensional convex hull with the largest value is taken as the target two-dimensional convex hull area. The main fracture surface is determined based on the target two-dimensional convex hull corresponding to the target two-dimensional convex hull area. The normal vector of the main fracture surface is calculated to obtain the fracture surface dip angle and fracture surface orientation. The fracture type is determined based on the fracture surface dip angle and fracture surface orientation. Fracture morphology analysis is achieved based on the microseismic data, the main fracture surface, the fracture type, and the fracture surface orientation. This application performs convex hull calculation and projected area calculation on the target microseismic data of the main fracture to determine the main fracture surface. The normal vector calculation of the main fracture surface yields the fracture surface dip angle and fracture surface orientation, thereby accurately determining the fracture type and orientation. This reduces computational load, improves computational efficiency, and ultimately achieves a comprehensive and accurate characterization and analysis of the true fracture morphology.
[0073] See Figure 7 As shown, an embodiment of the present invention discloses a crack morphology analysis device, which may specifically include:
[0074] The data classification module 11 is used to acquire microseismic data of hydraulic fracturing fractures and classify the microseismic data to obtain target microseismic data with fracture morphology as the main fracture.
[0075] The projection area calculation module 12 is used to perform convex hull calculation on the target microseismic data to obtain a three-dimensional convex hull, and to perform projection area calculation on the three-dimensional convex hull to obtain each two-dimensional convex hull and the corresponding two-dimensional convex hull area, and to take the area of the two-dimensional convex hull with the largest value as the target two-dimensional convex hull area.
[0076] The normal vector calculation module 13 is used to determine the main crack surface based on the target two-dimensional convex hull corresponding to the area of the target two-dimensional convex hull, calculate the normal vector of the main crack surface to obtain the crack surface inclination angle and crack surface orientation, and determine the crack type based on the crack surface inclination angle and crack surface orientation.
[0077] The crack morphology analysis module 14 is used to perform crack morphology analysis based on the microseismic data, the main crack surface, the crack type, and the crack surface orientation.
[0078] In this embodiment, microseismic data of hydraulic fracturing is acquired, and the microseismic data is classified to obtain target microseismic data with the main fracture morphology. Convex hull calculation is performed on the target microseismic data to obtain a three-dimensional convex hull. The projected area of the three-dimensional convex hull is calculated to obtain each two-dimensional convex hull and its corresponding area. The area of the two-dimensional convex hull with the largest value is taken as the target two-dimensional convex hull area. The main fracture surface is determined based on the target two-dimensional convex hull corresponding to the target two-dimensional convex hull area. The normal vector of the main fracture surface is calculated to obtain the fracture surface dip angle and fracture surface orientation. The fracture type is determined based on the fracture surface dip angle and fracture surface orientation. Fracture morphology analysis is achieved based on the microseismic data, the main fracture surface, the fracture type, and the fracture surface orientation. This application performs convex hull calculation and projected area calculation on the target microseismic data of the main fracture to determine the main fracture surface. The normal vector calculation of the main fracture surface yields the fracture surface dip angle and fracture surface orientation, thereby accurately determining the fracture type and orientation. This reduces computational load, improves computational efficiency, and ultimately achieves a comprehensive and accurate characterization and analysis of the true fracture morphology.
[0079] In some specific embodiments, the data classification module 11 may specifically include:
[0080] The data classification module is used to classify the microseismic data using a preset data classification algorithm to obtain target microseismic data with fracture morphology as the main fracture; wherein, the data classification algorithm includes a clustering analysis algorithm.
[0081] In some specific embodiments, the projected area calculation module 12 may specifically include:
[0082] The convex hull calculation module is used to perform convex hull calculation on the target microseismic data using the convex hull algorithm to obtain a three-dimensional convex hull.
[0083] The two-dimensional convex hull area calculation module is used to project the three-dimensional convex hull to obtain the projected two-dimensional convex hulls, and to calculate the area of each two-dimensional convex hull to obtain the corresponding two-dimensional convex hull area.
[0084] In some specific embodiments, the projected area calculation module 12 may specifically include:
[0085] The target two-dimensional convex hull determination module is used to determine the target two-dimensional convex hull corresponding to the area of the target two-dimensional convex hull;
[0086] The main crack surface determination module is used to take the plane at the center of the microseismic data points in the target two-dimensional convex hull as the main crack surface.
[0087] In some specific embodiments, the normal vector calculation module 13 may specifically include:
[0088] The normal direction determination module is used to determine the normal direction of the main crack surface;
[0089] The normal vector calculation module is used to calculate the normal vector of the main crack surface based on the normal direction to obtain the crack surface dip angle and crack surface orientation.
[0090] In some specific embodiments, the crack morphology analysis module 14 may specifically include:
[0091] The microseismic data projection module is used to project all the microseismic data onto the main fracture surface to obtain the projected main fracture surface.
[0092] The crack morphology analysis module is used to perform crack morphology analysis based on the projected main crack surface, the crack type, and the crack surface orientation.
[0093] In some specific embodiments, the crack morphology analysis module 14 may specifically include:
[0094] The crack morphology analysis module is used to perform crack morphology analysis based on the projected main crack surface, the crack type, and the crack surface orientation using a random fractal crack network characterization method.
[0095] Figure 8 This is a schematic diagram of an electronic device provided in an embodiment of this application. The electronic device 20 may specifically include: at least one processor 21, at least one memory 22, a power supply 23, a communication interface 24, an input / output interface 25, and a communication bus 26. The memory 22 stores a computer program, which is loaded and executed by the processor 21 to implement the relevant steps in the crack morphology analysis method performed by the electronic device disclosed in any of the foregoing embodiments.
[0096] In this embodiment, the power supply 23 is used to provide operating voltage for each hardware device on the electronic device 20; the communication interface 24 can create a data transmission channel between the electronic device 20 and external devices, and the communication protocol it follows can be any communication protocol applicable to the technical solution of this application, and is not specifically limited here; the input / output interface 25 is used to acquire external input data or output data to the outside world, and its specific interface type can be selected according to specific application needs, and is not specifically limited here.
[0097] In addition, the memory 22, as a carrier for resource storage, can be a read-only memory, random access memory, disk or optical disk, etc. The resources stored on it include operating system 221, computer program 222 and data 223, etc., and the storage method can be temporary storage or permanent storage.
[0098] The operating system 221 manages and controls the various hardware devices on the electronic device 20 and the computer program 222 to enable the processor 21 to perform calculations and processing on the data 223 in the memory 22. The operating system 221 can be Windows, Unix, Linux, etc. The computer program 222, in addition to including a computer program capable of performing the crack morphology analysis method executed by the electronic device 20 as disclosed in any of the foregoing embodiments, may further include computer programs capable of performing other specific tasks. The data 223 may include data received by the crack morphology analysis device from external devices, as well as data collected by its own input / output interface 25.
[0099] The steps of the methods or algorithms described in conjunction with the embodiments disclosed herein can be implemented directly by hardware, a software module executed by a processor, or a combination of both. The software module can be located in random access memory (RAM), main memory, read-only memory (ROM), electrically programmable ROM, electrically erasable programmable ROM, registers, hard disk, removable disk, CD-ROM, or any other form of storage medium known in the art.
[0100] Furthermore, this application also discloses a computer-readable storage medium storing a computer program. When the computer program is loaded and executed by a processor, it implements the crack morphology analysis method steps disclosed in any of the foregoing embodiments.
[0101] Finally, it should be noted that in this document, relational terms such as "first" and "second" are used only to distinguish one entity or operation from another, and do not necessarily require or imply any such actual relationship or order between these entities or operations. Furthermore, the terms "comprising," "including," or any other variations thereof are intended to cover non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements includes not only those elements but also other elements not expressly listed, or elements inherent to such a process, method, article, or apparatus. Without further limitations, an element defined by the phrase "comprising one..." does not exclude the presence of other identical elements in the process, method, article, or apparatus that includes said element.
[0102] The present invention provides a detailed description of a crack morphology analysis method, apparatus, device, and storage medium. Specific examples have been used to illustrate the principles and implementation methods of the present invention. The descriptions of the above embodiments are only for the purpose of helping to understand the method and core ideas of the present invention. At the same time, for those skilled in the art, there will be changes in the specific implementation methods and application scope based on the ideas of the present invention. Therefore, the content of this specification should not be construed as a limitation of the present invention.
Claims
1. A method for analyzing crack morphology, characterized in that, include: Acquire microseismic data of the hydraulic fracturing fractures, and classify the microseismic data to obtain target microseismic data with fracture morphology as the main fracture. The target microseismic data is subjected to convex hull calculation to obtain a three-dimensional convex hull. The projected area of the three-dimensional convex hull is calculated to obtain each two-dimensional convex hull and its corresponding area. The area of the two-dimensional convex hull with the largest value is taken as the target two-dimensional convex hull area. The main crack surface is determined based on the target two-dimensional convex hull area corresponding to the target two-dimensional convex hull area. The normal vector of the main crack surface is calculated to obtain the crack surface inclination angle and crack surface orientation. The crack type is determined based on the crack surface inclination angle and crack surface orientation. Crack morphology analysis is performed based on the microseismic data, the main crack surface, the crack type, and the crack surface orientation.
2. The crack morphology analysis method according to claim 1, characterized in that, The process of classifying the microseismic data to obtain target microseismic data with fracture morphology as the main fracture includes: The microseismic data is classified using a preset data classification algorithm to obtain target microseismic data with fracture morphology as the main fracture; wherein, the data classification algorithm includes a clustering analysis algorithm.
3. The crack morphology analysis method according to claim 1, characterized in that, The step of performing convex hull calculation on the target microseismic data to obtain a three-dimensional convex hull, and then calculating the projected area of the three-dimensional convex hull to obtain each two-dimensional convex hull and its corresponding area, includes: The convex hull algorithm is used to calculate the convex hull of the target microseismic data to obtain a three-dimensional convex hull; The three-dimensional convex hull is projected to obtain the projected two-dimensional convex hulls. The area of each two-dimensional convex hull is calculated to obtain the corresponding two-dimensional convex hull area.
4. The crack morphology analysis method according to claim 1, characterized in that, The step of determining the main crack surface based on the target two-dimensional convex hull corresponding to the area of the target two-dimensional convex hull includes: Determine the target two-dimensional convex hull corresponding to the area of the target two-dimensional convex hull; The plane at the center of the microseismic data points in the target two-dimensional convex hull is taken as the main crack surface.
5. The crack morphology analysis method according to claim 1, characterized in that, The step of calculating the normal vector of the main crack surface to obtain the crack surface dip angle and crack surface orientation includes: Determine the normal direction of the main crack surface; The normal vector of the main crack surface is calculated based on the normal direction to obtain the crack surface dip angle and crack surface orientation.
6. The crack morphology analysis method according to any one of claims 1 to 5, characterized in that, The fracture morphology analysis based on the microseismic data, the main fracture surface, the fracture type, and the fracture surface orientation includes: All the microseismic data are projected onto the main fracture surface to obtain the projected main fracture surface; Crack morphology analysis is achieved based on the projected main crack surface, the crack type, and the crack surface orientation.
7. The crack morphology analysis method according to claim 6, characterized in that, The crack morphology analysis based on the projected main crack surface, crack type, and crack surface orientation includes: A random fractal crack network characterization method is adopted, and crack morphology analysis is achieved based on the projected main crack surface, the crack type, and the crack surface orientation.
8. A crack morphology analysis device, characterized in that, include: The data classification module is used to acquire microseismic data of hydraulic fracturing fractures and classify the microseismic data to obtain target microseismic data with fracture morphology as the main fracture. The projected area calculation module is used to perform convex hull calculation on the target microseismic data to obtain a three-dimensional convex hull, and to perform projected area calculation on the three-dimensional convex hull to obtain each two-dimensional convex hull and its corresponding area. The area of the two-dimensional convex hull with the largest value is taken as the target two-dimensional convex hull area. The normal vector calculation module is used to determine the main crack surface based on the target two-dimensional convex hull corresponding to the area of the target two-dimensional convex hull, calculate the normal vector of the main crack surface to obtain the crack surface inclination angle and crack surface orientation, and determine the crack type based on the crack surface inclination angle and crack surface orientation. The crack morphology analysis module is used to perform crack morphology analysis based on the microseismic data, the main crack surface, the crack type, and the crack surface orientation.
9. An electronic device, characterized in that, include: Memory, used to store computer programs; A processor for executing the computer program to implement the crack morphology analysis method as described in any one of claims 1 to 7.
10. A computer-readable storage medium, characterized in that, Used to store computer programs; wherein, when the computer programs are executed by a processor, they implement the crack morphology analysis method as described in any one of claims 1 to 7.