Pre-drilling assessments of average porosity and permeability in the geothermal reservoirs of the Danish area

Abstract

Denmark constitutes a low-enthalpy geothermal area, and currently geothermal production takes places from two sandstone-rich formations: the Bunter Sandstone and the Gassum formations. These formations form major geothermal reservoirs in the Danish area, but exploration is associated with high geological uncertainty and information about reservoir permeability is difficult to obtain. Prediction of porosity and permeability prior to drilling is therefore essential in order to reduce risks. Geologically these two formations represent excellent examples of sandstone diversity, since they were deposited in a variety of environments during arid and humid climatic conditions. The study is based on geological and petrophysical data acquired in deep wells onshore Denmark, including conventional core analysis data and well-logs. A method for assessing and predicting the average porosity and permeability of geothermal prospects within the Danish area is presented. Firstly, a porosity-depth trend is established in order to predict porosity. Subsequently, in order to predict permeability, a porosity-permeability relation is established and then refined in steps. Both one basin-wide and one local permeability model are generated. Two porosity-depth models are established. It is shown that the average permeability of a geothermal prospect can be modelled (predicted) using a local permeability model, i.e. a model valid for a geological province including the prospect. The local permeability model is related to a general permeability model through a constant, and the general model thus acts as a template. The applied averaging technique reduces the scatter that is normally seen in a porosity-permeability plot including all raw core analysis measurements and thus narrows the uncertainty band attached to the average permeability estimate for a reservoir layer. A "best practice" technique for predicting average porosity and permeability of geothermal prospects on the basis of core analysis data and well-logs is suggested. The porosity is primarily related to depth, whereas the permeability also depends on porosity, mineralogy and grain size, which are controlled by the depositional environment. Our results indicate that porosity and permeability assessments should be based on averaged data and not raw conventional core analysis data. The uncertainty range of permeability values is significantly lower, when average values are used.

Kristensen et al. Geotherm Energy (2016) 4:6 DOI 10.1186/s4051701600486

Predrilling assessments ofaverage porosity andpermeability inthe geothermal reservoirs ofthe Danish area

Lars Kristensen1*, Morten Leth Hjuler1, Peter Frykman1, Mette Olivarius1,2, Rikke Weibel1, Lars Henrik Nielsen1 and Anders Mathiesen1

*Correspondence: lk@geus. dk

1 Geological Surveyof Denmark and Greenland (GEUS), ster Voldgade 10, 1350 Copenhagen, Denmark Full list of author information is available at the end of the article

2016 Kristensen et al. This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/

Web End =http://creativecommons.org/licenses/by/4.0/ ), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

Kristensen et al. Geotherm Energy (2016) 4:6

Background

Recent evaluations of seismic reection surveys and well data acquired during former geothermal and hydrocarbon exploration activities indicate that several sandstone-rich formations in the Norwegian-Danish Basin and the North German Basin contain substantial geothermal low-enthalpy resources (Mathiesen etal. 2009; Norden 2011; Kirsch etal. 2015). These basins are classic sedimentary basins characterized by long-term subsidence and inlling by sediments (Bertelsen 1980; Nielsen 2003; Bachmann etal. 2010; Lott et al. 2010). The geothermal potential of the Danish area has been addressed by GEUS in a regional assessment to the Danish Energy Agency and in a large number of customer reports prepared for district heating companies during the last 10years (e.g. Hjuler et al. 2014; Mathiesen et al. 2010a, b). Furthermore, GEUS has evaluated and assessed the geothermal potential of selected city areas in Denmark. These activities are reviewed in Mathiesen etal. (2010b) and Vosgerau etal. (in press).

In the Norwegian-Danish Basin and the North German Basin, the widely distributed Bunter Sandstone and Gassum formations constitute major geothermal reservoirs, but also formations with more local distribution, such as Skagerrak, Haldager Sand, Flyvbjerg and Frederikshavn formations, have geothermal potentials (Rgen et al. 2015; Mathiesen etal. 2010b, Nielsen etal. 2004). The geothermal potential is related directly to reservoir quality, which in the Danish onshore area traditionally is addressed by considering clay content, net sand thickness, porosity and permeability by means of wireline logs, core analysis data, well tests and seismic data (e.g. Mathiesen etal. 2013). The available seismic and well data point to thick reservoir formations and in addition, formation temperatures, porosities and permeabilities have proven to be sufficiently high for geo-thermal water production in large parts of Denmark.

Denmark is a low-enthalpy geothermal area with a temperature gradient of 2530C/ km, and only minor temperature anomalies are encountered in the subsurface (Balling etal. 1981; Mathiesen etal. 2013; Balling etal. 2014; Poulsen etal. submitted). The geo-thermal potential onshore Denmark has been mapped for the Bunter Sandstone (including Skagerrak), Gassum, Haldager Sand and Frederikshavn formations (Mathiesen etal. 2009, 2010b). The mapping considered distribution and estimated resources of these formations, and potential geothermal reservoirs were identied as sandstone layers having thicknesses greater than 25m in the depth interval 8003000m, corresponding to formation temperatures in the range 25100C. The permeability in deeper-seated reservoirs is considered too low for geothermal water production (Mathiesen etal. 2009, 2010b; Weibel et al. submitted). Presently, two plants produce from the Gassum formation; saline water (43C) is produced from a depth of 1250m at Thisted, whereas water of 48C is produced at Snderborg (depth 1200m). At a plant in Copenhagen (Margretheholm), saline water (74C) is produced from the Bunter Sandstone formation located at a depth of c. 2600m (Rgen etal. 2015). Each plant is congured with a production well and an injection well located about 1km from the producer, returns the cooled water.

A prognosis (pre-drilling assessment) of reservoir porosity and permeability for the three Danish geothermal plants was not carried out in a quantitative fashion, but the intention of this paper is to suggest a methodology that makes it possible to predict reservoir properties of new geothermal prospects, as the lack of accurate predictions is

Page 2 of 27

Kristensen et al. Geotherm Energy (2016) 4:6

regarded as the bottle-neck for new projects. In contrast, the depth, thickness and temperature can in most areas be estimated with an uncertainty of less than 10%. Our work thus aims at presenting a method for assessing the average porosity and permeability of geothermal prospects within the Danish area. The method is considered a best practice approach.

Especially in poorly explored regions, prediction of reservoir parameters along with an assessment of the geological development is needed prior to drilling. Pre-drilling assessments of uncertainties and general reservoir parameters (thickness, porosity, permeability and temperature) are crucial for estimating the geothermal potential in such areas. Seismic interpretation is essential in order to determine depth and thicknesses. Setting up a work programme for a geothermal project thus requires a careful integration of existing geological and geophysical information.

The present study is based on geological and petrophysical data acquired in deep wells onshore Denmark. The database is comprehensive and it contains widely distributed data, both vertically and geographically. However, the data density varies considerably; both closely spaced and sporadically distributed data form part of the database. A statistical approach to data analysis is not feasible when dealing with such a database with large variations in data density and distribution. As an alternative the data analysis is herein based on an empirical approach. Two porosity-depth trends along with a number of empirical porositypermeability relationships that are established using existing welllogs and core analysis data, and then used for reservoir characterization. These trends and relations are considered to have a predictive potential and may be used for characterizing a geothermal prospect, provided that the approximate depth to the anticipated geothermal reservoir is known. Our empirical porositypermeability trends are established on the basis of core analysis data available from wells drilled and cored during the period c. 19502010. The large time span implies that quality of core analysis data, i.e. porosity and permeability measurements on plug samples, varies considerably but generally the quality is fair to good. The data material indicates that the correspondence between core porosity and core permeability data is not perfect, but reasonable porositypermeability correlations may be obtained, despite the rather scattered data. As permeability can be predicted with least condence, the main emphasis is laid on improving permeability prediction in areas with sparse data coverage.

