Systems and Methods for Enhanced Natural Hydrogen Exploration

The present disclosure relates generally to hydrogen resources and, more particularly, to systems and methods for enhanced natural hydrogen exploration using fault data.

As concerns over climate change continue to increase, there is growing interest in mitigating the effects of industrial processes, such as cement and steel production, and combustion processes utilizing fossil fuels. Hydrogen-based energy is one approach that has been suggested for mitigating the effects of carbon dioxide and other greenhouse gases. Hydrogen-based energy is carbon free. Thus, hydrogen-based industrial operation may emit fewer greenhouse gases (CHGs) by mitigating the presence of carbon during industrial procedures. Ideally, formations targeted for hydrogen processing operations are effective to collect hydrogen in a manner that is economically efficient and environmentally sound. However, typical hydrogen reservoirs are expensive to access and environmentally costly to produce.

Although current techniques for producing hydrogen energy, and for hydrogen exploration in particular, are based on technological advancements made over many years, current hydrogen exploration technology may still be ineffective to achieve ideal hydrogen production results. For example, hydrogen production in a typical formation may be expensive and environmentally tenuous. As a result, implementation of hydrogen exploration technologies may be hindered, consequently reducing CGH mitigation. Accordingly, there is an impetus to improve current hydrogen exploration technology to improve the cost and environmental impact of hydrogen production, including, for example: reducing the environmental impact of hydrogen exploration procedures, increasing the efficiency in a variety of target hydrogen production formations, increasing the throughput of hydrogen production, decreasing the cost of hydrogen production, improving the characterization of a target formation for optimal hydrogen exploration, reducing the cost and inefficiency associated with hydrogen exploration, and the like.

Consequently, there exists a need for further improvements in hydrogen exploration technology to overcome the aforementioned technical challenges and other challenges not mentioned.

Various details of the present disclosure are hereinafter summarized to provide a basic understanding. This summary is not an exhaustive overview of the disclosure and is neither intended to identify certain elements of the disclosure, nor to delineate the scope thereof. Rather, the primary purpose of this summary is to present some concepts of the disclosure in a simplified form prior to the more detailed description that is presented hereinafter.

According to an embodiment consistent with the present disclosure, a method for hydrogen exploration may include one or more steps. The method may include obtaining one or more hydrogen prospects, the one or more hydrogen prospects selected based on one or more parameters; The method may include constructing one or more fracture networks based on the one or more hydrogen prospects, the one or more fracture networks each having a source rock area value. The method may include, based on the source rock area value for each of the one or more hydrogen prospects, generating a source rock exposure value for each of the one or more prospects. The method may include outputting a production prediction based on the source rock exposure value for each of the one or more hydrogen prospects.

In another embodiment, a system for hydrogen exploration may include a memory and one or more processors. In one example, the one or more processors configured to cause the apparatus to perform one or more steps. The one or more steps may include obtaining one or more hydrogen prospects, the one or more hydrogen prospects selected based on one or more parameters. The one or more steps may include constructing one or more fracture networks based on the one or more hydrogen prospects, the one or more fracture networks each having a source rock area value. The one or more steps may include, based on the source rock area value for each of the one or more hydrogen prospects, generating a source rock exposure value for each of the one or more prospects. The one or more steps may include outputting a production prediction based on the source rock exposure value for each of the one or more hydrogen prospects.

Any combinations of the various embodiments and implementations disclosed herein can be used in a further embodiment, consistent with the disclosure. These and other aspects and features can be appreciated from the following description of certain embodiments presented herein in accordance with the disclosure and the accompanying drawings and claims.

Embodiments of the present disclosure will now be described in detail with reference to the accompanying drawing figures. Like elements in the various figures may be denoted by like reference numerals. Further, in the following detailed description, specific details are set forth in order to provide a more thorough understanding of the claimed subject matter. However, it will be apparent to one of ordinary skill in the art that the embodiments disclosed herein may be practiced without these specific details, or with details that are not described herein in the interest of clarity. Thus in some instances, well-known features have not been described in detail to avoid unnecessarily complicating the description. Additionally, it will be apparent to one of ordinary skill in the art that the scale of the elements presented in the accompanying drawing figures may vary without departing from the scope of the present disclosure.

