Patent Application
20160161933
This disclosure relates generally to the field of automated stage design and in particular to methods and systems for automated fracture stage design for hydraulic fracturing.
A common technique utilized to maximize production from permeability challenged reservoirs is fracturing. Fracturing involves creating fractures in the formation along which gas, oil, brine, or other substances may migrate to the well. These fractures are generally created to increase reservoir contact area and to maximize production. The most common form of fracturing is referred to as hydraulic fracturing which often involves injecting a high pressure liquid mixed with sand and chemicals into the formation.
A wellbore that is being fractured is often divided into different stages, and all stages may not respond the same way during a hydraulic fracture operation. For example, certain stages along a wellbore have increased fracture efficiency and thus respond better to a fracturing operation than others. A stage with high efficiency can promote a complex and highly connected fracture network which can remain open during the production and have acceptable volume of hydrocarbons in place. One way to optimize a hydraulic fracturing operation is to target the stages that have high fracture efficiency and thus increase production while minimizing the costs of operations. Thus the process of selecting and designing different stages is an important part of increasing overall efficiency of the operation.
Current industry practice for stage selection and design generally involves using single or multiple logs to perform a screening analysis along a wellbore. This analysis is generally quick and results in calculating an indicator which is then used as a proxy to indicate areas with high fracture potential. This process requires a user to input a certain set of parameters for stage selection and design. Once the indicator is calculated, location, number, and size of the different stages are usually manually selected. This is an oversimplified analysis of fracture efficiency, time consuming, non-repeatable, and in no way efficient.
Thus, currently used methods of performing fracturing stage design are inefficient and have certain shortcomings. The following disclosure addresses these and other issues.
In one embodiment a non-transitory program storage device, readable by a processor is provided. The non-transitory program storage device includes instructions stored thereon to cause one or more processors to receive a plurality of input parameters, each input parameter relating to a wellbore in a reservoir formation, to calculate one or more efficiency indices for the wellbore based on the input parameters, wherein each of the one or more efficiency indices may vary along a length of the wellbore, and to select one or more stage locations for the wellbore, each stage having a particular length, where the selection is done by utilizing an optimization process used to select stage locations based at least in part on at least one of the one or more efficiency indices having a minimal total intra-stage variation along the length of the wellbore.
In another embodiment, a method for designing fracture network staging for a fracturing operation is provided. The method includes receiving a plurality of input parameters that each relates to a wellbore in a reservoir formation, calculating one or more efficiency indices for the wellbore based on the input parameters, where each of the one or more efficiency indices may vary along a length of the wellbore, and selecting one or more stage locations for the wellbore, each stage having a particular length, where the selection is done by utilizing an optimization process used to select stage locations based at least in part on at least one of the one or more efficiency indices having a minimal total intra-stage variation along the length of the wellbore.
In yet another embodiment, a system is provided. The system includes in one embodiment, a memory, a display device, and a processor operatively coupled to the memory and the display device and adapted to execute program code stored in the memory. The program code is executed to receive a plurality of input parameters, each input parameter relating to a wellbore in a reservoir formation, to calculate one or more efficiency indices for the wellbore based on the input parameters, wherein each of the one or more efficiency indices may vary along a length of the wellbore, and to select one or more stage locations for the wellbore, each stage having a particular length, where the selection is done by utilizing an optimization process used to select stage locations based at least in part on at least one of the one or more efficiency indices having a minimal total intra-stage variation along the length of the wellbore.
In the following description, for purposes of explanation, numerous specific details are set forth in order to provide a thorough understanding of the inventive concept. As part of this description, some of this disclosure's drawings represent structures and devices in block diagram form in order to avoid obscuring the invention. Reference in this disclosure to “one embodiment” or to “an embodiment” means that a particular feature, structure, or characteristic described in connection with the embodiment is included in at least one embodiment of the invention, and multiple references to “one embodiment” or “an embodiment” should not be understood as necessarily all referring to the same embodiment.
It will be appreciated that in the development of any actual implementation (as in any development project), numerous decisions must be made to achieve the developers' specific goals (e.g., compliance with system- and business-related constraints), and that these goals will vary from one implementation to another. It will also be appreciated that such development efforts might be complex and time-consuming, but would nevertheless be a routine undertaking for those of ordinary skill in the art of data processing having the benefit of this disclosure.
In permeability challenged reservoirs, one of the primary goals of a fracturing operation is to create an extensive and well-connected fracture network that can remain open during production. However, not all locations around a well respond in the same manner to fracturing. Some areas respond more effectively than others. To maximize production from these areas, the wellbore is generally divided into different stages so that the areas within each stage have approximately a homogeneous response. The different stages are then treated differently during the fracturing operation depending on their fracture efficiency.
In addition to geomechanical properties affecting fracture geometry along the length of a wellbore, reservoir properties may also affect production from a fracture network. For example, properties such as water saturation or porosity can affect the overall fracture efficiency. If water saturation is high or porosity is low in a reservoir, for example, the efficiency of the reservoir will be low which affects the overall efficiency of the fracturing operation. Thus, similar to geomechanical properties, reservoir properties may need to be taken into account to ensure they have uniform properties along the length of each stage.
A stage design operation would need to predict fracture efficiency along the length of the wellbore based on known attributes and to use the predicted efficiency to divide the wellbore into various stages. In the past this was done manually and through a non-integrated quick look analysis which could be time consuming and give inaccurate and non-repeatable results. One embodiment of the present invention provides an integrated and automated procedure for selecting and designing stages that is quick, efficient, highly accurate, and repeatable. The automated procedure employs an optimization approach, which mathematically formalizes the notion of stage design accuracy.
In one embodiment, the fracture analysis and design of fracture networks operation involves presenting a user with a user interface screen such as the user interface screen 200 of
In addition to reservoir attributes, user interface screen 200 includes a section 204 for inputting geomechanical attributes and a section 206 for inputting user-defined attributes. Geomechanical attributes may include one or more of facture density, closure stress gradient, formation anisotropy, UCS, Friction angle, Poisson's ratio, Young's Modulus, critically stressed fracture potential, fracture closure index, stress anisotropy, and fracture height index. Any other geomechanical attribute may also be included. As with the reservoir attributes 202, a user may select which attribute to include in the analysis and may input values for poor 208, better 210, best 212, and weight 214 for each attribute. Alternatively, the user may choose to use predetermined values. The user-defined attributes section 206 may include any other attributes the user chooses to include in the analyses. For example, user-defined attributes may include mineral content, gas content, fracture cement fill type, and the like. In one embodiment, the user can define and add as many user-defined attributes to the user-defined attributes section 206 as the user wishes. Again, once these attributes are chosen, the user may input values for poor 208, better 210, best 212, and weight 214 for each attribute or may choose to use predetermined values. It should be appreciated that the integrated stage design procedure is not bound to any particular weighting scheme. Rather different embodiments may choose to implement different weighting options e.g., linear or any non-linear form. In one embodiment, the user also has the option of uploading data from a file for the chosen attributes. Various methods known in the art can be used for uploading the data. The uploaded values may be displayed in table 220. As shown, table 220 displays values of each of the attributes at each measured depth. Based on these values, the integrated procedure for stage design may determine and predict how each area of the wellbore will respond to a fracturing operation.
In one embodiment, once the attributes and their weights are selected and data is uploaded, the user may select to run an analysis on the wellbore. The results of such an analysis may be displayed to the user in a variety of different forms. For example, they may be presented on a graph, in a text file or in a table format. In one embodiment, the results are shown in graphs such as the ones shown in screen 300 of
In one embodiment, screen 300 includes a graph 302 for displaying an absolute fracture efficiency index at each measured depth. The absolute fracture efficiency index is a value calculated based on cut-off values and weights to evaluate the fracturing efficiency of a wellbore as a whole. For example, if all the attribute values at a specific depth fall into the best region based on the cut-off values, this depth will be determined as having acceptable fracturing efficiency. The absolute fracture efficiency index could help a user determine whether or not a wellbore is suitable for fracturing.
The screen 300 provides a variety of optional features for a user to select how to view the graphs and what to view on each graph. For example, it provides an option for the user to select a smoothing degree 306 for the graph 302. The user may select a medium smoothing degree and press on the update button 308 to view the graph 302 in the form illustrated.
Screen 300 also includes a graph 310 for illustrating relative fracture efficiency index curves. The relative fracture efficiency index is calculated by normalizing the fracture efficiency index through the length of the wellbore based on a weighting scheme that takes into account the reservoir, geomechanical and user-defined attributes and their weights and calculate an index that predicts the efficiency of fracturing at each depth in the wellbore. The weighting scheme could be either linear or non-linear. In particular, the weighting scheme can allow the selection of which types of attributes to take into account. For example, the relative fracture efficiency index may be calculated based on only one of the three types of reservoir, geomechanical, or user-defined attributes. Graph 310 can display the relative fracture efficiency curve calculated based on each of these types of attributes to show the contribution of each part on the overall efficiency. In one embodiment, the user has the option of selecting which type of relative fracture efficiency curve (i.e., calculate based on reservoir, geomechanical, or user-defined attributes) to plot by using the dropdown box 316 and pressing on the update button 318. The user may also select a smoothing degree 314 for the graph 310.
In addition to the absolute fracture efficiency index and the relative reservoir or geomechanical index, screen 300 also provides a graph 320 for viewing the overall relative fracture efficiency index at each depth. As discussed above, the overall relative fracture efficiency index takes into account the weight of each selected attribute to calculate an index that predicts the relative efficiency of fracturing at each depth in the wellbore. The weighted combination of selected attributes is normalized along the well length to form the relative fracture efficiency index. This index may be used to show what places in the wellbore provide a better fracturing efficiency relative to others for each wellbore, even if the wellbore as a whole does not have a high absolute fracture efficiency index. In other words, the values of cut-offs are not used for calculating the relative efficiency index. Again, as with plot 302 and 310, a smoothing degree may be selected through the dropdown box 324 and by pressing the update button 326.
After the various fracture efficiency indices are calculated, the information can be used to automatically design stages for the wellbore. In one embodiment, this process may be started by a user using the graphical interface of
Using an automated stage design algorithm, the well can be divided into different stages with least total intra-stage variation of relative fracture efficiency index as shown in plot 320 by lines 312. These lines divide the wellbore along its length into different stages. The stages are automatically selected based on a combination of values for minimum and maximum stage length, number of stages and the relative fracture efficiency index for the wellbore at each depth. The automated stage design algorithm involves an optimization process used to select stage locations yielding a minimum total fracture efficiency index variation within each stage. The optimization goal is to keep the efficiency indices as uniform as possible within each stage. Integrating all of the staging parameters and input curves, the algorithm automatically chooses the dividing lines for the stages, thus reproducibly, accurately, and efficiently designing the stages. The plot not only shows how to stage the wellbore, it can be combined with absolute fracture efficiency index to illustrate the efficiency of the fracturing operation at each stage, thus allowing the user to determine how to treat each stage. Although shown here on a graph, it should be understood that the resulting stage design can be outputted and presented to the user in a variety of different manners known in the art.
The automatic stage design algorithm and procedure may provide an option to skip stages. This is illustrated in
The algorithm may also provide the option of selecting whether to take faults into account. In general, it is preferred not to fracture areas that are near a fault or other similar geological features. By selecting the Faults box 340, as shown in
Once the geomechanical, reservoir attributes, or user-defined attributes are determined (or received), weights are assigned to each of those attributes (block 410). In one embodiment, this is done by operation 400 based on predetermined criteria or specific hydrocarbon plays. For example, if a default setting is selected, the operation 400 may assign a predetermined weight for each given attribute. Alternatively, a user may be able to enter desired weights for each attribute that will be used in the calculations. In addition to weights, cutoffs are also assigned to each attributes (block 412). In one embodiment, cutoffs are values that determine the quality of the attribute. For example, water saturation above a certain threshold value may be designated as a poor value, while water saturation below another certain threshold value may be considered a best value. Values in between the two threshold values may be considered average values. In one embodiment, default cutoff values are assigned for each hydrocarbon play. Alternatively, a user may be able to enter the cutoff values for one or more of the attributes.
After the attributes are calculated and weights and cutoffs are assigned, the operation moves to block 416 of
Once the fracture efficiency indices for the wellbore have been calculated and it is determined that the wellbore should be fractured, an automated stage design algorithm may begin running (block 420). The automated stage design algorithm may use the values from a chosen fracture efficiency index curve (e.g., absolute or relative) to automatically select and design stage(s) for the wellbore yielding nearly uniform fracture characteristics along each stage. In one embodiment, the automated algorithm is designed so as to run on multiple index curves if the user selects this option. Instances wherein multiple index curves may be input to the algorithm include cases when each raw attribute curve is viewed as an individual indicator of fracture efficiency and when the user does not want any loss of information as a result of the aggregation of the raw curves into a single index curve. More generally, multiple index curves may be fed to the automation algorithm when various fracture efficiency indices are generated each using a different weighting/aggregation scheme. This is because each index may be viewed as an aggregate feature and by having a collection of aggregate features, loss of information resulting from each single aggregation may be reduced or avoided altogether. Running the stage automation algorithm may be done completely automatically or the user may have the option of selecting certain parameters in the design. Thus, at this point, the operation determines if there are any user inputs to the stage design algorithm (block 422). If it is determined that no user inputs should be used in stage design, then the operation is indicted as fully automated (block 424). In such a case, the automated operation calculates and outputs the stage lengths, number, and locations of the designed stages (block 438). The results may be shown to the user in a graphical form (e.g. on a graph), tabular manner, or any other method known in the art, before the operation ends (block 440).
When at block 422 it is determined that the operation is not fully automated and the user should input parameters, then operation 400 moves to determine if skipping stages is allowed (block 426). This may be decided by a user or may be predetermined by the operation. When skipping stages are allowed, the operation receives a threshold value for skipping (block 428) stages. This value, may, for example be a value for the average relative fracture efficiency index for a stage below which the stage should be skipped. After receiving the threshold value or when it is determined that skipping stages is not allowed, the operation moves to determine if a fault is detected (block 430). In one embodiment, this step may be optional or only performed when selected by a user. When a fault is detected, the operation moves to receive the fault's coordinates and characteristics (block 432). After receiving the fault coordinates or when no fault is detected, the operation moves to receive stage design constraint (block 434). This may include minimum and maximum length for the stages and the number of stages a user desires to have. Other stage related information may also be provided and taken into account. Once the data has been received, the operation calculates and outputs the stage lengths, number and location of the designed stage network (block 442) for the wellbore. The results may be shown to the user (block 436) in a graphical form (e.g. on a graph), may be presented in tabular manner, or any other method known in the art, before the operation ends (block 440). Thus, operation 400 as detailed in
The stage design algorithm provides full autonomy in the process of dividing a reservoir zone into fracture stages by employing an optimization-driven staging approach for real-time operational decision making, therefore assuring consistency and repeatability of the interpretation in addition to relieving the human analyst from a laborious, error-prone, and often-subjective interpretation.
In general, the stage design automation algorithm casts the stage design problem as that of discrete optimization and solves it. The essence of the optimization relies on the fact that any partitioning of the wellbore area into fracture stages given fracture efficiency indices can be objectively assessed in terms of a measure that quantifies how well the efficiency index curves are segmented into stages. The algorithm then seeks to optimize such measure. An objective measure for assessing the quality of any partitioning should reward partitions wherein individual stages are as homogeneous in terms of their efficiency indices as possible. It follows that by defining the homogeneity of the efficiency indices for any given stage, the total stage homogeneity may serve as an objective measure of the quality of any stage design (i.e., partitioning). If the collection of efficiency indices at any fixed depth are referred to as a sample then the efficiency index homogeneity for any stage may be defined as the average distance from all in-stage samples to the mean in-stage sample, assuming any particular distance function (e.g., Euclidean). Hence, the lower the average distance to the mean is, the more homogenous the stage is. Thus the algorithm should seek to minimize the sum of all intra-stage average distances.
In one embodiment, to minimize the objective measure of stage partitioning, the automation algorithm makes use of a local optimization module illustrated by operation 500 of
The locally driven objective function minimizer follows an iterative strategy. Given the initial stage design, it is possible to assess the effect of moving one stage boundary to the left or the right by a fixed step size on the value of the objective function. The differences in the objective value from all individual fixed-step left and right boundary moves may be computed. The combination of boundary moves that yield the largest sum of all individual decrements to the original objective value may be chosen. Once all chosen boundary moves are applied, the new objective value may be computed and the process may be iterated until convergence (i.e., when no boundary moves are possible). The process of selecting the combination of the boundary moves at each iteration may be cast as an integer programming problem. Theoretically, such an iterative process is convergent given that the objective value is lowered at every iteration and that there is a finite set of partitions (solutions) to search-through. It should be noted that this optimization process is subjective to the choice of the prefixed step size. It is therefore beneficial if the entire process is repeated for different step sizes. In one embodiment, an ordered collection of decreasing step sizes is predetermined and the optimization process is iteratively sequenced such that the terminal solution using one fixed step size feeds into the initial solution of the optimization process with the next lower step size and so on.
The local optimization module of the automated algorithm for stage partitioning is outlined step-by-step in operation 500 of
It should be noted that assuming every sample within the wellbore area is sampled at a constant depth interval, then every stage partition (i.e., solution) of the set containing all the samples within the wellbore area mathematically corresponds to a unique integer composition of the number of samples where the number of terms in the composition is equal to the number of stages and each term in the composition is within any imposed stage length constraints. The solution space is therefore on the order of all permissible integer compositions of the sample size having as many terms as the allowed number of stages. This optimization process may be viewed as a projected steepest-descent algorithm wherein the projection space assures that every solution is an integer composition satisfying the problem constraints and that the descent algorithm walks through a path of integer compositions in a decreasing order of objective value starting from the initial integer composition. This particular view of the solution space is helpful for constructing the initial solution sampling module. The step size in the projected steepest-descent algorithm trades off the accuracy resolution of the solution for the required number of iterations for convergence. By appropriately choosing the step size collection, it is possible to achieve both good time overhead and solution accuracy.
One embodiment of the initial solution sampling algorithm is outlined in operation 600 of
Sampling a constrained integer partition of the number of samples can be done by imposing a lexicographical (i.e., dictionary) ordering on the constrained integer partitions. The lexicographically-ordered constrained integer partitions may then be visualized via a tree structure wherein every node in the tree represents an integer term in a partition and every path from the tree root to a terminal (i.e., leaf) node represents one constrained integer partition. The sampling algorithm may then iteratively sample one term at a time from every level in the tree in a top-down fashion. No data structure is required for the sampling as it is possible to compute the range of all integer terms at any given tree level knowing the path from the root containing all sampled parent terms. Sampling from the tree of all constrained integer partitions may be performed according to any chosen discrete probability. Once a constrained integer partition is sampled randomly, a random permutation may be chosen so as to permute the generated constrained integer partition into a constrained integer composition.
The constrained integer composition sampling algorithm is equivalent to a random top-down traversal of the tree associated with the lexicographical ordering of the space of all respective constrained integer partitions. The random tree exploration is controlled by a probability distribution that is defined via the choice of the sequence of random sampling distributions at all iterations followed by the random permutation sampling distribution.
Accordingly, operation 600 begins (block 602) by receiving data (block 604) such as the efficiency index curves and stage constraint parameters. After data is received, the operation chooses a sample size M (block 606) and then then iteratively samples M random constrained integer compositions according to a chosen probability density function (block 608).
Once M random solutions (i.e., constrained integer compositions) are sampled, the operation computes M local optima associated with the M random solutions by running operation 500 M times (block 610). Then the best solution amongst the M computed local optima is identified (block 612) and returned at block 614.
The operations of
In deciding which stages to flag for skipping, a user-defined minimum threshold on the average stage index (or indices for multiple input curves) may be used in one embodiment. In such a case all stages having an average index (indices) lower than the minimum threshold are flagged for skipping. The threshold may alternatively be defined with any other static of interest such the maximum allowed variability within any stage. An alternative embodiment for flagging the stages to be skipped may utilize an outlier detection algorithm that automatically flags the stages that need to be skipped. In such an embodiment, an outlier may be defined as any stage with an efficiency index distribution of values that is lower by a specific threshold than those of the most other stages. An additional alternative embodiment for stage skipping flagging may be to employ a training framework wherein a training dataset of multiple depth staging configurations are collected each along with its target stages to be flagged. The training framework may learn a classification function able to flag stages to be skipped in any never before seen staging configuration.
In addition to the need to consider skipping some stages, it may also be necessary to take into account areas close to fault(s) around which fracturing should be avoided. As discussed before it is generally preferable to skip fracturing depth ranges that are affected by faults or similar geological areas. Accordingly, to take into account faults, the operations of
Generally faults, defined by their center locations and proximity ranges, divide the complete depth zone up into disjointed depth subzones. In other words, the continuum of the complete depth zone is disrupted by several episodes of fault-inflicted depth subzones. A sub-problem over a non-fault-inflicted subzone may be solved by the operations of
Thus, the forgoing solutions provide embodiments for partitioning a wellbore depth zone into multiple stages automatically, accurately, and efficiently while providing the ability to take into account skipping and faults.
In the foregoing description, for purposes of explanation, specific details are set forth in order to provide a thorough understanding of the disclosed embodiments. It will be apparent, however, to one skilled in the art that the disclosed embodiments may be practiced without these specific details. In other instances, structure and devices are shown in block diagram form in order to avoid obscuring the disclosed embodiments. References to numbers without subscripts or suffixes are understood to reference all instance of subscripts and suffixes corresponding to the referenced number. Moreover, the language used in this disclosure has been principally selected for readability and instructional purposes, and may not have been selected to delineate or circumscribe the inventive subject matter, resort to the claims being necessary to determine such inventive subject matter. Reference in the specification to “one embodiment” or to “an embodiment” means that a particular feature, structure, or characteristic described in connection with the embodiments is included in at least one disclosed embodiment, and multiple references to “one embodiment” or “an embodiment” should not be understood as necessarily all referring to the same embodiment.
It is also to be understood that the above description is intended to be illustrative, and not restrictive. For example, above-described embodiments may be used in combination with each other and illustrative process acts may be performed in an order different than discussed. Many other embodiments will be apparent to those of skill in the art upon reviewing the above description. The scope of the invention therefore should be determined with reference to the appended claims, along with the full scope of equivalents to which such claims are entitled. In the appended claims, terms “including” and “in which” are used as plain-English equivalents of the respective terms “comprising” and “wherein.”