With respect to geothermal exploration in the Danish onshore area, two issues are particularly important for the local district heating companies holding the geothermal exploration license areas: (1) the geological uncertainty prior to drilling the rst well needs to be thoroughly assessed and (2) the prognosis for permeability and transmissivity in a potential reservoir needs to be as well constrained as possible. In some cases it has been concluded that the geological uncertainty and the associated exploration risks are too high for justifying the drilling of the rst and costly geothermal exploration well.

Previous studies ofrelevance

Many workers have discussed the correlation between porosity and permeability measurements, and despite somewhat imperfect correlations, they succeeded in establishing either exponential relationships (e.g. Tiab and Donaldson 2004) or trends dened by power functions (e.g. Doyen 1988; Mavko and Nur 1997).

Page 3 of 27

Kristensen et al. Geotherm Energy (2016) 4:6

A number of methods for determining permeability from porosity are described in literature, including empirical approaches and various modelling techniques. Already in 1927, Kozeny proposed to predict permeability from porosity, geometry of the pores and the specic surface of the solids in contact with the uid. The Kozeny equation was later modied by Carman (1937, 1956) to become the KozenyCarman equation, which estimates the permeability of a porous media based on a grain size distribution. Block (1991) outlined a predictive approach based on the correlation between porosity and permeability in which the predictive applicability is constrained by the limits dened by a calibration dataset that includes various petrographic variables or parameters. Block (1991) even provided an equation in which the permeability is expressed as a function of grain size, sorting and rigid grain content. Evans etal. (1997) discussed a permeability prediction methodology that is based on a combination of empirical and geological modelling approaches. Originally the Kozeny equation was derived for clay-free sand with high porosity, and for that reason Walderhaug etal. (2012) suggested a modied Kozeny equation for predicting permeability in quartz-rich sandstones, even with high clay content. This modied Kozeny equation includes a parameter reecting the type of pore system, and usually this parameter is to be considered a constant for a specic sandstone unit. With respect to modelling, Vadapalli et al. (2014) present fractal and Monte Carlo simulation approaches for permeability prediction.

Porosity andpermeability assessments

Our porosity modelling is based on data from the Bunter Sandstone and Gassum formations. The permeability modelling is demonstrated by data from the Gassum formation.

Quantication of reservoir permeability, which is the single most critical factor for geothermal uid extraction, is complicated as very few insitu measurements are available and moreover, permeability logs are not available from wells in the Danish onshore area.

Determination of reservoir permeability requires good-quality well test data, i.e. tests based on sufficiently long test intervals and ow/build-up periods of long duration. However, such high-quality test data only exist for a few Danish onshore wells (e.g. at Stenlille and Tnder; Fig. 1), and well test data are generally not available. So as an alternative, it is suggested to assess the permeability from an analysis of porosity log data supplemented by information from porositypermeability relationships based on core analysis data. Core analysis data commonly include gas permeabilities measured at laboratory conditions, but conversion into corresponding reservoir uid permeabilities (at eld scale) is not a straightforward task. This conversion normally involves transformation of the core scale gas permeability into liquid permeability followed by upscaling from laboratory scale to eld scale. The correction from gas to liquid permeability is normally carried out as a Klinkenberg correction of laboratory measurements (Klinkenberg 1941). The Klinkenberg correction is quite important for low-permeability sandstones, where a factor of 0.5 may apply, but for high permeabilities (>100mD) the Klinkenberg correction is less pronounced and a factor in the order of 0.80.9 may apply (Tanikawa and Shimamoto 2009; Duan and Yang 2014). The use of core permeability values for upscaling to a full reservoir volume assumes that the sampling of the core samples is representative, which is not always the case. The arithmetic averaging used

Page 4 of 27

Kristensen et al. Geotherm Energy (2016) 4:6

Page 5 of 27

herein is based on the concept of horizontal, along layer ow, and occurring in a layered formation, which has been shown to be closely comparable to that from more advanced upscaling schemes (Kazemi etal. 2012).

The Danish onshore core analysis database is considerably larger than the number of well tests, and the purpose of analysing core permeability data is to provide an estimate of the average permeability that characterizes a particular layer, but also to address the uncertainty of a predicted permeability estimate. A single permeability value interpreted from well test data is, however, not directly comparable to an average permeability based on core data, but the two permeability assessments are comparable at a relative scale, on the assumption that the thickness of test interval corresponds to the length of the cored interval.

Geological setting

The Norwegian-Danish Basin contains a thick Upper PermianMesozoic succession of sedimentary rocks (Fig.1). The basin was formed in Late CarboniferousEarly Permian time by crustal stretching, succeeded by deposition of Rotliegendes coarse-grained clastic sediments and later by Zechstein salts (Stemmerik etal. 1987; Nielsen 2003; Vejbk 1989, 1997; Michelsen and Nielsen 1993; Pharaoh etal. 2010). The Permian period followed basinal subsidence governed primarily by thermal cooling and local faulting, and the Triassic sediments include a succession of sandstones, mudstones, carbonates and

Kristensen et al. Geotherm Energy (2016) 4:6

evaporites, whereas the Jurassic strata are dominated by mudstones and sandstones (Berthelsen 1980; Michelsen and Clausen 2002; Nielsen 2003). The Mesozoic deposition is inuenced by dierential uplift across the basin and occasionally, halokinesis has resulted in local deviations from the regional uplift trend (Japsen etal. 2007). The regional trend points to gradually increasing uplift toward the northeast. The Cretaceous sediments encompass limestones, mudstones, sandstones, marls and chalks. The Cenozoic succession comprises a series of primarily mudstones, sandstones and limestones.

Bunter sandstone formation

The Lower Triassic Bunter sandstone formation is widely distributed in the North German Basin and in parts of the Norwegian-Danish Basin (Fig.2) (Berthelsen 1980; Nielsen and Japsen 1991; Michelsen and Clausen 2002; Bachmann etal. 2010). In Denmark, this formation is especially relevant for geothermal exploration in the southern and eastern parts of the country due to sufficient thickness and generally good reservoir properties (Mathiesen etal. 2010a, b). Formation thickness varies from about 300m in the Margretheholm wells in the eastern part of the Norwegian-Danish Basin to 600700m in the central parts of the basin and to c. 200 m in the Tnder wells in the North German Basin (Fig.1). In general the formation thins towards the RingkbingFyn High, and it is locally absent on the High (Nielsen and Japsen 1991; Olivarius etal. 2015b). The reservoir quality of the formation varies considerably, as the lithology varies from sandstones with a subordinate amount of mudstones, to siltstones and mudstones with few sandstone layers (Berthelsen 1980; Clemmensen 1985). However, the reservoir quality of the thickest sandstone layers is generally good to excellent (Olivarius etal. 2015a). The Bunter Sandstone formation was deposited in an arid to semi-arid climate, and the sandstones represent deposition by uvial channel systems and eolian dunes, whereas the mudstones mainly were deposited in lakes and on ood plains (Clemmensen 1985). Eroded material was supplied from the RingkbingFyn High by alluvialuvial systems and by wind transport across the North German Basin (Olivarius etal. 2015b; Clausen and Pedersen 1999).The sandstone layers are generally continuous, especially the eolian sandstones, and consist mostly of very ne to medium-grained sand (Olivarius et al. 2015b). Similarly, eroded material from the Fennoscandian Shield was transported into the Danish Basin (Bertelsen 1980; Olivarius and Nielsen 2016). The sandstone deposits are ne to coarse grained and generally, grain size increases in a northeasterly direction, i.e. towards the anks of the basin. In the southern and central parts of the Danish Basin, the deposits belong to the Bunter Sandstone formation and to the Skagerrak formation in the northeastern part.

Gassum formation

The Upper TriassicLower Jurassic Gassum formation is present in most of the Norwegian-Danish Basin (Fig.2), except above large salt structures and on the RingkbingFyn High (Nielsen 2003). It is also present in the northern North German Basin, but here its patchy and shallow occurrence makes it less suitable for geothermal exploration. In general, the formation thickness varies from c. 30m to more than 300m (Michelsen etal. 2003). Nielsen (2003) described the geological development of the Gassum formation in detail: the formation consists of shallow marine, uvial and estuarine sandstones

