B value-based microseismic natural fracture separation fitting method and device

Abstract

The invention discloses a microseismic natural fracture separation fitting method and device based on a B value, which can be used in the field of petroleum and natural gas exploration and development, and the method comprises the steps: obtaining microseismic data of a target mining area, and extracting magnitude distribution features and B value features according to the microseismic data; according to the earthquake magnitude distribution characteristics and the B value characteristics, according to induction reasons of the microseismic events of the target mining area, separating the microseismic events, and grading the microseismic events; judging whether each level of microseismic event is confused or not, extracting a main fracture surface of the natural fracture according to a grading result of each level of microseismic event and a judgment result of whether confusion occurs or not, and fitting the fracture surface of each microseismic event to obtain a well periphery fracture model; and according to the fracture models around the wells, fracturing design is conducted on the horizontal wells corresponding to the fracture models around the wells. The method can effectively predict the spatial distribution characteristics of the fault and fracture-developed zone, and improves the exploration and development benefits of shale gas.

Description

TECHNICAL FIELD

[0001] The present application relates to the field of oil and gas exploration and development, and particularly relates to a microseismic natural fracture separation fitting method and device based on B value. BACKGROUND

[0002] This section is intended to provide background or context to the embodiments of the application. The description herein does not constitute admission that the information provided herein is prior art.

[0003] As a clean energy, the development and utilization of shale gas can help reduce dependence on external energy and improve energy security. Deep shale gas resources are abundant, and are the main battlefield for achieving sustained production of shale gas. The shale gas resource in a certain mining right area is 1510 billion cubic meters, and the proportion of deep shale gas is 84%. However, the geological engineering conditions of deep shale gas are extremely complex, are affected by multiple tectonic movements, and have complex deep structures and developed natural fractures, which are much more difficult to develop than expected. Especially affected by fracturing sleeve change and fracturing channeling, it is difficult to accurately predict the spatial distribution characteristics of fault and fracture development zones, which restricts the scale benefit development process, and there is no effective solution to this problem in the prior art. SUMMARY

[0004] The embodiments of the present application provide a microseismic natural fracture separation fitting method based on B value, which is used to effectively predict the spatial distribution characteristics of fault and fracture development zones and improve the benefit of shale gas exploration and development. The method comprises:

[0005] Obtaining microseismic data of a target mining area, extracting magnitude distribution characteristics and B value characteristics according to the microseismic data;

[0006] According to the magnitude distribution characteristics and the B value characteristics, separating each microseismic event according to the inducing reason of the microseismic event of the target mining area, and classifying each microseismic event according to the separation result;

[0007] Judging whether each level of microseismic event is confused or not, extracting a main fracture surface of a natural fracture according to the classification result of each level of microseismic event and the judgment result of whether the confusion occurs or not, and fitting the fracture surface of each microseismic event respectively;

[0008] According to the fitting result, depicting the wellbore fracture of each microseismic event to obtain a wellbore fracture model;

[0009] According to each wellbore fracture model, performing fracturing design on a horizontal well corresponding to each wellbore fracture model.

[0010] The embodiments of the present application also provide a microseismic natural fracture separation fitting device based on B value, which is used to effectively predict the spatial distribution characteristics of fault and fracture development zones and improve the benefit of shale gas exploration and development. The device comprises:

[0011] The data acquisition module is configured to acquire microseismic data of a target mining area, and extract magnitude distribution characteristics and B value characteristics from the microseismic data.

[0012] The separation module is configured to separate each microseismic event according to the induced causes of the microseismic events of the target mining area according to the magnitude distribution characteristics and the B value characteristics, and grade each microseismic event according to the separation result.

[0013] The fitting module is configured to determine whether each level of microseismic event is confused, extract a main fracture surface of a natural fracture according to the grading result of each level of microseismic event and the determination result of whether each level of microseismic event is confused, and fit the fracture surface of each microseismic event, respectively.

[0014] The fracture delineation module is configured to delineate the borehole annulus fractures of each microseismic event according to the fitting result to obtain a borehole annulus fracture model, and perform fracturing design on a horizontal well corresponding to each borehole annulus fracture model according to each borehole annulus fracture model.

[0015] The embodiment of the present application further provides a computer device, which comprises a memory, a processor and a computer program stored in the memory and capable of running on the processor, and the processor implements the above-mentioned B value-based microseismic natural fracture separation and fitting method when executing the computer program.

