A Methodology for Assessing the Favourability of Geopressured-Geothermal Systems in Sedimentary Basin Plays: A Case Study in Abruzzo (Italy)

We exploit the concept of the geothermal favourability, widely used for hydrothermal and EGS systems, to present an innovative methodology for assessing geopressured-geothermal resources occurring in terrigenous units in sedimentary basin plays. Geopressured-geothermal systems are an unconventional resource for power trigeneration exploiting three forms of energy from hydrocarbons, hydrothermal fluids, and well-head overpressure. This paper is intended to be a practical analytical framework for the systematic integration of the relevant data required to assess geopressured-geothermal resources. For this purpose, innovative parameters were also implemented in the methodology. The final result is the favourability map for identifying prospective areas to be further investigated for the appraisal of the geopressured-geothermal potential. We applied our methodology to the foredeep-foreland domains of the Apennines thrust belt in the Abruzzo region (central Italy). We analysed hundreds of deep hydrocarbon wells in order to create 3D geological and thermo-fluid dynamic models at a regional scale as well as to obtain information on the pressure regimes and on the chemistry of the system. The final favourability map for the Abruzzo case study is a first attempt at ranking these kinds of unconventional geothermal resources in a region that has been historically explored and exploited mostly for hydrocarbons.


Introduction
Geothermal exploration is a complex, time-consuming, and expensive activity. Integrating geological, geochemical, and geophysical data can speed up exploration stages. Very few reference studies are available regarding the assessment at regional scale of geopressured-geothermal resources, which industrial interest is recently increasing worldwide.
We present a new methodology to assess the favourability of geopressured-geothermal systems occurring in terrigenous units in sedimentary basin plays. The aim is to drive effective exploration to the most promising areas throughout a systematic data integration. The research was developed within the framework of the Geothermal Atlas of Southern Italy Project [1].
Geopressured-geothermal systems (hereafter also referred to as "geopressured") are an unconventional resource for power trigeneration. They exploit three forms of energy [2]: (i) chemical energy from the combustion of hydrocarbons, (ii) thermal energy from hydrothermal fluids, and (iii) kinetic energy from well-head overpressure due to abnormal geopressured regimes. This resource is of particular interest due to the possibility of improving the economic feasibility of an industrial geothermal project or of uneconomic/depleted hydrocarbon wells.
The USA focused on geopressured systems from the 1970s to the 1990s above all for industry and produced extensive knowledge of the geopressured system of the Gulf of Mexico (e.g., [3,4]). Such industrial interest has recently been renewed with projects aimed at demonstrating their commercial feasibility, as in Louisiana [5]. With regard to Europe, Hungary has been developing geopressured projects [6] while cogeneration plants for the direct use of geothermal heat and dissolved methane are in operation. Previous works have described the prospective factors for the assessment of geopressured systems (e.g., [7]), mainly aimed at the Gulf of Mexico. Italy has also been involved in some works, for example Alimonti and Gnoni [8] presented a study on the heat recovery from depleted wells in a well-known geopressured field in the Po Plain (Northern Italy).
In the recent scientific literature, the geothermal favourability concept has been widely proposed in order to study geothermal systems, focusing on hydrothermal and enhanced geothermal systems (EGS). Various approaches have been exploited in order to quantitatively integrate different sources of data, usually organized in information layers in a GIS environment [9][10][11][12][13][14][15][16][17][18]. The favourability consists in a semiquantitative data integration for identifying prospective areas to be further investigated for the appraisal of the geothermal potential.
Our methodology is based on the integration of layers of evidence by Index Overlay. Data analysis by the Index Overlay method has been widely proposed in literature (e.g., [9,16]) for the assessment of other types of geothermal systems. We provide a novel tool for assessing geopressured resources that considers specific prospective factors. We applied our methodology to the foredeep-foreland domains of the Apennines thrust belt in the Abruzzo region (central Italy). This is one of the first attempts to assess geopressured systems at the regional scale in Italy. Favourability maps were computed in order to assess the geopressured resources that would be suitable for power production.
The Abruzzo case study is also important due to the possibility of developing geothermal projects in a region belonging to the Adriatic petroleum province [19][20][21] and in a gas (methane)-prone area, characterized by low geothermal gradients.
We analysed hundreds of deep hydrocarbon wells, and we carried out geological modelling and coupled thermo-fluid dynamic numerical simulations. Information on the pressure regimes and on the chemistry of the system was also obtained. The final favourability map for the Abruzzo region was computed which provides a ranking of the most prospective areas.
Once having described the general approach, the article focuses on the application to the test site in Abruzzo in order to better explicate the various steps of the methodology. Finally, the obtained maps are analysed and discussed.