Page 6 of 27

Kristensen et al. Geotherm Energy (2016) 4:6

Page 7 of 27

interbedded with marine and lagoonal mudstones, and also siltstones and minor coal beds are present. During the Late Triassic to Early Jurassic, the sedimentation was inuenced by repeated sea-level uctuations and rivers transported eroded material from the Fennoscandian Shield into the Norwegian-Danish Basin. Several sandstone layers were deposited as regressive shoreface sands, but also as uvial sands deposited in river channels or estuaries. The Gassum formation was deposited in a humid climate, and the sandstone layers consist mainly of ne to medium-grained sand. The mudstones were

Kristensen et al. Geotherm Energy (2016) 4:6

primarily deposited in lagoons, lakes and oshore shelf areas. The reservoir quality of the sandstones is generally good to excellent due to the presence of rather continuous sandstone beds with low clay content (Weibel etal. submitted).

The inuence ofdiagenesis onreservoir quality

The porosity and permeability of the Bunter Sandstone and Gassum formations generally decrease with increasing burial depth, due to mechanical compaction and diagenetic alterations (Olivarius etal. 2015a, b; Weibel etal. submitted). Diagenesis (e.g. cementation) may even aect the sandstones at shallow burial, but otherwise mechanical compaction is dominant at shallow burial (Table1). Chemical compaction becomes more pronounced at depths greater than 23km owing to higher temperature and pressure. Quartz cement is considered the major porosity reducing element in quartz-dominated sandstones at depth greater than 2000m (Ehrenberg 1990). The sandstones of the Gassum formation are aected by pronounced quartz cementation at greater depths, but also the presence of carbonate cement (calcite, siderite and ankerite) and clays aects reservoir quality as well (Weibel et al. submitted). In the Bunter Sandstone formation quartz cementation is very limited due to shallow burial, so the reservoir quality is mainly inuenced by calcite, anhydrite and halite cements along with clay minerals (Olivarius et al. 2015a). This dierence in the diagenetic development of the Bunter Sandstone and Gassum formations may be attributed to the dierences in the depositional environments caused by the dierent climatic conditions during deposition, i.e. arid conditions versus humid conditions (see Olivarius et al. 2015a and Weibel et al. submitted).

In general, the permeability reduction with burial depth is more pronounced for ne-grained sandstones than for coarse grained, but also the content of detrital clay, sorting and other elements related to variations in the depositional environments and source area proximity aect permeability (Weibel et al. submitted; Olivarius et al. 2015a, b). The presence of diagenetic cements may result in substantial permeability reduction, as the cement reduces the size of the pore throats. A synopsis of the key factors aecting porosity and permeability compiled from Olivarius etal. (2015a) and Weibel etal. (submitted) is presented in Table1, covering both the Bunter Sandstone and Gassum formations. The tabulation summarizes the results of two diagenesis studies based on analyses of thin sections, scanning electron microscope images, cores and cuttings samples.

Methods

A 5-step procedure has been developed in order to improve predictions of reservoir parameters valid for geothermal prospects in areas with poor data coverage, i.e. areas with no wells, but with sufficient seismic data to determine approximate depths and thicknesses of potential reservoirs.

First step is to establish a regional model for porosity prediction.

In step 2, a permeability model is established.

In step 3, the permeability model is rened.

Step 4 includes the establishment of local permeability models.

Step 5 is implemented to reduce the uncertainty range of the permeability.

Page 8 of 27

Kristensen et al. Geotherm Energy (2016) 4:6

Page 9 of 27

Table 1 Factors inuencing porosity andpermeability inthe Bunter Sandstone andGassum formations

Factors Depth<2500m Depth>2500m

Depth (max. burial depth)

Bunter Sandstone Fm Mechanical compaction. Core data represent a narrow depth interval

Gassum Fm Mechanical compaction. The porosity reduction with depth is comparable to the porosity estimated from a mechanical compaction curve. Depth is an important factor, when dealing with permeabilities of shallowly buried sandstones

Bunter Sandstone Fm No core data Gassum Fm Chemical compaction. The porosity reduction with depth may be higher than indicated by a mechanical compaction curve. Permeability is not that depth dependent

Grain size Bunter Sandstone Fm Porosity reduction is highest for very ne-grained sandstones and less for coarser-grained sandstones. Increasing grain size leads to higher permeabilities

Gassum Fm Limited inuence on porosity. Increasing grain size leads to higher permeabilities

Gassum Fm The porosity reduction is highest for very ne-grained sandstones and less for coarser-grained sandstones

With respect to permeability, grain size has less inuence than at shallow depths. Increasing grain size still leads to higher permeabilities, however

Detrital clay Bunter Sandstone Fm Presence of inter-granular clay and clay clasts reduce porosity, but not substantially since much microporosity is present within the clays. Even small amounts of inter-granular clay reduce permeability considerably, whereas larger amounts of clay clasts are needed to produce a similar reduction in permeability

Gassum Fm Presence of inter-granular clay and clay clasts reduce porosity. Detrital clays and/or clay clasts are often present. Clays and clay laminae lower the permeability, since some of the pore throats are very narrow

Gassum Fm The amount of clay increases with depth, but it has only minor eect on porosity, since the clay grows on the expense of other minerals. In addition, clay clasts increase the eect of compaction

High amounts of detrital clays result in reduced permeability

Cement Bunter Sandstone Fm Pervasive carbonate, anhydrite or halite cement may reduce porosity signicantly, whereas patchy carbonate cement does not have a notable eect on porosity. Pervasive carbonate, anhydrite or halite cement occludes pores and thus prevents uid ow. However, most commonly the cement is patchy and has only limited eect on uid ow

Gassum Fm Siderite and calcite cement occasionally result in larger porosity reduction than mechanical compaction. Presence of siderite cement leads to a marked reduction in permeability

Gassum Fm Pronounced porosity reduction with depth due to the presence of quartz and/or ankerite cement

The permeability is markedly reduced where authigenic illite is present. Kaolinite has limited reducing eect on permeability. The permeability is primarily reduced by quartz and ankerite cement

Coatings Bunter Sandstone Fm Sandstones with thick iron-oxide/ hydroxide coatings are characterized by high permeability. The coatings may also preserve porosity

Gassum Fm Sandstones with chlorite coatings have high permeability, unless other cement types are present

Gassum Fm Chlorite coatings may preserve porosity and permeability

Overall, the porosity distribution depends on the depositional environment and the maximum burial depth. Similarly, the permeability is strongly related to the depositional environment, which controls the distribution of grain sizes, the abundance of detrital clays, and the amount/type of cementing minerals etc. Based on information from Olivarius etal. (2015a) and Weibel etal. (submitted)

Our method for permeability prediction is based on averaging core and log data in order to derive average values for potential reservoir layers. The basic data for developing the 5-step procedure are Danish data as described below.

Welllog data, core analyses andwell test data

Well-log data of variable quality have been acquired in hydrocarbon and geothermal exploration wells drilled in the Danish onshore area during the last c. 60 years. The porosity variation is determined from the interpretation of well-log data, and log interpretation results form the basis for calculating net sand properties. Both the total and

Kristensen et al. Geotherm Energy (2016) 4:6

the eective porosity of the reservoir units are considered. The total porosity (PHIT) is the total volume of inter-granular porosity plus clay bound water, whereas the eective porosity (PHIE) does not include the water bound by the clay particles. The core porosity data are directly comparable to the log-derived total porosity data. The term net sand or potential reservoir layer is herein dened as sandstone having minimum 15% porosity and a clay or shale content of less than 30%. These cut-os are adopted from hydrocarbon exploration practice and applied herein to ensure that only potential reservoir sands with sufficient storage capacity and permeability are considered. Data representing non-reservoir sandstone and mudstone intervals are thus removed prior to further data analysis.