Aspects of the present disclosure provide systems and methods for natural hydrogen exploration using fault data to identify favorable exploration targets. A natural hydrogen system is one in which component reagents for generated hydrogen may be derived from features within natural geological formations (e.g., rock formations). In certain aspects, the reactions may occur in place, whether at a certain depth within a formation or near the surface of a formation. In at least one embodiment, methods described herein may be applied during hydrogen production that targets direct production from the host rocks (e.g., white hydrogen or gold hydrogen). In at least one embodiment, methods described herein may be applied with assistance via a stimulant (e.g., orange hydrogen). Faults, fault-related fracture zones, as well as fracture zones not directly formed by faulting but affected by subsequent fault activity (henceforth all referred to as “faults” below for brevity) may serve as migration pathways in a tight source formation that would otherwise have low permeability. Faults can be sealing (e.g., cemented) or flow-enhancing. In at least one embodiment, active faulting may be identified as an agent for circulating water (e.g., natural water or injected water) that may facilitate hydrogen reactions. Active faulting may also be identified as a set of conduits for hydrogen migration towards a trap or a producer well. Active fault identification may be especially important for identifying subsurface fluid migration because, unlike sedimentary rock formations targeted for petroleum production, hydrogen source rocks may be crystalline and may have little intrinsic porosity and permeability that would allow fluid migration without fracture.

Earthquakes may be direct indicators of active faulting in recent geologic history. Earthquake data may be useful to screen for geological zones able to produce natural hydrogen. Methods provided herein may facilitate analysis of a geological fault population using earthquake data to determine conditions conducive to active flow. Methods provided herein may further allow analysis of natural hydrogen productivity potential due to geochemical hydration reactions. In many cases, hydrogen may be generated from pre-existing hydration of several source rocks within a formation. An example source rock may include a source rock containing Fe-rich (particularly Fe) minerals (e.g., olivine, pyroxene, spinel, mica, illite, chlorite, smectite, etc.), serpentinite, kimberlite, peridotite, pyroxenite, basalt/diabase, gabbro, granite (sensu lato), clay-rich sedimentary rock and rock assemblages collectively referred to as ophiolites. Ophiolites, for example, are commonly referenced in natural hydrogen literature because of the well-known hydration reaction called serpentinization. A hydration reaction requires the source rock to receive an influx of water. In many cases, water may be derived from the ocean, groundwater, modern aquifers that recharge from rainwater, from deep-seated paleo-aquifers, from injected water during engineering operations, and the like. In the hydrogen systems framework that is presented according to aspects of the present disclosure, active faults and associated fractures may serve as a conduit for water influx as well as a pathway for enhanced migration of generated hydrogen. Accordingly, the hydration reactions may only occur in certain geological locations and under favorable reaction conditions. Systems and methods provided herein may utilize fault data to identify natural hydrogen systems that occur under favorable conditions, and to predict corresponding hydrogen generation.

illustrate example formations that may be targets for hydrogen production.illustrates a schematic view of an example formation for natural hydrogen exploration. In, point 1 is a radiolysis reaction, point 2 is a serpentinization reaction, point 3 is a deep-seated fault, point 4 is a set of seeps, point 5 is a microbial consumption point, point 6 is a set of abiotic reactions, point 7 is production from a set of traps, point 8 is a direct natural hydrogen drilling, and point 9 is an enhanced hydrogen stimulation plant.illustrates a schematic view of an example formation for natural hydrogen exploration.

The formation ofrelies on faultsto transport hydrogen from a deep serpentinized source rock. Fault zonesand associated fractures in hydrogen source rocksmay be associated with natural hydrogen. Faultsand associated fractures may function as conduits that enhance fluid flow at the subsurface. In some cases, fault zonepresence has been correlated with oil and gas anomalies. In, water infiltration into iron-rich mantle rock, or to the source rock serpentinization zonemay be assumed but may not be constrained. Faultsmay enhance fluid migration into the source rock/and may facilitate gas migration to traps that are disposed above the source rock/. As depicted in, natural hydrogen may be produced in situ by direct drilling and production, or by using enhanced methods such as circulating water through a pair of injector-producer wellbores. However, even under enhanced conditions, the tight nature of the source rocks/may suggest the presence of fault enhanced migration pathways. In some cases, there may be a correlation with a deep fault zone and natural hydrogen concentration at the surface. Similarly, there may be a correlation between natural hydrogen concentration and fracture systems in certain ophiolites. Additionally, active fault zones may control fluid migration between deep source rocks and outlets at the surface.

