Construction method and device of breaking joint body crack density three-dimensional model and electronic equipment

Abstract

The invention relates to a method and device for constructing a three-dimensional fracture density model of a breaking joint body and electronic equipment, and belongs to the technical field of oil reservoir exploration and development, and the method comprises the steps: obtaining a field outcrop image of a target area; marking cracks on the field outcrop image of the target area, and obtaining the density of the cracks; obtaining a scatter diagram based on the density of the crack and the first distance between the crack and the fault; fitting data on the scatter diagram to obtain a fitting function; constructing a geologic model of the target area; calculating a second distance from a grid in the geologic model to a fault main shaft surface; and based on the fitting function and the second distance, obtaining a three-dimensional model of the fracture density of the breaking joint body in the target area. According to the method, the fracture density of the breaking joint body can be depicted with higher precision, so that the recovery rate of oil and gas is improved.

Description

Technical Field

[0001] The present invention relates to the technical field of reservoir exploration and development, and in particular to a method, device and electronic device for constructing a three-dimensional model of fracture density of a fracture body. Background Technique

[0002] Faults are one of the main factors for generating fractures. The differences in faults determine various parameters such as the scale, density, and aperture of fracture development, and further determine the development degree of its reservoir, which has an important impact on the migration, accumulation and recovery rate of oil and gas, and determines the scale of oil and gas reserves and economic benefits.

[0003] The quantitative analysis of fracture density of a fracture body is relatively difficult. Currently, there are mainly two methods. One is to use indirect data such as logging and seismic data and use the geological model of fracture development as a constraint to predict fracture density; the other is to analyze and infer the fracture density by directly using data such as outcrops and cores to analyze the distance from the fault. The former is limited by the resolution and can only depict large-scale fractures, and cannot accurately depict small-scale fractures; the latter believes that the fracture density decreases linearly with the distance from the fault, which is quite different from the actual situation of the fluctuating decrease of the outcrop fracture density in the field.

[0004] In summary, there is a lack of a modeling method in the prior art that can depict the fracture density of a fracture body with higher accuracy. Summary of the Invention

[0005] In view of this, it is necessary to provide a method, device and electronic device for constructing a three-dimensional model of fracture density of a fracture body to achieve the purpose of a modeling method for depicting the fracture density of a fracture body with higher accuracy.

[0006] To achieve the above object, the present invention provides a method for constructing a three-dimensional model of fracture density of a fracture body, including: Obtain the outcrop image of the target area in the field; Mark the fractures on the outcrop image of the target area in the field and obtain the density of the fractures; Based on the density of the fractures and the first distance from the fractures to the fault, obtain a scatter plot; Fit the data on the scatter plot to obtain a fitting function; Construct a geological model of the target area; Calculate the second distance from the grid in the geological model to the main fault plane; Based on the fitting function and the second distance, obtain a three-dimensional model of fracture density of the fracture body in the target area.

[0007] In a possible implementation manner, the fitting the data on the scatter plot to obtain a fitting function includes: Fitting the data on the scatter plot using data processing and fitting software to obtain a fitting function.

[0008] In one possible implementation, constructing the geological model of the target area includes: Constructing the geological model of the target area based on 3D geological model software.

[0009] In one possible implementation, the 3D geological model software includes one or more of PETREL, GOCAD, 3DMine, and PLAXIS 3D.

[0010] In one possible implementation, calculating the second distance from the grid in the geological model to the main fault plane includes: Using each main fault plane as the starting plane and a grid of a preset size as the basic unit; Calculating the distance from the grid in the model to the main fault plane based on the starting plane and the basic unit.

[0011] In one possible implementation, the expression of the fitting function is:

[0012] In the formula, represents the density of fractures, with the unit of m -1 ; represents the distance from the fault, with the unit of m; A represents the vibration amplitude of fracture development; k represents the change rate of the attenuation degree of fracture development; ω represents the vibration frequency of fracture development; b represents the reference value of fracture development.