Only a limited number of well test data are available for permeability assessments in the Danish area, since most of the data for this study originates from dry hydrocarbon exploration wells. These wells were not tested, but usually logged and sometimes also cored. Accordingly, the present study focuses on alternative ways of predicting permeability, namely a combined analysis of petrophysical well-log data and routine permeability measurements on core plugs, with the objective to provide an average gas permeability of a particular reservoir unit. For that reason a direct link between the gas permeability measured in the laboratory and the actual liquid permeability of the reservoir in the subsurface is needed for precise pre-drilling estimates of the expected performance of a given geothermal reservoir.

All available core analysis data from the Danish onshore area have been considered in order to construct a robust database of porosity and permeability data. In general, measurements were performed according to the API RP-40 standard (American Petroleum Institute), i.e. He-porosity was measured at unconned conditions and gas permeability was measured at a conning pressure of c. 2.8MPa (400psi), and at a mean nitrogen gas pressure of c. 0.15MPa (c. 1.5bar absolute). It is not always possible to obtain information about the test conditions, as the large number of analyses were performed over several decades and by various companies. These companies were, however, expected to follow the existing standards and it is thus assumed that measurements were conducted at conditions corresponding or comparable to the API RP-40 standard. Both Klinkenberg corrected and uncorrected gas permeabilities are included in the database, but primarily uncorrected permeabilities form part of the porosity permeability plots presented herein.

Gas permeabilities measured in the laboratory do not equal reservoir permeabilities, whereas well test data are considered to provide reservoir permeability estimates on the condition that the height (thickness) of the test interval is well-known. Strictly speaking well test data only supply a transmissivity measure (i.e. reservoir height multiplied by reservoir permeability). Permeability interpreted from well test data is commonly up to 2 times higher than the corresponding core permeability (e.g. Wolfgramm etal. 2008); for example presence of fractures in the reservoir rock usually results in permeability enhancement. The presence of fractures or overall inhomogeneities cannot be validated for geothermal reservoir rocks onshore Denmark, but our interpretations of well test data from 5 wells indicate a liquid permeability that is about 1.5 times higher than the core permeability (Fig.3). A prerequisite is that both core and well test permeabilities exist for the same interval.

Page 10 of 27

Kristensen et al. Geotherm Energy (2016) 4:6

Page 11 of 27

Prediction ofporosity

First step in the 5-step procedure is to establish a porosity-depth model. The basic principle behind the porosity model is an assumed porosity-depth relationship, which previously has been documented for other basins (e.g. Gluyas and Cade 1997).

The porositydepth model (step 1)

Seismic data provide information about depth and thickness of a geothermal prospect, but usually information about the average porosity must be modelled using porosity data from nearby wells. Porosities have primarily been interpreted from well-log data, and these data form the basis of establishing a correspondence between porosity and depth. The interpreted porosity log for each well is therefore averaged in order to determine the average porosity for specic reservoir intervals. A reservoir interval presumes a minimum porosity (herein>15%) and a maximum shale content (herein<30%) as described above. Therefore cut-os were applied prior to calculating average porosities, i.e. the calculated porosities are average net porosities. The eect of applying cut-os is that reservoir layers within each formation are identied and then assigned a porosity value on a well-to-well basis. One porosity-depth trend is observed for the Bunter Sandstone formation and another for the Gassum formation (Figs.4, 5). The depth scale has been modied due to dierential uplift across the Norwegian-Danish Basin (Japsen

Kristensen et al. Geotherm Energy (2016) 4:6

Page 12 of 27

and Bidstrup 1999; Japsen etal. 2007), and hence the standard depth scale is replaced by estimated maximum burial depth.

Prediction ofpermeability

In step 24 of the 5-step procedure a regional porositypermeability model is established (step 2), rened (step 3) and adapted to be applicable to local conditions (step 4). Step 5 suggests a method to narrow the permeability uncertainty range.

Kristensen et al. Geotherm Energy (2016) 4:6

Page 13 of 27

The initial permeability model (step 2)

Initially an empirical porositypermeability relationship is established for each geological formation on the basis of available conventional core analysis data from all relevant

Kristensen et al. Geotherm Energy (2016) 4:6

wells using power-law-based curve tting in line with the Kozeny law (Kozeny 1927). The resulting function used for calculating permeabilities is termed the initial permeability model and is illustrated in Fig.6 and expressed in Eq.1;

where kini is the initial permeability (in mD) and ~ is the porosity (in %). The log-derived porosity was used as input data, and in this way a synthetic permeability log is assigned to all study wells, including the uncored wells. The permeability was calculated both from the eective and the total porosity, but the permeability curve calculated from the eective porosity is considered the most reliable permeability estimate, since it takes into account that the actual reservoir permeability is reduced in shaly and clayey intervals. However, the total porosity approximates the eective porosity in clean or almost shale-free reservoir intervals, and the synthetic permeability curves calculated from the total and eective porosity, respectively, are not signicantly dierent in reservoirs with (very) low clay content.

Page 14 of 27

kini = 3.47 104 ~4.38,

(1)

Kristensen et al. Geotherm Energy (2016) 4:6

The general permeability model (step 3)

The next step is to average both the log-derived porosity curve and the log-derived permeability curve for the reservoir interval within each formation. Prior to this, shale and porosity cut-os are applied in order to exclude non-reservoir data originating from clayey and cemented sections. These data do not contribute to the reservoir performance, but would otherwise inuence the calculation of the average permeability if included. We assume a basinal setting with horizontal layering where the uid ow is parallel to layer boundaries. In our opinion such a ow regime together with the cut-os justify the use of arithmetic average for the operation.

A cross-plot between averaged porosity and permeability data points forms the basis of dening a new porositypermeability relationship, the general permeability model, (Fig.6) represented in Eq.2:

where kG is the general permeability (mD), and ~ is the porosity (%). The porosity determined from the porosity-depth model (step 1) is used as input data for permeability modelling. Equation2 is considered to be more representative of the full reservoir sections (Fig.6; bold line). The calculated values are based on well-log data covering the reservoir section, and not only the cored parts of the formation. The introduction of the general permeability model also means that a specic permeability estimate becomes slightly higher than that calculated from the initial permeability model. The few available eld test measurements conrm this observation.

The local permeability model (step 4)

The general permeability model represents the entire Danish onshore area and constitutes a template for constructing more rened local permeability models. Reservoir intervals of a specic formation in a local region may be characterized by deviating permeability distributions compared to the general trend (step 3). If local core permeability data are available, such data should therefore be used to calibrate the general permeability model to the local area. The use of the general permeability model in constructing local models is implemented with Eq.3:

where kL is the local permeability (in mD) and C is a constant controlled by the local permeability variations, and it may be determined from permeability measurements on core plugs from one or two local wells.

Core porosity and permeability data from selected wells drilled throughout the Danish onshore area have been averaged for each well and then compared to the general trend line (Fig. 7). The gure focuses on deviations from the general model, and it appears from the gure that permeability averages deviate by a factor of up to four as indicated by the upper and lower bounds, even though some data points t the trend line. However, this spread is caused by the use of regional data representing the entire Danish onshore area, which masks local and often more limited permeability variations. The uncertainty connected to a local dataset is less than indicated by the regional data

Page 15 of 27

kG = 4.43 104 ~4.36,

(2)

kL = C kG,

(3)

Kristensen et al. Geotherm Energy (2016) 4:6

Page 16 of 27

set due to e.g. facies uniformity, and local permeability models are therefore needed for delivering locally adjusted permeability predictions.

Permeability uncertainty range (step5)