Geopressured-Geothermal Systems: State of the Art and Prospective Factors
A geopressured-geothermal system is constituted by a hydrothermal reservoir with a higher pore pressure than the hydrostatic reference and contains dissolved gaseous hydrocarbons. The first patent for exploiting this resource was filed in 1966 [22]. The trigeneration of energy derives from the abnormal formation pressure, the occurrence of dissolved methane, and the heat transferred by the fluid. Due to various technical problems (e.g., the sustainability of the exploitation), geopressured systems are still considered unconventional, i.e., requiring technological development. Figure 1 shows the classic flow diagram for trigeneration as proposed by Hughes [23]. The literature contains a few studies regarding the exploration of this resource worldwide (e.g., [24]) and several works from related congresses, held in the USA, and technical reports (e.g., [7,[25][26][27][28]). The most studied geopressured systems are located onshore of the Gulf of Mexico, particularly in Texas and Louisiana. In the USA, following the oil crisis in the 1970s, the DOE (Department of Energy) funded two important research programs [27]. Numerous abandoned hydrocarbon wells were tested, and new wells were drilled in prospective areas. The most important productive well was the "Pleasant Bayou 2" located in Texas, a test site in which the technical feasibility of power production was demonstrated [29].
The concomitant occurrence of hydrothermal resources and hydrocarbons, in abnormal pressure regimes, makes the genesis of geopressured systems particular to specific geological environments. The play concept [30][31][32] is of primary importance in the favourability assessment.  Figure 1: Flow diagram of the hybrid power system for the trigeneration from geopressured-geothermal resources (modified from [23,91]).

Geofluids
The most suitable plays are mainly related, but not limited, to clastic depositional environments. Although sedimentary basin plays are not usually suitable for high-temperature hydrothermal resources, their geological features favour the occurrence of strata overpressure and hydrocarbon accumulations. Table 1 summarizes the main prospective factors for the characterization of geopressured systems in such plays.
In a sedimentary basin environment, compaction disequilibrium is the predominant process that leads to pore (or formation) overpressure, and this is related to subsidence, burial, and compaction of sediments. The overpressure gradient is directly proportional to the deposition rate and inversely proportional to the permeability [33]. These basin sectors that experience the fast deposition of interbedded clays and sands are favoured. Other processes that influence the formation of pressure are the thermal expansion of brine, the occurrence of hydrocarbons, and tectonics [33,34].
Loucks et al. [3] classified the Gulf of Mexico geopressured system into three main pressure regimes: (i) hydrostatic, (ii) "soft" geopressured, and (iii) "hard" geopressured. The latter, with gradients higher than 0.7 psi/ft (15.83 bar/100 m) and marked by a change in electrical resistivity, is also characterized by an increase in the thermal gradient at its top. The authors related this observation to the decrease in thermal conductivity as a consequence of higher porosity, following the theory proposed by Lewis and Rose [35], given a constant heat flow.
The sedimentary processes underlying such geopressured sedimentary basins also promote gaseous hydrocarbons, mainly methane, which occur in the formation as dissolved gas in brines or also as "free" exsolved gas in saturated conditions. The continuous burial of an enormous amount of organic matter is key to the formation of hydrocarbons. In geopressured systems, the methane exploited is mainly in solution in the hydrothermal fluid, and the gas/water ratio and the methane solubility are fundamental. The P/T and salinity conditions of the fluid strongly influence the solubility of methane; i.e., the amount of exploitable gas for a certain flow rate of the fluid. The methane solubility is directly proportional to the pressure and temperature and inversely proportional to the salinity of the fluid [36].
Geothermal exploration, in particular for power production, searches for the hottest hydrothermal resource and highest geothermal gradients. In the case of geopressured systems, the geological conditions that usually facilitate their development are unfavourable to high-geothermal gradients, such as those in volcanic or intrusive geothermal plays. For example, fast sedimentation environments are related to a decrease in the thermal gradient. The thermal regime is, however, a fundamental factor in deriving the economic value of the geopressured resource, since thermal energy plays a major role in trigeneration.
Permeability is another key factor. In geopressured systems, the lithology and the primary permeability play important roles, characterizing (i) the clayey seals, with a low permeability that favours the development of overpressure regimes, and (ii) the productive intervals, usually in sandy formations. The structural setting is a key factor in the evolution of such systems, for example, the possibility of characterizing permeable and nonpermeable faults as well as growth faults.

Methodology
GIS spatial analysis for mapping the geothermal potential or favourable areas has been extensively used following different approaches in order to speed up the exploration stages [9-18, 37, 38].
In cases of abundant data availability, the spatial layers in a GIS environment can be treated on the basis of statistical estimations. Otherwise, the processing of layers can rely on the knowledge of experts. The first case is known as the data-driven method, while the second is referred to as the knowledge-driven method [39].
In many cases, the information available makes it impossible to carry out a probabilistic analysis. We therefore applied a knowledge-driven method using the Index Overlay (IO) technique to combine geological, geophysical, and geochemical information. The IO is a simple and flexible way to linearly integrate different spatial information ensuring a common scale of value. The resulting map is obtained from (1), where F is the favourability for each pixel, W i is the weight for the ith map, and S ij is the score for the jth class of the ith map [39]: The workflow setup for the computation of the favourability map is organized into three stages, as shown in Figure 2. It starts with a preliminary geological analysis of the play. Indeed, the methodology focuses on the sedimentary terrigenous basin plays and the study area should be selected accordingly.
The second stage of computation requires the collection and processing of geological, well logs, geochemical, and geophysical datasets. The following thematic inputs are properly set up: (i-ii) depth of the top and base of the geopressured-geothermal reservoir, (iii) depth of the top and base of basin deposits, (iv) isobaths of the target temperature, (v-vi) temperature at the Earth's surface and at the top of the basement underlying the reservoir, (vii) digital elevation model. (viii) formation pressure, and (ix) fluid and gas geochemistry.
The thematic inputs are combined by means of GIS spatial analysis tools to obtain the layers of evidence.
The layers of evidence are the spatial representation of the main prospective parameters described in Section 2. The methodology includes the following layers: (i) the effective geopressured reservoir, (ii) the thermal regime, (iii) the pressure regime, (iv) the deposit thickness, and (v) the geochemistry.
The five layers of evidence are combined to produce the final favourability map for geopressured systems as a result of the last stage of the workflow. The layers are in turn scored and weighted following the Index Overlay method ((1)). The classification for each layer of evidence consists of identifying five ranges of score values (classes). The classes are scored from 1 to 5, "very low" (less favourable area) to "very high" (most promising area), respectively. In order to combine the layers of evidence, each was weighted with values whose sum is equal to 1 ( Table 2). The weights, classes, and scores were set based on generic features of terrigenous sedimentary basins.
In this work, all the computations were performed in the Open Source Quantum GIS environment exploiting SAGA and GRASS GIS tools [40]. The resolution of the combined maps is the same, with the grid nodes of each layer overlapping.

Effective Geopressured-Geothermal
Reservoir. The effective reservoir concept was initially proposed in Trumpy et al. [9] for hydrothermal conventional systems in carbonates. Here, we adapted the idea of "effective" reservoir in order to develop a new concept for geopressured systems in sedimentary basin plays: the geopressured-geothermal effective reservoir. This layer of evidence is intended to assess only that part of a geopressured reservoir with a temperature suitable for geothermal exploitation. In this paper, a threshold value of 90°C was set as target, considering the technology of a conventional power production plant. The threshold temperature is chosen accordingly to the aim of exploration (e.g., power production, and direct use of heat). Basically, the idea is to disregard those sectors of reservoir that are colder than the threshold temperature. The layer is computed by means of a layer intersection between the depth of the   4 Geofluids  5 Geofluids 90°C isotherm and the base of the reservoir. Where the isotherm is deeper than the base of the reservoir, i.e., no effective reservoir occurs, the corresponding areas are neglected from the computation and considered as not favourable. Conversely, if the 90°C isotherm rests above the base, an effective reservoir is identified and the depth of the top is recorded in this layer of evidence. The result is a grid layer of a ranked depth of the top of the effective geopressured-geothermal reservoir.
The ranking classes (score, S) are related to the depth to be drilled in order to reach the top of the effective reservoir: the shallower the top, the higher the favourability. The choice of limit values ( Table 2) was driven by economic considerations regarding geothermal drilling, based on worldwide studies [41,42]. For example, the lowest class is for a drilling depth of higher than 4500 meters. The effective geopressured reservoir is the layer with the highest weight in the Index Overlay calculus for the final map with a value of 0.4 by 1.
The concept is properly illustrated in Subsection 6.1, using the 3D geological model of the study area.

Thermal
Regime. The layer of evidence named the thermal regime includes information related to the variation in the thermal state of the underground. In our study, the thermal regime is parameterized by the thermal gradient (see Table 2). Thus, higher thermal gradients imply shallower high temperatures.
The method to compute the thermal gradient is dependent upon the kind of available datasets. The thermal gradient is ranked following the range presented in Table 2, whereas the weight assigned for this layer of evidence is 0.1 by 1.

Pressure Regime.
We propose to produce a map of the pressure gradient through geostatistical analyses of pressure data (e.g., drill steam tests) along wells. The construction of this layer of evidence requires a careful data analysis due to the variability of pressure conditions both vertically and horizontally. For this reason, we rank the gradient computed in the targeted geopressured interval: the larger the overpressure, the higher the favourability.
The limits of the classes (see Table 2) are in part based on the hydrostatic (1st class), soft (3rd), and hard (4th) geopressured regimes as proposed by Loucks et al. [3]. We also used a "light" class (2nd) and the most favourable class (5th) for values larger than 18.82 bar/100 m, i.e., the "near-lithostatic" class. The pressure regime, which plays a very important role in the favourability computation, has a weight of 0.3 by 1.

Deposit
Thickness. The deposit thickness takes into account the role of the compaction disequilibrium for the genesis of overpressure regimes. The importance of a larger thickness of the (effective) reservoir can influence the quality of the resource in terms not only of pressure but also of temperature. In fact, assuming that the thermal gradient is not driven by pure convection, higher temperatures can be reached at larger depths.
This layer of evidence is obtained by simply classifying the thickness of the deposits of the studied basin. The classes (Table 2) were set based on knowledge-driven considerations related to the thicknesses of basin deposits worldwide. The weight is 0.1 by 1.
3.5. Geochemistry. The decision to focus on terrigenous sedimentary basin plays is driven by the possible assumption that their formation waters are saturated in methane. This layer of evidence is aimed at ranking the study area according to the evidences of the occurrence of CH 4 in reservoir. Due to the vertical and horizontal variability of the CH 4 content, the spatial mapping of this layer is extremely difficult.
The score is simply classified according to the occurrence or not of clear indications of CH 4 -saturated and oversaturated water, 5th and 1st classes, respectively (see Table 2). The weight of the layer of evidence is 0.1.

Geological-Structural Setting of Eastern Abruzzo (Italy)
The study area is located in the central-eastern sector of the Italian Apennines (inset in Figure 3(a)), which experienced several deformation events in response to the late Neogene tectonic convergence between the European and African plates. The outcropping tectonic units (Figure 3(a)) derived from the deformation of both shallow-water limestones (carbonate platform domains) and deeper-water carbonate (slope and pelagic basin domains) successions deposited on the southern Neotethyan passive margin ( [43] and references therein). Since the very late Messinian, these successions were affected by eastward-directed fold and thrusting which led, initially, to the growth of the main Apenninic ridges (Gran Sasso and Majella) and later (Late Pliocene-Early Pleistocene) of the Apennines and peri-Adriatic foothills [43][44][45][46][47][48][49][50][51][52].
The contractional deformation [53][54][55] piled the Apennine thrust units onto the Adriatic foreland (Triassic-Early Miocene) consistently with an overall in-sequence regional propagating model [44,56,57] (Figure 3(b)). From the Late Pliocene-Early Pleistocene, the western Abruzzo chain area was affected by an extensional tectonic regime which was responsible for the formation of an articulated system of normal faults [58,59] and associated Quaternary continental basins (Figure 3(a)).
The study area ( Figure 3) also includes the Neogene-Quaternary Abruzzo foredeep and the adjacent Adriatic foreland [60][61][62]. From a structural point of view, the foredeep can be subdivided into two sectors. The western sector is characterized by outcrops of syn-orogenic turbiditic successions of the upper Laga (Late Messinian) and Cellino (Zanclean) Formations. The eastern sector is dominated, at the surface, by late-orogenic shelf clays, evolving upward to coastal sands of the Late Pliocene-Pleistocene age (Mutignano Fm [63,64]).
Active deformation characterizes the eastern foredeep sector, which since early Pleistocene times has been undergoing eastward overthrusting above the Adriatic foreland [65,66]. Instrumental seismicity and borehole breakouts (Figure 3(a)) show a near horizontal WSW-ENE 6 Geofluids  Figure 3: (a) Simplified geological-structural map of the study area with the main lithological units and tectonic structures. All information was derived from geological maps on a 1:100.000 [92] and 1:50.000 [93] scale and from [68]. Tectonic structures were mainly derived from CNR [76]. Key: (1) continental deposits (Quaternary), (2) marly and clayey deposits of the Mutignano Fm and equivalent units (Late Pliocene-Early Pleistocene), (3) arenaceous and pelitic deposits of the Laga (A) and Cellino Fms and equivalent deposits (Late Messinian-Early Pliocene), (4) slope-to-basin and basinal carbonatic deposits (Late Triassic-Miocene), (5) platform carbonatic deposits (Late Triassic-Miocene), (6) undifferentiated deposits pertaining to the allochthonous Molise Nappe (Upper Cretaceous-Upper Miocene), (7) main thrusts ((a) outcropping, (b) inferred or buried), (8) normal faults ((a) outcropping, (b) inferred or buried), (9) boundary of the allochthonous Molise Nappe, CoS = coastal anticline, (10) direction of the minimum horizontal stress (SHmin) referring to a selection of A and B quality borehole breakout data as reported in [94], and (11) P-axes from a compilation of focal mechanisms taken from the RCMT and TDMT catalogues [95] plus other focal solutions deduced from specific papers (e.g., [96]) for the Italian earthquakes with Mw > 4 0 occurring between 1968 and May 2016. BCS = Bellante-Cellino structure. The black dotted box includes the sector investigated for favourability assessment. (b) Interpretative geological section (trace in (a)) showing the main thrusts deforming the carbonatic and foredeep deposits of the outer Abruzzo region. The outcropping successions and their thickness were constrained with information from geological cartography [92,93]. The thicknesses of the Mezo-Cenozoic carbonatic deposits in the Adriatic foreland were extracted from [78], while most of the data concerning the depth of the base of Pliocene deposits come from the present study (see Subsection 5.1). Key: compressional deformation, on average characterized by a low seismic budget (~0.3 mm/year) [67].
The lower Pliocene succession (Cellino Fm and equivalent units) is particularly important for our study, since it is the possible overpressured target. Its stratigraphic range spans from the Sphaeroidinellopsis sp. to the G. puncticulata biozones [64,68]. Its succession is up to 2 km thick and consists of several alternations of poorly cemented arenaceous bodies and thick pelitic units. On the whole, this succession is sandwiched in between the post-evaporitic pelites of the topmost Laga Fm and the thick package of clays of the lower Mutignano Fm.
Thrust and folds, affecting the Laga and Cellino deposits, are well documented in the subsurface of the Abruzzo periadriatic area. Two main buried structural trends are defined in the literature, referred to as the "Bellante-Cellino structure" (BCS) and the "Coastal anticline" (CoS), from west to east (Figure 3(a); [46,[60][61][62]69]). Between the BCS and the CoS, the upper Pliocene sedimentary wedge reaches a maximum thickness of 2000 m and is slightly folded. The Pleistocene succession is only slightly deformed in an E-dipping, low-angle monocline. East of the leading edge of the Apennine outer thrust, the Adriatic foreland is characterized by a west-dipping low-angle regional monocline, affected by minor normal faulting, due to flexural retreat, on which the Plio-Pleistocene siliciclastic succession lies conformably above the Mezo-Cenozoic carbonates [70].
The foredeep-foreland system of the Abruzzo region belongs to the Adriatic petroleum province [19][20][21], where many exploration plays and productive oil and gas fields have been in operation. Different kinds of plays occur both in siliciclastic basinal and carbonate platform systems. Gas fields are present in the Plio-Pleistocene turbiditic sequences in channelized or deep-sea fan deposits [19,69].
With regards to the geothermal resources, the geological conditions of the study area do not favour the development of high temperature systems. Mezo-Cenozoic carbonates represent the regional-scale reservoir for hydrothermal resources. The heat flow map of Italy by Della Vedova et al. [71] shows values mainly in the range of 40-50 mW/m 2 in the study area, with a positive anomaly in the south-eastern offshore part (up to 80 mW/m 2 ). This pattern is expected in areas experiencing a very high sedimentation rate. The geothermal gradient values are in the range of 30-40°C/km regarding the onshore areas, as shown in the Italian geothermal ranking by Cataldi et al. [72].

Source of Data from the Abruzzo Case Study
The favourability computation of the geopressuredgeothermal system of Abruzzo was based on a critical review of a large dataset. The main focus was the analysis of about 200 deep hydrocarbon wells, extracted from the Italian National Geothermal Database (BDNG) [73]. For each well, we gathered, from the BDNG, detailed information sheets related to headings (e.g., well name, coordinates), fluid geochemistry, temperature, formation tests, lithology, stratigraphy and ages, drilling (e.g., deviation), and occurrence of formation fluids (water, hydrocarbons).
We have also checked the original master logs, provided by the Italian Ministry of Economic Development [74], in order to support the additional analysis of selected wells, e.g., homogenization of stratigraphic information, and assessment of the validity of the pressure measurements from formation tests (as described in the next subsections).
Beside the well logs and data from the scientific literature, we carried out a geological modelling and coupled thermo-fluid dynamic numerical simulations in order to provide a reliable temperature distribution at depth. Indeed, we used as input for building the layers of evidence both measured data and modelled data.

Geological Data and Modelling.
We reconstructed the base of the Pliocene foredeep deposits (hereafter BP) in the study area by exploiting and integrating datasets from the literature ( Figure 4) and the abovementioned BDNG. This involved identifying the lowermost position of the base of Pliocene foredeep deposits, disregarding the local interposition of tectonic slices which could have piled up and thickened the foredeep succession. These particular settings were especially investigated in the northern Abruzzo (offshore) and in the southern Abruzzo (onshore) sectors. In the former, the compressive deformation also affected the foredeep deposits (see cross section in Figure 3(b)) from the uppermost Early Pliocene to the Early Pleistocene (moving from west to east). In the latter, minor slices of the Molise Allochthonous units have been sandwiched between lower-middle Pliocene deposits [75]. In both cases, the lowermost contacts (both stratigraphic or tectonic) were considered.
We first extracted information on the depth of BP from the stratigraphy reported in the wells drilled for petroleum exploration [74] and organized in the BDNG (Figure 4(a)). More than 180 wells were found to clearly intercept this level. Most (about 170) reported the stratigraphic contact of the lower Pliocene successions (conformable or unconformable) above terms belonging to the different basinal (northern sector) or slope-to-basin and platform (southern sector) depositional areas.
In our reconstruction, we also considered isobaths of the BP already reported in the literature. Different studies were considered in the north-western and south-eastern offshore areas. In the north-western area, since no well clearly reached the BP in the foredeep depocentral areas, we exclusively adopted the isobaths reported by CNR [76] and shown in Figure 4(b). In the south-eastern offshore areas, we integrated the isobaths with those reported in Argnani et al. [77] ( Figure 4(c)) converting the available isochrones into isobaths using published average velocities for the Pliocene and Pleistocene units (e.g., [44]). A recent interpretation of the M13-14 CROP seismic line [78], which crosses the Adriatic offshore, was also considered (section trace in Figure 4(a)).
We achieved our final target by converting all data into feature datasets (point and lines) in order to process them into a common georeferenced framework in a GIS environment. Finally, we integrated all the collected data and converted them into point features which were interpolated 8 Geofluids through an inverse distance weighting (IDW) method, using local barriers represented by the main tectonic structures reported for the study area and/or hypothesized from lateral discontinuities observed in the in-well stratigraphy. A 3D surface with 1000 m pixel resolution was obtained and is shown in the new map of the base of Pliocene in the Abruzzo outer sector (Figure 4(d)).

Pressure
Data. The pressure conditions in the periadriatic FTB-foredeep-foreland system are well known, above all in the northern sector. Carlin and Dainelli [79] reviewed pressure well data in order to define the various pressure systems in the Adriatic foredeep. Abnormal pressure regimes occur mainly in the Pliocene basin deposits, whereas the underlying carbonates are mainly in hydrostatic conditions. The authors defined three pressure regions in the Plio-Pleistocene succession, from the innermost: (i) inner thrust, (ii) deformation front, and (iii) undeformed foredeep. The compaction disequilibrium is the major cause for overpressures, coupled with a minor contribution of tectonics. The top of the geopressured zone corresponds to the base of the Pleistocene, which led us to target the Pliocene basin deposits as also highlighted by the temperature analysis. Our study area was only partially covered by Carlin and Dainelli [79] in their analysis. Thus, a careful analysis of the pressure well data was performed here in order to estimate the pressure gradients along available wells and to obtain first-order information on the pressure conditions.
The dataset is composed of hundreds of drill stem tests (DSTs) and some repeat formation tests (RFTs). Each test 9 Geofluids is described with additional information such as the measured (MD) and true vertical (TVD) depths, lithology, stratigraphy, type of fluid, measured pressure, and above all the duration of all the phases of the tests. In fact, we considered the stabilized pressures measured after a greater shut-in time than the flowing periods. The tests were grouped according to the geological unit tested belonging to Mezo-Cenozoic carbonates or Plio-Pleistocene basin deposits.
In Figure 5, the complete dataset is summarized in a depth vs pressure plot. The pressure data are plotted and classified on the basis of the geological unit. The analysis highlights mainly hydrostatic pressure conditions for the carbonate basement as well as for the Pleistocene sediments, whereas abnormal pressure regimes occur in the Pliocene deposits, in some cases approaching lithostatic conditions (the epochs refer to the GSA scale v4.0 from [80]).
A pressure gradient was estimated for each well in the targeted Pliocene interval. We firstly selected those wells with more than one pressure value in the Pliocene succession in order to compute the pressure gradient by linear regression. In some sectors, where only one pressure value was available in the Pliocene succession along a well, we extrapolated the pressure at the top of Pliocene from the surface assuming a hydrostatic regime in the overlying Pleistocene deposits. We also considered six pressure profiles from Carlin and Dainelli [79] for a total amount of 24 deep hydrocarbon wells suitable for the study of the pressure regimes in the Pliocene succession. Table 3 summarizes the main information on the pressure data indicating the location of the wells (see the ID) and the pressure regimes in Figure 6.

Thermal Data and Modelling.
In order to evaluate the regional thermal structure of the study area, we solved the mathematical model for geothermal heat and mass transfer.
The bottom-hole temperatures (BHTs), measured in the deep hydrocarbon exploratory wells, were used as control points in order to compare the numerical results with the borehole data (Figure 7(b)). The thermal effects due to the drilling had been previously evaluated in order to correct the BHT data [81,82]. Although the corrected BHTs had a mean error of the order of ±10%, they were suitable for highlighting the main regional thermal structure.
The sets of partial differential equations which describe the principles of conservation of mass, momentum, and energy are approximated through the finite element method in the COMSOL Multiphysics environment [83,84]. The steady-state solution was evaluated within a numerical domain of 130 × 100 × 20 km 3 accounting for the whole regional thermal structure of the Abruzzo-Molise outer sector. The tetrahedral mesh-grid counts >3.6·10 6 nodes and the length of tetrahedron edges vary from a minimum of 100 m to a maximum of 1000 m, allowing the numerical mesh to fit the spatial variations of the boundary surfaces. The regional fluid circulation occurs within the permeable Mezo-Cenozoic carbonate units, and in the eastern Abruzzo geological framework, the Plio-Pleistocene foredeep deposits act as the cap-rock of the deep-seated low-temperature hydrothermal system, almost throughout the entire investigated area. The carbonate units outcropping in the Apennine chain represent the main recharge areas of the regional reservoir.
The geometrical model considers three main lithothermal units, from the top to the bottom: (i) the impervious sedimentary cover unit acting as a cap-rock, (ii) the tectonically thickened carbonate units hosting the main regional reservoir, and (iii) the basement unit. Each lithothermal unit is composed of different rocks with similar thermal and hydraulic properties. The rocks were treated as a homogeneous and downward anisotropic porous material. Mixing laws were applied to estimate the effective thermal and hydraulic properties of the rock-fluid system accounting for the in situ conditions [85]. In the cover and basement domains, a purely conductive heat transport takes place, while in the carbonate reservoir, the regional hydraulic gradient and the thermal convection affect the temperature distribution by mass and energy transport. The buried upper and basal boundaries of the reservoir are set impermeable to fluid flow, allowing only conductive heat transfer.
Accounting for the recharge zones, we applied a stress boundary condition where the reservoir units crop out (Figure 7(a)). As the reservoir is assumed to be fully saturated, the pressure on those boundaries is set as equal to  [73,74] in the studied area. Table  3: Pressure measurements in wells. The ID refers to the location in Figure 6. The type of pressure test is reported plus information on measured depth (MD) lithology, age, and the geological unit.     Table 3, which reports the main information (from [73,74]), whereas the wells marked with B are from Carlin and Dainelli [79]. The area in red is the study area for the favourability assessment. A good fit between the measured and simulated temperatures was achieved for a basal temperature of 350°C and an average permeability of the regional carbonate reservoir of 1.6·10 -16 m 2 .

ID
The observed geothermal gradients spanned from 20 to 50°C/km, and the mean value was 28°C/km. The highest and lowest values were observed in the mountainous zones in which the heat advection, controlled by the regional groundwater flow, modifies the conductive thermal structure. Regarding the thermal structure of the study area, we computed the depth (a.s.l.) of the 90°C isotherm (also referred to as z90). The isobaths show two minima (z90 > 5 km) in the south-western corner separated by an upwelling zone (z90 ≈ 2 km). In the mainland, far from the carbonate outcrops, z90 settles around 2.5-2.7 km and deepens toward the offshore sector to 3-3.4 km.
On the basis of the thermal and geological models, the sectors in which z90 is above the base of the Pliocene deposits represent the target areas (Figure 7(a)).

Geochemistry.
In our case study, the geochemical analyses on fluids sampled in wells were used to analyse known methane-prospective areas (corresponding to tectonic elements), whereas a geostatistical analysis was not possible.
Minissale et al. [86,87] clearly defined multiple sources of gas discharge manifestations in central Italy. Along a transect crossing the Northern Apennine, the inner part of the chain is mostly characterized by CO 2 -rich emissions, whereas CH 4 -rich emissions are predominant in the outer part of the Apennines along the Adriatic coast, where the study area is located. The source of CH 4 in the Adriatic sector is related to the sedimentary process in the Periadriatic foredeep. In fact, Italy is considered a gas (methane)-prone area, the formation and accumulation of which are strictly related to the tectono-sedimentary evolution and specific areas.
Mattavelli and Novelli [88] identified a narrow area, corresponding to the foredeep of the Apennine belt, where 77% of Italian natural gases were discovered (Figure 8). Here, the synsedimentary tectonics, the highly efficient turbidite systems coupled with subsidence, favoured the generation (and also trapping) of natural gases. Eighty-two per cent of natural gas in this area is biogenic, characterized by almost pure and isotopically light CH 4 with C 2 (ethane) and other components (<0.5%). The authors found strong evidence for in situ generation of gas, through bacterial or diagenetic processes. The gaseous hydrocarbons are mainly stored in Plio-Pleistocene sediments.
We focused on the gas chemistry of sampled fluids in wells. Several chemical analyses of major ions were available for the water samples with complementary notes on the  Figure 8: Location of wells with chemical analyses available on gas and with traces of exsolved gas (data from [74]). The oil and gas leases are from [74,97], whereas the area with a high occurrence of hydrocarbons is from Mattavelli and Novelli [88]. The data shown are only for Abruzzo. The area in red is the study area for the favourability assessment. occurrence of gas (generic) or smell of hydrocarbons. A few were coupled with chemical analyses of well gases.
The analyses of well gases sampled in Plio-Pleistocene deposits show an amount of CH 4 that is always higher than 97% and CO 2 lower than 1%. In addition, tens of wells record exsolved gas traces in the Plio-Pleistocene interval. The same analyses in the underlying carbonates show lower and highly variable amounts of CH 4 . These gases were sampled in relation to supposed hydrocarbon targets in the carbonates. Figure 8 shows the wells with measured CH 4 , as well as those with traces of exsolved gas.
An example is given by the well "Montebello di Bertona 1" in which well gases were analysed both in terrigenous sediments and in carbonate successions. The gas sampled in the interbedded shale and sand, Pliocene in age, is composed of 99% of CH 4 , whereas in the underlying carbonates, CH 4 abruptly decreases and CO 2 increases to over 50%.
The analyses led us to assume formation water saturated and oversaturated in methane, with a clear indication of the diffuse presence of CH 4 in the Plio-Pleistocene deposits, which is locally accumulated. Table 4 summarizes the well gas analyses. We also computed the methane solubility in water by using the empirical method by Price et al. [36] and taking as an example the conditions at depth along the well "Morro D'oro 1" from a test in water with dissolved gas. Considering a formation pressure of 660 bar at 3733 m, a temperature of 90°C, and water salinity of 28.6 g/l, the solubility of methane is 28.75 SCF/Bbl. The solubility we computed is close to the values reported for the geopressured system of the Gulf of Mexico [27,36].

Layers of Evidence
With regard to the Abruzzo case study, in the first stage, a preliminary analysis was carried out in order to identify the areas to be assessed, those strictly related to the concept of sedimentary basin plays. The study area corresponds to the foredeep-foreland domains of the Apennines thrust belt in Abruzzo. The targets are the geopressured resources hosted in the Plio-Pleistocene siliciclastic succession. Basically, the sectors in which the depth of the 90°C isotherm is above the base of the Pliocene deposits represent the target areas.
We then created the layers of evidence by combining different thematic inputs throughout the spatial analyses. These layers of evidence were related to the favourable factors and referred to (see details in Section 3): (i) the effective geopressured-geothermal reservoir, (ii) the thermal regime, (iii) the pressure regime, (iv) the deposit thickness, and the (v) geochemistry.
The five layers of evidence were combined by Index Overlay to produce the final favourability map as a result of the last stage of the workflow. The layers of evidence have been computed with the same spatial resolution, 1 × 1 km, and with the grid nodes of each layer overlapping.
6.1. Effective Geopressured-Geothermal Reservoir. The concept of the effective geopressured-geothermal reservoir is illustrated in Figure 9, using an example from the 3D model of the study area. In this paper, a threshold value of 90°C was set considering the geothermal power production purposes. The layer was built by applying a layer intersection between the depth of the 90°C isotherm and the bottom surface of Table 4: Well gas chemical analyses, in relation to the CH 4 and CO 2 percentage. The ID refers to the location in Figure 8.  16 Geofluids the Pliocene deposits. The result is a grid layer of a ranked depth of the top of the effective geopressured-geothermal reservoir ( Figure 10). Where the isotherm is deeper than the base of the Pliocene deposits, i.e., no effective reservoir occurs, the corresponding areas are neglected from the computation.
The ranking classes and the weight are listed in Table 2.
6.2. Thermal Regime. In our study, the thermal regime is parameterized by the thermal gradient. The thermal gradient was computed starting from the average air temperature, which varies on the basis of the topography elevation, and using the layer of the temperature at the top of the carbonate formations. By a raster computation, we obtained an average air temperature for each cell from the digital elevation model (DEM) layer applying a lapse rate of 0.0065 (°C/m) and 0.0045 (°C/m) for positive (i.e., above sea level) and negative (i.e., below sea level) elevations, respectively [89,90]. A further raster computation then enabled us to define the thermal gradient map (Figure 11). The thermal regime was ranked and weighted as shown in Table 2. 6.3. Pressure Regime. We obtained information of the pressure regions through geostatistical analyses on the pressure gradients computed in the wells reaching the Pliocene succession. A statistical analysis showed only a partial trend along the X and Y coordinates. After various attempts, we have chosen the resulting interpolation obtained with the Universal Kriging algorithm, which guaranteed the least root-mean-square error. The resulting grid layer is shown in Figure 12.
The rank is related to the amount of overpressure that occurs in the Pliocene sedimentary succession.

Deposit
Thickness. This layer of evidence was obtained by simply classifying the thickness of the basin deposits (i.e., the Pliocene bottom depth from the ground level) using the raster recode function (Figure 13). The classes and weights are listed in Table 2. 6.5. Geochemistry. This layer of evidence, for the Abruzzo case study, is essentially based on the ranking of the methane-prospective area, corresponding to the foredeep domain, proposed by Mattavelli and Novelli [88] ( Figure 14). We assigned the highest class (5th) to this area, assuming the occurrence of CH 4 -saturated and oversaturated waters in reservoir. The available geochemical dataset on fluids sampled in wells supports this assumption (see Section 5.4). The dataset was not complete enough to compute geostatistical analysis. The remaining part of the area, not included in this CH 4 -prospective domain, was ranked with the less favourable class (1st).

The Favourability Map of Abruzzo: Results and Discussion
The quantitative integration of data using the Index Overlay method (see Section 6) resulted in the favourability map of a geopressured-geothermal system for the foredeep-foreland basin play of Abruzzo, shown in Figure 15. The role of the prospective factors at regional scale (Table 1) has been addressed in the five layers of evidence. The geopressured effective reservoir takes into account the depositional environment, the temperature, the depth of the reservoir, and intrinsically the lithology. Pressure and temperature have been considered in the pressure and thermal regime layers of evidence, respectively. The deposit thickness takes into account the depositional environment, the deposition rate, the size of the reservoir, and the lithology. The geochemistry layer provides information on the occurrence of dissolved gaseous hydrocarbons.
A point that deserves a clarification is the use or not of the faults as permeability indicator. The main limit of its use in the favourability computation is the possible hydraulic behaviour of the faulted volume of rocks. A fault can act as a permeability barrier or can enhance the permeability. At regional scale, it is not accurate and reliable to assign a favourable value around a fault trace. At local (i.e., exploration licence) or well scale, the detailed geophysical and geological studies allow the understanding of the hydraulic feature of a fault. We decided to disregard this contribution in our computation at regional scale, but we consider the structural setting a key information to be addressed during a local exploration project.
The study area was mostly ranked from not to low/medium favourable, with one exception corresponding to a wide continuous prospective sector in the centre. The most favourable sector, with a rank up to the 4th class, extends for less than 1000 km 2 and runs parallel to the shoreline along a  18 Geofluids NW-SE direction, both offshore and onshore. The 5th class (very favourable) was not retrieved. The cells of the grid that have been ranked (from 1st to 5th class) are those where the effective reservoir was detected; otherwise, the cells were not considered favourable. We assigned the highest weight in the Index Overlay computation to the effective reservoir layer of evidence (0.4; see Table 2) because it takes into account (i) the minimum temperature requirement for the exploitation, (ii) the financial factors related to the depth of the drilling target, and (iii) the occurrence of a potential reservoir. Since the subsurface temperature is usually close or lower than the average continental reference in the study area, the depth of the top of the effective reservoir belongs mostly to the medium and low favourable classes with scores of 3 and 2, respectively. This means that in the best scenario, wells should reach a 2.5-3 km depth to exploit 90°C.
The second most weighted layer of evidence (0.3; see Table 2) is the pressure regime. Our study area is mostly characterized by high to very high favourable classes (4 and 5 scores) with hard geopressured gradients in the Pliocene succession, locally approaching the lithostatic reference. We consider this layer of evidence as a useful first-order approximation of pressure regions. However, due to the low number of pressure data, the spatial analysis cannot be used to evaluate a precise value at depth of the formation pressure. As expected, the most abnormal pressure regime occurs along the deformation front subparallel to the coast. The inner sector (westwards) shows hydrostatic to soft geopressured regimes, whereas the undeformed foredeep (eastwards) shows soft to hard geopressured regimes. The results of our analysis are in good agreement with the pressure regions identified by Carlin and Dainelli [79] for the northern sectors of the Adriatic Sea.

Geofluids
Although temperature is embedded in the concept of the effective reservoir, which considers the 90°C isotherm, we included a layer of evidence exclusively related to the thermal regime, with a weight of 0.1 by 1. We thus included information on the occurrence of positive thermal anomalies. As expected, most of the study area shows low values of geothermal gradient except for the southern part, where higher values belonging to the 3rd class were found. Here, the deep water rising toward the groundwater base level can affect the conductive thermal gradient within the impermeable cover.
The layer of evidence related to the deposit thickness shows higher values, up to the 4th class, only in a small sector located in the deformation front along the northern coastal zone. The weight assigned is the same as the thermal regime, i.e., 0.1 by 1. We chose a weight value lower than the effective reservoir and pressure regime because this layer is mainly related to geological conditions that can indirectly favour the presence of geopressured resources and less direct information is provided, such as the possible extent at depth of reservoirs. We stress that the deposit thickness is correlated with the deposition rate, but it is also controlled by tectonics in the study area.
The aim of this study was to assess the geopressured resources occurring in terrigenous sedimentary basin plays, as they favour the generation of overpressure and above all widespread dissolved gaseous hydrocarbons. The importance of methane is already embedded in this focus. Furthermore, we assumed that there are methane-saturated waters in the CH 4 -prospective area [88] corresponding to the foredeep domain, after a duly analysis, assigning the highest class (5th). In our case study, the geochemical dataset was too poor to rank the area accordingly and the assigned weight is 0.1 by 1.

Geofluids
The procedure can be improved with a quantitative ranking of the layer of evidence regarding the geochemistry (e.g., percentage of methane in free gases, methane solubility) by using geostatistical tools.
From a practical point of view, the map can be considered constrained where direct data along deep wells are available. For this reason, we provide the final favourability map coupled with the well locations. Moreover, to further assess and exploit the geopressured resources in the study area, it should be considered that the favourable part located in the offshore is limited by technical (even if the sea floor of the Adriatic Sea is really shallow) and nontechnical (e.g., legal regulations) barriers. ii) The weight of the first layer of evidence (W i=1 ) was perturbed in the range 0.1-0.9 while homogenously distributing the remaining weight to the other layers (being the sum of the weights equal to 1). The procedure is repeated for each layer of evidence. For each scenario, a favourability map (out n ) was computed by IO (see (1)). We stress that the classes of each layer of evidence are kept equal to the reference scenario (iii) The percentage change in the output (Δout) of each scenario is computed for each pixel with respect to the reference map (v) The sensitivity s for each scenario is computed pixel by pixel as the ratio between percentage changes in the output and weights, as follows: The results are summarized in Figure 16. We show the sensitivity maps evaluated for each layer of evidence with weights varying from 0.1 to 0.9 with a step of 0.2.
The comparison among the various sensitivity maps highlights a larger sensitivity to the geopressured effective reservoir and the pressure regime layers of evidence. The main variations of the sensitivity index mimic the spatial distribution of the pressure regime, highlighting a very important role of this layer. This behaviour is expected because among the most weighted layers, the pressure regime shows the largest spatial variability; on the contrary, the effective reservoir is almost homogenous.
Conversely, our reference map is slightly sensitive to the perturbation on the thermal regime, geochemistry, and the deposits thickness layers, except for a few restricted areas. Particularly, the sector in the southwestern corner of the map shows larger sensitivity. This response corresponds to an area at the contact of outcropping carbonates where very high thermal gradients are modelled. This computation appears to be a local artefact, due to the fact that the thickness of the impermeable cover is really short and the resulting thermal gradients were extremely high.