[0013] In one possible implementation, obtaining the 3D model of the fracture density of the fault fracture body based on the fitting function and the second distance includes: Substituting the input second distance into the fitting function to obtain the fracture density; Obtaining the 3D model of the fracture density of the fault fracture body in the target area based on the fracture density.

[0014] On the other hand, the present invention also provides a device for constructing a 3D model of the fracture density of a fault fracture body, including: An image acquisition module for acquiring field outcrop images of the target area; A fracture line density acquisition module for marking fractures on the field outcrop images of the target area and acquiring the density of the fractures; A scatter plot acquisition module for obtaining a scatter plot based on the density of the fractures and the first distance from the fractures to the fault. A fitting function acquisition module, configured to fit the data on the scatter plot to obtain a fitting function; A geological model acquisition module, configured to construct a geological model of the target area; A second distance acquisition module, configured to calculate a second distance from the grid in the geological model to the main fault plane; A fracture density three-dimensional model acquisition module, configured to obtain a fracture density three-dimensional model of the fault fracture body in the target area based on the fitting function and the second distance.

[0015] On the other hand, the present invention further provides an electronic device, including a memory and a processor, wherein, The memory is configured to store a program; The processor is coupled to the memory and configured to execute the program stored in the memory to implement the steps in a method for constructing a fracture density three-dimensional model of a fault fracture body as described in any one of the above implementation manners.

[0016] On the other hand, the present invention further provides a computer-readable storage medium, configured to store a computer-readable program or instruction, and when the program or instruction is executed by a processor, it can implement the steps in a method for constructing a fracture density three-dimensional model of a fault fracture body as described in any one of the above implementation manners.

[0017] The beneficial effects of the present invention are as follows: A method for constructing a fracture density three-dimensional model of a fault fracture body provided by the present invention first obtains a field outcrop image of a target area, then marks fractures on the outcrop image of the target area and obtains the density of the fractures. Further, based on the density of the fractures and the first distance from the fractures to the fault, a scatter plot is obtained. It is determined through the scatter plot that the fracture density shows a fluctuating characteristic as the distance from the fault increases. Then, the data on the scatter plot is further fitted to obtain a fitting function. Then, a geological model of the target area is constructed, and the second distance from the grid in the geological model to the main fault plane is calculated. Finally, based on the fitting function and the second distance, a fracture density three-dimensional model of the fault fracture body in the target area is obtained. The present invention quantitatively predicts the fracture density of the fault fracture body through a scatter plot, establishes a fracture density three-dimensional model of the fault fracture body, and quantitatively analyzes the fracture density through the scatter plot to more accurately characterize the fracture density of the fault fracture body. Description of the Drawings

[0018] Figure 1 It is a flowchart of a method for an embodiment of a method for constructing a fracture density three-dimensional model of a fault fracture body provided by the present invention; Figure 2 It is a flowchart of a method for a specific implementation in a method for constructing a fracture density three-dimensional model of a fault fracture body provided by the present invention; Figure 3The outcrop fracture density interpretation map in the method for constructing a three-dimensional model of the fracture density of a faulted seam body provided by the present invention; Figure 4 The outcrop fracture density scatter plot in the method for constructing a three-dimensional model of the fracture density of a faulted seam body provided by the present invention; Figure 5 The geological model of the target area in the method for constructing a three-dimensional model of the fracture density of a faulted seam body provided by the present invention; Figure 6 The three-dimensional model map of the outcrop fracture density of the target area in the method for constructing a three-dimensional model of the fracture density of a faulted seam body provided by the present invention; Figure 7 The schematic flowchart of an embodiment of the device for constructing a three-dimensional model of the fracture density of a faulted seam body provided by the present invention; Figure 8 The schematic structural diagram of an embodiment of the electronic device provided by the present invention. Detailed implementation manners

[0019] Next, the preferred embodiments of the present invention will be specifically described with reference to the accompanying drawings. The accompanying drawings form a part of this application and are used together with the embodiments of the present invention to explain the principle of the present invention, rather than to limit the scope of the present invention.

[0020] The present invention provides a method, device and electronic device for constructing a three-dimensional model of the fracture density of a faulted seam body, which will be described separately below.