Even a local, calibrated porositypermeability model is associated with an uncertainty range owing to the variations in lithology and diagenetic alterations within a limited local region (geological province). Such a restricted area may only be represented by one or two wells with wireline logs and no cores, and therefore it is difficult to determine a suitable uncertainty range related to the permeability estimate representative of any local area on the basis of such a limited database. In order to get an impression of the local variability, a selected area with good well coverage has been studied in detail. 20 wells have been drilled within a restricted area at the Stenlille gas storage facility located in the eastern part of the Norwegian-Danish Basin (Fig.1), and a large database comprising both core analysis data and well-logs exists (these data are all available from the GEUS archives). Ten of these wells are cored in the upper part of the Gassum formation, and evaluation of 3D seismic data, correlation of well-logs and interpretation of cores indicate that relatively uniform geological conditions prevailed in the area during deposition. When all unprocessed core porosity and core permeability data representing the upper part of the Gassum formation are plotted in one cross-plot, the data points show

Kristensen et al. Geotherm Energy (2016) 4:6

an expected substantial scattering (Fig.8). Hence the averaged core porosity and core permeability data are used to establish a local Stenlille model for this part of the Gassum formation and most importantly, also to assess the uncertainty range of the local porositypermeability relationship (Fig.9). The gure shows that the scatter is notably reduced when applying the averaging technique, but also that the distribution is displaced towards lower average permeabilities compared to the general model. This displacement could be due to dierences in grain size as discussed later. The key message from the gure is, however, that the high and low bounds of the uncertainty band are limited by factors of 2 and 0.5 of the model mid-line.

Results

Result ofporosity modelling

The result of the porosity modelling is the establishment of two regional porosity-depth trends as presented in Figs.4 and 5.

Result ofpermeability modelling

Applying the averaging technique for predicting permeability on the basis of core permeability data and well-log interpretations adds to the applicability of porositypermeability relations. A best practice method is suggested for predicting the average permeability of a potential geothermal reservoir. The average permeability of a geo-thermal prospect is modelled (predicted) using the closest local permeability model, i.e. a permeability model valid for a nearby geological province. This permeability value along with estimated net sand thickness are considered the key factors, when assessing the geothermal potential of a particular prospect. From our experience, the geo-thermal resource is of economic value if the reservoir transmissivity is greater than 10

Page 17 of 27

Kristensen et al. Geotherm Energy (2016) 4:6

Page 18 of 27

Darcy-metre. This assessment is in line with indicative values published in Seibt and Kellner (2003). The net sand thickness may be found from seismic isochores combined with extrapolation of net-to-gross ratios. Acquisition, interpretation and depth conversion of seismic data are thus pre-requisites for addressing the geothermal potential of a prospect. The determination of net sand thickness in a particular well is exemplied in Fig.10.

Porositypermeability plots based on conventional core analysis data are commonly scattered as shown and discussed above, but the use of the averaging technique presented herein reduces the scatter and narrows the width of the uncertainty band associated with a predicted average permeability value. It is thus recognized that the uncertainty of the average value is signicantly smaller than the uncertainty of a single measurement.

Application ofthe methodology

The use of the methodology is exemplied by the derivation of a local model as presented in Fig.11. The model is illustrated by mid-line, high and low bounds, dening an uncertainty envelope, and it is assumed that the model is a representative of the Gassum formation at the northern basin margin. Averaged core analysis data from wells located

Kristensen et al. Geotherm Energy (2016) 4:6

Page 19 of 27

in this area (Brglum-1 and Flyvbjerg-1; Fig.1) point to permeabilities higher than indicated by the general model. A local model based on Eq.3 was therefore established with the scope of predicting the average permeability of the Gassum formation in the BrglumFlyvbjerg areas, expressed by Eq.4:

where kL is the local permeability (in mD) and the constant (C=1.7) is determined from analysis of local core analysis data.

kL = 1.7 kG,

(4)

Kristensen et al. Geotherm Energy (2016) 4:6

Page 20 of 27

The conventional core analysis data from the Brglum-1 and Flyvbjerg-1 wells are also plotted after applying a porosity cut-o (Fig.11). The mid-line honours calculated averages of the core analysis data from the two wells. Data points plotting below the uncertainty envelope were, however, excluded prior to calculating averages, because these data points represent sandstones with large amounts of ne-grained material (shale, mud-stone and siltstone) not contributing to reservoir performance. The width of the uncertainty band is transferred from our Stenlille study, i.e. factors of 2 and 0.5 dene the high and low bounds. The uncertainty band does not address the spread in the actual core analysis data, but it points out an uncertainty range that is related to the average permeability for a typical reservoir, presuming that the reservoir sandstones are almost shale-free. Furthermore, the example demonstrates that it is difficult to narrow the uncertainty range on the basis of a limited set of conventional core analysis data.

Discussion

A number of porositypermeability trends have been analysed to obtain a better understanding of the relationship between porosity and permeability. It is recognized that the porosity and permeability distribution is aected by burial depth and grain size (Table1),

Kristensen et al. Geotherm Energy (2016) 4:6

but the inuence of depositional environment, sediment source area, transport distance and diagenesis far from well control is not yet fully understood and is therefore difficult to quantify. Models for estimating porosities and permeabilities have, nevertheless, been established, but are associated with uncertainty.

Prediction ofporosity

The formation brines of the potential reservoir sandstones in the Danish onshore area are characterized by hydrostatic pressure, so sandstone porosity is not aected by over-pressure. Knowledge of the pressure distribution across the Danish onshore area comes from analysis of well test data. Our study indicates, however, that the sandstone porosity more likely is related to geological factors such as burial depth, lithofacies, clay content and diagenesis. In addition, the porosity-depth trend is aected by dierential uplift and erosion; thus present-day depths were corrected for uplift using exhumation data presented in Japsen and Bidstrup (1999) and Japsen etal. (2007). For simple modelling purposes in undrilled areas, the porosity may be considered to be controlled by depth alone, provided that present-day depths are replaced by estimated maximum burial depths.

One porosity-depth trend is observed for the Bunter Sandstone formation and another for the Gassum formation (Figs.4, 5), despite the fact that local deviations from the general trends are observed. These deviations are presumably related to geological factors as outlined above. An empirical porosity-depth relationship by Gluyas and Cade (1997) is also plotted to compare the current porosity-depth trends with a suggested normal compaction curve. The Gluyas and Cade curve presumes no uplift and reects mechanical compaction in uncemented sandstones from hydrocarbon elds with high porosity and no overpressure. These general porosity-depth trends (cf. Figs.4, 5) are suggested for predicting the average porosity of potential reservoir sandstones in undrilled areas. A complete match between depth and porosity has not been fully achieved although maximum burial depths are applied, since the porosity also depends on factors other than depth, e.g. gain size (Table1). However, for modelling purposes the use of a porosity-depth relation is considered satisfactory for porosity prediction and assessment. Consequently, modelled porosities are associated with uncertainty, even if the reservoir depth is well-known prior to drilling. In specic areas where the porosity distribution is well-known from local well and core data (e.g. at Tnder; Fig.1), a local porosity-depth trend should be used.

Porosity ofthe Bunter sandstone formation

Porosities of the Bunter sandstone formation are generally lower than derived from the mechanical compaction curve published by Gluyas and Cade (1997), indicating that most sandstones of the Bunter Sandstone formation contain clays and diagenetic cements (Fig.4). The Bunter sandstone formation retains relatively high porosity at greater depths, as detrital clays contain microporosity and diagenetic iron-oxide/ hydroxide coatings seem to retard quartz cementation (Table1; Olivarius etal. 2015a). The Tnder area in the southernmost part of Denmark (Fig. 1) is a potential site for geothermal exploitation and here the Bunter Sandstone reservoir (situated at a depth of 15002000m) is characterized by average porosities (c. 22%) that are higher than indicated by the general trend line. The porosity development at Tnder is mainly due to

Page 21 of 27

Kristensen et al. Geotherm Energy (2016) 4:6