[0016] The embodiment of the present application further provides a computer readable storage medium, which stores a computer program, and the computer program implements the above-mentioned B value-based microseismic natural fracture separation and fitting method when executed by a processor.

[0017] The embodiment of the present application further provides a computer program product, which comprises a computer program, and the computer program implements the above-mentioned B value-based microseismic natural fracture separation and fitting method when executed by a processor.

[0018] In the embodiment of the present application, microseismic data of a target mining area is acquired, and magnitude distribution characteristics and B value characteristics are extracted according to the microseismic data; according to the magnitude distribution characteristics and the B value characteristics, each microseismic event is separated according to the inducing reasons of the microseismic events of the target mining area, and each microseismic event is classified according to the separation result; whether each microseismic event is confused is judged, and a main fracture surface of a natural fracture is extracted according to the classification result of each microseismic event and the judgment result of whether confusion occurs, and the fracture surfaces of each microseismic event are fitted respectively; the wellbore fractures of each microseismic event are depicted according to the fitting result, and a wellbore fracture model is obtained; and each wellbore fracture model is used to design fracturing of a horizontal well corresponding to the wellbore fracture model. In this way, the microseismic events are distinguished by the magnitude distribution characteristics and the B value characteristics of the microseismic data, the main fracture surface of the natural fracture is extracted, the wellbore fractures are finely depicted, and finally, the differential segmentation is implemented according to the fracture development characteristics on the basis of high-precision fracture prediction, the matrix modification degree is improved, the risk assessment is accurate, and the casing deformation risk is prevented. BRIEF DESCRIPTION OF DRAWINGS

[0019] In order to more clearly illustrate the technical solutions in the embodiments of the present application or the prior art, the following will briefly introduce the drawings needed to be used in the embodiments or the prior art description. Obviously, the drawings in the following description are only some embodiments of the present application, and other drawings can be obtained by those skilled in the art without any creative effort. In the drawings:

[0020] Figure 1 A flow chart of the microseismic natural fracture separation and fitting method based on B value provided in the embodiment of the present application;

[0021] Figure 2 A schematic diagram of the microseismic natural fracture separation and fitting method based on B value provided in the embodiment of the present application;

[0022] Figure 3 A schematic diagram of the natural fracture and artificial fracture separation result using the magnitude histogram method provided in the embodiment of the present application;

[0023] Figure 4 A schematic diagram of the natural fracture and artificial fracture separation result using the B value analysis method provided in the embodiment of the present application;

[0024] Figure 5 A schematic diagram of the reservoir matrix rupture type of the microseismic event provided in the embodiment of the present application;

[0025] Figure 6 A schematic diagram of the natural fracture extension type of the microseismic event provided in the embodiment of the present application;

[0026] Figure 7A schematic diagram of the microseismic event fault activation induced type provided in the embodiments of the present application;

[0027] Figure 8 A schematic diagram of the microseismic event point data interpretation provided in the embodiments of the present application;

[0028] Figure 9 A microseismic event point weighted three-dimensional plane fitting diagram provided in the embodiments of the present application;

[0029] Figure 10 A microseismic event point network reconstruction technology fitting fracture surface result diagram provided in the embodiments of the present application;

[0030] Figure 11 A microseismic constraint seismic fracture prediction diagram provided in the embodiments of the present application;

[0031] Figure 12 A schematic diagram of a microseismic natural fracture separation fitting device based on B value provided in the embodiments of the present application;

[0032] Figure 13 A structural block diagram of an electronic device provided in the embodiments of the present application. DETAILED DESCRIPTION

[0033] In order to make the objects, technical solutions and advantages of the embodiments of the present application clearer, further detailed description will be made to the embodiments of the present application in combination with the drawings. Herein, the schematic embodiments of the present application and the descriptions thereof are used to explain the present application, but not as a limitation to the present application.

[0034] The term “and / or” herein is merely used to describe an association relationship, which means that there can exist three relationships, for example, A and / or B can mean that A exists alone, A and B exist simultaneously, and B exists alone. In addition, the term “at least one” herein means any one of multiple or any combination of at least two of multiple, for example, including at least one of A, B and C can mean including any one or more elements selected from the set consisting of A, B and C.

