Abstract
Upfaulted ridges of Neoarchean crystalline basement rocks formed in the Faeroe-Shetland basin as a consequence of Mesozoic rift processes and are an active target for oil exploration. We carried out a comprehensive fault and fracture attribute study on the extensive exposures of geologically equivalent crystalline basement rocks onshore in NW Scotland (Lewisian Gneiss Complex) as an analogue for the offshore oil and gas reservoirs of the uplifted Rona Ridge basement high. Our analysis shows a power-law distribution for fracture sizes (aperture and length), with random to clustered spacing and high connectivity indices. Regional variations between the Scottish mainland and the Outer Hebrides are recognized that compare directly with variations observed along the Rona Ridge in the Faeroe-Shetland basin. Here we develop a model for the scaling properties of the fracture systems in which variations in the aperture attributes are a function of the depth of erosion beneath the top basement unconformity. More generally, the combination of size, spatial and connectivity attributes we found in these basement highs demonstrates that they can form highly effective, well-plumbed reservoir systems in their own right.
Supplementary material: Additional methods and results are available at: https://doi.org/10.6084/m9.figshare.c.5017139
Thematic collection: This article is part of the The Geology of Fractured Reservoirs collection available at: https://www.lyellcollection.org/cc/the-geology-of-fractured-reservoirs
The metamorphic basement rocks of the Lewisian Gneiss Complex may once have seemed an unlikely target for hydrocarbons, but a series of recent discoveries means that they are now a focus for exploration activity in the Faeroe-Shetland basin (Fig. 1). The delineation of the Clair and Lancaster fields, and associated prospects, confirms that there are significant oil accumulations in Neoarchean basement lithologies of similar age to the onshore Lewisian Gneiss Complex in NW Scotland. These crystalline basement ridges were uplifted and exposed at surface during Mesozoic rifting before being buried again during the Cenozoic Atlantic margin opening (Stoker et al. 2018). Given that permeability in basement reservoirs is predominantly fracture-controlled (e.g. Achtziger-Zupančič et al. 2017) and given the general uncertainty associated with fractured reservoirs systems (Nelson 1985), a renewed interest in studying analogue basement-hosted fracture systems is unsurprising.
Map of the NW UK continental shelf showing location of fields, prospects, top basement depth map offshore and onshore crystalline basement exposures.
Well-exposed outcrops of the Lewisian Gneiss Complex occur in the mainland of NW Scotland and in the Outer Hebrides, the latter being an elongate uplifted crustal block with similar dimensions to the Rona Ridge offshore, where significant hydrocarbon discoveries have been made (Fig. 1).
Although rare, producing basement reservoirs in a range of fractured igneous and metamorphic host rocks are known from 27 countries worldwide (Gutmanis 2009, Gutmanis et al. 2015). They form by conventional means with migration from a mature basinal source rock into a fractured reservoir trap and are contained by a low permeability top seal. Oil accumulation in the crystalline basement of the Clair Field has long been known (e.g. Coney et al. 1993), but recent discoveries in other parts of the Rona Ridge, where basement has been specifically targeted, include the Lancaster Field, and the Lincoln, Halifax and Whirlwind prospects (Slightam 2012; Trice 2014).
The Lewisian Gneiss Complex of NW Scotland has, over its c. 3.2 Ga history, formed part of an active accretionary margin, a collisional foreland, a rifted margin (at least twice) and most recently a passive margin, and therefore retains a record of several generations of both ductile and brittle deformation, metamorphism and fluid-flow events. This complex history has produced a highly heterogeneous array of lithologies, metamorphic grades and structural styles (e.g. Park 1970).
Here, we present an analysis of fracture attribute datasets collected from brittle structures exposed across the onshore Lewisian Complex. The comprehensive nature of the data compilation (some 100 individual datasets) enables us to identify correlations between the mainland Lewisian and the Clair basement, and the Hebrides exposures with the Lancaster Field. We then propose a simple model that accounts for the first-order differences in fracture attributes and their scaling that is linked to recent work on the geological nature and development of the fracture systems and their infills (Holdsworth et al. 2019, 2020a; Trice et al. 2019). This work has led to a new understanding of the significance of fissuring processes in enhancing the capability of uplifted rift blocks of fractured crystalline basement to host significant accumulations of hydrocarbons and also provides a general model for explaining fluid flow in other uplifted basement lithologies in similar settings below regional unconformities.
Geological setting
Location and regional structure
The Precambrian rocks of the Lewisian Gneiss Complex of NW Scotland form a fragment of the continental basement of Laurentia that was isolated from North America by the opening of the North Atlantic (Bridgwater et al. 1973). The rocks comprise trondjemitic, tonalitic and granodioritic orthogneisses, with subordinate units of metabasic-ultrabasic and granitic composition, together with local units of metasedimentary rock. The complex then underwent a long history of major, crustal-scale geological events during the Archean and Paleoproterozoic (see Wheeler et al. 2010 and references therein) and is divided into a number of tectonic regions or ‘terranes’ that are separated by mainly steeply dipping shear zones or faults.
Two different tectonic views exist concerning the early geological evolution of the Lewisian Gneiss Complex. The first, based on the classic geological mapping by Sutton and Watson (1951), suggests that much of the basement gneiss is a single piece of continental crust that shares a common early history. This model was rooted in the recognition of two fundamentally separate groups of tectonothermal events, one predating and one post-dating the intrusion of a regional swarm of NW–SE-trending mafic to ultramafic dykes known as the Scourie Dyke swarm (Sutton and Watson 1951). Areas where evidence for these early events is not preserved were thought to have undergone intense overprinting and reworking during later Paleoproterozoic events (‘Laxfordian’). A more recent alternative hypothesis, proposed by Friend and Kinny (2001) and Kinny et al. (2005), is founded in zircon geochronology and suggests that each terrane has different Archean age spectra. They view the Lewisian as a collage of lithologically and geochronologically distinct tectonic units or terranes bounded by regional shear zones that were assembled progressively during a series of Precambrian amalgamation episodes.
Neoarchean orthogneisses of broadly similar composition and age extend north of the Scottish mainland at least as far as the northernmost tip of Shetland (Holdsworth et al. 2018; Kinny et al. 2019). Equivalent units underlie much of the Faroe-Shetland basin and the c. 200 km Rona Ridge, as shown by analyses of basement-penetrating offshore cores (Fig. 1; see Ritchie et al. 2011). These basement rocks have protoliths and early amphibolite facies tectonothermal events of broadly the same age as those of the Lewisian gneisses (c. 2.8–2.7 Ga), but lack the Paleoproterozoic (Laxfordian) overprinting events (Holdsworth et al. 2018). These rocks are directly comparable to those of the North Atlantic Craton in Eastern Greenland and Canada, while the reworked rocks of the Lewisian Gneiss Complex in NW Scotland and the Hebrides are thought to be southeasterly equivalents of the Nagssugtoqidian gneisses of Eastern Greenland (Mason and Brewer 2004; Holdsworth et al. 2018).
Early metamorphic assemblages and structures, together with the Scourie dykes, are heterogeneously overprinted by Laxfordian reworking in parts of the Lewisian Gneiss Complex. These older features are only clearly preserved in certain areas of the mainland complex, most notably the ‘Central Region’ or Assynt Terrane (Fig. 2). The main phases of the Laxfordian deformation and metamorphism predominate in the Rhiconich and Gruinard terranes that lie to the north and south of the Assynt Terrane respectively (Fig. 2). The NW–SE strike-slip-dominated shear zones that form the terrane boundaries on the Scottish mainland – and another 1 km-wide structure in the centre of the Assynt Terrane known as the Canisp Shear Zone (Fig. 2) – are thought to have formed and perhaps initially juxtaposed the three terranes during an early (Inverian) event c. 2.4 Ga (Park et al. 2002). All were then reactivated during episodic Laxfordian shearing (c. 1.9–1.66 Ga), often with alternating shear senses (Park et al. 2002). The predominantly amphibolite facies granodioritic orthogneisses of the Outer Hebrides preserve a superficially similar relative chronology of structures and metamorphic assemblages as on the mainland.
Lineament interpretation for well-exposed parts of the mainland and Hebrides basement of NW Scotland. Outcrop fracture sample sites are labelled and shown in blue (Hebrides), light green (mainland – Rhiconich Terrane) and dark green (Assynt Terrane). Summary rose diagrams of fracture orientations for the mainland and Hebrides. Inset map shows Clair Field (outline of Clair first development phase in black line) with lineaments from Pless (2012). Underlying onshore geology from BGS 1:625,000 geology map. Main units include Neoarchaean with/without Paleoproterozoic (Laxfordian) overprint: A = intermediate to granitic gneiss (Lewisian), Paleoproterozoic: Z = felsic intrusive rocks, Zm = Mafic intrusive rocks, Zs = metasedimentary rocks, M = Moine metasediments, Mesoproterozoic: S = Stoer Gp, Neoproterozoic: T = Torridonian, CO = Cambro-Ordovician sedimentary rocks, OS = Ordovcian/Silurian alkaline syenite, F = fault rocks (mylonites, cataclasites and pseudotachylytes), PT = Permo-Triassic sedimentary rocks. Major structures are labelled – KLB F = Kinlochbervie Fault.
Faulting and fracturing history
Currently exposed levels of the Lewisian Gneiss Complex passed through the brittle-ductile transition at some point after c. 1.66 Ga and were close to the surface by c. 1.2 Ga, the depositional age of the unconformably overlying, unmetamorphosed Stoer Group on the Scottish mainland (Beacom et al. 2001; Holdsworth et al. 2020b). Unsurprisingly for rocks that preserve a record of brittle deformation processes that occurred across a range of crustal depths, the Lewisian displays a wide range of micro- to regional-scale brittle fractures. A broad spectrum of types is developed that are difficult to strictly separate using an arbitrary classification scheme (Pless 2012; Franklin 2013). These include the following:
Joints are predominantly Mode 1-type tensile fractures based on a general lack of observed offsets of pre-existing features such as compositional banding. They are typically closed and only become open due to the effects of weathering – either in the geological past or present day – or due to later tectonic processes (such as fissuring – see below). They occur on a variety of scales but, like many Mode 1 fractures developed in crystalline basement rocks worldwide (e.g. Wang et al. 2019), are commonly of large lateral extent both horizontally and vertically (Fig. 3a).
Typical basement fracture types and fills. (a) Closely spaced laterally and vertically extensive jointing in granitic gneiss Lewisian basement, Uyea, Shetland (see Kinny et al. 2019). Note that later Devonian-age dykes have exploited these well-developed joint systems. (b) Composite carbonate veins cutting mafic gneisses, Traigh Dhail Mor, Isle of Lewis (1.2 km SW of Dail Beag, see Fig. 2). Note large open vug (V). (c) Cross-sectional view of part of a c. 30 m-wide fissure filled with chaotic millimetre- to metre-sized angular clasts of basement, and possible red sediment locally cemented by carbonate. Age of fill uncertain, but note that the contact with the wall rock has been exploited by a Cenozoic basalt dyke, suggesting that the breccia is likely Mesozoic in age. Traigh Dhail Mor, Isle of Lewis. (d) Close-up view of crudely laminated nature of the fill at Traigh Dhail Mhor suggesting an element of water-lain deposition. (e) Fissure filled with chaotic collapse breccia where the matrix is cataclasite and pseudotachylyte, Canisp Shear Zone, Achmelvich. Note that in this case, the development of the dilational cavity is thought to be related to seismogenic slip events along the well-developed foliation in the wall rocks at depths >5 km (see Hardman 2019 for details). (f) Foliated multicoloured gouges and breccias from the core of the Seaforth Fault, a major N–S Mesozoic normal fault with kilometre-scale offsets that cuts the Isle of Lewis (Fig. 2; see Franklin 2013 for details).
Veins are dominantly millimetre- to metre-scale tensile or hybrid fractures filled with a variety of hydrothermal minerals including (commonly) quartz, epidote, carbonates (calcite, siderite), chlorite, K-feldspar (adularia), iron oxides and (less commonly) base metal sulphides, prehnite and a variety of zeolites (Fig. 3b). The majority are entirely occluded by their mineral fills, but in some cases, partial fills and vuggy textures are preserved. Like many basement terrains, veins completely filled with dark, aphanitic pseudotachylyte (friction melts) are well developed locally and are typically associated with fault zones formed relatively early in the brittle deformation history (e.g. Imber et al. 2001; Holdsworth et al. 2020b).
Fissures are millimetre- to decametre-scale dilational (predominantly Mode 1) fractures filled or partially filled with often complex, composite fills formed at a range of crustal depths. Many formed close to the surface in the geological past and are spatially associated with regional unconformities at the base of the Torridonian or Mesozoic cover sequences (e.g. Beacom et al. 1999; Jonk et al. 2004). Fills here include wall rock collapse breccia, hydrothermal minerals and fine-grained sediment, sometimes with a laminated structure and cement consistent with having been deposited by flowing water in subterranean open cavity systems (Fig. 3c and d). Deeper fissure fills include magma, i.e. Paleozoic to Cenozoic dykes, and wall rock collapse breccias mixed with friction melt and hydrothermal minerals (Fig. 3c and e).
Shear fractures range from simple ‘clean break’ brittle faults with sub-millimetre-scale offsets through to large complex fault zones with kilometre-scale offsets. Fault rocks typically begin to appear once displacement exceeds more than a few millimetres, and include early formed pseudotachylytes and cataclasites, breccias and gouges; all with associated hydrothermal mineral assemblages similar to those seen in associated vein systems. Fault rocks formed earlier in the brittle deformation history are generally cohesive and highly indurated while those formed later and nearer to the surface are typically incohesive and easily weathered. Polished fault surfaces with slickenlines or hydrothermal mineral slickenfibres – particularly of quartz, epidote, chlorite or carbonate – are widely preserved (Fig. 3f). Large-scale fault zones – such as the Seaforth Fault in Lewis (Fig. 2; Franklin 2013) are typified by the development of well-defined cores with foliated gouges (Fig. 3g) and broad, chaotically fractured damage zones (e.g. Pless et al. 2015). Some – but not all – show evidence of reactivation (e.g. Imber et al. 2001; Holdsworth et al. 2020b).
A long history of fracturing is recognized on the Scottish mainland with at least three main fault/fracture sets preserved in the foreland region west of the Caledonian Moine Thrust Zone (Fig. 2). Each is associated with different fault geometries, kinematics and fault rock assemblages. These are (from earliest to latest):
NW–SE ‘Assyntian’ or ‘Late Laxfordian’ sinistral fault arrays (Holdsworth et al. 2020b) which are most abundant as reactivation events in pre-existing NW–SE Laxfordian shear zones (e.g. Canisp Shear Zone, Fig. 2) and along the margins of pre-existing Scourie Dykes (Beacom et al. 2001; Pless 2012). These structures are associated with the development of cohesive cataclasites and pseudotachylytes. Their Mesoproterozoic (c. 1.55 Ga) age is constrained by Re–Os dating of associated copper sulphide mineralization in the Assynt Terrane (Holdsworth et al. 2020b) and they demonstrably predate deposition of the unconformably overlying Stoer Group c. 1.2 Ga. A generally north–NE-trending set of complex polymodal fracture arrays was thought by Beacom et al. (1999, 2001) to be associated with the Stoer Group age rifting, but more recent fieldwork and thin section analysis (Hardman 2019) have shown that these fractures are synchronous with the NW–SE structures. The NE–SW structures are commonly associated with dilation and collapse brecciation, pseudotachylyte injection and epidote mineralization that predates Stoer Group deposition (Holdsworth et al. 2020b).
Post-Torridonan (c. 1.04 Ga) faults, including isolated thrusts and strike-slip faults related to the Paleozoic Moine Thrust Zone; many of the NW–SE Late Laxfordian faults also show evidence of reactivation close to the thrust belt (Krabbendam and Leslie 2010; Pless 2012). Most are clean breaks or are associated with the development of well-cemented breccia and gouge. Multiple sets of microfractures and fills of mainly Paleozoic age may also be present (e.g. Laubach and Diaz-Tushmann 2009; Ellis et al. 2012).
Mesozoic age structures that are generally NE–SW and NW–SE-trending dip-slip and strike-slip fracture sets that are widely associated with incohesive gouges and carbonate mineralization (Laubach and Marshak 1987). Many of the fissure structures are thought to have formed during Mesozoic rifting events when the basement was at or close to surface. These structures are likely more widespread than has generally been assumed and they typically show little evidence for reactivation except along major faults (e.g. Coigach Fault; Roberts and Holdsworth 1999). Holford et al. (2010) showed that NW Scotland has experienced multiple episodes of Mesozoic and Cenozoic burial and exhumation associated with passive margin formation, rifting processes and inversion.
The early brittle faulting history of the Outer Hebrides is dominated by the development of the SE-dipping Outer Hebrides Fault Zone (OHFZ) that initially developed as a mylonitic shear zone, possibly of Laxfordian or Grenvillian age (Imber et al. 2001, 2002) (Fig. 2). It then experienced a series of reactivation events from the Neoproterozoic to the Mesozoic, but direct geological or radioisotopic evidence for the age of movements is sparse. The likely presence of Torridonian rocks within Minch Basin (Fig. 2) suggests that the OHFZ may have been active as a normal fault c. 1.04 Ga, but there is no clear record of this faulting yet recognized in outcrops. Onshore, post-mylonite deformation along the OHFZ was initially brittle and was associated with the development of pseudotachylyte-bearing fault veins and thick, SE–east-dipping pseudotachylyte-ultracataclasite crush zones all along the eastern margin of the Hebridean island chain (e.g. Sibson 1977, Fig. 2). The brittle faults and crush zones are overprinted by a network of macroscopically ductile, greenschist facies phyllonitic shear zones that developed due to the influx of hydrothermal fluids during top-to-the-NE sinistral strike-slip shearing along the OHFZ (Butler et al. 1995; Imber et al. 2001). These shear zones were themselves then reactivated during late Caledonian brittle-ductile top-to-the-east extensional deformation (Imber et al. 2001). The Permo-Triassic Stornoway Formation (Steel and Wilson 1975) was deposited in eastwardly prograding alluvial fans associated with normal fault scarps developed in hanging wall of the OHFZ (Fig. 2). The sedimentary sequence contains clasts of basement gneisses and OHFZ-derived fault rocks, suggesting that the northern Outer Hebrides was exhumed by the earliest Mesozoic era. The rocks of both the Lewisian Gneiss Complex (including the OHFZ) and the Stornoway Formation are cut by E–W, NW–SE and NE–SW fractures, some of which – together with generally NNW-trending Tertiary dykes, form prominent topographic lineaments (e.g. Loch Seaforth) (Fig. 2) (Franklin 2013). Faulting events are widely associated with the development of generally incohesive gouge, breccia and fissure fills with local widths of at least 30 m, but perhaps up to 100 m, together with extensive carbonate mineralization. The Tertiary dykes mostly cross-cut the faults and fault rocks which are assumed therefore to be of Mesozoic age, but some fault sets show evidence of significant post-dyke reactivation during the Cenozoic, notably a prominent set of east–west-trending structures which are also characterized by a later phase of milky carbonate mineralization (Franklin 2013).
Methodologies
Sampling of fault and fracture networks
The datasets reported in this study were mainly acquired using the 1D linear scanline method (Priest and Hudson 1981; Baecher 1983; McCaffrey and Johnston 1996; Ortega et al. 2006). This method allows a relatively simple characterization, albeit with known biases, of fracture sizes and intensities, and generally can be deployed at most field localities. The data (observations along a sample line) are closely analogous to logs from borehole or drill core taken from a prospect or reservoir. To gain 2D (map) information on the spatial and topological relationships within the fractured system, we also conducted 2D window sampling (Odling 1992; Mauldon et al. 2001; Rohrbaugh et al. 2002; Zeeb et al. 2013; Watkins et al. 2015; Sanderson and Nixon 2015), enabling access to connectivity estimates for the fracture array, which are a key input for modelling fluid flow.
For the linear scanlines, fracture orientations, lengths and apertures, together with composition and texture of fracture infills and fracture terminations on joints and other faults were recorded at measured intervals along the sample line. The start and end point of each transect was recorded using a hand-held GPS unit. Most of the fractures are filled, or partially filled, with minerals (mainly quartz, epidote or calcite) and, following Laubach (2003) and Ortega et al. (2006), the apertures measured in this study are the opening displacement where the scanline intersects the fracture including any fill, i.e. the ‘kinematic aperture’. This is equivalent to the fracture thickness of McCaffrey and Johnston (1996) and Massiot et al. (2015).
Fracture samples
The 1D datasets were collected in the field mainly from natural exposures of the Lewisian Gneiss Complex (well exposed in coastal settings), but also from road cuttings where natural fractures may be easily distinguished from those created by blasting (for detailed descriptions of the sample locations, see Pless 2012 and Franklin 2013). The locations of study sites across the mainland Scotland and Hebrides are shown in Figure 2 with full details of individual sample lines given in the supplementary tables. The initial studies focused on size (aperture, length), spatial characterization (orientation and spacing) and the topological characteristics of the fracture systems. Our database contains more than 100 individual datasets (48 aperture and 29 length samples and 27 topological estimates) chosen because they capture the fracture systems that formed from Proterozoic to Cenozoic times (see above). Details of the fracture samples, including location, host lithology, number of fractures, sample line length for 1D samples, area for 2D samples and types of structure intersected are given in the Supplementary file (Tables S1 and S2). To extend the analysis to other scales, the above-mentioned scanline methods were adapted and applied to aerial photographs and optical data (BGS NextMap data) to quantify fracture lengths in 1D (see lineaments shown on Fig. 2). These datasets were collected before we had fully appreciated the importance and extent of fissure formation in the basement, particularly in the Hebrides, nonetheless we think the study provides important baseline information.
Data from the Clair Field comprise fractures logged in wells 206/7a-2 and 206/8-8 that were drilled by Elf into crystalline basement gneisses of the Clair ridge (see Holdsworth et al. 2018; Fig. 2). Core samples were examined at the Iron Mountain core storage facility, Aberdeen, and a fracture analysis was conducted by Pless (2012) and in this study. The basement core slab samples from 206/7a-2 are in 10 m lengths at irregular intervals from measured depths of 2140 m to 2600 m (see Holdsworth et al. 2018 and S1).
At regional scales, a fracture interpretation of Clair 3D top basement seismic attribute maps was performed (see Pless 2012 and Fig. 2 inset). From this fracture map we were able to derive fracture length distributions along 1D sample lines across the maps and in 2D windows. An equivalent study of fracture lengths in 1D and 2D was carried out on the onshore lineament maps from the mainland and Hebrides (see Fig. 2). The lineament maps show density variations related to the amount of younger cover rocks or Quaternary material (Fig. 2) and so our length analyses were conducted on lines that cross, or windows that sample, regions with high density and thin cover. We carefully filtered the datasets to make sure that those with low numbers (c. n < 40) were omitted. We also checked that the datasets were collected and formatted in a comparable way as they have been assembled from a number of studies.
For the topology study, photographs from outcrops in the Assynt Terrane, Clair core 206/7a-2 supplemented by samples from the 205/21-1A from further along the Rona Ridge (Lancaster Field) collected at BGS core store form the basis for picking of nodes and branches. All node and branch picking was carried out manually to ensure the correct network topology was recorded.
Data analysis
In this study we assessed the distribution of fracture size (aperture and length) attributes, collected from the 1D sample lines as a primary characterization of the brittle deformation within the basement Lewisian Gneiss Complex. We collected fracture data from drill core samples from the basement of the Clair Field for comparative purposes. We only report the aperture data here as the fracture lengths are heavily censored by the dimensions of the drill core. However, we report regional-scale 2D length data from onshore Lewisian terranes from the lineament dataset derived from the optical data that is equivalent in scale to the offshore seismic attributes maps. This allows us to constrain further the upscaling of fracture attribute and compare the offshore basement fracture mapping to that performed on seismic attribute maps.
Fracture sizes
Fracture intensity plotted as cumulative distribution (population) plots enables an assessment of the distribution, spatial and scaling properties of the fracture population (i.e. the ratio of small to large fractures for a given sample line length). Fracture attribute distributions display three main types of statistical distribution (Gillespie et al. 1993; Bonnet et al. 2001; Zeeb et al. 2013): (1) exponential, random or Poisson distributions are characteristic of a system with a randomized variable; (2) log-normal distributions are generally produced by systems with a characteristic length scale, for example layer-bound jointing (Narr 1991 and Olson 2007); and (3) power-law distributions lack a characteristic length scale in the fracture growth process (Zeeb et al. 2013) (see Supplementary file S1). Although some fracture populations are better described by scale-limited laws, such as log-normal or exponential distributions, it is generally accepted that power-law distributions and fractal geometry provide a widely applied descriptive tool for fracture system characterization (e.g. Bonnet et al. 2001). Ideally, the best-fit power-law distribution should be constrained over several orders of magnitude (Walsh and Watterson 1993; McCaffrey and Johnston 1996). However, in practice this is typically very difficult to achieve at a given scale due to sampling limitations. Fracture sampling issues (e.g. censoring and truncation) are commonly encountered and can result in an incomplete description of the full population. For instance, when large fractures are incompletely sampled in a power-law population, the resulting plot can resemble a log-normal distribution. Following Ortega et al. (2006), Dichiarante et al. (2020) have shown how a multi-scale approach can be used to better constrain the scaling laws for fracture size attributes. As pointed out by Clauset et al. (2009), use of the maximum likelihood estimator (MLE) is preferred over a least square regression analyses (R2) for the fitting of power-law distributions. In this study we used MLE scripts developed by Rizzo et al. (2017) as used in the FracPaQ toolbox (Healy et al. 2017). In addition, we followed the Dichiarante et al. (2020) modification in which the MLE for power-law, exponential and log-normal fits are calculated on systematically truncated and censored datasets to find the optimum distribution parameters (see S1).
Fracture spatial organization
The spatial organization of fracture systems are a property of the orientation and clustering of fractures in 1D sample lines. For many years the Coefficient of Variation (Cv) – the standard deviation of all spaces between adjacent fractures divided by the mean spacing (Gillespie et al. 1993, 1999) – has been used as to describe clustering. These authors showed that a Cv > 1 reflected a clustered distribution and could be expected in non-layered rocks (like basement). A random or Poisson distribution gives an exponential (Cv = 1). Superimposition of power-law distributions can give a ‘Kolmogorov’ distribution (log-normal) (Cv < 1). Log-normal (and normal distributions) are also produced by ‘saturation’ models when fractures are produced in well-bedded sequences (Bai et al. 2000). The Cv values for 1D basement sample lines are reported in this study; however, we note that that there are issues with the sensitivity of this method and that the method does not take into account the spatial arrangement of structures or the scale of clustering (Marrett et al. 2018). The correlation analysis method subsequently developed by Marrett et al. (2018) is now the preferred method for analysing fracture spatial distributions and will be the subject of further work on the datasets collected in this study.
Fracture topology
While the 1D scanline data provide information about fractures as single entities and their distribution, 2D topology analyses consider fractures as part of a network and provide access to fracture connectivity assessment. The 2D analysis used here has been carried out on fracture maps at regional-scale (metre-decametre) DEM images and seismic attribute maps (see Fig. 2). At smaller scales (centimetre-metre) we carried out topological analysis on core samples and outcrops. We followed the methodology of Sanderson and Nixon (2015) in defining nodes and branches. ‘Nodes’ are defined as the point where a fracture terminates (I-type), abuts against/splays from another fracture (Y-type) or intersects (cross-cuts) another fracture (X-type). ‘Branches’ are the portions of a fracture confined between two nodes.
The number of nodes and branches for a given fracture network is strictly related, meaning that by knowing one of the two elements for the fracture network, it is possible to quantify all its components. NI, NY and NX are defined as the number of I-, Y- and X-type nodes and PI, PY and PX their relative proportions. Once the number of nodes and/or branches making up a fracture array are known, the connectivity can be visualized using a ternary plot of the component proportions or can be quantified by calculating the number of connections existing in the 2D map. In general, X-type nodes provide four times and Y-type nodes three times more connectivity than I-type nodes (Sanderson and Nixon 2015). An array dominated by I-nodes is isolated, while arrays dominated by Y- and X-type nodes are increasingly more connected.
Results
Fracture lineaments from the mainland (Assynt and Rhiconich) terranes show strong NE–SW and WNW–ESE trends (Pless 2012; Figure 2). In contrast in the Hebrides, the main lineament trend is NNW–SSE with a subordinate ENE–WSW trend (Franklin 2013; Figure 2). The lineament maps show density variations that particularly reflect the amount of Quaternary cover; for example, see southern and western Lewis compared with the northern region (Fig. 2). At the regional scale, there is no qualitative variation in density of lineaments in relation to major structures such as the OHFZ, Canisp Shear Zone or the Seaforth Fault (Fig. 2). We also see no systematic variation at this scale with the host lithological units (Fig. 2). Pless (2012) have conducted an analysis of fracture density maps which confirms the qualitative observations.
Aperture data
Figure 4 shows cumulative distribution plots for the aperture distributions for localities in Lewisian Complex gneisses on the mainland (20 sample lines), Hebrides (17 lines) and Clair basement core (12 lines). Details of the individual samples and the distribution fitting parameters are given in Table S1. For the mainland, there is high degree of variability, but the data span more than three orders of magnitude from 0.00005 m to 0.5 m (0.05–500 mm) in aperture (Fig. 4a). We note that some constant values appear in the plots at small sizes and are a rounding effect that occurs during the field acquisition. We generally remove repeated values, as recommended by Ortega et al. (2006), but the application of the Terzhagi true thickness correction tends to smear out these clusters of sub-millimetre values towards even smaller values. In terms of the fracture intensity or spacing (y axes), the data show about an order of magnitude spread from low strain (0.05 fractures of 10 mm size per metre) to high strain (1 fracture per metre) (Fig. 4a). For the Hebrides, data span nearly five orders of magnitude from about 0.05 mm to 1000 mm (Fig. 4b). Fracture intensity or spacing (y axes) vary by about an order of magnitude from 0.1 fracture per metre to a higher strain of about 2 per metre for 10 mm aperture fractures. For the Clair core datasets, aperture values range from 0.05 mm to 100 mm and the intensity values are less variable than the onshore datasets ranging from 0.5 to 1.2 per metre for 10 mm aperture fractures.
Fracture aperture intensity data for: (a) mainland Scotland; (b) Hebrides; and (c) Clair basement. The grey polygon highlights the same Fracture Intensity/Aperture space with a slope of −1 and is shown for comparison in each plot. Orange bars show comparative fracture intensity ranges for 10 mm aperture fractures as discussed in text. Data from locations that sample Mesozoic structures on the Hebrides include Garrabost, Memorial Cairn, Pabail, Seisadar and Tolstadh.
Aperture distribution data for all regions can be described by power-law scaling or log-normal distributions with greater than 95% confidence calculated using the MLE method with a slight preference for power-law distributions (Fig. 4 and Table S1). The sample lines (Garrabost, Memorial Cairn, Pabail and Seisadar) identified by Franklin (2013 and Supplementary file S1) are those taken across Mesozoic structures and tend to be those that display the highest absolute aperture values (Fig. 4).
The advantage of plotting many datasets together (Fig. 4) is that general trends emerge above variations displayed by individual samples (see Discussion below). One clear signal that emerges is that the power-law exponent is lower for samples from the Hebrides than for the mainland and the Clair basement. This can be seen qualitatively in Figure 4. For the mainland and Clair data (Fig. 4a), the values tend to lie along the grey-shaded reference area which has boundaries with a slope = −1 on the plot except at the lower and upper ranges where truncation and censoring effects are likely. For the Hebrides, the datasets clearly plot along a shallower slope line compared with the shaded reference area. To test this inference, we performed a significance test of the difference between the MLE power-law scaling exponents (individual slope with >95% confidence fits) for the two regions. An average power-law exponent for each region was calculated and the t test statistics confirm that the mainland (average slope α = 1.23, SD = 0.49) and Hebrides (average α = 0.74, SD = 0.26) conditions; t(37) = 4.15, P = 0.0002 are different. These results show that the Hebrides and mainland fractures show different scaling properties, and this implies that there are fewer small aperture fractures in the Hebrides relative to the largest fractures when compared with those seen in the mainland.
Length distributions
The fracture length distributions for faults and fractures from both onshore and offshore regions are presented as cumulative distribution plots of intensity v. length in Figure 5. Length data at smaller scales from outcrops (c. 0.1–10 m) are plotted alongside line samples across the top Clair basement (c. 0.5–50 km) (Fig. 5a). Details of individual samples and distribution fitting are given in Table S2. Again, most of the samples can be described by power-law or log-normal distributions with greater than 95% confidence with a slight preference for log-normal distributions. A general scaling relationship (power-law) from outcrop to regional scale is suggested (Fig. 5a).
Two measures of fracture length intensity scaling. (a) Fracture lengths intersected in 1D samples plotted on a multi-scale diagram from mainland outcrops and the Clair top-basement seismic attribute map. (b) The intensity of fractures per unit area (m) is shown for 2D length data from window samples taken across mainland, Hebrides and Clair seismic attribute and topographic maps.
The regional-scale 2D length data from both onshore lineament mapping and offshore top basement seismic attribute map (Fig. 2 inset) are shown in Figure 5b with details of distribution fitting given in Table S2. The data show good agreement between the onshore Lewisian and the Clair fields (similar intensity values and slopes) at fracture lengths 1–50 km (Fig. 5b). Below 0.5–1 km, the distributions show truncation effects (inflection points on the curves) that are dependent on the scale at which the fractures have been mapped and the level of exposure (onshore this is c. 500 m and offshore it is c. 1 km).
Spatial organization
The Cv values for the spaces between adjacent fractures for each sample line are shown on Figure 6 plotted against the overall fracture intensity for each of the sample line datasets assembled in this study. Plotting in this way enables us to compare Cv values and assess the spatial organization at different scales. The values show a range of behaviours from more uniform spacing (<1) to more clustered distributions (>1). There is a large amount of variation, but two overall observations may be suggested: (1) regional-scale data tend to be more uniform and outcrop data more clustered (e.g. compare Clair regional and core Cv values); and (2) the Hebrides data show a tendency for more clustered spacing distributions compared with the mainland (Assynt and Rhiconich) terranes. Franklin (2013) indicated that this effect is most pronounced at the outcrop scale (Fig. 6) and is likely due to the prominent influence of Mesozoic faulting in the Hebrides region.
Plot of Coefficient of variation (Cv) v. Fracture Intensity for outcrop, mesoscale (virtual model) and regional (lineament maps) datasets from mainland (Assynt and Rhiconich), Hebrides and Clair.
Topology results
The topology analyses were carried out on a range of onshore and offshore samples including drill core, outcrop images, seismic attribute and regional datasets. Figure 7 shows a summary of the topology values that have been obtained from Clair and other Rona Ridge (Lancaster) core and the Assynt Terrane (see Table S3 for full results). All basement samples show connected fracture networks with CB values >>1 which is the threshold CB (connections per branch) for a connected network. Most outcrop and core samples show a predominance of Y node-dominated fracture networks (Fig. 7).
Fracture topology results from Clair drill core samples, the greater Rona Ridge, and the Assynt Terrane (outcrops and regional lineament samples). Examples of the three scales sampled are shown: regional scale; outcrop scale; and core scale.
Discussion
This study, which reports the largest attribute dataset ever assembled for basement-hosted fractures, shows that the Scottish mainland exposures broadly show similar scaling and connectivity properties to the Clair basement and the greater Rona Ridge. Aperture scaling from all three areas (Hebrides, mainland and Clair) can be described by a power-law distribution when appropriate censoring and truncation of individual datasets are taken into account (Fig. 3 and Table S1). A number of individual datasets, which tend to be those with lower sample numbers, may be equally or slightly better described by log-normal distributions. Fracture length datasets from both onshore and offshore may be described by either power-law or log-normal distributions. Length distributions are known to be particularly prone to censoring and truncation (Odling et al. 1999). However, Odling et al. (1999) and Dichiarante et al. (2020) have shown that a multi-scale analysis can help to confirm that power-law scaling is an appropriate choice to model the fracture length distributions. In the present study, the basement fracture lengths sampled in 1D show a scaling relationship across eight orders of magnitude and the 2D sample windows show consistent and comparable length distributions between onshore and offshore datasets. Fractures onshore and offshore show similar spatial characteristics as demonstrated by the Cv values. The fracture topology analyses show similar levels of connectivity between onshore and offshore basement terranes. We note that the fracture networks at three scales (regional, outcrop and core) from kilometre to centimetre scale appear to be strongly Y-node dominated, which supports the conclusions that the networks are all well connected (Sanderson and Nixon 2015). Y-node dominated connectivity might be expected in relatively massive basement rocks which have multiple fracturing events in which large apertures form. Later formed fractures will tend to abut against the earlier fractures rather than cross-cut, hence Y-node development is favoured over X-node. The power-law fracture distributions are typical of massive crystalline rocks (e.g. Genter et al. 1997; Gillespie et al. 1999; Odling et al. 1999). The long history of brittle deformation and reactivation of structures within the Lewisian gneisses produced areas in which there are multiple fracture sets with power-law size distributions and good connectivity but, as has been noted previously, these attributes alone are not enough to make a viable fractured reservoir (Nelson 1985).
Our characterization demonstrates that certainly the onshore basement terranes provide a good first-order analogue for the offshore Clair basement and greater Rona Ridge. Importantly, however, our analysis has also shown that important differences do exist between the areas, e.g. the Hebrides has different aperture scaling to the mainland and Clair which we discuss in the following sections as it potentially provides further insight into what produces better reservoir potential in the basement gneisses. If it is accepted that our MLE analysis indicates a general power-law behaviour for the fracture aperture distributions, the large number of datasets collated in this study enables the overall scaling properties of the distributions to emerge. In most fracture studies there is generally high variability in scaling and fracture intensity between individual sample lines (e.g. see McCaffrey et al. 2003). In previous work, we have shown for basement lithologies at the outcrop scale that fracture distributions are affected by lithology and proximity to higher-order structures. Beacom et al. (2001) showed that fracture densities and clustering are higher in metasedimentary rocks compared with the more common intermediate to acidic gneisses. Pless et al. (2015) analysed a well-exposed basement outcrop in the Rhiconich Terrane and found that fracture density is higher within a 220 m envelope adjacent to the Kinlochbervie Fault (Fig. 2). The outcrop-scale datasets reported in this study are all deliberately taken from intermediate to felsic gneisses which minimizes significant variation caused by lithology. This lithology also dominates in the offshore basement (e.g. Holdsworth et al. 2018). The variation in fracture intensity of about an order of magnitude in the outcrop data for the mainland (Assynt and Rhiconich terranes) does include variation due to proximity to major structures (Figs 3 and 4).
An increase in fracture intensity with proximity to major structures explains the difference we see in the variability between the onshore and the offshore datasets. We find that the Clair core aperture dataset, of equivalent scale to the outcrop data, show similar power-law scaling to mainland Scotland with exponents in the range of 1–1.2. However, all of the Clair datasets plot in the higher fracture intensity range and do not show the lower intensity patterns displayed by the mainland. Specifically, the Clair data generally occupy the area defined by the grey box defined in Figure 4 whereas only the higher fracture intensity samples from the mainland do this – including those closer to major structures like the Kinlochbervie Fault (Figs 2 and 4). The Clair fracture intensity data have a much more limited spatial coverage compared with the mainland fracture sample lines in that they come from a single horizontal well that was drilled close to the top-basement interface; the Clair Ridge Fault. Holdsworth et al. (2019) also reported that the Clair core aperture distributions (the same datasets as plotted herein) show a systematic variation with highest fracture intensity in cores taken closest to the top basement interface. The above discussion and the findings of Holdsworth et al. (2019) show that variations in fracture intensity of about an order of magnitude in aperture distributions might be expected due to proximity to major structures. What this variation does not account for is the significant variation in scaling (slope of the lines) between the Hebrides aperture datasets and those of the mainland terranes and Clair. As we have shown in this study, the fracture apertures collected from the Hebrides, from both high- and low-intensity regions, show significantly lower scaling exponents (in range 0.5–0.8) compared with Clair or the mainland (1–1.2 (Fig. 4)). In simple terms, this means that in any sample we take from the Hebrides, we see more fractures with large aperture and relatively fewer with smaller apertures. Given that the fracture length distributions appear similar for all the datasets, we seek an explanation that can account for the presence of relatively more larger-aperture structures in the Hebrides. One explanation could be that the Hebrides has experienced more Mesozoic faulting, but there is no evidence from the data that the overall fracture intensities are higher here than on the mainland or at Clair. Using geological observations, we propose a simple conceptual model based on the development of fissures in basement blocks in the near-surface during the Mesozoic in order to account for the scaling differences we observe.
In recent related work, Holdsworth et al. (2019, 2020a) report structures and textures from offshore fracture fills that reveal the widespread development of steeply inclined to sub-vertical, rift-related tensile fissures in the basement lithologies of the Rona Ridge. They suggest that near-surface fissuring during rift-related faulting, as seen in modern rift systems – such as those exposed in Iceland (Kettermann et al. 2019) – allowed pervasive influx of clastic sediment fills from above and hydrothermal mineral fills from below. These partial sediment and vuggy mineral fills could act as natural props holding open fracture systems enabling long-term permeability pathways and facilitating hydrocarbon migration (Holdsworth et al. 2019, 2020a). In the following section, we explore whether this model might explain the different scaling properties that we see in onshore–offshore NW Scotland.
Our model is based on the following assumptions: (1) the fracture systems in the Lewisian basement and equivalents offshore are both cumulative products of multiple episodes of brittle deformation that produced shear, hybrid and tensile fractures – some of which display evidence for reactivation; (2) the basement was exposed at surface during its history for significant periods of time as indicated by the preservation of the basal Torridonian (c. 1.2 Ga Stoer, c. 1.04 Ga Torridan groups), Cambrian (c. 0.5 Ga) and Mesozoic (<0.3 Ga) unconformities; and (3) the basement experienced at least one (most likely several) phases of rifting while at surface that produced significant fissure-type fracturing with sediment and mineral infills (e.g. Beacom et al. 1999; Jonk et al. 2004; Holdsworth et al. 2019, 2020a). The model presented in Figure 8 shows a basement block with cover sediments (representing older sequences such as the Devonian–Carboniferous Clair Group, for example) that has been deformed by brittle deformation related to rifting. We know that many of the larger fractures onshore in the Hebrides (Franklin 2013) and offshore (Holdsworth et al. 2019, 2020a) are sediment filled, contain vuggy cavities in mineral fills, and show clear evidence for past fluid flow (mineralization) and even present-day fluid transport. These types of structures have been recorded in other settings where high-strength crystalline (e.g. Montenat et al. 1991) or carbonate rocks (e.g. Wright et al. 2009) are exposed at surface; sub-unconformity fissure fills and related structures are also widely recorded in active rift settings (e.g. Frenzel and Woodcock 2014; Kettermann et al. 2016, 2019; Koehn et al. 2019).
Conceptual model for fracture systems and their attributes developed in an uplifted basement block (see also Holdsworth et al. 2020a). Cartoon logs A’-A’ and B-B’ correspond to two hypothetical, horizontally deviated wells drilled through the block at different structural levels or through their onshore analogue equivalents exposed in outcrop.
Analogue modelling studies (e.g. van Gent et al. 2010; Holland et al. 2011) demonstrate that fissure structures which form open tensile fractures (with sediment infills) at surface, likely change character with depth transitioning through hybrid (shear tensile structures) to shear fractures at depth with a concomitant reduction in consistent fracture aperture. This variation in fissure/fault character with depth becomes important when considering the erosional level of the basement terranes of Scotland and the Rona Ridge at various times in their geological history (Fig. 8). We hypothesize that near-surface, large aperture tensile fractures with a more distributed deformation (a lower aperture exponent <1 and Cv <1) indicate a position near the top of a basement block. For example, a sample from Well A-A’ in Figure 8, or an onshore exposure located at an equivalent position. In contrast, where the faults and fractures intersected have more of a shear component with damage zones clustered around the larger fault structures (thus aperture exponents >1 and Cv > 1), it indicates that erosion levels are somewhat greater (Well B’-B’ in Fig. 8 or equivalent exposure). We suggest that less-eroded fault blocks represent the Hebridean basement terranes (and perhaps also the basement of Lancaster – see Holdsworth et al. 2020a) whereas the Clair basement and the mainland exposure represent more deeply eroded equivalents.
Further work is needed both on subsurface datasets and the onshore analogues to better constrain the speculative model proposed here. This study largely compiles datasets collected prior to our new understanding of the key role of fissuring in creating viable basement reservoirs. There is a need for new datasets that focus on the fissure structures to test this hypothesis, but at the moment it serves as a semi-quantitative predictor of the fracture attributes and hence also their fluid storage capacities and flow performance. Our model for the Rona Ridge, mainland and Clair basement fracture systems suggests a possible depth-dependent influence component on the basement fracture systems. While this is primarily due to a downward change in fissure and fault characteristics, it is the appreciation of the depth of erosion of the uplifted fault blocks in each of the rift episodes that is key to understanding the preserved fracture attributes and their influence on reservoir behaviour. Other factors that need to be explored include the effect on fracture attributes of the presence and thickness of cover sequence present during rifting, but our model provides a hypothesis that can be further tested. Fracture characterization of reservoir analogues can help to reduce uncertainties in the development of subsurface models that are created to determine drilling locations and quantifying the likely economic returns in terms of hydrocarbon production and resource in fractured basement fields such as Lancaster and Clair. However, we agree with Nelson (1985) when he said that ‘Finding fractures is not enough’. It is finding where the right type of fractures are preserved, in this case places where Mesozoic sub-unconformity fissures have formed and are preserved, that is key to a good reservoir in the offshore crystalline basement of NW Scotland.
Conclusions
One of the most extensive investigations of fault and fracture attributes collected from brittle structures in the onshore and offshore Lewisian Gneiss Complex rocks of NW Scotland shows that fracture sizes display power-law scaling of aperture and length attributes and are highly connected across a wide range of scales. The results show that the onshore fracture systems may be used as a good analogue for the basement reservoirs of the Rona Ridge and likely other fractured basement reservoirs worldwide. The high connectivity and size attribute scaling characteristics of the faults and fractures that may form in uplifted, crystalline basement rift blocks confirms that given the right geological history – notably the development and preservation of near-surface, rift-related fissure systems beneath unconformities – they may make good reservoir targets in their own right.
Acknowledgements
Casey Nixon is thanked for help with the topological analysis and Geospatial Research Ltd for an early version of Figure 1. We thank reviewer Dave Sanderson, Editor Simon Price and an anonymous reviewer for their insightful comments.
Author contributions
KM: conceptualization (lead), data curation (lead), formal analysis (lead), funding acquisition (supporting), investigation (lead), methodology (lead), supervision (supporting), writing – original draft (lead), writing – review & editing (lead); RH: conceptualization (supporting), funding acquisition (lead), investigation (supporting), project administration (lead), supervision (lead), writing – original draft (supporting), writing – review & editing (supporting); JP: data acquisition (supporting), data curation (supporting), formal analysis (supporting), investigation (supporting), writing – review & editing (supporting); BF: data acquisition (supporting), data curation (supporting), formal analysis (supporting), investigation (supporting), writing – review & editing (supporting); KH: visualization (supporting), writing – review & editing (supporting)
Funding
We thank the Clair Joint Venture (BP, Chevron Britain Limited, Chrysaor and Shell), Hurricane Exploration, Conoco Phillips, and NERC for supporting this work on basement fracture systems.
Data availability statement
The datasets generated during and/or analysed during the current study are not publicly available because the data are still to be analysed in further studies but are available from the corresponding author on reasonable request.
Scientific editing by Simon Price
- © 2020 The Author(s). Published by The Geological Society of London
This is an Open Access article distributed under the terms of the Creative Commons Attribution 4.0 License (http://creativecommons.org/licenses/by/4.0/)