In the process of natural hydrogen generation, a hydration reaction may occur in iron bearing rocks. The hydrogen reaction may be caused by circulation of ground water or aquifer sources. In some cases, the reaction may be represented by the simplified equation as follows:

In (1), the original rock may contain iron (II) or Fe. As this iron comes into contact with water, it oxidizes and becomes iron (III) or Fein the resultant rock. The hydrogen in the water (left side of equation) may simultaneously become reduced to free H(right side of equation). Iron (II) is the dominant driving force of this hydration reaction, and may facilitate hydrogen generation.

The compositional elements of the hydration reactions may be characteristics of a type of source rock. Embodiments of the present disclosure are described in terms of serpentinization of olivine bearing rocks such as peridotite, gabbro, and basalt, or as an assemblage of these rocks (termed ‘ophiolite’). However, it is noted that other relevant hydration reactions are possible. The mineral olivine may be formed in a solid-solution and can have large variations in chemical composition, ranging between an Mg-rich and an Fe-rich endmember. Because composition of olivine can vary, the stoichiometry of equation (1) may vary. An example reaction equation for serpentinization is provided as follows:

The product of mixing water with olivine (left side of equation 2) is serpentine, magnetite, brucite, and hydrogen (right side of equation 2). The chemical reaction that produces natural hydrogen may depend on at least one of the chemical composition of the rock, access to a water source, chemical composition of the water. The chemical reaction may occur within an optimal temperature range of about 160° C. or more to about 450° C., such as about 210° C. or more to about 400° C. or less, such as about 300° C. An example of the modeled reaction is illustrated incomparing compositional variations between peridotiteand pyroxenitesource rocks. Specifically,illustrates a predicted equilibrium of hydrogen concentrations for serpentinization of lithologies having distinct olivine to orthopyroxene mass ratios as a function of temperature.

In some cases, hydrogen can be generated by hydration reactions in banded iron formations. Banded iron formations (BIFs) are massive chemical deposits composed of alternating layers of chert and iron-rich minerals (e.g., hematite, magnetite, and siderite). An intermediate reaction of iron (II) and water in banded iron formations is provided for ferric oxides (e.g., hematite and magnetite) and iron carbonate (e.g., siderite) BIFs and may yield free hydrogen as follows:

The temperature ranges may define the depths at which the reactions may occur, which may be determined by the thermal gradient, a function of the thermal conductivity of the rock and regional heat flow data. Rocks containing radioactive materials that are undergoing radioactive decay may have higher thermal conductivity. Rocks adjacent to magmatic sources such as volcanic intrusions may have higher heat flow. An example of temperature information acquired from wellbore data is illustrated in. Specifically,illustrates temperature and thermal conductivity in welland well. The linesrepresent thermal conductivity, black circlesare corrected bottom hole temperature readings, and the black linesare a continuous temperature values.

For oceanic settings such as mid-ocean ridges, the presence of ultramafic rock combined with high heat flow and constant water supply may provide suitable conditions for serpentinization reactions. For onshore conditions, the appropriate source rock types may occur at the surface or at depth, which may be determined from well bore data (e.g., wireline logs, core, drilling data) or geophysical data (e.g., seismic, gravity, and magnetic data). For onshore conditions with low geothermal gradients, the reaction temperatures may be unlikely to occur at the surface or at shallow depths (e.g., as shown in). However, some onshore geologic locations may have higher geothermal gradients leading to higher temperatures encountered at shallower depths. Hydration reactions in hydrogen source rocks will have their own kinetic characteristics whereby generation rates are controlled by pressure and temperature. For example, in a common example of a serpentinite formation, temperature dependence may be significant while pressure dependence may be negligible for the main structural types (e.g., lizardite, chrysotile, polygonal serpentine, and antigorite), as illustrated in the example phase diagram of.specifically illustrates a phase diagram having hydrogen generating reactions on the left side, at temperatures less than about 450° C., such as about 400° C. Sub-vertical phase transition boundaries may indicate minimal pressure dependence. A may represent antigorite. B may represent brucite. C may represent chlorite. Chr may represent chrysostile. L may represent lizardite. F may represent forsterite. T may represent talc, and W may represent water.

