Click to view article in PDF format.
Structural Analysis of a Core on Fractured Carbonate Reservoir, Brazil: Implications for Exploration and Reservoir Modeling*
Nolan M. Dehler1, Délzio Machado Jr.1, Raphael Hatushika1, Leonardo Gomes1, Barbara Trzaskos1
Search and Discovery Article #40408 (2009)
Posted May 4, 2009
*Adapted from poster presentation at AAPG International Conference and Exhibition, Cape Town, South Africa, October 26-29, 2008.
1 Department of Exploration and Production, Petrobras, Rio de Janeiro, Brazil ([email protected])
Fractured Albian carbonate reservoirs, localized along eastern Brazilian marginal basins, have proven to contain important hydrocarbon reserves. The purpose of this paper is to discuss the structural analysis carried out on a selected oriented core, and its application for exploring and modeling fractured reservoirs. Attributes of fractures and faults measured on the core were correlated with depositional sequences, facies, and other petrophysical properties. We observe that the most deformed facies are mechanically soft, porous laminated mudstones. We suggest that this occurs due to flexural slip. We suggest that a clear, positive correlation between
fracture
density and some petrophysical parameters occurs when deformation results from wider wavelength structures in relation to thickness, like large folds.
Nevertheless, this correlation is not so obvious when restricted fault-related deformation occurs. Like sub-seismic faults, flexural slip zones are virtually undetectable by seismic imaging, and may connect different structural levels in the reservoir. The detachment horizons may also represent a decouple between distinct mechanical layers, and hence a boundary between zones with distinct stress axis orientation. It could be expected that in this case, the
fracture
systems could have distinct hydraulic properties from the decouple surface.
|
|
Carbonate hydrocarbon reservoirs play an important role in storage and production of worldwide reserves. These reservoirs are frequently fractured, due to the more fragile rheology of carbonate rocks during deformation in upper crustal conditions. These rocks are also affected earlier by diagenetic events (e.g. recurrent dissolution and cementation), which can significantly change the hydraulic and mechanical characteristics of the system with time (Olson et al., 2007; see also Hesthammer et al., 2002). Hence the distribution of fractures in these rocks, although influenced by other factors such as composition, grain size, layer thickness, strain rate, structural position, depends upon the kind of chemical transformations they have suffered throughout their diagenetic history, in relation to the timing of deformation (see Nelson, 2001). The relative higher chemical mobility of these rocks leads to an additional complexity to the usual structural and geomechanical modeling, as it considers the diagenesis as a controlling factor of hydraulic characteristics of the system (Laubach, 2003).
The distribution of fractures and faults on rocks is heterogeneous and interdependent at a certain extent (Price and Cosgrove, 1991; Shipton and Cowie, 2003; Ortega et al., 2006; Micarelli et al., 2006) because faults control a fraction of the fracturing amount in most areas. The detection of faults is also critical in modeling the hydraulic behavior of
The core presents many types of structures such as joints and shear fractures (faults), brecciated zones, tension gashes and stylolites. Most
The microstructural analysis evaluated the relations between the tectonic features, diagenesis and the migration of hydrocarbon, and determined the existence of a main cementing event followed by a partial dissolution and porosity production. Most of
The structural data obtained from the core analysis were presented as curves of Between 4049 m and 4100 m there are two grainstone intervals separated by a layers of packstone and wackestone. An interval between 4053.5 m and 4063.9 m was lost, probably due to the presence of very porous oolitical grainstones, similar to the beds containing oil found immediately below (Figure 1). From 4049 m packstone and wackestone intercalations occur again, with a tendency of thinning upwards up to 4035m. From this depth upward, fine grained facies are predominant (mudstones). The mudstones at this level are generally laminated. This facies corresponds to the base of the Middle Albian sequence. The remaining intervals below the laminated mudstones belong to the Lower Albian sequence (Figure 1).
The grainstones seem to contain fewer fractures than the other facies. This control is evident in the lower portion of the core below the lost interval. Most fractures are concentrated in the interval between 4070 m and 4080 m, composed of packstone, wackestone and mudstone intercalations. In the upper portion, a prevailing wackestone and mudstone interval was observed, which contains the highest concentration of fractures in the entire section. The whole interval is fractured with average density around 2 fr/m, except to the upper 5 meters, which do not contain fractures. The
Figure 1 shows that the fracturing tends to be more intense in the fine-grained facies, especially in laminated ones. However, some intervals do not follow this correlation. Concentrations of fractures in grainstones, observed at 4087 m and secondarily at 4052 m, can not be fully understood in a model in which
The gamma ray log is more influenced by the faciology and can be used to characterize the lithological types, which keep straight relation with the
Structural Model – Implication for
The results obtained with the inversion of fault/slick data, as well as the description of natural structures of the core, suggest a stress partitioning, or vertical stratification of stress in the structure, at least for some rift structures. In the overlapping sequence, relatively more anisotropic, the data from fault inversion suggests the presence of instantaneous maximum horizontal stretching (σ3), toward the SE. The σ2 is also horizontal and trends toward the SW. In the lower sequence, the fault and slick data suggest the maximum instantaneous stretching points to the NNE, while the intermediary axis points ESE, showing a shift between the minimal and intermediary axis of the vertical stress ellipsoid (Figure 2). For both inversions the maximum stress is vertical, showing an extensional
Two models were assumed as possible (considering only the geometrical characteristics), for the production of flexural shear and fracturing associated with a rift structure: the first model is considered as being formed by the generalized sub-horizontal extension, accommodated by normal faults, layer rotation and flexural slip (Ferril et al. 1998 – see also Wernicke and Burchfiel 1982) (Figure 3b). The second model is associated with the arching of the structure, geometrically similar to what occurs in folding due to flexural slip (Ramsay 1967, Wilson and Cosgrove 1982). The main geometrical consequence for the exploration and production is that the full development of this system is placed at the flank of the arched structure, concentrating in the rheological interfaces of the main mechanical units (Stephenson et al., 2007). This may produce more fractured corridors parallel to the bedding surface, which may connect the upper and lower structural portions of the zone of interest, as well as more
The results of the structural analysis are also consistent with the previously proposed model for rift tectonics I the basin (see Figure 4). According to Gangá & Machado Jr. (2007) and Machado Jr. et al. (2007), rifts have been displaced mainly toward the southeast, but northeast-southwest extension was also present, as shown by the geometry of normal faulting (Figure 4). Hence, the obtained inversion of faults and slick data are consistent with the finite geometry pattern of normal faulting.
The presented data suggest that the parameters of porosity composition, grain size, and thickness of layers are important in the modelling of T In the studied case, the existence of faults associated with the movement on bedding planes are virtually untraceable by seismic imaging independently of its throw. These structures can connect different levels of the reservoir if the relation between the fracturing and the diagenetic modifications is favorable. The possibility of connection seems low due to the intense cementation. The presence of decoupling surfaces between the main mechanical units may favor a shift in the stress axis along a section normal to the surface. In these cases, the orientation of fractures which are more favorable to fluid flow may vary abruptly from this surface. We would like to acknowledge geologist Fernanda Mourão de Brito from Petrobras research center (CENPES) for providing sedimentologic analysis and laboratory support, and also to Petrobras for permission to publish this paper. Ferril, D.A., A.P. Morris, S.M. Jones, A. Stamatokos, 1998, Extensional layer-parallel shear and normal faulting: Journal of Structural Geology, v. 20, p. 355-362. Ganga, A.O., and D.L. Machado, Jr., 2007, Padrões de fraturamento em carbonatos com base em análise de curvatura de horizontes, in VI Simpósio de Interpretação Exploratória, Vitória, Espírito Santo. Laubach, S.E., 2003, Practical approaches to identifying sealed and open fractures: AAPG Bulletin, v. 87, p. 561-579. Machado, D.L. Jr., M.V. Sant’Anna, M. Basseto, 2007, Análise de fraturas para a estrutura Jabuti, ring fence de Marlim Leste – Bacia de Campos, in VI Simpósio de Interpretação Exploratória, Vitória, Espírito Santo. Micarelli, L., A. Benedicto, C.A.J. Wibberley, 2006, Structural evolution and permeability of normal fault zones in highly porous carbonates: Journal of Structural Geology, v. 28, p. 1214-1227. Nelson, R.N., 2001, Geologic Analysis of Naturally Fractured Reservoirs: Gulf Professional Publishing, 2nd edition, 332 p.
Olson, J.E., S.E. Laubach, R.H. Lander, 2007, Combining diagenesis and mechanics to quantify
Ortega, O.J., R.A. Marrett, S.E. Laubach, 2006, A scale-independent approach to Price, N.J., and J.W. Cosgrove, 1991, Structural Analysis of Geological Structures, Cambridge University Press, 502 p. Ramsay, J., 1967, Folding and Fracturing of Rocks, McGraw-Hill, New York, 568 p. Shipton, Z.K., and P.A. Cowie, 2003, A conceptual model for the origin of fault damage zone structures in high-porosity sandstone: Journal of Structural Geology, v. 25, p. 333-344. Stephenson, B.J., A. Koopman, H. Hillgartener, H. McQuillan, S. Bourne, J.J. Noad, K. Rawsley, 2007, Structural and Stratigraphic controls on fold-related fracturing in the Zagros Mountains, Iran: Implications for reservoir development, in L. Lonergan , R.J.H. Jolly, K. Rawsley, and D.J. Sanderson (Eds.), Fractured Reservoirs: Geological Society of London, Special Publication No. 270, p. 1-22. Wernicke, B., and B.C. Burchfiel, 1982, Modes of extensional tectonics: Journal of Structural Geology, v. 4, p. 105–115. Wilson, G., J. Cosgrove, 1982, An introduction to small-scale geological structures, George Allen and Unwin, London, 128 p.
|
