This invention relates generally to the field of prospecting and more particularly to data processing. Specifically, the invention is a method for modeling subsurface volumes for one or more geophysical properties. The present invention relates generally to the field of hydrocarbon production, and more particularly to the effect of anisotropic formation strength on the drilling of a stable wellbore along a given locus. As a specific example, the present techniques involve a petrophysical method for utilizing geophysical wireline logs to predict the shear strength anisotropy of sub-surface fine-grained rock formations (particularly shales) without having to core and plug the sub-surface formation of interest to directly measure shear strength anisotropy in the laboratory. Predicted anisotropic strength functions are subsequently used in numerical models of the subsurface, including analytical and computational methods for optimizing mud weight design in high-angle/horizontal extended reach wells in which inclinations can approach bed-parallel.
In the oil and gas industry, modeling of the subsurface is typically utilized for visualization and to assist with analyzing the subsurface volume for potential locations for hydrocarbon resources. Accordingly, various methods exist for estimating the geophysical properties of the subsurface volume (e.g., information in the model domain) by analyzing the measurements from measurement equipment (e.g., information in the data domain). The measurements include some information of the geophysical properties that may be utilized to generate the subsurface model.
The interpretation of the subsurface volume may be complicated by the rock properties within the subsurface. In particular, interest in the anisotropic strength characteristics exhibited by fine-grained rock formations has resulted primarily from a desire to understand wellbore failure mechanisms associated with drilling through shales. Although elastic anisotropy in shales can impact wellbore failure through the redistribution of borehole stresses, analytical studies have shown that strength anisotropy has a more significant control. See e.g., Gazaniol, D., T. Forsans, M. J. F. Boisson and J.-M. Piau. 1995. Wellbore failure mechanisms in shales: Prediction and prevention. J. Pet. Tech. July 1995: 589-595.
This effect is particularly pronounced for extended reach drilling (ERD) where the borehole trajectory can achieve near-parallelism with weak bedding planes, as shown in
For optimally orientated positions around the borehole circumference, failure may occur via oblique loading of the planes of weakness (intra-laminar failure) before failure can occur within the bulk material cutting through the planes of weakness (trans-laminar failure), necessitating higher mud weights to reduce stress concentration on weak bedding planes, or else limiting wellbore inclination.
As the recovery of natural resources, such as hydrocarbons rely, in part, on a subsurface model, a need exists to enhance subsurface models of one or more geophysical properties. In particular, a technique for predicting the degree of shear strength anisotropy encountered in fine-grained rock formations (including shales) without recourse to costly and time consuming laboratory measurements is needed. The technique may preferably be suited for utilizing geophysical wireline logs as input data to the predictive technique for evaluating strength anisotropy. The technique calculates a lithology-dependent anisotropic failure criterion that can subsequently be input into analytical or computational methods for optimizing mud weight design in ERD wells. See, e.g., Ottesen, S., R. H. Zheng and R. C. McCann. 1999. Borehole stability assessment using quantitative risk analysis. In Proceedings of the IADC/SPE Drilling Conference, Amsterdam, 9-11 Mar. 1999 and Pande, G. N., G. Beer and J. R. Williams. 1990. Numerical Methods In Rock Mechanics. John Wiley & Sons, NY.
In one embodiment, a method for predicting strength anisotropy in subsurface formations along a wellpath comprising: receiving a stress model of the subsurface formation in the area of the wellpath; obtaining bulk strength properties of target rock; developing anisotropic failure criteria based on the bulk strength properties; combining the stress model with the anisotropic failure criteria to develop an anisotropic wellbore stability model; and using the anisotropic wellbore stability model to enhance hydrocarbon recovery.
In yet another embodiment, a method for predicting strength anisotropy in subsurface formations is described. The method comprises: compiling a database of core-based strength measurements over a range of orientations and confining pressures for a plurality of fine-grained rock formations; utilizing an automated surface fitting procedure to determine a common anisotropic failure criterion for all database rocks; developing predictive algorithms relating anisotropy parameters in failure criterion and equivalent bulk strength properties; generating subsurface bulk strength properties in strata of interest; utilizing the predictive algorithms and subsurface bulk strength properties to calculate equivalent anisotropic strength properties in subsurface strata of interest; and utilizing predicted anisotropic rock strength for subsurface strata of interest in numerical models of the subsurface.
In still yet another embodiment, a system associated with the drilling of a wellbore is described. The system comprises: a processor; memory coupled to the processor; an application stored in the memory and configured to: compile a database of core-based strength measurements over a range of orientations and confining pressures for a plurality of fine-grained rock formations; utilize an automated surface fitting procedure to determine a common anisotropic failure criterion for all database rocks; develop predictive algorithms relating anisotropy parameters in failure criterion and equivalent bulk strength properties; generate subsurface bulk strength properties in strata of interest; utilize the predictive algorithms and subsurface bulk strength properties to calculate equivalent anisotropic strength properties in subsurface strata of interest; utilize predicted anisotropic rock strength for subsurface strata of interest in numerical models of the subsurface to develop an anisotropic wellbore stability model; and a monitor coupled to the processor and configured to present aspects of the anisotropic wellbore stability model.
In one or more embodiments, the method or system may include various enhancements. For example, observed anisotropic strength values may be determined based on the database of core-based strength measurements; predicted strength values may be determined based on the predictive algorithms; and statistical tools may be used to evaluate a goodness of fit between the observed anisotropic strength values and the predicted strength values; and/or the predictive algorithms may be modified based on the goodness of fit. Further, a database of core-based strength measurements may be modified based on the goodness of fit.
The foregoing and other advantages of the present disclosure may become apparent upon reviewing the following detailed description and drawings of non-limiting examples of embodiments.
In the following detailed description section, the specific embodiments of the present disclosure are described in connection with preferred embodiments. However, to the extent that the following description is specific to a particular embodiment or a particular use of the present disclosure, this is intended to be for exemplary purposes only and simply provides a description of the exemplary embodiments. Accordingly, the disclosure is not limited to the specific embodiments described below, but rather, it includes all alternatives, modifications, and equivalents falling within the true spirit and scope of the appended claims.
Various terms as used herein are defined below. To the extent a term used in a claim is not defined below, it should be given the broadest definition persons in the pertinent art have given that term as reflected in at least one printed publication or issued patent. Further, various terms and reference characters that are the same are utilized in the equations and associated description within the disclosure.
The present disclosure describes methods and systems to enhance subsurface models. The present techniques relates generally to the field of hydrocarbon production, and more particularly to the effect of anisotropic formation strength on the drilling of a stable wellbore along a given locus. To enhance these production operations, a petrophysical method may be based on geophysical wireline logs to predict the shear strength anisotropy of sub-surface fine-grained rock formations (e.g., shales) without having to core and plug the sub-surface formation of interest to directly measure shear strength anisotropy in the laboratory. Then, the predicted anisotropic strength functions are used in numerical models of the subsurface, including analytical and computational methods for optimizing mud weight design in high-angle/horizontal extended reach wells in which inclinations approach may be bed-parallel.
In conventional techniques, the direct measurement of anisotropy may be performed using core samples. Knowledge of shale anisotropy has been significantly advanced by the geophysics community through the analysis of seismic wave propagation (e.g., incorporating anisotropy into migration algorithms to provide for proper positioning of reflectors, [6]) such that the majority of experimental and theoretical investigations have focused on determining elastic anisotropy in shale formations. See, e.g., Cholach, P. Y. and D. R. Schmitt 2003. Intrinsic anisotropy of shales. In Proceedings of the Annual Meeting of the Society of Exploration Geophysicists, Dallas, 26-31 Oct. 2003 and Sondergeld, C. H. and C. S. Rai. 2011. Elastic anisotropy of shale. The Leading Edge March 2011. By comparison with this relatively advanced understanding of elastic anisotropy, research into shale strength anisotropy has received much less attention. This observation may partly reflect the fact that many more test samples are required to quantify strength as opposed to elastic anisotropy, which is discussed further below. Nevertheless, numerous experimental studies have sought to investigate anisotropic strength of rocks with strong planar geologic fabrics including: slate; schist and gneiss; shale and mudstone; siltstone; sandstone; limestone; and coal. See, e.g., McLamore, R. T. and K. E. Gray. 1967. A strength criterion for anisotropic rocks based upon experimental observations. In Proceedings of the Annual Meeting of the American Institute of Mining, Metallurgical and Petroleum Engineers, Los Angeles, 19-23 Feb. 1967; Donath, F. A. 1964. Strength variations and deformational behavior in anisotropic rock. In State of Stress in the Earth's Crust, ed. W. R. Judd, 281-297; Attewell, P. B. and M. R. Sandford. 1974. Intrinsic shear strength of a brittle, anisotropic rock—I, Experimental and mechanical interpretation. Int. J. Rock Mech. Min. Sci. 11: 423-430; Cho, J.-W., H. Kim, K.-B. Min and S. Jeon. 2011. An experimental study on deformation and strength anisotropy of transversely isotropic rocks in Korea. In Proceedings of the 12th International Congress on Rock Mechanics, Beijing, 18-21 Oct. 2011; Nasseri, M. H. B., K. S. Rao and T. Ramamurthy. 2003. Anisotropic strength and deformational behavior of Himalayan schists. Int. J. Rock Mech. Min. Sci. 40(1): 3-23; Shea, W. T. and A. K. Kronenberg. 1993. Strength and anisotropy of foliated rocks with varied mica contents. J. Struct. Geol. 15(9): 1097-1121; Niandou, H., J. F. Shao, J. P. Henry and D. Fourmaintraux. 1997. Laboratory investigation of the mechanical behavior of Tournemire shale. Int. J. Rock Mech. Min. Sci. 34(1): 3-16; Chenevert, M. E. and C. Gatlin. 1965. Mechanical anisotropies of laminated sedimentary rocks. Soc. Petrol. Eng. J. 5: 67-77; Ibanez, W. D. and A. K. Kronenberg. 1993. Experimental deformation of shale: Mechanical properties and microstructural indicators of mechanisms. Int. J. Rock Mech. Min. Sci. & Geomech. Abstr. 30(7): 723-734; Ajalloeian, R. and G. R. Lashkaripour. 2000. Strength anisotropies in mudrocks. Bull. Eng. Geol. Env. 59: 195-199; Attewell, P. B. and I. W. Farmer. 1976. Principles of Engineering Geology. 1st ed. London: Chapman & Hall; Whittles, D. N., E. Yasar, D. J. Reddish and P. W. Lloyd. 2002. Anisotropic strength and stiffness properties of some UK Coal Measure siltstones. Q. J. Eng. Geol. Hydro. 35: 155-166; McGill, G. E. and J. A. Raney. 1970. Experimental study of faulting in an anisotropic, inhomogeneous dolomitic limestone. Geol. Soc. America Bull. 81: 2949-2958; and Pomeroy, C. D., D. W. Hobbs and A. Mahmoud. 1971. The effect of weakness-plane orientation on the fracture of Barnsley Hards by triaxial compression. Int. J. Rock Mech. Min. Sci. 8: 227-238.
Conventional techniques for quantifying both elastic and strength anisotropy involve direct measurement of mechanical properties utilizing orientated core samples. The directional character of the deformational properties of anisotropic media is commonly analyzed in terms of the theory of elasticity using the 3D generalized form of Hooke's law. Expressed via the compliance matrix, the rock then has twenty-one independent elastic constants. See, e.g., Jaeger, J. C., N. G. W. Cook and R. W. Zimmerman. Fundamentals of Rock Mechanics. 4th Edn. Blackwell Publishing Ltd, MA. However, for most practical cases, anisotropic rocks are usually modeled as transversely isotropic media, implying that at each point in the rock there is an axis of rotational symmetry and that the rock has isotropic properties in the plane normal to that axis, in which case only five independent elastic constants are required to fully describe elastic deformability. As shown in diagram 200 of
While plasticity (e.g., deformation beyond the limit of purely elastic strain), theories for anisotropic geologic materials have also been explored, but by comparison with elastic anisotropy it is less clear how such theories should be utilized to define rational material properties testing procedures. See e.g., Pariseau, W. G. 1968. Plasticity theory for anisotropic rocks and soil. In Proceedings 10th US Symposium on Rock Mechanics, Austin, 20-22 May 1968. As shown in the diagram 250 of
Geographically-diverse examples of anisotropic rock strength models derived from direct core measurement being used to optimize drilling programs include: Gulf of Mexico; Canada; Venezuela and Colombia; UK continental shelf; and the Middle East. See, e.g., Brehm, A., C. Ward, D. Bradford and G. Riddle. 2006. Optimizing a deepwater subsalt drilling program by evaluating anisotropic rock strength effects on wellbore stability and near-wellbore stress effects on the fracture gradient. In Proceedings of the IADC/SPE Drilling Conference, Miami, 21-23 Feb. 2006; Lang, J., S. Li and J. Zhang. 2011. Wellbore stability modeling and real-time surveillance for deepwater drilling to weak bedding planes and depleted reservoirs. In Proceedings of the IADC/SPE Drilling Conference, Amsterdam, 1-3 Mar. 2011; Gallant, C., J. Zhang, C. A. Wolfe, J. Freeman, T. Al-Bazali and M. Reese. 2007. Wellbore stability considerations for drilling high-angle wells through finely laminated shale: A case study from Terra Nova. In Proceedings of the SPE Annual Technical Conference, Anaheim, 20-22 Nov. 2007; Aadnoy, B., G. Hareland, A. Kustamsi, T. de Freitas and J. Hayes. 2009. Borehole failure related to bedding plane. In Proceedings 43rd US Rock Mechanics Symposium, Asheville, 28 Jun.-1 Jul. 2009; Willson, S. M., N. C. Last, M. D. Zoback and D. Moos. 1999. Drilling in South America: A wellbore stability approach for complex geologic conditions. 1999. In Proceedings of the SPE Latin American and Caribbean Petroleum Engineering Conference, Caracas, 21-23 Apr. 1999; Narayanasamy, R., D. Barr and A. Milne. 2010. Wellbore-instability predictions within the Cretaceous mudstones, Clair field, west of Shetlands. SPE Drilling & Completion December 2010: 518-529; and Grini, M., S. Ottesen, H. S. Al Junaibi, J. Jorden, D. A. Boyd, W. Martin, B. Crawford, B. Alramahi and S. Gosavi. 2012. Extended reach drilling optimization in layer A shale through geomechanical characterization. In Proceedings of the SPE Abu Dhabi International Petroleum Exhibition & Conference, Abu Dhabi, 11-14 Nov. 2012.
In addition to the costly and time consuming nature of deriving shale strength anisotropy from direct core measurement, another distinct disadvantage relates to the potential for inducing damage in such low permeability, fine-grained rocks due to rapid exhumation of core from great depth to the surface. See, e.g., Holt, R. M. 1994. Effects of coring on petrophysical measurements. In Proceedings of the Society of Core Analysts Annual Technical Conference, Stavangar, 1994. The associated permanent microstructural alteration (resulting from diverse mechanisms including in situ stress and pore pressure release, temperature reduction, desiccation, exposure to non-native drilling fluids and poor handling and preservation techniques) can lead to spurious measurement of shale strength anisotropy unrepresentative of that associated with subsurface conditions. See, e.g., Santarelli, F. J. 1994. Similarities and differences between wellbore stability and sand production. In North Sea Oil and Gas Reservoirs—III, 30 Nov.-2 Dec., eds. Aasen, J. O., E. Berg, A. T. Buller, O. Hjelmeland, R. M. Holt, J. Kleppe and O. Torsxter, 331-338. Dordrecht: Kluwer.
Experimentally determined strength data, σ1=f(σ3, β), may then be used to calibrate various anisotropic failure criteria. Numerous such criteria have been proposed to capture the deformation behavior of anisotropic rocks under different loading conditions. Duveau et al summarize many of these criteria and categorize them into distinct groups. See, e.g., Duveau, G., J. F. Shao and J. P. Henry, 1998, Assessment of some failure criteria for strongly anisotropic rocks. Mech. Cohes-frict. Mater. 3: 1-26. The first group includes discontinuous criteria such as the single plane of weakness (SPW) model. In criteria of this form, there is a distinction between failure on a weak plane (such as a bedding or lamination plane) and failure of the bulk material (where the failure surface cuts through such planes of weakness). With this theory the strength of the rock is uniform with respect to orientation if the failure does not occur on a weak plane and variable if it does. Results for this simplest example of planar anisotropy can be based on the well-known two-dimensional theory for sliding on a plane of weakness whose normal makes an orientation angle β with the maximum principal stress, σ1-direction. See, e.g., Jaeger. If the plane of weakness has strength properties τ0w is the cohesive strength and φw is the angle of internal friction, the criterion for failure in the plane of weakness is given by the well-known linear Coulomb criterion. See, e.g., Coulomb, C. A. 1773. Application des regles de maxima et minima a quelques problemes de statique relatifs a 1′ architecture. [application of the rules of maxima and minima to some problems of statics related to architecture]. Acad. Roy. Sci. Mem. Math. Phys. 7: 343-382. The linear Coulomb criterion is shown by the following Equation (1).
|τ|=τ0w+σn tan φw Equation (1)
where the shear stress τ and normal stress σn on the weakness plane are determined by the maximum and minimum principal stresses (σ1 and σ3 respectively) and the orientation angle β according to the following Equations 2 and 3:
Substituting Equation (2) and Equation (3) into Equation (1) enables determination of the maximum principal stress difference (σ1−σ3) sustainable by the plane of weakness before failure, for a given orientation angle β: This is provided in Equation (4) below.
(σ1−σ3)=2(τ0w+σ3 tan φw)/((1−tan φw cot β)sin 2β) Equation (4)
As given by Equation (4), the stress difference at failure (σ1−σ3) tends to infinity as β→½ π and as β→φw. Between these values, failure of the plane of weakness is possible and has a minimum value at βminw given as given by Equation (5):
This is the orientation most favorable for failure. However, for unfavorable values of orientation angle β, failure may take place through the bulk material on a plane which cuts across the plane of weakness. If the Coulomb criterion is again utilized to describe bulk failure with material properties τ0b is the cohesive strength in the bulk material and φb is the angle of internal friction in the bulk material, failure occurs based on the same criterion as the weak plane (Equation (1)), but with the maximum principal stress difference (σ1−σ3) sustainable by the bulk material before failure now given by Equation (6):
The induced shear failure plane within the bulk material is orientated such that its outward normal forms an acute angle βb with the maximum principal stress direction where:
The SPW model therefore predicts that, for a given value of σ3, the magnitude of σ1 required to induce failure somewhere within the anisotropic rock is equal to the lesser of the two values as defined by Equation (4) and Equation (6). If, for a given orientation angle β, the value given by Equation (4) is less than that given by Equation (6), failure occurs along the plane of weakness. Conversely, if the value given by Equation (4) is greater than that given by Equation (6), the induced shear failure plane within the intact rock cuts through the plane of weakness at an orientation as given by Equation (7).
It is evident from Equation (4) that the variation of (σ1−σ3) with orientation angle β has the same form for each of the values of σ3, but with a different multiplying factor. The nature of this variation is illustrated in the graph 300 of
In the graph 300, the concave-upwards portions of each curve, such as portions 306a, 306b, 306c and 306d, correspond to the criterion given by Equation (4) and represent β-dependent failure along the plane of weakness. The horizontal portions of each curve are derived from Equation (6) and represent O-independent failure within the intact rock along a plane defined by Equation (7).
While the advantages of the SPW model are that it is theoretically-based and fully-described using only four material parameters (τ0b, φb, τ0w, φw) the considerable disadvantage of this criterion is illustrated schematically in
In this flow diagram 400, various blocks indicate differentiating criteria to classify the types of anisotropy. For example, the types of anisotropy are described in block 402, which divide into a cleavage type in block 404 and a bedding type in block 412. The cleavage type may exhibit a maximum strength at orientation angle β equal to 0°, which is utilized to divide the types into a U Type in block 407 (e.g., associated with profile 417) and an Undulatory Type in block 410 (e.g., associated with profile 418). The bedding type may include exhibiting maximum strength at orientation angle β equal to 90°, as shown in block 414. Then, the shoulder type may be determined as shown in block 416 (e.g., associated with profile 419).
These classifications may be provided as a profile, such as profiles 417, 418 and 419, which compare the strength versus orientation for each profile. Although some experimental results (particularly for slates) show general agreement with the SPW model, as shown in profile 417, many anisotropic sedimentary rocks exhibit pronounced asymmetry in the shape of their strength σ1 versus orientation angle β profiles which cannot be satisfactorily modeled using the SPW criterion, as shown in profiles 418 and 419.
A second category of anisotropic failure criteria distinguished by Duveau includes those models which describe a continuous variation in strength with orientation such that no explicit distinction is made between failure along planes of weakness and failure of the bulk material. In one such formulation, Jaeger generalized the well-known Coulomb theory for shear failure to the case in which cohesive strength τ0 was not constant but varied with orientation angle β. See, e.g., Jaeger, J. C. 1960. Shear failure of anisotropic rocks. Geol. Mag. 97: 65-72. The equation is shown below in Equation (8):
τ0=A*−B*cos 2(β−βmin) Equation (8)
where A* and B* are fitting parameters and βmin is the orientation angle at which cohesive strength exhibits a functional minimum. In this variable cohesive strength (VCS) model, the associated internal friction angle φ is not directly specified, but is assumed to be independent of direction. Because the VCS model uses a periodic (cosine) function to fit cohesive strength data and Coulomb theory predicts that βmin=45+φ/2 (Equation (7)), such that βmin≧45°, it follows that the magnitude of τ0 at the orientation angle β=0° is greater than that at the orientation angle β=90°. This is a constraint placed by the model that does not represent the observed variation in strength for many anisotropic rock types. McLamore et al. expanded upon the VCS model by postulating additional periodic variation in friction angle with orientation and allowing both cohesion and friction angle to have dependence on a cosine function raised to a power n, as shown in Equation (9):
(τ0,φ)=A1,2−β1,2[cos 2(β−βmin)]n Equation (9)
where A1 and B1 are τ0 and φ-specific empirical constants that describe the behavior of τ0 and φ over the range 0°≦β≦βmin and A2 and B2 the equivalent empirical constants for the range βmin≦β≦90°. The terms A1,2 and B1,2 represent A1 or A2 and B1 or B2 for different equations of this terms. It is therefore evident that Equation (9) requires a total of twelve empirically determined constants to fully describe strength anisotropy. Both the SPW and VCS anisotropic failure criteria detailed above can only be determined for a specific subsurface lithology of interest through orientated triaxial compression testing of recovered core, to directly measure σ1 equal f(σ3, β) data that can subsequently be used to calibrate each model. Although the modified VCS model of McLamore et al. offer increased congruency with experimental data, utilization of this criterion for describing strength anisotropy in a subsurface lithology of interest for which no core was available for geomechanics testing, involves having to predict the requisite twelve empirical constants which poses a considerable challenge.
With regard to predictive techniques, as rock strength is directly quantified from laboratory testing, core recovery issues dictate that many geomechanical problems have to be addressed in the absence of samples. Accordingly, there exists a plethora of empirical equations in the literature for relating mechanical strength properties to other physical properties potentially derivable from geophysical wireline logs, which can be utilized when no core is available for testing. Mechanical properties evaluation is primarily based on application of the most appropriate of these pre-determined predictors to logs derived from the subsurface strata of interest, calibrated where possible against limited core measurements often from offset wells.
Most empirical predictors are explicit to particular rock types (e.g., sandstone, carbonate or shale) and the majority relates rock strength to either acoustic velocities, the equivalent dynamic elastic moduli, or porosity. Predicted strength parameters are predominantly those that describe the linear Coulomb-type failure criterion (as noted in Coulomb), specifically unconfined compressive strength (UCS), cohesive strength τ0 and internal friction angle φ. Comprehensive summaries of numerous such empirical equations sourced from the literature that relate rock strength parameters to other physical properties are given in Change and Khakar. See, e.g., Chang, C., M. D. Zoback and A. Khaksar. 2006. Empirical relations between rock strength and physical properties in sedimentary rocks. J. Pet. Sci. & Eng. 51: 223-237 and Khaksar, A., P. G. Taylor, Z. Fang, T. Kayes, A. Salazar and K. Rahman. 2009. Rock strength from core and logs: where we stand and ways to go. Proceedings of the SPE EUROPEC/EAGE Annual Conference & Exhibition, Amsterdam, 8-11 Jun. 2009.
By way of an example for shales, several references describe robust negative power law relationships between UCS and porosity, while others describe positive correlations between UCS and compressional sonic velocity. However, as strength is controlled by complex rock fabric on the microscale, oftentimes core-derived rock strength data can show considerable variation as a function of single input parameters, such as velocity, elastic modulus or porosity, such that most empirical predictors are not sufficiently generic to fit all published strength versus physical property data. See, e.g., Lashkaripour, G. R., and M. B. Dusseault. 1993. A statistical study on shale properties: relationship among principal shale properties. In Proceedings of the Conference on Probabilistic Methods in Geotechnical Engineering, Canberra, 10-12 February; Horsrud, P. 2001. Estimating mechanical properties of shale from empirical correlations. SPE Drilling & Completion, 16: 68-73; Lal, M. 1999. Shale stability: drilling fluid interaction and shale strength. Proceedings of the SPE Latin American and Caribbean Petroleum Engineering Conference, Caracas, 21-23 Apr. 1999; and Ingram, G. M. and J. L. Urai. 1999. Top-seal leakage through faults and fractures: the role of mudrock properties. In, Muds and Mudstones: Physical and Fluid Flow Properties. Geological Society London, Special Publications, 158: 125-135.
Multi-variate algorithms can offer improved congruency with experimental data, for example Coulomb strength parameters in sandstone-to-shale siliciclastic rocks can be predicted from algorithms relating τ0 and φ to porosity and lithology, the latter defined by total clay content which can be directly quantified by the X-ray diffraction XRD technique or derived from petrophysical analysis of wireline logs. See, e.g., United States Patent Application No. 2011/0015907, Publication Date: Jan. 20, 2011. Another shale-specific strength predictor correlates shale strength to surface area as derived from the dielectric constant measurement (DCM) technique. See, e.g., Leung, P. K. and Steiger, R. P. 1992. Dielectric constant measurements: a new, rapid method to characterize shale at the wellsite. Proceedings of the IADC/SPE Drilling Conference, New Orleans, 18-21 Feb. 1992. Shale surface area represents a measure of the total fraction of swelling clays and DCM can be readily performed on drilling cuttings. An extensive shale strength database acquired through consolidated undrained triaxial testing may be compiled to develop the predictive regression which relates shale strength (octahedral shear stress at failure) to mean effective stress and shale surface area. See, e.g., Steiger, R. P. and P. K. Leung. 1992. Quantitative determination of the mechanical properties of shales. SPE Drilling & Completion, 7: 181-185.
Implicit in the above predictive techniques is that strength properties are derived from databases of laboratory measurements in which core testing is conducted on samples with any planar geologic fabric orientated perpendicular to the direction of maximum compressive stress (e.g., at orientation angle β=0°), such that predicted mechanical properties relate specifically to bulk (trans-laminar) strength in which induced failure planes cut across bedding laminations, so that these techniques offer no predictive capability for quantifying strength anisotropy associated with the presence of such weakness planes.
Compared with the relatively numerous bulk strength predictive techniques that have been developed via calibration with laboratory data as described above, techniques for predicting strength anisotropy without recourse to core testing are relatively few.
As an example, Willson et al proposed a formulation for describing unconfined compressive strength (UCS) variation relative to bedding orientation, in which the level of strength anisotropy is characterized via two constants in Equations (10) and (11):
k
1=UCS∥/UCS⊥ Equation (10)
k
2=UCSmin/UCSmax Equation (11)
In Equation (10) UCS∥ is the strength with bedding parallel to the σ1-direction (β=90°) and UCS⊥ is the strength with bedding perpendicular to the σ1-direction (e.g., orientation angle β=0°) and in Equation (11) UCSmin is the minimum strength for any orientation and UCSmax the maximum strength for any orientation. To derive a smoothly varying strength profile with orientation angle the following assumptions were made: UCS⊥=UCSmax; UCSmin occurs at orientation angle β=45°. Based on these assumptions the following orientation-dependent unconfined compressive strength variation was proposed:
UCS(β)=UCSmax(cos φ+k1 sin β)(1−sin β cos β)[1−2 sin β cos β(1−4k2/√{square root over (2)}(1+k1))] Equation (12)
Equation (12) describes strength anisotropy from knowledge of only three direct measurements (UCS∥, UCS⊥ and UCSmin). While UCS⊥ is potentially predictable from empirical correlations with geophysical wireline measurements (for example porosity or acoustic velocities, see, e.g., Chang) some recovered core from the subsurface strata of interest may be available for direct laboratory measurement of UCS∥ and UCSmin. In addition, knowledge of unconfined compressive strength (σ3=0) alone is insufficient for determining anisotropic strength at elevated confining pressures (e.g., σ3>0), which requires additional information regarding the magnitude of the associated internal friction angle φ. Another disadvantage relates to the form of this orientation-dependent unconfined compressive strength variation, which due to the imposed assumptions, is constrained to always predict UCSmin at orientation angle β=45° and UCS∥ at orientation angle (β=90°) less than UCS⊥ at orientation angle (β=0°) for all fine-grained rock types, which is not in agreement with many experimental observations, as noted above in the discussion associated with
Smith and Goldman describe a technique for assaying the compressive strength of rock, primarily based upon a database of laboratory measurements from which strength versus porosity correlations are extracted. See, e.g., U.S. Pat. No. 5,767,399. While the technique mainly addresses unconfined compressive strength, the present techniques expand this to accommodate other factors that impact rock strength, such as confining stress, bedding plane orientation and temperature. Due to the predictive methodology underlying the technique, ideally calibration involving laboratory testing of pure lithologies (e.g., clean sandstone and/or shale) recovered from the well of interest is required. However, lithologically similar samples (perhaps from offset wells) may suffice. The technique has limitations when used to predict strength in mixed lithologies, such as shaly sandstones or sandy shales, as it relies on a weighted average of the pure end member components which can lead to significant predictive uncertainty. In the present techniques, the strength versus orientation relationship is represented by a curve corresponding to an equation of the form in Equation (13):
σco=Se[σumax+Δσmax(σ3/σ3max)β](1−comax)+(1−Se)[σumin+Δσmin(σ3/σ3max)β](1−comin) Equation (13)
where, for 0°<β≦βc:
The advantage of the present techniques is that knowing the general form of Equation (13) as well as the general form of its representation as a curve (a conjunction of portions of two different sine waves) provides characterization of the full-range of experimentally observed orientation-dependent compressive strength variations, including asymmetric profiles, as shown schematically in profiles 418 and 419 of
In the present techniques, the strength prediction may be limited to correlations using porosity alone as data input, and precludes the use of other log-based strength correlations that may be more appropriate for a given lithology or region (generally field-, geologic age- or depth-specific) such as acoustic velocity and DCM predictors. Disadvantages to utilizing existing porosity-based shale strength predictors are detailed below.
Direct porosity measurement in fine-grained rock formations, such as shales, is difficult to accomplish using traditional core analysis methods (such as, Mercury injection, Dean Stark solvent extraction or Boyle's law Helium porosimetry) due to the large impact of coring-induced microcracking and fluid invasion effects. To overcome this, typical analysis is performed on crushed rock samples to eliminate coring-induced artifacts while maintaining representation of pore structure. However, measured porosity can still exhibit considerable dependency on the measurement technique due to difficulties in discriminating free and bound water. See, e.g., Handwerger, D. A., R. Suarez-Rivera, K. I. Vaughn and J. F. Keller. 2011. Improved petrophysical core measurements on tight shale reservoirs using retort and crushed samples. In Proceedings of the SPE Annual Technical Conference and Exhibition, Denver, 30 Oct.-2 Nov. 2011. Also, petrophysical techniques for deriving porosity from geophysical wireline logs are typically only calibrated against core measurements derived from the reservoir section, so that log-based porosity in intervening shale sections may be prone to considerable error without additional advanced interpretation.
The present techniques involve predicting the degree of shear strength anisotropy encountered in subsurface fine-grained rock formations (particularly shales) without recourse to coring the subsurface strata of interest and performing laboratory strength measurements. The present techniques use a pre-established database of strength measurements derived from orientated core testing to develop predictive algorithms relating parameters describing the degree of strength anisotropy associated with diverse geologic planes of weakness (such as bedding, lamination or cleavage planes) to bulk strength properties easily derived from geophysical wireline logs using a variety of established petrophysical techniques.
With regard to planar geologic fabrics in sedimentary rocks, the terms “bedding” and “lamination” refer to stratification thicker than 1 centimeter (cm) and thinner than 1 cm, respectively. However, both may be defined by changes in sediment grain-size, color or mineralogical composition, such that the terms may be used interchangeably to imply planes of weakness in sedimentary rocks resulting from changes in the pattern of sedimentation. Similarly, with regard to foliated textures in metamorphic rocks, the terms “slaty cleavage”, “schistosity” and “gneissic mineral banding” refer to specific fabrics associated with increasing metamorphic grade (pressure and temperature conditions). However, they represent the preferential alignment of minerals within a directional stress field, such that the terms may be used interchangeably to imply planes of weakness in metamorphic rocks resulting from preferential alignment of platy minerals such as chlorite, muscovite and biotite, or the segregation banding of granular and platy mineral concentrations.
As long as the subsurface strata of interest are broadly lithologically similar to a particular fine-grained rock type characterized in the database, then the pre-established predictive algorithms can be utilized to generate subsurface anisotropic strength parameters, without direct incorporation of orientated core measurements from the subsurface strata of interest within the database. For example, if the pre-established database contains mainly orientated strength measurements conducted on sedimentary shales, then it cannot be used to generate subsurface anisotropic strength parameters in metamorphic slates, and vice versa.
The present techniques may involve obtaining relationships between planes of weakness and bulk material strength for a plurality of rock samples; obtaining bulk strength of target rock (e.g., well logs, pre-established techniques, etc.); calculating the planes of weakness for the target rock (strength anisotropy); and using the results to enhance models, which may be used for wellbore stability analysis or enhancing recovery.
The method for predicting strength anisotropy in fine-grained subsurface formations may include various steps. For example, the method may include compiling a database of core-based strength measurements over a range of orientations and confining pressures for a plurality of fine-grained rock formations; utilizing an automated surface fitting procedure to determine a common 3D anisotropic failure criterion for all database rocks; developing predictive algorithms relating anisotropy parameters in 3D failure criterion and equivalent bulk strength properties; generating subsurface bulk strength properties in strata of interest using established petrophysical techniques; utilizing pre-established predictive algorithms and subsurface bulk strength properties to calculate equivalent anisotropic strength properties in subsurface strata of interest; utilizing predicted anisotropic rock strength for subsurface strata of interest in numerical models of the subsurface.
As an example, a workflow is shown in the flow chart 500 of
The present techniques provide various enhancements over other methods. For example, the present techniques predict the degree of shear strength anisotropy encountered in subsurface fine-grained rock formations (particularly shales) without recourse to coring the subsurface strata of interest and performing laboratory strength measurements. Also, the present techniques predict anisotropic strength functions that can subsequently be used in numerical models of the subsurface, including analytical and computational methods for optimizing mud weight design in high-angle/horizontal extended reach wells in which inclinations can approach bed-parallel. The present techniques also utilize bulk strength properties to predict anisotropic failure criteria without any restrictions on the source of the bulk strength properties (directly measured on core, or estimated from geophysical wireline logs or cuttings/cavings analyses). Further, the techniques are not restricted to any one specific anisotropic failure criterion (predicts any discontinuous or continuous model) or anisotropic rock type (predicts any fine-grained rock with either a planar sedimentary fabric or foliated metamorphic texture). Moreover, the present techniques incorporate standard analysis of variance statistical tools in the predictive workflow. Further still, the techniques provide a simple index for comparing the magnitude of strength anisotropy between diverse fine-grained rock types. To further describe the present techniques, various aspects are described further in
As noted above, the flow chart 500 in
Strength data included in the database for each fine-grained rock formation includes maximum axial compressive strength σ1 measured over a range of confining pressure σ3 magnitudes at different orientation angles β (defined as the angle between the normal to the planar layering and the direction of maximum compressive stress σ1). Preferably, for each fine-grained rock formation, a series of orientated, triaxially compressive strength tests should be conducted for at least seven β orientations in 15° increments from the vertical (orientation angle β=0°) to the horizontal (orientation angle β=90°) direction. The minimum number of orientation angles tested should be three (orientation angles β=0°, 45 or 60°, 90°). Preferably, for each fine-grained rock formation, orientated, triaxially compressive strength tests should be conducted for at least four constant confining pressure σ3 magnitudes for each orientation angle β. The minimum number of confining pressures tested should be three. Therefore, preferably, each fine-grained rock formation in the database should have twenty-eight orientated triaxially compressive strength tests, and for a fine-grained rock formation to be included in the database, it should have a minimum of at least nine such associated tests.
A schematic illustration of a database is shown in
As laboratory measurement of compressive failure strength is particularly difficult to accomplish successfully in fine-grained rocks due to their associated ultra-low permeabilities, proper methods should be used to quality control the experimental data before inclusion in the database. As one of the objectives is to predict rock strength anisotropy for wellbore stability applications associated with extended reach drilling through shale formations, 81% of all orientated triaxial compression tests in the illustrative database were conducted on fine-grained sedimentary rocks, of which 65% can be described as generic shales. The additional anisotropic sedimentary rock types and metamorphic slates are included for comparative purposes.
In block 504 of the predictive workflow in
τ0=UCS/2√{square root over (k)} Equation (14)
τ0(β)=A′·β3+B′·β2+C′·β+τ0|β=0 Equation (15)
tan φ=(k−1)/(2√{square root over (k)}) Equation (16)
τ(β)=χ·τ0|β=0+σn·tan
where:
χ=τ0(β)τ0|β=0=Aβ3+Bβ2+Cβ+1 Equation (18)
The anisotropic failure criterion of the present techniques has several distinct enhancements over the conventional approaches. The automated surface fitting procedure tends to minimize variability in measured strength not related to experimental variables (σ3 and β), but rather the result of between-sample material inhomogeneity. Further, the functional minimum in cohesive strength versus orientation at βmin as shown in
In block 506, the database fine-grained rock formations are used to develop empirical relationships between the anisotropy parameters describing the 3D failure criterion for each formation and their equivalent bulk strength properties. In this manner, the bulk strength properties refer generally to those properties associated with trans-laminar failure in which failure occurs within the bulk material by cutting across the planes of weakness, and specifically to those properties describing bulk failure at β=0° in which the direction of maximum applied compressive stress σ1 is coincident with the normal to the planar layering. This block 506 therefore utilizes the mechanical properties database to generate predictive algorithms that statistically relate the empirical constants A, B, C and
To quantify the magnitude of strength anisotropy, a cohesive strength reduction factor (CSR) representing the normalized strength difference associated with shearing across (trans-laminar) versus shear within (intra-laminar) the plane of anisotropy is calculated.
This factor may be represented by Equation (21):
such that, as noted in Equation (22):
τ0|β=βminmodel=τ0|β=0core·(1−CSR) Equation (22)
For CSR as defined in Equation (21) trans-laminar strength can be derived from conventional triaxial compression testing (β=0°) with loading perpendicular to bedding) whereas intra-laminar strength is derived from the modeled functional minimum at βmin as shown in
It is apparent that the slates exhibit the greatest CSR values, presumably due to the crystallographic preferred orientation of their constituent platy minerals. However, all other lithologies with planar sedimentary fabrics exhibit a systematic trend that may be fit using an asymmetric peak function for τ0|β=0≦45 MPa of the form of Equation (23):
where x=τ0|β=0, y=CSR and α0, α1 and α2 are empirical constants. As τ0|β=0 can be derived from geophysical wireline logs using a variety of predictive techniques (e.g. Chang and Khaksar).
The asymmetric peak function in
A,B,C=α
3+α4(τ0|β=0)+α5(CSR) Equation (24)
where α3, α4 and α5 are A, B and C and lithology (Type I, II and III) dependent empirical constants pre-determined through multi-variate regression of the database lithologies. A nomograph (graphical representation of Equation (24)) is shown in
For the database lithologies, the mean friction angle,
As φ|β=0 can be derived from geophysical wireline logs using a variety of predictive techniques (e.g. Chang and Khaksar),
In block 508, following development of empirical relationships between anisotropy parameters describing the selected 3D failure criterion and the equivalent bulk strength properties as detailed above, an additional step in the predictive workflow involves utilizing the well-known analysis of variance (ANOVA) collection of statistical tools to determine if the predictive algorithms so generated are statistically useful, that is unlikely to have occurred by chance alone. Goodness of fit statistics can be used to quantify discrepancy between observed values and the values expected under the model in question. For example, as a fit becomes better (e.g., as discrepancy between the observed and modeled values is lessened), the coefficient of determination approaches unity, the standard error decreases toward zero and the F-statistic tends towards infinity. If statistical analyses indicate a poor model fit to experimental observations, then the predictive workflow (as shown in
Alternatively (or in addition) to adjusting the database of core measurements as detailed above, the anisotropic failure criterion used to quantify orientation-dependent strength variation for all fine-grained rock formations can also be adjusted to increase statistical significance. For example, if ANOVA testing indicates that the predictive algorithms relating the polynomial VCS parameters A, B, C and
For example, the weakness plane cohesive strength, τ0w defined in Equation (1) for the discontinuous SPW model is positively correlated with bulk cohesive strength τ0|β=0 for all fine-grained database rocks with planar sedimentary fabrics.
In addition, statistically significant correlations between equivalent material parameters from different failure criteria can also be developed for predictive purposes. For example, the weakness plane angle of internal friction φw defined in Equation (1) for the discontinuous SPW model is positively correlated with the mean friction angle,
In yet another embodiment it is evident that, if statistically significant correlations can be generated for predicting the material parameters of more than one anisotropic failure criteria, then these failure criteria can be combined to increase congruency with experimental data. For example, the strength predicted by the continuous VCS model can be combined with the strength predicted by the discontinuous SPW model, where the resultant strength is the lesser of that predicted by the two models.
Further, in one or more embodiments, the present techniques may be utilized to produce hydrocarbons from a subsurface region. This subsurface region may be represented by a anisotropic rock strength model along an actual or proposed wellbore trajectory, which may be generated in accordance with the techniques noted above. For example, a method may include constructing a subsurface model via one or more of the various embodiments above. Then, the subsurface model may be utilized in hydrocarbon management activities to determine locations of resources, access the resources and produce hydrocarbons in a manner know to those skilled in the art. Disclosed aspects may be used in hydrocarbon management activities. As used herein, “hydrocarbon management” or “managing hydrocarbons” includes hydrocarbon extraction, hydrocarbon production, hydrocarbon exploration, identifying potential hydrocarbon resources, identifying well locations, determining well injection and/or extraction rates, identifying reservoir connectivity, acquiring, disposing of and/or abandoning hydrocarbon resources, reviewing prior hydrocarbon management decisions, and any other hydrocarbon-related acts or activities. The term “hydrocarbon management” is also used for the injection or storage of hydrocarbons or CO2, for example the sequestration of CO2, such as reservoir evaluation, development planning, and reservoir management. In one embodiment, the disclosed methodologies and techniques may be used to extract hydrocarbons from a subsurface region. In such an embodiment, geophysical information is utilized to define an anisotropic rock strength volume. Based at least in part on the anisotropic rock strength volume, hydrocarbon extraction may then be conducted to remove hydrocarbons from the subsurface region, which may be accomplished by drilling a well using oil drilling equipment. The equipment and techniques used to drill a well and/or extract the hydrocarbons are well known by those skilled in the relevant art. Other hydrocarbon extraction activities and, more generally, other hydrocarbon management activities, may be performed according to known principles.
Persons skilled in the technical field will readily recognize that in practical applications of the disclosed methodology, it is partially performed on a computer, typically a suitably programmed digital computer. Further, some portions of the detailed descriptions which follow are presented in terms of procedures, steps, logic blocks, processing and other symbolic representations of operations on data bits within a computer memory. These descriptions and representations are the means used by those skilled in the data processing arts to most effectively convey the substance of their work to others skilled in the art. In the present application, a procedure, step, logic block, process, or the like, is conceived to be a self-consistent sequence of steps or instructions leading to a desired result. The steps are those requiring physical manipulations of physical quantities. Usually, although not necessarily, these quantities take the form of electrical or magnetic signals capable of being stored, transferred, combined, compared, and otherwise manipulated in a computer system.
It should be borne in mind, however, that all of these and similar terms are to be associated with the appropriate physical quantities and are merely convenient labels applied to these quantities. Unless specifically stated otherwise as apparent from the following discussions, it is appreciated that throughout the present application, discussions utilizing the terms such as “processing” or “computing”, “calculating”, “determining”, “displaying”, “copying,” “producing,” “storing,” “adding,” “applying,” “executing,” “maintaining,” “updating,” “creating,” “constructing” “generating” or the like, refer to the action and processes of a computer system, or similar electronic computing device, that manipulates and transforms data represented as physical (electronic) quantities within the computer system's registers and memories into other data similarly represented as physical quantities within the computer system memories or registers or other such information storage, transmission or display devices.
Embodiments of the present techniques also relate to an apparatus for performing the operations herein. This apparatus may be specially constructed for the required purposes, or it may comprise a general-purpose computer selectively activated or reconfigured by a computer program stored in the computer (e.g., one or more sets of instructions). Such a computer program may be stored in a computer readable medium. A computer-readable medium includes any mechanism for storing or transmitting information in a form readable by a machine (e.g., a computer). For example, but not limited to, a computer-readable (e.g., machine-readable) medium includes a machine (e.g., a computer) readable storage medium (e.g., read only memory (“ROM”), random access memory (“RAM”), magnetic disk storage media, optical storage media, flash memory devices, etc.), and a machine (e.g., computer) readable transmission medium (electrical, optical, acoustical or other form of propagated signals (e.g., carrier waves, infrared signals, digital signals, etc.)).
Furthermore, as will be apparent to one of ordinary skill in the relevant art, the modules, features, attributes, methodologies, and other aspects of the invention can be implemented as software, hardware, firmware or any combination of the three. Of course, wherever a component of the present invention is implemented as software, the component can be implemented as a standalone program, as part of a larger program, as a plurality of separate programs, as a statically or dynamically linked library, as a kernel loadable module, as a device driver, and/or in every and any other way known now or in the future to those of skill in the art of computer programming. Additionally, the present techniques are in no way limited to implementation in any specific operating system or environment.
As an example, a computer system may be utilized and configured to implement one or more of the present aspects. The computer system may include a plurality of processors; memory in communication with the processors; and a set of instructions stored on the memory and accessible by the processors, wherein the set of instructions, when executed, are configured to: compile a database of core-based strength measurements over a range of orientations and confining pressures for a plurality of fine-grained rock formations; utilize an automated surface fitting procedure to determine a common anisotropic failure criterion for all database rocks; develop predictive algorithms relating anisotropy parameters in failure criterion and equivalent bulk strength properties; generate subsurface bulk strength properties in strata of interest; utilize the predictive algorithms and subsurface bulk strength properties to calculate equivalent anisotropic strength properties in subsurface strata of interest; utilize predicted anisotropic rock strength for subsurface strata of interest in numerical models of the subsurface. In certain embodiments, the set of instructions may perform the different aspects in the methods noted above or the algorithm noted above.
As an example, the techniques discussed herein may be implemented on a computing device, such as that shown in
The CPU 2002 may execute the various logical instructions according to the disclosed techniques. For example, the CPU 2002 may execute machine-level instructions for performing processing according to the exemplary operational flow described above. As a specific example, the CPU 2002 may execute machine-level instructions for performing the methods of
The computer system 2000 may also include random access memory (RAM) 2006, which may be SRAM, DRAM, SDRAM, or the like. The computer system 2000 may include read-only memory (ROM) 2008 which may be PROM, EPROM, EEPROM, or the like. The RAM 2006 and the ROM 2008 hold user and system data and programs, as is well known in the art. The programs may include code stored on the RAM 2006 that may be used for modeling geologic properties with homogenized mixed finite elements, in accordance with embodiments of the present techniques.
The computer system 2000 may also include an input/output (I/O) adapter 2010, a communications adapter 2022, a user interface adapter 2024, and a display adapter 2018. The I/O adapter 2010, user interface adapter 2024, and/or communications adapter 2022 may, in certain embodiments, enable a user to interact with computer system 2000 to input information. Further, the computer system 2000 may also include a graphical processing unit (GPU(s)) to enhance the graphical processing of the computer system 2000.
The I/O adapter 2010 may connect the bus 2004 to storage device(s) 2012, such as one or more of hard drive, compact disc (CD) drive, floppy disk drive, tape drive, flash drives, USB connected storage, etc. to computer system 2000. The storage devices may be used when RAM 2006 is insufficient for the memory requirements associated with storing data for operations of embodiments of the present techniques. For example, the storage device 2012 of computer system 2000 may be used for storing such information as computational meshes, intermediate results and combined data sets, and/or other data used or generated in accordance with embodiments of the present techniques.
The communications adapter 2022 is adapted to couple the computer system 2000 to a network (not shown), which may enable information to be input to and/or output from the system 2000 via the network, for example, the Internet or other wide-area network, a local-area network, a public or private switched telephone network, a wireless network, or any combination of the foregoing. The user interface adapter 2024 couples user input devices, such as a keyboard 2028, a pointing device 2026, and a microphone (not shown) and/or output devices, such as speaker(s) (not shown) to computer system 2000. The display driver 2016 and display adapter 2018 are driven by the CPU 2002 to control the display on the display device 2020, for example, to display information pertaining to a target area under analysis, such as displaying a generated representation of the computational mesh, the reservoir, or the target area, according to certain embodiments.
The present techniques are not limited to the architecture of the computer system 2000 shown in
One or more of the following embodiments in the following paragraphs may be utilized with the processes, apparatus, and systems, provided above, to prepare a model and/or be utilized to produce hydrocarbons:
A1. A method for predicting strength anisotropy in subsurface formations along a wellpath comprising: receiving a stress model of the subsurface formation in the area of the wellpath; obtaining bulk strength properties of target rock; developing anisotropic failure criteria based on the bulk strength properties; combining the stress model with the anisotropic failure criteria to develop an anisotropic wellbore stability model; and using the anisotropic wellbore stability model to enhance hydrocarbon recovery.
A2. The method of paragraph A1, the obtaining bulk strength of target rock is determined from the analysis of well logs.
A3. The method of paragraph A1 or A2, wherein the enhance hydrocarbon recovery comprises using the determined anisotropic failure criteria to model wellbore stability along the wellpath.
A4. The method recited in any of paragraphs A1 to A3, wherein determining the anisotropic failure criteria is determined without coring the subsurface strata of interest for the target rock and without performing laboratory strength measurements of the target rock.
A5. The method recited in any of paragraphs A1 to A4, wherein the developing anisotropic failure criteria based on the bulk strength properties comprises using a pre-established database of strength measurements derived from orientated core testing to develop predictive algorithms relating parameters describing the degree of strength anisotropy associated with diverse geologic planes of weakness to bulk strength properties.
B1. A method for predicting strength anisotropy in subsurface formations comprising: compiling a database of core-based strength measurements over a range of orientations and confining pressures for a plurality of fine-grained rock formations; utilizing an automated surface fitting procedure to determine a common anisotropic failure criterion for all database rocks; developing predictive algorithms relating anisotropy parameters in failure criterion and equivalent bulk strength properties; generating subsurface bulk strength properties in strata of interest; utilizing the predictive algorithms and subsurface bulk strength properties to calculate equivalent anisotropic strength properties in subsurface strata of interest; utilizing predicted anisotropic rock strength for subsurface strata of interest in numerical models of the subsurface.
B2. The method of paragraph B1 further comprising: determining observed anisotropic strength values based on the database of core-based strength measurements; determining predicted strength values based on the predictive algorithms; and using statistical tools to evaluate a goodness of fit between the observed anisotropic strength values and the predicted strength values.
B3. The method of paragraph B2 further comprising modifying the predictive algorithms based on the goodness of fit.
B4. The method of any one of paragraphs B3 to B4 further comprising modifying the database of core-based strength measurements based on the goodness of fit.
C1. A system associated with the drilling of a wellbore comprising: a processor; memory coupled to the processor; an application stored in the memory and configured to: compile a database of core-based strength measurements over a range of orientations and confining pressures for a plurality of fine-grained rock formations; utilize an automated surface fitting procedure to determine a common anisotropic failure criterion for all database rocks; develop predictive algorithms relating anisotropy parameters in failure criterion and equivalent bulk strength properties; generate subsurface bulk strength properties in strata of interest; utilize the predictive algorithms and subsurface bulk strength properties to calculate equivalent anisotropic strength properties in subsurface strata of interest; utilize predicted anisotropic rock strength for subsurface strata of interest in numerical models of the subsurface to develop an anisotropic wellbore stability model; and a monitor coupled to the processor and configured to present aspects of the anisotropic wellbore stability model.
C2. The system of paragraph C1 wherein the application is further configured to: determine observed anisotropic strength values based on the database of core-based strength measurements; determine predicted strength values based on the predictive algorithms; and use statistical tools to evaluate a goodness of fit between the observed anisotropic strength values and the predicted strength values.
C3. The system of paragraph C2 wherein the application is further configured to modify the predictive algorithms based on the goodness of fit.
C4. The system of paragraph C2 or C3 wherein the application is further configured to modify the database of core-based strength measurements based on the goodness of fit.
D1. A method for predicting strength anisotropy in subsurface formations along a wellpath comprising: obtaining relationships between planes of weakness and bulk material strength for a plurality of rock samples; obtaining bulk strength of target rock (e.g., well logs, pre-established techniques, etc.); calculating the planes of weakness for the target rock (strength anisotropy); and using the results to enhance models, which may be used for wellbore stability analysis and/or enhancing recovery.
E1. A method for predicting strength anisotropy in subsurface formations along a wellpath comprising: compiling a database of orientated triaxially compressive strength measurements for a plurality of fine-grained rock formations; utilizing an automated surface fitting procedure to generate a common 3D anisotropic failure criterion for each fine-grained rock formation; developing predictive algorithms relating the anisotropic strength parameters describing the 3D anisotropic failure criterion for each fine-grained rock formation and their equivalent bulk strength properties; statistical tools may be utilized to evaluate the goodness of fit between the observed and predicted strength values; in the event the statistical tool determines that no good fit is shown, the database of rock measurements and/or the predictive algorithm may be adjusted; in the event the statistical tool determines that a good fit is established, the process continues by generating bulk strength properties in the subsurface strata of interest using a variety of established petrophysical techniques; utilizing the predictive algorithms pre-established from the database and the bulk strength properties in the subsurface strata of interest to determine the equivalent anisotropic strength properties in the subsurface strata of interest; utilizing the predicted anisotropic rock strength properties of the subsurface strata of interest in numerical models of the subsurface, including analytical and computational methods for optimizing mud weight design in ERD wells.
F1. The method of any one of paragraph A1 to A5, B1 to B4, D1 or E1, wherein one or more of the steps are stored as instructions in memory and are executed by the processor accessing the memory to execute the instructions.
It should be understood that the preceding is merely a detailed description of specific embodiments of the invention and that numerous changes, modifications, and alternatives to the disclosed embodiments can be made in accordance with the disclosure here without departing from the scope of the invention. The preceding description, therefore, is not meant to limit the scope of the invention. Rather, the scope of the invention is to be determined only by the appended claims and their equivalents. It is also contemplated that structures and features embodied in the present examples can be altered, rearranged, substituted, deleted, duplicated, combined, or added to each other. The articles “the”, “a” and “an” are not necessarily limited to mean only one, but rather are inclusive and open ended so as to include, optionally, multiple such elements.
This application claims the benefit of U.S. Provisional Patent Application 61/777,977, filed Mar. 12, 2013, entitled PETROPHYSICAL METHOD FOR PREDICTING SHEAR STRENGTH ANISOTROPY IN FINE-GRAINED ROCK FORMATIONS, and 61/663,368, filed Jun. 22, 2012, entitled PETROPHYSICAL METHOD FOR PREDICTING SHEAR STRENGTH ANISOTROPY IN FINE-GRAINED ROCK FORMATIONS, the entirety of each are incorporated by reference herein.
Number | Date | Country | |
---|---|---|---|
61777977 | Mar 2013 | US | |
61663368 | Jun 2012 | US |