Aspects of the present disclosure provide methods and systems for utilizing earthquake data from seismic data repositories to determine active hydrogen resources. In at least one embodiment, methods described herein may utilize properties of faults within active seismic zones to estimate natural hydrogen generation potential resulting from water-rock interactions. The methods presented herein may assume hydrogen generation by hydration reactions related to Fe-bearing (ferrous) source rocks. In at least one embodiment, the natural hydrogen generation potential may represent a production prediction for use in hydrogen resource exploitation, either from natural hydrogen or from stimulated hydrogen production. The production prediction estimate described in at least one embodiment may be used in decisions regarding field development and may be used to define various specifications for the field operations associated with hydrogen extraction. In at least one embodiment, the production prediction may be a volume value (e.g., a static value) or a flux value (e.g., a dynamic value) of hydrogen resulting from natural generation as described below.

illustrates an example flow diagram of a method for predicting natural hydrogen production using fault parameters data, which may be performed according to at least one embodiment of the present disclosure. The one or more blocks ofmay be performed by a computer system or device (e.g., the computer systems and devices described below with respect to), including by a computer processor and a network communication interface. At step, the system or device obtains one or more hydrogen prospects, the one or more hydrogen prospects selected based on one or more parameters. In at least one embodiment, the one or more hydrogen prospects are selected based on at least one of rock characteristics of the one or more hydrogen prospects, fault characteristics of the one or more hydrogen prospects, and earthquake catalog data of the one or more hydrogen prospects. In at least one embodiment, the one or more hydrogen prospects represent fractured subsurface reservoirs connected to active fault systems.

At step, the system or device construct one or more fracture networks based on the one or more hydrogen prospects, the one or more fracture networks each having source rock area value. In at least one embodiment, constructing the one or more fracture networks based on the one or more hydrogen prospects includes defining a set of subsurface fault parameters for each of the one or more hydrogen prospects, the set of subsurface fault parameters characterizing an active seismic zone, and delineating, based on a comparison of the set of subsurface fault parameters and the active seismic zone, an active fault volume. In at least one embodiment, constructing the one or more fracture networks based on the one or more hydrogen prospects includes identifying a set of connections between the active fault volume and at least one water source, and identifying, based on the set of connections, water source flux. In at least one embodiment, constructing the one or more fracture networks based on the one or more hydrogen prospects includes forming a source rock exposure surface, and generating an area of the source rock exposure surface based on a framework model of the active fault volume. In at least one embodiment, forming the source rock exposure surface includes comparing the active fault volume to a source rock volume, and identifying exposure points on the active fault volume to generate the source rock exposure surface. In at least one embodiment, the framework model is generated using at least a statistical distribution of earthquakes. In at least one embodiment, the framework model is generated using a discrete fracture network (DFN) model. In at least one embodiment, the DFN model is constructed based on geological data. In at least one embodiment, the framework model is generated using equivalent volume data.

At step, the system or device, based on the source rock area for each of the one or more hydrogen prospects, generates a source rock exposure value for each of the one or more prospects. At step, the system or device outputs a production prediction based on the source rock exposure value for each of the one or more hydrogen prospects. In at least one embodiment, the production prediction estimates a natural hydrogen gas volume accumulated in the one or more hydrogen prospects. In at least one embodiment, outputting a production prediction based on the source rock exposure value for each of the one or more hydrogen prospects includes calculating the production prediction based on at least one of the area of the source rock exposure surface and the water source flux. In an optional step, the system or device may extract, based on the production prediction, natural hydrogen from a fractured subsurface reservoir represented by at least one of the one or more hydrogen prospects.