a dominance of aeolian deposition that favoured generation of very well-sorted and clay-free sandstones, but also the presence of nitrogen gas contributes to porosity preservation during burial cf. the Tnder-3 and -4 well completion reports (Mrsk 1981; Dong 1983). The Kegns-1 and Kvrs-1 wells, also located in the southern part of Denmark, have encountered the Bunter Sandstone formation at a similar depth (15002000m), but here the porosities (c. 17%) are somewhat lower than indicated by the general trend for the Bunter Sandstone formation. Cutting descriptions reported in the well completion reports (Mrsk 1985; Texaco 1985) indicate that this porosity dierence most likely is attributed to the presence of more ne-grained sandstones with higher clay content at Kegns and Kvrs.

Porosity ofthe Gassum formation

The porosities are generally higher for the Gassum formation than for the Bunter Sandstone formation when considering shallow depths (<c. 3000m), but the reservoir intervals of the Gassum formation show a steeper porositydepth gradient. The overall porosity deterioration with depth is about ve porosity units (%) per 1000m (Fig.5), and this trend is presumably a result of mechanical compaction combined with the eect of diagenetic alterations (Table1). As an example, the reservoir sandstones at the Stenlille gas storage facility are characterized by porosities (c. 30%) that are higher than indicated by the general trend line (see Fig. 1 for location). This porosity development is probably due to dominance of well-sorted shoreface sandstones. A maximum burial depth of c. 2100m is estimated for the Gassum formation at Stenlille, and the present-day depth of c. 1500m is due to uplift. In most wells, the porosity deterioration with depth corresponds to the general trend irrespective of correction for uplift (Figs.1, 4). In the Brglum-1 well, in contrast, the porosity of the reservoir interval is extraordinarily high (33 %) compared to burial depth, primarily due to the occurrence of unconsoli-dated sandstone layers in this interval. The reservoir intervals of the Gassum formation likewise exhibit porosities that locally are lower than modelled from the porosity-depth trend, presumably because of variations in diagenetic development (e.g. the Voldum-1 well in which the porosity is c. 18%).

Prediction ofpermeability

Reservoir performance, including ow rates, is closely related to porosity and permeability, and information about permeability is essential for determining reservoir transmissivity unless well test data are available. Correct prediction of sandstone permeability is a challenge for potential geothermal reservoirs, since uid ow depends on a number of factors. Compaction processes, precipitation of cement and presence of detrital clay lead to smaller pore throat sizes, thus lowering the resulting reservoir permeability.

Usually there is a relatively clear correlation between core porosity and core permeability data, but this correlation is somewhat ambiguous due to scattered data resulting from dierences in the depositional environment, clay content and diagenetic development. The variability in the core analysis data means that a perfect correlation between core porosity and core permeability data cannot be obtained. Despite these uncertainties, this computed porositypermeability relation forms the basis of calculating a log-based permeability curve for each well using the log-derived porosity as input data.

Page 22 of 27

Kristensen et al. Geotherm Energy (2016) 4:6

The observed variation in permeability for a xed porosity value (Fig.7) is most likely related to variations in the original depositional environment when the sandstone layers were formed, expressing itself in dierences in diagenetic development, lithology, grain size distribution and clay content (Table1). These variations emphasize the need of incorporating geological data, preferably both sedimentological and palaeogeographical data, and also the need for establishing local porositypermeability models representative of local geological provinces characterized by relatively uniform geological development. Such local permeability models are herein suggested for predicting the average permeability of geothermal prospects. In fact, access to local permeability data from wells located close to the prospect in question is essential for using our method for predicting permeability.

In addition to the porosity-depth model, it could be relevant to set up a general, basin-wide permeability-depth model. However, such a model will be rather uncertain, because depth is not the single most important factor for assessing the permeability. Sandstone porosity and grain size are also responsible for permeability development (cf. the Kozeny law). It is, nevertheless, suggested to link permeability to depth in 2 steps:

1. Depth is converted into porosity using a porosity-depth model.2. Porosity is converted into permeability using a porositypermeability relationship dened by a local permeability model. A local permeability model accounts for differences in grain size and diagenetic development (Table1). When combining these two items, the permeability of the Gassum formation may be expressed as function of depth (Eq.5):

where kL is the local permeability (mD), kG is the general permeability (mD), C is a constant, ~ is the porosity (%) and Z is the maximum burial depth (m). The porosity-depth equation for the Gassum formation (Fig.5) is also applied. Inserting C=1 in Eq.5 will give an initial assessment of the permeable conditions at depth.

The methods presented herein are developed on the basis of data from Danish geo-thermal reservoirs. The methodology for predicting average porosity and permeability is, however, considered applicable to similar sandstone reservoirs in other settings, but it has to be proved with additional data.

Permeability uncertainty range based onlocal eld data

The width of the uncertainty band related to average permeability is addressed by analysing core analysis data from the local Stenlille eld with several cored and logged wells located closely together. A petrophysical evaluation of the Gassum formation in the Stenlille-1 well is illustrated in Fig.10. The character of the Gassum formation at Stenlille is considered representative of geothermal reservoirs found in the formation onshore Denmark. The Gassum formation consists of a number of sandstone layers interbedded with clay-rich intervals. Hence the Gassum formation is assigned a lithology column along with porosity and permeability curves as shown in Fig.10. The porosity is

Page 23 of 27

kL = C kG = C 4.43 104 ~4.36 = C 4.43 104 (37 0.0054 Z)4.36,

(5)

Kristensen et al. Geotherm Energy (2016) 4:6

interpreted from well-log data, whereas the permeability is calculated from the porositypermeability relationship shown in Fig.8.

A local Stenlille permeability model was initially established on the basis of averaged core porosity and permeability data from 10 Stenlille wells cored in the upper part of the Gassum formation, corresponding to the main gas storage reservoir (Fig.9). In addition, the core analysis data are relevant for calibrating the log-based porosity and permeability interpretations (Fig.10). Less core material exists from the lower part of the Gassum formation and the geological conditions that prevailed during deposition of the lower part dier from those of the upper part (Nielsen etal. 1989; Nielsen 2003). With respect to permeability, the core analyses show that the overall Stenlille permeability is approximately a factor 2 less than the permeability estimated from the general, basin-scale model (Fig.9). The lower average permeability values observed at Stenlille are presumably caused by dominance of ner-grained sandstone compared to the sandstones dening the general permeability model. The presence of ne-grained reservoir sandstones at Stenlille is evidenced by sedimentological core log descriptions of the Stenlille-1 core (Nielsen etal. 1989). The permeabilities of these sandstones (Fig.8) are generally lower than commonly seen in Gassum formation sandstones deposited elsewhere, and the lower permeability level at Stenlille may therefore be linked to grain size.

The averaged core analysis data from the Stenlille wells suggest that the uncertainty range can be expressed by multipliers of 2 and 0.5 (Fig.9). This procedure involves multiplying and dividing the local Stenlille permeability trend by a factor of 2 in order to determine an upper and lower bound that delineate the data points, i.e. an envelope ranging from the lowest to the highest average permeability values. As the core permeability data do not follow a normal distribution, the uncertainty cannot be described by conventional standard deviation and instead, we suggest applying this uncertainty envelope. We assume the Stenlille dataset to be relevant for assessing the uncertainty range for a typical Gassum geothermal reservoir at any location in Denmark.

Conclusions

The major outcome of this study is the development of a Best practice technique for predicting porosity and permeability in sandstone reservoirs located in areas with poor data coverage. A 5-step procedure has been developed for assessing porosity and permeability in the Danish area:

1. Establishing a regional porosity-depth model.2. Establishing an initial permeability model based on a conventional porositypermeability plot using core analysis data from all cored study wells.

3. Establishing a general permeability model where both porosity and permeability are replaced by averaged log-derived data.

4. Using the general permeability model as a template for establishing local permeability models.

5. Using the Stenlille model to describe the uncertainty range of a permeability estimate obtained from a local permeability model.

Page 24 of 27

Kristensen et al. Geotherm Energy (2016) 4:6

Furthermore, the following conclusions can be provided based on our present knowledge:

Application of the 5-step technique using averaged porosity and permeability data signicantly reduces the scattering of data points that is normally seen in conventional porositypermeability plots. This observation is related to the fact that the uncertainty of an average value is signicantly lower than the uncertainty of a single measurement.

The available database is comprehensive and it comprises both widely and sporadically distributed data, meaning that a statistical approach to data analysis is not feasible. Instead an empirical approach has been used.

The average porosity of a particular sandstone reservoir in a geothermal prospect is herein considered to be related primarily to depth. Models for porosity prediction are established for the Bunter sandstone and the Gassum formations.

A permeability model that is based on averaging both log-derived porosities and permeabilities is introduced (general permeability model). The porosity determined from the porosity-depth model is used as input data for the permeability modelling.

We consider local permeability models suitable for permeability prediction. Local permeability models account for variations in burial depth, grain size and diagenetic alterations.

A permeability-depth modelling method is suggested. The permeability is, however, related to factors other than depth such as grain size, detrital clay content and diagenetic development, including presence of cement and authigenic clays (Table1).

We have compared permeabilities interpreted from well test data with core permeability measurements whenever possible (5 wells). Seemingly, core permeability data measured in the laboratory resemble the reservoir permeability, provided that an appropriate correction or upscaling factor is applied.

The geothermal potential of a particular prospect should be assessed using modelled porosity and permeability values combined with thicknesses derived from seismic interpretations. Thus an assessment of the transmissivity is essential, and in this context modelled permeabilities and estimated net sand thicknesses are key input parameters.

The uncertainty range related to the average permeability of a geothermal prospect is addressed using a comprehensive local eld dataset from the Stenlille gas storage facility. Local eld data are geographically and geologically constrained and therefore suitable for analysing a local uncertainty range connected to a specic site. Our data analysis suggests that the uncertainty range related to the average permeability can be expressed by multipliers of 2 and 0.5.

Authors contributions

LK carried out the petrophysical evaluation, performed the data analysis, participated in the construction of the porositypermeability models and drafted the manuscript. MLH participated in the log analyses and contributed considerably to the drafting of the manuscript. PF contributed to the development and design of the permeability prediction model. MO and RW performed the diagenetic investigations, dealing with the Bunter Sandstone formation (MO) and the Gassum formation (RW). LHN and AM participated in the design of the study and contributed to the regional interpretations of the results and to the drafting of the manuscript. All authors have read and commented early drafts and have approved the nal manuscript. All authors read and approved the nal manuscript.

Author details

1 Geological Survey of Denmark and Greenland (GEUS), ster Voldgade 10, 1350 Copenhagen, Denmark. 2 Department of Geoscience, Aarhus University, Hegh-Guldbergsgade 2, 8000 Aarhus, Denmark.

Page 25 of 27

Kristensen et al. Geotherm Energy (2016) 4:6

Acknowledgements

This paper is published with permission of the Geological Survey of Denmark and Greenland and is an outcome of the project The geothermal energy potential in Denmarkreservoir properties, temperature distribution and models for utilization under the Sustainable Energy and Environment programme founded by the Danish Council for Strategic Research. We would like to thank professor Kim Esbensen (GEUS) for valuable discussions and for constructive comments on the manuscript text. Stefan Slberg (GEUS) improved the artwork considerably.

Competing interests

The authors declare that they have no competing interests.

Received: 1 November 2015 Accepted: 21 March 2016

References

Bachmann GH, Geluk MC, Warrington G, Becker-Roman A, Beutler G, Hagdorn H, Hounslow MW, Nitsch E, Rhling H-G,

Simon T, Szulc A. Triassic. In: Doornenbal H, Stevenson A, editors. Petroleum geological atlas of the Southern Permian Basin area. Rotterdam: European Association of Geoscientists and Engineers Publications; 2010. p. 14973. Balling N, Kristensen JI, Breiner N, Poulsen KD, Rasmussen R, Saxov S. Geothermal measurements and subsurface temperature modelling in Denmark. GeoSkrifter. 1981;16:172.

Balling N, Poulsen SE, Bording TS, Fuchs S, Mathiesen A. Surface temperatures in Denmarkmeasurements and modelling. Abstract and poster presentation. Copenhagen: Energi-og miljkonference; 2014.

Berthelsen F. Lithostratigraphy and depositional history of the Danish Triassic. Geological Survey of Denmark; 1980. Series

B, No 4.

Block S. Empirical prediction of porosity and permeability in sandstones. Am Assoc Pet Geol Bull. 1991;75:114560. Carman PC. Fluid ow through granular beds. Trans Inst Chem Eng. 1937;15:15066.

Carman PC. Flow of gasses through porous media. London: Butterworth Scientic Publications; 1956.

Clausen OR, Pedersen PK. Late Triassic evolution of the southern margin of the Ringkbing-Fyn High. Mar Pet Geol.

1999;16:65365.

Clemmensen LB. Desert sand plain and sabkha deposits from the Bunter Sandstone formation (L. Triassic) at the northern margin of the German Basin. Geol Rundsch. 1985;74:51936.

Dong. Tnder-4 gas exploration well. Well summary. Dansk Olie og Naturgas A/S; 1983.

Doyen PM. Permeability, conductivity and pore geometry of sandstone. J Geophys Res. 1988;93:772940.

Duan Q, Yang X. Experimental studies on gas and water permeability of fault rocks from the rupture of the 2008 Wenchuan earthquake, China. Sci China Earth Sci. 2014;57(11):282534.

Ehrenberg SN. Relationship between diagenesis and reservoir quality in sandstones of the Garn formation, Haltenbanken, Mid-Norway Continental Shelf. Am Assoc Pet Geol Bull. 1990;74:153858.

Evans J, Cade C, Bryant S. A geological approach to permeability prediction in clastic reservoirs. In: Kupecz JA, Gluyas J,

Block S, editors. Reservoir quality prediction in sandstones and carbonates AAPG Memoir 69; 1997. p. 91101. Gluyas J, Cade CA. Prediction of porosity in compacted sands. In: Kupecz JA, Gluyas J, Block S, editors. Reservoir quality prediction in sandstones and carbonates. AAPG Memoir 69; 1997. p. 1928.

Hjuler ML, Erlstrm M, Lindstrm S, Nielsen LH, Kristensen L, Mathiesen A, Bidstrup T. Extended evaluation of possible geothermal reservoirs in the Helsingr area including geological data from Helsingborg and resund. Contribution to an evaluation of the geothermal potential. Geological Survey of Denmark and Greenland Report 2014/29; 2014. Japsen P, Bidstrup T. Quantication of late Cenozoic erosion in Denmark based on sonic data and basin modelling. Bull

Geol Soc Den. 1999;46:7999.

Japsen P, Green PF, Nielsen LH, Rasmussen ES, Bidstrup T. Mesozoic-Cenozoic exhumation events in the eastern North

Sea Basin: a multi-disciplinary study based on palaeothermal, palaeoburial, stratigraphic and seismic data. Basin Res. 2007;19:45190.

Kazemi A, Shaikhina D, Pickup G, Corbett P. Comparison of upscaling methods in a heterogeneous carbonate model. SPE

154499. EAGE Annual Conference and Exhibition, Copenhagen, June 47; 2012.

Kirsch R, Balling N, Fuchs S, Hese F, Hjuler ML, Kristensen L, Mathiesen A, Nielsen CM, Nielsen LH, Oermann P, Poulsen

NE, Rabbel W. GeoPower: Varme fra undergrunden til SnderjyllandSchleswigErdwrme fr die Region Sdjt-landSchleswig. Technical Report. Schriftenreihe: LLUR SHGeologie und Boden; 21. Landesamt fr Landwirtschaft, Umwelt und lndliche Rume des Schleswig-Holstein; 2015.

Klinkenberg LJ. The permeability of porous media to liquids and gases, drilling and production practice, American Petroleum Inst.; 1941. p. 20013.

