Patent Application
20230095763
Fractures typically present naturally in rocks and often exhibit high conductive pathways for fluids. Hydraulic properties of the natural fractures are essential parameters for modeling fluid flow and transport in fractured formations. Current technologies mostly focus on capturing static geometric properties of the natural fractures, such as strike and dip angle. However, dynamic hydraulic properties, such as fracture permeability and hydraulic aperture, cannot be calculated directly because they depend on fracture conductivity which requires special and costly measurements. In addition to these, due to heterogeneous nature of the natural fractures, for example, non-uniform apertures, roughness, and tortuosity, it is often challenging to estimate the dynamic hydraulic properties of the natural fractures. However, estimating the dynamic hydraulic properties of the natural fractures is vital to assess fluid flow behavior in subsurface rock formation. When modeling hydrocarbon reservoirs, fractures conductivity, which is captured by the fracture permeability and effective hydraulic aperture, is important to evaluate recovery performance of different recovery schemes, including waterflood and enhanced oil recovery.
This summary is provided to introduce a selection of concepts that are further described below in the detailed description. This summary is not intended to identify key or essential features of the claimed subject matter, nor is it intended to be used as an aid in limiting the scope of the claimed subject matter.
In one aspect, embodiments disclosed herein relate to a method for fracture dynamic hydraulic properties estimation and reservoir simulation that includes obtaining, by a computer processor, a first set of high-resolution images of a first fracture. The method includes obtaining, by the computer processor and a first model, a first set of fracture detections based on the first set of high-resolution images. The method includes generating, by the computer processor, a plurality of numerical calculations based on the first set of fracture detections of the first fracture. The method includes generating, by the computer processor, a second model based on the plurality of numerical calculations and the first set of fracture detections. The method includes obtaining, by the computer processor, a second set of high-resolution images of a second fracture of a new reservoir. The method includes generating, by the computer processor using the first model, a second set of fracture detections of the second fracture. The method includes generating, by the computer processor using the second model, dynamic hydraulic estimations of the second fracture. The method includes generating, by the computer processor and a third model, a three-dimensional (3D) reservoir simulation of the new reservoir based on the second set of high-resolution images and the dynamic hydraulic estimations of the second fracture. The method further includes determining, by the computer processor and using the dynamic hydraulic estimations of the second fracture and the 3D reservoir simulation, a plurality of recovery schemes for the new reservoir.
In one aspect, embodiments relate to a system for fracture dynamic hydraulic properties estimation and reservoir simulation. The system includes a plurality sets of high-resolution images for a plurality fractures. The system includes a fracture manager comprising a computer processor. The fracture manager obtains a first set of high-resolution images of a first fracture. The fracture manager obtains, using a first model, a first set of fracture detections based on the first set of high-resolution images. The fracture manager generates a plurality of numerical calculations based on the first set of fracture detections of the first fracture. The fracture manager generates a second model based on the plurality of numerical calculations and the first set of fracture detections. The fracture manager obtains a second set of high-resolution images of a second fracture of a new reservoir. The fracture manager generates, using the first model, a second set of fracture detections of the second fracture. The fracture manager generates using the second model, dynamic hydraulic estimations of the second fracture. The fracture manager generates, using a third model, a three-dimensional (3D) reservoir simulation of the new reservoir based on the second set of high-resolution images and the dynamic hydraulic estimations of the second fracture. The fracture manager further determines, using the dynamic hydraulic estimations of the second fracture and the 3D reservoir simulation, a plurality of recovery schemes for the new reservoir.
In one aspect, embodiments relate to A non-transitory computer readable medium storing instructions. The instructions obtain a first set of high-resolution images of a first fracture. The instructions obtain a first set of fracture detections based on the first set of high-resolution images. The instructions generate a plurality of numerical calculations based on the first set of fracture detections of the first fracture. The instructions generate a second model based on the plurality of numerical calculations and the first set of fracture detections. The instructions obtain a second set of high-resolution images of a second fracture of a new reservoir. The instructions generate a second set of fracture detections of the second fracture. The instructions generate dynamic hydraulic estimations of the second fracture. The instructions generate a three-dimensional (3D) reservoir simulation of the new reservoir based on the second set of high-resolution images and the dynamic hydraulic estimations of the second fracture. The instructions further determine a plurality of recovery schemes for the new reservoir.
Other aspects and advantages of the claimed subject matter will be apparent from the following description and the appended claims.
Specific embodiments of the disclosed technology will now be described in detail with reference to the accompanying figures. Like elements in the various figures are denoted by like reference numerals for consistency.
In the following detailed description of embodiments of the disclosure, numerous specific details are set forth in order to provide a more thorough understanding of the disclosure. However, it will be apparent to one of ordinary skill in the art that the disclosure may be practiced without these specific details. In other instances, well-known features have not been described in detail to avoid unnecessarily complicating the description.
Throughout the application, ordinal numbers (e.g., first, second, third, etc.) may be used as an adjective for an element (i.e., any noun in the application). The use of ordinal numbers is not to imply or create any particular ordering of the elements nor to limit any element to being only a single element unless expressly disclosed, such as using the terms “before”, “after”, “single”, and other such terminology. Rather, the use of ordinal numbers is to distinguish between the elements. By way of an example, a first element is distinct from a second element, and the first element may encompass more than one element and succeed (or precede) the second element in an ordering of elements.
As hydrocarbon discoveries are declining significantly, there is a compelling need in the gas and oil industry to enhance oil recovery from maturing fields and to develop challenging resources, such as low-quality rock reservoirs, naturally-fractured reservoirs, and shales. These reservoirs are multiple continuum where fluid flow is governed by a complex interplay between fractures and matrix systems. Modeling such fractures is crucial to understand matrix-fracture interactions and recovery mechanisms.
A fracture, also referred as a natural fracture, is a crack or surface of breakage within rock. In conventional hydrocarbon reservoirs, understanding the fractures is crucial to recovery the oil and gas behavior. For example, fractures can enhance rock permeability and water injectivity. Thus, fractures are induced mechanically in reservoirs to boost hydrocarbon flow. On the other hand, fractures may lead to early water breakthrough and poor sweep. In unconventional reservoirs, fractures are the primary conduits for flow and production.
High-resolution wellbore images are commonly used to characterize fractured reservoirs. In particular, high-resolution borehole images can reduce uncertainties when realizing heterogeneous rock matrix and embedded discrete fracture network (DFN). In particular, a DFN refers to a computational model that represents the geometrical properties of each individual fracture, and the topological relationships between individual fractures. Further, reservoir modeling requires both static geometric properties (such as fracture density, orientation, dip, and strike), and dynamic hydraulic properties (such as fracture permeability and effective aperture). However, existing techniques focus on detecting the static geometric properties. None of the existing techniques discusses estimating the dynamic hydraulic properties by using high-resolution wellbore images. In addition, traditional techniques and methods require physical measurements and simulations to determine the fracture properties. Due to complexity and expensive cost, these measurements and simulations can only be done on a limited scale.
This disclosure provides an automatic procedure that estimates the dynamic hydraulic properties of a fracture by using high-resolution wellbore images. The disclosure provides an efficient and cost-effective procedure that incorporate big data of wells, cores, and outcrops based on machine-learning (ML) models. Embodiments disclosed herein also provide workflows that aim at capturing fracture data from different sources and integrate a reservoir simulation model accordingly. In particular, the provided workflows can be integrated with existing workflows.
In general, embodiments of the disclosure include a system and a method that automatically process fracture detections and dynamic hydraulic properties estimation from high-resolution image. The present disclosure relates to a method of collecting high-resolution images, detecting fractures from the collected images, and estimating dynamic hydraulic properties from the detected fractures. More specifically, a new workflow based on machine learning techniques aims to automate the process of fracture detection from high-resolution wellbore images, and other images of rock cores and outcrops, and estimation of the hydraulic properties of natural fractures. In some embodiments, the method may utilize a first ML model to generate a plurality of detected fractures. Further, the method may utilize the plurality of detected fractures to generate a second ML model that estimates dynamic hydraulic properties. Additionally, the method may include generating dynamic hydraulic properties of a fracture utilizing the second ML model. Moreover, the method may include updating the second ML model based on newly collected high-resolution images and fracture data. Furthermore, the method may utilize the collected high-resolution images and the estimated dynamic hydraulic properties to construct reservoir simulation model.
Turning to
Specifically, the wellbore images (251) may refer to high-resolution images of wellbore interceptions. For example, the wellbore images (251) may reflect wellbore interceptions of the one or more wells (111-115) in
Keeping with
In some embodiments, the fracture manager (200) may include a graphical user interface (GUI) (e.g., graphical user interface (GUI) (205)) that receives instructions and/or inputs from users. More specifically, the users may enter various types of instructions and/or inputs via the GUI (205) to start certain actions, such as calculating, evaluating, selecting, and/or updating data or parameters.
In some embodiments, the fracture manager (200) may include a data controller (e.g., data controller (210)). The data controller (210) may be software and/or hardware implemented on any suitable computing device, and may include functionalities for collecting various data from the high-resolution images source (250) and processing the collected data. For example, the data controller (210) may collect the wellbore images (251), the rock core images (252), and the outcrop images (253) from the high-resolution image resource (250). The data controller may include data processors (e.g., data processors (215)) and data storage (e.g., data storage (216)). Specifically, the data processors (215) process the data collected from the high-resolution images source (250) as well as the data stored in the data storage (216). The data storage (216) may store the various data collected from the high-resolution images source (250), and other data and parameters for the other sections and functionalities of the fracture manager (200).
Keeping with
Further, the fracture manager (200) may include a hydraulic property estimation model (e.g., hydraulic property estimation model (230)) that utilizes one or more deep-learning algorithms (e.g., deep-learning algorithm (235)) to estimate dynamic hydraulic properties of a fracture (e.g., dynamic hydraulic property estimations (270)). In particular, the hydraulic property estimation model estimates permeabilities (e.g., fracture permeabilities (271)) and hydraulic apertures (e.g., hydraulic apertures (272)) of the detected fractures. In particular, the fracture permeabilities may be associated with the secondary porosity created by an open fractures. In many reservoirs, the fracture permeabilities can be a major controlling factor of flow of fluids. Further, the fracture permeabilities are critical parameters in reservoir simulations and ultimately impact future determinations based on the simulation results.
Further, the hydraulic apertures may refer to effective perpendicular width of the open fracture. More details about the hydraulic property estimation model (230) and the deep-learning algorithm (235) are provided below in
Continuing with
The 3D reservoir simulations are highly related to reservoir development. Specifically, generated to evaluate the amount of hydrocarbon in a reservoir, such as original oil in place and recoverable reserves, which can be used to evaluate profitability of the reservoir and to design future development plans of the reservoir. For example, reservoir simulation models can be used to predict field performance and ultimate recovery for various field development scenarios, such as well orientation, hydraulic fracturing, enhanced oil recovery, and etc., to evaluate effects of different operational conditions recovery and to compare economics of different recovery methods. More specifically, the 3D reservoir simulations are used to determine, for example, how to develop a field to maximize economic recovery, the best enhanced recovery scheme for the reservoir, the best completion scheme for wells, and reservoir production source.
In some embodiments, the high-resolution image source (250) and the fracture manger (200) may be implemented on the same computing device, or different computing systems connected by a network. In some embodiments, the high-resolution image source (250), the fracture manger (200) and/or other elements, including but not limited to network elements, user equipment, user devices, servers, and/or network storage devices may be implemented on computing systems similar to the computing system (600) shown and described in
While
Turning to
In Step 1, high-resolution images are obtained from wellbores interceptions (300). Wellbore interceptions illustrate where a wellbore that is drilled intersects a fracture in the formation, as shown in 300. The high-resolution images include at least wellbore images, rock core images, and outcrop images. For example, the high-resolution images may refer to the wellbore images (251), rock core images (252) and the outcrop images (253) obtained by the data controller (220) as illustrated in
In Step 2, fractures are detected from the wellbore images utilizing a fracture detection model. The fractures may be detected from the fracture detection model (220) and may be represented by the fracture detections (260) as described in
In Step 3, dynamic hydraulic properties are generated based on the detected fractures from Step 2 and a hydraulic property estimation model. For example, the dynamic hydraulic properties may be generated by the hydraulic property estimation model (230) and may be represented by the dynamic hydraulic property estimations (270) as illustrated in
In Step 4, a 3D reservoir simulation (308) is generated. Specifically, the 3D reservoir simulation is generated based on the high-resolution images from Step 1, the fracture detections from Step 2, and the dynamic hydraulic property estimations from Step 3. The 3D reservoir simulation may be generated by the reservoir simulation model (240) and may be represented by the 3D reservoir simulations (280) as explained in
Those skilled in the art will appreciate that the process of
Turning to
In
ML models include supervised ML models and unsupervised ML models. More specifically, supervised ML models include classification, regression models, etc. Unsupervised ML models include, for example, clustering models. Deep-learning algorithms are a part of ML algorithms based on artificial neural networks with representation learning. For example, the deep-learning algorithm may run data through multiple layers of neural network algorithms, each of which passes a simplified representation of the data to the next layer. More specifically, each artificial neural network consists a plurality of neurons that are staked next to each other and organized in layers. Each neuron may receive various inputs, multiplies the inputs by weights, sums them up, and applies a non-linear function.
In some embodiments, the fracture detection model (420) uses ML algorithms, such as the deep-learning algorithm (425), to identify and recognize fractures from the wellbore images (410) using an image segmentation technique, and outputs the identified and recognized fractures as the fracture detections (460). More specifically, the deep-learning algorithm (425) may utilize training based on pixel labeling of grey-scale images. Further, the deep-learning algorithm (425) may also utilize other image segmentation algorithms, such as Random Forest, which is based on probabilities of pixels.
In some embodiments, the deep-learning algorithm (425) may utilize U-Net procedure to generate the fracture detections (460). Specifically, the U-Net algorithm comprises two paths, a contracting path (also known as a convolutional path), and an expanding path (also known as a de-convolutional path). More specifically, the contracting path includes several cycles where each cycle has multiple repeating convolution operators (3×3), followed by a Rectified Linear Unit (ReLU), followed by a pool operator (2×2). The de-convolution path includes several cycles where each cycle has multiple repeating convolution operators (3×3), followed by a ReLU. Layer at the output end includes one convolution and one softmax.
However, the deep-learning algorithm (425) is not limited to this, and may be one or more neural network architectures, such as artificial neural networks (ANN), convolutional neural networks (CNN), recurrent neural networks, general adversarial neural networks, deep belief networks, autoencoders, etc.
Keeping with
In some embodiments, the trained hydraulic property estimation model (430) takes the fracture detections (460) from the fracture detection model (420) as input, and estimates hydraulic property estimations (e.g., dynamic hydraulic property estimations (470)) corresponding to the fracture detections (460) as output. In particular, the dynamic hydraulic property estimations (470) include the fracture permeabilities and the hydraulic apertures are outputted.
In some embodiments, the wellbore images (410) may be collected from a new fracture of an existing or a new reservoir. By utilizing the fracture detection model (420) and the trained hydraulic property estimation model (430), dynamic hydraulic property estimations (470) corresponding to the newly collected wellbore images are generated.
While
Turning to
In Block 510, a plurality of high-resolution wellbore images are obtained. For example, the obtained high-resolution wellbore images may be represented by the wellbore images (410) as shown in
In Block 520, a plurality of fractures are detected and identified based on the wellbore images and a fracture detection model. For example, the fracture detection model may be the fracture detection model (420) as illustrated in
In Block 530, a plurality of numerical calculations are generated as training data based on the detected fracture detections. For example, Block 530 may refer to the first part of generating the hydraulic property estimation model (430), and the generated calculations may be the training data (430) as illustrated in
More specifically, is velocity vector, ρ is density, p is pressure, and μ is viscosity.
ρ(·∇
)=−∇p+μ∇2
∇·
=0 (1)
Continuing with
Specifically, as shown in
More specifically, high-fidelity discrete wellbore images including the fractures are pre-processed to generate input data for the a CNN model. Similar to those described in
To accelerates the feature extraction and training-validating processes, a filling-in process is conducted in Steps 2-3 to mark the zone bounded by the fracture profile, followed by an image-coarsening process. In particular, during the image-coarsening process, a 50×1000 pixel-resolution is applied. This resolution is a compromise between accuracy and efficiency, which is sufficient for characterizing the roughed fractures without significant loss of accuracy.
Furthermore, the region of interest (rectangle zone in Steps 2-3 of
In particular,
In Block 550, the trained hydraulic property estimation model from Block 540 is verified. In some embodiments, the verification utilizes the numerical calculations from Block 530. In particular, the calculations used for this verification may or may not be same as the calculations used for training the hydraulic property estimation model. More specifically, in order to obtain an accurate, efficient, and stable hydraulic property estimation model, relative error D of the generated model may range from −22% to 22%, with an essential normal distribution crossing the line of D=0 and most points located within D=±10%. In particular, accuracy of the generated model highly depends on accuracy of the training data that is used to train the model. In some embodiments, high-quality data from simulation and lab measurements are acquired to validate the hydraulic property estimation model. For example, a trained and verified hydraulic property estimation model may be represented by the hydraulic property estimation model (230) as shown in
In Block 560, a determination as to whether the trained hydraulic property estimation model is accurate. If the trained model is determined as accurate, the flowchart moves to Block 570. If the answer to Block 560 is no, the flowchart loops back to Block 510, to update the trained model with newly collected wellbore images, and the process repeats until the trained model is determined as accurate.
In Block 570, a new plurality of high-resolution wellbore images of a new fracture are collected. These new images may be collected from a reservoir that is same as or different from the reservoir, from which the high-resolution wellbore images in Block 510 are collected.
In Block 580, dynamic hydraulic properties for the new fracture are generated using the fracture detection model and the trained hydraulic property estimation model.
In Block 590, a 3D reservoir simulation for the new fracture is generated using the new plurality of high-resolution wellbore images.
In Block 595, the hydraulic properties and the 3D reservoir simulation for the new fracture are used to determine performance of a plurality of recovery schemes designed for a reservoir that includes the new fracture. By obtaining permeability distribution of different fractures and rocks based on the 3D reservoir simulation, mathematical equations for fluid transport in each block of the 3D reservoir simulation can be solved. These mathematical equations are further used to predict produced hydrocarbon at the simulated reservoir and surface conditions, for example, reservoir and surface temperature, pressure, and water saturation.
Those skilled in the art will appreciate that the process of
Turning to
In one or more embodiments, for example, the input device 620 may be coupled to a receiver and a transmitter used for exchanging communication with one or more peripherals connected to the network system 630. The receiver may receive information relating to hydraulic properties as described in
Further, one or more elements of the computing system 600 may be located at a remote location and be connected to the other elements over the network system 630. The network system 630 may be a cloud-based interface performing processing at a remote location from the well site and connected to the other elements over a network. In this case, the computing system 600 may be connected through a remote connection established using a 5G connection, such as protocols established in Release 15 and subsequent releases of the 3GPP/New Radio (NR) standards.
The computing system in
While
Although only a few example embodiments have been described in detail above, those skilled in the art will readily appreciate that many modifications are possible in the example embodiments without materially departing from this invention. Accordingly, all such modifications are intended to be included within the scope of this disclosure as defined in the following claims. In the claims, means- plus-function clauses are intended to cover the structures described herein as performing the recited function and not only structural equivalents, but also equivalent structures. Thus, although a nail and a screw may not be structural equivalents in that a nail employs a cylindrical surface to secure wooden parts together, whereas a screw employs a helical surface, in the environment of fastening wooden parts, a nail and a screw may be equivalent structures. It is the express intention of the applicant not to invoke 35 U.S.C. § 112, paragraph 6 for any limitations of any of the claims herein, except for those in which the claim expressly uses the words ‘means for’ together with an associated function.