illustrates an example flow diagram of a method for predicting natural hydrogen production using fault parameters data, which may be performed according to at least one embodiment of the present disclosure. The one or more blocks ofmay be performed by a computer system or device (e.g., the computer systems and devices described below with respect to), including by a computer processor and a network communication interface. At step, a system or device performs a regional exploration screening for natural hydrogen play potential using geological maps, source rock data, regional fault maps, and other relevant screening data. At step, the system or device estimates the preliminary hydrogen production potential using analog data or analytical calculations. At step, the system or device uses the estimations of stepto determine if the potential hydrogen play may warrant further study. At step, the system or device reads a set of subsurface fault parameters characterizing an active seismic zone and delineates an active fault volume based on the subsurface fault parameters relationship with the active seismic zone. At step, the system or device determines whether the active fault volume connects to a water source to access a water source flux, and to determine if the active fault volume is exposed to a source rock to form a source rock exposure.

At step, the system or device calculates the hydrogen generation potential for a field scenario. In at least one embodiment, an accurate calculation may be generated by computing the total surface area of the source rock that is exposed to water over time in order to solve the hydrothermal hydration reaction described above. The total surface area may be derived from framework models of the fault network via at least one of the following: (1) Using the statistical distribution combining the theoretical magnitude to length relationship (e.g., as described with respect to) and real data (e.g., as described with respect to); (2) Building a discrete fracture network (DFN) model conditioned to the minimum available geologic data (e.g., as described with respect to); (3) Building a DFN model conditioned to more robust geological data by combiningwith geologic descriptions such as those describe with respect to, which may be derived from maps, three-dimensional (3D) seismic data, well data, and the like; (4) Building an equivalent volume data that incorporates a natural fracture network based on other forms of fault (e.g., side-wall leak-off, finite width fault core, fault splays, finite width damage zone) comprising an equivalent reaction zone (e.g., as described with respect to). Where available, supporting information from microseismic monitoring may be used to further constrain the predicted extent of faults at the depth of interest and in contact with source rock interval. Using one of the above derivation methods (1)-(4) may allow for constrained reaction predictions.

At step, the system or device calculates a production prediction of the natural hydrogen gas accumulated in a subsurface reservoir connected to the active fault volume based on the area of the source rock exposure and on the water source flux. In at least one embodiment, the computed hydrogen productivity in the cumulative equivalent area (e.g., as described with respect to) or cumulative equivalent volume (e.g., as described with respect to) may be determined using reaction equations such as equations (1)-(5). In one non-limiting example, the curves illustrated inmay show hydrogen generation in mmol/kg, and may be predicted using computer modeling techniques and adapted and scaled for large reaction volumes. In at least one embodiment, the chemical reactions may be solved computationally and linked to dynamic flow construction modules. In at least one embodiment, the reactions and fluid flux may be computed independently and passed as input/output values. In at least one embodiment, the reactions may be determined from laboratory tests and treated empirically. In at least one embodiment, the essential workflow includes determination of the surface area of active faults, which may be represented as cumulative equivalent area or cumulative equivalent volume, and then to solve the hydration reactions given the hydrothermal environmental conditions. A simplified version of this reaction concept may be depicted and described with respect to.

are example models of natural hydrogen reservoirs occurring in active seismic zones. Such hydrogen reservoirs may be targeted in the method of.illustrate the relationship between active faults, inactive faults, water circulation, ultramafic source rocks, serpentinization reactions, and hydrogen migration paths. In the scenarios presented in, the cause of the earthquakes may be unknown. However, active portions of the fault zone resulting from the earthquake may preferentially intake a higher volume of water than impermeable or low permeability host rock that may surround the fault zone.depicts an example fluid flow localized within a narrow fault zone. Water-rock interactions may be limited to a narrow zone surrounding the main fault. Water may enter the fault zone by deep connectivity with one of the following (a) meteoric water at the ground surface; (b) present day buried aquifer; (c) a paleo-aquifer that collected water during a previous geologic period including but not limited to periods of glacial melting.depicts a similar example fault zone that is connected to a fracture network on one side of the fault. Activation of the fault zone may allow a greater-than-average volume of water to enter the host rock volume through connected fracture flow.reflects a minimum volume scenario.reflects a moderate volume scenario based on asymmetric fracture connectivity. Other scenarios may be envisioned with different patterns of fault and fracture connectivity, according to certain aspects of the present disclosure.