[0021] Figure 1 The schematic flowchart of an embodiment of the method for constructing a three-dimensional model of the fracture density of a faulted seam body provided by the present invention. As Figure 1 shown, the method for constructing a three-dimensional model of the fracture density of a faulted seam body includes: S101. Obtain the outcrop image of the target area; S102. Mark the fractures on the outcrop image of the target area and obtain the fracture line density of the fractures; S103. Based on the fracture line density and the first distance from the fracture to the fault, obtain a scatter plot; S104. Fit the data on the scatter plot to obtain a fitting function; S105. Construct the geological model of the target area; S106. Calculate the second distance from the grid in the geological model to the main fault plane; S107. Based on the fitting function and the second distance, obtain the three-dimensional model of the fracture density of the faulted seam body in the target area.

[0022] Compared with the prior art, a method for constructing a three-dimensional model of fracture density of a fracture body provided in this embodiment first obtains a field outcrop image of a target area, then marks fractures on the field outcrop image of the target area and obtains the density of the fractures. Further, based on the density of the fractures and the first distance from the fractures to the fault, a scatter plot is obtained. It is determined through the scatter plot that the fracture density shows a fluctuating characteristic as the distance from the fault increases. Then, the data on the scatter plot is further fitted to obtain a fitting function. Then, a geological model of the target area is constructed, and the second distance from the grid in the geological model to the main fault plane is calculated. Finally, based on the fitting function and the second distance, a three-dimensional model of fracture density of the fracture body in the target area is obtained. Through the quantitative prediction of the fracture density of the fracture body by the scatter plot, the present invention establishes a three-dimensional model of fracture density of the fracture body, and quantitatively analyzes the fracture density through the scatter plot, so as to more accurately depict the fracture density of the fracture body.

[0023] It should be noted that the method of the present invention is carried out on a computer terminal device.

[0024] It should be noted that there are two display forms of the fracture body. One is the natural outcrop in the field; the other is the fracture body reservoir formed by underground fractures.

[0025] In a specific embodiment of the present invention, as Figure 2 shown, the method for constructing a three-dimensional model of fracture density of the fracture body specifically includes the following steps: Step 1, construct a high-precision fracture body field outcrop in the target area by using a drone and a high-precision three-dimensional laser scanner, and carry out fracture interpretation and mark the fractures.

[0026] Step 2, count the density of the fractures and make a scatter plot of the fracture density and the distance from the fault.

[0027] Step 3, according to the fact that the fracture density obtained from the fracture density scatter plot in Step 2 shows a fluctuating decreasing characteristic of every 20m - 30m as the distance from the fault increases, determine the preferred unilateral exponential decay cosine function and parameters to ensure its consistency with the outcrop fracture density distribution; Step 4, establish a work area in the target area, import the layer data and fault data of the target area, and use the traditional corner point grid modeling method to jointly establish a structural geological model of the fault and the layer grid.

[0028] Step 5, through the established fault model, take each main fault plane as the starting plane for distance calculation, and take a 2m 2m 1m grid as the basic unit to calculate the distance from each grid in the model to the main fault plane; Step 6: Based on the one-sided exponential decay cosine function and parameters optimized in Step 3, and on the constructed model established in Step 4 and the distances from each grid to the fault calculated in Step 5, establish a three-dimensional model of the fracture density of the fracture body.

[0029] In some embodiments of the present invention, in Step S101, obtaining the outcrop image of the target area specifically includes: Construct a high-precision outcrop of the fracture body in the target area by using an unmanned aerial vehicle and a high-precision three-dimensional laser scanner, and conduct fracture interpretation and mark the fractures.

[0030] In some embodiments of the present invention, in Step S102, mark the fractures on the outcrop image of the target area and obtain the fracture line density of the fractures. As Figure 3 shown, mark the fractures on the outcrop image of the target area to obtain the fracture line density of the fractures.

