# Prospectivity analysis for underground hydrogen storage, Taranaki basin, Aotearoa New Zealand: A multi-criteria decision-making approach

**Authors:** Karen E. Higgs, Dominic P. Strogen, Andrew Nicol, David Dempsey, Kerry Leith, Kari Bassett, Catherine Reid, Edward Yates, Matthew Parker, Alan Bischoff, Ludmila Adam, Michael Rowe  
**Published in:** International Journal of Hydrogen Energy (2024)  
**DOI (version of record):** https://doi.org/10.1016/j.ijhydene.2024.05.098  
**This copy:** Published version. Text extracted from the PDF at https://ddempsey.github.io/papers/2024-higgs-prospectivity-analysis-underground-hydrogen-storage/higgs2024-prospectivity-analysis-underground.pdf; equations, tables and figure placement may be imperfect, so cite the DOI.  
**Notice:** Published version, reproduced under the CC BY 4.0 licence. © The Author(s). Version of record: https://doi.org/10.1016/j.ijhydene.2024.05.098  
**Licence:** CC BY 4.0 http://creativecommons.org/licenses/by/4.0/  

---
International Journal of Hydrogen Energy 71 (2024) 1468–1485 

Contents lists available at ScienceDirect 

International Journal of Hydrogen Energy 

journal homepage: www.elsevier.com/locate/he 

Prospectivity analysis for underground hydrogen storage, Taranaki basin, Aotearoa New Zealand: A multi-criteria decision-making approach 

Karen E. Higgs<sup>a,b,c,*</sup> , Dominic P. Strogen<sup>c</sup> , Andrew Nicol<sup>b</sup> , David Dempsey<sup>b</sup> , Kerry Leith<sup>c</sup> , Kari Bassett<sup>b</sup> , Catherine Reid<sup>b</sup> , Edward Yates<sup>b</sup> , Matthew Parker<sup>b</sup> , Alan Bischoff<sup>d</sup> , Ludmila Adam<sup>e</sup> , Michael Rowe<sup>e</sup> 

a _Independent, The Old Church Rooms, Bryn-y-maen, Colwyn Bay, United Kingdom_ b _University of Canterbury, Private Bag 4800, Canterbury, 8041, New Zealand_ c _GNS Science, 1 Fairway Drive, Avalon, Lower Hutt, 5040, New Zealand_ d _Geological Survey of Finland, Vuorimiehentie 5, P.O. Box 96, FI-02151, Espoo, Finland_ e _University of Auckland, 20 Symonds Street, Auckland, 1010, New Zealand_ 

A R T I C L E I N F O 

A B S T R A C T 

Handling Editor: Dr Mehran Rezaei 

_Keywords:_ Geological hydrogen storage Renewable energy storage Depleted hydrocarbon fields Porous media Site evaluation 

Seasonal underground hydrogen storage (UHS) in porous media provides an as yet untested method for storing surplus renewable energy and balancing our energy demands. This study investigates the technical suitability for UHS in depleted hydrocarbon fields and one deep aquifer site in Taranaki Basin, Aotearoa New Zealand. Prospective sites are assessed using a decision tree approach, providing a “fast-track” method for identifying potential sites, and a decision matrix approach for ranking optimal sites. Based on expert elicitation, the most important factors to consider are storage capacity, reservoir depth, and parameters that affect hydrogen injectivity/withdrawal and containment. Results from both approaches suggest that Paleogene reservoirs from gas (or gas cap) fields provide the best option for demonstrating UHS in Aotearoa New Zealand, and that the country’s projected 2050 hydrogen storage demand could be exceeded by developing one or two high ranking sites. Lower priority is assigned to heterolithic and typically finer grained, labile and, clay-rich Miocene oil reservoirs, and to deep aquifers that have no proven hydrocarbon containment. 

# **1. Introduction** 

Energy production using hydrogen gas (H2) for combustion and electrochemical fuel cells emits no greenhouse gases and, as such, is expected to play an important role in a future global carbon-neutral economy [1–4]. While hydrogen can be produced by different industrial processes, it is water electrolysis that produces “green hydrogen” by using renewable electricity [5–7]. The temporary storage of excess renewable energy, as green hydrogen, has been proposed to mitigate the fluctuating supply of electricity from renewable sources, while balancing energy supply and demand [2,8–12]. Wind and solar power are two renewable sources set to provide significant energy in the future, yet their output is heavily dependent on short-term climatic fluctuations. Consequently, underground hydrogen storage (UHS) is considered an attractive large-scale energy storage option for renewable spill [2,3,5, 7,12]. 

The concept of UHS is to convert surplus green electricity into H2 and inject it underground into porous geological formations (such as depleted hydrocarbon reservoirs or deep aquifers) or salt caverns/ engineered cavities. Caverns and cavities may provide a relatively small, peak load storage solution (e.g., 10s–250 GWh range [7,13]), while UHS in geological formations can provide sufficient volume to fulfil seasonal storage requirements (e.g., TWh range [14–16]). To maximise the yield of renewable energy systems, H2 would be stored during periods of high surplus and low energy demand, and then be extracted as dictated by increasing demand, thus supporting management of the energy system. 

UHS in porous formations requires a suitable geological structure in the form of a trap (structural, stratigraphic or combination), a porous and permeable reservoir to store H2, and an impermeable caprock (top seal) to prevent upwards leakage. In addition, UHS requires a cushion gas, which would be stored below the working gas, primarily to maintain the required operating pressure and deliverability rate [9,11,17,18] 

- Corresponding author. Independent, The Old Church Rooms, Bryn-y-maen, Colwyn Bay, United Kingdom. _E-mail address:_ k.higgs@reservoirgeosolutions.co.uk (K.E. Higgs). 

https://doi.org/10.1016/j.ijhydene.2024.05.098 

Received 20 December 2023; Received in revised form 8 March 2024; Accepted 7 May 2024 

Available online 30 May 2024 

0360-3199/© 2024 The Authors. Published by Elsevier Ltd on behalf of Hydrogen Energy Publications LLC. This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/). 

_International Journal of Hydrogen Energy 71 (2024) 1468–1485_ 

_K.E. Higgs et al._ 

**Fig. 1.** Simplified UHS scenario. H2 is injected and stored as the working gas zone in a permeable reservoir, trapped in a structural closure, and overlain by an impermeable seal (caprock). Injected H2 will displace _in situ_ pore fluid due to its lower density. The producible working gas zone is underlain by cushion gas, which will undergo alternate compression and expansion during injection and withdrawal cycles respectively to maintain pressure and deliverability. 

(Fig. 1). The presence of a cushion gas may also help to prevent water from entering the storage compartment [9] and minimise H2 losses [18]. Various gases can serve as the cushion gas, with a concentration gradient expected between H2 and a denser gas (e.g., N2 or CH4), and the degree of mixing from advection, diffusion, buoyancy and other mobilization mechanisms dependent on cushion gas composition and site-specific storage operations (e.g., injection/withdrawal rate, reservoir properties etc.) [6,19]. 

The success of UHS will be dependent on locating suitable sites that satisfy a wide range of technical (geological) requirements, together with sociocultural, economic, political, and regulatory criteria [3,5,10, 20,21]. To date, there are no operating industrial-scale UHS sites for high-purity H2 in porous media, hence there are few established regulations and only a theoretical understanding of the social and economic impacts. Additionally, little is known about how H2 will behave in the subsurface [22–24], and much more fundamental work on hydrogen systems is required to evaluate UHS feasibility [12,23]. 

The objective of our study is to undertake prospectivity analysis for UHS in Aotearoa New Zealand (A-NZ) as preparatory work to further UHS characterisation and research. First we introduce some previous work on underground gas storage (section 2), highlighting specific aspects and/or issues related to UHS in other parts of the world. The potential for UHS in Taranaki Basin (A-NZ) is then briefly described (section 3) and includes a summary of the A-NZ energy market and regional geology. This is followed by a description of methodologies (section 4) and results (section 5) for two different multi-criteria decision-making approaches to site evaluation. A comparison of the two approaches and summary of outcomes is presented in the discussion (sections 6.1 and 6.2 respectively). Our study provides an example for ranking UHS prospects based on current knowledge and using a comprehensive set of technical parameters. However, there remain some outstanding technical uncertainties with UHS technology, and further scientific research will be necessary to fully understand the requirements for successful deployment of UHS. Technical research directions for A-NZ that address some of the uncertainties are presented in the discussion section (section 6.3). Enviro-socio-cultural impacts unique to A-NZ are being investigated in parallel, and together will help 

generate a workflow for optimal UHS site selection that screens and ranks potential sites and provides an input to risk assessment. 

# **2. Previous work** 

Many studies have been undertaken in recent years that discuss and characterise the potential for underground storage of gases in subsurface geological formations. The oil and gas industry has decades-long experience with underground natural gas storage (UGS) to buffer variable demands, many of which focus on depleted hydrocarbon reservoirs [3,4 21]. More recent work has explored long-term (permanent) underground storage of CO2, which has the potential to play a significant role in the mitigation of global greenhouse gas emissions through carbon capture and storage (CCS). Both UGS and CCS studies include site screening, evaluation of storage capacity, best practices, monitoring, potential reactivity, and provide valuable insights for site evaluation in a general sense [25,26]. However, the storage of H2 requires specific geological and operational conditions [9,23,27] and a more tailored approach must be developed for assessing UHS sites that consider the physico-chemical composition of H2, and H2 storage goals and processes. 

The importance of hydrogen in the energy transition has gained significant attention in the last decade. UHS facilities have been developed in salt caverns to contain H2 for industrial uses, and this technology is now well established [9,13,14,28]. However, experience with UHS in porous geological formations is limited to storage of gas mixtures (town gas) [6,11] and pilot projects in Argentina and Austria [11]. Preliminary UHS screening studies have been undertaken in several countries, including Australia [22], Brazil [29], Canada [28,30], China [31], Germany [32], Japan [33], Netherlands [34], Poland [20,35], Portugal [36], Spain [8], United Kingdom [4,9,37,38], and USA [14], and two larger projects (HyUSpre and Hystories) have recently completed an evaluation of storage capacities for porous reservoirs in multiple European countries [39–41]. The concept of UHS in Aotearoa New Zealand has also been considered, both from a geological perspective [24,42], and including cultural, environmental, commercial, and infrastructural aspects [42]. 

The feasibility of UHS in terms of hydrogen injection/withdrawal, transportation, storage, and economic viability has been the focus of several studies [8,14,43–46]. There are also a number of recent publications that provide comprehensive reviews on the technical requirements for UHS and/or potential issues that still need to be resolved [2,3,6,11,12,22,23]. These reviews have identified some key challenges for UHS as: the potential for H2 leakage through faults and/or top seal; potential contamination and/or changes to rock properties related to geochemical reactions and microbial growth; operational issues, including the interaction of stored H2 with residual fluids and/or cushion gas under reservoir conditions. 

Recent investigations into H2 flow dynamics and properties look at the behaviour of H2 in response to its physico-chemical properties, such as high diffusivity, low viscosity, and high mobility [21–23]. These investigations comprise numerical simulations of UHS systems, looking at plume migration, injection/withdrawal cycles, recovery factors, and interphase interaction [18,19,39,40,47–53]. Input data for modelling rock/H2/brine/petroleum systems is limited (e.g., wettability, interfacial tension, capillary pressure data) but has recently been gaining attention ([27,54], and references therein). These studies generally support the feasibility of UHS, and suggest minimal H2 losses through caprocks of depleted gas reservoirs [52]. 

Significant work has also recently focused on the potential for chemical reactions [32,55–62] and microbial interactions [39,40,63,64] within UHS sites. While reactivity of silicates with H2 is considered to be minimal over the timeframes for storage (months to years), there are several studies to suggest reactions are likely between H2 and carbonates, sulphates and sulphides. These reactions could result in changes to petrophysical rock properties of reservoir or seal [32,57,65–69], and 

1469 

_International Journal of Hydrogen Energy 71 (2024) 1468–1485_ 

_K.E. Higgs et al._ 

with the formation of methane and/or hydrogen sulphide contaminating the stored gas [6,40,55]. Microbial growth is expected to be significant for low-salinity, low-temperature reservoirs, which could significantly impact UHS through H2 losses, clogging and corrosion [64]. 

At present there is very limited published work that focuses on selecting suitable sites for UHS [35,39–41,53,70,71]. Studies with weighted criteria suggest that of their assessed parameters, the most important are flow capacity and reservoir depth [71], overburden lithology [35], maximum well deliverability rate and cushion gas requirement [70], or readiness of the development (time-to-market) and risk induced by microbiological activity [40]. Research to date suggests that the success of UHS will be dependent on site-specific geological and operational conditions based on reliable site characterisation, and on identified technical parameter requirements specific to UHS. The reader is referred to supplementary data file 1, which summarises our current understanding of the main technical considerations for UHS based on many of the recent studies. 

# **3. Potential for UHS, Taranaki basin, Aotearoa New Zealand** 

# _3.1. Aotearoa New Zealand energy market_ 

A-NZ’s energy market provides a small economy by global standards. Renewables, including hydro, geothermal and increasingly wind, dominate the electricity generation portfolio, with 87.1% of total generation from renewable sources in 2022 [72]. However, seasonal shortfall of electricity in A-NZ is currently accommodated by natural gas and/or coal, and other sectors are also largely reliant on hydrocarbons (e.g., industry, agriculture and transport) [72]. Data from the Ministry of Business, Innovation and Employment (MBIE) shows that total energy demand for 2022 was 150 TWh (543 PJ) of which 100 TWh (360 PJ) was sourced from hydrocarbons (oil, natural gas and coal). In a bid to decarbonise the economy the New Zealand Government has set a target of a 100% renewable electricity sector by 2030 and committed to achieving net-zero by 2050 [42,73]. Green hydrogen is being investigated as a replacement for hydrocarbons through its use in downstream industries (e.g., fertilizer production at Kapuni) and transport (e.g., Hiringa Energy hydrogen refuelling network) [24,42]. Hydrogen also has potential in A-NZ for export (e.g., methanol or urea), and for mitigating the fluctuating energy demand by low-carbon energy storage/withdrawal (e.g., UHS). 