Kozeny J. ber kapillare Leitung des Wassers im Boden. Akademi der Wissenschaften in Wein, Sitzungsberichte, Abt. 2a,

136; 1927. p. 271306.

Lott GK, Wong TE, Dusar M, Andsbjerg J, Mnning E, Feldman-Olszewska A, Verreussel RMCH. Jurassic. In: Doornenbal H,

Stevenson A, editors. Petroleum geological atlas of the Southern Permian Basin area. Rotterdam: European Association of Geoscientists and Engineers Publications; 2010. p. 17593.

Mathiesen A, Kristensen, L, Bidstrup T, Nielsen LH. Vurdering af det geotermiske potentiale i Danmark. Geological Survey of Denmark and Greenland Report 2009/59; 2009.

Mathiesen A, Kristensen, L, Bidstrup T, Nielsen LH. Evaluation of possible geothermal reservoirs in the Tnder area. Contribution to an evaluation of the geothermal potential. Geological survey of Denmark and Greenland Report 2010/90 (condential); 2010a.

Mathiesen A, Nielsen LH, Bidstrup T. Identifying potential geothermal reservoirs in Denmark. Geol Surv Denmark Greenland Bull. 2010;20:1922.

Page 26 of 27

Kristensen et al. Geotherm Energy (2016) 4:6

Mathiesen A, Kristensen, L, Nielsen CM, Weibel R, Hjuler ML, Rgen B, Mahler A, Nielsen LH. Assessment of sedimentary geothermal aquifer parameters in Denmark with focus on transmissivity. In: Proceedings European geothermal congress 2013. Pisa; 2013.

Mavko G, Nur A. The eect of a percolation threshold in the Kozeny-Carman relation. Society of exploration geophysicists. Geophysics. 1997;62:14802.

Michelsen O, Nielsen LH. Structural development of the Fennoscandian Border Zone, oshore Denmark. Mar Pet Geol.

1993;10:12434.

Michelsen O, Clausen OR. Detailed stratigraphic subdivision and regional correlation of the southern Danish Triassic succession. Mar Pet Geol. 2002;19:56387.

Michelsen O, Nielsen LH, Johannessen P, Andsbjerg J, Surlyk F. Jurassic lithostratigraphy and stratigraphic development onshore and oshore Denmark. In: Ineson JR, Surlyk F, editors. The Jurassic of Denmark and Greenland. Geological Survey of Denmark and Greenland Bulletin 1; 2003. p. 147216.

Mrsk. Tnder-3 Completion Report. Dansk Boreselskab A/S; 1981.

Mrsk. Kvrs-1 Completion Report. Mrks Olie og Gas A/S; 1985.

Nielsen LH, Larsen F, Frandsen N. Upper TriassicLower Jurassic tidal deposits of the Gassum formation on Sjlland,

Denmark. Geological Survey of Denmark; 1989. Series A, No 23.

Nielsen LH, Japsen P. Deep wells in Denmark 19311990. Lithostratigraphic subdivision. Geological Survey of Denmark.

DGU; 1991. Series A, No. 31.

Nielsen LH. Late TriassicJurassic development of the Danish Basin and Fennoscandian Border Zone, Southern Scandinavia. In: Ineson JR, Surlyk F, editors. The Jurassic of Denmark and Greenland. Geological Survey of Denmark and Greenland Bulletin 1; 2003. p. 459526.

Nielsen LH, Mathiesen A, Bidstrup T. Geothermal energy in Denmark. Review of survey activities 2003. Geol Surv Denmark Greenl Bull. 2004;4:1720.

Norden B, editor. Geothermal energy utilization in low-enthalpy sedimentary environments. Scientic Technical Report

11/06. German Research Centre for Geosciences (GFZ), Potsdam; 2011.

Olivarius M, Weibel R, Hjuler ML, Kristensen L, Mathiesen A, Nielsen LH, Kjller C. Diagenetic eects on porositypermeability relationships in red beds of the Lower Triassic Bunter Sandstone formation in the North German Basin. Sed Geol. 2015a;321:13953.

Olivarius M, Weibel R, Friis H, Boldreel LO, Keulen N, Thomsen TB. Provenance of the Lower Triassic Bunter Sandstone formation: implications for distribution and architecture of aeolian vs. uvial reservoirs in the North German Basin. Basin Res. 2015b. doi:http://dx.doi.org/10.1111/bre.12140

Web End =10.1111/bre.12140 .

Olivarius M, Nielsen LH. Triassic paleogeography of the greater eastern Norwegian-Danish Basin: constrains from provenance analysis of the Skagerrak formation. Mar Pet Geol. 2016;69:16882.

Pharaoh TC, Dusar M, Geluk MC, Kochel F, Krawezyk C, Krzywiec P, Schech-Wenderoth M, Thybo H, Vejbk O, Diederik van Wees J. Tectonic evolution. In: Doornenbal H, Stevenson A, editors. Petroleum geological atlas of the Southern Permian Basin area. Rotterdam: European Association of Geoscientists and Engineers Publications; 2010. p. 2557. Poulsen SE, Balling N, Bording TS, Mathiesen, A, Nielsen SB (submitted) Regional-scale subsurface temperature modelling with inverse methodology: application to Danish sedimentary basins.

Rgen B, Ditlefsen C, Vangkilde-Pedersen T, Nielsen LH, Mahler A. Geothermal energy use, 2015 Country Update for

Denmark. Proceedings World Geothermal Congress 2015, Melbourne; 2015.

Seibt P, Kellner T. Practical experience in the reinjection of cooled thermal waters back into sandstone reservoirs. Geothermics. 2003;32:73341.

Stemmerik L, Frykman P, Christensen OW, Stentoft N. The Zechstein carbonates of southern Jylland. In: Brooks J, Glennie,

KW, editors. Petroleum Geology of North West Europe; 1987. p. 36574.

Texaco. Kegns-1 Completion Report Well 5410/5-1. Texaco; 1985.

Tiab D, Donaldson EC. Petrophysics. Theory and practice of measuring reservoir rock and uid transport properties.

Amsterdam: Elsevier; 2004.

Tanikawa W, Shimamoto T. Comparison of Klinkenberg-corrected gas permeability and water permeability in sedimentary rocks. Int J Rock Mech Min Sci. 2009;46(2):22938.

Vadapalli U, Srivastava N, Vedanti N, Dimri VP. Estimation of permeability of a sandstone reservoir by a fractal and Monte

Carlo simulation approach: a case study. Nonlin Process Geophys. 2014;21:918.

Vejbk OV. Eects of asthenospheric heat ow in basin modellering exemplied with the Danish Basin. Earth Planet Sci

Lett. 1989;95:97114.

Vejbk OV. Dybe strukturer i danske sedimentre bassiner. Geologisk Tidsskrift. 1997;4:131.

Vosgerau, H, Mathiesen, A, Andersen, MS, Boldreel, LO, Hjuler, ML, Kamla, E, Kristensen L, Laier, T, Pedersen, CB, Pjetursson, B, Nielsen, LH. A new WebGIS portal with geological and geophysical data for exploration of deep geothermal energy. Review of survey activities 2015. Geol Surv Denmark Greenl Bull; (in press).

Walderhaug O, Eliassen A, Aase NE. Prediction of permeability in quartz-rich sandstones Examples from the Norwegian shelf and the Fontainebleau sandstone. J Sediment Res. 2012;82:899912.

Weibel R, Olivarius M, Kristensen L, Friis H, Hjuler ML, Kjller C, Mathiesen A, Nielsen LH. Predicting permeability of low enthalpy geothermal reservoirs: a case study from the Upper TriassicLower Jurassic Gassum formation, Norwegian-Danish Basin; (submitted).

Wolfgramm M, Rauppach K, Seibt P. Reservoir-geological characterization of Mesozoic sandstones in the North German

Basin by petrophysical and petrograpical data. Geol Wiss. 2008;36(45):24965.

Page 27 of 27