[0031] As Figure 3 shown, the fracture density in this target area has obvious regularity, that is, the right side is closer to the fault, and the fracture density is larger; the left side is farther from the fault, and the fracture density is smaller, and the fracture density shows a fluctuating decreasing characteristic with every increase of 20m - 30m in the distance from the fault. And based on this, it is divided into a fracture zone with a larger fracture density and a fracture zone with a smaller fracture density. Therefore, according to the scatter plot, it is determined that the fracture density shows a fluctuating decreasing characteristic with every increase of 20m - 30m in the distance from the fault, and the one-sided exponential decay cosine function and parameters are optimized to ensure that it is consistent with the outcrop fracture density distribution. In some embodiments of the present invention, the fitting of the data on the scatter plot to obtain a fitting function includes: Based on data processing and fitting software, fit the data on the scatter plot to obtain a fitting function.

[0032] In some embodiments of the present invention, Figure 4 is the scatter plot of the outcrop fracture density according to the public example. As Figure 4 shown, taking the measured fracture density as the standard, according to the actual situation that only the fracture zone or the fracture zone with a larger fracture density contains exploitable oil and gas during the oil and gas exploration process, select the fracture densities of the first 11 points to optimize different functions. And by comparing the original fitting function with the optimized fitting function this time, it is found that the R 2 of the optimized function is 0.63, while the R 2 of the original optimized function is only 0.42. Thus, it is determined that the one-sided exponential decay cosine function is used as the fitting function for the relationship between the fracture density of the fracture body and the distance from the fault. The expression of the fitting function is:

[0033] In the formula, represents the density of the fractures, and the unit ism -1 ; represents the distance from the fault, with the unit of m; A represents the vibration amplitude of fracture development; k represents the change rate of attenuation degree of fracture development; ω represents the vibration frequency of fracture development; b represents the reference value of fracture development.

[0034] In some embodiments of the present invention, the construction of the geological model of the target area includes: Constructing the geological model of the target area based on 3D geological model software.

[0035] Specifically, establishing the work area of the target area, importing the layer data and fault data of the target area, and using the traditional corner point grid modeling method to jointly establish the structural geological model of the fault and the layer grid, as Figure 5 shown as the geological model.

[0036] In some embodiments of the present invention, the 3D geological model software includes one or more of PETREL, GOCAD, 3DMine, and PLAXIS 3D.

[0037] In some embodiments of the present invention, the calculation of the second distance from the grid in the geological model to the main fault plane includes: Taking each main fault plane as the starting plane and using a grid with a preset size as the basic unit; Based on the starting plane and the basic unit, calculating the distance from the grid in the model to the main fault plane.

[0038] Specifically, based on the geological model, taking each main fault plane as the starting plane for distance calculation and using a 2m * 2m * 1m grid as the basic unit, calculating the distance from each grid in the model to the main fault plane.

[0039] In some embodiments of the present invention, as Figure 6 shown as the 3D model of fracture density of the fault - fracture body in the target area, the obtaining of the 3D model of fracture density of the fault - fracture body in the target area based on the fitting function and the second distance includes: Substituting the input second distance into the fitting function to obtain the fracture density; Based on the fracture density, obtaining the 3D model of fracture density of the fault - fracture body in the target area.

[0040] In order to better implement a method for constructing a 3D model of fracture density of a fault - fracture body in an embodiment of the present invention, on the basis of a method for constructing a 3D model of fracture density of a fault - fracture body, correspondingly, as Figure 7As shown in the figure, an embodiment of the present invention further provides a device for constructing a three-dimensional model of the fracture density of a fracture body. A device 700 for constructing a three-dimensional model of the fracture density of a fracture body includes: An image acquisition module 701, configured to acquire an outcrop image of a target area in the wild; A fracture line density acquisition module 702, configured to mark fractures on the outcrop image of the target area in the wild and acquire the density of the fractures; A scatter plot acquisition module 703, configured to obtain a scatter plot based on the density of the fractures and the first distance from the fractures to the fault; A fitting function acquisition module 704, configured to fit the data on the scatter plot to obtain a fitting function; A geological model acquisition module 705, configured to construct a geological model of the target area; A second distance acquisition module 706, configured to calculate the second distance from the grid in the geological model to the main fault plane; A three-dimensional fracture density model acquisition module 707, configured to obtain a three-dimensional model of the fracture density of the fracture body in the target area based on the fitting function and the second distance.