[0035] In the description of the present specification, "include", "including", "have", "has", and the like are open terms, that is, mean including but not limited to. The description referring to the terms "one embodiment", "one specific embodiment", "some embodiments", "for example", and the like means that the specific features, structures, or characteristics described in connection with the embodiment or example are included in at least one embodiment or example of the present application. In the present specification, the illustrative description of the above terms does not necessarily refer to the same embodiment or example. Moreover, the specific features, structures, or characteristics described can be combined in any one or more embodiments or examples in a suitable manner. The order of steps involved in each embodiment is used to illustrate the implementation of the present application, and the order of steps is not limited and can be appropriately adjusted as needed.

[0036] The natural fracture development feature is an important index for adjusting the fracturing operation parameters, and is very important for development plan optimization. The microseismic fracture testing technology principle originates from the monitoring of natural earthquakes. In the hydraulic fracturing well, the formation generates a fracture, the formation around the fracture is stretched or dislocated, and the elastic wave energy is radiated outward, which is accepted by the receiver located beside the fracturing well. The space position of the receiver is obtained by processing the received signal, which represents the position of the fracture. However, the received information contains two kinds of fracture information, that is, the natural fracture and the artificial reconstruction fracture data. If the natural fracture development condition near the well group is obtained through processing, the potential channel communicated with the adjacent well before fracturing can be predicted, and corresponding measures are taken to adjust the fracturing parameters and operation scheme, so that the fracturing casing deformation and fracturing channeling can be effectively prevented, the development risk is reduced, the fracturing effect is improved, and the shale gas is effectively developed.

[0037] The deep shale reservoir develops multi-scale natural fractures. The prediction of the underground fracture by using the conventional seismic data is affected by the resolution, the prediction scale is large, and the demand of the fracturing engineering cannot be met. However, the microseismic data can be used as prior information to constrain the fracture density distribution and development direction predicted by using the seismic data, so that the relatively accurate fracture space distribution is obtained, which can provide effective guidance for the design of the horizontal well group and the optimization of the fracturing parameters of the block, so that the casing deformation and fracturing channeling risk can be greatly reduced. To realize the constraint of the microseismic data on the conventional seismic prediction fracture, the natural fracture separation fitting needs to be realized from the microseismic data.

[0038] Based on this, the embodiment of the present application provides a microseismic natural fracture separation fitting method based on B value, as shown in Figure 1 , which comprises the following steps:

[0039] Step 101: Obtain the microseismic data of the target mining area, and extract the magnitude distribution characteristics and B value characteristics according to the microseismic data;

[0040] Step 102: according to the magnitude distribution characteristics and the B value characteristics, separating each microseismic event according to the inducing reason of the microseismic event of the target mining area, and grading each microseismic event according to the separation result;

[0041] Step 103: judging whether each level of microseismic event is confused, extracting the main fracture surface of the natural fracture according to the grading result of each level of microseismic event and the judgment result of whether confusion occurs, and fitting the fracture surface of each microseismic event respectively;

[0042] Step 104: according to the fitting result, depicting the wellbore fracture of each microseismic event to obtain a wellbore fracture model;

[0043] Step 105: according to each wellbore fracture model, carrying out fracturing design on the horizontal well corresponding to each wellbore fracture model.

[0044] The microseismic natural fracture separation and fitting method based on B value provided by the embodiment of the application distinguishes the microseismic events induced by natural fractures from the microseismic events induced by the fracturing reservoir matrix by the magnitude distribution characteristics and the B value characteristics of the microseismic data, extracts the main fracture surface of the natural fracture by using the cluster analysis and the fault surface fitting technology, realizes the fine depiction of the wellbore fracture, and finally realizes the differential segmentation according to the fracture development characteristics on the basis of the high-precision fracture prediction, improves the matrix modification degree, accurately evaluates the risk, and prevents the casing deformation risk.

[0045] In an embodiment, the inducing reasons of the microseismic events include the microseismic events induced by natural fractures and the microseismic events induced by the fracturing reservoir matrix.

[0046] In an embodiment, according to the magnitude distribution characteristics and the B value characteristics, separating each microseismic event according to the inducing reason of the microseismic event of the target mining area includes:

[0047] According to the magnitude distribution characteristics, separating the microseismic events induced by natural fractures and the microseismic events induced by the fracturing reservoir matrix by the histogram method;

[0048] According to the B value characteristics, separating the microseismic events induced by natural fractures and the microseismic events induced by the fracturing reservoir matrix by the B value analysis method.

[0049] In an embodiment, the grading result of the microseismic event includes:

[0050] The microseismic events induced by the fracture of the reservoir matrix, the microseismic events induced by the fracture propagation of the fracturing through the natural fracture, and the microseismic events induced by the fault activation in the hydraulic fracturing process.

[0051] In one embodiment, based on the classification results of microseismic events at various levels and the judgment results of whether confusion has occurred, the main fracture surface of natural fractures is extracted, and the fracture surface of each microseismic event is fitted separately, including:

[0052] For microseismic events induced by reservoir matrix rupture, fracture imaging is performed using the fractured section as the basic interpretation unit.

[0053] Cluster analysis was performed on microseismic events induced by fracture propagation through natural fractures and microseismic events induced by fault activation during hydraulic fracturing. Fracture clusters were used as the basic interpretive units for fracture imaging.

[0054] For microseismic events that are not confused, a weighted three-dimensional plane fitting is performed based on the spatial location and magnitude of each fracture segment as a natural fracture seismic event to obtain the main fracture surface;

[0055] For microseismic events that are confused, cluster analysis and weighted three-dimensional plane fitting are performed on the microseismic events to obtain the main fault surface.

[0056] In practice, microseismic events induced by natural fractures are distinguished from those induced by fracturing the reservoir matrix by analyzing the magnitude distribution and B-value characteristics of microseismic data. Based on the separation results, microseismic events are classified into three types. Depending on the type and whether microseismic event points are confused, cluster analysis and fault plane fitting techniques are used to extract the main fracture surface of natural fractures. Different fitting techniques are employed based on whether microseismic event points are confused, achieving a detailed characterization of wellbore fractures. Based on the fracture model, differentiated fracturing designs are implemented for specific horizontal wells, effectively reducing casing deformation and hydraulic channeling, and significantly improving the efficiency of shale gas exploration and development.

[0057] For example, through the processing and analysis of microseismic data in the Zigong East 3D area of ​​the Sichuan shale gas region, it was found that the fracture prediction results using the step-by-step fracture identification technology showed that the fractures in the Zigong East block were mainly distributed in a northeast direction. However, there was a significant difference between the microseismic monitoring and prediction results from the 215H1 platform. After constraining with microseismic information, the fracture distribution characteristics and fracturing results tended to be consistent. Based on the fracture orientation model of the study area, the spatial distribution characteristics of fault and fracture development zones can be effectively predicted. At the same time, based on the fracture model, differentiated fracturing designs can be carried out for specific horizontal wells, which can reduce casing deformation and pressure channeling, and significantly improve the efficiency of shale gas exploration and development.

[0058] like Figure 2 As shown in the flowchart of the microseismic natural fracture separation and fitting technology based on the B-value, the main core contents include:

[0059] (1) Natural fracture and artificial fracture separation technology: the number of microseismic events and their magnitude have a power law relationship. Studies show that the standard induced fractures that occur in standard shale often correspond to a B value of about 2. If the events that occur during a period of construction correspond to existing fractures or faults, the B value will quickly drop to about 1. The frequency-magnitude relationship of any seismic sequence can be written as:

[0060]

[0061] where N M represents the cumulative number of earthquakes or events with a magnitude greater than or equal to M. A and B are constants. A is the intercept, and B represents the slope of the straight line.

[0062] Separation of microseismic event points of natural fractures and artificial fractures by magnitude mainly includes two methods:

[0063] Histogram method: using numerical mutation point extraction, as shown in Figure 3 .

[0064] B value analysis method: according to the slope of the magnitude-frequency logarithm, as shown in Figure 4 .

[0065] (2) Based on the magnitude and B value to distinguish three types of microseismic events:

[0066] Reservoir matrix fracture type: magnitude -3 to -1, B value >1.5, fracturing microseismic events are affected by small-scale natural fractures, as shown in Figure 5 .

[0067] Natural fracture extension type: magnitude -1 to 2, B value 0.5 to 1.5, large-scale natural fractures have a strong influence, and there is no fault activation, as shown in Figure 6 .

[0068] Fault activation induced type: magnitude 1 to 2 or higher; B value <0.5, fault is easy to activate, as shown in Figure 7 .

[0069] (3) Natural fracture fracture surface fitting technology, as shown in the technical process, mainly includes two cases: Figure 8