In some cases, water source flux may be present if the fault zone penetrates a known source of water (e.g., an aquifer at a certain stratigraphic depth, ground water when a fault extends to the surface). In at least one embodiment of the present disclosure, a water source is at or near the surface, the presence and amount of water flux may be determined using ground water recharge data such as information about rainfall or other hydrogeologic data. Source rock exposure is present if the fault zone penetrates the source rock geological formation as presented in. The source rock exposure may be defined on the structural rupture area of the active faults.

In at least one embodiment, earthquake data may be incorporated into systems and methods described herein to provide information on the location, depth, and frequency of active faulting. Enhanced water circulation may be present in seismogenic terranes hydraulically connecting rocks from basement to surface along active fault zones. Earthquake data may be historic and/or real-time earthquake data. An example of earthquake data is shown in. The data may contain dates, times, geospatial coordinates, estimated depth(s) (shown in km), event magnitude, other seismological attributes, or some combination therein. Some events, such as the eventin, include strike information of the activated fault plane and dip information of the activated fault plane (e.g., depicted as ‘beachball’ plot). In some cases, earthquake data may indicate that earthquake magnitude is related to fault size. For example, as illustrated in, the rupture area may reflect a portion of the total fault surface given length, L, and width, W, that slips by some distance, D.illustrates a cross-section showing a fault area and the values that are used to compute the seismic moment. In some examples, fault patch may slip according to a non-linear slip distribution and may have an irregular shape. In such non-linear cases, the size of the fault as defined by a maximum length dimension may be approximated by assuming a square, circular, or elliptical patch size.

illustrates an example relationship between earthquake magnitude, moment (energy release), slip magnitude, and fault size. In one example, the black linehighlights an example fault of magnitudewith several centimeters of slip on a fault having a length of about 0.5 km or more to about 4.5 km or less (e.g., about 1 km to about 4 km). The relationships illustrated inmay provide estimation parameters for upper and lower bounds of the size of a given fault using earthquake magnitude data. Accordingly, using a larger population of earthquakes, it may be possible to collect statistical information on the range of earthquake magnitudes and corresponding fault lengths. In one example, a map-based method for linking mapped fault lengths with potential earthquake magnitude may be utilized to compute earthquake magnitudes derived from surface fault lengths (FLEM) as follows:

In some cases, pre-existing and verified empirical relationships may also be used to compute fault area from an earthquake catalog data. For example,illustrates an example relationship between strike-slip earthquake data and seismic moment. In at least one embodiment, the cumulative fault length or area for active faults can therefore be constrained using statistical information from earthquake data, like the information from. In at least one embodiment, mechanical models may also be utilized to validate assumptions about active fault trends or seismogenic locations. For example, fault and fracture stability computations may be employed to understand the relationship between tectonic stresses, faults, and borehole data in a variety of hydrogen formations. In at least one embodiment, if map data is available, cumulative fault length (or area) may be compared to ensure an active fault length for a region is less than the total mapped faults. In at least one embodiment, fault area may be computed by assuming a height to length ratio for rectangular or square fault patches or using a circular or elliptical fault patch shape assumption.

In some cases, fault statistics may be obtained from earthquake catalog data. Example catalog data is shown in.shows an example cumulative number of earthquakes as a function of earthquake magnitude. The magnitude of completeness Mc (e.g., about 2.70) for this data set may indicate the minimum earthquake size that is accurately detected by the seismic monitoring network. Linemarking the slope of the curve reflects the b-value (e.g., about 0.94±0.05), which is a calculation of the scaling ratio between small and large earthquakes. In at least one embodiment, using the relationships into convert from magnitude to fault length, it may be possible to determine the statistical size distribution for this active fault population.is a histogram of the earthquake data of,showing the statistical distribution of hypocenter depth. Information about the ranges of earthquake depth is necessary to determine if there is a nearby source of water to infiltrate the fault zones and host the hydration reaction process. Earthquake depth may also be related to the geothermal gradient, which determines whether active faults exposed to an influx of water can carry out the hydration reaction within the optimal temperature ranges.