[0041] The device 700 for constructing a three-dimensional model of the fracture density of a fracture body provided in the above embodiment can implement the technical solutions described in the above embodiment of the method for constructing a three-dimensional model of the fracture density of a fracture body. The specific implementation principles of the above modules or units can be referred to the corresponding content in the above embodiment of the method for constructing a three-dimensional model of the fracture density of a fracture body, and will not be elaborated here.

[0042] As Figure 8 shown in the figure, the present invention also correspondingly provides an electronic device 800. The electronic device 800 includes a processor 801, a memory 802, and a display 803. Figure 8 Only some components of the electronic device 800 are shown, but it should be understood that it is not required to implement all the shown components, and more or fewer components can be alternatively implemented.

[0043] In some embodiments, the processor 801 may be a central processing unit (CPU), a microprocessor, or other data processing chips, and is configured to run the program code stored in the memory 802 or process data, such as the method for constructing a three-dimensional model of the fracture density of a fracture body in the present invention.

[0044] In some embodiments, the processor 801 may be a single server or a server group. The server group may be centralized or distributed. In some embodiments, the processor 801 may be local or remote. In some embodiments, the processor 801 may be implemented on a cloud platform. In some embodiments, the cloud platform may include a private cloud, a public cloud, a hybrid cloud, a community cloud, a distributed cloud, an internal cloud, a multi-cloud, etc., or any combination thereof.

[0045] The memory 802 may be an internal storage unit of the electronic device 800 in some embodiments, such as the hard disk or memory of the electronic device 800. The memory 802 may also be an external storage device of the electronic device 800 in some other embodiments, such as a plug-in hard disk equipped on the electronic device 800, a Smart Media Card (SMC), a Secure Digital (SD) card, a Flash Card, etc.

[0046] Furthermore, the memory 802 may include both the internal storage unit and the external storage device of the electronic device 800. The memory 802 is used to store the application software installed in the electronic device 800 and various types of data.

[0047] The display 803 may be an LED display, a liquid crystal display, a touch liquid crystal display, and an OLED (Organic Light-Emitting Diode) toucher, etc. in some embodiments. The display 803 is used to display the information in the electronic device 800 and to display a visual user interface. The components 801 - 803 of the electronic device 800 communicate with each other through a system bus.

[0048] In one embodiment, when the processor 801 executes a program for constructing a three-dimensional model of the fracture density of a faulted body in the memory 802, the following steps may be implemented: Obtain field outcrop images of the target area; Mark fractures on the outcrop images of the target area and obtain the density of the fractures; Based on the density of the fractures and the first distance from the fractures to the fault, obtain a scatter plot; Fit the data on the scatter plot to obtain a fitting function; Construct a geological model of the target area; Calculate the second distance from the grid in the geological model to the main fault plane; Based on the fitting function and the second distance, obtain a three-dimensional model of the fracture density of the faulted body in the target area.

[0049] It should be understood that when the processor 801 executes the program for constructing a three-dimensional model of the crack density of a crack body in the memory 802, in addition to the above functions, other functions can also be realized. For specific details, reference can be made to the description of the corresponding method embodiments above.

[0050] Furthermore, the type of the electronic device 800 mentioned in the embodiments of the present invention is not specifically limited. The electronic device 800 can be a mobile phone, a tablet computer, a personal digital assistant (PDA), a wearable device, a laptop computer, or other portable electronic devices. Exemplary embodiments of the portable electronic device include, but are not limited to, portable electronic devices running IOS, android, microsoft, or other operating systems. The above-mentioned portable electronic devices can also be other portable electronic devices, such as a laptop computer with a touch-sensitive surface (such as a touch panel). It should also be understood that in some other embodiments of the present invention, the electronic device 800 may not be a portable electronic device, but a desktop computer with a touch-sensitive surface (such as a touch panel).

[0051] Those skilled in the art can understand that all or part of the processes of implementing the above method embodiments can be completed by instructing relevant hardware through a computer program, and the program can be stored in a computer-readable storage medium. Among them, the computer-readable storage medium is a magnetic disk, an optical disk, a read-only memory, or a random access memory, etc.