[0070] One is the microseismic event point without confusion, according to the normalized spatial position and magnitude of each fracturing section micro natural fracture seismic event point, weighted three-dimensional plane fitting is carried out to obtain the main fracture surface, and the length, area, inclination degree and other occurrence information of the main fracture surface of each fracturing section are analyzed, as shown in Figure 9 .

[0071] ​Another is the microseismic event point confusion, first microseismic event point clustering analysis, and then weighted three-dimensional plane fitting. Often in actual situation, there may be multiple groups of different fracture groups in the fracture network in the fracturing section, therefore, only simply extracting the main fault plane cannot effectively fracture the fracture network structure in the section; the application adopts a robust three-dimensional fracture network reconstruction method (RFM3D) such as Figure 10 as shown.

[0072] Figure 11 As shown, there is a significant difference between the microseismic monitoring and prediction results of the 215H1 platform, and the fracture distribution characteristics and the fracturing results tend to be consistent after the microseismic information constraint. According to the fracture azimuth model of the research area, the spatial distribution characteristics of the fault and the fracture development zone can be effectively predicted. At the same time, based on high-precision fracture prediction, according to the fracture development characteristics, differential segmentation is implemented, high-strength multi-cluster fracturing process, multi-cycle staged regulation and control fracturing process and deep penetration fracturing process and other differential fracturing are adopted, and the matrix fullness and casing deformation and channeling prevention and control are organically integrated.

[0073] The embodiment of the application also provides a B-value-based microseismic natural fracture separation fitting device, as described in the following embodiment. Since the principle of solving the problem of the device is similar to that of the B-value-based microseismic natural fracture separation fitting method, the implementation of the device can be referred to the implementation of the method, and the repeated parts will not be described here.

[0074] Figure 12 A schematic diagram of the B-value-based microseismic natural fracture separation fitting device provided in the embodiment of the application is shown in Figure 12 , and the device comprises:

[0075] The data acquisition module 1201 is configured to acquire microseismic data of a target mining area, and extract magnitude distribution characteristics and B-value characteristics according to the microseismic data;

[0076] The separation module 1202 is configured to separate each microseismic event according to the magnitude distribution characteristics and the B-value characteristics according to the induced reasons of the microseismic events of the target mining area, and grade each microseismic event according to the separation result;

[0077] The fitting module 1203 is configured to judge whether each level of microseismic event is confused, extract a natural fracture main fault plane according to the grading result of each level of microseismic event and the judgment result of whether the confusion occurs, and fit the fault plane of each microseismic event respectively;

[0078] The crack characterization module 1204 is configured to characterize the wellbore cracks of each microseismic event according to the fitting result to obtain a wellbore crack model; and perform fracturing design on the horizontal well corresponding to each wellbore crack model according to the wellbore crack model.

[0079] In an embodiment, the microseismic event inducing reasons include natural crack induced microseismic events and fracturing reservoir matrix induced microseismic events.

[0080] In an embodiment, the separation module 1202 is specifically configured to:

[0081] According to the magnitude distribution characteristics, the natural crack induced microseismic events and the fracturing reservoir matrix induced microseismic events are separated by a histogram method;

[0082] According to the B value characteristics, the natural crack induced microseismic events and the fracturing reservoir matrix induced microseismic events are separated by a B value analysis method.

[0083] In an embodiment, the microseismic event classification result includes:

[0084] Reservoir matrix rupture induced microseismic events, fracturing crack propagation induced microseismic events through natural cracks, and microseismic events induced by fault activation in the hydraulic fracturing process.

[0085] In an embodiment, the fitting module 1203 is specifically configured to:

[0086] For the reservoir matrix rupture induced microseismic events, the fracture imaging is performed by taking the fracturing section as a basic interpretation unit;

[0087] For the fracturing crack propagation induced microseismic events through natural cracks and the microseismic events induced by fault activation in the hydraulic fracturing process, the cluster analysis is performed, and the fracture imaging is performed by taking the fracture cluster as a basic interpretation unit;

[0088] For the microseismic events without confusion, the weighted three-dimensional plane fitting is performed according to the spatial position and magnitude of each fracturing section for the natural crack seismic event to obtain a main fracture surface;

[0089] For the microseismic events with confusion, the cluster analysis and the weighted three-dimensional plane fitting are performed on the microseismic events to obtain a main fracture surface.