Yergin [74] suggested that hydrogen could account for at least 10% of the global energy system by 2050. For A-NZ, this would amount to ~15 TWh of hydrogen according to present energy demand, likely increasing by 2050. In A-NZ’s North Island, there is a natural gas pipeline network and one underground storage facility, Ahuroa, with capacity for about 10% of annual gas use. The network owner, FirstGas Group, have proposed that hydrogen could replace natural gas in most sectors of the country by 2050, and help to decarbonise other parts of the economy, such as transport [75]. FirstGas project the hydrogen demand in 2050 to be 42 TWh and have announced plans to decarbonise the gas pipeline network with the aim of blending green gas (including H2 and biogas) into the natural gas network from 2030 [24,75]. If we assume that the 2050 annual hydrogen demand for A-NZ will be in the range 15–42 TWh, and that storage requirements are similar to the present, that is, 10% of the total energy demand, then this amounts to 1.5 to 4.2 TWh of useable UHS (not including cushion gas requirements). 

Salt caverns are currently the only proven method for storing and recovering high-purity H2, yet salt deposits do not occur in A-NZ. There are, however, numerous Cretaceous–Cenozoic sedimentary basins, and UHS in porous formations is considered the most viable large-scale hydrogen storage option for the country [24]. Taranaki Basin is the only sedimentary basin in A-NZ with commercial petroleum production, and has been extensively studied, providing numerous open-file datasets that could be utilised in the investigation of UHS (e.g., core, seismic, well logs, production and pressure data etc.). It therefore represents an 

excellent region to study the technical requirements for UHS and to develop best practices for geostorage prospectivity analysis. The existing petroleum infrastructure and growing renewable (notably wind) electricity resources makes Taranaki well positioned to emerge as a major hydrogen hub in A-NZ and globally. 

# _3.2. Taranaki basin geology_ 

Taranaki Basin contains up to 11 km of mid-Cretaceous–Cenozoic strata and is bound along its eastern margin by the buried, crustal-scale Taranaki Fault System [76,77] (Fig. 2A). Basin development was initiated in the mid-Cretaceous associated with rifting prior to breakup of the eastern margin of Gondwana and formation of the Tasman Sea [76,78, 79]. Deposition of syn-rift sediments was followed by a prolonged period of relative tectonic quiescence, with a post-rift passive margin persisting through the late Paleocene–Eocene [76,80]. Subduction and arc magmatism started at the end of the Eocene, and contraction accelerated through the Miocene creating trapping folds in the basin [81,82]. 

The eastern and southern parts of Taranaki Basin are within the active plate boundary deformation zone, and have been affected by complex overprinting of extensional and contractional tectonic processes [76] (Eastern Mobile Belt on Fig. 2A). Many of the anticlinal structures that have been targeted for petroleum exploration occur in this region, and include the Tarata Thrust Zone beneath eastern Taranaki Peninsula [76,77]. An active volcano (Mt Taranaki on Fig. 2A) is located on the Taranaki Peninsula, west of the Tarata Thrust Zone. By comparison, the western parts of Taranaki Basin have remained relatively undeformed since the latest Paleocene [79] (Western Stable Platform on Fig. 2A), and deposition is largely characterised by progradational sedimentation on a regionally subsiding seafloor [76]. 

The stratigraphy of Taranaki Basin is illustrated in Fig. 2B, which shows the distribution of reservoir and sealing formations relative to basin history. Proven commercial reservoir intervals range from Paleocene to Pliocene age and consist of a wide range of sandstone facies (fluvial, marginal–shallow marine, and deep marine), characterised by variable reservoir mineralogy, quality, and diagenetic alteration. Sandstone composition and quality is strongly related to sediment provenance and basin history, with older Paleogene reservoirs generally coarser-grained, and comprising a relatively mature quartzo-feldspathic composition compared to the younger, finer-grained and immature, active margin, Neogene reservoirs [83]. The Tikorangi Limestone is the only proven carbonate reservoir in the Taranaki Basin, and is characterised by extremely poor matrix permeability with reservoir quality reliant on fracture porosity. Proven sealing caprocks are primarily Paleocene–Eocene offshore silty mudstones, Oligocene deep marine calcareous mudstones/limestones (where unfractured), and Miocene–Pliocene deep marine mudstones [76]. 

# _3.3. UHS options in Taranaki basin_ 

# _3.3.1. Depleted oil and gas reservoirs_ 

Depleted oil and gas reservoirs are porous formations where the pore space contains residual hydrocarbons as well as an aqueous phase, following production. Depleted reservoirs generally have well-identified geological structures, proven reservoir quality, a proven caprock, and existing infrastructure, and therefore provide good options for upscaling hydrogen geostorage to an industrial scale. These reservoirs are usually supported with an aquifer at the bottom or side(s), from which the pore volume may have been variably filled, depending on the drive mechanism for the particular reservoir. Residual gas within a depleted gas reservoir might be beneficial for UHS if it provides pressure support as the cushion gas [4,6,37], or a disadvantage if mixing with the remaining natural gas reduces the H2 purity [11,51]. Generally, residual liquid hydrocarbons (oil/condensate) in a potential UHS site is considered to be more problematic than natural gas due to low recovery factors, the possibility of chemical reactions, dissolution of H2, and reduction of H2 

1470 

_K.E. Higgs et al.                                                                                                                                                                                                                                 International Journal of Hydrogen Energy 71 (2024) 1468–1485_ 

_International Journal of Hydrogen Energy 71 (2024) 1468–1485_ 

**Fig. 2. A)** Taranaki Basin map showing the main structural elements. **B)** Lithostratigraphic nomenclature, stratigraphic succession, and petroleum fields in Taranaki Basin relative to basin history; modified from King and Thrasher [76] and Strogen et al. [79]. 

purity [11,23]. 

Taranaki is an oil and gas/condensate province, with over 30 petroleum fields representing potential UHS sites and representing over 750 TWh of hypothetical potential H2 storage capacity (Fig. 3). Most commercial discoveries are located onshore in a zone immediately west of the Taranaki Fault (Eastern Mobile Belt on Fig. 2A), with fewer discoveries located in offshore western parts of the basin (including the giant Maui Field with 4465 PJ/1350 PJ gas/oil reserves [84]; Fig. 4). Gas/condensate has been produced from the largest structures and oldest reservoir intervals (Paleogene), which represent the rift to passive margin and early contractional stages of basin history (Fig. 2B); oil is typically more prevalent in younger active margin reservoirs (Neogene), often within smaller structures and at shallower burial depths. Published data by MBIE [84] suggest that although many fields are still in production, the lifetime for these fields is limited and most operations are likely to cease production between 2030 and 2050. Early identification of potential oil and gas fields due for decommissioning in the next 10–20 years is essential if the projected annual hydrogen demand is to be realised and decarbonisation goals are to be met. 

# _3.3.2. Deep aquifers_ 

Aquifers are porous formations where the pore space is filled by fresh or saline water; geostorage is generally limited to saline aquifers due to the risk of contamination of potable aquifers. UHS in a deep saline aquifer requires the same conditions as a depleted oil or gas reservoir, namely good reservoir properties, and a trap and caprock (top seal) to prevent leakage. As H2 is injected into an aquifer, some of the water will be displaced (downward or sideways) due to the density difference and immiscibility between gas and liquid, which increases the pore fluid 

pressure [12]. This is technically more challenging than storage in depleted hydrocarbon reservoirs [47] and would also require high reservoir permeability; Tarkowski et al. [21] suggest a minimum permeability of 100 mD for gas storage in water-bearing rocks. UHS in a deep aquifer would also require more cushion gas compared to depleted hydrocarbon reservoirs (e.g., 80% of the total gas volume [14,21] compared to ~30–50% for depleted fields [22]), which would increase installation costs and decrease the effective H2 storage volume. However, the absence of hydrocarbons would have the benefit of being able to maintain H2 purity (subject to possible mineral and/or microbial effects) where H2 could serve as the cushion gas. 

There are many potential sites for UHS in deep aquifers in Taranaki Basin, and the total storage capacity of deep aquifers is understood to be much greater than that of the oil and gas reservoirs. However, the main issue with these sites is the more limited geological data required for assessing site suitability, together with the lack of proven trap, reservoir and caprock. Aquifers will therefore be significantly more expensive to appraise and develop than depleted oil and gas reservoirs due to the uncertain geological constraints and separation from natural gas infrastructure. As such, deep aquifers are considered a lower priority for demonstrating technical feasibility of UHS in A-NZ and should be explored once UHS in geological formations has been proven. 

# **4. Methodology** 

UHS site evaluation has been undertaken on all significant past and present depleted/depleting fields in Taranaki Basin, one undeveloped (non-commercial) gas discovery, and one deep aquifer site (Fig. 4). Geological suitability of a storage site is dependent on many factors that 

1471 

_K.E. Higgs et al.                                                                                                                                                                                                                                 International Journal of Hydrogen Energy 71 (2024) 1468–1485_ 

_International Journal of Hydrogen Energy 71 (2024) 1468–1485_ 

**Fig. 3.** Recoverable reserves for **(A)** gas/condensate and **(C)** oil fields in Taranaki Basin (logarithmic scales), with estimated equivalent H2 storage ( **B & D** respectively). Storage capacity for all reservoirs in single or combined fields is based on 2P total recoverable hydrocarbon reserves published in MBIE Energy in New Zealand (2021; 2017 for Tui, Pateke, Amokura; 2018 for Moturoa) https://www.mbie.govt.nz/building-and-energy/energy-and-natural-resources/, and FirstGas [75] for Ahuroa. Refer to supplementary files and Table 6 for information on calculation of estimated H2 storage capacity. 

have been discussed in several peer-reviewed publications [3,6,11,12, 21–23,36]. Here, we have approached prospectivity analysis using criteria-driven processes based on the current understanding of technical requirements for UHS. We introduce two methodologies for integrating the assessments of an expert panel on these criteria and their relative importance. The panel was composed of ten geoscientists with differing backgrounds and expertise including sedimentology, reservoir geology, geoenergy, petroleum geology, structural geology, geophysics, reservoir engineering, engineering geology, paleontology, volcanology, and geochemistry. Key parameters, criteria and weightings, were first assessed by the panel to provide a consistent and objective approach across all evaluated sites. Subsequent consultations included a series of dedicated discussions focused on reviewing sites and scoring values for the identified parameters. 

Eighteen technical parameters were identified for evaluation, which are broadly subdivided into the following groups: (1) site basics, (2) reservoir injectivity, (3) H2 containment, (4) H2 reactivity/contamination, and (5) data availability/quality. Some parameters require specific (quantifiable) conditions, whilst others are qualitative and are evaluated based on whether conditions are favourable for installing a UHS site. The main technical considerations for UHS and how they have been applied in our evaluation of potential UHS sites in Taranaki Basin are summarised in Table 1. Detailed information on the rationale for parameters and their criteria is provided in supplementary data file 1, with site details, associated parameter results and data references in supplementary data file 2. 

Potential H2 storage capacity is one of the key parameters used in our site evaluation and has been estimated from documented recoverable petroleum reserves by using a hydrogen energy conversion factor for natural gas [22] or an equivalent hydrogen mass for oil volume based on specific reservoir temperatures and pressures [86]. The storage capacity of very large prospects has been assessed for potential smaller compartments using full field reserves data and/or an evaluation of publicly available legacy data (e.g., individual reservoir reserve estimates, production data, reservoir models etc.). Similarly, legacy datasets have 

been reviewed to assess the potential capacity of specific sites where recoverable reserves are documented for combined fields, or where there is no reserves information (undeveloped sites). These methods and results are presented in supplementary data files 1 and 2, with detail on H2 storage capacity estimates for large prospects requiring smaller compartments presented in supplementary data file 3 (see also Table 6). 

The methodology applied is heavily dependent on data quality and availability, and we have specifically used parameters and criteria that are quickly and easily assessed from publicly available datasets (e.g., Ref. [94]), and do not require detailed technical field appraisal. Both approaches to prospectivity analysis can be adapted for different sedimentary basins based on geological characteristics of the region, and project-specific requirements. 

# _4.1. Decision tree approach to site evaluation_ 

A binary yes-or-no decision tree technique was used to evaluate the various site options in Taranaki. A series of binary decisions was first considered, based on criteria and values of geological evaluation parameters considered applicable for UHS (Table 1). A hierarchical model was subsequently constructed, whereby the most fundamental decisions were posed first. 

The tree starts at the root node, storage size, which we propose is an important factor to consider for UHS in A-NZ. Property values have been assigned to each node, and in most cases the decisions relate to a yes/no answer for criteria that we consider pertinent for UHS. Each branch of the decision tree generates a new node; square leaf nodes indicate another decision to be made, while circle leaf nodes indicate a chance event or unknown outcome. Questions posed by the square leaf nodes and their thresholds, where quantitative, are shown on Fig. 5 (also refer to supplementary data file 1). Four types of circle leaf node have been assigned as follows: red = site not considered further, orange = unlikely to be economic, yellow = additional information required, green = preferred option. “Additional information” has been recommended where the datasets and/or the knowledge is deemed inadequate for 

1472 

_International Journal of Hydrogen Energy 71 (2024) 1468–1485_ 

_K.E. Higgs et al.                                                                                                                                                                                                                                 International Journal of Hydrogen Energy 71 (2024) 1468–1485_ 

**Fig. 4.** Map showing location of 40 assessed sites for UHS in Taranaki Basin, including all significant past and present petroleum fields; site 40 is an uncommercial gas discovery, and site 23 has no proven hydrocarbons and represents a deep aquifer site. Petroleum field locations from NZP&M. 

assessment. “Unlikely to be economic” has been applied to those sites where the geological features are considered unfavourable or challenging (e.g., storage size is too small, or burial depth is too deep). However, all options would require a full economic assessment after completing site characterisation and before proceeding with UHS. 

This approach efficiently evaluates the feasibility of potential sites and their respective fulfilment of the objective according to quantitative or qualitative and deterministic criteria. Each decision (square leaf node) has been assigned a number, starting at one for the root node, to compare potential UHS site positions on the tree. 

# _4.2. Matrix approach to site evaluation_ 

A criteria-driven matrix approach has historically been used by the petroleum industry for volumetric assessments and prospect analysis (e. g., Refs. [95,96]), and more recently to identify and prioritise sites for CCS [26,97,98]. We have applied a similar weighted decision matrix approach to potential UHS sites in Taranaki to compare results with the decision tree and generate a numerical ranking of sites. An initial site 

screening was first applied, based on selected parameter threshold criteria that we consider minimum requirements for the successful deployment of UHS in A-NZ. These are: H2 capacity 1–20 TWh (working gas + cushion gas), top reservoir depth _>_ 1 km, and mean reservoir permeability _>_ 5 mD (refer Table 1). 