[0052] The above is only a preferred specific embodiment of the present invention, but the protection scope of the present invention is not limited thereto. Any changes or substitutions that can be easily thought of by those skilled in the art within the technical scope disclosed by the present invention should be covered by the protection scope of the present invention.

Claims

1. A method for constructing a three-dimensional model of fracture density of a fracture body, characterized in that: include: Acquire field outcrop images of target areas; Mark the fractures on the field outcrop images of the target area and obtain the density of the fractures; Obtaining a scatter plot based on the density of the cracks and a first distance between the cracks and the fault; Fitting the data on the scatter plot to obtain a fitting function; Constructing a geological model of the target area; Calculating a second distance from a grid in the geological model to a main axis of the fault; Based on the fitting function and the second distance, a three-dimensional model of crack density of the fault body in the target area is obtained.

2. The method for constructing a three-dimensional model of fracture density of a fracture body according to claim 1, characterized in that: The step of fitting the data on the scatter plot to obtain a fitting function includes: The data on the scatter plot are fitted based on data processing and fitting software to obtain a fitting function.

3. The method for constructing a three-dimensional model of fracture density of a fracture body according to claim 1, characterized in that: The step of constructing a geological model of the target area includes: Construct a geological model of the target area based on three-dimensional geological model software.

4. The method for constructing a three-dimensional model of fracture density of a fracture body according to claim 3, characterized in that: The three-dimensional geological model software includes one or more of PETREL, GOCAD, 3DMine and PLAXIS 3D.

5. The method for constructing a three-dimensional model of fracture density of a fracture body according to claim 1, characterized in that: The calculating a second distance from a grid in the geological model to a main axis plane of the fault comprises: Each fault principal axis plane is taken as the starting plane, and a grid of a preset size is taken as the basic unit; Based on the starting surface and the basic unit, the distance from the grid in the model to the main axis of the fault is calculated.

6. The method for constructing a three-dimensional model of fracture density of a fracture body according to claim 1, characterized in that: The expression of the fitting function is: In the formula, Represents the density of cracks, in units of m -1 ; Indicates the distance from the fault, in m; A Indicates the vibration amplitude of crack development; k Indicates the rate of change of fracture development attenuation degree; ω Indicates the vibration frequency of crack development; b Indicates the baseline value of crack development.

7. The method for constructing a three-dimensional model of fracture density of a fracture body according to claim 1, characterized in that: The step of obtaining a three-dimensional model of crack density of a fracture body in a target area based on the fitting function and the second distance includes: Substituting the input second distance into the fitting function to obtain the crack density; Based on the crack density, a three-dimensional model of the crack density of the fault body in the target area is obtained.

8. A device for constructing a three-dimensional model of fracture density of a fracture body, characterized in that: include: An image acquisition module is used to acquire field outcrop images of target areas; A fracture line density acquisition module is used to mark fractures on the field outcrop image of the target area and obtain the density of the fractures; A scatter plot acquisition module, configured to obtain a scatter plot based on the density of the cracks and a first distance between the cracks and the fault; A fitting function acquisition module, used for fitting the data on the scatter plot to obtain a fitting function; A geological model acquisition module is used to construct a geological model of the target area; A second distance acquisition module, used for calculating a second distance from a grid in the geological model to a main axis plane of a fault; The crack density three-dimensional model acquisition module is used to obtain the crack density three-dimensional model of the fault body in the target area based on the fitting function and the second distance.

9. An electronic device, characterized in that: comprising a memory and a processor, wherein: The memory is used to store programs; The processor is coupled to the memory and is used to execute the program stored in the memory to implement the steps in the method for constructing a three-dimensional model of fracture density of a fracture body as described in any one of claims 1 to 7 above.

10. A computer-readable storage medium, characterized in that: Used to store computer-readable programs or instructions, which, when executed by a processor, can implement the steps in the method for constructing a three-dimensional model of fracture density of a fracture body as described in any one of claims 1 to 7.

Images