[0090] Based on the foregoing inventive concept, as shown in Figure 13 The application further provides a computer device 1300, which comprises a memory 1310, a processor 1320, and a computer program 1330 stored in the memory 1310 and capable of running on the processor 1320, wherein the processor 1320 implements the foregoing B value based microseismic natural crack separation fitting method when executing the computer program 1330.

[0091] The embodiment of the present application further provides a computer readable storage medium, which stores a computer program, and the computer program is executed by a processor to realize the B-value-based microseismic natural fracture separation fitting method.

[0092] The embodiment of the present application further provides a computer program product, which comprises a computer program, and the computer program is executed by a processor to realize the B-value-based microseismic natural fracture separation fitting method.

[0093] To sum up, in the embodiment of the present application, the microseismic data of a target mining area is acquired, the magnitude distribution characteristics and B-value characteristics are extracted according to the microseismic data, according to the magnitude distribution characteristics and B-value characteristics, each microseismic event is separated according to the inducing reasons of the microseismic events of the target mining area, and each microseismic event is classified according to the separation result, whether the microseismic events are confused is judged, the main fracture surface of the natural fracture is extracted according to the classification result of each microseismic event and the judgment result of whether the microseismic events are confused, and the fracture surfaces of each microseismic event are fitted respectively, the wellbore fractures of each microseismic event are depicted according to the fitting result, and the wellbore fracture model is obtained, and the horizontal wells corresponding to each wellbore fracture model are fractured according to each wellbore fracture model. In this way, the microseismic events are distinguished through the magnitude distribution characteristics and B-value characteristics of the microseismic data, the main fracture surface of the natural fracture is extracted, the wellbore fractures are finely depicted, and finally, the differential segmentation is implemented according to the fracture development characteristics on the basis of high-precision fracture prediction, the matrix modification degree is improved, the risk evaluation is accurate, and the casing deformation risk is prevented.

[0094] Those skilled in the art will understand that the embodiments of the present application can be provided as methods, systems, or computer program products. Therefore, the present application can adopt a completely hardware embodiment, a completely software embodiment, or an embodiment combining software and hardware aspects. Moreover, the present application can adopt the form of a computer program product implemented on one or more computer usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer usable program codes.

[0095] The present application is described with reference to flowcharts and / or block diagrams of the methods, devices (systems), and computer program products according to the embodiments of the present application. It should be understood that each flow and / or block in the flowcharts and / or block diagrams, and the combination of the flows and / or blocks in the flowcharts and / or block diagrams can be implemented by computer program instructions. These computer program instructions can be provided to the processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing devices to produce a machine, so that the instructions executed by the processor of the computer or other programmable data processing devices produce a device that implements the flowcharts and / or block diagrams. Figure 1one or more processes and / or blocks Figure 1 an apparatus for performing the functions specified in the flowchart or multiple flows and / or blocks.

[0096] These computer program instructions can also be stored in a computer readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer readable memory produce an article of manufacture including instructions which implement the flowchart or multiple flows and / or blocks. Figure 1 one or more processes and / or blocks Figure 1 an apparatus for performing the functions specified in the flowchart or multiple flows and / or blocks.

[0097] These computer program instructions can also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the flowchart or multiple flows and / or blocks. Figure 1 one or more processes and / or blocks Figure 1 an apparatus for performing the functions specified in the flowchart or multiple flows and / or blocks.

[0098] The above-described specific embodiments, the purpose, technical solutions and beneficial effects of the present application are further described in detail, it should be understood that the above-described is only the specific embodiments of the present application, and is not used to limit the protection scope of the present application, any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present application should be included in the protection scope of the present application.

Claims

1. A B-value based microseismic natural fracture separation fitting method, characterized in that, The method comprises the following steps: obtaining microseismic data of a target mining area, and extracting magnitude distribution characteristics and B value characteristics from the microseismic data; separating each microseismic event according to the induced causes of the microseismic events of the target mining area according to the magnitude distribution characteristics and the B value characteristics, and grading each microseismic event according to the separation result; determining whether each level of microseismic event is confused, and fitting the fracture surface of each microseismic event according to the grading result of each level of microseismic event and the determination result of whether the microseismic event is confused to extract the main fracture surface of the natural fracture; drawing each wellbore crack according to the fitting result to obtain a wellbore crack model; designing fracturing of each wellbore corresponding to each wellbore crack model according to each wellbore crack model.

2. The method of claim 1, wherein, The induced causes of the microseismic events include natural fracture induced microseismic events and fracturing reservoir matrix induced microseismic events.