The minimum and maximum storage capacity criteria are based on predicted energy storage requirements for A-NZ (see section 3.1). The minimum depth criterion (1 km) is also specific to A-NZ, whereby shallow reservoirs ( _<_ 1 km) are very poorly consolidated and have typically experienced drilling problems during hydrocarbon production. It is therefore expected that repeated injection and withdrawal cycles, as required for UHS, will result in major issues related to reservoir and caprock degradation. While permeability of the A-NZ reservoirs is highly variable, several fields produce natural gas from tight reservoirs (e.g., Mangahewa, Turangi), where the average permeability is low ( _<_ 1–5 mD). A fairly low average permeability of 5 mD was selected as the minimum criterion to exclude most of these tight gas reservoirs, which generally require fracture stimulation, while acknowledging that a lower permeability is expected for producing H2 compared to natural 

1473 

_K.E. Higgs et al.                                                                                                                                                                                                                                 International Journal of Hydrogen Energy 71 (2024) 1468–1485_ 

**Table 1** 

Technical parameters and the assessment methods applied to UHS prospectivity analysis. The main technical considerations for each parameter and how they have been applied to Taranaki are summarised, with details provided in supplementary data file 1. 

|Technical Group|Technical<br>Parameter|Technical Considerations|Application in Taranaki|Assessment<br>Method|
|---|---|---|---|---|
|**Site Basics**|**1) Structural**<br>**style - trap**|• Defined as structural, stratigraphic, or combined;<br>• Trap complexity increases risk of closure&storage<br>capacity;<br>• Fault rock increases trap uncertainty;<br>• A stratigraphic component requires a well-<br>constrained geological model.|Dominated by structural traps formed by<br>compressional tectonics; good–poor:<br>• Simple 4-way dip anticline;<br>• 3-way dip anticline with fault seal;<br>• Anticline with fault compartmentalisation;<br>• Overthrust anticline;<br>•Multiple thrust slivers.|Qualitative:<br>expert elicitation|
||**2) H2 storage**<br>**capacity**|• Seasonal H2storage cycles should have the capacity<br>to last several months [4,9]<br>• This requires larger storage capacity than tanks and<br>lt  [124685]|<br>• H2storage capacity is estimated from recoverable<br>petroleum reserves using gas to H2recovery factor<br>[9,22], or oil volume to H2factor based on local<br>diti[86]|Quantitative<br>(TWh calculated)|
|||sa caverns,, <br>• However, limiting the size of the storage container<br>will have economic and operational advantages<br>• We suggest optimal storage capacity should be<br>based on site-specific energy requirements|conons;<br>• Sites or compartments with**capacities of 1**–**20 TWh**<br>**H2 are considered for UHS in Taranaki**(working<br>+ cushion gas);<br>• New Zealand UHS requirements can be met over the<br>short term by 1–3 sites providing combined storage<br>f 510 TWh H||
||**3) Reservoir**<br>**porosity**|• Reservoir porosity is a measure of free space in a<br>rock that can potentially be used to store gas[21];<br>• Total porosity is routinely measured on samples<br>from hydrocarbon wells;<br>• Total porosity∕=effective porosity and is not a direct<br>f  i|o– 2.<br>• Core porosity data is provided in well reports;<br>• Effective porosity calculations are not available for<br>all sites;<br>• It is assumed that reservoirs with higher (average)<br>total porosity are favourable for UHS.|Quantitative<br>(average<br>measured %)|
||**4) Storage depth**|measure o storage capacty.<br>• Deeper sites (with higher T&P) provide more H2<br>storage capacity than shallower sites[22,87], but<br>this is offset by the cost of drilling deep[26].<br>Previous authors have suggested:<br>• UHS should be in the range 500–2000 m[2,36];<br>• Offshore UHS should be_>_1.5 km to ensure H2<br>densities of 10 kg m<sup>−3 </sup>[87].|• Shallow reservoirs in Taranaki comprise<br>poorly–unconsolidated sands and are unsuitable for<br>H2storage;**minimum depth to top reservoir of 1**<br>**km has been applied;**<br>• **Optimal depth is considered to be 1.5**–**3.5 km**for<br>UHS in Taranaki but with many good reservoirs in<br>the better consolidated and more deeply buried<br>Paleoene strata|Quantitative (km<br>to top reservoir)|
|**Reservoir**<br>**Injectivity**|**5) Reservoir**<br>**permeability**|• Reservoir permeability is the ability of a porous<br>medium to transmit fluid when saturated by the<br>fluid [88];<br>• Permeability is routinely measured on samples from<br>hydrocarbon wells;<br>• Higher permeability= higher injectivity [87];<br>• Low viscosity of H2will result in enhanced mobility<br>compared to natural gas or CO2[6].|g .<br>• Permeability is provided in well reports;<br>• It is assumed that reservoirs with higher<br>permeability are favourable for UHS;<br>• **A minimum (average) permeability of 5 mD has**<br>**been applied;**<br>• No differentiation has been made at this stage for<br>depleted hydrocarbon fields and deep aquifers.|Quantitative<br>(average<br>measured mD)|
||**6) Gas**<br>**production rates**|<br>• Gas production rate represents the rate per unit time<br>that natural gas is produced;<br>• It is generally proportional to reservoir<br>permeability;<br>• H2can be withdrawn and injected at higher volume<br>rates than natural gas[16].|• Maximum historical yearly gas production rates and<br>maximum 2021 rates are provided by MBIE[84];<br>• It is assumed that higher gas production rates are<br>favourable for UHS.|Quantitative<br>(max yearly TJ/<br>day)|
||**7) Stratigraphic**<br>**heterogeneity**|<br>• Stratigraphic heterogeneity refers to lithological<br>complexity (e.g., bed thickness);<br>• High heterogeneity can result in injectivity and<br>operational issues (e.g., pressure compartments);<br>• Need to consider depositional facies/lithology&<br>faults with clay gouge.|• A preliminary qualitative assessment has been based<br>on available well reports;<br>• Quality and quantity of pressure, facies and<br>geometry data is highly variable;<br>• Sites with low heterogeneity are favoured over<br>complex sites with numerous pressure<br>comartments|Qualitative:<br>expert elicitation|
||**8) Multiple**<br>**production pools**|• Multiple production pools refers to several small<br>accumulations produced under the same field name;<br>• Multiple (small) production pools will result in<br>lower injectivity due to reservoir<br>compartmentalisation.|p.<br>In Taranaki multiple (small) pools include:<br>• Small pools in the same structural closure but within<br>different formations;<br>• Pools at different levels of one formation;<br>• The same reservoir interval but within structurally<br>dit t|Qualitative:<br>expert elicitation|
||**9) Potential**<br>**formation**<br>**damage**|• Formation damage during drilling can cause<br>production problems and impair injectivity by<br>reducing permeability in the near-wellbore region<br>[89];<br>• Need to consider how the formation will perform<br>over the multiple injection and production cycles<br>id f UHS|scree raps.<br>• A preliminary qualitative assessment of potential<br>formation damage has been based on well reports<br>detailing reservoir texture and composition;<br>• In Taranaki the main considerations are the<br>potential for sand production, fines migration and<br>presence of swelling clays.|Qualitative:<br>expert elicitation|
||**10) Aquifer risk**|requre or .<br>• A strong aquifer in a depleted petroleum field will be<br>detrimental to H2injectivity[21];<br>• Water production in UHS simulations occurs more<br>quickly in the water zone compared to oil or gas<br>zones[51];<br>• Fields using depletion drive recovery are favourable<br>to aquifer support recovery [90].|• A preliminary qualitative assessment of aquifer risk<br>has been based on well reports;<br>• However, public data on aquifer support and water<br>influx are not easily accessible;<br>• Aquifer risk has been qualitatively ranked, with a<br>high uncertainty associated with this parameter.|Qualitative:<br>expert elicitation|

( _continued on next page_ ) 

1474 

_K.E. Higgs et al.                                                                                                                                                                                                                                 International Journal of Hydrogen Energy 71 (2024) 1468–1485_ 

_International Journal of Hydrogen Energy 71 (2024) 1468–1485_ 

## **Table 1** ( _continued_ ) 