Conclusions
This paper is intended to be a practical analytical framework for the systematic integration of the relevant data required to  Figure 14: Geochemistry layer of evidence. 22 Geofluids assess the geopressured-geothermal resources. The approach described can be considered valid and applicable at a global scale as the whole procedure is based on generic features of terrigenous sedimentary basins and was further tested on the Abruzzo case study.
The main results of this study can be summarized as follows: (i) We propose a novel integrated methodology aimed at assessing and mapping the favourability of geopressured resources in sedimentary basin plays worldwide. We used a consolidated system of weighting and scoring in a GIS environment, the Index Overlay, to analyse various prospective parameters.
(ii) We identified five layers of evidence that describe the prospective factors needed for assessing such unconventional resources: (i) the effective geopressured-geothermal reservoir, (ii) the thermal regime, (iii) the pressure regime, (iv) the deposit thickness, and (v) the geochemistry. The effective geopressured-geothermal reservoir is a concept that can be of help for the industrial exploration. The depth of its top is a key parameter, easy to map.
(iii) Regarding the Abruzzo case study, an updated dataset of subsurface data was organized and important outcomes are now available such as the 3D geological model, including an up-to-date base of Pliocene deposits, and the study on the pressure and thermal regimes.
(iv) The final favourability map for the Abruzzo case study is a first attempt at ranking these kinds of unconventional geothermal resources in a region that has been historically explored and exploited  Figure 15: Favourability map of the geopressured-geothermal system for Abruzzo. 23 Geofluids only for hydrocarbons. In order to provide a quantitative estimation of the resources required for industrial projects, detailed exploration activities at the local scale (research permits) are required.

Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.   Figure 16: Sensitivity analysis for the favourability map of the Abruzzo region. The outcomes refer to the computation of the sensitivity by perturbing weights of a selected layer of evidence while distributing homogeneously the remaining part of the weight to the other layers (by sum equal to 1). W 1 = 0 1, W 2 = 0 3, W 3 = 0 5, W 4 = 0 7, and W 5 = 0 9. The contour of the value ±0.15 is shown. The missing outcomes are those with the perturbed weight of the selected layer of evidence equal to the reference weight. 24 Geofluids