illustrates a discrete fracture network representation of resolved active faults represented as equivalent area representations. In at least one embodiment, the exposed surface area for each discrete fault is proportionate to potential hydrogen reactivity. Outputshows one possible realization of a fault plane with connected fractures generated from the active faults.illustrates a discrete fracture network representation of resolved active faults represented as equivalent volume representations. In at least one embodiment, the exposed volume is proportionate to potential hydrogen reactivity. The volume may be determined by the area of each fault and any connected fractures that facilitate leak-off into or water influx from the host rock.illustrates an example concept for natural hydrogen production in faulted basement source rock.

The methods described herein may be implemented on a computer system. In at least embodiment, the methods may be implemented on computer systemof. As shown in, the exploration and production computer systemincludes an analysis tool, a data repositoryfor storing input data, intermediate data, and resultant outputs of the analysis tool, and a user interface. The computer systemis connected to a networkand may access earthquake data from the earthquake data repository. The data repositorymay include a repositoryfor subsurface fault parameters, a repositoryfor active seismic zones, a repositoryfor active fault volumes, a repositoryfor water source fluxes, a repositoryfor source rock exposure areas, a repositoryfor framework models, and a repositoryfor production predictions. The analysis toolmay include a modulefor active seismic zone identification, a modulefor active fault volume building, a modulefor water source connection identification, a modulefor source rock exposure identification, a modulefor source rock exposure area calculation, and a modulefor production prediction calculation.

The computer systemmay be part of, integrated with, or independent from any of the features described with respect to,, and.

is an example of a block diagram of a systemfor natural hydrogen exploration using fault data. The systemcan be implemented using one or more modules, shown in block form in the drawings. The one or more modules can be in software or hardware form, or a combination thereof. In some examples, the systemcan be implemented as machine readable instructions for execution on one or more computing platforms(referred to as a computing platform herein), as shown in. The computing platformcan include one or more computing devices selected from, for example, a desktop computer, a server, a controller, a blade, a mobile phone, a tablet, a laptop, a personal digital assistant (PDA), and the like.

The computing platformcan include a processorand a memory. By way of example, the memorycan be implemented, for example, as a non-transitory computer storage medium, such as volatile memory (e.g., random access memory), non-volatile memory (e.g., a hard disk drive, a solid-state drive, a flash memory, or the like), or a combination thereof. The processorcan be implemented, for example, as one or more processor cores. The memorycan store machine-readable instructions that can be retrieved and executed by the processorto implement the system. Each of the processorand the memorycan be implemented on a similar or a different computing platform. The computing platformcan be implemented in a cloud computing environment (for example, as disclosed herein) and thus on a cloud infrastructure. In such a situation, features of the computing platformcan be representative of a single instance of hardware or multiple instances of hardware executing across the multiple of instances (e.g., distributed) of hardware (e.g., computers, routers, memory, processors, or a combination thereof). Alternatively, the computing platformcan be implemented on a single dedicated server or workstation.

In view of the structural and functional features described above, example methods will be better appreciated with reference to. While, for purposes of simplicity of explanation, the example methods ofare shown and described as executing serially, it is to be understood and appreciated that the present examples are not limited by the illustrated order, as some actions could in other examples occur in different orders, multiple times and/or concurrently from that shown and described herein. Moreover, it is not necessary that all described actions be performed to implement the methods, and conversely, some actions may be performed that are omitted from the description.