|Technical Group|Technical<br>Parameter|Technical Considerations|Application in Taranaki|Assessment<br>Method|
|---|---|---|---|---|
|**H2 Containment**|**11) Hydrocarbon**<br>**phase**<br>**12) Top seal risk**|• The composition of the depleted, trapped<br>hydrocarbon will have specific physico-chemical<br>properties different to H2[19–24];<br>• Gas systems generate higher buoyant pressures than<br>oil systems suggesting that a depleted gas reservoir<br>may be more likely to hold back H2than a depleted<br>oil reservoir;<br>• H2is more diffusive, with lower viscosity and<br>density, and is hence potentially more mobile than<br>methane or CO2[6,21];<br>• H2has a lower solubility in water than other gases,<br>reducing risk of dissolution into the seal[21].<br>• Top seal is dependent on lithology/facies<br>(composition, continuity, quality) and wells;<br>• Seal quality is based on capillary pressure&pore<br>size distribution measurements[93];<br>• A given caprock is likely to retain at least as large a<br>vertical column of H2as natural gas[22,87]and H2<br>losses through caprocks are likely to be low[22,65];<br>• Further work is required to assess the impact of H2<br>on seal properties[2223]|• Depleted gas reservoirs are favoured over depleted<br>oil reservoirs on the premise that caprocks proven to<br>contain natural gas are more likely to contain H2<br>than caprocks for liquid hydrocarbons;<br>• A gas cap in a mixed phase reservoir is considered<br>more favourable than an oil only reservoir;<br>• Many of deeper (older) reservoir intervals in<br>Taranaki contain gas; shallow (younger) reservoirs<br>are oil prone (updip oil spill [91,92]);<br>• Deep aquifer sites have no proven containment.<br>• Effective hydrocarbon seals are proven in the<br>depleted fields in Taranaki;<br>• A preliminary qualitative assessment of top seal<br>capacity&continuity has been based on available<br>well reports and an understanding of facies trends;<br>• The number of well penetrations into the reservoir<br>interval has been included in the assessment for<br>potential top seal risk.|Oil, gas/cond, 2-<br>phase, deep<br>aquifer<br>Qualitative:<br>expert elicitation|
||**13) Fault seal**<br>**risk**|,.<br>• Fault seal risk is difficult to determine and requires<br>evaluation of the fault shale gouge ratio and<br>buoyancy pressures [93];<br>• An individual fault can locally provide a seal (e.g.,<br>pressure compartments) or a leakage pathway;<br>• Fault and seal integrity can be compromised by<br>drilling activity.|• Detailed fault analysis in Taranaki is generally<br>sparce and inconsistent, and a preliminary<br>qualitative assessment has been based on available<br>well reports;<br>• It is assumed that sites with extensive faulting<br>through the caprock and/or with a fault-seal<br>component to the trap are less favourable for UHS.<br>i|Qualitative:<br>expert elicitation|
|**H2 Reactivity/**<br>**Contamination**<sup>**a**</sup>|**14) H2 **–**CO2**<br>**15) H2 **–**coal**<br>**16) H2 **–**oil**<br>**17) H2 **–<br>**carbonate**|• Work to date suggests that there is likely to be little<br>reaction between siliciclastic minerals and H2over<br>UHS time-scales[32,57–61];<br>• However, the potential for biochemical reactions is<br>considered a major uncertainty[63,65]with H2<br>providing the main source of energy for<br>carbon-based minerals and hydrocarbons[21,22,32,<br>60,65,68], sulphates and sulphides[12,23,57,61,<br>65];<br>• These reactions can lead to H2losses and<br>contamination[3,12,23,55,63,65], and/or changes<br>in rock properties[56];<br>• Further work is required to assess potential<br>reactivity under reservoir conditions[57,59,63].|• A qualitative, site-specific assessment has been made<br>for reservoir/seal composition.<br>In Taranaki:<br>• Reservoir and seal rocks are dominated by<br>siliciclastic minerals;<br>• CO2content is highly variable;<br>• Coals/carbonaceous matter is relatively common in<br>Paleocene-Eocene reservoirs;<br>• Carbonate cements are relatively common in<br>Oligocene and Miocene strata;<br>• Sulphates are generally not present;<br>• Pyrite is generally minor (≤1%);<br>• It is assumed that reservoir and seal rocks with low<br>CO2, coal, carbonate and oil are favourable for UHS.|Qualitative:<br>expert elicitation|
|**Data**|**18) Data**<br>**availability/**<br>**quality**|<br>• Data availability, vintage and quality of datasets are<br>of paramount importance for UHS prospectivity<br>analysis.|<br>• In New Zealand there is a requirement to lodge data<br>with MBIE, which becomes publicly available after 5<br>years;<br>• Potential sites have been qualitatively assessed for<br>data availability and quality.|Qualitative:<br>expert elicitation|

> a reactivity of H2-sulphide/sulphate would need to be addressed in the site characterisation work. 

gas (due to the lower viscosity of H2) [6]. An effective trap and top seal are two of the most important geological factors for UHS, but these parameters were not included at the preliminary screening stage given that UHS prospectivity has focused on developed fields with proven hydrocarbon accumulations (hence inferred trap and caprock). However, a minimum top seal thickness and/or maximum seal permeability would be recommended for screening of deep aquifers where trap and containment are largely unknown. 

Sites that met the matrix-qualifying criteria have been assessed and ranked based on the geological parameters outlined in Table 1. For these sites, a value from 1 (poorest) to 3 (best) was assigned to each parameter. All parameters were also given a weighting from 1 (lowest) to 5 (highest), reflecting our understanding of the relative importance of different attributes for UHS. Confidence scores (high, mid, low) are provided for both parameter values (by site) and weightings. The rank score for each potential storage site was calculated as a value between 0 and 1 by summing the normalised, weighted parameter values (normalised criterion value x normalised weighting). As such, the score reflects a combined assessment of the overall technical suitability for UHS at each site. The data confidence and weighting confidence scores are not factored into this technical score but provide a measure of the 

uncertainties relating to individual sites and parameters. 

This weighted decision matrix approach is subject to the biases of the authors and is very much dependent upon our current understanding of the most favourable parameters for UHS. The method used was preferred over a paired comparison matrix that assesses relative importance of two factors [4,26] because of the considerable uncertainties that remain regarding H2 behaviour in the subsurface [6,22, 23]. For example, there are high uncertainties with different H2 reactivity parameters and how they will individually impact H2 purity (Table 1), and we consider a pairwise comparison of these parameters would be more susceptible to author bias than the broad low–high (1–5) weighting scheme we have applied. We have undertaken some sensitivity analysis to address the uncertainty with parameter weightings, which assesses how rank score changes with different weightings. Sensitivity analysis was concentrated on parameter weightings that have the greatest uncertainties, and in this way, we were able to investigate how bias might influence the results. Weighting factors are liable to be changed and adapted as our understanding of the importance of specific parameters and criteria for UHS evolves. 

1475 

_K.E. Higgs et al.                                                                                                                                                                                                                                 International Journal of Hydrogen Energy 71 (2024) 1468–1485_ 

_International Journal of Hydrogen Energy 71 (2024) 1468–1485_ 

**Fig. 5.** Summary decision tree for evaluating the technical suitability of UHS sites. Square nodes represent questions with numbers reflecting the tree position. Solid line branches represent a “yes” response and dashed branches represent “no”. The position of sites is shown on the tree by the site name; reservoir interval in parentheses where there are multiple reservoir intervals. 

# **5. Site prospectivity** 

A total of 40 sites have been assessed in this work, which are referred to by their site name and reservoir interval (site, reservoir; Fig. 4). They comprise gas fields, oil-bearing fields with gas caps, oil fields, one deep aquifer stratigraphically located above a commercial field (Kapuni, Moki), and one site with proven, but undeveloped (non-commercial) gas (Karewa, Mangaa). Sites without commercial production have some subsurface geological data, but are poorly characterised compared to the depleted oil and gas reservoirs. Site details, parameter results and data references are provided in supplementary data file 2. 

# _5.1. Decision tree results_ 

The summary decision tree is presented as Fig. 5, which illustrates the position of all 40 potential sites. Sites that are positioned high on the tree structure (e.g., nodes 2–3) are considered to have poor geological characteristics for UHS. They include sites ending on a red node (fail), or an orange node (unlikely to be economic), based on whether they meet the minimum and maximum size criteria (node 1), depth criteria (nodes 2), and minimum permeability criterion (nodes 3). 

To progress down the main part of the tree structure the UHS site must have a proven oil or gas accumulation (node 4), which demonstrates a working trap-reservoir-caprock. If there is no proven accumulation the site will stop at node 4 and require additional information (yellow node). In the current study we have only included one deep aquifer where there are no proven hydrocarbons (Kapuni, Moki). However, any other deep aquifer that meets the storage size, permeability and depth criteria would finish at this node. 

Sites that present further down the tree structure are considered more suitable for demonstrating UHS, with two branches relating to 1) natural gas reservoirs, including gas caps, and 2) oil only reservoirs. Natural gas reservoirs occupy the most favourable tree branch for UHS. Parameters that are considered along this branch comprise structural style, data availability/quality, number of production pools, potential 

for fault or top seal leakage, and potential reactivity. This represents a simplified set of parameters compared to the matrix approach (section 5.2). All sites that proceed along the oil only reservoirs branch are considered to have potential to leak gas through the seal, finishing at node 6, and requiring additional information (yellow node; e.g., Cheal). However, a greater understanding on H2 mobility is required to assess the competence of seals for UHS, and the position of oil fields may change depending on results from future research. 

Table 2 summarises the decision tree results, which lists nine eliminated sites (red nodes), two candidates with favourable geological characteristics (green nodes), and the relative position of other sites. The top two candidates are the natural gas storage site Ahuroa (Ahuroa, Tariki), and the offshore Maui A field (Maui A, Kaimiro). However, it is noted that the Ahuroa tree position is dependent upon provision of confidential data at site characterisation stage. Other potential UHS sites (finishing on nodes 7–9) are all Paleogene depleted gas or gas-cap reservoirs, which have ended on a yellow node for potential reactivity, H2 leakage (caprock or fault), or data availability. 

The position of sites on the decision tree is partly controlled by the parameter order (i.e. order that questions are posed). For example, if a question posed relatively high in the tree structure is switched with a question posed low in the tree structure (e.g., permeability switched with potential reactivity), then the resultant node position of sites would be different, although the node colour would be the same (e.g., yellow for additional information required, or green for preferred option). Site position on the decision tree is therefore dependent upon our current understanding of the most favourable parameters for UHS (e.g., sufficient permeability is considered more favourable for demonstrating UHS than the absence of potential reactants). 

# _5.2. Matrix results_ 

A list of criteria and relative weightings for the 18 technical parameters used in site ranking is presented in Table 3. The quality and quantity of data used to assess potential storage sites has a significant 

1476 

_K.E. Higgs et al.                                                                                                                                                                                                                                 International Journal of Hydrogen Energy 71 (2024) 1468–1485_ 

_International Journal of Hydrogen Energy 71 (2024) 1468–1485_ 

**Table 2** 

Technical suitability of sites for UHS in Taranaki Basin based on a decision tree approach. 

|#|Site Name|Reservoir Interval|Deci|sion Tree Node Position and Results||
|---|---|---|---|---|---|
|**1**|Ahuroa<sup>a</sup>|Tariki|**10**|**Green Node¼ Geologically Preferred Ca**|**ndidates**|
|**2**|Maui A<br>i|Kaimiro|**9**|||
|**3–5 **|Kapuni N field|Mangahewa (K1A)|**9**|**Yellow Node 9**, potential reactivity|**Geologically Possible Candidates: Additional information**|
||Maui B|Kaimiro|**9**||**required**|
||McKee<sup>a</sup>|McKee|**9**|**Yellow Node 9**, potential leakage||
|**6–9 **|Kupe|Farewell|**8**|**Yellow Node 8**, potential leakage (fault||
||Pohokura|Mangahewa (U1)|**8**|or seal)||
||Maari|Mangahewa|**8**|||
||Karewa|Mangaa|**8**|||
|**10**|Tariki|Tariki|**7**|**Yellow Node 7**, data needed to address||
|||||uncertainties||
|**11–16**|Maui B|Farewell|**6**|**Yellow Node 6**, oil reservoir with no gas|**Lower UHS potential based on current understanding: Additional**|
||Manaia|Mangahewa|**6**|containment|**information required**|
||Maari|Moki|**6**|||
||Kaimiro|Mount Messenger|**6**|||
||Ngatoro|Mount Messenger|**6**|||
||Cheal and Cheal|Urenui&Mount Mess|**6**|||
||E|||||
|**17**|Kapuni|Moki|**4**|**Yellow Node 4**, no proven containment||
|**18–31**|Tui area<sup>b</sup>|Farewell|**3**|**Orange Node 3**, deep reservoir|**Unlikely to be economic**|
||Maui A|Mangahewa|**2**|**Orange 2**, container size too large||
||Maui B|Mangahewa|**2**|||
||Rimu|Tariki, Rimu Lst,|**2**|**Orange Node 2**, container size too small||
|||Basement||||
||Copper Moki|Mt Mess and Moki|**2**|||
||Supplejack|Mt Mess and Moki|**2**|||
||Surrey|Mount Messenger|**2**|||
||Sidewinder|Mt Mess and Moki|**2**|||
||Windsor|Mount Messenger|**2**|||
||Goldie|Mount Messenger|**2**|||
||Moturoa|Matemateaonga|**2**|||
||Kahili|Tariki|**2**|||
|**32–40**|Kauri|Kauri|**3**|**Red Node 3**, failed on permeability|**Site not considered further**|
||Waihapa/|Tikorangi|**3**|||
||Ngaere|||||
||Mangahewa|Mangahewa|**3**|**Red Node 3**, failed on depth and||
||Ohanga/Onaero|Mangahewa/McKee|**3**|permeability||
||<br>Kowhai|<br>Mangahewa|**3**|||
||Turangi|Mangahewa|**3**|||
||Manutahi|Matemateaonga|**2**|**Red Node 2**, failed on container size and||
||Cardiff|Mangahewa|**2**|depth||
||Radnor|Mangahewa|**2**|||

> a site position dependent upon confidential data. 

> b Tui area comprises three separate fields (Tui, Pateke, Amokura). 

## **Table 3** 

Criteria and relative weightings applied to parameters assessed for UHS prospectivity analysis, Taranaki Basin. Weighting values from 1 (low) to 5 (high). Normalised weighting = weighting value/sum of weighting values. 

|Technical Group|Technical Parameter|Criteria and Va|lue||Weighting|Weighting<br>i|Normalised|
|---|---|---|---|---|---|---|---|
|||Poor - 1|Mid - 2|Good - 3||Confidence|Weighting|
|**Site Basics**|1) Structural style|Complex|Mid|Simple|**3**|High|0.056|
||2) Storage capacity (TWh)|15–20|1–2, 8–15|3–8|**4**|Mid|0.074|
||3) Reservoir porosity (av. %)|_<_10|10–15|_>_15|**2**|High|0.037|
||4) Storage depth (km)|1–1.5,_>_3.5|2.5–3.5|1.5–2.5|**3**|Mid|0.056|
|**Reservoir**|5) Reservoir permeability (av.|5–10|10–100|_>_100|**4**|High|0.074|
|**Injectivity**|mD)|||||||
||6) Stratigraphic heterogeneity|High|Mid|Low|**1**|Mid|0.019|
||7) Multiple production pools|Many pools|2–3 contacts|Single<br>accumulation|**3**|Mid|0.056|
||8) Gas production rates (TJ/d)|_<_5<br>i|5–50|_>_50|**2**|Mid|0.037|
||9) Potential formation damage<br>l|Confirmed<br>i|Suspected|No indication|**2**|Mid|0.037|
||10) Water influx|i<br>Confirmed|Suspected|No indication|**1**|Low|0.019|
|**H2 Containment**|l<br>11) Hydrocarbon phase|i<br>Oil reservoir<br>i|Gas cap|Gas reservoir|**5**|Mid|0.093|
||12) Fault seal risk|Confirmed<br>i|Suspected|No indication|**5**|High|0.093|
||13) Top seal risk|i<br>Confirmed|Suspected|No indication|**5**|High|0.093|
|**H2 Reactivity**|14) H2–CO2reactivity|i<br>_>_20% CO2|3–20% CO2|_<_3% CO2|**3**|Low|0.056|
||15) H2-coal reactivity|Common coal|Rare coal|No coal|**2**|Low|0.037|
||16) H2-oil interaction|Common oil|Rare oil|No oil|**3**|Low|0.056|
||17) H2-carbonate reactivity|Common<br>carb.|Locally<br>carb.|Rare carb.|**1**|Low|0.019|
|**Data**|18) Data Availability/Quality|Poor|Fair|Good|**5**|High|0.093|

1477 

_International Journal of Hydrogen Energy 71 (2024) 1468–1485_ 

_K.E. Higgs et al._ 

influence on the results, and hence has been given the highest weighting (parameter 18, weighting 5). Highest weighting values have also been assigned to parameters that relate to containment, including risk of leakage through top seal or fault seal (parameters 12 and 13), and hydrocarbon phase (whereby a natural gas accumulation is considered more likely to contain H2 compared to an oil reservoir with no proven gas cap; parameter 11). A high weighting (weighting 4) has been applied to reservoir permeability (parameter 5), which is essential for injectivity and to storage capacity. Lower weightings have been placed on parameters where there remains a high uncertainty with respect to their impact on hydrogen storage, and/or where there is limited information in the database (making comparisons difficult and resulting in high uncertainty). It is anticipated that future research should help to resolve some of these uncertainties, which may result in a revision to the parameter weightings. Additionally, in cases where the data are sparse or of poor quality, collection of new data may ultimately change the site scores. 

Site ranking results are presented in Table 4. The top-ranking sites all comprise Paleogene depleted gas or gas-cap reservoirs, with the lowest ranking fields representing younger (Miocene) reservoirs and oil reservoirs from older (Paleogene) accumulations. Typically, the Miocene reservoirs rank low due to their small size, reservoir complexity/compartmentalised petroleum pools (which increases risk associated with reservoir injectivity), and their uncertainty for trapping gas (generally oil is the primary reservoir fluid). Additionally, these Miocene reservoirs are often poorly lithified, which, with the fine grain size, labile and clayrich composition, makes them prone to reservoir degradation that could worsen over time due to the multiple injection/withdrawal cycles required for UHS. 

The lowest matrix score is from the deep aquifer site (Kapuni, Moki) that overlies the commercial Kapuni Field at Mangahewa level. This very low score is primarily due to a significant risk of containment, and uncertainties associated with lack of publicly available data/interpretations. The second site with no commercial discovery (Karewa, Mangaa) is ranked at number 11, above the Miocene plays and some Paleogene oil sites. Karewa has some good quality data, including a 3D seismic survey and single borehole, with data availability and quality better than expected for a deep aquifer site. Data and interpretations are, however, much less detailed than those for a developed field. Critical factors that hinder UHS development of this site relate to the uncertainty with containment (fault and top seal risk), and reservoir quality (heterogeneity, production rates, formation damage). A list of disqualified sites from the screening stage is provided in supplementary data file 2, with size being the main basis for ruling out 

potential sites. Many of the sites that are deemed to be too small for anticipated UHS requirements are onshore Miocene oil fields (e.g., Surrey). Several of the stratigraphically older reservoirs produce from deep/low permeability and locally small reservoirs, which would make UHS challenging (e.g., Cardiff) and have also been excluded. Only two sites have been eliminated due to storage capacities that are too large, and where a smaller compartment was not considered possible (Maui A and Maui B, Mangahewa; refer to supplementary data file 3). 

# _5.3. Sensitivity analysis_ 

There is a very small range in matrix scores for the top 10 sites (0.806 to 0.725, Table 4) and site order could change with input from different experts. We have therefore undertaken sensitivity analysis, primarily to investigate the impact of changing parameter weightings where there remains significant uncertainty with the weighting and/or the impact of the parameter on H2 storage (e.g., H2 reactivity). We have consistently applied the highest weightings to containment parameters, permeability, storage capacity, and data availability/quality. 

A ranking comparison for some of these sensitivity analyses is provided in Table 5. In all cases the top 10 sites are dominated by Paleogene depleted gas or gas-cap reservoirs, with lower ranks associated with the younger Miocene reservoirs. However, there is some variability in the actual rank position for individual sites, which is dependent on the importance placed on different variables for the successful deployment of UHS. 

The top three candidates are variably the natural gas storage site Ahuroa (Ahuroa, Tariki), the offshore Maui A field (Maui A, Kaimiro), the onshore Kapuni field (Kapuni, K1A Mangahewa), and the on/ offshore Pohokura field (Pohokura, U1 Mangahewa). In most cases, Ahuroa is the preferred candidate, but its position is partly dependent upon high-quality, confidential (UGS) data. An adjacent field (Tariki, Tariki) is broadly analogous to Ahuroa in terms of its structural style, reservoir and seal but has a lower rank score. This result is primarily due to poor data availability/quality at Tariki, and the site rank is liable to change with acquisition of new datasets. 

Sensitivity analysis shows that results are impacted by the weightings chosen by the panel and will need to be updated to integrate learnings from new studies as they become available. Based on current knowledge, the ranking matrix suggests that there are multiple sites in Taranaki that are likely to have broadly suitable geological characteristics for UHS. 

**Table 4** 

Ranking of sites in Taranaki Basin based on a matrix approach where the site number from Fig. 3 is shown in parentheses. Rank score = sum of weighted parameter values (normalised criterion value x normalised weighting). The ranking matrix addresses geological parameters required for UHS; economic, social/cultural and regulatory issues would need to be separately evaluated for each potential site before proceeding with any pilot project. 

|#|Site Name (Site Number)|Onshore/Offshore|Hydrocarbon Phase|Reservoir Age|Reservoir Interval|Seal Interval|Rank Score|
|---|---|---|---|---|---|---|---|
|**1**|Ahuroa (17)|Onshore|Gas/cond|Oligocene|Tariki|Otaraoa|**0.806**|
|**2**|Maui A (33)|Offshore|Gas/cond|Eocene|Kaimiro|Turi|**0.796**|
|**3**|Kapuni (24)|Onshore|Gas/cond|Eocene|Mangahewa (K1A)|Turi/Otaraoa|**0.784**|
|**4**|Pohokura (1)|Offshore|Gas/cond|Eocene|Mangahewa (U1)|Turi|**0.765**|
|**5**|McKee (4)|Onshore|Gas cap to oil leg|Eocene|McKee|Otaraoa|**0.762**|
|**6**|Kupe (28)|Offshore|Gas leg with oil rim|Paleocene|Farewell|Otaraoa|**0.756**|
|**7**|Tariki (16)|Onshore|Gas/cond|Oligocene|Tariki|Otaraoa|**0.756**|
|**8**|Maui B (35)|Offshore|Gas cap to oil leg|Eocene|Kaimiro|Turi|**0.753**|
|**9**|Tui area<sup>3</sup>(37)|Offshore|Oil|Paleocene|Farewell|Turi|**0.753**|
|**10**|Maari (30)|Offshore|Gas cap to oil leg|Eocene|Mangahewa|Turi/Otaraoa|**0.725**|
|**11**|Karewa (40)|Offshore|Uncommercial gas|Pliocene|Mangaa|Manganui|**0.713**|
|**12**|Maui B (36)|Offshore|Oil|Paleocene|Farewell|Turi|**0.701**|
|**13**|Manaia (31)|Offshore|Oil|Eocene|Mangahewa|Turi/Otaraoa|**0.617**|
|**14**|Maari (29)|Offshore|Oil, minor gas|Miocene|Moki|Manganui|**0.611**|
|**15**|Kaimiro (8)|Onshore|Oil, minor gas|Miocene|Mount Messenger|Manganui|**0.599**|
|**16**|Ngatoro (12)|Onshore|Oil, minor gas|Miocene|Mount Messenger|Manganui|**0.599**|
|**17**|Cheal/Cheal E (20)|Onshore|Oil and gas|Miocene|Urenui&Mount Mess|Manganui|**0.574**|
|**18**|Kapuni (23)|Onshore|Deep aquifer|Miocene|Moki|Manganui|**0.525**|

1478 

_K.E. Higgs et al.                                                                                                                                                                                                                                 International Journal of Hydrogen Energy 71 (2024) 1468–1485_ 

_International Journal of Hydrogen Energy 71 (2024) 1468–1485_ 

**Table 5** 

Comparison of site rank based on different weightings applied to selected parameters, where the rank score is shown in parentheses. **A)** Original weightings as presented in Table 3. **B)** Reactivity parameter weightings increased for CO2 (4), coal (3), oil (4) and carbonate (3). **C)** Reactivity parameter weightings lowered for CO2 (2), coal (1), oil (2). **D)** Reservoir parameter weightings increased for stratigraphic heterogeneity (3), formation damage (4), structural complexity (4) and aquifer risk (3). Top 10 sites are Paleogene reservoirs for all scenarios and lowest ranked sites are Miocene reservoirs (italics). 

|#|A) Original Weightings|B) Reactivity–higher weighting|C) Reactivity–lower weighting|D) Reservoir Parameters|
|---|---|---|---|---|
|**1**|Ahuroa, Tariki (0.806)|Ahuroa, Tariki (0.811)|Maui A, Kaimiro (0.804)|Ahuroa, Tariki (0.831)|
|**2**|Maui A, Kaimiro (0.796)|Maui A, Kaimiro (0.797)|Kapuni, K1A Mangahewa (0.797)|Kapuni, K1A Mangahewa (0.820)|
|**3**|Kapuni, K1A Mangahewa (0.784)|Pohokura, U1 Mangahewa (0.774)|Ahuroa, Tariki (0.794)|Maui A, Kaimiro (809)|
|**4**|Pohokura, U1 Mangahewa (0.765)|Tariki, Tariki (0.766)|Pohokura, U1 Mangahewa (0.765)|Pohokura, U1 Mangahewa (0.803)|
|**5**|McKee, McKee (0.762)|McKee, McKee (0.763)|Maui B, Kaimiro (0.765)|Kupe, Farewell (792)|
|**6**|Kupe, Farewell (0.756)|Tui area, Farewell (0.763)|Kupe, Farewell (0.761)|Tariki, Tariki (0.787)|
|**7**|Tariki, Tariki (0.756)|Kapuni, K1A Mangahewa (0.757)|McKee, McKee (0.758)|McKee, McKee (0.781)|
|**8**|Maui B, Kaimiro (0.753)|Maui B, Kaimiro (0.751)|Tui area, Farewell (0.752)|Maui B, Kaimiro (0.770)|
|**9**|Tui area, Farewell (0.753)|Kupe, Farewell (0.737)|Tariki, Tariki (0.742)|Tui area, Farewell (0.770)|
|**10**|Maari, Mangahewa (0.725)|Maari, Mangahewa (0.726)|Maari, Mangahewa (0.722)|Maari, Mangahewa (0.757)|
|**11**|Karewa, Mangaa (0.713)|Karewa, Mangaa (0.720)|Karewa, Mangaa (0.703)|Maui B, Farewell (0.732)|
|**12**|Maui B, Farewell (0.701)|Maui B, Farewell (0.709)|Maui B, Farewell (0.703)|Karewa, Mangaa (0.730)|
|**13**|Manaia, Mangahewa (0.617)|_Maari, Moki (0.621)_|Manaia, Mangahewa (0.618)|Manaia, Mangahewa (0.664)|
|**14**|_Maari, Moki (0.611)_|Manaia, Mangahewa (0.619)|_Maari, Moki (0.601)_|_Maari, Moki (0.628)_|
|**15**|_Kaimiro, Mt Mess (0.599)_|_Kaimiro, Mt Mess (0.605)_|_Kaimiro, Mt Mess (0.595)_|_Kaimiro, Mt Mess (0.596)_|
|**16**|_Ngatoro, Mt Mess (0.599)_|_Ngatoro, Mt Mess (0.605)_|_Ngatoro, Mt Mess (0.595)_|_Ngatoro, Mt Mess (0.596)_|
|**17**|_Cheal&Cheal E, Uren_+ _MM (0.574)_|_Cheal& Cheal E, Uren_+_MM (0.588)_|_Cheal&Cheal E, Uren_+ _MM (0.562)_|_Cheal&Cheal E, Uren_+ _MM (0.574)_|
|**18**|_Kapuni, Moki (0.525)_|_Kapuni, Moki (0.537)_|_Kapuni, Moki (0.503)_|_Kapuni, Moki (0.546)_|

# **6. Discussion** 

# _6.1. Comparison of approaches_ 

Results are compared for the two different multi-criteria decisionmaking approaches used in UHS prospectivity analysis. The decision tree is a “fast-track” method that provides an evaluation without the need for assessing other sites. Conversely, the matrix approach is a more in-depth assessment and requires parameters for all sites prior to evaluation. 

The top two sites for validating the technical feasibility of UHS in A–NZ are the same for both approaches, demonstrating the consistent criteria applied by the expert panel. Most other sites that fall within the top ten options using the decision tree (failing at node 7 or above, Table 2) are also ranked in the top ten on the matrix. The one exception is Karewa, an undeveloped gas discovery, which occurs in the top ten of the decision tree based on proven containment (gas); this site is ranked slightly lower in the matrix (position eleven) because of the uncertainty with the size of the accumulation and the amount of hydrocarbons that might have been leaked. The relatively simple approach of the tree structure does not allow for these more complex questions to be addressed in the analysis. 

Generally, the criteria used in both approaches have been consistent and it is therefore not surprising that results for both methods are similar. Small variations in the results partly reflect differences in the application of reservoir depth criteria, and are also partly due to the combination of certain parameters in the decision tree approach. 

The decision tree has included a minimum and maximum depth query, which occurs high up in the tree structure and provides the user with a visual prompt that deeply buried reservoirs are unlikely to be economic. We have used a maximum depth of 3.5 km, but this will vary for different geographical locations, and will be dependent on geological, geographical, social, regulatory and economic factors. A maximum depth has not been applied in the early screening stage of the matrix to prevent exclusion of good quality, deep reservoir options. This results in a poorer position on the decision tree for deep reservoirs (e.g., Tui area) compared to the matrix. 

Several parameters have been combined on the decision tree, which simplifies the prospectivity assessment and does not accommodate the relative importance of different parameters. For example, a site containing many potential reactants (e.g., CO2, coal, residual oil, carbonate) is not distinguished on the tree from a site that only has potential reactivity with coal, and the relative impact of H2–CO2, H2–coal, H2–oil, 

and H2–carbonate reactions is not addressed. 

Notwithstanding small variations in results for the two approaches, both suggest that based on our current understanding, the most geologically prospective UHS sites in Taranaki Basin are Paleogene reservoirs (rift to passive margin and earliest contractional deposits, Fig. 2B), and that younger reservoirs (active margin deposits) are typically poorer candidates. The one exception to this is the Karewa site (undeveloped, sub-commercial Pliocene reservoir), which is characterised by a simpler structural trap (4-way dip closure) compared to Miocene geostorage prospects, and has a small gas accumulation, compared to oil in most Miocene reservoirs. 

Overall, the two sites without commercial hydrocarbons display similar positions for both decision tree and matrix approaches (Karewa, Mangaa positions 6–9 and 11 respectively, and Kapuni, Moki positions 17 and 18 respectively). Deep aquifer sites are downgraded in the matrix due to the absence of proven containment (hydrocarbon phase parameter = 0 where none proven, and top seal risk), with low parameter values typically also put on data availability/quality. If hydrocarbons are unproven, these sites finish early on the decision tree, with advice for more information required (yellow node). While this effectively deals with the Kapuni, Moki deep aquifer (stopping at node 4), the Karewa, Mangaa option has branched further down the tree to finish in a similar position to commercial gas reservoirs (i.e., Pohokura and Kupe, node 8; Fig. 5). These three sites are identified as all having the potential for some breach, but the presence of large gas columns at Pohokura and Kupe are not distinguished from the sub-commercial gas at Karewa. 

The strength of the decision tree approach to prospectivity is simplicity, and the ability to apply it without the need for other site information. The decision tree visually demonstrates cause-and-effect relationships, providing a simplified view of complex issues and helps to clarify the choices, risks, objectives and gains. Additional information can be documented on the tree, highlighting any technical issues related to the failure node. 

In comparison, the matrix approach provides a numerical rank score that would be beneficial for the large-scale screening of numerous storage sites. The matrix results do not present the user with an understanding of potential UHS issues. However, a visual comparison of site scores can be made using spider diagrams, the shape of which illustrates the strengths and weaknesses for different sites. Spider plots comparing the top five and bottom five ranking sites in Taranaki show that young Miocene reservoirs are characterised by spiky plots, indicative of their highly variable parameter values (Fig. 6). Most Miocene reservoirs score high values for porosity, coal/CO2 and data availability, but low values 

1479 

_International Journal of Hydrogen Energy 71 (2024) 1468–1485_ 

**Table 6** 

Estimated H2 storage capacity (working gas + cushion gas) for qualifying sites in Taranaki Basin based on recoverable 2P 2021 petroleum reserves from MBIE [84] unless otherwise stated. Natural gas to H2 storage capacity uses the energy conversion factor of 0.27 [22], and oil to H2 storage capacity converts oil volume to equivalent hydrogen mass based on specific reservoir temperatures and pressures [86]. Equivalent H2 storage is recorded using the gas/condensate leg (g), the oil leg (o), or both (g + o); GDT = gas down-to. Illustrations, input data and references used for estimating H2 storage capacity in large prospects is provided in supplementary data file 3. 

|Site Name|Reservoir Interval|Hydrocarbon|Gas reserves|Oil reserves|Equiv. H2|Comment|
|---|---|---|---|---|---|---|
|||Phase|(TWh)|(TWh)|storage<br>(TWh)||
|**Large prospe**|**cts requiring smaller c**|**ompartments, deta**|**ils in italics**|||i|
|**Maui (A& **|Mangahewa,|Gas/cond, oil|1280|375|369 (g+ o)|Full field reserves/storage for 3 reservoirs, 2 sites.|
|**B)**|Kaimiro, Farewell||||||
|**_Maui A_**|_Kaimiro_|_Gas/cond_|_ND_|_ND_|_21.6 (g)_|_Petroleum reserves for individual reservoirs based on legacy reserve_<br>_i_|
|**_Maui B_**|_Kaimiro_|_Gas cap oil leg_|_ND_|_ND_|_17.8 (g)_|_estimates; seismic mapping from open file datasets used to demonstrate_|
|**_Maui B_**|_Farewell_|_Oil_|_ND_|_ND_|_<22 (o)_|_i_<br>_smaller closures possible in structure (GDT); Maui A, B at Mangahewa level_<br>_are considered too big._<br>i|
|**Pohokura**|Mangahewa<br>(U1–U4)|Gas/cond|376|90|108 (g+ o)|Full field reserves/storage for 4 intervals (U1–U4).|
|**_Pohokura_**|_Mangahewa (U1)_|_Gas/cond_|_ND_|_ND_|_47.8 (g)_|_U1 reserves based on operator reservoir model; smaller closure possible in_<br>_structure (GDT)._<br>i|
|**Kapuni**|Mangahewa (K1A-<br>K1E)|Gas/cond|333|108|99 (g+o)|Full field reserves/storage for several intervals (K1A-K1E).|
|**_Kapuni_**|_Mangahewa (K1A)_|_Gas/cond_|_ND_|_ND_|_19.7 (g)_|_K1A reserves based on operator production data; smaller closure possible in_<br>_structure (GDT)._<br>i|
|**Kupe**|Farewell|Gas leg oil rim|162|34|46 (g+o)|Full field reserves/storage.|
|**_Kupe_**|_Farewell_|_Gas/cond_|_ND_|_ND_|_36 (g)_|_Storage in gas leg with smaller closure possible in structure (GDT)._<br>i|
|**McKee**|McKee|Gas cap oil leg|71.7|79.8|23 (g+o)|Full field reserves/storage.|
|**_McKee_**|_McKee_|_Gas/cond_|_ND_|_ND_|_19.4 (g)_|i<br>_Storage in gas leg with smaller closure possible in separate fault blocks._|
|**H2 storage in**|**gas/condensate prosp**|**ects**|||||
|**Ahuroa**|Tariki|Gas/cond|5|N/A|1.4 (g)|Reserves data from natural gas storage (FirstGas).|
|**Tariki**|Tariki|Gas/cond|17.9|N/A|4.8 (g)||
|**H2 storage in**|**oil prospects± gas**|||||i|
|**Tui area**|Farewell|Oil|N/A|62.6|3.7 m (o)|Combined reserves data for 3 fields (Tui-Amokura-Pateke); est. 1–2<br>TWh H2storage each.<br>i|
|**Maari**|Moki|Oil, minor gas|16.1|80.9|3.1 (g+ o)|Combined field reserves data for Maari and Manaia with separate|
|**Maari**|Mangahewa|Gas cap oil leg|||2.9 (g+ o)|i<br>estimations for equivalent H2storage based on legacy reserves|
|**Manaia**|Mangahewa|Oil, minor gas|||1.2 (g+ o)|estimates.|
|**Kaimiro/**|<br>Mount Messenger|<br>Oil, minor gas|18.6|16|<br>5.6 (g+ o)|Several small pools with combined reserves data; Kaimiro and Ngatoro|
|**Ngatoro**||||||represent largest pools.|
|**Cheal/**|Urenui&MM|Oil and gas|2.8|8.9|1.1 (g+ o)||
|**Cheal E**|||||||
|**H2 storage in**|**undeveloped sites**||||||
|**Karewa**|Mangaa|Non-com gas|(43.1)|ND|11.6 (g)|Contingent gas reserves from MBIE[84].|
|**Kapuni**|Moki|Aquifer|(14.6)|ND|3.9|Unrisked reserves from operator estimates.|

for oil and some injectivity and containment parameters. Once parameters for UHS are better understood and parameter weightings have been revised to reflect this, bias in the decision-making process is expected to reduce, and we suggest that the matrix approach could provide the user with a robust method of ranking multiple sites. A paired comparison matrix [35] may be applicable in future prospectivity analysis to further reduce potential biases of an expert panel, although that process also relies on expert opinion. Notably, there remains significant uncertainty with oil reservoirs, particularly relating to seal competence (for H2), contamination issues (e.g., potential for methane production), and high residual fluids (see section 6.3). Pending further research, this could effectively disqualify these reservoirs as potential UHS sites. Furthermore, it is recommended that deep aquifer options, in addition to depleted reservoirs, need to be explored further once UHS technology in porous media has been proven. Given their specific requirements and profiles, it is suggested that different matrix parameters and criteria might be applicable for the three storage classes (i.e., gas, oil, and aquifer), and that future prospectivity analysis could result in three separate ranking profiles. 

# _6.2. Summary of outcomes_ 

We have assessed and ranked potential subsurface reservoirs in A-NZ for hydrogen storage based on publicly accessible data and geological interpretations, criteria-driven processes, and current understanding of 

the technical requirements for UHS. We suggest that storage sites in A- NZ should comprise a volume that is suitable for the size of the market (~0.5–10 TWh working gas), be deep enough to ensure sufficient H2 densities and to prevent excessive production issues ( _>_ 1 km), and have adequate reservoir quality for injectivity/withdrawal (average permeability _>_ 5 mD). Site-screening criteria that disqualified the largest number of potential UHS sites was minimum size, but this could be revised to a lower value if a storage development project was proven economic at a smaller scale. Once site-screening conditions have been satisfied, the most important parameters for suitable storage sites are thought to be related to H2 containment and data availability/quality. However, the matrix ranking of sites is highly sensitive to weightings assigned to different parameters, and further work is required to determine the relative importance of different factors for UHS (see section 6.3). 

Depleted hydrocarbon reservoirs are considered to represent the best options for developing commercial UHS projects, with depleted gas (or gas cap) reservoirs favoured over oil reservoirs, and Paleogene reservoirs favoured over younger reservoirs. Notably, some of the Miocene sites have top reservoir depths between 1 and 1.5 km and their inclusion as potential UHS sites is dependent upon the minimum depth criterion ( _>_ 1 km). Some of these sites would be eliminated if the minimum depth was revised to 1.5 km, making it consistent with the recommendation of Hassanpouryouzband et al. [87] to maximise H2 storage potential. Relatively low ranking of the single deep aquifer option suggests that 

1480 

_K.E. Higgs et al.                                                                                                                                                                                                                                 International Journal of Hydrogen Energy 71 (2024) 1468–1485_ 

_International Journal of Hydrogen Energy 71 (2024) 1468–1485_ 

**Fig. 6.** Spider plots comparing parameter values used in the decision matrix approach to site selection. **A-E)** top five sites, deeper Paleogene reservoirs; **F-J)** bottom five sites, shallow Miocene reservoirs. Values from least favourable for UHS (0 at centre of plots) to most favourable (3 at outside edge of plots). Refer to Table 3 for parameter details and criteria. 

1481 

_International Journal of Hydrogen Energy 71 (2024) 1468–1485_ 

_K.E. Higgs et al._ 

this type of play is not advisable in the short-term, and aquifer sites will inevitably require significant technical work programmes prior to consideration. 

Results from our prospectivity analysis are sensitive to capacity estimates, with reasonable estimates determined from recoverable reserves data and/or an evaluation of publicly available legacy data. Estimated H2 storage capacity and a summary of input data used for qualifying sites is shown in Table 6, with higher uncertainties associated with capacity estimates derived from oil reservoirs, undeveloped sites, or combined fields. The maximum storage for Ahuroa (Ahuroa, Tariki) plus Maui A (Maui A, Kaimiro), which represent the top two candidates for demonstrating technical suitability of UHS in Taranaki, is 23 TWh H2 (supplementary data file 2, Table 6). However, capacity estimates used in the current study assume total H2 storage (i.e., working gas + cushion gas). Since the cushion gas does not participate in injection/withdrawal cycles it needs to be excluded from the more detailed site characterisation work. The amount of cushion gas will be site-specific and may vary from ~20% to 80% of the total gas storage volume. It will also be dependent on gas composition. Assuming a 50% cushion gas requirement, total estimated H2 storage capacity across the top two sites exceeds the predicted 2050 domestic demand of ~5 TWh H2 storage (refer section 3.1). 

It has been pointed out by several authors that a conflict of interest may arise for sites with potential for underground storage of natural gas, H2, and/or CO2 [21,99]. In the case of A-NZ, the Ahuroa Field is currently being used for natural gas storage and, perhaps unsurprisingly, it also has favourable geological characteristics for UHS (Tables 2 and 4). While a change of use from UGS to UHS would need to make commercial sense to those with a stake in the field, it would conform with the government directive to achieve a net-zero carbon economy by 2050 [73]. Some other storage sites in A-NZ are very large (e.g., Maui A, Mangahewa; Kapuni, K3E Mangahewa), and much larger than would be required for seasonal UHS. As such, it might be more appropriate for these very large sites to be used for CCS or, if reservoirs are compartmentalised, for simultaneous use with both H2 and CO2 storage. These operational issues (including potential gas leakage, contamination, and well placement) would need to be resolved as part of site characterisation work. 

The current study has only considered the technical suitability of a site for UHS, yet economic, sociocultural, and regulatory factors can make a site unsuitable for development. Conversion of current infrastructure for hydrogen gas would be significant, and additionally, in A- NZ, the pore space rights and ownership, access and permission would need to be carefully assessed and discussed with respect to the indigenous population. Alongside improving our technical understanding of UHS systems, criteria must therefore be developed to evaluate regulatory, political, socio-economic, and environmental characteristics. 

# _6.3. Research directions_ 

Previous work has shown that while there is significant experience with UGS and CCS technology, there is limited understanding of the requirements for UHS. The specificity of H2 requires significant research and testing before UHS technology is implemented on a commercial scale [21]. This is also highlighted by our study, with both approaches to prospectivity analysis dependent on, and restricted by, the current knowledge base. Criteria used in the current work will need to be revised in light of additional research into H2 behaviour, in particular how it affects the seal and its potential reactivity with native pore fluid and mineral phases. Ultimately this will lead to a revision of applied criteria and weightings, thereby reducing risk in the selection process and providing the most technically suitable sites for UHS technology. 

Some of the key uncertainties in A-NZ site evaluation that need addressing are listed below: 

- H2 reactivity: experimental work to investigate the impact of mineralogy and microbes on H2 reactivity, losses and contamination. Kapuni, K1A Mangahewa is a potential site to investigate the effect of CO2, coals and carbonates on UHS. Sulphide abundance in A-NZ rocks is uncertain, but work is recommended to address what amount of sulphide will render a site unsuitable for UHS; sulphates have not been documented in Taranaki but may need to be addressed in other areas. 

- H2 mobility through seal: experimental and modelling work to investigate the sealing requirements for H2 properties, and the effect of multiple injection and withdrawal cycles where fluctuations in the reservoir pressure may result in induced effects on seal and reservoir. Examples may include seals with facies heterogeneities (e.g., Pohokura), fault-seal (e.g., Kupe), highly faulted structures (e.g., McKee), and/or an oil-only reservoir (e.g., Maui B, Farewell). 

- Cushion gas composition: experimental work and dynamic modelling with gas mixtures (e.g., H2–CH4, H2–CO2) should help us to assess gas storage site conversions, predict the degree of gas mixing (contamination), and make recommendations for cushion gas com- 

positions, which are likely to be site-specific. 

- Volume of residual hydrocarbons: experimental and modelling work to investigate the impact of residual hydrocarbons on injectivity/ withdrawal. Examples would include fields where significant oil is likely to remain in place following hydrocarbon production, or the uncommercial site Karewa, which has unknown gas reserves that would need displacing or removing prior to hydrogen injection. If results conclude that residual petroleum will be significantly detrimental to hydrogen purity, we suggest that an additional parameter of original hydrocarbon in place should be included in the prospectivity analysis. Time to petroleum depletion may also need to be considered for site prospectivity in other areas. 

- Study of aquifer recharge: modelling the displacement of brine relative to reservoir heterogeneity and transport properties of H2 and how it affects injection and withdrawal rates and scales. Good examples would be Maui A and B, Kaimiro, which have active aquifer charge in high-permeability zones, or Ahuroa, which has recorded aquifer ingress when operating at low pressures. 

- Minimum depth criterion: experimental and modelling work to determine the effect of multiple injection and withdrawal cycles on poorly consolidated Miocene reservoirs with depths 1–1.5 km. 

- Induced seismicity: an investigation into the historical seismicity of Taranaki with fault mapping and geomechanical studies to predict the degree of induced seismicity likely to result from H2 injection. Knowledge on managing and derisking seismicity due to fluid injection would be drawn from CCS studies [100,101]. For Taranaki, a regional network of seismometers (GeoNet) provides a unique opportunity for studying historical seismicity, and if UHS is implemented, this could be used with a local seismic array for seismicity monitoring. 

The geological feasibility of a site represents just one step in the process of UHS site selection, and there is a need to integrate results from the technical study with social, economic, political, and regulatory criteria before an operator will consider moving towards a pilot project. It is therefore recommended that similar appraisal work is undertaken with the facility operators and local community with any future site assessment also considering storage capacity in relation to economics and location of renewable energy production. 

# **7. Conclusions** 

Potential hydrogen storage sites have been assessed to determine the most promising locations for demonstrating technical feasibility of underground hydrogen storage (UHS) in Taranaki Basin, Aotearoa New Zealand. Two different methods have been undertaken. The decision tree provides a visual and simplified approach to site evaluation, 

1482 

_International Journal of Hydrogen Energy 71 (2024) 1468–1485_ 

_K.E. Higgs et al._ 

without the need for input from other sites. As such it can be quickly applied and highlights potential issues that require further investigation. The weighted decision matrix approach uses a criteria-driven workflow and scoring system that provides a more robust method of ranking sites once parameters are better understood. Both methodologies can easily be adapted with revised criteria for regions that have different technical challenges. 

We consider that the most important parameters for assessing UHS feasibility are storage capacity, reservoir depth, and parameters that affect injectivity and containment. Provisional recommended criteria for UHS in Taranaki Basin are H2 capacity (working gas + cushion gas) in the range 1–20 TWh, minimum top reservoir depth of 1 km, and minimum average permeability of 5 mD. Depleted oil and gas reservoirs are considered the best options for demonstrating UHS technology given their established datasets, proven reservoirs, seals and traps, and existing infrastructure. Depleted gas reservoirs (and gas caps overlying oil legs) are favoured over depleted oil reservoirs on the premise that caprocks proven to contain natural gas are more likely to contain H2 than caprocks for liquid hydrocarbons. 

Individual UHS sites in Taranaki Basin offer a wide range of H2 storage capacities (≪1 TWh to ≫20 TWh). Based on our assessment, we suggest that commercial development of one or two depleted fields (or compartments thereof) would provide sufficient storage capacity for UHS development in A-NZ by 2050 ( _>_ 5 TWh). Highest ranking sites are all from Paleogene reservoirs, which are typically fine-to coarse-grained, quartzo-feldspathic sandstones that were deposited during the rift to passive margin and earliest contractional stages of basin history. The lowest ranking sites are mostly from the finer grained, labile, clay-rich, small and typically compartmentalised Miocene reservoirs that were deposited in an active margin and are relatively poorly lithified. 

Prospectivity analysis has highlighted several parameters that need to be better understood before UHS sites can be identified. Our studies suggest that the main risks requiring further research include: a) potential losses and contamination of stored gas associated with H2 reactivity (with CO2, coal etc.), and b) potential leakage of H2 through seals due to the higher diffusivity, lower viscosity and lower density of hydrogen compared to methane. We aim to select sites for characterisation that can address some of these outstanding risks, with subsequent revision of parameters, weightings, and criteria to reflect hydrogen storage requirements. Ultimately this work will help in developing a consistent methodology for UHS prospectivity analysis to ensure safe, secure, and economic deployment of hydrogen storage that can be applied globally. 

# **Declaration of competing interest** 

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper. 

# **Acknowledgements** 

Funding for this work was received from the Ministry of Business, Innovation and Employment (MBIE) to the University of Canterbury under the 2022 Endeavour Fund Research Programmes (grant number UOCX2207). We would like to acknowledge New Zealand Petroleum & Minerals (NZP&M), Ministry of Business, Innovation and Employment (MBIE) for their open-file data repository https://geodata.nzpam.govt. nz/, and FirstGas Ltd. for their support of the project. We acknowledge all colleagues who have contributed data to this study, and especially thank Mac Beggs for his helpful review of this paper. 

# **Appendix A. Supplementary data** 

Supplementary data to this article can be found online at https://doi. org/10.1016/j.ijhydene.2024.05.098. 

# **References** 

- [1] Schmidt O, Melchior S, Hawkes A, Staffell I. Projecting the future levelized cost of electricity storage technologies. Joule 2019;3:81–100. https://doi.org/10.1016/j. joule.2018.12.008. 

- [2] Olabi AG, Bahri AS, Abdelghafar AA, Baroutaji A, Sayed ET, Alami AH, Rezk H, Abdelkareem MA. Large-vscale hydrogen production and storage technologies: current status and future directions. Int J Hydrogen Energy 2021;46:23498–528. https://doi.org/10.1016/j.ijhydene.2020.10.110. 

- [3] Tarkowski R, Uliasz-Misiak B. Towards underground hydrogen storage: a review of barriers. Renew Sustain Energy Rev 2022;162:112451. https://doi.org/ 10.1016/j.rser.2022.112451. 

- [4] Peecock A, Edlmann K, Mouli-Castillo J, Martinez-Felipe A, McKenna R. Mapping hydrogen storage capacities of UK offshore hydrocarbon fields and exploring potential synergies with offshore wind. In: Miocic JM, Heinemann N, Edlmann K, Alcalde J, Schultz RA, editors. Enabling secure subsurface storage in future energy systems: geological Society, vol. 528. London: Special Publications; 2023. https://doi.org/10.1144/SP528-2022-40. 

- [5] IRENA. Green hydrogen cost reduction: scaling up electrolysers to meet the 1.5<sup>◦</sup> C climate goal. Int Renew Energy Agency 2020. https://www.irena.org/pu blications/2020/Dec/Green-hydrogen-cost-reduction. 

- [6] Heinemann N, Alcalde J, Miocic JM, Hangx SJT, Kallmeyer J, OstertagHenning C, Hassanpouryouzband A, Thaysen EM, Strobel GJ, SchmidtHattenberger C, Edlmann K, Wilkinson M, Bentham M, Haszeldine RS, Carbonell R, Rudloff A. Enabling large scale hydrogen storage in porous media – the scientific challenges. Energy Environ Sci 2021;14:853–64. https://doi.org/ 10.1039/D0EE03536J. 

- [7] Langmi HW, Engelbrecht N, Modisha PM, Bessarabov D. Hydrogen storage. In: Smolinka T, Garche J, editors. Electrochemical power sources: Fundamentals, systems, and applications; hydrogen production by water electrolysis. Elsevier; 2022. p. 455–85. https://doi.org/10.1016/B978-0-12-819424-9.00006-9. 

- [8] Simon J, Ferriz AM, Correas LC. HyUnder – hydrogen underground storage at large scale: case study Spain. Energy Proc 2015;73:136–44. https://doi.org/ 10.1016/j.egypro.2015.07.661. 

- [9] Heinemann N, Booth MG, Haszeldine RS, Wilkinson M, Scafidi J, Edlmann K. Hydrogen storage in porous geological formations - onshore play opportunities in the Midland Valley. Int J Hydrogen Energy 2018;43:20861–74. https://doi.org/ 10.1016/j.ijhydene.2018.09.149. 

- [10] Chapman A, Itaoka K, Hirose K, Davidson FT, Nagasawa K, Lloyd AC. A review of four case studies assessing the potential for hydrogen penetration of the future energy system. Int J Hydrogen Energy 2019;44:6371–82. https://doi.org/ 10.1016/j.ijhydene.2019.01.168. 

- [11] Zivar D, Kumar S, Foroozesh J. Underground hydrogen storage: a comprehensive review. Int J Hydrogen Energy 2021;46:23436–62. https://doi.org/10.1016/j. ijhydene.2020.08.138. 

- [12] Muhammed NS, Haq B, Al-Shehri D, Al-Ahmed A, Rahman MM, Zaman E. A review on underground hydrogen storage: Insight into geological sites, influencing factors and future outlook. Energy Rep 2022;8:461–99. https://doi. org/10.1016/j.egyr.2021.12.002. 

- [13] IEA Hydrogen TCP-Task 42. Underground Hydrogen Storage: technology monitor report. https://www.ieahydrogen.org/task/task-42-underground-hydrogen-s torage/; 2023. 

- [14] Lord AS, Kobos PH, Borns DJ. Geologic storage of hydrogen: scaling up to meet city transportation demands. Int J Hydrogen Energy 2014;39:15570–82. https:// doi.org/10.1016/j.ijhydene.2014.07.121. 

- [15] Bauer S, Pfeiffer T, Boockmeyer A, Dahmke A, Beyer Ch. Quantifying induced effects of subsurface renewable energy storage. Energy Proc 2015;76:633–41. https://doi.org/10.1016/j.egypro.2015.07.885. 

- [16] Juez-Larr´e J, Machado CG, Groeneberg RM, Belfroid SSPC, Yousefim SH. A detailed comparative performance study of underground storage of natural gas and hydrogen in The Netherlands. Int J Hydrogen Energy 2023;74:28843–68. https://doi.org/10.1016/j.ijhydene.2023.03.347. 

- [17] Heinemann N, Scafidi J, Pickup G, Thaysen EM, Hassanpouryouzband A, Wilkinson M, et al. Hydrogen storage in saline aquifers: the role of cushion gas for injection and production. Int J Hydrogen Energy 2021;46:39284–96. https://doi. org/10.1016/j.ijhydene.2021.09.174. 

- [18] Kanaani M, Sedaee B, Asadian-Pakfar M. Role of cushion gas on underground hydrogen storage in depleted oil reservoirs. J Energy Storage 2022;45:103783. https://doi.org/10.1016/j.est.2021.103783. 

- [19] Abdellatif M, Hashemi M, Azizmohammadi S. Large-scale underground hydrogen storage: Integrated modeling of reservoir-wellbore system. Int J Hydrogen Energy 2023;48:19160–71. https://doi.org/10.1016/j.ijhydene.2023.01.227. 

- [20] Tarkowski R. Perspectives of using the geological subsurface for hydrogen storage in Poland. Int J Hydrogen Energy 2017;42:347–55. https://doi.org/10.1016/j. ijhydene.2016.10.136. 

- [21] Tarkowski R, Uliasz-Misiak B, Tarkowski P. Storage of hydrogen, natural gas, and carbon dioxide – geological and legal conditions. Int J Hydrogen Energy 2021;46: 20010–22. https://doi.org/10.1016/j.ijhydene.2021.03.131. 

- [22] Ennis-King J, Michael K, Strand J, Sander R, Green C. Underground storage of hydrogen: mapping out the options for Australia. Future Fuels CRC 2021. http s://www.futurefuelscrc.com/project/underground-storage-of-hydrogen-mappin g-out-the-options-for-australia-rp1-1-04/. 

- [23] Pan B, Yin X, Ju Y, Iglauer S. Underground hydrogen storage: influencing parameters and future outlook. Adv Colloid Interface Sci 2021;294:102473. https://doi.org/10.1016/j.cis.2021.102473. 

1483 

_K.E. Higgs et al.                                                                                                                                                                                                                                 International Journal of Hydrogen Energy 71 (2024) 1468–1485_ 

- [24] Nicol A, Dempsey D, Yates E, Higgs K, Beggs M, Adam L. Underground hydrogen storage in the Taranaki region, New Zealand. Rep Firstgas 2022. https://hdl.hand le.net/10092/106667. 

- [25] Ma J, Kühn M, Nakaten N. Power-to-gas based subsurface energy storage: a review. Renew Sustain Energy Rev 2018;97:478–96. https://doi.org/10.1016/j. rser.2018.08.056. 

- [26] Callas C, Saltzer SD, Davis JS, Hashemi SS, Kovscek AR, Okoroafor ER, Wen G, Zoback MD, Benson SM. Criteria and workflow for selecting depleted hydrocarbon reservoirs for carbon storage. Appl Energy 2022;324. https://doi. org/10.1016/j.apenergy.2022.119668. 

- [27] Esfandyari H, Sarmadivaleh M, Esmaeilzadeh F, Ali M, Iglauer S, Keshavarz A. Experimental evaluation of rock mineralogy on hydrogen-wettability: implications for hydrogen geostorage. J Energy Storage 2022;52:104866. https:// doi.org/10.1016/j.est.2022.104866. 

- [28] Lemieux A, Shkarupin A, Sharp K. Geologic feasibility of underground hydrogen storage in Canada. Int J Hydrogen Energy 2020;45:32243–59. https://doi.org/ 10.1016/j.ijhydene.2020.08.244. 

- [29] Ciotta M, Tassinari C, Guilherme L, Zacharias L, van de Zwaan B, Peyerl D. Hydrogen storage in depleted offshore gas fields in Brazil: potential and implications for energy security. Int J Hydrogen Energy 2023;48:39967–80. https://doi.org/10.1016/j.ijhydene.2023.08.209. 

- [30] Lemieux A, Sharp K, Shkarupin A. Preliminary assessment of underground hydrogen storage sites in Ontario, Canada. Int J Hydrogen Energy 2019;44: 15193–204. https://doi.org/10.1016/j.ijhydene.2019.04.113. 

- [31] Bai M, Song K, Sun Y, He M, Li Y, Sun J. An overview of hydrogen underground storage technology and prospects in China. J Petrol Sci Eng 2014;124:132–6. https://doi.org/10.1016/j.petrol.2014.09.037. 

- [32] Henkel S, Pudlo D, Gaupp R. Research sites of the H2STORE project and the relevance of lithological variations for hydrogen storage at depths. Energy Proc 2013;40:25–33. https://doi.org/10.1016/j.egypro.2013.08.004. 

- [33] Safari A, Zeng L, Nguele R, Sugai Y, Sarmadivaleh M. Review on using the depleted gas reservoirs for the underground H2 storage: a case study in Niigata prefecture, Japan. Int J Hydrogen Energy 2023;48:10579–602. https://doi.org/ 10.1016/j.ijhydene.2022.12.108. 

- [34] Netherlands Enterprise Agency. The effects of hydrogen injection in natural gas networks for the Dutch underground storages. In: Final report. Commissioned by the ministry of economic affairs; 2017. Publication number: RVO-079-1701/RPDUZA. 

- [35] Lewandowska-Smierzchalska J, Tarkowski R, Uliasz-Misiak B. Screening and ranking framework for underground hydrogen storage site selection in Poland. Int J Hydrogen Energy 2018;43:4401–14. https://doi.org/10.1016/j. ijhydene.2018.01.089. 

- [36] Carneiro JF, Matos CR, van Gessel S. Opportunities for large-scale energy storage in geological formations in mainland Portugal. Renew Sustain Energy Rev 2019; 99:201–11. https://doi.org/10.1016/j.rser.2018.09.036. 

- [37] Mouli-Castillo J, Heinemann N, Edlmann K. Mapping geological hydrogen storage capacity and regional heating demands: an applied UK case study. Appl Energy 2021;283:116348. https://doi.org/10.1016/j.apenergy.2020.116348. 

- [38] Scafidi J, Wilkinson M, Gilfillan SMV, Heinemann N, Haszeldine RS. A quantitative assessment of the hydrogen storage capacity of the UK continental shelf. Int J Hydrogen Energy 2021;46:8629–39. https://doi.org/10.1016/j. ijhydene.2020.12.106. 

- [39] HyUSpre. 2023. https://www.hyuspre.eu/. [Accessed 5 March 2024]. 

- [40] Hystories. 2023. https://hystories.eu/. [Accessed 5 March 2024]. 

- [41] Le Gallo Y. Assessment of hydrogen storage capacity in porous media at the European scale. Int J Hydrogen Energy 2024;58:400–8. https://doi.org/10.1016/ j.ijhydene.2024.01.158. 

- [42] Yates E, Bischoff A, Beggs M, Jackson N. Green hydrogen geostorage in Aotearoa – New Zealand. Elemental group for Ministry of business. Innovat Employment 2021. https://doi.org/10.13140/RG.2.2.36303.82083. 

- [43] Taylor JB, Alderson JEA, Kalyanam KM, Lyle AB, Phyllips LA. Technical and economic assessment of methods for the storage of large quantities of hydrogen. Int J Hydrogen Energy 1986;11:5–22. https://doi.org/10.1016/0360-3199(86) 90104-7. 

- [44] Anderson PD, Norgaard PH, Olesen MH, Tanner AN. Final report of the roads2hycom project: fuel cells and hydrogen in a sustainable energy economy. https://orbit.dtu.dk/en/publications/final-report-of-the-roads2hycom-projectfuel-cells-and-hydrogen-i; 2009. 

- [45] Meier K. Hydrogen production with sea water electrolysis using Norwegian offshore wind energy potentials. Int J Energy Environ Eng 2014;104. https://doi. org/10.1007/s40095-014-0104-6. 

- [46] Edlmann K, Haszeldine S, Heinemann N, Hassanpouryouzband A, Thaysen E. Hystorpor: unlocking the UK’s hydrogen storage. https://blogs.ed.ac.uk/hystorp or/2021/06/07/hystorpor-unlocking-the-uks-hydrogen-storage/. [Accessed 26 January 2024]. 

- [47] Feldmann F, Hagemann B, Ganzer L, Panfilov M. Numerical simulation of hydrodynamic and gas mixing processes in underground hydrogen storages. Environ Earth Sci 2016;75:1165. https://doi.org/10.1007/s12665-016-5948-z. 

- [48] Sainz-Garcia A, Abarca E, Rubi V, Grandia F. Assessment of feasible strategies for seasonal underground hydrogen storage in a saline aquifer. Int J Hydrogen Energy 2017;42:16657–66. https://doi.org/10.1016/j.ijhydene.2017.05.076. 

- [49] Pfeiffer WT, Bauer S. Comparing simulations of hydrogen storage in a sandstone formation using heterogeneous and homogenous flow property models. Petrol Geosci 2018;25:325–36. https://doi.org/10.1144/petgeo2018-101. 

- [50] Lubon K, Tarkowski R. Numerical simulation of hydrogen injection and withdrawal to and from a deep aquifer in NW Poland. Int J Hydrogen Energy 2020;45:2068–83. https://doi.org/10.1016/j.ijhydene.2019.11.055. 

- [51] Lysyy M, Fernø M, Ersland G. Seasonal hydrogen storage in a depleted oil and gas field. Int J Hydrogen Energy 2021;46:25160–74. https://doi.org/10.1016/j. ijhydene.2021.05.030. 

- [52] Huang T, Moridis GJ, Blasingame TA, Abdulkader AM, Yan B. Compositional reservoir simulation of underground hydrogen storage in depleted gas reservoirs. Int J Hydrogen Energy 2023;48:36035–50. https://doi.org/10.1016/j. ijhydene.2023.05.355. 

- [53] Okoroafor ER, Saltzer SD, Kovscek AR. Toward underground hydrogen storage in porous media: reservoir engineering insights. Int J Hydrogen Energy 2022: 33781–802. https://doi.org/10.1016/j.ijhydene.2022.07.239. 

- [54] Muhammed NS, Haq B, Abdullah D, Shehri A. Hydrogen storage in depleted gas reservoirs using nitrogen cushion gas: a contact angle and surface tension study. Int J Hydrogen Energy 2023;48:38782–807. https://doi.org/10.1016/j. ijhydene.2023.06.208. 

- [55] Truche L, Jodin-Caumon MC, Lerouge C, Berger G, Mosser-Ruck R, Giffaut E, et al. Sulphide mineral reactions in clay-rich rock induced by high hydrogen pressure. Application to disturbed or natural settings up to 250 ◦C and 30 bar. Chem Geol 2013;351:217–28. https://doi.org/10.1016/j.chemgeo.2013.05.025. 

- [56] Henkel S, Pudlo D, Werner L, Enzmann F, Reitenbach V, Albrecht D. Mineral reactions in the geological underground induced by H2 and CO2 injections. Energy Proc 2014;63:8026–35. https://doi.org/10.1016/j.egypro.2014.11.839. 

- [57] Hemme C, Van Berk W. Hydrogeochemical modelling to identify potential risks of underground storage in depleted gas fields. Appl Sci 2018;2282. https://doi.org/ 10.3390/app8112282. 

- [58] Yekta AE, Pichavant M, Audigane P. Evaluation of geochemical reactivity of hydrogen in sandstone: application to geological storage. Appl Geochem 2018;95: 182–94. https://doi.org/10.1016/j.apgeochem.2018.05.021. 

- [59] Rezaee R. Assessment of natural hydrocarbon systems in Western Australia. Int J Hydrogen Energy 2021;46:33068–77. https://doi.org/10.1016/j. ijhydene.2021.07.149. 

- [60] Bo Z, Zeng L, Chen Y, Xie Q. Geochemical reactions-induced hydrogen loss during underground hydrogen storage in sandstone reservoirs. Int J Hydrogen Energy 2021;46:19998–20009. https://doi.org/10.1016/j.ijhydene.2021.03.116. 

- [61] Labus K, Tarkowski R. Modeling hydrogen – rock – brine interactions for the Jurassic reservoir and caprocks from Polish lowlands. Int J Hydrogen Energy 2022;47:10947–62. https://doi.org/10.1016/j.ijhydene.2022.01.134. 

- [62] Ziemianski PP, Derkowski A. Structural and textural control of high-pressure hydrogen adsorption on expandable and non-expandable clay minerals in geologic conditions. Int J Hydrogen Energy 2022;47:28794–805. https://doi.org/ 10.1016/j.ijhydene.2022.06.204. 

- [63] Dopffel N, Jansen S, Gerritse J. Microbial side effects of underground hydrogen storage – knowledge gaps, risks and opportunities for successful implementation. Int J Hydrogen Energy 2021;46:8594–606. https://doi.org/10.1016/j. ijhydene.2020.12.058. 

- [64] Thaysen EM, McMahon S, Strobel GJ, Butler IB, Ngwenya BT, Heinemann N, Wilkinson M, Hassanpouryouzband A, McDermott CI, Edlmann K. Estimating microbial growth and hydrogen consumption in hydrogen storage in porous media. Renew Sustain Energy Rev 2021;151:111481. https://doi.org/10.1016/j. rser.2021.111481. 

- [65] Amid A, Mignard D, Wilkinson M. Seasonal storage of hydrogen in a depleted natural gas reservoir. Int J Hydrogen Energy 2016;41:5549–58. https://doi.org/ 10.1016/j.ijhydene.2016.02.036. 

- [66] Flesch S, Pudlo D, Albrecht D, Jacob A, Enzmann F. Hydrogen underground storage – petrographic and petrophysical variations in reservoir sandstones from laboratory experiments under simulated reservoir conditions. Int J Hydrogen Energy 2018;43:20822–35. https://doi.org/10.1016/j.ijhydene.2018.09.112. 

- [67] Hassannayebi N, Azizmohammadi S, De Lucia M, Ott H. Underground hydrogen storage: application of geochemical modelling in a case study in the Molasse Basin, Upper Austria. Environ Earth Sci 2019;78:1–13. https://doi.org/10.1007/ s12665-019-8184-5. 

- [68] Shi Z, Jessen K, Tsotsis TT. Impacts of the subsurface storage of natural gas and hydrogen mixtures. Int J Hydrogen Energy 2020;45:8757–73. https://doi.org/ 10.1016/j.ijhydene.2020.01.044. 

- [69] Bo Z, Zeng L, Chen Y, Xie Q. Geochemical reactions-induced hydrogen loss during underground hydrogen storage in sandstone reservoirs. Int J Hydrogen Energy 2021;46:19998–20009. https://doi.org/10.1016/j.ijhydene.2021.03.116. 

- [70] Harati S, Gomari SR, Ramegowda M, Pak T. Multi-criteria site selection workflow for geological storage of hydrogen in depleted gas fields: a case for the UK. Int J Hydrogen Energy 2024;51:143–57. https://doi.org/10.1016/j. ijhydene.2023.10.345. 

- [71] Safari A, Sugai Y, Sarmadivaleh M, Imai M. Screening and ranking Japanese gas fields for underground H2 storage potential: impact of the reservoir drive mechanism. J Energy Storage 2023;70:107679. https://doi.org/10.1016/j. est.2023.107679. 

- [72] Ministry of Business. Innovation & employment. Electricity statistics. 2022. htt ps://www.mbie.govt.nz/building-and-energy/energy-and-natural-resource s/energy-statistics-and-modelling/energy-statistics/electricity-statistics/. 

- [73] Ministry of Business, Innovation & Employment. A vision for hydrogen in New Zealand: green paper. https://www.mbie.govt.nz/building-and-energy/energyand-natural-resources/energy-strategies-for-new-zealand/hydrogen-in-new-ze aland/a-vision-for-hydrogen-in-new-zealand/; 2019. 

- [74] Yergin D. The new map: energy, climate, and the clash of nations. Penguin Press; 2020. ISBN 978-1594206436. 

1484 

_International Journal of Hydrogen Energy 71 (2024) 1468–1485_ 

_K.E. Higgs et al._ 

- [75] Firstgas. Bringing zero carbon gas to Aotearoa; Hydrogen feasibility study - summary report. https://cms.firstgas.co.nz/assets/Uploads/Documents/FirstgasGroup_Hydrogen-Feasibility-Study.pdf; 2021. 

- [76] King PR, Thrasher GP. Cretaceous-Cenozoic geology and petroleum systems of the Taranaki Basin, New Zealand, vol. 13. Institute of Geological and Nuclear Sciences Monograph; 1996. 

- [77] Stagpoole V, Nicol A. Regional structure and kinematic history of a large subduction back thrust: Taranaki Fault, New Zealand. J Geophys Res 2008;113: B01403. https://doi.org/10.1029/2007JB005170. 

- [78] Thrasher GP. Tectonics of the Taranaki rift. New Zealand oil exploration conference proceedings. Wellington: Ministry of Commerce; 1990. p. 124–33. 

- [79] Strogen DP, Seebeck H, Nicol A, King PR. Two-phase cretaceous-paleocene rifting in the Taranaki Basin region, New Zealand; implications for Gondwana breakup. J Geol Soc 2017;174:929–46. https://doi.org/10.1144/jgs2016-160. 

- [80] Strogen DP. (compiler). Updated paleogeographic maps for the Taranaki Basin and surrounds. GNS Sci Rep 2010/53 2011. 

- [81] King PR, Thrasher GP. Post-Eocene development of the Taranaki Basin, New Zealand: convergent overprint of a passive margin. In: Watkins JS, Zhiqiang F, McMillen KJ, editors. Geology and geophysics of continental margins, vol. 53. American Association of Petroleum Geologists Memoir; 1992. p. 93–118. https:// doi.org/10.1306/M53552C7. 

- [82] Bull S, Nicol A, Strogen D, Kroeger KF, Seebeck HS. Tectonic controls on Miocene sedimentation in the Southern Taranaki Basin and implications for New Zealand plate boundary deformation. Basin Res 2019;31:253–73. https://doi.org/ 10.1111/bre.12319. 

- [83] Higgs KE, King PR. Sandstone provenance and sediment dispersal in a complex tectonic setting: Taranaki Basin, New Zealand. Sediment Geol 2018;372:112–39. https://doi.org/10.1016/j.sedgeo.2018.05.004. 

- [84] Ministry of Business. Innovation & employment. Petrol Reserv Data 2022. https ://www.mbie.govt.nz/building-and-energy/energy-and-natural-resources/energ 

   - y-statistics-and-modelling/energy-statistics/petroleum-reserves-data/. dataset. 

- [85] Aftab A, Hassanpouryouzband A, Xie Q, Machuca LL, Sarmadivaleh M. Toward a fundamental understanding of geological hydrogen storage. Ind Eng Chem Res 2022;61:3233–53. https://doi.org/10.1021/acs.iecr.1c04380. 

- [86] Sykes R, Andrew SM, Joanne C. The New Zealand petroleum PVT database. GNS Sci Data Ser 8a 2011. https://data.gns.cri.nz/metadata/srv/eng/catalog.search #/metadata/098b2a66-b9f8-4e65-8ff8-d1e6f26ed9fe. dataset. 

- [87] Hassanpouryouzband A, Joonaki E, Edlmann K, Haszeldine RS. Offshore geological storage of hydrogen: is this our best option to achieve net-zero? ACS Energy Lett 2021;6:2181–6. https://doi.org/10.1021/acsenergylett.1c00845. 

- [88] Bear J. Dynamics of fluids in porous media. New York: Dover; 1988. ISBN 9780486656755. 

- [89] Faergestad E, editor. Formation damage: the defining series. Oilfield Review. Schlumberger; 2016. 

- [90] Wang J, Fu J, Xie J, Wang J. Quantitative characterisation of gas loss and numerical simulations of underground gas storage based on gas displacement experiments performed with systems of small-core devices connected in series. J Nat Gas Sci Eng 2020;81:103495. https://doi.org/10.1016/j. jngse.2020.103495. 

- [91] Gussow WC. Differential entrapment of oil and gas: a fundamental principle. AAPG (Am Assoc Pet Geol) Bull 1954;38:816–53. https://doi.org/10.1306/ 5CEADF11-16BB-11D7-8645000102C1865D. 

- [92] Schowalter TM. Mechanics of secondary hydrocarbon migration and entrapment. AAPG (Am Assoc Pet Geol) Bull 1979;63:723–60. https://doi.org/10.1306/ 2F9182CA-16CE-11D7-8645000102C1865D. 

- [93] Bailey WR, Underschultz J, Dewhurst DN, Kovack G, Mildren S, Raven M. Multidisciplinary approach to fault and top seal appraisal; Pyrenees-Macedon oil and gas fields, Exmouth sub-basin, Australian northwest shelf. Mar Petrol Geol 2006; 23:241–59. https://doi.org/10.1016/j.marpetgeo.2005.08.004. 

- [94] New Zealand Petroleum and Minerals Geodata Catalogue, https://geodata.nzpam .govt.nz/; [accessed 04 March 2024]. 

- [95] Davies GG, Whiteside MW, Young MS. An integrated approach to prospect evaluation. Paper presented at the SPE Offshore Europe. United Kingdom: Aberdeen; September 1991. https://doi.org/10.2118/23157-MS. 

- [96] Ward S, Adams M. PEP51150 prospectivity review report. NZP&M petroleum report 5616. 2018. https://geodata.nzpam.govt.nz/report/pr5616. 

- [97] Alcalde J, Heinemann N, James A, Bond CE, Ghanbari S, Mackay EJ, Haszeldine S, Faulkner DR, Worden RH, Allen MJ. A criteria-driven approach to the CO2 storage site selection of East Mey for the acorn project in the North Sea. Mar Petrol Geol 2021;133:105309. https://doi.org/10.1016/j. marpetgeo.2021.105309. 

- [98] Kim TW, Callas C, Saltzer SD, Kovscek AR. Assessment of oil and gas fields in California as potential CO2 storage sites. Int J Greenh Gas Control 2022;114: 103579. https://doi.org/10.1016/j.ijggc.2022.103579. 

- [99] Matos CR, Carneiro JF, Silva PP. Overview of large-scale underground energy storage technologies for integration of renewable energies and criteria for reservoir identification. J Energy Storage 2019;21:241–58. https://doi.org/ 10.1016/j.est.2018.11.023. 

- [100] Krevor S, De Coninck H, Gasda SE, Ghaleigh NS, de Gooyert V, Hajibeygi H, Juanes R, Neufeld J, Roberts JJ, Swennenhuis F. Subsurface carbon dioxide and hydrogen storage for a sustainable energy future. Nat Rev Earth Environ 2023;4 (2):102–18. https://doi.org/10.1038/s43017-022-00376-8. 

- [101] Lee KK, Ellsworth WL, Giardini D, Townend J, Ge S, Shimamoto T, Yeo IW, Kang TS, Rhie J, Sheen DH, Chang C. Managing injection-induced seismic risks. Science 2019;364:730–2. https://doi.org/10.1126/science.aax1878. 

1485 