3. The method of claim 2, wherein, The separating each microseismic event according to the induced causes of the microseismic events of the target mining area according to the magnitude distribution characteristics and the B value characteristics comprises: separating the natural fracture induced microseismic events and the fracturing reservoir matrix induced microseismic events by histogram method according to the magnitude distribution characteristics; separating the natural fracture induced microseismic events and the fracturing reservoir matrix induced microseismic events by B value analysis method according to the B value characteristics.

4. The method of claim 1, wherein, The grading result of the microseismic events comprises: reservoir matrix rupture induced microseismic events, fracturing through natural fracture induced microseismic events, and microseismic events induced by fault activation in the process of hydraulic fracturing.

5. The method of claim 4, wherein, The fitting the fracture surface of each microseismic event according to the grading result of each level of microseismic event and the determination result of whether the microseismic event is confused to extract the main fracture surface of the natural fracture comprises: for the reservoir matrix rupture induced microseismic events, taking the fracturing section as a basic interpretation unit to perform crack imaging; for the fracturing through natural fracture induced microseismic events and the microseismic events induced by fault activation in the process of hydraulic fracturing, performing cluster analysis to take the fracture cluster as a basic interpretation unit to perform crack imaging; for the microseismic events without confusion, performing weighted three-dimensional plane fitting according to the spatial position and magnitude of each fracturing section to obtain the main fracture surface; for the confused microseismic events, performing cluster analysis and weighted three-dimensional plane fitting on the microseismic events to obtain the main fracture surface.

6. A B-value based microseismic natural fracture separation fitting device, characterized in that, The method comprises the following steps: a data acquisition module is configured to obtain microseismic data of a target mining area, and extract magnitude distribution characteristics and B value characteristics from the microseismic data; a separation module is configured to separate each microseismic event according to the induced causes of the microseismic events of the target mining area according to the magnitude distribution characteristics and the B value characteristics, and grade each microseismic event according to the separation result; a fitting module is configured to determine whether each level of microseismic event is confused, and fit the fracture surface of each microseismic event according to the grading result of each level of microseismic event and the determination result of whether the microseismic event is confused to extract the main fracture surface of the natural fracture. The fracture delineation module is configured to delineate the wellbore fractures of each microseismic event according to the fitting result to obtain a wellbore fracture model; and to perform fracturing design on the horizontal well corresponding to each wellbore fracture model according to the wellbore fracture model.

7. The apparatus of claim 6, wherein, The microseismic event inducing reasons include natural fracture induced microseismic event and fracturing reservoir matrix induced microseismic event.

8. The apparatus of claim 7, wherein, The separation module is specifically configured to: According to the magnitude distribution characteristics, the natural fracture induced microseismic event and the fracturing reservoir matrix induced microseismic event are separated by a histogram method; According to the B value characteristics, the natural fracture induced microseismic event and the fracturing reservoir matrix induced microseismic event are separated by a B value analysis method.

9. The apparatus of claim 6, wherein, The microseismic event classification result includes: Reservoir matrix rupture induced microseismic event, fracturing through natural fracture induced microseismic event and fault activation induced microseismic event in the hydraulic fracturing process.

10. The apparatus of claim 9, wherein, The fitting module is specifically configured to: For the reservoir matrix rupture induced microseismic event, the fracture imaging is performed by taking the fracturing section as a basic interpretation unit; For the fracturing through natural fracture induced microseismic event and the fault activation induced microseismic event in the hydraulic fracturing process, the cluster analysis is performed, and the fracture imaging is performed by taking the fracture cluster as a basic interpretation unit; For the microseismic event without confusion, the weighted three-dimensional plane fitting is performed according to the spatial position and magnitude of the natural fracture event of each fracturing section to obtain a main fracture surface; For the confused microseismic event, the cluster analysis and the weighted three-dimensional plane fitting are performed on the microseismic event to obtain a main fracture surface.

11. A computer device comprising a memory, a processor, and a computer program stored on the memory and executable on the processor, characterized in that, The processor executes the computer program to realize the method of any one of claims 1 to 5.

12. A computer-readable storage medium, characterized in that, The computer readable storage medium stores a computer program, and the computer program is executed by the processor to realize the method of any one of claims 1 to 5.

13. A computer program product, characterised in that, The computer program product includes a computer program, and the computer program is executed by the processor to realize the method of any one of claims 1 to 5.

Images