In view of the foregoing structural and functional description, those skilled in the art will appreciate that portions of the embodiments may be embodied as a method, data processing system, or computer program product. Accordingly, these portions of the present embodiments may take the form of an entirely hardware embodiment, an entirely software embodiment, or an embodiment combining software and hardware, such as shown and described with respect to the computer system of. Furthermore, portions of the embodiments may be a computer program product on a computer-readable storage medium having computer readable program code on the medium. Any non-transitory, tangible storage media possessing structure may be utilized including, but not limited to, static and dynamic storage devices, volatile and non-volatile memories, hard disks, optical storage devices, and magnetic storage devices, but excludes any medium that is not eligible for patent protection under 35 U.S.C. § 111 (such as a propagating electrical or electromagnetic signals per se). As an example and not by way of limitation, computer-readable storage media may include a semiconductor-based circuit or device or other IC (such, as for example, a field-programmable gate array (FPGA) or an ASIC), a hard disk, an HDD, a hybrid hard drive (HHD), an optical disc, an optical disc drive (ODD), a magneto-optical disc, a magneto-optical drive, a floppy disk, a floppy disk drive (FDD), magnetic tape, a holographic storage medium, a solid-state drive (SSD), a RAM-drive, a SECURE DIGITAL card, a SECURE DIGITAL drive, or another suitable computer-readable storage medium or a combination of two or more of these, where appropriate. A computer-readable non-transitory storage medium may be volatile, nonvolatile, or a combination of volatile and non-volatile, as appropriate.

Certain embodiments have also been described herein with reference to block illustrations of methods, systems, and computer program products. It will be understood that blocks and/or combinations of blocks in the illustrations, as well as methods or steps or acts or processes described herein, can be implemented by a computer program comprising a routine of set instructions stored in a machine-readable storage medium as described herein. These instructions may be provided to one or more processors of a general purpose computer, special purpose computer, or other programmable data processing apparatus (or a combination of devices and circuits) to produce a machine, such that the instructions of the machine, when executed by the processor, implement the functions specified in the block or blocks, or in the acts, steps, methods and processes described herein.

These processor-executable instructions may also be stored in computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory result in an article of manufacture including instructions which implement the function specified. The computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to realize a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in flowchart blocks that may be described herein.

In this regard,illustrates one example of a computer systemthat can be employed to execute one or more embodiments of the present disclosure. Computer systemcan be implemented on one or more general purpose networked computer systems, embedded computer systems, routers, switches, server devices, client devices, various intermediate devices/nodes or standalone computer systems. Additionally, computer systemcan be implemented on various mobile clients such as, for example, a personal digital assistant (PDA), laptop computer, pager, and the like, provided it includes sufficient processing capabilities.

Computer systemincludes processing unit, system memory, and system busthat couples various system components, including the system memory, to processing unit. System memorycan include volatile (e.g. RAM, DRAM, SDRAM, Double Data Rate (DDR) RAM, etc.) and non-volatile (e.g. Flash, NAND, etc.) memory. Dual microprocessors and other multi-processor architectures also can be used as processing unit. System busmay be any of several types of bus structure including a memory bus or memory controller, a peripheral bus, and a local bus using any of a variety of bus architectures. System memoryincludes read only memory (ROM)and random access memory (RAM). A basic input/output system (BIOS)can reside in ROMcontaining the basic routines that help to transfer information among elements within computer system.

Computer systemcan include a hard disk drive, magnetic disk drive, e.g., to read from or write to removable disk, and an optical disk drive, e.g., for reading CD-ROM diskor to read from or write to other optical media. Hard disk drive, magnetic disk drive, and optical disk driveare connected to system busby a hard disk drive interface, a magnetic disk drive interface, and an optical drive interface, respectively. The drives and associated computer-readable media provide nonvolatile storage of data, data structures, and computer-executable instructions for computer system. Although the description of computer-readable media above refers to a hard disk, a removable magnetic disk and a CD, other types of media that are readable by a computer, such as magnetic cassettes, flash memory cards, digital video disks and the like, in a variety of forms, may also be used in the operating environment; further, any such media may contain computer-executable instructions for implementing one or more parts of embodiments shown and described herein.

A number of program modules may be stored in drives and RAM, including operating system, one or more application programs, other program modules, and program data. In some examples, the application programscan include a modulefor active seismic zone identification, a modulefor active fault volume building, a modulefor water source connection identification, a modulefor source rock exposure identification, a modulefor source rock exposure area calculation, and a modulefor production prediction calculation, and the program datacan include a repositoryfor subsurface fault parameters, a repositoryfor active seismic zones, a repositoryfor active fault volumes, a repositoryfor water source fluxes, a repositoryfor source rock exposure areas, a repositoryfor framework models, and a repositoryfor production predictions. The application programsand program datacan include functions and methods programmed to enhance natural hydrogen exploration using fault data, such as shown and described herein.